dsf.modelling.hod module

Halo occupation distribution power spectrum models.

dsf.modelling.hod.hod_galaxy_bias(cosmo, *, a, concentration=None, **hod_kwargs)[source]

Calculate the occupation-weighted linear galaxy bias.

The effective bias is computed by averaging the halo bias over the central and satellite occupations, weighted by the halo mass function.

Parameters:

  • cosmo (Cosmology) – CCL cosmology.

  • a (float) – Scale factor.

  • concentration (Concentration | None) – Optional CCL concentration model. Defaults to Duffy08. The model must be compatible with the 200c halo mass definition.

  • **hod_kwargs (Any) – HOD profile parameters.

Returns:

Effective linear galaxy bias.

Raises:

ValueError – If the scale factor is invalid or the galaxy normalization is nonfinite or nonpositive.

Return type:

float

dsf.modelling.hod.pk2d_hod(cosmo, *, k_array, a_array, concentration=None, **hod_kwargs)[source]

Construct standard HOD power spectra.

Parameters:

  • cosmo (Cosmology) – CCL cosmology.

  • k_array (NDArray[float64]) – Wavenumber grid.

  • a_array (NDArray[float64]) – Scale-factor grid.

  • concentration (Concentration | None) – Optional CCL concentration model. Defaults to Duffy08. The model must be compatible with the 200c halo mass definition.

  • **hod_kwargs (Any) – HOD profile parameters.

Returns:

Dictionary containing pk_deltasigma, pk_gm, and pk_gg.

Return type:

dict[str, Pk2D]