
We compute the sphericalsky weaklensing power spectrum of the shear and
convergence. We discuss various approximations, such as flatsky, and first
and second order Limber equations for the projection. We find that the impact
of adopting these approximations is negligible when constraining cosmological
parameters from current weak lensing surveys. This is demonstrated using data
from the CanadaFranceHawaii Telescope Lensing Survey (CFHTLenS). We find that
the reported tension with Planck Cosmic Microwave Background (CMB) temperature
anisotropy results cannot be alleviated.
For future largescale surveys with unprecedented precision, we show that the
spherical secondorder Limber approximation will provide sufficient accuracy.
In this case, the cosmicshear power spectrum is shown to be in agreement with
the full projection at the subpercent level for l > 3, with the corresponding
errors an order of magnitude below cosmic variance for all l. When computing
the twopoint shear correlation function, we show that the flatsky fast Hankel
transformation results in errors below two percent compared to the full
spherical transformation. In the spirit of reproducible research, our numerical
implementation of all approximations and the full projection are publicly
available within the package nicaea at
http://www.cosmostat.org/software/nicaea.

Weaklensing peak counts provide a straightforward way to constrain cosmology
by linking local maxima of the lensing signal to the mass function. Recent
applications to data have already been numerous and fruitful. However, the
importance of understanding and dealing with systematics increases as data
quality reaches an unprecedented level. One of the sources of systematics is
the convergenceshear inversion. This effect, inevitable when carrying out a
convergence field from observations, is usually neglected by theoretical peak
models. Thus, it could have an impact on cosmological results. In this paper,
we study the bias from neglecting (mismodeling) the inversion. Our tests show
a small but nonnegligible bias. The cosmological dependence of this bias seems
to be related to the parameter Sigma_8. When this bias propagates to the
parameter estimation, we discovered that constraint contours involving the dark
energy equation of state can differ by 2sigma. Such an effect can be even
larger for future highprecision surveys and we argue that the inversion should
be properly modeled for theoretical peak models.

We present a study of the dependencies of shear bias on simulation (input)
and measured (output) parameters, noise, pointspread function anisotropy,
pixel size, and the model bias coming from two different and independent galaxy
shape estimators. We used simulated images from Galsim based on the GREAT3
controlspaceconstant branch, and we measured shear bias from a modelfitting
method (gFIT) and a momentbased method (KaiserSquiresBroadhurst). We show
the bias dependencies found on input and output parameters for both methods,
and we identify the main dependencies and causes. Most of the results are
consistent between the two estimators, an interesting result given the
differences of the methods. We also find important dependences on orientation
and morphology properties such as flux, size, and ellipticity. We show that
noise and pixelization play an important role in the bias dependencies on the
output properties and galaxy orientation. We show some examples of model bias
that produce a bias dependence on the Sersic index n as well as a different
shear bias between galaxies consisting of a single Sersic profile and galaxies
with a disc and a bulge. We also see an important coupling between several
properties on the bias dependences. Because of this, we need to study several
measured properties simultaneously in order to properly understand the nature
of shear bias. This paper serves as a first step towards a companion paper that
describes a machine learning approach to modelling shear bias as a complex
function of many observed properties.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations (BAO) and the
growth of structure through redshiftspace distortions with a widearea galaxy
and quasar redshift survey. To trace the underlying dark matter distribution,
spectroscopic targets will be selected in four classes from imaging data. We
will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to
even higher redshift, DESI will target bright [O II] emission line galaxies up
to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying
dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the
Ly$\alpha$ forest absorption features in their spectra, which will be used to
trace the distribution of neutral hydrogen. When moonlight prevents efficient
observations of the faint targets of the baseline survey, DESI will conduct a
magnitudelimited Bright Galaxy Survey comprising approximately 10 million
galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy
and quasar redshifts will be obtained to measure the BAO feature and determine
the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations and the growth
of structure through redshiftspace distortions with a widearea galaxy and
quasar redshift survey. The DESI instrument is a roboticallyactuated,
fiberfed spectrograph capable of taking up to 5,000 simultaneous spectra over
a wavelength range from 360 nm to 980 nm. The fibers feed ten threearm
spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500,
depending on wavelength. The DESI instrument will be used to conduct a
fiveyear survey designed to cover 14,000 deg$^2$. This powerful instrument
will be installed at prime focus on the 4m Mayall telescope in Kitt Peak,
Arizona, along with a new optical corrector, which will provide a threedegree
diameter field of view. The DESI collaboration will also deliver a
spectroscopic pipeline and data management system to reduce and archive all
data for eventual public use.

