measure_cat_3x2pt_pcls module
Measure the 3x2pt Pseudo-Cl power spectra from the mock catalogues generated by the catalogue_sim package. This code reads in a file describing the n(z) bin boundaries to measure the 3x2pt tomographic power spectra. At present (13/06/24) this file must be produced by the create_nz_boundaries.py script within the pcl_measurement package. It is in the roadmap to have this file read in e.g. from an external user-defined file.
The power spectra are measured to a resolution and ell range as specified by the set_variables_3x2pt_measurement.ini config file. The outputs are the ‘raw’ 3x2pt Pseudo-Cl power spectra and noise power spectra - before these may have scale cuts or bandpower binning applied later in the pipeline
- pcl_measurement.measure_cat_3x2pt_pcls.execute_pcl_measurement()
Execute the script to measure the 3x2pt Pseudo-Cl power spectrum. First create a config dictionary of parameters, then create 3x2pt maps, then measure the power spectra.
Returns
Saves 3x2pt power spectra (and noise contributions) to disk
- pcl_measurement.measure_cat_3x2pt_pcls.maps_from_cats(config_dict, iter_no)
From the simulated catalogues, first create the 3x2pt maps, i.e. convergence, y1, y2, galaxy overdensity and store in a dictionary to then extract and measure the power spectra
Parameters
- config_dict(dict)
Dictionary containing the measurement parameters (e.g. redshift range, Healpix resolution) to extract information from the catalogues and create maps. Created by measure_pcls_config()
- iter_no(int)
Index denoting realisation number of the simulation
Returns
Dictionary (cut_maps_dic) containing the tomographic maps corresponding to each field of the 3x2pt data vector
- pcl_measurement.measure_cat_3x2pt_pcls.measure_00_pcls(lmin_out, lmax_out, cut_maps_dic, spectra_type, bin_i, bin_j, measured_pcl_save_dir)
Function to measure the Pseudo Cl power spectrum between two spin-0 fields
Parameters
- lmin_out(float)
Minimum ell that measured power spectra are stored and saved to disk
- lmax_out(float)
Maximum ell that measured power spectra are stored and saved to disk
- cut_maps_dic(dict)
Dictionary of tomographic 3x2pt maps to measure power spectra from - generated by maps_from_cats()
- spectra_type(str)
Spectra type that is measured. For 00, must be either ‘gal-gal’ for the galaxy overdensity, or ‘TT’ for the convergence component of the shear
- bin_i(float)
Tomographic bin id number of the first field
- bin_j(float)
Tomographic bin id number of the second field
- measured_pcl_save_dir(str)
Location on disk to store the measured Pseudo-Cls
Returns
Saves a file tree of the measured 00 3x2pt Pseudo-Cl power spectra
- pcl_measurement.measure_cat_3x2pt_pcls.measure_02_pcls(lmin_out, lmax_out, cut_maps_dic, spectra_type, bin_i, bin_j, measured_pcl_save_dir)
Function to measure the Pseudo Cl power spectrum between a spin-0 and spin-2 field
Parameters
- lmin_out(float)
Minimum ell that measured power spectra are stored and saved to disk
- lmax_out(float)
Maximum ell that measured power spectra are stored and saved to disk
- cut_maps_dic(dict)
Dictionary of tomographic 3x2pt maps to measure power spectra from - generated by maps_from_cats()
- spectra_type(str)
Spectra type that is measured. For 02, must be either ‘TE’ or ‘TB’ for the cosmic shear, or ‘gal_E’ or ‘gal_B’ for the cross correlation between galaxy overdensity and cosmic shear
- bin_i(float)
Tomographic bin id number of the first field
- bin_j(float)
Tomographic bin id number of the second field
- measured_pcl_save_dir(str)
Location on disk to store the measured Pseudo-Cls
Returns
Saves a file tree of the measured 02 3x2pt Pseudo-Cl power spectra
- pcl_measurement.measure_cat_3x2pt_pcls.measure_22_pcls(lmin_out, lmax_out, cut_maps_dic, spectra_type, bin_i, bin_j, measured_pcl_save_dir)
Function to measure the Pseudo Cl power spectrum between two spin-2 fields
Parameters
- lmin_out(float)
Minimum ell that measured power spectra are stored and saved to disk
- lmax_out(float)
Maximum ell that measured power spectra are stored and saved to disk
- cut_maps_dic(dict)
Dictionary of tomographic 3x2pt maps to measure power spectra from - generated by maps_from_cats()
- spectra_type(str)
Spectra type that is measured. For 22, must be either ‘EE’, ‘EB’, ‘BE’ or ‘BB’ for the cosmic shear
- bin_i(float)
Tomographic bin id number of the first field
- bin_j(float)
Tomographic bin id number of the second field
- measured_pcl_save_dir(str)
Location on disk to store the measured Pseudo-Cls
Returns
Saves a file tree of the measured 22 3x2pt Pseudo-Cl power spectra
- pcl_measurement.measure_cat_3x2pt_pcls.measure_pcls_config(pipeline_variables_path)
Set up a config dictionary of parameters that control how the 3x2pt Pseudo-Cls are measured - will be read in from set_variables_3x2pt_measurement.ini config file
Parameters
- pipeline_variables_path(str)
Path to location of pipeline variables file (‘set_variables_3x2pt_measurement.ini’)
Returns
Dictionary of Pseudo-Cl measurement parameters