UNIONS: The impact of systematic errors on weak-lensing peak counts

 

Authors: E. Ayçoberry, V. Ajani, A. Guinot, M. Kilbinger, V. Pettorino, S. Farrens, J.-L. Starck, R. Gavazzi, M. Hudson
Journal: A&A
Year: 2022
DOI:  
Download: ADS | arXiv


Abstract

Context. The Ultraviolet Near-Infrared Optical Northern Survey (UNIONS) is an ongoing deep photometric multi-band survey of the Northern sky. As part of UNIONS, the Canada-France Imaging Survey (CFIS) provides r-band data which we use to study weak-lensing peak counts for cosmological inference.
Aims. We assess systematic effects for weak-lensing peak counts and their impact on cosmological parameters for the UNIONS survey. In particular, we present results on local calibration, metacalibration shear bias, baryonic feedback, the source galaxy redshift estimate, intrinsic alignment, and the cluster member dilution.

Methods. For each uncertainty and systematic effect, we describe our mitigation scheme and the impact on cosmological parameter constraints. We obtain constraints on cosmological parameters from MCMC using CFIS data and MassiveNuS N-body simulations as a model for peak counts statistics.
Results. Depending on the calibration (local versus global, and the inclusion of the residual multiplicative shear bias), the mean matter density parameter Ωm can shift up to −0.024 (−0.5σ). We also see that including baryonic corrections can shift Ωm by +0.027 (+0.5σ) with respect to the DM-only simulations. Reducing the impact of the intrinsic alignment and cluster member dilution through signal-to-noise cuts can lead to a shift in Ωm of +0.027 (+0.5σ). Finally, with a mean redshift uncertainty of ∆z ̄ = 0.03, we see that the shift of Ωm (+0.001 which corresponds to +0.02σ) is not significant.

Conclusions. This paper investigates for the first time with UNIONS weak-lensing data and peak counts the impact of systematic effects. The value of Ωm is the most impacted and can shift up to ∼ 0.03 which corresponds to 0.5σ depending on the choices for each systematics. We expect constraints to become more reliable with future (larger) data catalogues, for which the current pipeline will provide a starting point. The code used to obtain the results is available in the following Github repository.

ShapePipe: a new shape measurement pipeline and weak-lensing application to UNIONS/CFIS data

 

Authors: A. Guinot, M. Kilbinger, S. Farrens, A. Peel, A. Pujol, M. Schmitz, J.-L. Starck, T. Erben, R. Gavazzi, S. Gwyn, M. Hudson,  H. Hiledebrandt, T. Liaudat , et. al
Journal: A&A
Year: 2022
DOI:  
Download: ADS | arXiv


Abstract

UNIONS is an ongoing collaboration that will provide the largest deep photometric survey of the Northern sky in four optical bands to date. As part of this collaboration, CFIS is taking r-band data with an average seeing of 0.65 arcsec, which is complete to magnitude 24.5 and thus ideal for weak-lensing studies. We perform the first weak-lensing analysis of CFIS r-band data over an area spanning 1700 deg2 of the sky. We create a catalogue with measured shapes for 40 million galaxies, corresponding to an effective density of 6.8 galaxies per square arcminute, and demonstrate a low level of systematic biases. This work serves as the basis for further cosmological studies using the full UNIONS survey of 4800 deg2 when completed. Here we present ShapePipe, a newly developed weak-lensing pipeline. This pipeline makes use of state-of-the-art methods such as Ngmix for accurate galaxy shape measurement. Shear calibration is performed with metacalibration. We carry out extensive validation tests on the Point Spread Function (PSF), and on the galaxy shapes. In addition, we create realistic image simulations to validate the estimated shear. We quantify the PSF model accuracy and show that the level of systematics is low as measured by the PSF residuals. Their effect on the shear two-point correlation function is sub-dominant compared to the cosmological contribution on angular scales <100 arcmin. The additive shear bias is below 5x104, and the residual multiplicative shear bias is at most 103 as measured on image simulations. Using COSEBIs we show that there are no significant B-modes present in second-order shear statistics. We present convergence maps and see clear correlations of the E-mode with known cluster positions. We measure the stacked tangential shear profile around Planck clusters at a significance higher than 4σ.

ShapePipe: A modular weak-lensing processing and analysis pipeline

 