Peak statistics in weak lensing maps access the nonGaussian information
contained in the largescale distribution of matter in the Universe. They are
therefore a promising complement to twopoint and higherorder statistics to
constrain our cosmological models. To prepare for the highprecision data of
nextgeneration surveys, we assess the constraining power of peak counts in a
simulated Euclidlike survey on the cosmological parameters
$\Omega_\mathrm{m}$, $\sigma_8$, and $w_0^\mathrm{de}$. In particular, we study
how the Camelus modela fast stochastic algorithm for predicting peakscan be
applied to such large surveys. We measure the peak count abundance in a mock
shear catalogue of ~5,000 sq. deg. using a multiscale mass map filtering
technique. We then constrain the parameters of the mock survey using Camelus
combined with approximate Bayesian computation (ABC). We find that peak
statistics yield a tight but significantly biased constraint in the
$\sigma_8$$\Omega_\mathrm{m}$ plane, indicating the need to better understand
and control the model's systematics. We calibrate the model to remove the bias
and compare results to those from the twopoint correlation functions (2PCF)
measured on the same field. In this case, we find the derived parameter
$\Sigma_8=\sigma_8(\Omega_\mathrm{m}/0.27)^\alpha=0.76_{0.03}^{+0.02}$ with
$\alpha=0.65$ for peaks, while for 2PCF the value is
$\Sigma_8=0.76_{0.01}^{+0.02}$ with $\alpha=0.70$. We therefore see comparable
constraining power between the two probes, and the offset of their
$\sigma_8$$\Omega_\mathrm{m}$ degeneracy directions suggests that a combined
analysis would yield tighter constraints than either measure alone. As
expected, $w_0^\mathrm{de}$ cannot be well constrained without a tomographic
analysis, but its degeneracy directions with the other two varied parameters
are still clear for both peaks and 2PCF. (abridged)

This is the third in a series of papers that develop a new and flexible model
to predict weaklensing (WL) peak counts, which have been shown to be a very
valuable nonGaussian probe of cosmology. In this paper, we compare the
cosmological information extracted from WL peak counts using different
filtering techniques of the galaxy shear data, including linear filtering with
a Gaussian and two compensated filters (the starlet wavelet and the aperture
mass), and the nonlinear filtering method MRLens. We present improvements to
our model that account for realistic survey conditions, which are masks,
sheartoconvergence transformations, and nonconstant noise. We create
simulated peak counts from our stochastic model, from which we obtain
constraints on the matter density $\Omega_\mathrm{m}$, the power spectrum
normalisation $\sigma_8$, and the darkenergy parameter $w_0$. We use two
methods for parameter inference, a copula likelihood, and approximate Bayesian
computation (ABC). We measure the contour width in the
$\Omega_\mathrm{m}$$\sigma_8$ degeneracy direction and the figure of merit to
compare parameter constraints from different filtering techniques. We find that
starlet filtering outperforms the Gaussian kernel, and that including peak
counts from different smoothing scales helps to lift parameter degeneracies.
Peak counts from different smoothing scales with a compensated filter show very
little crosscorrelation, and adding information from different scales can
therefore strongly enhance the available information. Measuring peak counts
separately from different scales yields tighter constraints than using a
combined peak histogram from a single map that includes multiscale information.
Our results suggest that a compensated filter function with counts included
separately from different smoothing scales yields the tightest constraints on
cosmological parameters from WL peaks.

