def radio_lumfn(L, params, z=0.): """ Calculate radio luminosity function for a given set of model parameters. """ # Number density as a function of sfr, dn/dlog(sfr) sfr = L * 5.52e-29 # erg/s/Hz, Bell (2003), Eq. 6 dndlogsfr_sfms, dndlogsfr_pass = g.sfr_fn(hm, sfr, z=z, params=params) phi = dndlogsfr_sfms #+ dndlogsfr_pass # FIXME: Remove passive galaxies return phi
def radio_lumfn(L, _params): """ Calculate radio luminosity function for a given set of model parameters. """ # Number density as a function of sfr, dn/dlog(sfr) sfr = L * 5.52e-29 # erg/s/Hz, Bell (2003), Eq. 6 dndlogsfr_sfms, dndlogsfr_pass = g.sfr_fn(hm, sfr, z=0., params=_params) #phi = dndlogsfr_sfms #+ dndlogsfr_pass return dndlogsfr_sfms, dndlogsfr_pass