Authors: S. Farrens, A. Guinot, M. Kilbinger, T. Liaudat , L. Baumont, X. Jimenez, A. Peel , A. Pujol , M. Schmitz, J.-L. Starck, and A. Z. Vitorelli
Journal: A&A
Year: 2022
DOI: 10.1051/0004-6361/202243970
Download: ADS | arXiv


Abstract

We present the first public release of ShapePipe, an open-source and modular weak-lensing measurement, analysis, and validation pipeline written in Python. We describe the design of the software and justify the choices made. We provide a brief description of all the modules currently available and summarise how the pipeline has been applied to real Ultraviolet Near-Infrared Optical Northern Survey data. Finally, we mention plans for future applications and development. The code and accompanying documentation are publicly available on GitHub.

AMICO galaxy clusters in KiDS-DR3: measurement of the halo bias and power spectrum normalization from a stacked weak lensing analysis

 

Authors: L. Ingoglia, G. Covone, M. Sereno, ..., S. Farrens, et al.
Journal: MNRAS
Year: 2022
DOI: 10.1093/mnras/stac046
Download: ADS | arXiv


Abstract

Galaxy clusters are biased tracers of the underlying matter density field. At very large radii beyond about 10 Mpc/\textit{h}, the shear profile shows evidence of a second-halo term. This is related to the correlated matter distribution around galaxy clusters and proportional to the so-called halo bias. We present an observational analysis of the halo bias-mass relation based on the AMICO galaxy cluster catalog, comprising around 7000 candidates detected in the third release of the KiDS survey. We split the cluster sample into 14 redshift-richness bins and derive the halo bias and the virial mass in each bin by means of a stacked weak lensing analysis. The observed halo bias-mass relation and the theoretical predictions based on the \\Lambda\CDM standard cosmological model show an agreement within \2\sigma\. The mean measurements of bias and mass over the full catalog give \M_{200c} = (4.9 \pm 0.3) \times 10^{13} M_{\odot}/\textit{h}\ and \b_h \sigma_8^2 = 1.2 \pm 0.1\. With the additional prior of a bias-mass relation from numerical simulations, we constrain the normalization of the power spectrum with a fixed matter density \\Omega_m = 0.3\, finding \\sigma_8 = 0.63 \pm 0.10\.

Deep transfer learning for blended source identification in galaxy survey data

 

Authors: S. Farrens, A. Lacan, A. Guinot, A. Z. Vitorelli
Journal: A&A
Year: 2022
DOI: 10.1051/0004-6361/202141166
Download: ADS | arXiv


Abstract

We present BlendHunter, a proof-of-concept deep-transfer-learning-based approach for the automated and robust identification of blended sources in galaxy survey data. We take the VGG-16 network with pre-trained convolutional layers and train the fully connected layers on parametric models of COSMOS images. We test the efficacy of the transfer learning by taking the weights learned on the parametric models and using them to identify blends in more realistic Canada-France Imaging Survey (CFIS)-like images. We compare the performance of this method to SEP (a Python implementation of SExtractor) as a function of noise level and the separation between sources. We find that BlendHunter outperforms SEP by ~ 15% in terms of classification accuracy for close blends (< 10 pixel separation between sources) regardless of the noise level used for training. Additionally, the method provides consistent results to SEP for distant blends (>10 pixel separation between sources) provided the network is trained on data with noise that has a relatively close standard deviation to that of the target images. The code and data have been made publicly available to ensure the reproducibility of the results.

PySAP: Python Sparse Data Analysis Package for Multidisciplinary Image Processing

 

Authors: S. Farrens, A. Grigis, L. El Gueddari, Z. Ramzi, Chaithya G. R., S. Starck, B. Sarthou, H. Cherkaoui, P.Ciuciu, J-L. Starck
Journal: Astronomy and Computing
Year: 2020
DOI: 10.1016/j.ascom.2020.100402
Download: ADS | arXiv


Abstract

We present the open-source image processing software package PySAP (Python Sparse data Analysis Package) developed for the COmpressed Sensing for Magnetic resonance Imaging and Cosmology (COSMIC) project. This package provides a set of flexible tools that can be applied to a variety of compressed sensing and image reconstruction problems in various research domains. In particular, PySAP offers fast wavelet transforms and a range of integrated optimisation algorithms. In this paper we present the features available in PySAP and provide practical demonstrations on astrophysical and magnetic resonance imaging data.


Code

PySAP Code


Euclid preparation: VI. Verifying the Performance of Cosmic Shear Experiments

Euclid preparation: VI. Verifying the Performance of Cosmic Shear Experiments

