Cluster Abundance Distribution

Cluster Abundance Distribution — FIXME

Synopsis

struct              NcClusterAbundance;
struct              NcClusterAbundanceClass;
                    NcClusterAbundanceDataBin;
                    NcClusterAbundanceDataBinM;
typedef             NcClusterAbundanceDataBinZ;
                    NcClusterAbundanceDataP;
gdouble             (*NcClusterAbundanceN)              (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMl,
                                                         gdouble lnMu,
                                                         gdouble zl,
                                                         gdouble zu);
enum                NcClusterAbundanceOpt;
gdouble             (*NcClusterAbundanced2NdzdlnM)      (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnM,
                                                         gdouble z);
gdouble             (*NcClusterAbundancedNdlnM)         (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnM,
                                                         gdouble zl,
                                                         gdouble zu);
gdouble             (*NcClusterAbundancedNdz)           (NcClusterAbundance *cad,
                                                         NcHICosmo *godel,
                                                         gdouble lnMl,
                                                         gdouble lnMu,
                                                         gdouble z);
void                nc_bias_mean_prepare                (NcClusterAbundance *cad,
                                                         NcHICosmo *model);
gdouble             nc_bias_mean_val                    (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMl,
                                                         gdouble lnMu,
                                                         gdouble z);
gdouble             nc_ca_mean_bias                     (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnM,
                                                         gdouble z);
gdouble             nc_ca_mean_bias_Mobs_denominator    (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble z);
gdouble             nc_ca_mean_bias_Mobs_numerator      (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble z);
gdouble             nc_ca_mean_bias_denominator         (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnM,
                                                         gdouble z);
gdouble             nc_ca_mean_bias_numerator           (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnM,
                                                         gdouble z);
gdouble             nc_cluster_abundance_N_Mobs_local_selection
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs_i,
                                                         gdouble lnMobs_f,
                                                         gdouble z_i,
                                                         gdouble z_f);
gdouble             nc_cluster_abundance_N_photoz_Mobs_local_selection
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs_i,
                                                         gdouble lnMobs_f,
                                                         gdouble zp_i,
                                                         gdouble zp_f);
#define             nc_cluster_abundance_N_val          (cad,
                                                         cp,
                                                         lnMl,
                                                         lnMu,
                                                         zl,
                                                         zu)
void                nc_cluster_abundance_and_options    (NcClusterAbundance *cad,
                                                         NcClusterAbundanceOpt opt);
void                nc_cluster_abundance_bin_realization
                                                        (GArray *zr,
                                                         gsl_histogram **h);
NcClusterAbundance * nc_cluster_abundance_copy          (NcClusterAbundance *cad);
gdouble             nc_cluster_abundance_d2NdzdlnM_Mobs (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble z);
gdouble             nc_cluster_abundance_d2NdzdlnM_Mobs_local_selection
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble z);
gdouble             nc_cluster_abundance_d2NdzdlnM_photoz
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnM,
                                                         gdouble z_phot);
gdouble             nc_cluster_abundance_d2NdzdlnM_photoz_Mobs
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble zp);
gdouble             nc_cluster_abundance_d2NdzdlnM_photoz_Mobs_local_selection
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble zp);
gdouble             nc_cluster_abundance_d2NdzdlnM_purity_val
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble z);
#define             nc_cluster_abundance_d2NdzdlnM_val  (cad,
                                                         cp,
                                                         lnM,
                                                         z)
#define             nc_cluster_abundance_dNdlnM_val     (cad,
                                                         cp,
                                                         lnM,
                                                         zl,
                                                         zu)
gdouble             nc_cluster_abundance_dNdz_Mobs_local_selection
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs_i,
                                                         gdouble lnMobs_f,
                                                         gdouble z);
gdouble             nc_cluster_abundance_dNdz_photoz_Mobs_local_selection
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs_i,
                                                         gdouble lnMobs_f,
                                                         gdouble z);
gdouble             nc_cluster_abundance_dNdz_purity_val
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMl,
                                                         gdouble lnMu,
                                                         gdouble z);
#define             nc_cluster_abundance_dNdz_val       (cad,
                                                         cp,
                                                         lnMl,
                                                         lnMu,
                                                         z)
