product of the bethe lattice dos, fermi distribution

def bethe_fermi(energy, quasipart, shift, hopping, beta):
    """product of the bethe lattice dos, fermi distribution"""
    return fermi_dist(quasipart * energy - shift, beta) \
            * bethe_lattice(energy, hopping)