The XXL survey is the largest survey carried out by XMMNewton. Covering an
area of 50deg$^2$, the survey contains $\sim450$ galaxy clusters out to a
redshift $\sim$2 and to an Xray flux limit of
$\sim5\times10^{15}erg\,s^{1}cm^{2}$. This paper is part of the first
release of XXL results focussed on the bright cluster sample. We investigate
the scaling relation between weaklensing mass and Xray temperature for the
brightest clusters in XXL. The scaling relation is used to estimate the mass of
all 100 clusters in XXL100GC. Based on a subsample of 38 objects that lie
within the intersection of the northern XXL field and the publicly available
CFHTLenS catalog, we derive the $M_{WL}$ of each system with careful
considerations of the systematics. The clusters lie at $0.1<z<0.6$ and span a
range of $ T\simeq15keV$. We combine our sample with 58 clusters from the
literature, increasing the range out to 10keV. To date, this is the largest
sample of clusters with $M_{WL}$ measurements that has been used to study the
masstemperature relation. The fit ($M\propto T^b$) to the XXL clusters returns
a slope $b=1.78^{+0.37}_{0.32}$ and intrinsic scatter $\sigma_{\ln
MT}\simeq0.53$; the scatter is dominated by disturbed clusters. The fit to the
combined sample of 96 clusters is in tension with selfsimilarity,
$b=1.67\pm0.12$ and $\sigma_{\ln MT}\simeq0.41$. Overall our results
demonstrate the feasibility of groundbased weaklensing scaling relation
studies down to cool systems of $\sim1keV$ temperature and highlight that the
current data and samples are a limit to our statistical precision. As such we
are unable to determine whether the validity of hydrostatic equilibrium is a
function of halo mass. An enlarged sample of cool systems, deeper weaklensing
data, and robust modelling of the selection function will help to explore these
issues further.

Peak counts have been shown to be an excellent tool to extract the
nonGaussian part of the weak lensing signal. Recently, we developped a fast
stochastic forward model to predict weaklensing peak counts. Our model is able
to reconstruct the underlying distribution of observables for analyses. In this
work, we explore and compare various strategies for constraining parameter
using our model, focusing on the matter density $\Omega_\mathrm{m}$ and the
density fluctuation amplitude $\sigma_8$. First, we examine the impact from the
cosmological dependency of covariances (CDC). Second, we perform the analysis
with the copula likelihood, a technique which makes a weaker assumption
compared to the Gaussian likelihood. Third, direct, nonanalytic parameter
estimations are applied using the full information of the distribution. Fourth,
we obtain constraints with approximate Bayesian computation (ABC), an
efficient, robust, and likelihoodfree algorithm based on acceptreject
sampling. We find that neglecting the CDC effect enlarges parameter contours by
22%, and that the covariancevarying copula likelihood is a very good
approximation to the true likelihood. The direct techniques work well in spite
of noisier contours. Concerning ABC, the iterative process converges quickly to
a posterior distribution that is in an excellent agreement with results from
our other analyses. The time cost for ABC is reduced by two orders of
magnitude. The stochastic nature of our weaklensing peak count model allows us
to use various techniques that approach the true underlying probability
distribution of observables, without making simplifying assumptions. Our work
can be generalized to other observables where forward simulations provide
samples of the underlying distribution.

