Darkenergy constraints and correlations with systematics from CFHTLS weak lensing, SNLS supernovae Ia and WMAP5^{}
M. Kilbinger^{1}  K. Benabed^{1}  J. Guy^{2}  P. Astier^{2}  I. Tereno^{1,3}  L. Fu^{4,1,5}  D. Wraith^{6,1}  J. Coupon^{1}  Y. Mellier^{1}  C. Balland^{2}  F. R. Bouchet^{1}  T. Hamana^{7}  D. Hardin^{2}  H. J. McCracken^{1}  R. Pain^{2}  N. Regnault^{2}  M. Schultheis^{8}  H. Yahagi^{9}
1  Institut d'Astrophysique de Paris, CNRS UMR 7095 & UPMC, 98bis
boulevard Arago, 75014 Paris, France
2 
LPNHE, CNRSIN2P3 and Universités Paris VI & VII,
4 place Jussieu, 75252 Paris Cedex 05, France
3 
ArgelanderInstitut für Astronomie, Universität Bonn, Auf dem
Hügel 71, 53121 Bonn, Germany
4 
INAF  Osservatorio Astronomico di Capodimonte, via Moiariello
16, 80131 Napoli, Italy
5 
Shanghai Key Lab for Astrophysics, Shanghai Normal University,
Shanghai 200234, PR China
6 
CEREMADE, Université Paris Dauphine, 75775 Paris Cedex 16, France
7 
National Astronomical Observatory of Japan, Mitaka, Tokyo
1818588, Japan
8 
Observatoire de Besançon, 41bis avenue de l'Observatoire, BP
1615, 25010 Besançon Cedex, France
9 
Research Institute for Information Technology, University of Kyushu
6101 Hakozaki, Higashiward, Fukuoka 8128581, Japan
Received 28 October 2008 / Accepted 18 January 2009
Abstract
Aims. We combine measurements of weak gravitational lensing from the CFHTLSWide survey, supernovae Ia from CFHT SNLS and CMB anisotropies from WMAP5 to obtain joint constraints on cosmological parameters, in particular, the darkenergy equationofstate parameter w. We assess the influence of systematics in the data on the results and look for possible correlations with cosmological parameters.
Methods. We implemented an MCMC algorithm to sample the parameter space of a flat CDM model with a darkenergy component of constant w. Systematics in the data are parametrised and included in the analysis. We determine the influence of photometric calibration of SNIa data on cosmological results by calculating the response of the distance modulus to photometric zeropoint variations. The weak lensing data set is tested for anomalous fieldtofield variations and a systematic shape measurement bias for highredshift galaxies.
Results. Ignoring photometric uncertainties for SNLS biases cosmological parameters by at most 20% of the statistical errors, using supernovae alone; the parameter uncertainties are underestimated by 10%. The weaklensing fieldtofield variance between 1 deg^{2}MegaCam pointings is 515% higher than predicted from Nbody simulations. We find no bias in the lensing signal at high redshift, within the framework of a simple model, and marginalising over cosmological parameters. Assuming a systematic underestimation of the lensing signal, the normalisation
increases by up to 8%. Combining all three probes we obtain
0.10 < 1+w < 0.06 at 68% confidence (
0.18<1+w<0.12 at 95%), including systematic errors. Our results are therefore consistent with the cosmological constant .
Systematics in the data increase the error bars by up to 35%; the bestfit values change by less than 0.15.
Key words: cosmology: observations  cosmology: cosmological parameters  methods: statistical
1 Introduction
The CanadaFranceHawaiiTelescope Legacy Survey (CFHTLS^{}) has yielded interesting constraints on cosmological parameters using different probes, for example weak gravitational lensing and supernovae of type Ia. In this paper we combine two recent measurements from the CFHTLS: the cosmic shear data (Fu et al. 2008, hereafter F08) and supernova (SN) Ia data from the Supernova Legacy Survey (SNLS, Astier et al. 2006, hereafter A06). We complement these with CMB measurements from the Wilkinson Microwave Anisotropy Probe (WMAP5, Hinshaw et al. 2009, hereafter H09).
There are hundreds of typeIa supernovae for which highquality observations are available. Many surveys spanning wide redshift ranges and using different telescopes and search strategies can be accessed. A recent compilation of various surveys and the derived cosmological constraints can be found in Kowalski et al. (2008). Despite this impressive number of available objects we have chosen to restrict ourselves to the SNLS sample. The resulting constraints will suffer from a greater statistical uncertainty; however, the systematic errors are better understood by using a single, homogeneous survey. By providing imaging with a single telescope and camera for example, SNLS allows a common photometric calibration strategy for the whole survey.
The lensing data presented in F08 might suffer from systematic errors in the shear calibration and shape measurement. In this work, we performed two tests for the presence of systematics. First, we compared the variance between individual MegaCam pointings to simulations. Second, we quantified a potential systematic underestimation of the lensing signal at high redshifts. The influence on cosmological parameters was estimated. Due to the unknown origin of those systematics, the analysis must remain less rigorous than for the SNIa case. It is by no means comprehensive since we only quantify the effect of those systematics and do not investigate their origin, which is beyond the scope of this paper.
The analysis of the fiveyear WMAP data has yielded impressive results on many cosmological parameters from a great number of models (Dunkley et al. 2009; Komatsu et al. 2009; Hinshaw et al. 2009). Thanks to the rich features in the angular power spectrum, many cosmological parameters can be determined with high precision. However, degeneracies between parameters remain, in particular, for darkenergy models. In order to lift these degeneracies, measurements of CMB anisotropies have to be complemented with other probes.
The WMAP5 team has combined their data with other probes which are sensitive only to the geometry of the Universe, i.e. SNIa and baryonic acoustic oscillations (BAO). Other teams have included probes of structure growth, like counts of Xray clusters (Mantz et al. 2008), SDSS and Lyforest (Xia et al. 2008). In this work, we combine WMAP5 with weak gravitational lensing which is sensitive to both the geometry and the growth of structure. Probing both domains will allow future surveys to distinguish between dark energy and modified gravity as a possible cause for the present acceleration of the Universe (e.g. Jain & Zhang 2008). Such a test will be feasible with upcoming and proposed surveys such as KIDS^{}, DES^{}, LSST^{}, JDEM^{} and Euclid^{}.
The firstyear WMAP data has already been combined with weak lensing, using the RCS survey (Hoekstra et al. 2002), to improve constraints on and (Contaldi et al. 2003). In F08, cosmic shear has been supplemented with the thirdyear WMAP data. In this work we extend the latter, simple analysis by dropping the (not well motivated) priors on some parameters for lensing ( ).
The ``concordance'' flat CDM model of cosmology provides an excellent fit to WMAP5 and most other probes of the geometry and largescale structure of the Universe. It only contains six free parameters, , , , , h and (or functions thereof). This model assumes a cosmological constant as the cause of the observed accelerated expansion of the Universe today. The cosmological constant has yet evaded all plausible physical explanations of its nature and origin and, further, brings in problems of finetuning and coincidence. A signature of some ``dark energy'' beyond a simple cosmological constant might be an equationofstate (eos) parameter which deviates from the vacuum energy value of 1. We therefore extend the concordance sixparameter model by including the darkenergy eos parameter w. Although this wCDM model with a constant w is not better motivated physically than the cosmological constant, a significant observed deviation from w=1 will definitely be an indication for new physics. Moreover, the data at present are not good enough to constrain more dark energy parameters in a general way. This is only feasible for very specific models, e.g. quintessence (see Schimd et al. 2007) or early dark energy (Bartelmann et al. 2006; Francis et al. 2008), and will be subject of a future work using the CFHTLS data (Tereno et al., in prep.).
One goal of this paper is to focus on systematic errors, nuisance parameters and their interplay with cosmological parameters. For the SNIa data, apart from the usual lightcurve parameters, we take into account photometric calibration errors. A joint Bayesian analysis including systematic and cosmological parameters is done and correlations are revealed. This will be mandatory for future surveys with statistical errors which will be smaller by several orders of magnitude as compared to today. To decrease systematics further and further is very challenging, and technical limitations might set a barrier to this endeavour. Therefore, it is important to quantify the effect of systematics and nuisance factors on cosmological constraints.
This paper is organised as follows. Section 2 contains a brief description of the data together with the likelihood used later in the analysis. It also discusses sources of systematics in the data and their correlations with cosmological parameters. In Sect. 3 we define the cosmological model tested in this work, introduce our implementation of the Monte Carlo Markov Chain (MCMC) technique and present the cosmological results of the analysis. We conclude with a discussion in Sect. 4 and an outlook in Sect. 5.
2 Data and method
2.1 CFHTLSWide cosmic shear
We use the cosmic shear results from the CFHTLSWide year data release (T0003, Fu et al. 2008). On 57 square degrees (35 sq. deg. effective area), about 2 10^{6} galaxies with i_{AB}magnitudes between 21.5 and 24.5 were imaged. The data, reduction analyses and shear pipeline are described in detail in F08. We use the aperturemass dispersion (Schneider et al. 1998) measured between 2 and 230 arcmin. Due to the compensated nature of its filter, this secondorder measure is least sensitive to largescale systmatics in the data. The source redshift distribution is obtained by using the CFHTLSDeep p(z) (Ilbert et al. 2006) and by rescaling it according to the Wide i_{AB} magnitue distribution and weaklensing galaxy weights.
To ascertain that the quality and reliability of the shear measurements are sufficient for this work, we perform further tests of the data beyond what has been done in F08. This is addressed in Sects. 2.1.3 and 2.1.2 where we assess the importance of potential systematics for the current data, and estimate their influence on inferred cosmological parameters.
As in F08 the loglikelihood is
,
where the
is modelled as
The predicted aperturemass dispersion given by a model parameter vector is fitted to measured at angular scales . The covariance C of is the one used in F08 and contains shape noise, (nonGaussian) cosmic variance and the residual Bmode. Those parts of the covariance which depend on the shear correlation (mixed and cosmic variance terms) are calculated using a theoretical model of the largescale structure and therefore depend on cosmological parameters. We ignore this dependence and keep the covariance constant, corresponding to the fiducial cosmology as in F08. This biases the posterior confidence regions, but the effect is weak over the region of parameter space permitted by CMB and lensing, see Eifler et al. (2008) for a detailed discussion. Moreover, we can drop the term in the loglikelihood because its parameterdependence manifests itself only for very small survey areas (Kilbinger & Munshi 2006).
As usual, the following relation between the aperturemass dispersion and the weak lensing power spectrum holds (Schneider et al. 1998),
(2) 
The lensing power spectrum is a projection of the 3d matterdensity power spectrum , weighted by the redshift distribution (Kaiser 1992),
(3) 
For the nonlinear evolution of the power spectrum, the fitting formula of Smith et al. (2003) is used. Although tested for CDM models, it provides reasonable good fits to wCDM cosmologies as well (Ma 2007). The accuracy of any nonlinear fitting function is limited; when using the ansatz of Peacock & Dodds (1996) instead of Smith et al. (2003), the resulting bestfit differs by 2% (F08).
We parametrise the redshift distribution using the function
which we fit to the obtained redshift histogram. The corresponding is
Here, n_{i} is the normalised number of galaxies in the redshift bin and p(z_{i}) the fitting function, evaluated at the redshift bin centre. The uncertainty of n_{i} contains Poisson noise, photoz uncertainty and cosmic variance, as described in F08 and Benjamin et al. (2007). The sum in Eq. (5) extends over the redshift range . In this range, cosmic variance is the dominant uncertainty. As in F08 we neglect the crosscorrelation between different bins.
2.1.1 Systematics in the lensing data
In the following two subsections we address the question of potential residual systematics in the cosmic shear data. We estimate the influence of some systematics on cosmological constraints. The shear catalogue used here has been extensively tested in F08. Internal consistency checks have been performed involving the comparison between two data reduction pipelines, residual Bmodes and crosscorrelations between galaxies and stars. In addition, the shape measurement pipeline has been calibrated with the STEP1 and STEP2 simulations (Heymans et al. 2006; Massey et al. 2007), and a bias of less than 3% has been determined.
Despite that, there are indications of remaining systematics which do not manifest themselves in the merged shear catalogue. For example, there are variations of the shear signal between individual MegaCam pointings which vanish on average. Moreover, with the help of photometric redshifts covering parts of the survey, a problem with the redshiftscaling of the shear correlation becomes apparent.
2.1.2 Variance between CFHTLS pointings
In F08 it was shown the systematics are small globally. The Bmode is consistent with zero on most scales (note however a significant detection around 100 arcmin). The level of systematics is low on average, e.g. the residual crosscorrelation between uncorrected star ellipticities and galaxy shapes after PSF correction. There might however be systematics present on individual fields. Even if they vanish on average they are an additional noise source and increase the measurement error bars.
Here, we address the question whether individual pointings show anomalous variations which might be due to incorrect shape measurements or PSF correction. The fieldtofield variance of the lensing signal is compared with numerical simulations.
We measure the shear aperturemass dispersion on all 57 fields, each one corresponding to a 1 deg^{2}MegaCam image, and calculate the fieldtofield variance. We use Nbody simulations (Yahagi & Yoshii 2001; Yahagi 2005) through which we shoot light rays to obtain shear maps (Hamana & Mellier 2001) using a redshift distribution corresponding to the mean n(z) estimated for the observations. The underlying model is a flat WMAP3cosmology ( ).
An important issue are the correlations between CFHTLSWide pointings due to largescale structure. Unfortunately, we cannot simulate the exact survey geometry of the CFHTLSWide T0003 survey which extends up to 8 degrees angular separation, since our simulations are only in size. We therefore use several (independent) raytracing simulations to cover the observed fields; three for W1, one for W2 and four for W3. This results in lesscorrelated fields and a smaller cosmic variance compared to the observations. We therefore expect the simulations to slightly underestimate the observed fieldtofield variance. The galaxies on each simulated pointing are distributed homogeneously, we do not take into account masking. We do however simulate the varying number density between fields by applying the observed number densities to the simulations.
From Fig. 1 we see that the observed Emode variance of is 1020% higher than the predicted value, on scales between 1 and 30 arcmin. The Bmode variance is in very good agreement, with a tendency to be slightly lower than predicted. From the fact that the Bmode variance is of the same order as the Emode, one sees that the considered angular scales are shapenoise and not cosmicvariance dominated. The powerlaw shape of the variance is another indication for this. The lowlevel oscillations on scales on the order of the correlation length are probably due to noise.
The observed tophat fieldtofield variance is substantially greater than the one from simulations, in particular at large angular scales. This might be due to its higher sensitivity to residuals on large scales compared to the aperturemass. Note however that because of the high correlations between different angular scales for tophat, the significance of any discrepancy is hard to interpret. Due to the small number of simulations we cannot attempt an error estimate on the fieldtofield variance.
Note that the variance for the merged catalogue used in F08, and this work for the cosmological constraints, is roughly a factor of four less than the one from individual fields.
The greater observed fieldtofield fluctuations could be a sign for uncorrected residual systematics. We have not included in the simulations the varying mean redshift due to the different numbers of exposures for each field in the final stacks. The variance in limiting magnitude between fields is 0.25. With the empirical law between limiting magnitude and mean redshift from van Waerbeke et al. (2006) this translates into a zvariance of 0.03. Using the approximate relation from linear theory we find an additional expected variance of about 4.5%. We conclude that the observed Emode fieldtofield variance is higher than predicted by at most 515%. It is difficult to assess the influence of this additional error on the complete galaxy catalogue used here. First, we do not know how scales larger than 30 arcmin are affected. Second, the shear correlation used to constrain cosmology is calculated with many more galaxy pairs than in this fieldtofield analysis, with a large number of pairs stemming from different MegaCam pointings.
Figure 1: Variance of from the numerical simulations (solid lines) and CFHTLS data (dashed). The rms used for the merged catalogue (from Fu et al. 2008), which we use for the cosmological constraints, is the dotted line. 

