Visual description of the multi-tracer products that can be simulated with CoLoRe. The upper plot shows one of the beams used internally by CoLoRe for domain decomposition, with the redshift and angular coordinates of 3 different DESI samples, and maps of the density, radial velocity, and lensing shear constructed from sources at two different redshifts. The lower plot shows the density skewers calculated for three arbitrary DESI quasars contained in the same beam, as a function of comoving distance.Abstract
We present \texttt{CoLoRe}, a public software package designed to efficiently generate synthetic realizations of multiple cosmological surveys. \texttt{CoLoRe} can simulate the growth of structure with different levels of accuracy, with the current implementation supporting lognormal fields as well as first- and second-order Lagrangian perturbation theory. The code simulates the density field on an all-sky light cone up to a chosen maximum redshift and uses it to generate a variety of two- and three-dimensional observables, including galaxy positions and velocities, gravitational lensing maps (shear, magnification, and convergence), the integrated Sachs–Wolfe effect, line-intensity mapping, and line-of-sight skewers for simulations of the Lyman-$\alpha$ forest. We test the accuracy of the simulated maps against analytical theoretical predictions and demonstrate the performance of the code through a multi-survey simulation that includes DESI galaxies and quasars, LSST galaxies and weak lensing, and SKA intensity mapping and radio galaxies. We expect \texttt{CoLoRe} to be particularly useful for studies aimed at characterizing the impact of systematics in multi-experiment analyses, quantifying the covariance between different data sets, and testing cross-correlation pipelines for upcoming cosmological surveys.