UNIONS shear analyses

Usage	Development	Release

---

Author: Martin Kilbinger
Email: martin.kilbinger@cea.fr
Year: 2022

---

(Copied from the Pyralid template. See pyraliddemo for a demo package created with the Pyralid template.)

Contents

1. Requirements
2. AGN-GGL
3. Management
4. Deployment
5. CosmoStat

Requirements

The python packages treecorr, àstropy, statsmodel and cs_util need to be installed.

AGN galaxy-galaxy lensing

Input data

First, you need the AGN catalogue, optically selected in the SDSS survey. This data file is part of this repository, data/agn_ggl/SDSS_SMBH_202206.txt. Second, the weak-lensing catalogue used here is the (internally) release v1.0, produced by either the ShapePipe or Lensfit weak-lensing pipeline. E.g. use the file unions_shapepipe_2022_v1.0.fits.

Split input catalogues in black-hole mass bins

Run

scripts/split_sample.py -v -n 2

to create equi-populated mass sub-samples.

Compute GGL

Ignoring weights, use

scripts/compute_ng_binned_samples.py --input_path_fg data_mass_sub/SDSS_SMBH_202206_0_n_split_2.fits --input_path_bg unions_shapepipe_2022_v1.0.fits --key_ra_fg ra --key_dec_fg dec -v --out_path data_mass_sub/ggl_agn_0.txt

to compute the weak-lensing tangential shear of UNIONS galaxies around the lower-mass bin (#0). Repeat for the high-mass bin (#1).