We present weak lensing constraints on the ellipticity of galaxyscale matter
haloes and the galaxyhalo misalignment. Using data from the
CanadaFranceHawaii Telescope Lensing Survey (CFHTLenS), we measure the
weightedaverage ratio of the aligned projected ellipticity components of
galaxy matter haloes and their embedded galaxies, $f_\mathrm{h}$, split by
galaxy type. We then compare our observations to measurements taken from the
Millennium Simulation, assuming different models of galaxyhalo misalignment.
Using the Millennium Simulation we verify that the statistical estimator used
removes contamination from cosmic shear. We also detect an additional signal in
the simulation, which we interpret as the impact of intrinsic shapeshear
alignments between the lenses and their largescale structure environment.
These alignments are likely to have caused some of the previous observational
constraints on $f_\mathrm{h}$ to be biased high. From CFHTLenS we find
$f_\mathrm{h}=0.04 \pm 0.25$ for earlytype galaxies, which is consistent with
current models for the galaxyhalo misalignment predicting $f_\mathrm{h}\simeq
0.20$. For latetype galaxies we measure $f_\mathrm{h}=0.69_{0.36}^{+0.37}$
from CFHTLenS. This can be compared to the simulated results which yield
$f_\mathrm{h}\simeq 0.02$ for misaligned latetype models.

Dark energy may be the first sign of new fundamental physics in the Universe,
taking either a physical form or revealing a correction to Einsteinian gravity.
Weak gravitational lensing and galaxy peculiar velocities provide complementary
probes of General Relativity, and in combination allow us to test modified
theories of gravity in a unique way. We perform such an analysis by combining
measurements of cosmic shear tomography from the CanadaFrance Hawaii Telescope
Lensing Survey (CFHTLenS) with the growth of structure from the WiggleZ Dark
Energy Survey and the SixdegreeField Galaxy Survey (6dFGS), producing the
strongest existing joint constraints on the metric potentials that describe
general theories of gravity. For scaleindependent modifications to the metric
potentials which evolve linearly with the effective dark energy density, we
find presentday cosmological deviations in the Newtonian potential and
curvature potential from the prediction of General Relativity to be (Delta
Psi)/Psi = 0.05 \pm 0.25 and (Delta Phi)/Phi = 0.05 \pm 0.3 respectively (68
per cent CL).

Cosmic shear is the distortion of images of distant galaxies due to weak
gravitational lensing by the largescale structure in the Universe. Such images
are coherently deformed by the tidal field of matter inhomogeneities along the
line of sight. By measuring galaxy shape correlations, we can study the
properties and evolution of structure on large scales as well as the geometry
of the Universe. Thus, cosmic shear has become a powerful probe into the nature
of dark matter and the origin of the current accelerated expansion of the
Universe. Over the last years, cosmic shear has evolved into a reliable and
robust cosmological probe, providing measurements of the expansion history of
the Universe and the growth of its structure.
We review here the principles of weak gravitational lensing and show how
cosmic shear is interpreted in a cosmological context. Then we give an overview
of weaklensing measurements, and present the main observational cosmicshear
results since it was discovered 15 years ago, as well as the implications for
cosmology. We then conclude with an outlook on the various future surveys and
missions, for which cosmic shear is one of the main science drivers, and
discuss promising new weak cosmological lensing techniques for future
observations.

We present first results from the third GRavitational lEnsing Accuracy
Testing (GREAT3) challenge, the third in a sequence of challenges for testing
methods of inferring weak gravitational lensing shear distortions from
simulated galaxy images. GREAT3 was divided into experiments to test three
specific questions, and included simulated space and groundbased data with
constant or cosmologicallyvarying shear fields. The simplest (control)
experiment included parametric galaxies with a realistic distribution of
signaltonoise, size, and ellipticity, and a complex point spread function
(PSF). The other experiments tested the additional impact of realistic galaxy
morphology, multiple exposure imaging, and the uncertainty about a
spatiallyvarying PSF; the last two questions will be explored in Paper II. The
24 participating teams competed to estimate lensing shears to within systematic
error tolerances for upcoming StageIV dark energy surveys, making 1525
submissions overall. GREAT3 saw considerable variety and innovation in the
types of methods applied. Several teams now meet or exceed the targets in many
of the tests conducted (to within the statistical errors). We conclude that the
presence of realistic galaxy morphology in simulations changes shear
calibration biases by $\sim 1$ per cent for a wide range of methods. Other
effects such as truncation biases due to finite galaxy postage stamps, and the
impact of galaxy type as measured by the S\'{e}rsic index, are quantified for
the first time. Our results generalize previous studies regarding sensitivities
to galaxy size and signaltonoise, and to PSF properties such as seeing and
defocus. Almost all methods' results support the simple model in which additive
shear biases depend linearly on PSF ellipticity.

