hmf.core._kernels.fits.behroozi_modify_dndm

hmf.core._kernels.fits.behroozi_modify_dndm(m_msun: ArrayLike, dndm: ArrayLike, z: ArrayLike, ngtm: ArrayLike, h: ArrayLike) → ndarray[tuple[Any, ...], dtype[float64]][source]

The Behroozi et al. (2013) App. G correction to the Tinker (2008) mass function.

\[\frac{dn}{dM} = \theta\frac{dn_{\rm T08}}{dM} - n_{\rm T08}(>M)\frac{d\theta}{dM}, \qquad \theta = 10^{\alpha(M/M_\star)^\gamma},\]

with \(\alpha = 0.144/(1 + \exp[14.79(a - 0.213)])\), \(\gamma = 0.5/(1 + \exp(6.5a))\), \(a = 1/(1+z)\) and \(M_\star = 10^{11.5}\,M_\odot\). This is the derivative of \(n(>M) = \theta\, n_{\rm T08}(>M)\).

Parameters:
  • m_msun – Halo masses in Msun (not Msun/h: the pivot mass is in Msun).

  • dndm – The Tinker (2008) mass function, per unit mass in Msun/h (canonical h^4 / (Msun Mpc^3)).

  • z – Redshift.

  • ngtm – The cumulative Tinker (2008) mass function n(>M), in h^3 / Mpc^3.

  • h – The dimensionless Hubble parameter, used to take dtheta/dM per unit Msun/h, consistently with dndm.

Returns:

The corrected dn/dM, in the units of dndm. Non-finite inputs give non-finite outputs: nothing is replaced by a sentinel.

Return type:

numpy.ndarray