def METHOD_NAME(xyzfile, charge, spin):
    '''
    Assemble a molecule object from xyzfile, charge and spin.
    '''
    wholetext = "%s %s\n" % (charge, spin)
    if os.path.isfile(xyzfile):
        with open(xyzfile, "r") as fo:
            natoms = int(fo.readline().split()[0])
            fo.readline()
            for ii in range(natoms):
                wholetext += fo.readline()
    wholetext += "\nsymmetry c1\nnoreorient\nnocom\n"
    mol = psi4.geometry("""%s""" % wholetext)
    return mol