Issue 
A&A
Volume 627, July 2019



Article Number  A137  
Number of page(s)  17  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201834629  
Published online  12 July 2019 
Testing gravity with galaxygalaxy lensing and redshiftspace distortions using CFHTStripe 82, CFHTLenS, and BOSS CMASS datasets^{⋆}
^{1}
AixMarseille Univ., CNRS, CNES, LAM, Marseille, France
email: eric.jullo@lam.fr
^{2}
AixMarseille Univ., CNRS/IN2P3, CPPM, Marseille, France
^{3}
Dipartimento di Fisica e Scienze della Terra, Università degli Studi di Ferrara, Via Saragat 1, 44122 Ferrara, Italy
^{4}
INAF – Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, Via Gobetti 93/3, 40129 Bologna, Italy
^{5}
Dipartimento di Fisica e Astronomia, Alma Mater Studiorum Università di Bologna, Via Gobetti 93/2, 40129 Bologna, Italy
^{6}
INFN – Sezione di Bologna, Viale Berti Pichat 6/2, 40127 Bologna, Italy
^{7}
MaxPlanckInstitut für extraterrestrische Physik, Giessenbachstrasse 1, 85748 Garching bei München, Germany
^{8}
Shanghai Astronomical Observatory (SHAO), Nandan Road 80, Shanghai 200030, PR China
^{9}
Centro Brasileiro de Pesquisas Físicas, Rio de Janeiro, RJ 22290180, Brasil
^{10}
Institute of Physics, Laboratory of Astrophysics, Ecole Polytechnique Fédérale de Lausanne (EPFL), Observatoire de Sauverny, 1290 Versoix, Switzerland
^{11}
Instituto de Astrofísica de Andalucía (CSIC), Glorieta de la Astronomía, 18080 Granada, Spain
^{12}
Departamento de Física Teórica, Módulo 15, Universidad Autónoma de Madrid, 28049 Madrid, Spain
^{13}
Centro de Investigación Avanzada en Física Fundamental (CIAFF), Universidad Autónoma de Madrid, 28049 Madrid, Spain
^{14}
LeibnizInstitut für Astrophysik (AIP), An der Sternwarte 16, 14482 Potsdam, Germany
Received:
12
November
2018
Accepted:
2
May
2019
The combination of galaxygalaxy lensing (GGL) and redshift space distortion of galaxy clustering (RSD) is a privileged technique to test general relativity predictions and break degeneracies between the growth rate of structure parameter f and the amplitude of the linear power spectrum σ_{8}. We performed a joint GGL and RSD analysis on 250 sq. deg using shape catalogues from CFHTLenS and CFHTStripe 82 and spectroscopic redshifts from the BOSS CMASS sample. We adjusted a model that includes nonlinear biasing, RSD, and Alcock–Paczynski effects. We used an Nbody simulation supplemented by an abundance matching prescription for CMASS galaxies to build a set of overlapping lensing and clustering mocks. Together with additional spectroscopic data, this helps us to quantify and correct several systematic errors, such as photometric redshifts. We find f(z = 0.57) = 0.95 ± 0.23, σ_{8}(z = 0.57) = 0.55 ± 0.07 and Ω_{m} = 0.31 ± 0.08, in agreement with Planck cosmological results 2018. We also estimate the probe of gravity E_{G} = 0.43 ± 0.10, in agreement with ΛCDM−GR predictions of E_{G} = 0.40. This analysis reveals that RSD efficiently decreases the GGL uncertainty on Ω_{m} by a factor of 4 and by 30% on σ_{8}. We make our mock catalogues available on the Skies and Universe database.
Key words: cosmological parameters / cosmology: observations / largescale structure of Universe
The catalogues are available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/627/A137 and at http://www.skiesanduniverses.org
© E. Jullo et al. 2019
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Since its inception, general relativity theory (GR) has been constantly tested, starting with observations in the solar system and in our Galaxy (see e.g. Damour 2000). At cosmological scales, the advent of wide field survey experiments currently yields very high precision measurements in both the early and late ages of the universe. A Universe dominated by cold dark matter and a cosmological constant in the context of GR (hereafter ΛCDM−GR model) reproduces all these observations with very high accuracy and for this reason, the model is often referred to as the standard reference model.
However, some slight tensions are emerging between predictions based on the cosmic microwave background measurements from the Planck mission at redshift z = 1089 and measurements at redshifts z < 1 obtained from galaxy clustering or gravitational lensing. In particular with Planck, the amplitude of the matter power spectrum σ_{8} is larger and the Hubble constant H_{0} is smaller than what is estimated at redshifts z < 1 at about 2σ confidence level (C.L.; Planck Collaboration XIII 2016; Beutler et al. 2014; Alam et al. 2017; Hildebrandt et al. 2017; DES Collaboration 2018). Although systematic errors in the analyses can explain a significant fraction of these discrepancies, they might nonetheless suggest some issues with our understanding and modelling of the expansion of the Universe, or of the largescale structure formation probed by galaxy clustering and gravitational lensing.
The common approach to test ΛCDM−GR at cosmological scales is either to measure the expansion history H(z) of the Universe (e.g. Betoule et al. 2014; Alam et al. 2017; Magaña et al. 2015), or to measure the growth of structures traced by the velocity or density fields in redshift space (e.g. de la Torre et al. 2013; Tully et al. 2016; Martinet et al. 2018). In this paper, we combine galaxygalaxy lensing (GGL) and RSD to test both aspects simultaneously at redshift z = 0.57. The amplitude of GGL measurements is sensitive to H(z) and the density field, whereas RSD probes the growth of structure f(z) through galaxy peculiar velocities. The combination of these two observables has demonstrated its effectiveness at isolating the independent effects of the growth rate of structure f(z), the amplitude of the matter power spectrum σ_{8}, and the dark energy equation of state parameter w involved in H(z) calculation (Simpson et al. 2013; de la Torre et al. 2017; Joudaki et al. 2017).
Zhang et al. (2007) proposed an alternative method to test deviations to GR. Assuming small scalar perturbations around the Friedmann–Lemaître–Robertson–Walker metric (FLRW) in the conformal Newtonian gauge ds^{2} = −a(τ)^{2}[1 + 2Ψ]dτ^{2} + a(τ)^{2}[1 − 2Φ]dx^{2}, where a is a scale factor, τ is the conformal time, and x are comoving coordinates, these authors proposed a statistics E_{G} sensitive to the gravitational slip between the two gravitational potentials Φ and Ψ as follows:
where all quantities are estimated at the redshift of interest. Reyes et al. (2010) proposed an associated observational estimator E_{G} = Υ_{gm}/βΥ_{gg} (see below in Sect. 6.4), which converges to ⟨E_{G}⟩ in the largescale limit where the galaxy bias b and the distortion parameter β = f/b converge to constant values. The smallscale filtered galaxymatter crosscorrelation Υ_{gm} probed with GGL is sensitive to both b and ∇^{2}(Φ − Ψ) since photons traverse equal quantity of space and time. The galaxyvelocity crosscorrelation βΥ_{gg} probed with RSD is sensitive to galaxy bias and the Newtonian potential Ψ. In GR and in absence of anisotropic stress, Φ = −Ψ so lensing is sensitive of 2∇^{2}Φ. In the linear regime, the Poisson equation relates the potential to the matter density contrast δ, such that . This estimator therefore converges to in the standard model.
In their seminal paper, Zhang et al. (2007) predicted deviations from GR with four alternative models: ΛCDM, flat (Dvali et al. 2000, hereafter DGP), f(R) gravity (Bean et al. 2007), and TeVeS/MOND (Bekenstein 2004). Apart from the TeVeS/MOND model, which introduces a wavelength difference between dynamical and lensing powerspectra, all other models add at most 10% deviations compared to GR predictions. Leonard et al. (2015) reached similar conclusions with other models based on the empirical extension of the Poisson equations with the commonly used parameters Σ(a) and μ(a) (Amendola et al. 2008; Ferreira & Skordis 2010). Most importantly, these authors found that details of the analysis (e.g. integration length along the line of sight for projected estimators) could mimic deviations similar to those predicted with alternative models of gravity, thus the need for a careful study of these biases. In any case with 20%–30% precision, current datasets are not yet at the level of accuracy required to observe these deviations, and unsurprisingly no deviation to GR predictions has been detected so far (Blake et al. 2016; Pullen et al. 2016; de la Torre et al. 2017; Alam et al. 2017; Amon et al. 2018).
Nowadays, cosmological analyses require measurements with exquisite control of systematic errors, at all levels from data acquisition to cosmological model inference. The wide range of expertise needed to reach the requirements is demonstrated by the size of the ongoing and forthcoming cosmological experiments such as Dark Energy Survey (Dark Energy Survey Collaboration 2005), the Kilo Degree Survey (Hildebrandt et al. 2017), the HyperSuprime Cam survey (Aihara et al. 2018), the extended Baryonic Oscillation Sky Survey (Dawson et al. 2013), the Prime Focus Spectrograph project (Sugai et al. 2012), the Dark Energy Survey Instrument project (DESI Collaboration 2016a,b), the Large Scale Synoptic Telescope (LSST Dark Energy Science Collaboration 2012), and the Euclid mission Laureijs et al. (2011).
In this paper, we extend the Leauthaud et al. (2017) analysis (hereafter L17), by adding RSD measurements of CMASS galaxies from the Baryon acoustic Oscillation Spectroscopic Survey (BOSS) to GGL measurements in the CFHTStripe 82 and CFHTLS fields. Thanks to refined simulations, we precisely quantify systematic errors, and thus manage to reconcile real and simulated measurements of clustering and lensing. The work presented builds on the theoretical model and joint RSD and GGL analysis developed in de la Torre et al. (2017, hereafter DLT17).
The outline of the paper is as follows. First we present our galaxy bias model, and its inclusion in standard clustering and lensing estimators. Next, we present our datasets and measurement estimators. Our tests on simulations are presented in Sect. 5, and our estimates of the cosmological parameters in Sect. 6. Finally, we present our measurement of E_{G} and conclude. Systematic errors are discussed in the Appendix. Unless otherwise mentioned, we express the GGL projected densities Σ_{gm} and distances in comoving coordinates. We assume the fiducial ΛCDM−GR cosmology with flat universe, Ω_{m} = 0.31, h = 0.6777, Ω_{b} = 0.048, σ_{8} = 0.82 (Planck Collaboration XIII 2016).
2. Method
In the following, we compute the RSD twopoint galaxy correlation functions in configuration space. We decomposed the threedimensional galaxy separation vector s into polar (s, μ) or Cartesian (r_{p}, π) coordinates in the frame defined by the line of sight and the normal to it, where s is the norm of s, μ is the cosine of the angle between s and the line of sight, π and r_{p} are the projections of s on the line of sight and its normal, respectively. In the flatsky approximation, the transformation between Cartesian and polar coordinates is μ = π/s, (Fisher et al. 1994). Conversely, the GGL formalism is defined in real space, where the separation vector r is decomposed into Cartesian coordinates (R, χ), where χ and R are the projections of r on the line of sight and its normal, respectively. In GGL, the radial window function of integration is hundreds of h^{−1} Mpc, and the effects of RSD can safely be neglected (Baldauf et al. 2010). Hereafter, we assume that R in the model corresponds to r_{p} in the observations.
2.1. Galaxy bias model
In this work, we want to measure the growth rate f and amplitude of the matter power spectrum σ_{8} with GGL and galaxyclustering measurements. These measurements are not typically estimated at the same scale. While the GGL signal typically emerges in the range of transverse distances 0.1 < r_{p} < 20 h^{−1} Mpc, the galaxy clustering signal rises between 10 < r_{p} < 100 h^{−1} Mpc. To maximize the overlap between these two observables in the nonlinear regime, we adopted the fourth order perturbation model in the initial density field as proposed by McDonald & Roy (2009). Assuming homogeneity and isotropy in the density field, they derived the following expression for the halomatter power spectrum:
where P_{δδ} and P_{lin} represent the nonlinear and linear matter power spectra respectively; P_{b2, δ} and P_{bs2, δ} are the oneloop power spectra between the density field δ, its derivative and the variance of the tidal tensor field s(x). The term includes various third order terms of the galaxy bias model (see McDonald & Roy 2009, for more details). Assuming coevolution between the halo and matter density fields, and the bias being purely local in Lagrangian space at initial conditions, Baldauf et al. (2012) computed the second order halo density field in both Eulerian and Lagrangian space and found the relation b_{s2} = −4/7(b_{1} − 1). Under the same assumptions as above to compute b_{s2}, Saito et al. (2014) obtained the relation b_{3nl} = 32/315(b_{1} − 1). The analytical expressions for all these terms are given in Appendix A of DLT17.
2.2. Galaxygalaxy lensing model
The measured GGL differential excess surface density is defined as
where the mean projected surface density can be read as
and Σ_{gm}(R) is the projected surface density defined as a function of the galaxymatter crosscorrelation function (Dvornik et al. 2018)
where the mean matter density is constant in comoving coordinates. The galaxymatter crosscorrelation function ξ_{gm} is obtained from the Fourier transform of the galaxymatter power spectrum P_{gm}(k) defined above.
In practice, we used an FFTLOG unbiased Hankel transform with parameter in logarithmic space to perform the Fourier transform^{1}. We truncated the power spectrum at k_{min} = 10^{−5} and k_{max} = 1000 to minimize cutoff aliasing during the FFT operation, and we splineinterpolated the resulting correlation function to obtain the desired binning.
2.3. Redshift space distortions model
In this work, we used the Taruya et al. (2010) model to describe the RSD effect. In the ideal case in which galaxies are perfect tracers of the matter density field, this model takes the form
where θ is the divergence of the velocity field defined as θ = −∇⋅v/(aHf). The values P_{δδ}, P_{θθ}, and P_{δθ} are the nonlinear matter density, velocity divergence, and density–velocity divergence power spectra, respectively; C_{A}(k, μ, f) and C_{B}(k, μ, f) terms derive from the general anisotropic power spectrum of matter and their expressions are given in Taruya et al. (2010) and de la Torre & Guzzo (2012).
The damping function D(kμσ_{v}) essentially (but not only) describes the Fingers of God effect on the twopoint correlation function, and we modelled it as a Lorentzian damping in Fourier space, i.e.
where σ_{v} represents an effective pairwise velocity dispersion that we fitted for and then treated as a nuisance parameter.
This model can be generalized to the case of biased tracers, by including our bias model. Hence, we obtain (Beutler et al. 2014; GilMarín et al. 2014)
where,
In the above equations P_{b2, δ}, P_{bs2, δ}, P_{b2, b2}, P_{b2, bs2}, P_{bs2, bs2} and are oneloop integrals, of which analytical expressions can be found in Appendix A of DLT17. We computed the linear matter power spectrum P_{lin} using the CLASS Bolzmann code (Lesgourgues 2011), and the nonlinear matter power spectrum P_{δδ} using the semianalytic prescriptions HALOFIT (Smith et al. 2003; Takahashi et al. 2012). To predict the velocity spectra P_{θθ} and P_{δθ}, we use the nearly universal fitting functions from Bel et al. (2019), already used in DLT17 and Pezzotta et al. (2017). These are built such that they converge to P_{lin} at large scales, but reproduce nonlinearities at small scales. Pezzotta et al. (2017) highlighted that adding a redshift dependency with σ_{8}(z) such that
and
was helping. The coefficients (p_{1} = 1.906, p_{2} = 2.163, p_{3} = 2.972, p_{4} = 2034) were deduced from a fit to measurements performed on the DEMNUni simulations (dark energy and massive neutrino universe). These two fitting functions are accurate within 5% to the measurements in simulations and appear to be insensitive to the presence of neutrinos (Carbone et al. 2016). The overall degree of nonlinearity in these terms is therefore solely controlled by σ_{8}(z), which is left free when fitting the model to observations. Although these fitting functions possibly duplicate a fraction of the highorder modes included in the perturbation theory model above, we demonstrate in DLT17 and in Sect. 5 below that it does not bias significantly our cosmological estimates given data uncertainties.
Finally, we obtain the multipole moments of the anisotropic correlation functions in configuration space
where j_{ℓ}(x) is the spherical Bessel function and is the anisotropic powerspectrum multipole moment of order ℓ defined as
where L_{ℓ}(x) are the Legendre polynomial of order ℓ.
At linear scales, f and σ_{8} are degenerate, but extending to nonlinear scales with the Taruya et al. (2010) model, , , and appear in the calculation of the correction terms C_{A} and C_{B}, and hence help break the degeneracy. Accordingly, in our model (f, b_{1}, b_{2}, σ_{v}, σ_{8}) are treated as separate parameters in the fit (de la Torre & Guzzo 2012).
2.4. Spectroscopic redshift uncertainties
It is worth mentioning that redshift errors can potentially affect the anisotropic RSD signal. They have the same effect as galaxy random motions in virialized objects. We model the redshift errors by multiplying the anisotropic power spectrum by the Fourier transform of a Gaussian damping function of the form
such that our predicted signal can be finally written as
Bolton et al. (2012) measured the error on the estimated spectroscopic velocities, thanks to multiple observations of the same CMASS galaxies, and found approximately δ_{v} = 32 km s^{−1}, which translates to σ_{z} = 0.797 h^{−1} Mpc in comoving distances at redshift z = 0.57 with our fiducial cosmology. This effect is therefore negligible, but we included it to have a cleaner estimate of σ_{v}.
2.5. Suppressing smallscale modelling uncertainties
Although considered as sufficient for galaxyclustering analysis, we find that our weaklensing (WL) model deviates from our measurements with simulated data at scales r_{p} ∼ 3 h^{−1} Mpc (see Fig. 7 in DLT17). In order to damp the contribution of any signal below a given cutoff radius R_{0}, we computed the annular differential excess surface density (ASAD) estimator from the data (Baldauf et al. 2010). For the lensing observable ΔΣ_{gm}(r_{p}), it is given by
and for the galaxy clustering
These two estimators become useful to estimate E_{G} in the following. We derive the projected correlation w_{p}(r_{p}) from the projection of the multipole decomposition of the correlation function in redshift space
The α_{2ℓ} coefficients are given in Baldauf et al. (2010)
We integrate along the line of sight up to π_{max} = 40 h^{−1} Mpc to match the integration length used with the data (see the estimators Sect. 4.2). According to Singh et al. (2019), they found consistent results whether they use π_{max} = 50 h^{−1} Mpc or 100 h^{−1} Mpc. Given the low number CMASS galaxies in this analysis, we set π_{max} = 40 h^{−1} Mpc to minimize the noise.
The ASAD can also be predicted from theory. For the lensing part, Υ_{gm}(r_{p}, R_{0}) is obtained by filtering the crosscorrelation function ξ_{gm}(r)
with the window function W_{Υ}(x, r_{p}, R_{0}) (Baldauf et al. 2010) defined as
where Θ(x) is the Heaviside step function. In a similar manner, we computed Υ_{gg}(r_{p}, R_{0}) by simply replacing ξ_{gm}(r) by ξ_{gg}(r) in Eq. (23). We included the RSD effect in the calculation of Υ_{gg}(r_{p}). In both cases, we integrated in logarithmic scale up to r_{max} = 100 h^{−1} Mpc.
We note that we do not include intrinsic alignment in our modelling. This choice is motivated by the marginal constraints obtained in Joudaki et al. (2017) on the amplitude of this effect , with smallscale cut on γ_{t} at θ > 12 arcmin. Since we applied the smallscale Υ filter, we anticipate very little constraint on this parameter as well, at a significant additional computing cost.
2.6. Alcock–Paczynski effect
We may mention that additional distortions can occur in the correlation functions owing to possible differences between the true and the fiducial cosmological models used to compute the distances. This effect was first identified by Alcock & Paczynski (1979; hereafter AP) as a means to constrain the cosmological model. However these distortions degenerate with the RSD effect and considerably limit the constraining power of the AP effect (Ballinger et al. 1996; Matsubara & Suto 1996). Fortunately, the scaledependence of the AP and RSD effects differ and thus allow us this degeneracy to break (Seo & Eisenstein 2003; Blake et al. 2011; Chuang & Wang 2012).
In this work, we adopted the AP model proposed by Xu et al. (2013). The isotropic and anisotropic distortions are expressed respectively as
where quantities computed with the fiducial cosmology as denoted with primes. Those parameters modify the transverse and the radial distances such that
Given these distortions, the observed redshiftspace monopole and quadrupole expressed in configuration space become
The GGL estimator becomes
3. Data
In our GGL analysis, the lenses are the CMASS galaxies and the sources are galaxies in the CFHTLens and CFHTStripe 82 WL catalogues. Lenses have spectroscopic redshifts and sources have photometric redshifts. For each lens, we can then discard all uncorrelated foreground sources and use the background sources to estimate the lensing signal. The final GGL measurement is the average of the signals for each lens.
3.1. Weaklensing datasets
3.1.1. CFHTLens catalogue
In 2013, the CFHTLenS team released a public WL catalogue covering an area of 154 sq. deg in four wide fields (W1, W2, W3, and W4; Erben et al. 2013; Heymans et al. 2012). So far, the depth of the input CFHT Legacy Survey imaging is unrivaled, with a 5σ point source limiting magnitude i_{AB} ∼ 25.5. The LENSFIT algorithm is used to measure the shape of every object detected with i_{AB} < 24.7. Then, we selected galaxies with good shape measurements (FITCLASS = 0 and WEIGHT > 3).
We obtained photometric redshifts from five optical band photometry u, g, r, i, z and reach a precision of about 5% up to z ∼ 1 (Hildebrandt et al. 2012). Galaxygalaxy lensing measurements can be significantly biased by inaccurate photometric redshifts (Nakajima et al. 2012). We computed the photometric redshift bias estimator ⟨b_{z}⟩, based on spectroscopic and photometric catalogues matched in position, and averaged over the CMASS redshift distribution (see appendix details). Since the spectroscopic calibration sample is significantly shallower than the photometric sample, we discarded galaxies fainter than the 90% completeness limit of the spectroscopic sample (see below), i.e. we only kept galaxies brighter than i_{AB} < 24. After this selection, we obtained ⟨b_{z}⟩= + 0.003 ± 0.003, ⟨b_{z}⟩= − 0.014 ± 0.004 and ⟨b_{z}⟩= + 0.022 ± 0.003 in fields W1, W3, and W4, respectively. We discarded field W2 because it only contains 200 CMASS galaxies on its northern edge.
Our final catalogue contains 3.5 millions galaxies over an effective area of about 127 sq. deg. The galaxy density^{2} is n_{eff} = 7.0 galaxies arcmin^{−2}. The median redshift is z_{med} = 0.70.
3.1.2. CFHTStripe 82 catalogue
The CFHTStripe 82 survey (CS82; Moraes et al. 2014) is an iband imaging survey containing 173 tiles (PIs: J.P. Kneib, A. Leauthaud, M. Makler, L. Van Waerbeke). It covers about 160 sq. deg of the Sloan Digital Sky Survey (SDSS, Gunn et al. 2006) stripe 82 region, with a 5σ pointsource magnitude limit i_{AB} ∼ 24.1, and a mean seeing of 0.6″. The effective area is 129.2 sq. deg after masking out bright stars and other image artefacts (L17). We used a new version 3.0 of the shape catalogue that has shapes measured with LENSFIT down to magnitude i_{AB} < 24.7. This new version benefits from internal calibration in LENSFIT based on image simulations inherited from the CFHTLenS project. Shape measurements are accurate at the 2% level, without relying on any additional linear correction. In addition, this new catalogue contains about 40% more galaxies, mostly because of a better handling of galaxy deblending and instrument artefacts in LENSFIT (L. van Waerbeke, priv. comm.).
Photometric redshifts in the original version of the catalogue (Bundy et al. 2015) were computed with BPZ (Benítez 2000) using ugriz from the Stripe 82 coadds (Annis et al. 2014) and UJHK from UKIDSS. We use nearestneighbour interpolation in sky coordinates, i magnitude, and g − r, r − i, i − z colour space to get photometric redshifts for the new galaxies. We verified that the redshift distribution is unchanged. We applied the same procedure as in the CFHTLS fields to estimate the bias due to photometric redshifts in our GGL measurements. However, given the relatively shallow spectroscopic survey coverage of the CS82 field compared to CFHTLS fields (90% completeness reached at i_{AB} = 22.5), we are forced to select galaxies only down to i_{AB} < 22.5. Although this cut is quite severe, it allow us to confidently model and correct photometric redshift bias in this field. The lack of deeper spectroscopic information prevents us from exploiting the complete WL catalogue. For i_{AB} < 22.5 sources and CMASS lenses, we find a bias b_{z} = −0.028 ± 0.006. In contrast to L17, we apply no cut based on the ODD quality flag because we find it has no impact on our lensing measurements given our stringent cut in magnitude. Our final catalogue contains 2.2 million galaxies. The galaxy density is n_{eff} = 4.7 galaxies arcmin^{−2}. The median redshift is z_{median} = 0.53.
3.2. Spectroscopic dataset: BOSS CMASS sample
The BOSS spectroscopic survey (Eisenstein et al. 2011) is a programme of the SDSS project. The constant (stellar) mass (CMASS) galaxy sample is one of the galaxy samples observed in this survey. This sample consists of galaxies selected with SDSS photometry, such that they lie in the redshift range 0.43 < z < 0.7 and represent a sample of galaxies approximately volumelimited in stellar mass (Reid et al. 2016). Early clustering analysis found that CMASS galaxies lie in massive haloes and have a mean halo mass of 2.6 × 10^{13} h^{−1} M_{⊙}, a largescale bias of b ∼ 2.0, and a satellite fraction of 10% (White et al. 2011).
We used the public DR12v5 version of the CMASS catalogue (Alam et al. 2015). The galaxy surface density is about 100 deg^{−2} (Reid et al. 2016). We only considered CMASS overlapping with our four lensing fields, i.e. covering an area 250 sq. deg. Our catalogue of lenses contains 28 039 CMASS galaxies, distributed as reported in Table 1. The redshift distribution of CMASS galaxies compared to CS82 and CFHTLens lensing sources is shown in Fig. 1.
Number of CMASS galaxies per field, effective lensing area after masking and number of WL sources.
Fig. 1.
Redshift distribution of CMASS galaxies (blue) compared to CS82 (orange) and CFHTLens (green) source distributions, after WL selection has been performed. WL n(z) are based on photometric redshifts (see text for details). 
In spite of a careful photometric selection, the observed CMASS galaxy sample remains contaminated by various observational effects (Ross et al. 2012). We took these into account by applying the galaxy weights w_{g} = w_{star}w_{see}(w_{zf} + w_{cp} − 1) as defined in Ross et al. (2017). We also included the Feldman et al. (1994, hereafter FKP) weights with the parameter P_{0} = 20 000 h^{−3} Mpc^{3} (Ross et al. 2012), such that the noise in the power spectrum is minimum at the baryon acoustic oscillation scale k = 0.1 h Mpc^{−1}. Although not optimal for our study focussed on smallscale clustering, this value of P_{0} allows for a consistent comparison with previous measurements. For consistency, we take the same value of P_{0} for our mock catalogues and data. Finally, we used the DR12v5random0 catalogues trimmed to the regions overlapping with WL data.
4. Measurement estimators
4.1. Galaxygalaxy lensing estimation
We computed ΔΣ_{gm}(r_{p}) using the estimator
where r_{p} is the comoving transverse distance between the lens and the source at redshifts z_{l} and z_{s}, respectively. The subscript “r” denotes the random catalogue of lensing objects. Our number of random objects N_{r} is ten times the number of lenses N_{l}. Their redshift distribution n(z) is that from CMASS galaxies (Nuza et al. 2013). The subtraction of the random signal decreases the variance at large scales (Singh et al. 2017; Shirasaki et al. 2017). The value ϵ_{+} represents the tangential component of a source ellipticity around a lens. The weight is the product of the shape measurement weight w_{s} from LENSFIT and the critical density. This inverse variance scheme downweights pairs which are close in redshift (Mandelbaum et al. 2006). The critical lensing density Σ_{cr}(z_{l}, z_{s}) in comoving units is defined as
where D_{S}, D_{LS}, D_{L} are the observersource, lenssource, and observerlens angular diameter distances^{3} We used the bestfit estimate of the photometric redshift to compute the distances, instead of the full probability distribution, as suggested in Blake et al. (2016). However, our approach described below and based on full raytracing simulations consistently takes this simplification into account.
4.2. Anisotropic galaxy clustering estimation
We computed the twopoint galaxy correlation function in the polar and Cartesian coordinate systems. The anisotropy in the signal is due to the RSD effect we are after. The estimator is the same in each coordinate system and is defined as
where (x, y)=(s, μ) or (r_{p}, π). GG, GR and RR are the normalized number of pairs between galaxygalaxy, galaxyrandom, and random–random, respectively, at a given separation.
We compressed the information contained in ξ(s, μ) by projecting it on the Legendre polynomials using the expressions for the correlationfunction multipole moments
where L_{ℓ} is the Legendre polynomial of order ℓ. We used the monopole and quadrupole ℓ = (0, 2) only because the higher order nonnull multipoles are too noisy.
We also computed the projected correlation function w_{p}(r_{p}) by projecting ξ(r_{p}, π) along the line of sight such that
where we find the optimal value π_{max} = 40 h^{−1} Mpc to minimize the noise due to the limited number of pairs in our fields.
4.3. Joint lensing and clustering likelihood
We performed a maximumlikelihood analysis to derive the cosmological parameters from the GGL and RSD measurements. In each field i, we measured the data vector and we computed the likelihood function per field ℒ^{i} such that
where m is the model prediction, and is the precision matrix estimated from the simulations.
Our four fields are statistically uncorrelated, and therefore the global likelihood is just the product of the individual likelihoods for each field, i.e.
Field W4 partly overlaps with field S82, but this overlapping represents < 6% of the total area. In addition, CFHTLens catalogue used for W4 goes deeper than CS82 catalogue used for S82, thus decreasing further the correlation between the two fields.
5. Simulations
5.1. Light cones and lensing mock catalogues
In order to accurately estimate largescale variance and possibly unveil new systematic errors, we produce light cones with the same geometry as the observed fields. We used the BigMultidark Nbody simulation, as it appears to be a good compromise between particle resolution and cosmological volume (m_{p} = 2.5 × 10^{10}h^{−1} M_{⊙}, L_{box} = 2.5 h^{−1} Gpc, Planck cosmology with h = 0.6777; Klypin et al. 2016). Following the approach described in Giocoli et al. (2016), we simulated four fields, CS82, W1, W3, and W4, which have light cones extending up to redshift z = 2.3 for the CFHTLS fields and z = 2 for CS82. We computed lensing properties, such as deflected positions, shear, and convergence, by ray tracing through 25 lens planes separated by 161 h^{−1} Mpc comoving (Giocoli et al. 2016) using the GLAMER^{4} code. The spatial resolution of the lensing maps is 6 arcsec.
5.1.1. Lensing properties
We simulated lensing catalogues of sources including survey mask, intrinsic shape, and photometric redshift noises. For survey mask, we placed the source galaxies are the location of observed sources. Thus, we naturally reproduced the footprint and the holes around bright stars and other artefacts of the real WL catalogue. Effects due to the intrinsic clustering of sources in projection were also included. We got the shear properties for each source by spatially interpolating the values from the shear maps computed with GLAMER. For sources intrinsic ellipticities, we randomly drew observed ellipticites ϵ^{obs} from the WL catalogue, that we multiply by a random orientation ϕ^{int}, such that and .
5.1.2. Photometric redshifts
To simulate photometric redshifts with catastrophic failures, we designed a method related to that described in Lima et al. (2008), also referred as the direct calibration method (DIR) in the Kilo Degree Survey (KIDS, Hildebrandt et al. 2017). We start by estimating the true redshift distribution n_{true}(z) for our CFHTLens (i_{AB} < 24) and CS82 (i_{AB} < 22.5) WL catalogues from our spectroscopic calibration sample described in Appendix A. In practice, we computed the histograms of the WL and spectroscopic (ZP) catalogues in the magnitudecolor space (i, g − r, r − i, i − z), that we limited to the region ([18; 25],[−1; 3],[−1; 3],[−1; 3]). We applied the same binning for both catalogues. For each bin of coordinates m, we derive the weights W(m)=N_{WL}(m)/N_{ZS}(m), where N is the number of sources per bin. We assume all sources in a bin have the same weight. Finally, we obtain the true distribution in redshift bin i with the following sum:
A drawback of this approach is that if spectroscopic selection does not cover part of the redshift range, then it truncates n_{true}(z). However, we see in the following that the coverage is sufficient for our purpose.
Then, we compute the joint probability P(z_{BPZ}, z_{spec}) for each field, as shown in Fig. 2. We observe that the spectroscopic redshift completeness at z < 0.5 in field W1 and W4 is very low because most of the redshifts come from the CMASS sample. Fortunately, this has little impact on our simulation of photometric redshift noise because our analysis focusses on the crosscorrelation of CMASS galaxies with lensing sources at z > 0.5. We also observe that the scatter in the z_{BPZ} of CS82 field is almost twice as large as in field W3, and differs between the three CFHTLens fields. This justifies our fieldbyfield treatment of the photometric redshift noise. Finally, we assign photometric redshifts to the simulated sources by randomly drawing a photometric redshift from P(z_{BPZ}, z_{spec}), where we assume the spectroscopic redshift z_{spec} is the true redshift that was assigned at the beginning of the procedure.
Fig. 2.
Probability distribution of having a photometric redshift with BPZ and a spectroscopic redshift for each field. Contours are given for 1, 2, and 3σ C.L. 
5.2. Spectroscopic CMASS mock catalogues
We adopted the subhalo abundance matching (SHAM) procedure described in RodríguezTorres et al. (2016) to produce CMASS mock catalogues. Starting from the ROCKSTAR public catalogues (Behroozi et al. 2013)^{5}, we computed a scattered peak velocity , where 𝒩 is the normal distribution, and σ_{SHAM} = 0.31. We also simulated the CMASS incompleteness in stellar mass and redshift, based on the stellar mass function (SMF) from the Portsmouth SEDFIT DR12 stellar mass catalogue with Kroupa initial mass function (Maraston et al. 2013). We binned the catalogue in 12 redshift intervals between 0.43 < z < 0.7 and in 18 stellar mass bins between 10.5 < log_{10}(M * /M_{⊙}) < 12.3. Thus, we obtained a tabulated SMF that we can interpolate in stellar mass and redshift. Finally from cumulative stellar mass and halo mass functions, we constructed a number density matching such that . Since different cosmologies were assumed in the Portsmouth catalogue and in the BigMultidark simulations, h = 0.73 and h = 0.6777 respectively, we renormalized the stellar masses to the BigMultidark cosmology. As shown in Fig. 3, our number densities for each of the four fields are in good agreement with the measurements from Anderson et al. (2012).
Fig. 3.
Left panel: number density of CMASS mock galaxies for our 4 simulated fields. Limits of our analysis are indicated with blue dashed lines. Measurements from Anderson et al. (2012, A12) are in grey. Right panel: CMASS CMF for the 4 fields reproducing the observed incompleteness. The mock catalogue is complete at high mass in agreement with the model proposed in RodríguezTorres et al. (2016) in grey. 
We also include the effect of peculiar velocities by summing together in redshiftspace the halo position r_{c} and the peculiar velocity vector v in real space using the relation , where is the lineofsight unit vector, a is the scale factor, and H(z_{real}) is the Hubble parameter at redshift z_{real}, the redshift corresponding to r_{c}. Finally, we masked the borders of the square simulated fields W3 and W4 to reproduce their complex geometry, and we computed the FKP weights, assuming the same P_{0} = 20 000 h^{−3} Mpc^{3} as in the data. Since data have been corrected for fiber collision, redshift failure, stellar density, and seeing, we did not simulate these effects.
5.3. Bias due to photometric redshift noise
We computed successively the lensing signal for catalogues with and without photometric redshift noise, and compare the measurements in Fig. 4. We find that the large photometric scatter observed in field S82 (Fig. 2) seems to result in a bias of about 10% in the lensing signal at scales R < 10 h^{−1} Mpc, whereas the CFHTLens fields seem insignificantly affected. We argue that this might explain the discrepancy highlighted in L17 between lensing measurements obtained with real and mock data. Indeed, in the following, we show that our lensing measurements with mock data contaminated by photometric redshift noise are in agreement with real data.
Fig. 4.
Comparison of lensing measurements performed on simulated catalogues affected and not affected by photometric redshifts noise and ΔΣ, respectively. The grey shaded areas correspond to the uncertainties on the mean value obtained by resampling the multiple noises in different light cones. 
5.4. Bias from smallscale modelling
We used the simulation to quantify the bias in the estimation of the cosmological parameters f and Ω_{m} due to our model prediction of the small scales. Successively, we cut data points of ξ_{0} and ξ_{2} at scales s_{min} = 11.2, 14.1, and 17.8 h^{−1} Mpc, and Υ_{gm} at scales R_{0} = 1.0 and 1.5 h^{−1} Mpc. Overall, we find that the values s_{min} = 17.8 h^{−1} Mpc and R_{0} = 1 h^{−1} Mpc provide the best compromise between systematic bias and statistical precision as can be seen in Fig. 5.
Fig. 5.
Bias between recovered parameters f and Ω_{m} relative to the values used in the mocks as functions of the minimum scale of the multipoles s_{min}, and the cutoff radius of the Υ(R, R_{0}) lensing estimator. Values of s_{min} = 17.8 h^{−1} Mpc and R_{0} = 1.0 h^{−1} Mpc provide unbiased estimates of f and Ω_{m}. All measurements were performed without tapering smoothing in the covariance matrices. 
5.5. Covariance matrices
To obtain an unbiased estimate of the precision matrices, we need minimal errors in the covariance matrices and therefore a large number of mock catalogues. Noise in the covariance matrices increases the errors on the model parameter estimation (see e.g. Taylor & Joachimi 2014). Unfortunately, we were limited by the size of our simulation box L = 2 h^{−1} Gpc. Escoffier et al. (2016) proposed a method to increase the number of mocks, based on Jackknife resampling of the mock catalogues (see Table 2). Following their prescription, we split each catalogue into N_{JK} spatial subregions and measured the clustering and lensing signals in each Jackknife subsample using estimators given in Eqs. (34) and (32). The covariance matrix for each mock catalogue is then
Properties of the simulated fields in terms of independent mock catalogue, random resampling of lensing shape noise and photometric redshifts per catalogue, and number of subregions for Jackknife resampling.
where the mean vector is obtained from the Jackknife samples
In addition, given our limited number of independent mock catalogue N_{m}, we increase their number for lensing by resampling N_{r} times the observed lensing ellipticity distribution function, and the photometric redshifts distribution. We find this strategy to efficiently improve the accuracy of the covariance matrix for the lensing, especially at small scales. The final covariance matrix is therefore obtained by averaging the Jackknife covariance matrices
Finally, we computed the precision matrix
Escoffier et al. (2016) have shown that this expression provides an unbiased estimate of the true precision matrix.
In spite of our resampling strategy, our covariance matrices are still noisy. Therefore, we adopted the tapering method proposed by Paz & Sánchez (2015) to damp the noise by a filter function beyond a given tapering scale T_{p}. This technique is based on the assumption that correlation between pairs of data points far apart is negligible and little information is lost by treating these points as being independent. Although very efficient, it is commonly accepted that this method might inadvertently remove nonGaussian terms (Paz & Sánchez 2015). However this effect is beyond the scope of this analysis given our data and the range of scales investigated in this work. In Fig. 6, we observe that large tapering yields errors similar to no tapering. In contrast, small tapering zeros all offdiagonal terms, and can also lead to overestimated errors. We find the errors on f and Ω_{m} to reach a minimum value at the tapering scale T_{p} ∼ 12 h^{−1} Mpc. We adopted this scale in the rest of this analysis. We should note that all measurements were performed with s_{min} = 14.1 h^{−1} Mpc and R_{0} = 1.5 h^{−1} Mpc. However, we repeated some measurements with our final setup (s_{min} = 17.8 h^{−1} Mpc and R_{0} = 1.0 h^{−1} Mpc) and find that these parameters have almost no impact on the tapering scale behaviour. The covariance and precision matrices obtained before and after tapering at this scale are shown in Fig. 7. We can observe that the noise in the offdiagonal terms is significantly reduced after tapering. This is particularly obvious between clustering and lensing, which cover very different range of scales.
Fig. 6.
Variation of the relative errors on the parameters f and Ω_{m}, as a function of the smooth scale T_{p} in the covariance matrices. There is no improvement below T_{p} = 12 h^{−1}. 
Fig. 7.
Matrices of covariance (top panel) and precision (bottom panel) in logarithmic color scale for the 4 fields used in this analysis. In both panels, the upper triangular part of the matrices represents the case without tapering, while the lower part represents the case with tapering T_{p} = 12 h^{−1} Mpc. Noise between far apart scales is significantly decreased and the errors on the model parameters converge to a minimum. 
6. Cosmological results
The quality tests and errors assessment that we performed with the simulations give us confidence that our dataset can lead to reliable cosmological constraints.
6.1. Galaxyclustering and galaxygalaxy lensing measurements
In Figs. 8 and 9, we show our RSD and GGL measurements, along with our theoretical predictions, assuming the fiducial parameters of the simulation, and a constant linear bias b_{1} = 1.8. We find a good agreement within 1σ C.L. between mocks, data, and theoretical predictions for all fields. We notice that the quadrupole of the correlation function measurement in the field W3 is lower than the 1σ C.L., and that the GGL measurement in the field W4 is lower than 1σ C.L. at scales R < 1 h^{−1} Mpc. For field W3, we found that setting σ_{8}(z = 0.57) = 0.9 and b_{1} = 1.5 could reconcile predictions with measurements, thus suggesting a sample variance effect. These values are within the 3σ C.L. of the RSDonly fit of the data (see Fig. 10). For field W4, we attribute the discrepancy to our poor modelling of baryonic or lensing effects at small scales, which average out too slowly in the data to reproduce the simulated darkmatter only profile. Nonetheless, the overall good agreement gives us confidence that we can proceed with the cosmological analysis.
Fig. 8.
Monopole (red) and quadrupole (blue) measurements with mock catalogues (shaded region), real data (solid lines), and theoretical predictions with a linear bias parameter b_{1} = 1.8 (dashed lines). Black dots represent prereconstruction measurements with the full DR12v5 CMASS sample from Cuesta et al. (2016). 
Fig. 9.
Filtered Υ_{gm} and nonfiltered ΔΣ GGL measurements with mocks (shaded regions), ΔΣ and Υ data (blue and orange points respectively), and theory with a linear bias parameter b_{1} = 1.8 (dashed line). Black dots in S82 panel represent ΔΣ measurements from L16, and Υ_{gm} measurements from Alam et al. (2017) in CFHTLS panels. 
6.2. Growth of structure and background constraints
We estimate the cosmological parameters f, σ_{8}, and Ω_{m} by combining ξ_{0}, ξ_{2}, and Υ_{gm} measurements. The power of this combination to break the degeneracy between f and σ_{8} has already been demonstrated (see e.g. DLT17, Joudaki et al. 2017). In this paper, we move one step further by estimating Ω_{m} as well from the data. Figure 10 shows the independent lensing, clustering, and combined constraints on these parameters. Bestfit values and 1σ error estimates are reported in Table 3. A corner plot with all the parameters involved in the fit is reported in Fig. A.2. On the one hand, we find that GGL alone constrains Ω_{m} at 45% and σ_{8} at 22%. It provides no constraint on the structure growth rate f. On the other hand, RSD also constrains σ_{8} at 20% but leaves Ω_{m} completely unconstrained as expected from the model. When used in combination, GGL and RSD measurements yield 12% precision constraint on σ_{8}, i.e. almost as if the two datasets were independent. In fact, Fig. 10 shows that the wellknown WL degeneracy between Ω_{m} and σ_{8} intersects almost perpendicularly with the constraint on σ_{8} from RSD.
Fig. 10.
Improvement on estimating Ω_{m} and σ_{8} when combining RSD and WL measurements. The blue and black curves are, respectively, obtained with WL and RSD constraints only. The orange curves are obtained with the combination of WL and RSD. Contours are given at 1, 2, and 3σ C.L. Vertical lines indicate Planck TT,TE,EE+lowE 2018 results. 
Bestfit and derived parameters obtained by fitting the RSD only, GGL only, and their combination.
In Fig. 11, we present our estimate of the growth rate f, and compare to other measurements. In spite of having a wider area, we obtain a constraint similar to that found in DLT17 with VIPERS. Clearly, the number of RSD tracers determines the precision. In both analysis, we have about 28,000 galaxies in the range 0.5 < z < 0.7. Regarding WL, the number densities of background sources at z > 0.7 in both analysis are similar. We have n_{eff} = 3.45 in CFHTLS fields and n_{eff} = 2.33 in CS82 and CFHTLS fields combined.
Fig. 11.
Growth rate f as a function of redshift compared to recent measurements. The black line and surrounding grey shared area indicate the Planck TT,TE,EE+lowE 2018 mean and 1σ uncertainty predictions for ΛCDM−GR flat model. 
We also compare our results with analyses performed on the full CMASS sample. Singh et al. (2019) performed a joint analysis with Planck cosmic microwave background lensing and SDSS galaxy lensing and obtained three times tighter constraints than ours. Their results are in agreement with ours at the 1σ C.L.
Finally, combining CMASS power spectrum and bispectrum, GilMarín et al. (2017) also obtained very competitive constraints at redshift z = 0.57 in agreement with ours. These two estimates find a tension on f with Planck predictions at z = 0.57. Interestingly, this tension was also observed in other RSD analysis with the CMASS sample, but not with the LOWZ sample (e.g. Alam et al. 2017; Beutler et al. 2014).
6.3. Comparison with other measurements
From MCMC, we can derive new parameter constraints, defined as a combination of single parameters. In particular, we look at the quantity , very common in gravitational lensing analyses. We find S_{8} = 0.72 ± 0.08, which is in agreement with the value estimated in L17, but 2 − 3σ smaller than the cosmic microwave background measurements S_{8} = 0.832 ± 0.013 (Planck Collaboration VI 2018). Similarly, our estimate of σ_{8} = 0.73 ± 0.08 is 2 − 3σ smaller than the measurement σ_{8} = 0.8111 ± 0.0060 from the Planck collaboration 2018. Our results are also in agreement with KIDS shear peaks statistics S_{8} = 0.75 ± 0.059 (Martinet et al. 2018; Shan et al. 2018), KIDS tomographic WL S_{8} = 0.745 ± 0.039 (Hildebrandt et al. 2017), and DES cosmological constraints from WL and clustering . We note that our fit only performed with RSD measurements yield an estimate of σ_{8} = 0.78 ± 0.26, in better agreement with Planck estimates.
The linear galaxy bias parameter is known to be degenerate with the cosmological parameters Ω_{m} and σ_{8}. In our fitting procedure, we assume a uniform prior on b_{1} between 1 and 3, which largely encompasses the expected value for the CMASS sample. In their clustering analyses, GilMarín et al. (2017) found b_{1}σ_{8}(z = 0.57) = 1.237 ± 0.011, and Chuang et al. (2013) found b_{1}σ_{8}(z = 0.57) = 1.18 ± 0.14. We find b_{1}σ_{8}(z = 0.57) = 1.256 ± 0.097 in full agreement with these previous measurements. Marginalizing over σ_{8}, we find b_{1} = 2.33 ± 0.33, in agreement with White et al. (2011) and subsequent analyses (e.g. Ho et al. 2012; Nuza et al. 2013; RodríguezTorres et al. 2016).
Our model also contains a second order biasing term, but our estimated value b_{2} = −0.04 ± 0.53 is not sufficient to discuss the nonlinearity of the CMASS sample. We note that GilMarín et al. (2017) found b_{2} = 0.606 ± 0.069, which is in agreement with us.
Finally, we also include Alcock–Paczynski effect in our model, but found no significant constraint given the data, α = 1.01 ± 0.05 and ϵ = 0.00 ± 0.05. We note that no significant constraint could either be obtained by GilMarín et al. (2017) with the full CMASS DR12 sample.
To conclude, we demonstrated the effectiveness of combining RSD and GGL to break the degeneracies between the amplitude of the largescale structure fluctuations σ_{8} and their growth rate f at redshift z = 0.57. We also found that the constraints on the cosmic matter density Ω_{m}, usually derived with WL, could be significantly improved by combining with RSD. Given the data, our measurements are still in agreement with Planck predictions.
6.4. Measuring E_{G}
To corroborate the information obtained with the analysis in the previous section and probe any deviation to ΛCDM−GR predictions, we estimate E_{G}, as defined in Reyes et al. (2010). The E_{G} estimator is function of projected scale r_{p}, and is defined as (Zhang et al. 2007)
This estimator is particularly interesting because it apparently just relies on observations. However, we show in the following that this might not be the case in practice.
Indeed, the E_{G} estimator suffers from a few downsides. First, this estimator relies on a previous determination of β. However, statistical and systematic error propagation into E_{G} error might be awkward, unless proper correction terms and covariance matrices are determined from ad hoc mock catalogues of lensing and clustering. Although seldom the case in the past, this is becoming more and more common (Blake et al. 2016; Amon et al. 2018; Singh et al. 2019).
Second, it is assumed that galaxy bias is linear, scaleindependent, and the galaxy density field is fully correlated to the underlying matter density field, i.e. the crosscorrelation factor r_{cc} = 1. Of course, these assumptions hold in the linear regime, but the scale at which they break depends on the galaxy sample. Using CMASS mock catalogues, several authors have shown that they hold in the range 5 < r_{p} < 60 h^{−1} Mpc (Baldauf et al. 2010; White et al. 2011; Amon et al. 2018; Singh et al. 2019). This depends on the requested precision on the model though, and recent works have proposed to take nonlinearity and other effects into account with normalizing functions derived from simulations (Alam et al. 2017; Singh et al. 2019). The multiplication of these correction terms nonetheless tend to reveal the limitation of the E_{G} estimator.
Marta Pinho et al. (2018) have noted that E_{G} depends not only on gravity but also on the background (e.g. quantified with the matter density Ω_{m0} in ΛCDM). Although it is always possible to predict E_{G} for different cosmological models (see e.g. Zhang et al. 2007, in which predictions are computed for ΛCDM, Flat DGP, f(R) gravity, TeVeS/MOND), a discrepancy with the observations therefore does not specifically point to a failure of GR, but can also be attributed to the background. In this respect, these authors claim that an estimator such as η, based on independent estimates of fσ_{8}(z), H(z), E_{G} might be more appropriate. To our point of view, adjusting an actual model including modified gravity parameters might be as helpful.
In spite of these limitations, E_{G} has become quite popular recently, mostly because of the advent of wide field imaging and spectroscopic surveys. This estimator has been measured several times, but no significant deviation from ΛCDM−GR has been found so far. In particular with the CMASS sample at redshift z = 0.57, Amon et al. (2018) found E_{G} = 0.26 ± 0.08, Blake et al. (2016) found E_{G} = 0.30 ± 0.07, Pullen et al. (2016) found E_{G} = 0.24 ± 0.06, Alam et al. (2017) found E_{G} = 0.42 ± 0.06, and Singh et al. (2019) found E_{G} = 0.39 ± 0.05. The dispersion in the estimates reveal that the method is probably not fully mature yet, and deserves further investigation, in particular regarding the observational biases such as photometric redshifts.
Figure 12 shows our measurements of E_{G} as a function of scale. We estimate β = 0.41 ± 0.15 from our fit to the RSD measurements only. Although it is difficult to compare our work with other works because authors use different models, this value is larger but statistically consistent with that found with the full CMASS sample β = 0.34 ± 0.02 from Amon et al. (2018). For each bin of E_{G}(r_{p}), we add in quadrature the errors on the ratio derived from the data and the error on β = f/b_{1} derived from the fit, with the following chain rule formula:
Fig. 12.
Measurement of E_{G} with combined constraints in the fields CFHTStripe 82 and CFHTLs W1, W3, and W4. The horizontal black line indicates the Planck TT,TE,EE+lowE 2018 prediction. We note that CFHTS82 data help shrink the error bars by about 30%. CFHTLens points have been shifted rightwards for clarity. 
Using the MCMC samples from the fit of the GGL and RSD measurements, we use our model to reconstruct the ratio U^{i}. We also determine the correlation coefficients ρ(β, U^{i})∼0.3, i.e. β and U^{i} are significantly correlated.
We average in the scale range 10 < r_{p} < 60h^{−1} Mpc, and we find E_{G} = 0.48 ± 0.15 for CFHTLens field only, and E_{G} = 0.43 ± 0.11 for CFHTLens and CS82 fields combined, i.e. a 30% improvement in precision for a 100% increase in area. In the average E_{G} calculation, we consider the full covariance matrix between the E_{G} points estimated from our simulations in Sect. 5.1. We note finally that our current precision does not justify applying scaledependent bias, redshift weighting, or integration window corrections since their effect is less than 5% at the scales we consider (see Alam et al. 2017; Singh et al. 2019).
To put our measurement in context, we collected the E_{G} measurements at different redshifts from the literature in Fig. 13. Overall, we observe a trend of E_{G} values lower than predicted by Planck 2018. In the appendix, we forward model the E_{G} signal based on the MCMC samples output from the joint fit of the GGL and RSD measurements on mocks. Figure A.1 shows that the probability distribution function of the E_{G} estimator is skewed towards low values. Taking its mean value then necessarily leads to a biasedlow estimation of E_{G}(r_{p}). This result confirms the previous claim from Alam et al. (2017), and might also explain why so many E_{G} measurements are below the Planck 2018 predictions.
Fig. 13.
Measurements of E_{G} at various redshifts. Effective redshifts of the measurements have been slightly shifted for clarity. 
7. Conclusions
Understanding the current acceleration of the expansion of the Universe is one of the major goal of cosmology today. The combination of GGL and RSD is a remarkable avenue to distinguish the effect of gravity due to largescale structures, and the effect of some scalar field on the background expansion rate.
In this work, we have demonstrated the power of this combination applied to the wellstudied CMASS galaxy sample at the effective redshift z = 0.57. Using a comprehensive set of lensing and galaxy mock catalogues, we investigated several sources of systematic biases and determined the confidence limits for our datasets. In particular, we found that thanks to spectroscopic data, we could correct the bias due to photometric redshift uncertainty for galaxies brighter than i_{AB} < 22.5, and i_{AB} < 24 in our CFHTS82 and CFHTLens WL catalogues, respectively. These conservative magnitude cuts allow us to match our GGL measurements in the CFHTS82 and CFHTLS fields, although at the cost of drastically reducing the number of WL sources. This conclusive remark highlights the crucial need of spectroscopic redshifts to calibrate the photometric redshift faint galaxies.
Building on this encouraging result, we pursue a cosmological analysis of the combined dataset. Thanks to the joint GGL and RSD constraints, we efficiently break the degeneracy between galaxy bias b_{1}, matter density Ω_{m}, matter power spectrum amplitude σ_{8}, and the structure growth rate f at z = 0.57. We find astrophysical CMASS parameters and cosmological parameters in agreement with measurements previously obtained by other authors (White et al. 2011; Beutler et al. 2014; Chuang et al. 2013; GilMarín et al. 2017; Joudaki et al. 2017) and with Planck 2018 predictions in the frame of the ΛCDM−GR model.
Finally, we combine GGL and RSD measurements to estimate E_{G}. By averaging in the range of scales 10 < r_{p} < 60 h^{−1} Mpc, we find E_{G}(z = 0.57) = 0.43 ± 0.11, which is in perfect agreement with Planck 2018 prediction E_{G} = 0.40. Also, we use our mocks to characterize the statistical properties of E_{G}, and find that it has an asymmetric probability distribution, which tends to underestimate its mean value. This might explain part of the low values found in previous analysis. We also find that the reconstructed value of E_{G} = Ω_{m0}/f derived from the fit of the GGL and RSD measurements results in a value with smaller errors bars than that obtained directly from the data. More importantly, this value naturally includes the crosscorrelation terms between β and Υ_{gg}.
Back in 2012, Gaztañaga et al. (2012) was already advocating that overlapping lensing and spectroscopic surveys were 100 times more constraining on the dark energy equation of state and cosmic growth history parameter γ. Although it might not be the cleanest way to test gravity, the recent progress in estimating E_{G} at different redshifts with different tracers comes as a confirmation. In the future, wider imaging and spectroscopic surveys will result in very tight constraints on cosmological parameters. In contrast, it will probably take us more time to fully profit from smaller but deeper imaging surveys. Deep imaging surveys are helpful for many reasons, but also introduce additional systematic errors on the lensing side, in particular with respect to blending (HarnoisDéraps et al. 2018; Euclid Collaboration 2019). Nonetheless, both strategies lead to very exciting perspectives regarding our understanding of the dark sector.
We use the definition (Heymans et al. 2012), where w_{i} is a galaxy weight and Ω is the opening angle.
The factor (1 + z_{s})^{2} is missed in Eq. (10) of de la Torre et al. (2017), but was properly taken into account in the calculations.
Gravitational Lensing with Adaptive Mesh Refinement (Metcalf & Petkova 2014).
Acknowledgments
Based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, at the CanadaFranceHawaii Telescope (CFHT), which is operated by the National Research Council (NRC) of Canada, the Institut National des Science de l’Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii. The Brazilian partnership on CFHT is managed by the Laboratório Nacional de Astrofísica (LNA). We appreciate the support of the Laboratório Interinstitucional de eAstronomia (LIneA). We thank the CFHTLenS team. Funding for SDSSIII has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the US Department of Energy Office of Science. The SDSSIII web site is http://www.sdss3.org/. SDSSIII is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSSIII Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofisica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, Max Planck Institute for Extraterrestrial Physics, New Mexico State University, New York University, Ohio State University, Pennsylvania State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University. The BigMDPL simulation has been performed on the SuperMUC supercomputer at the LeibnizRechenzentrum (LRZ) in Munich, using the computing resources awarded to the PRACE project number 2012060963. We thank the Red Española de Supercomputación for granting us computing time in the Marenostrum Supercomputer at the BSCCNS where part of the analyses presented in this paper have been performed. We appreciate the support of the OCEVU Labex (Grant N^{o} ANR11LABX0060) and the A*MIDEX project (Grant N^{o} ANR11IDEX000102) funded by the Investissements d’Avenir French government programme managed by the ANR. We also acknowledge support from the ANR eBOSS project (ANR16CE310021) of the French National Research Agency. CG acknowledges support from Centre National d’Etudes Spatiales, Italian Ministry of Foreign Affairs and International Cooperation Directorate General for Country Promotion (Project “Crack the lens?”), from the agreement ASI n.I/023/12/0 “Attività relative alla fase B2/C per la missione Euclid”, and from the Italian Ministry for Education, University and Research (MIUR) through the SIR individual grant SIMCODE (project number RBSI14P4IH). JPK acknowledges support from the “Cosmology with 3DMaps of the Universe” SNF grant #175751. GY acknowledges financial support from MINECO/FEDER under project grant AYA201563810P
References
 Aihara, H., Arimoto, N., Armstrong, R., et al. 2018, PASJ, 70, S4 [NASA ADS] [Google Scholar]
 Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS, 219, 12 [NASA ADS] [CrossRef] [Google Scholar]
 Alam, S., Miyatake, H., More, S., Ho, S., & Mandelbaum, R. 2017, MNRAS, 465, 4853 [NASA ADS] [CrossRef] [Google Scholar]
 Alcock, C., & Paczynski, B. 1979, Nature, 281, 358 [NASA ADS] [CrossRef] [Google Scholar]
 Amendola, L., Kunz, M., & Sapone, D. 2008, JCAP, 2008, 013 [Google Scholar]
 Amendola, L., Kunz, M., Motta, M., Saltas, I. D., & Sawicki, I. 2013, Phys. Rev. D, 87, 023501 [NASA ADS] [CrossRef] [Google Scholar]
 Amon, A., Blake, C., Heymans, C., et al. 2018, MNRAS, 479, 3422 [NASA ADS] [CrossRef] [Google Scholar]
 Anderson, L., Aubourg, E., Bailey, S., et al. 2012, MNRAS, 427, 3435 [Google Scholar]
 Annis, J., SoaresSantos, M., Strauss, M. A., et al. 2014, ApJ, 794, 120 [NASA ADS] [CrossRef] [Google Scholar]
 Baldauf, T., Smith, R. E., Seljak, U., & Mandelbaum, R. 2010, Phys. Rev. D, 81, 063531 [NASA ADS] [CrossRef] [Google Scholar]
 Baldauf, T., Seljak, U., Desjacques, V., & McDonald, P. 2012, Phys. Rev. D, 86, 083540 [NASA ADS] [CrossRef] [Google Scholar]
 Ballinger, W. E., Peacock, J. A., & Heavens, A. F. 1996, MNRAS, 282, 877 [NASA ADS] [CrossRef] [Google Scholar]
 Bean, R., Bernat, D., Pogosian, L., Silvestri, A., & Trodden, M. 2007, Phys. Rev. D, 75, 064020 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Behroozi, P. S., Wechsler, R. H., & Wu, H.Y. 2013, ApJ, 762, 109 [NASA ADS] [CrossRef] [Google Scholar]
 Bekenstein, J. D. 2004, Phys. Rev. D, 70, 083509 [NASA ADS] [CrossRef] [Google Scholar]
 Bel, J., Pezzotta, A., Carbone, C., Sefusatti, E., & Guzzo, L. 2019, A&A, 622, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Benítez, N. 2000, ApJ, 536, 571 [NASA ADS] [CrossRef] [Google Scholar]
 Betoule, M., Kessler, R., Guy, J., et al. 2014, A&A, 568, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Beutler, F., Saito, S., Seo, H.J., et al. 2014, MNRAS, 443, 1065 [NASA ADS] [CrossRef] [Google Scholar]
 Blake, C., Glazebrook, K., Davis, T. M., et al. 2011, MNRAS, 418, 1725 [NASA ADS] [CrossRef] [Google Scholar]
 Blake, C., Joudaki, S., Heymans, C., et al. 2016, MNRAS, 456, 2806 [NASA ADS] [CrossRef] [Google Scholar]
 Bolton, A. S., Schlegel, D. J., Aubourg, É., et al. 2012, AJ, 144, 144 [NASA ADS] [CrossRef] [Google Scholar]
 Bundy, K., Leauthaud, A., Saito, S., et al. 2015, ApJS, 221, 15 [NASA ADS] [CrossRef] [Google Scholar]
 Carbone, C., Petkova, M., & Dolag, K. 2016, JCAP, 7, 034 [CrossRef] [Google Scholar]
 Chuang, C.H., & Wang, Y. 2012, MNRAS, 426, 226 [NASA ADS] [CrossRef] [Google Scholar]
 Chuang, C.H., Prada, F., Cuesta, A. J., et al. 2013, MNRAS, 433, 3559 [NASA ADS] [CrossRef] [Google Scholar]
 Coil, A. L., Blanton, M. R., Burles, S. M., et al. 2011, ApJ, 741, 8 [NASA ADS] [CrossRef] [Google Scholar]
 Cuesta, A. J., VargasMagaña, M., Beutler, F., et al. 2016, MNRAS, 457, 1770 [NASA ADS] [CrossRef] [Google Scholar]
 Damour, T. 2000, Nucl. Phys. B Proc. Suppl., 80, 41 [NASA ADS] [CrossRef] [Google Scholar]
 Dark Energy Survey Collaboration 2005, ArXiv eprints [arXiv:astroph/0510346] [Google Scholar]
 Dawson, K. S., Schlegel, D. J., Ahn, C. P., et al. 2013, AJ, 145, 10 [Google Scholar]
 de la Torre, S., & Guzzo, L. 2012, MNRAS, 427, 327 [NASA ADS] [CrossRef] [Google Scholar]
 de la Torre, S., Guzzo, L., Peacock, J. A., et al. 2013, A&A, 557, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 de la Torre, S., Jullo, E., Giocoli, C., et al. 2017, A&A, 608, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 DES Collaboration (Abbott, T. M. C., et al.) 2018, Phys. Rev. D, 98, 043526 [Google Scholar]
 DESI Collaboration (Aghamousa, A., et al.) 2016a, ArXiv eprints [arXiv:1611.00036] [Google Scholar]
 DESI Collaboration (Aghamousa, A., et al.) 2016b, ArXiv eprints [arXiv:1611.00037] [Google Scholar]
 Dvali, G., Gabadadze, G., & Porrati, M. 2000, Phys. Lett. B, 485, 208 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Dvornik, A., Hoekstra, H., Kuijken, K., et al. 2018, MNRAS, 479, 1240 [NASA ADS] [CrossRef] [Google Scholar]
 Eisenstein, D. J., Weinberg, D. H., Agol, E., et al. 2011, AJ, 142, 72 [Google Scholar]
 Erben, T., Hildebrandt, H., Miller, L., et al. 2013, MNRAS, 433, 2545 [NASA ADS] [CrossRef] [Google Scholar]
 Escoffier, S., Cousinou, M. C., Tilquin, A., et al. 2016, ArXiv eprints [arXiv:1606.00233] [Google Scholar]
 Euclid Collaboration (Martinet, N., et al.) 2019, A&A, 627, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Feldman, H. A., Kaiser, N., & Peacock, J. A. 1994, ApJ, 426, 23 [NASA ADS] [CrossRef] [Google Scholar]
 Ferreira, P. G., & Skordis, C. 2010, Phys. Rev. D, 81, 104020 [NASA ADS] [CrossRef] [Google Scholar]
 Fisher, K. B., Davis, M., Strauss, M. A., Yahil, A., & Huchra, J. P. 1994, MNRAS, 267, 927 [NASA ADS] [CrossRef] [Google Scholar]
 Garilli, B., Le Fèvre, O., Guzzo, L., et al. 2008, A&A, 486, 683 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gaztañaga, E., Eriksen, M., Crocce, M., et al. 2012, MNRAS, 422, 2904 [NASA ADS] [CrossRef] [Google Scholar]
 GilMarín, H., Wagner, C., Noreña, J., Verde, L., & Percival, W. 2014, JCAP, 12, 029 [NASA ADS] [CrossRef] [Google Scholar]
 GilMarín, H., Percival, W. J., Verde, L., et al. 2017, MNRAS, 465, 1757 [NASA ADS] [CrossRef] [Google Scholar]
 Giocoli, C., Jullo, E., Metcalf, R. B., et al. 2016, MNRAS, 461, 209 [NASA ADS] [CrossRef] [Google Scholar]
 Gunn, J. E., Siegmund, W. A., Mannery, E. J., et al. 2006, AJ, 131, 2332 [NASA ADS] [CrossRef] [Google Scholar]
 Guzzo, L., Scodeggio, M., Garilli, B., et al. 2014, A&A, 566, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 HarnoisDéraps, J., Amon, A., Choi, A., et al. 2018, MNRAS, 481, 1337 [NASA ADS] [CrossRef] [Google Scholar]
 Heymans, C., Van Waerbeke, L., Miller, L., et al. 2012, MNRAS, 427, 146 [NASA ADS] [CrossRef] [Google Scholar]
 Hildebrandt, H., Erben, T., Kuijken, K., et al. 2012, MNRAS, 421, 2355 [NASA ADS] [CrossRef] [Google Scholar]
 Hildebrandt, H., Viola, M., Heymans, C., et al. 2017, MNRAS, 465, 1454 [Google Scholar]
 Ho, S., Cuesta, A., Seo, H.J., et al. 2012, ApJ, 761, 14 [NASA ADS] [CrossRef] [Google Scholar]
 Joudaki, S., Mead, A., Blake, C., et al. 2017, MNRAS, 471, 1259 [NASA ADS] [CrossRef] [Google Scholar]
 Klypin, A., Yepes, G., Gottlöber, S., Prada, F., & Heß, S. 2016, MNRAS, 457, 4340 [NASA ADS] [CrossRef] [Google Scholar]
 Laureijs, R., Amiaux, J., Arduini, S., et al. 2011, ArXiv eprints [arXiv:1110.3193] [Google Scholar]
 Leauthaud, A., Saito, S., Hilbert, S., et al. 2017, MNRAS, 467, 3024 [NASA ADS] [CrossRef] [Google Scholar]
 Leonard, C. D., Ferreira, P. G., & Heymans, C. 2015, JCAP, 2015, 051 [Google Scholar]
 Lesgourgues, J. 2011, ArXiv eprints [arXiv:1104.2932] [Google Scholar]
 Lima, M., Cunha, C. E., Oyaizu, H., et al. 2008, MNRAS, 390, 118 [NASA ADS] [CrossRef] [Google Scholar]
 LSST Dark Energy Science Collaboration 2012, ArXiv eprints [arXiv:1211.0310] [Google Scholar]
 Magaña, J., Motta, V., Cárdenas, V. H., Verdugo, T., & Jullo, E. 2015, ApJ, 813, 69 [NASA ADS] [CrossRef] [Google Scholar]
 Mandelbaum, R., Hirata, C. M., Broderick, T., Seljak, U., & Brinkmann, J. 2006, MNRAS, 370, 1008 [NASA ADS] [CrossRef] [Google Scholar]
 Mandelbaum, R., Seljak, U., Hirata, C. M., et al. 2008, MNRAS, 386, 781 [NASA ADS] [CrossRef] [Google Scholar]
 Maraston, C., Pforr, J., Henriques, B. M., et al. 2013, MNRAS, 435, 2764 [NASA ADS] [CrossRef] [Google Scholar]
 Marta Pinho, A., Casas, S., & Amendola, L. 2018, JCAP, 11, 027 [NASA ADS] [CrossRef] [Google Scholar]
 Martinet, N., Schneider, P., Hildebrandt, H., et al. 2018, MNRAS, 474, 712 [NASA ADS] [CrossRef] [Google Scholar]
 Matsubara, T., & Suto, Y. 1996, ApJ, 470, L1 [NASA ADS] [CrossRef] [Google Scholar]
 McDonald, P., & Roy, A. 2009, JCAP, 8, 020 [NASA ADS] [CrossRef] [Google Scholar]
 Metcalf, R. B., & Petkova, M. 2014, MNRAS, 445, 1942 [NASA ADS] [CrossRef] [Google Scholar]
 Moraes, B., Kneib, J.P., Leauthaud, A., et al. 2014, Rev. Mex. Astron. Astrofis. Conf. Ser., 44, 202 [NASA ADS] [Google Scholar]
 Nakajima, R., Mandelbaum, R., Seljak, U., et al. 2012, MNRAS, 420, 3240 [NASA ADS] [Google Scholar]
 Newman, J. A., Cooper, M. C., Davis, M., et al. 2013, ApJS, 208, 5 [Google Scholar]
 Nuza, S. E., Sánchez, A. G., Prada, F., et al. 2013, MNRAS, 432, 743 [NASA ADS] [CrossRef] [Google Scholar]
 Paz, D. J., & Sánchez, A. G. 2015, MNRAS, 454, 4326 [NASA ADS] [CrossRef] [Google Scholar]
 Pezzotta, A., de la Torre, S., Bel, J., et al. 2017, A&A, 604, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VI. 2018, ArXiv eprints [arXiv:1807.06209] [Google Scholar]
 Pullen, A. R., Alam, S., He, S., & Ho, S. 2016, MNRAS, 460, 4098 [NASA ADS] [CrossRef] [Google Scholar]
 Reid, B., Ho, S., Padmanabhan, N., et al. 2016, MNRAS, 455, 1553 [NASA ADS] [CrossRef] [Google Scholar]
 Reyes, R., Mandelbaum, R., Seljak, U., et al. 2010, Nature, 464, 256 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 RodríguezTorres, S. A., Chuang, C.H., Prada, F., et al. 2016, MNRAS, 460, 1173 [NASA ADS] [CrossRef] [Google Scholar]
 Ross, A. J., Percival, W. J., Sánchez, A. G., et al. 2012, MNRAS, 424, 564 [NASA ADS] [CrossRef] [Google Scholar]
 Ross, A. J., Beutler, F., Chuang, C.H., et al. 2017, MNRAS, 464, 1168 [NASA ADS] [CrossRef] [Google Scholar]
 Saito, S., Baldauf, T., Vlah, Z., et al. 2014, Phys. Rev. D, 90, 123522 [NASA ADS] [CrossRef] [Google Scholar]
 Seo, H.J., & Eisenstein, D. J. 2003, ApJ, 598, 720 [NASA ADS] [CrossRef] [Google Scholar]
 Shan, H., Liu, X., Hildebrandt, H., et al. 2018, MNRAS, 474, 1116 [NASA ADS] [CrossRef] [Google Scholar]
 Shirasaki, M., Takada, M., Miyatake, H., et al. 2017, MNRAS, 470, 3476 [NASA ADS] [CrossRef] [Google Scholar]
 Simpson, F., Heymans, C., Parkinson, D., et al. 2013, MNRAS, 429, 2249 [NASA ADS] [CrossRef] [Google Scholar]
 Singh, S., Mandelbaum, R., Seljak, U., Slosar, A., & Vazquez Gonzalez, J. 2017, MNRAS, 471, 3827 [NASA ADS] [CrossRef] [Google Scholar]
 Singh, S., Alam, S., Mandelbaum, R., et al. 2019, MNRAS, 482, 785 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, R. E., Peacock, J. A., Jenkins, A., et al. 2003, MNRAS, 341, 1311 [NASA ADS] [CrossRef] [Google Scholar]
 Sugai, H., Karoji, H., Takato, N., et al. 2012, in Groundbased and Airborne Instrumentation for Astronomy IV, SPIE Conf. Ser., 8446, 84460Y [CrossRef] [Google Scholar]
 Takahashi, R., Sato, M., Nishimichi, T., Taruya, A., & Oguri, M. 2012, ApJ, 761, 152 [NASA ADS] [CrossRef] [Google Scholar]
 Taruya, A., Nishimichi, T., & Saito, S. 2010, Phys. Rev. D, 82, 063522 [NASA ADS] [CrossRef] [Google Scholar]
 Taylor, A., & Joachimi, B. 2014, MNRAS, 442, 2728 [NASA ADS] [CrossRef] [Google Scholar]
 Tully, R. B., Courtois, H. M., & Sorce, J. G. 2016, AJ, 152, 50 [NASA ADS] [CrossRef] [Google Scholar]
 White, M., Blanton, M., Bolton, A., et al. 2011, ApJ, 728, 126 [NASA ADS] [CrossRef] [Google Scholar]
 Xu, X., Cuesta, A. J., Padmanabhan, N., Eisenstein, D. J., & McBride, C. K. 2013, MNRAS, 431, 2834 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, P., Liguori, M., Bean, R., & Dodelson, S. 2007, Phys. Rev. Lett., 99, 141302 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
