Method

NumCosmoGalaxySDObsRedshiftGausscompute_equal_area_photoz_bins

Declaration [src]

NcmVector*
nc_galaxy_sd_obs_redshift_gauss_compute_equal_area_photoz_bins (
  NcGalaxySDObsRedshiftGauss* gsdorgauss,
  guint n_bins,
  gdouble zp_max
)

Description [src]

Computes equal-area photometric redshift bin edges.

This function computes bin edges with equal areas under the convolved photo-z distribution, suitable for weak lensing source samples where equal-area binning in photo-z space is desired.

The algorithm:

  1. Computes the convolved photo-z distribution $P(z_p) = \int P(z) \times \mathrm{Gauss}(z_p|z,\sigma_z) \, \mathrm{d}z$
  2. Builds the cumulative distribution function (CDF)
  3. Inverts the CDF to find n_bins+1 edges with equal integrated probability.

Parameters

n_bins

Type: guint

Number of bins to create.

zp_max

Type: gdouble

Maximum photometric redshift to consider.

Return value

Type: NcmVector

A NcmVector with n_bins+1 photo-z bin edges.

The caller of the method takes ownership of the returned data, and is responsible for freeing it.