void                nc_cluster_abundance_free           (NcClusterAbundance *cad);
gdouble             nc_cluster_abundance_get_lnM_sigma0 (NcClusterAbundance *cad);
gdouble             nc_cluster_abundance_get_lnMf       (NcClusterAbundance *cad);
gdouble             nc_cluster_abundance_get_lnMi       (NcClusterAbundance *cad);
NcClusterPhotoz *   nc_cluster_abundance_get_photoz     (NcClusterAbundance *cad);
gdouble             nc_cluster_abundance_get_zf         (NcClusterAbundance *cad);
gdouble             nc_cluster_abundance_get_zi         (NcClusterAbundance *cad);
NcClusterAbundance * nc_cluster_abundance_new           (NcClusterAbundanceOpt opt,
                                                         NcMassFunction *mfp,
                                                         NcHaloBiasFunc *mbiasf,
                                                         gdouble zi,
                                                         gdouble zf,
                                                         gdouble lnMi,
                                                         gdouble lnMf,
                                                         NcClusterPhotoz *photoz,
                                                         gdouble lnM_sigma0);
void                nc_cluster_abundance_prepare        (NcClusterAbundance *cad,
                                                         NcHICosmo *model);
void                nc_cluster_abundance_prepare_inv_dNdlnM_z
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble z);
void                nc_cluster_abundance_prepare_inv_dNdlnM_z_no_obs
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble z);
void                nc_cluster_abundance_prepare_inv_dNdz
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model);
void                nc_cluster_abundance_prepare_inv_dNdz_no_obs
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model);
GPtrArray *         nc_cluster_abundance_realizations_read_from_file
                                                        (gchar *file_realization,
                                                         gint n_realizations);
void                nc_cluster_abundance_realizations_save_to_file
                                                        (GPtrArray *realizations,
                                                         gchar *filename);
void                nc_cluster_abundance_set_lnM_sigma0 (NcClusterAbundance *cad,
                                                         const gdouble lnM_sigma0);
void                nc_cluster_abundance_set_lnMf       (NcClusterAbundance *cad,
                                                         const gdouble lnMf);
void                nc_cluster_abundance_set_lnMi       (NcClusterAbundance *cad,
                                                         const gdouble lnMi);
void                nc_cluster_abundance_set_options    (NcClusterAbundance *cad,
                                                         NcClusterAbundanceOpt opt);
void                nc_cluster_abundance_set_photoz     (NcClusterAbundance *cad,
                                                         NcClusterPhotoz *photoz);
void                nc_cluster_abundance_set_zf         (NcClusterAbundance *cad,
                                                         const gdouble zf);
void                nc_cluster_abundance_set_zi         (NcClusterAbundance *cad,
                                                         const gdouble zi);

Object Hierarchy

  GObject
   +----NcClusterAbundance

Properties

  "lnM-sigma0"               gdouble               : Read / Write
  "mass-function"            NcMassFunction*       : Read / Write / Construct Only
  "maximum-mass"             gdouble               : Read / Write
  "maximum-redshift"         gdouble               : Read / Write
  "mean-bias"                NcHaloBiasFunc*       : Read / Write / Construct Only
  "minimum-mass"             gdouble               : Read / Write
  "minimum-redshift"         gdouble               : Read / Write
  "options"                  NcClusterAbundanceOpt  : Read / Write / Construct Only
  "photoz"                   NcClusterPhotoz*      : Read / Write / Construct Only

Description

FIXME

Details

struct NcClusterAbundance

struct NcClusterAbundance;

struct NcClusterAbundanceClass

struct NcClusterAbundanceClass {
};

NcClusterAbundanceDataBin

typedef struct _NcClusterAbundanceDataBin NcClusterAbundanceDataBin;

NcClusterAbundanceDataBinM

typedef struct _NcClusterAbundanceDataBinM NcClusterAbundanceDataBinM;

NcClusterAbundanceDataBinZ

typedef struct _NcClusterAbundanceDataBinz NcClusterAbundanceDataBinZ;

NcClusterAbundanceDataP

typedef struct _NcClusterAbundanceDataP NcClusterAbundanceDataP;

NcClusterAbundanceN ()

gdouble             (*NcClusterAbundanceN)              (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMl,
                                                         gdouble lnMu,
                                                         gdouble zl,
                                                         gdouble zu);

