Galaxy--point spread function correlations as a probe of weak-lensing systematics with UNIONS data

Weak gravitational lensing requires precise measurements of galaxy shapes and therefore accurate knowledge of the point spread function (PSF) model. The latter can be a source of systematics that affect the shear two-point correlation function. A key aspect of weak-lensing analysis is the forecasting of the systematics due to the PSF. Correlation functions of galaxies and the PSF, the so-called ̊ho and τ statistics, are used to evaluate the level of systematics coming from the PSF model and PSF corrections and contributing to the two-point correlation function used to perform cosmological inference. Our goal is to introduce a fast and simple method to estimate this level of systematics and to assess its agreement with state-of-the-art approaches. We introduce a new way to estimate the covariance matrix of τ statistics using analytical expressions. The covariance allows us to estimate parameters directly related to the level of systematics associated with the PSF and provides us with a tool to validate the PSF model used in a weak-lensing analysis. We applied these methods to data from the Ultraviolet Near-Infrared Optical Northern Survey (UNIONS). We show that semi-analytical covariance yields results comparable to those obtained by using covariances obtained from simulations or jackknife resampling. The approach requires less computation time and is therefore well suited to rapid comparison of the systematic level obtained from different catalogues. We also show how one can break degeneracies between parameters with a redefinition of the τ statistics. The methods developed in this work will be useful tools in the analysis of current weak-lensing data but also of Stage IV surveys such as Euclid LSST, and Roman . They provide fast and accurate diagnostics on PSF systematics that are crucial in the context of cosmic shear studies.