Open with DEXTER 
2.1.3 Systematics in the redshiftscaling of the shear signal
Preliminary analysis of weak lensing tomography of the CFHTLSWide shows that beyond a redshift of about unity, the cosmic shear signal does not increase as expected, but instead is systematically underestimated (Fu 2008). For this analysis we use the fourth CFHTLS data release (T0004) which provides fiveband photometry on 35 square degrees. Photometric redshifts for each galaxy have been obtained using the templatefitting code Le Phare (Coupon et al. 2008). The quality of the new photoz's is compatible to those used in F08 (T0003 data release, Ilbert et al. 2006). An improvment has been obtained due to additional and larger spectroscopic samples. We take a subset of T0004 consisting of a contiguous area of 19 square degrees. The measured galaxy shapes are the same as in F08.
As an illustration, we show the tophat shear variance measured at 5 and 25 arcmin for various redshift bins, see Table 1. The measured values, corresponding to 19 square degrees, are compared to predictions using a flat CDM model with and . For the highest bin with redshifts above 0.95, the measured values are clearly inconsistent with the predictions.
Table 1: Shear tophat variance at 5' and 25' and for five redshift ranges. The error bars are Poisson noise.
We have to assume that this anomalous redshiftscaling is also present in the shear catalogue used in this paper for the cosmology constraints. The redshifts for this preliminary analysis are taken from the CFHTLS Wide and are therefore slightly less accurate than the T0003Deep ones, used to infer the F08 redshift distribution. However, the former do not suffer from cosmic variance. The calibration with spectroscopic redshifts assures their reliability for . The depth of the F08 shear galaxies is only half a magnitude higher, and those faint galaxies are downweighted in the weak lensing analysis. We therefore assume that the problem is mainly due to the shape measurement and not to the photometric redshifts. One reason for this could be a shear calibration bias which depends on galaxy properties that are a function of redshift. A bias in the shape measurement which is a function of galaxy size, magnitude or galaxy signaltonoise will affect low and highz galaxies differently.
For this work we choose a simple toy model to parametrise this potential underestimation of the redshiftscaling of the shear signal, see Lesgourgues et al. (2007) for a similar approach. This decrease of the lensing signal corresponds to an effective lensing efficiency which is lower than expected at high redshifts. To model this we multiply the redshift distribution n(z) with a constant c_{0}>0, for redshifts z>z_{0}. We chose z_{0} to be 1.0. c_{0}=1 is the unbiased case with no degradation of the shear signal. c_{0}=0 is the pessimistic case where there is no shear signal for z>z_{0}. We allow values greater than unity corresponding to an overestimation of the shear signal. Note that we use the original, unaltered redshift distribution when fitting the n(z) histogram (Eq. (5)).
Using a model with fixed and and marginalising over and redshift parameters we measure a value of c_{0} consistent with unity but with large error margins, c_{0} = 1.1 (see also Fig. 2, right panel). At 95% confidence nearly the whole range of c_{0} is permitted. Ignoring the highz calibration error leads to smaller error bars for other parameters, in particular (left panel of Fig. 2). For a fixed we respectively get and 0.05 with and without modelling the systematic error. Because we allow the calibration to be smaller and greater than unity there is no significant bias on . If we restrict c_{0} to values smaller or equal 1, thus using the information about the underestimation of the shear signal as prior, we obtain a bestfit of 0.83 ( ). Thus, might be underestimated by which is about one sigma, provided our simple model is correct. From Fig. 3 (left panel) we see the potential shear measurement systematics shifts the likelihood and introduces a tail for high .
We repeat the above analysis with z_{0} = 0.8 and 1.2 and find only small shifts for the powerspectrum normalisation. Writing , we find if c_{0} is varied between 0 and 2, and for c_{0}<1, corresponding in both cases to a 2% change for in the considered range for z_{0}.
Interestingly, we see no correlation between c_{0} and any other individual parameter when the whole posterior is considered. However, c_{0} is of course correlated with the combination of and which determines the shear amplitude. When fixing we see a strong correlation with , as expected (Fig. 3, right panel).
Figure 2: Left panel: 2d marginals (68% and 95%) for and for the two cases of including the highz calibration bias (solid lines) and ignoring it (dashed curves), respectively. Right panel: 1d marginal likelihood for c_{0}. The vertical lines indicate mean (dashed) and 68% confidence intervals (dotted). 