Weaklensing peak counts has been shown to be a powerful tool for cosmology.
It provides nonGaussian information of large scale structures, complementary
to second order statistics. We propose a new flexible method to predict weak
lensing peak counts, which can be adapted to realistic scenarios, such as a
real source distribution, intrinsic galaxy alignment, mask effects, photo$z$
errors from surveys, etc. The new model is also suitable for applying the
tomography technique and nonlinear filters. A probabilistic approach to model
peak counts is presented. First, we sample halos from a mass function. Second,
we assign them NFW profiles. Third, we place those halos randomly on the field
of view. The creation of these "fast simulations" requires much less computing
time compared to $N$body runs. Then, we perform raytracing through these fast
simulation boxes and select peaks from weaklensing maps to predict peak number
counts. The computation is achieved by our \textsc{Camelus} algorithm, which we
make available at http://www.cosmostat.org/software/camelus/ . We compare our
results to $N$body simulations to validate our model. We find that our
approach is in good agreement with full $N$body runs. We show that the lensing
signal dominates shape noise and Poisson noise for peaks with SNR between 4 and
6. Also, counts from the same SNR range are sensitive to $\Omega_\mathrm{m}$
and $\sigma_8$. We show how our model can discriminate between various
combinations of those two parameters. In summary, we offer a powerful tool to
study weak lensing peaks. The potential of our forward model is its high
flexibility, making the use of peak counts under realistic survey conditions
feasible.

Using a sample of 25683 starforming and 2821 passive galaxies at $z\sim2$,
selected in the COSMOS field following the BzK color criterion, we study the
hosting halo mass and environment of galaxies as a function of their physical
properties. Spitzer and Herschel provide accurate SFR estimates for starburst
galaxies. We measure the auto and crosscorrelation functions of various
galaxy subsamples and infer the properties of their hosting halos using both
an HOD model and the linear bias at large scale. We find that passive and
starforming galaxies obey a similarly rising relation between the halo and
stellar mass. The mean host halo mass of star forming galaxies increases with
the star formation rate between 30 and 200 M$_\odot$.yr$^{1}$, but flattens
for higher values, except if we select only mainsequence galaxies. This
reflects the expected transition from a regime of secular coevolution of the
halos and the galaxies to a regime of episodic starburst. We find similar large
scale biases for mainsequence, passive, and starburst galaxies at equal
stellar mass, suggesting that these populations live in halos of the same mass.
We detect an excess of clustering on small scales for passive galaxies and
showed, by measuring the largescale bias of close pairs, that this excess is
caused by a small fraction ($\sim16%$) of passive galaxies being hosted by
massive halos ($\sim 3 \times 10^{13}$ M$_\odot$) as satellites. Finally,
extrapolating the growth of halos hosting the z$\sim$2 population, we show that
M$_\star \sim 10^{10}$ M$_\odot$ galaxies at z$\sim$2 will evolve, on average,
into massive (M$_\star \sim 10^{11}$ M$_\odot$), field galaxies in the local
Universe and M$_\star \sim 10^{11}$ M$_\odot$ galaxies at z=2 into local,
massive, group galaxies. The most massive mainsequence galaxies and close
pairs of massive, passive galaxies end up in today's clusters.