Appendix A: Weaklensing systematics tests
Masking. In order to assess the impact of missing tiles and largescale masking (e.g. due to very bright stars), we compute the density of CS82 galaxies on a grid with pixel size ∼1 deg. Then, we randomly draw mock galaxies in the field such that the overall redshift distribution and total number of sources matches observations. Finally, we downsample this catalogue according the density fluctuations attributed to masking. We find that masking increases the statistical noise in the GGL measurement by about 20% at all scales. However we could not identify any obvious systematic bias related to masking.
Photometric redshifts bias.Mandelbaum et al. (2008) and Nakajima et al. (2012) proposed an alternative method to estimate the bias introduced by photometric redshifts on GGL measurements. These authors proposed to estimate the bias b_{z}(z_{lens}) between photometric redshifts and spectroscopic redshifts ΔΣ measurements,
The summation is performed over the subset of source galaxies with both spectroscopic and photometric redshifts. We adapted the original expression from Mandelbaum et al. (2008) such that the inverse critical densities converges to zero when the source redshift becomes smaller than z_{lens}. The value w_{j} is the weight on source galaxy j in the lensing catalogue. In order to estimate the effective bias on our GGL measurements with CMASS galaxies, we need to integrate the bias function b_{z}(z_{lens}) over the CMASS redshift distribution p(z_{lens}) such that
where the weight on each lens place is correcting for the fact that the number of sources involved in a given aperture in physical coordinates includes more objects at lower than at higher redshifts. We bootstrapped our catalogues to estimate the uncertainties on our bias estimates.
For this measurement, we used VVDS (i_{AB} < 22.5, Garilli et al. 2008), DEEP2 (R_{AB} < 24.1, Newman et al. 2013), PRIMUS (i_{AB} < 23.5, Coil et al. 2011), VIPERS (i_{AB} < 22.5, Guzzo et al. 2014), and SDSSDR13 spectroscopic redshifts, which we matched to our lensing sources in our four fields. On Stripe 82, we obtained b_{z} = −0.028 ± 0.006, b_{z} = −0.131 ± 0.004 and b_{z} = −0.082 ± 0.004 for BPZ, neural network or LePhare codes, respectively. With CFHTLens, we obtained b_{z} = +0.003 ± 0.003, b_{z} = −0.014 ± 0.004 and b_{z} = +0.022 ± 0.003 on fields W1, W3, and W4, respectively.
We also adapted our code to assess the improvement obtained by using the photometric redshift probability of each source galaxy p(z_{s}) instead of maximumlikelihood values. The critical densities then become
On field Stripe 82, we found that using the full probability p(z_{s}) halved the bias obtained with LePhare code to b_{z} = −0.031 ± 0.005. Nonetheless, using the bestfit redshifts provided by BPZ still yields the smallest bias.
Catastrophic photometric redshifts. In order to assess the impact of catastrophic redshifts on the lensing measurements, we computed the twodimensional probability p(z_{LP}z_{ANNZ}) of obtaining a photometric redshift with Le Phare given a photometric redshift obtained with neural network. Assuming this later to be the true redshift, we degraded the true redshifts in our mocks to reproduce the catastrophic outlier effects. Overall, we found that catastrophic redshifts could bias the lensing signal by about b_{z} = +0.03. This is in agreement with our estimations above with spectroscopic redshifts and the estimates found in Leauthaud et al. (2017).
Asymmetric posterior on. E_{G} It is very typical that observational estimators obtained from a ratio of observables have asymmetric probability distribution function. Indeed using our simulations, we found that E_{G} is systematically lower than ΛCDM−GR predictions, with a long tail towards larger values of E_{G}, as shown in Fig. A.1. When applying the usual estimator on our mocks, we also find mean values smaller than expected, although still in statistical agreement. We note finally that 1σ constraints are tighter with the fit than with the usual E_{G} estimator. Given that Amendola et al. (2013) have demonstrated that E_{G} is probably not as effective as anticipated in some particular cosmological models, and Alam et al. (2017) and Singh et al. (2019) have started to apply bias corrections to this estimator based on mock catalogues built with fiducial cosmological models, we think it might be as efficient and clean to fit the correlation functions, and derive E_{G} by marginalizing over the model parameters, or simply compare Ω_{m0} and f to their ΛCDM+GR predictions.
Fig. A.1.
Recovered signal in the mocks, when (grey shared area with 1σ and 2σ C.L. contours) and (orange 1σ shaded area) are computed from the MCMC samples, and when is directly estimated from the mocks (red data points with 1σ error bars). In this latter case, we take β = 0.84/2.13 = 0.39, as obtained from a previous fit of our model to the RSDonly measurements in the mocks. Both definitions are in agreement with the value of Ω_{m}/f computed using the Planck cosmology 2018 of the simulation (black dashed line). Measurements are performed with R_{0} = 1.0 h^{−1} Mpc, s_{min} = 17.8 h^{−1} Mpc and no tapering. 
Fig. A.2.
Contours at 1, 2, and 3σ C.L. of all our model parameters estimated with the RSD only, GGL only, and their combination. In all cases, we set s_{min} = 17.8 h^{−1} Mpc and R_{0} = 1.0 h^{−1} Mpc (see text for details). 
All Tables
Number of CMASS galaxies per field, effective lensing area after masking and number of WL sources.
Properties of the simulated fields in terms of independent mock catalogue, random resampling of lensing shape noise and photometric redshifts per catalogue, and number of subregions for Jackknife resampling.
Bestfit and derived parameters obtained by fitting the RSD only, GGL only, and their combination.
All Figures
Fig. 1.
Redshift distribution of CMASS galaxies (blue) compared to CS82 (orange) and CFHTLens (green) source distributions, after WL selection has been performed. WL n(z) are based on photometric redshifts (see text for details). 