Open with DEXTER 
Figure 3: Left panel: 1d marginal for for , the vertical lines are as in Fig. 2. Right panel: correlation between c_{0} and (68% and 95% confidence levels) for . The vertical lines indicate mean (dashed) and 68% confidence intervals (dotted). 

Open with DEXTER 
To summarise the results on weak lensing systematics, we note that this study is far from complete. Other tests focusing on the PSF correction have been made in F08, and extensive studies going into much more detail as presented here will be published soon (van Waerbeke et al., in prep.).
2.2 SNLS supernovae Ia
The supernova data set corresponds to the first data release of the SNLS, see A06. We use 105 supernovae in total, 71 of which were observed with the CFHTLS. Their redshift range is between 0.25 and 1. The remaining 44 objects are nearby supernovae, with 0.015<z<0.13. In Sect. 2.2.1 we describe the ingredients of the standard likelihood analysis, ignoring systematics, which is similar to A06 with slight modifications. Section 2.2.2 introduces the photometric calibration parameters and defines the updated likelihood function including systematics.
2.2.1 Standard likelihood analysis
We use the results from the SNIa lightcurve fits of A06. For each supernova the fit yields the restframe Bband magnitude m_{B}^{*}, the shape or stretch parameter s, and the colour c. These
quantities are combined in the estimator of the distance modulus for the object as follows:
The universal absolute SNIa magnitude is ; and are universal linear response factors to stretch and colour, respectively.
The error needed for the likelihood includes: (1) the measured uncertainties of the lightcurve parameters m_{B, i}^{*}, s_{i}, c_{i} and their correlations. The corresponding error contribution to can be written as with the parameter vector and the covariance W_{2} of the data vector (m_{B,i}^{*}, s_{i}, c_{i}). (2) An additional redshift uncertainty due to peculiar velocities corresponding to km s^{1}. (3) An intrinsic scatter in absolute magnitude of .
Those three errors are added in quadrature to the loglikelihood expression, which is
The linear dependence on the Hubble constant of the luminosity distance is taken out as an additional summand in Eq. (7) and integrated into the absolute luminosity. From now on we will use the parameter .
The error term from the lightcurve fit, , depends on the stretch and colour responsitivies and . A06 kept those parameters fixed during the minimisation and updated them iteratively in consecutive minimisation runs. This leads to a bias in those parameters of 715%. Instead, we include the dependence of and in the denominator in Eq. (7) to obtain unbiased bestfit results for those parameters which are therefore larger than those cited in A06, see Table 4.
We do not take into account the term involving the covariance determinant in the loglikelihood. Even though in this case this term is not a constant because of the parameterdependent variance, we verified that the effect on cosmological parameters is very small. Moreover, it is in accordance to most other SNIa analyses (e.g. Kowalski et al. 2008; Astier et al. 2006).
2.2.2 Extended analysis with systematics
We extend the analysis of A06 by including the response of the distance modulus to a photometric zeropoint shift in each of the seven filters
(g, r, i, z, U, B, V) and in the Vega (BR) colour in the Landolt system. The reference photometric zeropoints are those used in A06, which were estimated without using the Hubble diagram. If
denotes the vector of those eight zeropoint magnitude shifts, the linear response
is the change in distance
modulus for the supernova for a small change of the
zeropoint,
By changing the zeropoints and redoing the lightcurve fit, we obtain the values of for each supernova. We include this change to the distance modulus linearly in the model, which results in an additional term in the likelihood.
We assume that the zeropoint magnitude shift parameters are uncorrelated variables, since the data in different optical bands have been reduced independently. The noise is the result from a number of different source and can therefore be well approximated to be Gaussian. By definition they mean of is zero. The rms is taken to be for all filters except for z ( ) where we assume the rms to be . The numerical values of those uncertainties are taken from A06 (Sect. 4.1) and correspond to the limits of reproducibility of the photometric calibration. The corresponding (diagonal) covariance matrix is W_{1}.
This prior information is multiplied to the SNIalikelihood in the form of a multivariate Gaussian
likelihood. With the additional term, the loglikelihood expression corresponding to the extended analysis is then
Now, the parameter vector contains the zeropoint parameters . The corresponding optical bands are indicated in Table 2.
Table 2: List of SNIa systematic parameters and their symbols.
2.2.3 Systematics for SNIa
Unlike systematics in weak lensing shape measurements which are difficult to model, the observationrelated systematics for SN are more easily parametrised. In the next section we take into account errors in the estimated distance modulus due to uncertainties in the photometric calibration. The results of the MCMC analysis for the SNIa internal parameters are given in Table 4.
2.2.4 Bias due to systematics
We compare the case of ignoring the systematic errors, using Eq. (7) as the loglikelihood, with the case of fully taking into account the systematics according to the loglikelihood (9). As can be seen in Fig. 4, ignoring the zeropoint errors leads to an asymmetric decrease of the error bars. The constraints get tighter, mainly along the direction of constant luminosity distance, which is the parameterdegeneracy direction. The error bars decrease by about 10%, see Table 3. The bias on parameter means is small, between 10% and 20% of the statistical uncertainty. The bias on the intrinsic SNIa parameters ( ) is even smaller, not more than a few percent of the statistical uncertainty. For a fixed , the absolute biases on both the eos parameter (for ) and the darkenergy density (for w=1) are smaller than for the marginalised case, but remain to be about a tenth of the statistical error.
Figure 4: Confidence contours (68%, 95.5%, 99.7%) for full treatment of errors due to photometry zeropoints (blue solid lines) and ignoring those errors (green dashed curves). These two cases correspond to the two cases in Table 3. The left panel corresponds to w=1, the right panel is for a flat wCDM Universe. The constraints are obtained using SNIa alone. 