enum NcClusterAbundanceOpt

typedef enum {
  NC_CLUSTER_ABUNDANCE_NONE         = 0,
  NC_CLUSTER_ABUNDANCE_BINMASS      = 1 << 0,
  NC_CLUSTER_ABUNDANCE_BINZ         = 1 << 1,
  NC_CLUSTER_ABUNDANCE_POISSON      = 1 << 2,
  NC_CLUSTER_ABUNDANCE_PHOTOZ       = 1 << 3,
  NC_CLUSTER_ABUNDANCE_MOBS         = 1 << 4,
  NC_CLUSTER_ABUNDANCE_MOBS_LOCAL   = 1 << 5,
  NC_CLUSTER_ABUNDANCE_COMPLETENESS = 1 << 6,
  NC_CLUSTER_ABUNDANCE_PURITY       = 1 << 7,
  NC_CLUSTER_ABUNDANCE_REAL_ZM      = 1 << 8,
  NC_CLUSTER_ABUNDANCE_OBS_ZM       = 1 << 9,
  /* NC_MEAN_BIAS                   = 1 << 5, FIXME */
  /* NC_MEAN_BIAS_MOBS              = 1 << 6, FIXME */
} NcClusterAbundanceOpt;

FIXME

NC_CLUSTER_ABUNDANCE_NONE

without observational uncertainties

NC_CLUSTER_ABUNDANCE_BINMASS

FIXME

NC_CLUSTER_ABUNDANCE_BINZ

FIXME

NC_CLUSTER_ABUNDANCE_POISSON

FIXME

NC_CLUSTER_ABUNDANCE_PHOTOZ

include photometric redshift distribution (Gaussian)

NC_CLUSTER_ABUNDANCE_MOBS

include mass-observable distribution (log-normal)

NC_CLUSTER_ABUNDANCE_MOBS_LOCAL

include mass-observable distribution (log-normal) with standard deviation as a function of mass and redshift

NC_CLUSTER_ABUNDANCE_COMPLETENESS

include completeness

NC_CLUSTER_ABUNDANCE_PURITY

include purity

NC_CLUSTER_ABUNDANCE_REAL_ZM

FIXME

NC_CLUSTER_ABUNDANCE_OBS_ZM

FIXME

NcClusterAbundanced2NdzdlnM ()

gdouble             (*NcClusterAbundanced2NdzdlnM)      (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnM,
                                                         gdouble z);

NcClusterAbundancedNdlnM ()

gdouble             (*NcClusterAbundancedNdlnM)         (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnM,
                                                         gdouble zl,
                                                         gdouble zu);

NcClusterAbundancedNdz ()

gdouble             (*NcClusterAbundancedNdz)           (NcClusterAbundance *cad,
                                                         NcHICosmo *godel,
                                                         gdouble lnMl,
                                                         gdouble lnMu,
                                                         gdouble z);

nc_bias_mean_prepare ()

void                nc_bias_mean_prepare                (NcClusterAbundance *cad,
                                                         NcHICosmo *model);

nc_bias_mean_val ()

gdouble             nc_bias_mean_val                    (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMl,
                                                         gdouble lnMu,
                                                         gdouble z);

nc_ca_mean_bias ()

gdouble             nc_ca_mean_bias                     (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnM,
                                                         gdouble z);

nc_ca_mean_bias_Mobs_denominator ()

gdouble             nc_ca_mean_bias_Mobs_denominator    (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble z);

nc_ca_mean_bias_Mobs_numerator ()

gdouble             nc_ca_mean_bias_Mobs_numerator      (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble z);

nc_ca_mean_bias_denominator ()

gdouble             nc_ca_mean_bias_denominator         (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnM,
                                                         gdouble z);

nc_ca_mean_bias_numerator ()

gdouble             nc_ca_mean_bias_numerator           (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnM,
                                                         gdouble z);

nc_cluster_abundance_N_Mobs_local_selection ()

gdouble             nc_cluster_abundance_N_Mobs_local_selection
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs_i,
                                                         gdouble lnMobs_f,
                                                         gdouble z_i,
                                                         gdouble z_f);