Higherorder, nonGaussian aspects of the largescale structure carry
valuable information on structure formation and cosmology, which is
complementary to secondorder statistics. In this work we measure second and
thirdorder weaklensing aperturemass moments from CFHTLenS and combine those
with CMB anisotropy probes. The third moment is measured with a significance of
$2\sigma$. The combined constraint on $\Sigma_8 = \sigma_8 (\Omega_{\rm
m}/0.27)^\alpha$ is improved by 10%, in comparison to the secondorder only,
and the allowed ranges for $\Omega_{\rm m}$ and $\sigma_8$ are substantially
reduced. Including general triangles of the lensing bispectrum yields tighter
constraints compared to probing mainly equilateral triangles. Second and
thirdorder CFHTLenS lensing measurements improve Planck CMB constraints on
$\Omega_{\rm m}$ and $\sigma_8$ by 26% for flat $\Lambda$CDM. For a model with
free curvature, the joint CFHTLenSPlanck result is $\Omega_{\rm m} = 0.28 \pm
0.02$ (68% confidence), which is an improvement of 43% compared to Planck
alone. We test how our results are potentially subject to three astrophysical
sources of contamination: sourcelens clustering, the intrinsic alignment of
galaxy shapes, and baryonic effects. We explore future limitations of the
cosmological use of thirdorder weak lensing, such as the nonlinear model and
the Gaussianity of the likelihood function.

We present the observational evidence of the 2halo term in the stacked shear
profile of a sample of about 1200 optically selected galaxy clusters based on
imaging data and the public shear catalog from the CFHTLenS. We find that the
halo bias, a measure of the correlated distribution of matter around galaxy
clusters, has amplitude and correlation with galaxy cluster mass in very good
agreement with the predictions based on the LCDM standard cosmological model.
The massconcentration relation is flat but higher than theoretical
predictions. We also confirm the close scaling relation between the optical
richness of galaxy clusters and their mass.

We present a quantitative analysis of the largest contiguous maps of
projected mass density obtained from gravitational lensing shear. We use data
from the 154 deg2 covered by the CanadaFranceHawaii Telescope Lensing Survey.
Our study is the first attempt to quantitatively characterize the scientific
value of lensing maps, which could serve in the future as a complementary
approach to the study of the dark universe with gravitational lensing. We show
that mass maps contain unique cosmological information beyond that of
traditional twopoints statistical analysis techniques. Using a series of
numerical simulations, we first show how, reproducing the CFHTLenS observing
conditions, gravitational lensing inversion provides a reliable estimate of the
projected matter distribution of large scale structure. We validate our
analysis by quantifying the robustness of the maps with various statistical
estimators. We then apply the same process to the CFHTLenS data. We find that
the 2points correlation function of the projected mass is consistent with the
cosmological analysis performed on the shear correlation function discussed in
the CFHTLenS companion papers. The maps also lead to a significant measurement
of the third order moment of the projected mass, which is in agreement with
analytic predictions, and to a marginal detection of the fourth order moment.
Tests for residual systematics are found to be consistent with zero for the
statistical estimators we used. A new approach for the comparison of the
reconstructed mass map to that predicted from the galaxy distribution reveals
the existence of giant voids in the dark matter maps as large as 3 degrees on
the sky. Our analysis shows that lensing mass maps can be used for new
techniques such as peak statistics and the morphological analysis of the
projected dark matter distribution.

We present a finelybinned tomographic weak lensing analysis of the
CanadaFranceHawaii Telescope Lensing Survey, CFHTLenS, mitigating
contamination to the signal from the presence of intrinsic galaxy alignments
via the simultaneous fit of a cosmological model and an intrinsic alignment
model. CFHTLenS spans 154 square degrees in five optical bands, with accurate
shear and photometric redshifts for a galaxy sample with a median redshift of
zm =0.70. We estimate the 21 sets of cosmic shear correlation functions
associated with six redshift bins, each spanning the angular range of
1.5<theta<35 arcmin. We combine this CFHTLenS data with auxiliary cosmological
probes: the cosmic microwave background with data from WMAP7, baryon acoustic
oscillations with data from BOSS, and a prior on the Hubble constant from the
HST distance ladder. This leads to constraints on the normalisation of the
matter power spectrum sigma_8 = 0.799 +/ 0.015 and the matter density
parameter Omega_m = 0.271 +/ 0.010 for a flat Lambda CDM cosmology. For a flat
wCDM cosmology we constrain the dark energy equation of state parameter w =
1.02 +/ 0.09. We also provide constraints for curved Lambda CDM and wCDM
cosmologies. We find the intrinsic alignment contamination to be galaxytype
dependent with a significant intrinsic alignment signal found for earlytype
galaxies, in contrast to the latetype galaxy sample for which the intrinsic
alignment signal is found to be consistent with zero.