Open with DEXTER 
Table 3: Parameter means and 68%confidence intervals when, respectively, ignoring and taking into account the systematic errors in the form of photometric zeropoint fluctuations, see Sect. 2.2.4.
As mentioned in Sect. 2.2, we obtain unbiased bestfit values for the stretch and colour response parameters, and , respectively (see Table 4). These differ by about 15% from the (biased) values given in A06. The absolute magnitude M is consistent with A06.
The parameters describing the zeropoint shifts ( to ) are all consistent with zero (Table 4). Except for they have zero mean and rms of about . The influence of the derivative (Eq. (8)) on the second term of the likelihood (Eq. (9)) is small in comparison with the first term. Nevertheless, correlations with other parameters are introduced as is discussed in the next section. The mean of the zband zeropoint shifts is negative (although not significantly so) and its variance smaller than the expected value of . It is also this parameter which shows the highest correlation with cosmological parameters.
Table 4: Mean and 68% errors for SNIa internal parameters, for a flat wCDM model. The values of are in units of 0.01 mag.
2.2.5 Correlation between systematic errors and cosmological parameters
In most cases the zeropoint shifts are uncorrelated with other parameters. However, some pairs of one or more nuisance parameters show correlations, most notably, the absolute SNIa magnitude M is correlated both with the B and Vband uncertainty and , respectively, see Fig. 5. This correlation can be explained by looking at the restframe colour parameter, , which contains the BV colour excess at the time of the Bband maximum. Redefining the B or Vband zeropoints causes a systematic change in c which is compensated by a corresponding change in M. The slope of the measured correlation is about 1.5. The correlation between M and is much weaker and has a smaller slope of 0.3 to 0.4. The correlation between M and the difference is therefore on the order of the bestfit value of , which is expected from Eq. (6). The difference in correlation can be explained by the fact that the observed magnitudes of each object translate into different restframe bands for the lightcurve fit depending on the object's redshift.
Further, the zband zeropoint offset ( ) shows a correlation with both and w (Fig. 6). We find and . The latter value can be compared to A06 (Table 5) who varied the zeropoints to infer the influence on cosmological parameters; our value is of the same order of magnitude as the one from A06 but has opposite sign. For a fixed we obtain , about half of the value cited in A06 (Table 5).
Figure 5: Correlations between the universal absolute SNIa magnitude M and the Bband ( , left panel) and Vband ( , right panel), respectively. The zeropoint shift parameters are given in units of 0.01 mag. 

Open with DEXTER 
Figure 6: Correlations between the zband zeropoint offset ( ) and ( left panel) respectively w ( right panel). 