This function computes the number of clusters with mass between [lnMobs_i, lnMobs_f] in redshift bins, i.e., N_i = /f$ \int_z_i^\z_f dz \int_{lnMobs_i}^{lnMobs_f} d\lnMobs 1/p(lnMobs, z) \int_0^\infty d\ln M \frac{d^2N(\ln M, z)}{dzd\ln M} * c(\lnM, z) * P(\ln M^{obs}|\ln M, z) /f$.

cad :

a NcClusterAbundance.

model :

a NcHICosmo.

lnMobs_i :

logarithm base e of the observable minimum mass threshold.

lnMobs_f :

logarithm base e of the observable maximum mass threshold.

z_i :

redshift, minimum threshold of a ith bin.

z_f :

redshift, maximum threshold of a ith bin.

Returns :

a gdouble which represents the number of clusters per redshift bin taking into account the completeness and purity.

nc_cluster_abundance_N_photoz_Mobs_local_selection ()

gdouble             nc_cluster_abundance_N_photoz_Mobs_local_selection
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs_i,
                                                         gdouble lnMobs_f,
                                                         gdouble zp_i,
                                                         gdouble zp_f);

This function computes the number of clusters with mass between [lnMobs_i, lnMobs_f] in redshift bins, i.e., N_i = /f$ \int_z_i^\z_f dz \int_{lnMobs_i}^{lnMobs_f} d\lnMobs 1/p(lnMobs, z) \int_0^\infty d\ln M \frac{d^2N(\ln M, z)}{dzd\ln M} * c(\lnM, z) * P(\ln M^{obs}|\ln M, z) * P(zphot | z)/f$.

cad :

a NcClusterAbundance.

model :

a NcmModel.

lnMobs_i :

logarithm base e of the observable minimum mass threshold.

lnMobs_f :

logarithm base e of the observable maximum mass threshold.

zp_i :

photometric redshift, minimum threshold of a ith bin.

zp_f :

photometric redshift, maximum threshold of a ith bin.

Returns :

a gdouble which represents the number of clusters per redshift bin taking into account the completeness and purity.

nc_cluster_abundance_N_val()

#define nc_cluster_abundance_N_val(cad,cp,lnMl,lnMu,zl,zu) (cad)->N_val(cad,cp,lnMl,lnMu,zl,zu)

nc_cluster_abundance_and_options ()

void                nc_cluster_abundance_and_options    (NcClusterAbundance *cad,
                                                         NcClusterAbundanceOpt opt);

nc_cluster_abundance_bin_realization ()

void                nc_cluster_abundance_bin_realization
                                                        (GArray *zr,
                                                         gsl_histogram **h);

FIXME

zr :

FIXME

h :

FIXME

nc_cluster_abundance_copy ()

NcClusterAbundance * nc_cluster_abundance_copy          (NcClusterAbundance *cad);

Duplicates the NcClusterAbundance object setting the same values of the original propertities.

cad :

a NcClusterAbundance.

Returns :

A new NcClusterAbundance. [transfer full]

nc_cluster_abundance_d2NdzdlnM_Mobs ()

gdouble             nc_cluster_abundance_d2NdzdlnM_Mobs (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble z);

This function computes /f$ \int_{\ln M^{obs} - 7\sigma_{\ln M}}^{\ln M^{obs} + 7\sigma_{\ln M}} d\ln M \, \frac{d^2N}{dzdlnM} * P(\ln M^{obs}|\ln M) /f$. The integral limits were determined requiring a precision to five decimal places.

cad :

a NcClusterAbundance.

model :

a NcHICosmo.

lnMobs :

the logarithm base e of the "observable" mass (gdouble).

z :

redshift (gdouble).

Returns :

a gdouble which corresponds to /f$ \int_{\ln M^{obs} - 7\sigma_{\ln M}}^{\ln M^{obs} + 7\sigma_{\ln M}} d\ln M \, \frac{d^2N}{dzdlnM} * P(\ln M^{obs}|\ln M) /f$.

nc_cluster_abundance_d2NdzdlnM_Mobs_local_selection ()

gdouble             nc_cluster_abundance_d2NdzdlnM_Mobs_local_selection
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble z);