In the text 
Fig. 2.
Probability distribution of having a photometric redshift with BPZ and a spectroscopic redshift for each field. Contours are given for 1, 2, and 3σ C.L. 

In the text 
Fig. 3.
Left panel: number density of CMASS mock galaxies for our 4 simulated fields. Limits of our analysis are indicated with blue dashed lines. Measurements from Anderson et al. (2012, A12) are in grey. Right panel: CMASS CMF for the 4 fields reproducing the observed incompleteness. The mock catalogue is complete at high mass in agreement with the model proposed in RodríguezTorres et al. (2016) in grey. 

In the text 
Fig. 4.
Comparison of lensing measurements performed on simulated catalogues affected and not affected by photometric redshifts noise and ΔΣ, respectively. The grey shaded areas correspond to the uncertainties on the mean value obtained by resampling the multiple noises in different light cones. 

In the text 
Fig. 5.
Bias between recovered parameters f and Ω_{m} relative to the values used in the mocks as functions of the minimum scale of the multipoles s_{min}, and the cutoff radius of the Υ(R, R_{0}) lensing estimator. Values of s_{min} = 17.8 h^{−1} Mpc and R_{0} = 1.0 h^{−1} Mpc provide unbiased estimates of f and Ω_{m}. All measurements were performed without tapering smoothing in the covariance matrices. 

