The concept of blind analysis, a key procedure to remove the human-based systematic error called confirmation bias, has long been an integral part of data analysis in many research areas. In cosmology, blind analysis is recently making its entrance, as the field progresses into a fully fledged high-precision science. The credibility, reliability and robustness of results from future sky-surveys will dramatically increase if the effect of confirmation bias is kept under control by using an appropriate blinding procedure. Here, we present a catalog-level blinding scheme for galaxy clustering data apt to be used in future spectroscopic galaxy surveys. We shift the individual galaxy positions along the line of sight based on 1) a geometric shift mimicking the Alcock-Paczynski effect and 2) a perturbative shift akin to redshift-space distortions. This procedure has several advantages. After combining the two steps above, it is almost impossible to accidentally unblind. The procedure induces a shift in cosmological parameters without changing the galaxies' angular positions, hence without interfering with the effects of angular systematics. Since the method is applied at catalog level, there is no need to adopt several blinding schemes tuned to different summary statistics, likelihood choices or types of analyses. By testing the method on mock catalogs and the BOSS DR12 catalog we demonstrate its performance in blinding galaxy clustering data for relevant cosmological parameters sensitive to the background expansion rate and the growth rate of structures.

We present a tomographic cosmic shear analysis of the Kilo-Degree Survey (KiDS) combined with the VISTA Kilo-Degree Infrared Galaxy Survey (VIKING). This is the first time that a full optical to near-infrared data set has been used for a wide-field cosmological weak lensing experiment. This unprecedented data, spanning $450~$deg$^2$, allows us to improve significantly the estimation of photometric redshifts, such that we are able to include robustly higher-redshift sources for the lensing measurement, and - most importantly - solidify our knowledge of the redshift distributions of the sources. Based on a flat $\Lambda$CDM model we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.737_{-0.036}^{+0.040}$ in a blind analysis from cosmic shear alone. The tension between KiDS cosmic shear and the Planck-Legacy CMB measurements remains in this systematically more robust analysis, with $S_8$ differing by $2.3\sigma$. This result is insensitive to changes in the priors on nuisance parameters for intrinsic alignment, baryon feedback, and neutrino mass. KiDS shear measurements are calibrated with a new, more realistic set of image simulations and no significant B-modes are detected in the survey, indicating that systematic errors are under control. When calibrating our redshift distributions by assuming the 30-band COSMOS-2015 photometric redshifts are correct (following the Dark Energy Survey and the Hyper Suprime-Cam Survey), we find the tension with Planck is alleviated. The robust determination of source redshift distributions remains one of the most challenging aspects for future cosmic shear surveys.

We present constraints on Horndeski gravity from a combined analysis of cosmic shear, galaxy-galaxy lensing and galaxy clustering from $450\,\mathrm{deg}^2$ of the Kilo-Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey. The Horndeski class of dark energy/modified gravity models includes the majority of universally coupled extensions to $\Lambda$CDM with one scalar field in addition to the metric. We study the functions of time that fully describe the evolution of linear perturbations in Horndeski gravity. Our results are compatible throughout with a $\Lambda$CDM model. By imposing gravitational wave constraints, we fix the tensor speed excess to zero and consider a subset of models including e.g. quintessence and $f(R)$ theories. Assuming proportionality of the Horndeski functions $\alpha_B$ and $\alpha_M$ (kinetic braiding and the Planck mass run rate, respectively) to the dark energy density fraction $\Omega_{\mathrm{DE}}(a) = 1 - \Omega_{\mathrm{m}}(a)$, we find for the proportionality coefficients $\hat{\alpha}_B = 0.20_{-0.33}^{+0.20} \,$ and $\, \hat{\alpha}_M = 0.25_{-0.29}^{+0.19}$. Our value of $S_8 \equiv \sigma_8 \sqrt{\Omega_{\mathrm{m}}/0.3}$ is in better agreement with the $Planck$ estimate when measured in the enlarged Horndeski parameter space than in a pure $\Lambda$CDM scenario. In our joint three-probe analysis we report a downward shift of the $S_8$ best fit value from the $Planck$ measurement of $\Delta S_8 = 0.016_{-0.046}^{+0.048}$ in Horndeski gravity, compared to $\Delta S_8 = 0.059_{-0.039}^{+0.040}$ in $\Lambda$CDM. Our constraints are robust to the modelling uncertainty of the non-linear matter power spectrum in Horndeski gravity. Our likelihood code for multi-probe analysis in both $\Lambda$CDM and Horndeski gravity is publicly available at http://github.com/alessiospuriomancini/KiDSHorndeski .