Authors: Euclid Collaboration, P. Paykari, ..., S. Farrens, M. Kilbinger, V. Pettorino, S. Pires, J.-L. Starck, F. Sureau, et al.
Journal: Astronomy and Astrophysics
Year: 2020
DOI:  10.1051/0004-6361/201936980
Download:

ADS | arXiv

 


Abstract

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

Euclid preparation. V. Predicted yield of redshift 7 < z < 9 quasars from the wide survey

Euclid preparation: V. Predicted yield of redshift 7

Authors: Euclid Collaboration, R. Barnett, ..., S. Farrens, M. Kilbinger, V. Pettorino, F. Sureau, et al.
Journal: Astronomy and Astrophysics
Year: 2019
DOI:  10.1051/0004-6361/201936427
Download:

ADS | arXiv

 


Abstract

We provide predictions of the yield of 7<z<9 quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; Φ) with redshift, Φ∝10k(z−6), k=−0.72, and a further steeper rate of decline, k=−0.92; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to JAB∼23. Quasars at z>8 may be selected from Euclid OYJH photometry alone, but selection over the redshift interval 7<z<8 is greatly improved by the addition of z-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond z=6. For the case that the decline of the QLF accelerates beyond z=6, with k=−0.92, Euclid should nevertheless find over 100 quasars with 7.0<z<7.5, and ∼25 quasars beyond the current record of z=7.5, including ∼8 beyond z=8.0. The first Euclid quasars at z>7.5 should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, 7<z<8, M1450<−25, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at JAB∼23. The precision with which k can be determined over 7<z<8 depends on the value of k, but assuming k=−0.72 it can be measured to a 1 sigma uncertainty of 0.07.

Euclid preparation III. Galaxy cluster detection in the wide photometric survey, performance and algorithm selection

 

Authors: Euclid Collaboration, R. Adam, ..., S. Farrens, et al.
Journal: A&A
Year: 2019
DOI: 10.1051/0004-6361/201935088
Download: ADS | arXiv


Abstract

Galaxy cluster counts in bins of mass and redshift have been shown to be a competitive probe to test cosmological models. This method requires an efficient blind detection of clusters from surveys with a well-known selection function and robust mass estimates. The Euclid wide survey will cover 15000 deg2 of the sky in the optical and near-infrared bands, down to magnitude 24 in the H-band. The resulting data will make it possible to detect a large number of galaxy clusters spanning a wide-range of masses up to redshift ∼2. This paper presents the final results of the Euclid Cluster Finder Challenge (CFC). The objective of these challenges was to select the cluster detection algorithms that best meet the requirements of the Euclid mission. The final CFC included six independent detection algorithms, based on different techniques, such as photometric redshift tomography, optimal filtering, hierarchical approach, wavelet and friend-of-friends algorithms. These algorithms were blindly applied to a mock galaxy catalog with representative Euclid-like properties. The relative performance of the algorithms was assessed by matching the resulting detections to known clusters in the simulations. Several matching procedures were tested, thus making it possible to estimate the associated systematic effects on completeness to <3%. All the tested algorithms are very competitive in terms of performance, with three of them reaching >80% completeness for a mean purity of 80% down to masses of 1014 M⊙ and up to redshift z=2. Based on these results, two algorithms were selected to be implemented in the Euclid pipeline, the AMICO code, based on matched filtering, and the PZWav code, based on an adaptive wavelet approach.

A Distributed Learning Architecture for Scientific Imaging Problems

 

Authors: A. Panousopoulou, S. Farrens, K. Fotiadou, A. Woiselle, G. Tsagkatakis, J-L. Starck,  P. Tsakalides
Journal: arXiv
Year: 2018
Download: ADS | arXiv


Abstract

Current trends in scientific imaging are challenged by the emerging need of integrating sophisticated machine learning with Big Data analytics platforms. This work proposes an in-memory distributed learning architecture for enabling sophisticated learning and optimization techniques on scientific imaging problems, which are characterized by the combination of variant information from different origins. We apply the resulting, Spark-compliant, architecture on two emerging use cases from the scientific imaging domain, namely: (a) the space variant deconvolution of galaxy imaging surveys (astrophysics), (b) the super-resolution based on coupled dictionary training (remote sensing). We conduct evaluation studies considering relevant datasets, and the results report at least 60\% improvement in time response against the conventional computing solutions. Ultimately, the offered discussion provides useful practical insights on the impact of key Spark tuning parameters on the speedup achieved, and the memory/disk footprint.