Skip to content

Latest commit

 

History

History
48 lines (37 loc) · 2.01 KB

README.md

File metadata and controls

48 lines (37 loc) · 2.01 KB

UNIONS Cosmological Inference Pipeline

by Lisa Goh and Sacha Guerrini, CEA Paris-Saclay

This folder contains the files neccessary to run the cosmological inference pipeline on the UNIONS galaxy catalogues.

Requirements

To run the pipeline, one would need to have installed CosmoSIS and CosmoCov. To PSF leakage parameters, the fork of cosmosis-standard-library of Sacha Guerrini has to be used.

To Run

The inference pipeline is now orchestrated through Python. Run the main Snakemake workflow from the parent directory:

snakemake -j<jobs> inference_fiducial

This will automatically execute all steps:

  1. Calculate 2PCF ($\xi_{pm}$) via cosmo_val.py
  2. Compute covariance matrices using CosmoCov
  3. Prepare CosmoSIS data (FITS) via cosmosis_fitting.py
  4. Run CosmoSIS inference

For standalone FITS data preparation (real-space inputs plus optional pseudo-$C_\ell$ data), you can also use the Python script directly:

python scripts/cosmosis_fitting.py \
  --cosmosis-root "catalog_version_config" \
  --data-dir "/path/to/output/chains" \
  --nz-file "/path/to/nz_file.txt" \
  --out-file "/path/to/output.fits" \
  --xi "/path/to/xi_plus.fits" "/path/to/xi_minus.fits" \
  --cov-xi "/path/to/covariance.txt" \
  --use-rho-tau \
  --rho-stats "/path/to/rho_stats.fits" \
  --tau-stats "/path/to/tau_stats.fits" \
  --cov-tau "/path/to/cov_tau.npy" \
  --cl-file "/path/to/pseudo_cl.fits" \
  --cov-cl "/path/to/pseudo_cl_cov.fits"

You can view all available options with:

python scripts/cosmosis_fitting.py --help

Ensure the pseudo-$C_\ell$ spectra (pseudo_cl_*.fits) and their covariance (pseudo_cl_cov_*.fits) produced by cosmo_val.py exist for the requested catalog version (or mock seed) before running the standalone command.

This is the pipeline used to derive cosmological constraints with cosmic shear data from the UNIONS v1.4 catalogue.