This function computes \int_0^\infty d\ln M \frac{d^2N(\ln M, z)}{dzd\ln M} * c(\lnM, z) * P(\ln M^{obs}|\ln M, z). The sample/data has specific redshift and mass ranges, therefore we have to test if zp and lnMobs are inside these ranges. We have also to do the integral in pieces with limits corresponding to the lnM range of the histogram setting the completeness /f$ c(\ln M, z) /f$ and /f$ \sigma_{\lnM} /f$ spcecific of each bin. We determined the integral limits /f$ [\ln M^{obs} - 7\sigma_{\ln M}, \ln M^{obs} + 7\sigma_{\ln M}] /f$ to obtain a precision to five decimal places.

It is worth emphasizing that purity is not included, since it is a function of observable mass and photometric redshift. See nc_cluster_abundance_d2NdzdlnM_purity_val, nc_cluster_abundance_dNdz_purity_val and nc_cluster_abundance_N_purity_val.

cad :

a NcClusterAbundance.

model :

a NcHICosmo.

lnMobs :

logarithm (base e) of the observable mass

z :

redshift

Returns :

a gdouble which represents /f$ \frac{d^2N(\ln M^{obs}, z)}{dzd\lnM} /f$ taking into account the completeness.

nc_cluster_abundance_d2NdzdlnM_photoz ()

gdouble             nc_cluster_abundance_d2NdzdlnM_photoz
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnM,
                                                         gdouble z_phot);

This function computes /f$ \int_{z_{phot} - 10\sigma_{phot}}^{z_{phot} + 10\sigma_{phot}} dz \, \frac{d^2N}{dzdlnM} * P(z^{photo}|z) /f$. The integral limits were determined requiring a precision to five decimal places.

cad :

a NcClusterAbundance.

model :

a NcHICosmo.

lnM :

the logarithm base e of the mass (gdouble).

z_phot :

the photometric redshift (gdouble).

Returns :

a gdouble which corresponds to /f$ \int_{z_{phot} - 10\sigma_{phot}}^{z_{phot} + 10\sigma_{phot}} dz \, \frac{d^2N}{dzdlnM} * P(z^{photo}|z) /f$.

nc_cluster_abundance_d2NdzdlnM_photoz_Mobs ()

gdouble             nc_cluster_abundance_d2NdzdlnM_photoz_Mobs
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble zp);

This function computes /f$ \int_0^\infty dz \int_0^\infty d\ln M \frac{d^2N(\ln M, z)}{dzd\ln M} * P(z^{phot}|z) * P(\ln M^{obs}|\ln M, z) /f$. We studied the convergence of this integral to optimize this function. We verified that it converges to 5 decimal places at the redshift interval /f$ [z^{phot} - 10\sigma^{phot}, z^{phot} + 10\sigma^{phot}] /f$ and the mass interval /f$ [\ln M^{obs} - 7\sigma_{\ln M}, \ln M^{obs} + 7\sigma_{\ln M}] /f$.

cad :

a NcClusterAbundance.

model :

a NcHICosmo.

lnMobs :

logarithm (base e) of the observable mass.

zp :

photometric redshift.

Returns :

a gdouble which represents /f$ \frac{d^2N(\ln M^{obs}, z^{phot})}{dzd\lnM} /f$.

nc_cluster_abundance_d2NdzdlnM_photoz_Mobs_local_selection ()

gdouble             nc_cluster_abundance_d2NdzdlnM_photoz_Mobs_local_selection
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble zp);

This function computes /f$ \int_0^\infty dz \int_0^\infty d\ln M \frac{d^2N(\ln M, z)}{dzd\ln M} * c(\lnM, z) * P(\ln M^{obs}|\ln M, z) * P(z^{photo}, z) /f$. Analogously to our previous work, we determine the integration limits considering the maximum value of /f$ \sigma_{\ln M} /f$. Therefore, the mass interval is /f$ [\ln M^{obs} - 7\sigma_{\ln M}, \ln M^{obs} + 7\sigma_{\ln M}] /f$ and the redshift interval is /f$ [z^{photo} - 10\sigma, z^{photo} + 10\sigma] /f$. Purity is included later because it only depends on observables mass and redshift.

cad :

a NcClusterAbundance.

model :

a NcHICosmo.

lnMobs :

the logarithm (base e) of the observable mass (gdouble).

zp :

the photometric redshift (gdouble).

Returns :