In the text 
Fig. 6.
Variation of the relative errors on the parameters f and Ω_{m}, as a function of the smooth scale T_{p} in the covariance matrices. There is no improvement below T_{p} = 12 h^{−1}. 

In the text 
Fig. 7.
Matrices of covariance (top panel) and precision (bottom panel) in logarithmic color scale for the 4 fields used in this analysis. In both panels, the upper triangular part of the matrices represents the case without tapering, while the lower part represents the case with tapering T_{p} = 12 h^{−1} Mpc. Noise between far apart scales is significantly decreased and the errors on the model parameters converge to a minimum. 

In the text 
Fig. 8.
Monopole (red) and quadrupole (blue) measurements with mock catalogues (shaded region), real data (solid lines), and theoretical predictions with a linear bias parameter b_{1} = 1.8 (dashed lines). Black dots represent prereconstruction measurements with the full DR12v5 CMASS sample from Cuesta et al. (2016). 

In the text 
Fig. 9.
Filtered Υ_{gm} and nonfiltered ΔΣ GGL measurements with mocks (shaded regions), ΔΣ and Υ data (blue and orange points respectively), and theory with a linear bias parameter b_{1} = 1.8 (dashed line). Black dots in S82 panel represent ΔΣ measurements from L16, and Υ_{gm} measurements from Alam et al. (2017) in CFHTLS panels. 