Open with DEXTER 
2.2.6 Further sources of systematics
A06 discussed further sources of systematic uncertainties which can affect cosmological results from SNIa. Some of those uncertainties could be parametrised and invoked in a joint analysis with cosmology as presented in this paper. We state a few more examples here and leave a thorough analysis for future work.
For example, the difference between measured restframe Uband magnitude and the one ``predicted'' from the lightcurve fit using only two bands can be included in the likelihood function.
In addition, the Malmquist bias can be modelled as a function of redshift. The inclusion of this bias will be necessary when redshiftdependent darkenergy models are to be tested. This requires accurate knowledge of the sample and of selection effects.
Furthermore, weak gravitational lensing of distant SNIa can be modelled using the darkmatter power spectrum, either from (nonlinear) theoretical prescriptions (Ménard & Dalal 2005) or numerical simulations (Martel & Premadi 2008).
2.3 WMAP5 CMB anisotropies
To calculate CMB temperature and polarisation power and crossspectra we use the publicly available package CAMB^{} (Lewis et al. 2000). The likelihood is evaluated using the public WMAP5 code^{} (Dunkley et al. 2009). We include all ``standard'' components which are the same for three and fiveyear. Those are the low TT and TE/EE/BB spectra, the high TT and TE spectra, and the pointsource TT correction. The low ( for TT, for polarisation) likelihoods are calculated using Gibbs sampling (Dunkley et al. 2009; Page et al. 2007). The high sampling uses pseudo according to Hinshaw et al. (2007).
In contrast to the published WMAP5 results, we do not include corrections due to SZ. The SZ amplitude is unconstrained by WMAP5 and is not degenerate with other parameters (Dunkley et al. 2009). A recent analysis found no biases for WMAP5 when ignoring the contribution of the thermal SZ effect from clusters (Taburet et al. 2009).
3 Parameter analysis and cosmological constraints
3.1 Cosmological model and parameters
We assume a flat darkenergy cold dark matter cosmology (wCDM) with the parameter vector . Clustering of dark energy is not taken into account. For CMB, instead of the normalisation defined at a scale of 8 Mpc/h, , the parameter which is sampled is , the curvature perturbations amplitude at the pivot scale . As stated in Sect. 2.1 we use the Smith et al. (2003) fitting formula to model the nonlinear power spectrum needed for weak lensing. We do not include a massive neutrino component to the massenergy tensor and assume as the effective number of massless neutrinos (Mangano et al. 2002), which is the preferred value for WMAP5 (Dunkley et al. 2009). Tereno et al. (2008) obtained constraints on the neutrino mass using various probes including CFHTLSWide weak lensing, see also Gong et al. (2008) and Ichiki et al. (2009). Additional parameters are (a,b,c,c_{0}) from lensing and from SN Ia. These parameters are described in the respective subsections of Sect. 2; including systematics we sample a total of 22 parameters.
3.2 Sampling the parameter space with MCMC
We use an adaptive MetropolisHastings algorithm (Metropolis et al. 1953; Hastings 1970) to generate Monte Carlo Markov Chains (MCMC) as sample of the posterior. The Fisher matrix evaluated at the maximumlikelihood (ML) parameter serves as initial multivariate Gaussian proposal. The ML point is estimated by a conjugategradient search (Press et al. 1992). For lensing and SNIa the typical chain length is 100 000 with an acceptance rate of about 1520%. Every 1000 steps the chain covariance is updated and replaces the previous proposal. The proposal variance is multiplied by a factor which is optimal for a Gaussian posterior (Dunkley et al. 2005; Hanson & Cunningham 1998).
For WMAP5 on the other hand, we choose a different strategy. With acceptance rates as stated above we found the posterior not well sampled and as a consequence the parameter uncertainties largely underestimated. We refrain from updating the proposal but kept the initial guess, which is larger than the Fisher matrix, for the whole duration of the chain. This choice results in small acceptance rates of 3% to 10% but guarantees an unbiased sampling of the posterior.
The experiments are combined by multiplying the corresponding likelihoods or, equivalently, by summing up the loglikelihoods. The confidence intervals are obtained by creating histograms of the parameter vectors of the chain and estimating regions with 68%, 95% and 99.7% of the density.
In particular for WMAP5 the calculation of the converged chain is very timeconsuming and takes a few days to a week on a fast multicore machine. To overcome this and other problems related to MCMC, such as the difficulty to assess convergence and the high correlation of a Markov chain, a new method called Population Monte Carlo (PMC, Cappé et al. 2007,2004) has been developed which is based on adaptive importance sampling. In a companion paper, we test this method with simulations and apply it to cosmology posteriors (Wraith et al. 2009).
3.3 Combined constraints from lensing and CMB
Figure 7: 68% and 95% confidence levels for lensing only, assuming wCDM (solid red lines), CDM (dashed red), and CMB only, for wCDM (solid blue), CDM (dashed blue curves). 

Open with DEXTER 
The main contribution from cosmic shear (without tomographic information), in combination with CMB anisotropy measurements, are constraints on the normalisation (Contaldi et al. 2003; Tereno et al. 2005). For a CDM model, the degeneracy direction between and is orthogonal between lensing and CMB. Lensing is sensitive to the amount of structure and therefore, a high normalisation has to be countered by a smaller amount of the total matter. On the other hand, an overall increase of the CMB angular power spectrum by a higher normalisation has to be compensated by more matter to dampen the peaks.
This picture changes if the dark energy equationofstate is a free parameter. Whereas the lensing contours for and do not broaden much, the CMB constraints increase dramatically. Most importantly, they increase along the lensingdirection of degeneracy, see Fig. 7. The result is that adding lensing to CMB data will not improve the constraints in the plane as much as it would do in the framework of a cosmological constant model.
Lensing without tomography is not very sensitive to w, thus leaving it a free parameter has a minor influence on the other parameters. For CMB however, strong correlations enter through the angular diameter distance which CMB measures only at z=1100. Those degeneracies (see Fig. 9) can be broken by adding an additional distance measurement at low redshift as a lever, e.g. using BAO or SNIa, (Blake & Glazebrook 2003; Frieman et al. 2003; Seo & Eisenstein 2003).
3.4 Combined constraints from lensing, SNIa and CMB
Table 5: CMB, lensing and SNIa in various combinations. The mean and 68% marginals are given. For the first four cases systematics are ignored, the last column includes all systematics, from both lensing and supernovae (see Sect. 2.2.4).
Figure 8: 68% and 95% confidence levels for the three parameters which are affected most by systematics ( ). Solid (dashed) contours correspond to the case of ignoring (including) systematics. 

Open with DEXTER 
Figure 9: 68% and 95% confidence levels for CMB (orange, solid lines), CMB+lensing (green dashed), CMB+SNIa (magenta dotted) and CMB+lensing+SNIa (blue dashdotted). Systematics are ignored in this plot. 

