Introduction

From a OU-SHE/OU-PHZ catalogue lensMCCatalog2DMassWL (or momentsMLCatalog2DMassWL, or regaussCatalog2DMassWL, or kSBCatalog2DMassWL), 2D-MASS-WL will provide the following output products:

• E/B modes for small patches, noisy and/or denoised, with their associated SNR maps,

• E/B modes for a single cluster, noisy and/or denoised, with associated SNR maps, also collected in a single tar for all sufficiently massive clusters in a catalogue

• MultiPatch E/B modes for the “full” sky, with their associated SNR maps, computed via decomposition of the sky into multiple patches for which convergence maps are computed,

• Spherical E/B modes for the “full” sky, with their associated SNR maps; an alternative is based on Spherical E/B convergence maps <2DMassIntroConvergenceSpherical>,

• Peak catalogues for statistical studies.

E/B modes and SNR maps for tests (small patches)

To test the production of convergence maps for small fields (for object studies such as cluster), an input FITS catalogue (e.g lensMCCatalog2DMassWL) collecting information from OU-SHE/OU-PHZ is needed, as well as a set of parameters twoDMassParams2DConvergencePatch (for convergence maps) specifying the location of the patches and running parameters. The noisy/denoised convergence maps and SNR products are then obtained after projection of the input region, located on the sphere, and collected in twoDMassConvergencePatch.

The flowchart is the following:

E/B modes and SNR maps for cluster studies

A catalogue clustercatalog2DMassWL of optical clusters with centre position, radius and richness is used as input as well as a set of parameters twoDMassParams2DConvergenceClusters and an input FITS catalogue (e.g. lensMCCatalog2DMassWL). The output product is a tar file twoDMassConvergenceClusters containing FITS data arrays with as HDUs the noisy/denoised convergence maps, a map containing the number of galaxies per pixel and SNR maps for a patch centered on each detected (and sufficiently massive) cluster. This tar is associated to the name of the catalogue processed to ease further joint processing. A single product for a given cluster twoDMassConvergenceSingleCluster can also be obtained (with same FITS structure as collected in the tar file)

Multi-Patch E/B convergence and SNR maps

For large fields processed on the sphere (for statistical studies, to compute the catalogue of peaks), a first option is to perform multiple local Cartesian projections on the sky and estimate E/B modes on these patches.

An input FITS catalogue (e.g. lensMCCatalog2DMassWL) collecting information from OU-SHE/OU-PHZ is needed, as well as a set of parameters twoDMassParams2DConvergencePatchesToSphere defining in particular the patch characteristics, output HEALPix resolution and other running parameters. Potentially a footprint mask VMPZIDHEALPixFootprintMask is also needed to define the decomposition of the sky. The output twoDMassConvergencePatchesToSphere collects as HEALPix maps the noisy/denoised convergence, number of galaxies per pixel, Monte Carlo realisations and SNR maps.

The input and output are described in:

Spherical E/B convergence and SNR maps

For large fields processed on the sphere (for statistical studies, to compute the catalogue of peaks), another option is to fully process the catalogue on the sphere. An input FITS catalogue (e.g. lensMCCatalog2DMassWL) collecting information from OU-SHE/OU-PHZ is needed, as well as a set of parameters twoDMassParams2DConvergenceSphere specifying the HEALPix resolution and other running parameters.

The product twoDMassConvergenceSphere contains as products computed on the sphere (HEALPix maps) the noisy/denoised convergence maps, the number of galaxies per pixel, Monte Carlo realisations and SNR maps.

Peak Catalogs

Finally, peak catalogues twoDMassPeakCatalog can be estimated for statistical studies. This is performed by detecting peaks either from aperture mass maps computed on the input catalogue (e.g. lensMCCatalog2DMassWL) with parameters twoDMassParamsPeakCatalogMassAperture2D or from the wavelet decomposition of the convergence maps using parameters twoDMassParamsPeakCatalogConvergence, along with SNR maps. A FITS keyword indicate which estimator has been used.