Calibration and Splitting

These stages deal with calibrating shear and splitting up catalogs into sub-catalogs.

class txpipe.lens_selector.TXLensCatalogSplitter(args, comm=None, aliases=None)[source]

Split a lens catalog file into a new file with separate bins

Splitting up like this helps reduce memory usage in TreeCorr later

Inputs:

  • lens_tomography_catalog_unweighted: TomographyCatalog

  • photometry_catalog: PhotometryCatalog

  • fiducial_cosmology: FiducialCosmology

  • lens_photoz_pdfs: HDFFile

Outputs:

  • binned_lens_catalog_unweighted: HDFFile

Parallel: Yes - MPI

Configuration
  • initial_size: (int) Default=100000.
  • chunk_rows: (int) Default=100000.
  • extra_cols: (list) Default=[''].
  • redshift_column: (str) Default=zmean.
class txpipe.lens_selector.TXTruthLensCatalogSplitter(args, comm=None, aliases=None)[source]

Split a lens catalog file into a new file with separate bins with true redshifts.

The redshifts are used to calculate the comoving distances.

Inputs:

  • lens_tomography_catalog_unweighted: TomographyCatalog

  • photometry_catalog: PhotometryCatalog

  • fiducial_cosmology: FiducialCosmology

Outputs:

  • binned_lens_catalog_unweighted: HDFFile

Parallel: Yes - MPI

Configuration
  • initial_size: (int) Default=100000.
  • chunk_rows: (int) Default=100000.
  • extra_cols: (list) Default=[''].
  • redshift_column: (str) Default=redshift_true.
class txpipe.lens_selector.TXExternalLensCatalogSplitter(args, comm=None, aliases=None)[source]

Split an external lens catalog into bins

Implemented as a subclass of TXLensCatalogSplitter, and changes only file names.

Inputs:

  • lens_tomography_catalog_unweighted: TomographyCatalog

  • lens_catalog: HDFFile

  • fiducial_cosmology: FiducialCosmology

Outputs:

  • binned_lens_catalog_unweighted: HDFFile

Parallel: Yes - MPI

Configuration
  • initial_size: (int) Default=100000.
  • chunk_rows: (int) Default=100000.
  • extra_cols: (list) Default=[''].
  • redshift_column: (str) Default=redshift.
class txpipe.lens_selector.TXTruthLensCatalogSplitterWeighted(args, comm=None, aliases=None)[source]

Split a lens catalog file into a new file with separate bins with true redshifts.

The redshifts are used to calculate the comoving distances.

Inputs:

  • lens_tomography_catalog: TomographyCatalog

  • photometry_catalog: PhotometryCatalog

  • fiducial_cosmology: FiducialCosmology

Outputs:

  • binned_lens_catalog: HDFFile

Parallel: Yes - MPI

Configuration
  • initial_size: (int) Default=100000.
  • chunk_rows: (int) Default=100000.
  • extra_cols: (list) Default=[''].
  • redshift_column: (str) Default=redshift_true.
class txpipe.twopoint_null_tests.TXStarCatalogSplitter(args, comm=None, aliases=None)[source]

Split a star catalog into bright and dim stars

Inputs:

  • star_catalog: HDFFile

Outputs:

  • binned_star_catalog: HDFFile

Parallel: No - Serial

Configuration
  • chunk_rows: (int) Default=100000.
  • initial_size: (int) Default=100000.
class txpipe.calibrate.TXShearCalibration(args, comm=None, aliases=None)[source]

Split the shear catalog into calibrated bins

This class runs after source selection has been done, because the final calibration factor can only be estimated once we have read the entire catalog and chosen tomographic bins (since it is an ensemble average of cal factors).

Once that stage has run and computed both the tomographic bin for each sample and the calibration factors, this stage takes the full catalog and splits it into one HDF5 group per bin.

We are not (yet) saving per-patch catalogs for TreeCorr here. We might want to do that later.

Inputs:

  • shear_catalog: ShearCatalog

  • shear_tomography_catalog: TomographyCatalog

Outputs:

  • binned_shear_catalog: ShearCatalog

Parallel: Yes - MPI

Configuration
  • use_true_shear: (bool) Default=False.
  • chunk_rows: (int) Default=100000.
  • subtract_mean_shear: (bool) Default=True.
  • extra_cols: (list) Default=[''].
  • shear_catalog_type: (str) Default=.
  • shear_prefix: (str) Default=.