Open with DEXTER 
We will now discuss the main results of this paper. Joint constraints using lensing, SNIa and CMB are compared for the two cases with and without taking systematics into account. In the former, both lensing and supernovae systematics are included. The results are given in Table 5, Figs. 8 and 9. We show the bestfit angular and 3d power spectra in Fig. 10. A CAMB parameter file with our bestfit parameter values is available for download^{}.
With the current data, SNIa is more efficient than lensing in helping decrease uncertainties. CMB+SNIa gives nearly as tight constraints as CMB+SNIa+lensing for most parameters. The reason is that SNIa data show different degeneracy directions, in particular for the pair (see right panel of Fig. 4). This helps to pin down w and thus, the main degeneracy for CMB is largely lifted. As stated in the previous section, this is not the case for lensing  lensing without tomography cannot constrain w. The consequence is that even parameters to which SNIa is not sensitive, e.g. , are very accurately determined for the combination CMB+SNIa.
Nevertheless, lensing improves constraints from WMAP5 substantially. Some CMBrelated neardegeneracies which arise in the wCDM model are partially lifted in combination with lensing. In particular for ( ), there is a large gain when lensing data is added.
The effect of systematics on the parameters means and errors can be assessed by comparing the last two columns of Table 5. The shift of the bestfit values is less than 15% of the statistical error in all cases. Including systematics in the analysis increases the error bars by 1035%, where , w and are affected most. Varying the lensing redshift for the bias z_{0} from 1.0 to 0.8 and 1.2 changes the results by less then one percent.
Figure 10: The combined (lensing+SNIa+CMB) bestfit z=0 matter power spectrum ( left panel) and angular power spectrum ( right panel). The shaded region corresponds to the 68% uncertainty including all systematics. Note that this error does not correspond to the uncertainty due to cosmic variance. 