We present cosmological constraints from 2D weak gravitational lensing by the
largescale structure in the CanadaFrance Hawaii Telescope Lensing Survey
(CFHTLenS) which spans 154 square degrees in five optical bands. Using accurate
photometric redshifts and measured shapes for 4.2 million galaxies between
redshifts of 0.2 and 1.3, we compute the 2D cosmic shear correlation function
over angular scales ranging between 0.8 and 350 arcmin. Using nonlinear models
of the darkmatter power spectrum, we constrain cosmological parameters by
exploring the parameter space with Population Monte Carlo sampling. The best
constraints from lensing alone are obtained for the smallscale
densityfluctuations amplitude sigma_8 scaled with the total matter density
Omega_m. For a flat LambdaCDM model we obtain sigma_8(Omega_m/0.27)^0.6 =
0.79+0.03.
We combine the CFHTLenS data with WMAP7, BOSS and an HST distanceladder
prior on the Hubble constant to get joint constraints. For a flat LambdaCDM
model, we find Omega_m = 0.283+0.010 and sigma_8 = 0.813+0.014. In the case
of a curved wCDM universe, we obtain Omega_m = 0.27+0.03, sigma_8 =
0.83+0.04, w_0 = 1.10+0.15 and Omega_K = 0.006+0.0060.004.
We calculate the Bayesian evidence to compare flat and curved LambdaCDM and
darkenergy CDM models. From the combination of all four probes, we find models
with curvature to be at moderately disfavoured with respect to the flat case. A
simple darkenergy model is indistinguishable from LambdaCDM. Our results
therefore do not necessitate any deviations from the standard cosmological
model.

We present the public release of the Bayesian sampling algorithm for
cosmology, CosmoPMC (Cosmology Population Monte Carlo). CosmoPMC explores the
parameter space of various cosmological probes, and also provides a robust
estimate of the Bayesian evidence. CosmoPMC is based on an adaptive importance
sampling method called Population Monte Carlo (PMC). Various cosmology
likelihood modules are implemented, and new modules can be added easily. The
importancesampling algorithm is written in C, and fully parallelised using the
Message Passing Interface (MPI). Due to very little overhead, the wallclock
time required for sampling scales approximately with the number of CPUs. The
CosmoPMC package contains postprocessing and plotting programs, and in
addition a MonteCarlo Markov chain (MCMC) algorithm. The sampling engine is
implemented in the library pmclib, and can be used independently. The software
is available for download at http://www.cosmopmc.info.

The CanadaFranceHawaii Telescope Lensing Survey (CFHTLenS) comprises deep
multicolour (u*g'r'i'z') photometry spanning 154 square degrees, with accurate
photometric redshifts and shape measurements. We demonstrate that the redshift
probability distribution function summed over galaxies provides an accurate
representation of the galaxy redshift distribution accounting for random and
catastrophic errors for galaxies with best fitting photometric redshifts z_p <
1.3.
We present cosmological constraints using tomographic weak gravitational
lensing by largescale structure. We use two broad redshift bins 0.5 < z_p <=
0.85 and 0.85 < z_p <= 1.3 free of intrinsic alignment contamination, and
measure the shear correlation function on angular scales in the range ~140
arcmin. We show that the problematic redshift scaling of the shear signal,
found in previous CFHTLS data analyses, does not afflict the CFHTLenS data. For
a flat LambdaCDM model and a fixed matter density Omega_m=0.27, we find the
normalisation of the matter power spectrum sigma_8=0.771 \pm 0.041. When
combined with cosmic microwave background data (WMAP7), baryon acoustic
oscillation data (BOSS), and a prior on the Hubble constant from the HST
distance ladder, we find that CFHTLenS improves the precision of the fully
marginalised parameter estimates by an average factor of 1.52. Combining our
results with the above cosmological probes, we find Omega_m=0.2762 \pm 0.0074
and sigma_8=0.802 \pm 0.013.