a gdouble which represents /f$ \frac{d^2N(\ln M^{obs}, z^{photo})}{dzd\lnM} /f$ taking into account the completeness.

nc_cluster_abundance_d2NdzdlnM_purity_val ()

gdouble             nc_cluster_abundance_d2NdzdlnM_purity_val
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs,
                                                         gdouble z);

nc_cluster_abundance_d2NdzdlnM_val()

#define nc_cluster_abundance_d2NdzdlnM_val(cad,cp,lnM,z) (cad)->d2NdzdlnM_val(cad,cp,lnM,z)

nc_cluster_abundance_dNdlnM_val()

#define nc_cluster_abundance_dNdlnM_val(cad,cp,lnM,zl,zu) (cad)->dNdlnM_val(cad,cp,lnM,zl,zu)

nc_cluster_abundance_dNdz_Mobs_local_selection ()

gdouble             nc_cluster_abundance_dNdz_Mobs_local_selection
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs_i,
                                                         gdouble lnMobs_f,
                                                         gdouble z);

This function computes \int_lnMi^lnMf d lnMobs \frac{d^2N(\ln M, z)}{dzd\ln Mobs} / p(lnMobs, z). The sample/data has specific redshift and mass ranges, therefore we have to test if z and lnMobs are inside these ranges. We have also to do the integral in pieces with limits corresponding to the lnM range of the histogram setting the specific purity /f$ p(\ln Mobs, z) /f$ of each bin.

cad :

a NcClusterAbundance.

model :

a NcHICosmo.

lnMobs_i :

logarithm base e of the observable minimum mass threshold.

lnMobs_f :

logarithm base e of the observable maximum mass threshold.

z :

redshift.

Returns :

a gdouble which represents /f$ \frac{dN(\ln M^{obs}, z)}{dz} /f$ taking into account the completeness, purity and P(lnMobs | lnM) with varying standard deviation.

nc_cluster_abundance_dNdz_photoz_Mobs_local_selection ()

gdouble             nc_cluster_abundance_dNdz_photoz_Mobs_local_selection
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMobs_i,
                                                         gdouble lnMobs_f,
                                                         gdouble z);

This function computes \int_lnMi^lnMf d lnMobs \frac{d^2N(\ln M, z)}{dzd\ln Mobs} / p(lnMobs, z). The sample/data has specific redshift and mass ranges, therefore we have to test if z and lnMobs are inside these ranges. We have also to do the integral in pieces with limits corresponding to the lnM range of the histogram setting the specific purity /f$ p(\ln Mobs, z) /f$ of each bin.

cad :

a NcClusterAbundance.

model :

a NcHICosmo.

lnMobs_i :

logarithm base e of the observable minimum mass threshold.

lnMobs_f :

logarithm base e of the observable maximum mass threshold.

z :

redshift.

Returns :

a gdouble which represents /f$ \frac{dN(\ln M^{obs}, z)}{dz} /f$ taking into account the completeness, purity and P(lnMobs | lnM) with varying standard deviation.

nc_cluster_abundance_dNdz_purity_val ()

gdouble             nc_cluster_abundance_dNdz_purity_val
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble lnMl,
                                                         gdouble lnMu,
                                                         gdouble z);

nc_cluster_abundance_dNdz_val()

#define nc_cluster_abundance_dNdz_val(cad,cp,lnMl,lnMu,z) (cad)->dNdz_val(cad,cp,lnMl,lnMu,z)

nc_cluster_abundance_free ()

void                nc_cluster_abundance_free           (NcClusterAbundance *cad);

Atomically decrements the reference count of cad by one. If the reference count drops to 0, all memory allocated by cad is released.

cad :

a NcClusterAbundance.

nc_cluster_abundance_get_lnM_sigma0 ()

gdouble             nc_cluster_abundance_get_lnM_sigma0 (NcClusterAbundance *cad);

Gets the value of the "lnM-sigma0" property.

cad :

a NcClusterAbundance.

Returns :

the value of "lnM-sigma0" property.

nc_cluster_abundance_get_lnMf ()

gdouble             nc_cluster_abundance_get_lnMf       (NcClusterAbundance *cad);

Gets the value of the "maximum-mass" property.

cad :

a NcClusterAbundance.

Returns :

the value of "maximum-mass" property.

nc_cluster_abundance_get_lnMi ()

