onepower.pk.Spectra.fg_fit

Spectra.fg_fit(mass, mb, fstar, beta=2)

Compute the gas fraction for a general baryonic feedback model, 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