We present the CanadaFranceHawaii Telescope Lensing Survey (CFHTLenS) that
accurately determines a weak gravitational lensing signal from the full 154
square degrees of deep multicolour data obtained by the CFHT Legacy Survey.
Weak gravitational lensing by largescale structure is widely recognised as one
of the most powerful but technically challenging probes of cosmology. We
outline the CFHTLenS analysis pipeline, describing how and why every step of
the chain from the raw pixel data to the lensing shear and photometric redshift
measurement has been revised and improved compared to previous analyses of a
subset of the same data. We present a novel method to identify data which
contributes a nonnegligible contamination to our sample and quantify the
required level of calibration for the survey. Through a series of
cosmologyinsensitive tests we demonstrate the robustness of the resulting
cosmic shear signal, presenting a scienceready shear and photometric redshift
catalogue for future exploitation.

We present a comparison of Fisher matrix forecasts for cosmological probes
with Monte Carlo Markov Chain (MCMC) posterior likelihood estimation methods.
We analyse the performance of future Dark Energy Task Force (DETF) stageIII
and stage IV darkenergy surveys using supernovae, baryon acoustic
oscillations and weak lensing as probes. We concentrate in particular on the
darkenergy equation of state parameters w0 and wa. For purely geometrical
probes, and especially when marginalising over wa, we find considerable
disagreement between the two methods, since in this case the Fisher matrix can
not reproduce the highly nonelliptical shape of the likelihood function. More
quantitatively, the Fisher method underestimates the marginalized errors for
purely geometrical probes between 30%70%. For cases including structure
formation such as weak lensing, we find that the posterior probability contours
from the Fisher matrix estimation are in good agreement with the MCMC contours
and the forecasted errors only changing on the 5% level. We then explore
nonlinear transformations resulting in physicallymotivated parameters and
investigate whether these parameterisations exhibit a Gaussian behaviour. We
conclude that for the purely geometrical probes and, more generally, in cases
where it is not known whether the likelihood is close to Gaussian, the Fisher
matrix is not the appropriate tool to produce reliable forecasts.

Extending previous studies, we derive generic predictions for lower order
cumulants and their correlators for individual tomographic bins as well as
between two different bins. We derive the corresponding one and twopoint
joint probability distribution function for the tomographic convergence maps
from different bins as a function of angular smoothing scale. The modelling of
weak lensing statistics is obtained by adopting a detailed prescription for the
underlying density contrast. In this paper we concentrate on the convergence
field $\kappa$ and use tophat filter; though the techniques presented can
readily be extended to model the PDF of shear components or to include other
windows such as the compensated filter. The functional form for the underlying
PDF and bias is modelled in terms of the nonlinear or the quasilinear form
depending on the smoothing angular scale. Results from other semianalytical
models e.g. the lognormal distribution are also presented. Introducing a
reduced convergence for individual bins, we are able to show that the
tomographic PDFs and bias for each bin sample the same functional form of the
underlying PDF of density contrast but with varying variance. The joint
probability distribution of the convergence maps that correspond to two
different tomographic bins can be constructed from individual tomographic PDF
and bias. We study their dependence on cosmological parameters for source
distributions corresponding to the realistic surveys such as LSST and DES. We
briefly outline how photometric redshift information can be incorporated in our
computation of cumulants, cumulant correlators and the PDFs. Connection of our
results to the full 3D calculations is elucidated. Analytical results for
inclusion of realistic noise and finite survey size are presented in detail.