Open with DEXTER 
4 Discussion
In this paper we combine three different cosmological probes to test a possible deviation from a cosmological constant: (1) weak gravitational lensing as a probe of structure formation and geometry in the redshift range of about 0.2 to 0.8. (2) Supernovae Ia as standard(isable) candles up to redshift unity. (3) CMB anisotropies supplying a wealth of information of the recombination era ( ) and, to a lesser extend, of the Universe up to the reionisation epoch ( ). The data sets are (1) CFHTLSWide T0003 (see F08); (2) SNLS firstyear (A06) and (3) WMAP fiveyear (H09), respectively.
We test models in the context of a flat CDM cosmology with a dark or vacuumenergy component with free but constant equationofstate parameter . This corresponds to the simplest extension of the ``vanilla'' CDM model which goes beyond a cosmological constant.
The joint constraints including the modelled systematics in the data sets are w = 1.02^{+0.08 +0.14}_{0.08 0.16} (68% and 95% confidence, respectively). Without taking the systematics into account, the result is w = 1.01 0.06 0.12, representing 25% smaller error bars. With the current data there is no evidence for a dynamical darkenergy component not being the cosmological constant.
Two potential sources of bias in the thirdyear CFHTLSWide lensing data are scrutinised. One, the measured variations of the shear signal (aperturemass dispersion ) between MegaCam pointings are compared to Nbody simulations. We estimate the measured fluctuations to be higher than expected by not more than about 5% to 15% on scales below 30 arcmin. Whereas this might be a hint of systematics in the data, it is not straightforward to model and to assess its effect on cosmology. The second issue are systematics in the shape measurements which seem to lead to an underestimation of the lensing signal at high redshift. We devised a very simple model of this potential systematics by multiplying the lensing efficiency above z=1 with a constant c_{0}>0, to mimic the effect of a decreased measured shear. Marginalising over cosmological parameters, using weak lensing alone, yields c_{0} = 1.1 0.6. All probes combined do not constrain c_{0} much better, we find c_{0} = 1.1 0.5. Restricting ourselves to c_{0}<1, implying an underestimation of the lensing signal, increases by about 8% for a fixed , which is roughly equal to the statistical error. Therefore, in the framework of this simple model, the value of might be biased by 8% if this effect is ignored (as it has been in F08).
Combining all probes we assess the influence of systematics on the cosmological results. The contribution of systematics to the total error budget ranges between 5% and 40% (Table 5). The parameters which are affected most by systematics are and w for which this contribution is greater than 20%. Those are the parameters for which weak lensing adds significantly to the joint constraints.
There are indications of more unaccounted systematics in the lensing data. Discussions about the origin of those systematics and ways to remove them are addressed elsewhere (van Waerbeke et al., in prep.). Our findings strengthen the confidence in cosmological results from cosmic shear, presented in this work and earlier, using the same data, in F08.
Our constraints are slightly tighter than Seljak et al. (2005) who obtained w = 0.99^{+0.09 + 0.16}_{0.09 0.20} using WMAP1, SNIa and, as probes of structure formation, SDSS galaxy correlations and Ly forest clustering. Similar constraints were quoted in Wang & Mukherjee (2006), w = 1.00^{+0.08 + 0.16}_{0.08  0.17} from WMAP3, SNLS and SDSS. Our results are comparable to the ones stated by Komatsu et al. (2009), 0.11 < 1+w < 0.14 (95%) stemming from WMAP5, BAO and SNIa. Consistency is also achieved with the constraints from Mantz et al. (2008) who combined WMAP5 with the Xray cluster mass function, cluster baryon fraction and SNIa to get w = 1.02 0.06 (68.3%).
Earlier results including weak gravitational lensing used less wide and/or more shallow data. Using the firstyear release of the CFHTLS (T0001, Hoekstra et al. 2006), and with rather tight priors on other cosmological parameters, an upper bound of was derived. Jarvis et al. (2006) combined weak lensing from the 75 deg^{2} CTIO survey with SNIa and CMB and obtained . The error estimate for the latter includes statistical uncertainty, PSF systematics and shape measurement calibration biases.
This paper puts particular emphasis on the treatment of systematics in the data. Known observationrelated systematics are parametrised if possible and included in the analysis together with cosmological parameters. This allows us to directly quantify the influence of systematics on cosmology and to find possible correlations. For SNIa, we calculate the response of the distance modulus to fluctuations of the photometry zeropoints. Ignoring these error sources leads to an underestimation of the parameter errors by 10%. The parameter values itself are biased by a fraction of 10% to 20% of the statistical uncertainty.
5 Outlook
For future, highprecision experiments it will be of great importance to understand and control systematic effects. In particular for exotic models, such as timevarying darkenergy or modifications of GR, more and more subtle influences of systematics mimicking a signal have to be excluded. This calls for combined analyses of cosmological and systematics parameters.
Weak lensing serves as an important and independent probe of . The measurement of the lensing skewness will allow us to lift degeneracies between parameters, in particular between and (Kilbinger & Schneider 2005). Power and bispectrum tomography can be used to constrain the timeevolution of the darkenergy eos parameter (Takada & Jain 2004). In addition, as a complementary and nearly independent probe of the nonGaussian, highdensity regime of the largescale structure, weak lensing cluster counts can improve darkenergy parameter constraints (Takada & Bridle 2007).
For lensing, systematic effects which have to be taken into account for future analyses include measurement errors (galaxy shapes, PSF correction, photoz's), astrophysical sources (intrinsic alignments, source clustering) and theoretical uncertainties (nonlinear and baryonic physics). The modelling of all those effects, as suggested recently by Bernstein (2008), leads to a huge number of parameters, on the order of several dozens to hundreds, depending on the number of redshift bins. In case of simultaneous determination of cosmic shear and intrinsic alignments a high number of redshift bins is required (Bridle & King 2007; Joachimi & Schneider 2008).
Already with the present data, more general models of dark energy can be constrained. A wide variety of such models are tested with recent SNIa data (Rubin et al. 2008). Because dark energy and modified gravity can influence distances in a nondistinguishable way, it is important to include measures of the growth of structure. A survey like the CFHTLSWide reaches out into the linear regime and is less prone to smallscale uncertainties due to baryonic and nonlinear physics, which makes such a survey an excellent probe of the present acceleration of the Universe.
Acknowledgements
We acknowledge the use of the Legacy Archive for Microwave Background Data Analysis (LAMBDA). Support for LAMBDA is provided by the NASA Office of Space Science. We thank the TERAPIX group for support and computational facilities. We acknowledge the CFHTLS lensing systematics collaboration for uncovering the weaklensing systematics that are parametrised and marginalised over in this analysis. We thank O. Cappé, J.F. Cardoso and O. Doré for helpful discussions, and E. Bertin, S. Prunet and B. Rowe for valuable comments on the manuscript. We also would like to thank the anonymous referee for helpful suggestions which improved the paper. MK is supported by the CNRS ANR ``ECOSSTAT'', contract number ANR05BLAN028304. IT and LF acknowledge the support of the European Commission Programme 6th framework, Marie Curie Training and Research Network ``DUEL'', contract number MRTNCT2006036133.
References
 Astier, P., Guy, J., Regnault, N., et al. 2006, A&A, 447, 31 [NASA ADS] [CrossRef] [EDP Sciences] (A06) (In the text)
 Bartelmann, M., Doran, M., & Wetterich, C. 2006, A&A, 454, 27 [NASA ADS] [CrossRef] [EDP Sciences]
 Benjamin, J., Heymans, C., Semboloni, E., et al. 2007, MNRAS, 381, 702 [NASA ADS] [CrossRef] (In the text)
 Bernstein, G. M. 2008, ApJ, submitted [arXiv:0808.3400] (In the text)
 Blake, C., & Glazebrook, K. 2003, ApJ, 594, 665 [NASA ADS] [CrossRef]
 Bridle, S., & King, L. 2007, New J. Phys., 9, 444 [CrossRef]
 Cappé, O., Guillin, A., Marin, J.M., & Robert, C. 2004, J. Comput. Graph. Statist., 13, 907 [CrossRef]
 Cappé, O., Douc, R., Guillin, A., Marin, J.M., & Robert, C. 2007, in Statistics and Computing, in press [arXiv:0710.4242]
 Contaldi, C. R., Hoekstra, H., & Lewis, A. 2003, Phys. Rev. Lett., 90, 221303 [NASA ADS] [CrossRef]/1 (In the text)
 Coupon, J., Ilbert, O., Kilbinger, M., et al. 2008, A&A, submitted[arXiv:0811.3326] (In the text)
 Dunkley, J., Bucher, M., Ferreira, P. G., Moodley, K., & Skordis, C. 2005, MNRAS, 356, 925 [NASA ADS] [CrossRef]
 Dunkley, J., Komatsu, E., Nolta, M. R., et al. 2009, ApJS, 180, 306 [NASA ADS] [CrossRef]
 Eifler, T., Schneider, P., & Hartlap, J. 2008, A&A, submitted[arXiv:0810.4254] (In the text)
 Francis, M. J., Lewis, G. F., & Linder, E. V. 2008, MNRAS, accepted[arXiv:0808.2840]
 Frieman, J. A., Huterer, D., Linder, E. V., & Turner, M. S. 2003, Phys. Rev. D, 67, 083505 [NASA ADS] [CrossRef]
 Fu, L. 2008, Ph.D. Thesis, Université Paris 6 (In the text)
 Fu, L., Semboloni, E., Hoekstra, H., et al. 2008, A&A, 479, 9 [NASA ADS] [CrossRef] [EDP Sciences] (F08) (In the text)
 Gong, Y., Zhang, T.J., Lan, T., & Chen, X.L. 2008 [arXiv:0810.3572] (In the text)
 Hamana, T., & Mellier, Y. 2001, MNRAS, 327, 169 [NASA ADS] [CrossRef] (In the text)
 Hanson, K. M., & Cunningham, G. S. 1998, in Medical Imaging: Image Processing, 371
 Hastings, W. K. 1970, Biometrika, 57, 97 [CrossRef]
 Heymans, C., Van Waerbeke, L., Bacon, D., et al. 2006, MNRAS, 368, 1323 [NASA ADS] [CrossRef]
 Hinshaw, G., Nolta, M. R., Bennett, C. L., et al. 2007, ApJS, 170, 288 [NASA ADS] [CrossRef] (In the text)
 Hinshaw, G., Weiland, J. L., Hill, R. S., et al. 2009, ApJS, 180, 225 [NASA ADS] [CrossRef] (H09) (In the text)
 Hoekstra, H., Yee, H. K. C., & Gladders, M. D. 2002, ApJ, 577, 595 [NASA ADS] [CrossRef] (In the text)
 Hoekstra, H., Mellier, Y., van Waerbeke, L., et al. 2006, ApJ, 457, 116 [NASA ADS] [CrossRef] (In the text)
 Ichiki, K., Takada, M., & Takahashi, T. 2009, Phys. Rev. D, 79, 023520 [NASA ADS] [CrossRef] (In the text)
 Ilbert, O., Arnouts, S., McCracken, H. J., et al. 2006, A&A, 457, 841 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Jain, B., & Zhang, P. 2008, Phys. Rev. D, 78, 063503 [NASA ADS] [CrossRef] (In the text)
 Jarvis, M., Jain, B., Bernstein, G., & Dolney, D. 2006, ApJ, 644, 71 [NASA ADS] [CrossRef] (In the text)
 Joachimi, B., & Schneider, P. 2008, A&A, 488, 829 [NASA ADS] [CrossRef] [EDP Sciences]
 Kaiser, N. 1992, ApJ, 388, 272 [NASA ADS] [CrossRef] (In the text)
 Kilbinger, M., & Schneider, P. 2005, A&A, 442, 69 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Kilbinger, M., & Munshi, D. 2006, MNRAS, 366, 983 [NASA ADS] (In the text)
 Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, ApJS, 180, 330 [NASA ADS] [CrossRef]
 Kowalski, M., Rubin, D., Aldering, G., et al. 2008, ApJ, 686, 749 [NASA ADS] [CrossRef] (In the text)
 Lesgourgues, J., Viel, M., Haehnelt, M. G., & Massey, R. 2007, J. Cosmol. AstroPart. Phys., 11, 8 [NASA ADS] [CrossRef] (In the text)
 Lewis, A., Challinor, A., & Lasenby, A. 2000, ApJ, 538, 473 [NASA ADS] [CrossRef] (In the text)
 Ma, Z. 2007, ApJ, 665, 887 [NASA ADS] [CrossRef] (In the text)
 Mangano, G., Miele, G., Pastor, S., & Peloso, M. 2002, Phys. Lett. B, 534, 8 [NASA ADS] [CrossRef] (In the text)
 Mantz, A., Allen, S. W., Ebeling, H., & Rapetti, D. 2008, MNRAS, 387, 1179 [NASA ADS] [CrossRef] (In the text)
 Martel, H., & Premadi, P. 2008, ApJ, 673, 657 [NASA ADS] [CrossRef] (In the text)
 Massey, R., Heymans, C., Bergé, J., et al. 2007, MNRAS, 376, 13 [NASA ADS] [CrossRef]
 Ménard, B., & Dalal, N. 2005, MNRAS, 358, 101 [NASA ADS] [CrossRef] (In the text)
 Metropolis, N., Rosenbluth, A., Rosenbluth, M., Teller, A., & Teller, E. 1953, J. Chem. Phys., 21, 1087 [NASA ADS] [CrossRef]
 Page, L., Hinshaw, G., Komatsu, E., et al. 2007, ApJS, 170, 335 [NASA ADS] [CrossRef]
 Peacock, J. A., & Dodds, S. J. 1996, MNRAS, 280, L19 [NASA ADS] (In the text)
 Press, W. H., Teukolsky, S. A., Flannery, B. P., & Vetterling, W. T. 1992, Numerical Recipes in C (Cambridge University Press) (In the text)
 Rubin, D., Linder, E. V., Kowalski, M., et al. 2008, ArXiv eprints (In the text)
 Schimd, C., Tereno, I., Uzan, J.P., et al. 2007, A&A, 463, 405 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Schneider, P., van Waerbeke, L., Jain, B., & Kruse, G. 1998, MNRAS, 296, 873 [NASA ADS] [CrossRef] (In the text)
 Seljak, U., Makarov, A., McDonald, P., et al. 2005, Phys. Rev. D, 71, 103515 [NASA ADS] [CrossRef] (In the text)
 Seo, H.J., & Eisenstein, D. J. 2003, ApJ, 598, 720 [NASA ADS] [CrossRef]
 Smith, R. E., Peacock, J. A., Jenkins, A., et al. 2003, MNRAS, 341, 1311 [NASA ADS] [CrossRef] (In the text)
 Taburet, N., Aghanim, N., Douspis, M., & Langer, M. 2009, MNRAS, 392, 1153 [NASA ADS] [CrossRef] (In the text)
 Takada, M., & Bridle, S. 2007, New J. Phys., 9, 446 [CrossRef] (In the text)
 Takada, M., & Jain, B. 2004, MNRAS, 348, 897 [NASA ADS] [CrossRef] (In the text)
 Tereno, I., Doré, O., van Waerbeke, L., & Mellier, Y. 2005, A&A, 429, 383 [NASA ADS] [CrossRef] [EDP Sciences]
 Tereno, I., Schimd, C., Uzan, J.P., et al. 2008, A&A, accepted[arXiv:0810.0555] (In the text)
 van Waerbeke, L., White, M., Hoekstra, H., & Heymans, C. 2006, Astrop. Phys., 26, 91 [NASA ADS] [CrossRef] (In the text)
 Wang, Y., & Mukherjee, P. 2006, ApJ, 650, 1 [NASA ADS] [CrossRef] (In the text)
 Wraith, D., Kilbinger, M., Benabed, K., Robert, C., et al. 2009, Phys. Rev. D, submitted (In the text)
 Xia, J.Q., Li, H., Zhao, G.B., & Zhang, X. 2008, Phys. Rev. D, 78, 083524 [NASA ADS] [CrossRef] (In the text)
 Yahagi, H. 2005, PASJ, 57, 779 [NASA ADS]
 Yahagi, H., & Yoshii, Y. 2001, ApJ, 558, 463 [NASA ADS] [CrossRef]
Footnotes
 ... WMAP5^{}
 Based on observations obtained with M EGAP RIME/M EGAC AM, 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 Sciences de l'Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii. This work is based in part on data products produced at T ERAPIX and the Canadian Astronomy Data Centre as part of the CanadaFranceHawaii Telescope Legacy Survey, a collaborative project of NRC and CNRS.
 ... (CFHTLS^{}
 http://www.cfht.hawaii.edu/Science/CFHTLS
 ... KIDS^{}
 http://www.astrowise.org/projects/KIDS
 ... DES^{}
 https://www.darkenergysurvey.org
 ... LSST^{}
 www.lsst.org
 ... JDEM^{}
 http://universe.nasa.gov/program/probes/jdem.html
 ... Euclid^{}
 http://sci.esa.int/sciencee/www/area/index.cfm?fareaid=102
 ... CAMB^{}
 http://camb.info
 ... WMAP5 code^{}
 http://lambda.gsfc.nasa.gov
 ... download^{}
 http://www2.iap.fr/users/kilbinge/params.ini
All Tables
Table 1: Shear tophat variance at 5' and 25' and for five redshift ranges. The error bars are Poisson noise.
Table 2: List of SNIa systematic parameters and their symbols.
Table 3: Parameter means and 68%confidence intervals when, respectively, ignoring and taking into account the systematic errors in the form of photometric zeropoint fluctuations, see Sect. 2.2.4.
Table 4: Mean and 68% errors for SNIa internal parameters, for a flat wCDM model. The values of are in units of 0.01 mag.
Table 5: CMB, lensing and SNIa in various combinations. The mean and 68% marginals are given. For the first four cases systematics are ignored, the last column includes all systematics, from both lensing and supernovae (see Sect. 2.2.4).
All Figures
Figure 1: Variance of from the numerical simulations (solid lines) and CFHTLS data (dashed). The rms used for the merged catalogue (from Fu et al. 2008), which we use for the cosmological constraints, is the dotted line. 

Open with DEXTER  
In the text 
Figure 2: Left panel: 2d marginals (68% and 95%) for and for the two cases of including the highz calibration bias (solid lines) and ignoring it (dashed curves), respectively. Right panel: 1d marginal likelihood for c_{0}. The vertical lines indicate mean (dashed) and 68% confidence intervals (dotted). 

Open with DEXTER  
In the text 
Figure 3: Left panel: 1d marginal for for , the vertical lines are as in Fig. 2. Right panel: correlation between c_{0} and (68% and 95% confidence levels) for . The vertical lines indicate mean (dashed) and 68% confidence intervals (dotted). 

Open with DEXTER  
In the text 
Figure 4: Confidence contours (68%, 95.5%, 99.7%) for full treatment of errors due to photometry zeropoints (blue solid lines) and ignoring those errors (green dashed curves). These two cases correspond to the two cases in Table 3. The left panel corresponds to w=1, the right panel is for a flat wCDM Universe. The constraints are obtained using SNIa alone. 

Open with DEXTER  
In the text 
Figure 5: Correlations between the universal absolute SNIa magnitude M and the Bband ( , left panel) and Vband ( , right panel), respectively. The zeropoint shift parameters are given in units of 0.01 mag. 

Open with DEXTER  
In the text 
Figure 6: Correlations between the zband zeropoint offset ( ) and ( left panel) respectively w ( right panel). 

Open with DEXTER  
In the text 
Figure 7: 68% and 95% confidence levels for lensing only, assuming wCDM (solid red lines), CDM (dashed red), and CMB only, for wCDM (solid blue), CDM (dashed blue curves). 

Open with DEXTER  
In the text 
Figure 8: 68% and 95% confidence levels for the three parameters which are affected most by systematics ( ). Solid (dashed) contours correspond to the case of ignoring (including) systematics. 

Open with DEXTER  
In the text 
Figure 9: 68% and 95% confidence levels for CMB (orange, solid lines), CMB+lensing (green dashed), CMB+SNIa (magenta dotted) and CMB+lensing+SNIa (blue dashdotted). Systematics are ignored in this plot. 

Open with DEXTER  
In the text 
Figure 10: The combined (lensing+SNIa+CMB) bestfit z=0 matter power spectrum ( left panel) and angular power spectrum ( right panel). The shaded region corresponds to the 68% uncertainty including all systematics. Note that this error does not correspond to the uncertainty due to cosmic variance. 

Open with DEXTER  
In the text 
Copyright ESO 2009