gdouble             nc_cluster_abundance_get_lnMi       (NcClusterAbundance *cad);

Gets the value of the "minimum-mass" property.

cad :

a NcClusterAbundance.

Returns :

the value of "minimum-mass" property.

nc_cluster_abundance_get_photoz ()

NcClusterPhotoz *   nc_cluster_abundance_get_photoz     (NcClusterAbundance *cad);

Gets the value of the "photoz" property.

cad :

a NcClusterAbundance.

Returns :

the value of "photoz" property. [transfer none]

nc_cluster_abundance_get_zf ()

gdouble             nc_cluster_abundance_get_zf         (NcClusterAbundance *cad);

Gets the value of the "maximum-redshift" property.

cad :

a NcClusterAbundance.

Returns :

the value of "maximum-redshift" property.

nc_cluster_abundance_get_zi ()

gdouble             nc_cluster_abundance_get_zi         (NcClusterAbundance *cad);

Gets the value of the "minimum-redshift" property.

cad :

a NcClusterAbundance.

Returns :

the value of "minimum-redshift" property.

nc_cluster_abundance_new ()

NcClusterAbundance * nc_cluster_abundance_new           (NcClusterAbundanceOpt opt,
                                                         NcMassFunction *mfp,
                                                         NcHaloBiasFunc *mbiasf,
                                                         gdouble zi,
                                                         gdouble zf,
                                                         gdouble lnMi,
                                                         gdouble lnMf,
                                                         NcClusterPhotoz *photoz,
                                                         gdouble lnM_sigma0);

This function allocates memory for a new NcClusterAbundance object and sets its properties to the values from the input arguments.

Returns: A new NcClusterAbundance.

opt :

a NcClusterAbundanceOpt.

mfp :

a NcMassFunction.

mbiasf :

a NcHaloBiasFunc. [allow-none]

zi :

minimum redshift.

zf :

maximum redshift.

lnMi :

minimum logarithm base e of mass.

lnMf :

maximum logarithm base e of mass.

photoz :

a NcClusterPhotoz.

lnM_sigma0 :

FIXME

nc_cluster_abundance_prepare ()

void                nc_cluster_abundance_prepare        (NcClusterAbundance *cad,
                                                         NcHICosmo *model);

nc_cluster_abundance_prepare_inv_dNdlnM_z ()

void                nc_cluster_abundance_prepare_inv_dNdlnM_z
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble z);

This function prepares a spline where the x array corresponds to the value of \int_lnM0 ^lnM1 d2N/dzdlnM dM/ \int_lnMi^lnMf dN/dz dM given a redshift z and the y array contains the values of logarithms base e of the mass. It is used to generate a sample of lnM values.

cad :

a NcClusterAbundance.

model :

a NcHICosmo.

z :

redshift.

nc_cluster_abundance_prepare_inv_dNdlnM_z_no_obs ()

void                nc_cluster_abundance_prepare_inv_dNdlnM_z_no_obs
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model,
                                                         gdouble z);

This function prepares a spline where the x array corresponds to the value of \int_lnM0 ^lnM1 d2N/dzdlnM dM/ \int_lnMi^lnMf dN/dz dM given a redshift z and the y array contains the values of logarithms base e of the mass. It is used to generate a sample of lnM values.

cad :

a NcClusterAbundance.

model :

a NcHICosmo.

z :

redshift.

nc_cluster_abundance_prepare_inv_dNdz ()

void                nc_cluster_abundance_prepare_inv_dNdz
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model);

This function prepares a bidimensional spline...

cad :

a NcClusterAbundance.

model :

a NcHICosmo.

nc_cluster_abundance_prepare_inv_dNdz_no_obs ()

void                nc_cluster_abundance_prepare_inv_dNdz_no_obs
                                                        (NcClusterAbundance *cad,
                                                         NcHICosmo *model);

This function prepares a bidimensional spline...

cad :

a NcClusterAbundance.

model :

a NcHICosmo.

nc_cluster_abundance_realizations_read_from_file ()

GPtrArray *         nc_cluster_abundance_realizations_read_from_file
                                                        (gchar *file_realization,
                                                         gint n_realizations);

GPtrArray *realizations is an array of array with the z values of all realizations. To complete...

file_realization :

