onepower.pk.Spectra.fg

Spectra.fg(mass, z_vec, fstar, beta=2)

Compute the gas fraction, eq 24 of Mead et al. 2021.

\(f_{\rm g}(M) = [\Omega_{\rm b}/\Omega_{\rm m} - f_{\star}] (M/M_{\rm b})^{\beta}/ (1 + (M/M_{\rm b})^{\beta})\)

where \(f_{\rm g}\) is the halo gas fraction, the pre-factor in parenthesis is the available gas reservoir, while \(M_{\rm b}\) > 0` and \(\beta > 0\) are fitted parameters. Haloes of \(M >> M_{\rm b}\) are unaffected while those of \(M < M_{\rm b}\) have lost more than half of their gas

theta_agn = log10_TAGN - 7.8 table 4 of 2009.01858, units of M_sun/h