In the text 
Fig. 10.
Improvement on estimating Ω_{m} and σ_{8} when combining RSD and WL measurements. The blue and black curves are, respectively, obtained with WL and RSD constraints only. The orange curves are obtained with the combination of WL and RSD. Contours are given at 1, 2, and 3σ C.L. Vertical lines indicate Planck TT,TE,EE+lowE 2018 results. 

In the text 
Fig. 11.
Growth rate f as a function of redshift compared to recent measurements. The black line and surrounding grey shared area indicate the Planck TT,TE,EE+lowE 2018 mean and 1σ uncertainty predictions for ΛCDM−GR flat model. 

In the text 
Fig. 12.
Measurement of E_{G} with combined constraints in the fields CFHTStripe 82 and CFHTLs W1, W3, and W4. The horizontal black line indicates the Planck TT,TE,EE+lowE 2018 prediction. We note that CFHTS82 data help shrink the error bars by about 30%. CFHTLens points have been shifted rightwards for clarity. 

In the text 
Fig. 13.
Measurements of E_{G} at various redshifts. Effective redshifts of the measurements have been slightly shifted for clarity. 

In the text 
Fig. A.1.
Recovered signal in the mocks, when (grey shared area with 1σ and 2σ C.L. contours) and (orange 1σ shaded area) are computed from the MCMC samples, and when is directly estimated from the mocks (red data points with 1σ error bars). In this latter case, we take β = 0.84/2.13 = 0.39, as obtained from a previous fit of our model to the RSDonly measurements in the mocks. Both definitions are in agreement with the value of Ω_{m}/f computed using the Planck cosmology 2018 of the simulation (black dashed line). Measurements are performed with R_{0} = 1.0 h^{−1} Mpc, s_{min} = 17.8 h^{−1} Mpc and no tapering. 

In the text 
Fig. A.2.
Contours at 1, 2, and 3σ C.L. of all our model parameters estimated with the RSD only, GGL only, and their combination. In all cases, we set s_{min} = 17.8 h^{−1} Mpc and R_{0} = 1.0 h^{−1} Mpc (see text for details). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.