is the file's name which contains the redshift values of the clusters obtained with random Poisson generator. [array length=n_realizations]

n_realizations :

is the number of realizations generated.

Returns :

FIXME. [transfer full]

nc_cluster_abundance_realizations_save_to_file ()

void                nc_cluster_abundance_realizations_save_to_file
                                                        (GPtrArray *realizations,
                                                         gchar *filename);

nc_cluster_abundance_set_lnM_sigma0 ()

void                nc_cluster_abundance_set_lnM_sigma0 (NcClusterAbundance *cad,
                                                         const gdouble lnM_sigma0);

Sets the value lnM_sigma0 to the "lnM-sigma0" property.

cad :

a NcClusterAbundance.

lnM_sigma0 :

value of "lnM-sigma0".

nc_cluster_abundance_set_lnMf ()

void                nc_cluster_abundance_set_lnMf       (NcClusterAbundance *cad,
                                                         const gdouble lnMf);

Sets the value lnMf to the "maximum-mass" property.

cad :

a NcClusterAbundance.

lnMf :

value of "maximum-mass".

nc_cluster_abundance_set_lnMi ()

void                nc_cluster_abundance_set_lnMi       (NcClusterAbundance *cad,
                                                         const gdouble lnMi);

Sets the value lnMi to the "minimum-mass" property.

cad :

a NcClusterAbundance.

lnMi :

value of "minimum-mass".

nc_cluster_abundance_set_options ()

void                nc_cluster_abundance_set_options    (NcClusterAbundance *cad,
                                                         NcClusterAbundanceOpt opt);

nc_cluster_abundance_set_photoz ()

void                nc_cluster_abundance_set_photoz     (NcClusterAbundance *cad,
                                                         NcClusterPhotoz *photoz);

Sets the value photoz to the "photoz" property.

cad :

a NcClusterAbundance.

photoz :

value of "photoz".

nc_cluster_abundance_set_zf ()

void                nc_cluster_abundance_set_zf         (NcClusterAbundance *cad,
                                                         const gdouble zf);

Sets the value zf to the "maximum-redshift" property.

cad :

a NcClusterAbundance.

zf :

value of "maximum-redshift".

nc_cluster_abundance_set_zi ()

void                nc_cluster_abundance_set_zi         (NcClusterAbundance *cad,
                                                         const gdouble zi);

Sets the value zi to the "minimum-redshift" property.

cad :

a NcClusterAbundance.

zi :

value of "minimum-redshift".

Property Details

The "lnM-sigma0" property

  "lnM-sigma0"               gdouble               : Read / Write

Satandard deviation factor of a log-normal mass-observable distribution: sd = sigma0.

Default value: 0


The "mass-function" property

  "mass-function"            NcMassFunction*       : Read / Write / Construct Only

Mass Function.


The "maximum-mass" property

  "maximum-mass"             gdouble               : Read / Write

Maximum mass value is restricted to [11.0 * ln10, 16.0 * ln10] in units of h^{-1} M_sun.

Allowed values: [25.3284,36.8414]

Default value: 25.3284


The "maximum-redshift" property

  "maximum-redshift"         gdouble               : Read / Write

Maximum redshift value is restricted to [_NC_CLUSTER_ABUNDANCE_MIN_Z, G_MAXDOUBLE].

Allowed values: >= 1e-06

Default value: 1e-06


The "mean-bias" property

  "mean-bias"                NcHaloBiasFunc*       : Read / Write / Construct Only

Mean Halo Bias Function.


The "minimum-mass" property

  "minimum-mass"             gdouble               : Read / Write

Minimum mass value is restricted to [11.0 * ln10, 16.0 * ln10] in units of h^{-1} M_sun.

Allowed values: [25.3284,36.8414]

Default value: 25.3284


The "minimum-redshift" property

  "minimum-redshift"         gdouble               : Read / Write

Minimum redshift value is restricted to [_NC_CLUSTER_ABUNDANCE_MIN_Z, G_MAXDOUBLE].

Allowed values: >= 1e-06

Default value: 1e-06


The "options" property

  "options"                  NcClusterAbundanceOpt  : Read / Write / Construct Only

Options.


The "photoz" property

  "photoz"                   NcClusterPhotoz*      : Read / Write / Construct Only

Photoz.