massfunctions.py 文件源码

def quenchedfrac_zfourge(M, z):
    par_q1, par_q2 = zfourgeparams(z, type='quiescent')
    x_q1 = 10.**(M-par_q1[1])
    dn_q1 = np.log10(np.log(10)*np.exp(-1.*x_q1)*x_q1*(10.**par_q1[3]*x_q1**(par_q1[2]) + 10.**(par_q1[5])*x_q1**(par_q1[4])))
    x_q2 = 10.**(M-par_q2[1])
    dn_q2 = np.log10(np.log(10)*np.exp(-1.*x_q2)*x_q2*(10.**par_q2[3]*x_q2**(par_q2[2]) + 10.**(par_q2[5])*x_q2**(par_q2[4])))

    par_sf1, par_sf2 = zfourgeparams(z, type='star-forming')
    x_sf1 = 10.**(M-par_sf1[1])
    dn_sf1 = np.log10(np.log(10)*np.exp(-1.*x_sf1)*x_sf1*(10.**par_sf1[3]*x_sf1**(par_sf1[2]) + 10.**(par_sf1[5])*x_sf1**(par_sf1[4])))
    x_sf2 = 10.**(M-par_sf2[1])
    dn_sf2 = np.log10(np.log(10)*np.exp(-1.*x_sf2)*x_sf2*(10.**par_sf2[3]*x_sf2**(par_sf2[2]) + 10.**(par_sf2[5])*x_sf2**(par_sf2[4])))

    fq1 = 10.**dn_q1/(10.**dn_q1+10.**dn_sf1)
    fq2 = 10.**dn_q2/(10.**dn_q2+10.**dn_sf2)

    return (fq1*(par_q2[0]-z)+fq2*(z-par_q1[0]))/(par_q2[0]-par_q1[0])

# ------ OBSOLETE, left in for backwards-compatibility ------ #