Library of functions to perform halo abundance matching¶
This class handles a light cone and assigns galaxies to the halos
Library of function to create halo catalogs matched to a density.
-
class
HaloSelection.MultiDarkMock(hdu, area, mockOutput_dir, mockName, zmin, zmax, nGal_Deg2)[source]¶ Parameters: - hdu – hdu of the lightcone
- area – area in deg2
- mockOutput_dir – directory where to output mocks
- mockName – name of the file where to save the mock
- zmin – minimum redshift array defining the bins
- zmax – maximum redshift array defining the bins
- nGal_Deg2 – number of galaxies per square degrees
-
compare_clustering_data_mock(w_data, xi_data, theta_min_chi2=-2.3, theta_max_chi2=-1.5, w_bins=15, s_min_chi2=0.5, s_max_chi2=1.2, s_bins=10)[source]¶ Compares the clustering of the mock catalog and the clustering of the data. :param w_data: angular clustering from the data [x, y, yErr]. :param xi_data: monopole clustering from the data [x, y, yErr].
-
get_distrib_QTY(colN, z1, z2)[source]¶ Computes the cumulative histogram of a column for halos in the range z1, z2. :param colN: name of the column you want to take the histogram. :param z1: minimum redshift :param z2: maximum redshift
-
get_distrib_QTY_cen(colN, z1, z2)[source]¶ Computes the cumulative histogram of a column for central halos in the range z1, z2. :param colN: name of the column you want to take the histogram. :param z1: minimum redshift :param z2: maximum redshift
-
get_distrib_QTY_sat(colN, z1, z2)[source]¶ Computes the cumulative histogram of a column for satellite halos in the range z1, z2. :param colN: name of the column you want to take the histogram. :param z1: minimum redshift :param z2: maximum redshift
-
initialize()[source]¶ Initializes the procedure by putting into memroy arrays of central and satellites.
-
make_GaussianFsat_catalog(colN, means, scatters, fsats)[source]¶ Creates lists of ids of centrals and satellite galaxies corresponding to the density given in the n(z). :param colN: name of the column you construct the catalog with :param means: means of the Gaussians, array the same length of the redshift bin :param scatters: scatters of the Gaussians, array the same length of the redshift bin :param fsats: fractions of satellite, array the same length of the redshift bin
-
make_Gaussian_catalog(colN, means, scatters)[source]¶ Creates lists of ids of centrals and satellite galaxies corresponding to the density given in the n(z). :param colN: name of the column you construct the catalog with :param means: means of the Gaussians, array the same length of the redshift bin :param scatters: scatters of the Gaussians, array the same length of the redshift bin
-
make_LogNorm_catalog(colN, means, scatters)[source]¶ Creates lists of ids of centrals and satellite galaxies corresponding to the density given in the n(z). :param colN: name of the column you construct the catalog with :param means: means of the Gaussians, array the same length of the redshift bin :param scatters: scatters of the Gaussians, array the same length of the redshift bin
-
make_shamIncomplete_catalog(colN, incompletenessFactor)[source]¶ Creates lists of ids of halos corresponding to the density given in the n(z). For every bin of redshift, it gets the distribution o fthe column of interest and matches to the density of galaxies in the NZ given. Then provides a column of ids extracted from teh lightcone. :param colN: name of the column you wish to work on for the sham.
-
make_shamMAX_catalog(colN, maxQTY)[source]¶ Creates lists of ids of centrals and satellite galaxies corresponding to the density given in the n(z). For every bin of redshift, it gets the distribution o fthe column of interest and matches to the density of galaxies in the NZ given. Then provides a column of ids extracted from teh lightcone. :param colN: name of the column you wish to work on for the sham.
-
make_sham_catalog(colN='mvir')[source]¶ Creates lists of ids of halos corresponding to the density given in the n(z). For every bin of redshift, it gets the distribution o fthe column of interest and matches to the density of galaxies in the NZ given. Then provides a column of ids extracted from teh lightcone. :param colN: name of the column you wish to work on for the sham.
-
select_Gaussian(meanQTY, scatterQTY, nGal_perbin, IDhz, QTY)[source]¶ Creates lists of ids of centrals and satellite galaxies corresponding to the density given in the n(z) and to a gaussian distribution . For every bin of redshift, it gets the distribution o fthe column of interest and matches to the density of galaxies in the NZ given. Then provides a column of ids extracted from teh lightcone. :param colN: name of the column you wish to work on for the sham. :param meanQTY: mean of the distribution :param scatterQTY: scatter of the distribution :param nGal_perbin: total number of galaxies in this bins to mock :param IDhz: IDs of the halos in this bin :param QTY: array of the column to do the match on, mass, velocity, ...
-
select_GaussianFsat(meanQTY, scatterQTY, fsat, nGal_perbin, IDhz_c, QTY_c, IDhz_s, QTY_s)[source]¶ Extracts the ids of halos to create a mock with a Gaussian distribution. :param colN: name of the column you wish to work on for the sham. :param meanQTY: mean of the distribution :param scatterQTY: scatter of the distribution :param fsat: fraction of satellite in this bin :param nGal_perbin: total number of galaxies in this bins to mock :param IDhz_c: IDs of the central halos in this bin :param QTY_c: column to do the match on, mass, velocity, ... for the central halos :param IDhz_s: IDs of the satellite halos in this bin :param QTY_s: column to do the match on, mass, velocity, ... for the satellite halos
-
select_LogNorm(meanQTY, scatterQTY, nGal_perbin, IDhz, QTY, nn, bb)[source]¶ Creates lists of ids of centrals and satellite galaxies corresponding to the density given in the n(z) and to a gaussian distribution . For every bin of redshift, it gets the distribution o fthe column of interest and matches to the density of galaxies in the NZ given. Then provides a column of ids extracted from teh lightcone. :param colN: name of the column you wish to work on for the sham. :param meanQTY: mean of the distribution :param scatterQTY: scatter of the distribution :param nGal_perbin: total number of galaxies in this bins to mock :param IDhz: IDs of the halos in this bin :param QTY: array of the column to do the match on, mass, velocity, ...
-
select_sham(nGal_perbin, IDhz, QTY, nn, bb)[source]¶ Returns the ids corresponding to a given density. :param nGal_perbin: number of galaxies to be selected :param IDhz: parent ID distribution :param QTY: quantity to select halos on :param nn: cumulative distribution o QTY :param bb: bins of the cumulative distribution
-
select_shamIncomplete(incompFactor, nGal_perbin, IDhz, QTY, nn, bb)[source]¶ Returns the ids corresponding to a given density and an incompleteness factor. :param nGal_perbin: number of galaxies to be selected :param IDhz: parent ID distribution :param QTY: quantity to select halos on :param nn: cumulative distribution o QTY :param bb: bins of the cumulative distribution :param incompFactor: incompleteness factor compared to the max of QTY : max(QTY)/incompFactor will be set as the max of the extracted distribution.
-
select_shamMAX(QTY_max, nGal_perbin, IDhz, QTY, nn, bb)[source]¶ Returns the ids corresponding to a given density with a maximum in the QTY. :param nGal_perbin: number of galaxies to be selected :param IDhz: parent ID distribution :param QTY: quantity to select halos on :param nn: cumulative distribution o QTY :param bb: bins of the cumulative distribution :param QTY_max: the max of QTY is set to QTY_max.
-
writeClusteringParamFile(type, decade='')[source]¶ Writes the clustering commands that command the CUTE code, see Alonso et al. 2012 https://arxiv.org/abs/1210.1833 :param type: monopole or angular or ... :param decade: string suffix that is appended if you study different scales (decades) _d1, _d2, _d3 are used for the angular clustering.