Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. The Monte Carlo Control Loops (MCCL) is a forward modelling method designed to tackle this problem. It relies on the Ultra Fast Image Generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. We present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum $C_{\ell}$ and the redshift distribution $n(z)$ of the background galaxy sample. The method includes maps of the systematic sources, Point Spread Function (PSF), an Approximate Bayesian Computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and the fast estimation of the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a non-tomographic setting, we derive a set of $C_\ell$ and $n(z)$ curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of $\Omega_m$, $\sigma_8$, and intrinsic alignment amplitude $A_{\rm{IA}}$, defined as $S_8D_{\rm{IA}} = \sigma_8(\Omega_m/0.3)^{0.5}D_{\rm{IA}}$, where $D_{\rm{IA}}=1-0.11(A_{\rm{IA}}-1)$. We find $S_8D_{\rm{IA}}=0.895^{+0.054}_{-0.039}$, where systematics are at the level of roughly 60\% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys.

We present a structural and morphological catalogue for 45 million objects selected from the first year of data from the Dark Energy Survey (DES). Single Sersic fits and non-parametric measurements are produced for g, r and i filters. The parameters from the best-fitting Sersic model (total magnitude, half-light radius, Sersic index, axis ratio and position angle) are measured with Galfit; the non-parametric coefficients (concentration, asymmetry, clumpiness, Gini, M20) are provided using the Zurich Estimator of Structural Types (ZEST+). To study the statistical uncertainties, we consider a sample of state-of-the-art image simulations with a realistic distribution in the input parameter space and then process and analyse them as we do with real data: this enables us to quantify the observational biases due to PSF blurring and magnitude effects and correct the measurements as a function of magnitude, galaxy size, Sersic index (concentration for the analysis of the non-parametric measurements) and ellipticity. We present the largest structural catalogue to date: we find that accurate and complete measurements for all the structural parameters are typically obtained for galaxies with SExtractor MAG AUTO I < 21. Indeed, the parameters in the filters i and r can be overall well recovered up to MAG AUTO < 21.5, corresponding to a fitting completeness of ~90% below this threshold, for a total of 25 million galaxies. The combination of parametric and non-parametric structural measurements makes this catalogue an important instrument to explore and understand how galaxies form and evolve. The catalogue described in this paper will be publicly released alongside the Dark Energy Survey collaboration Y1 cosmology data products at the following URL: https://des.ncsa.illinois.edu/releases/y1a1/gold/morphology.

Determining the redshift distribution $n(z)$ of galaxy samples is essential for several cosmological probes including weak lensing. For imaging surveys, this is usually done using photometric redshifts estimated on an object-by-object basis. We present a new approach for directly measuring the global $n(z)$ of cosmological galaxy samples, including uncertainties, using forward modeling. Our method relies on image simulations produced using UFig (Ultra Fast Image Generator) and on ABC (Approximate Bayesian Computation) within the $MCCL$ (Monte-Carlo Control Loops) framework. The galaxy population is modeled using parametric forms for the luminosity functions, spectral energy distributions, sizes and radial profiles of both blue and red galaxies. We apply exactly the same analysis to the real data and to the simulated images, which also include instrumental and observational effects. By adjusting the parameters of the simulations, we derive a set of acceptable models that are statistically consistent with the data. We then apply the same cuts to the simulations that were used to construct the target galaxy sample in the real data. The redshifts of the galaxies in the resulting simulated samples yield a set of $n(z)$ distributions for the acceptable models. We demonstrate the method by determining $n(z)$ for a cosmic shear like galaxy sample from the 4-band Subaru Suprime-Cam data in the COSMOS field. We also complement this imaging data with a spectroscopic calibration sample from the VVDS survey. We compare our resulting posterior $n(z)$ distributions to the one derived from photometric redshifts estimated using 36 photometric bands in COSMOS and find good agreement. This offers good prospects for applying our approach to current and future large imaging surveys.

Weak Gravitational Lensing is a powerful probe of the dark sector of the Universe. One of the main challenges for this technique is the treatment of systematics in the measurement of cosmic shear from galaxy shapes. In an earlier work, Refregier & Amara (2014) have proposed the Monte Carlo Control Loops (MCCL) to overcome these effects using a forward modeling approach. We focus here on one of the control loops in this method, the task of which is the calibration of the shear measurement. For this purpose, we first consider the requirements on the shear systematics for a given survey and propagate them to different systematics terms. We use two one-point statistics to calibrate the shear measurement and six further one-point statistics as diagnostics. We also propagate the systematics levels that we estimate from the one-point functions to the two-point functions for the different systematic error sources. This allows us to assess the consistency between the systematics levels measured in different ways. To test the method, we construct synthetic sky surveys with an area of 1,700 deg$^2$. With some simplifying assumptions, we are able to meet the requirements on the shear calibration for this survey configuration. Furthermore, we account for the total residual shear systematics in terms of the contributing sources. We discuss how this MCCL framework can be applied to current and future weak lensing surveys.

We present the first constraints on cosmology from the Dark Energy Survey (DES), using weak lensing measurements from the preliminary Science Verification (SV) data. We use 139 square degrees of SV data, which is less than 3\% of the full DES survey area. Using cosmic shear 2-point measurements over three redshift bins we find $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5} = 0.81 \pm 0.06$ (68\% confidence), after marginalising over 7 systematics parameters and 3 other cosmological parameters. We examine the robustness of our results to the choice of data vector and systematics assumed, and find them to be stable. About $20$\% of our error bar comes from marginalising over shear and photometric redshift calibration uncertainties. The current state-of-the-art cosmic shear measurements from CFHTLenS are mildly discrepant with the cosmological constraints from Planck CMB data; our results are consistent with both datasets. Our uncertainties are $\sim$30\% larger than those from CFHTLenS when we carry out a comparable analysis of the two datasets, which we attribute largely to the lower number density of our shear catalogue. We investigate constraints on dark energy and find that, with this small fraction of the full survey, the DES SV constraints make negligible impact on the Planck constraints. The moderate disagreement between the CFHTLenS and Planck values of $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ is present regardless of the value of $w$.

Cosmic shear is sensitive to fluctuations in the cosmological matter density field, including on small physical scales, where matter clustering is affected by baryonic physics in galaxies and galaxy clusters, such as star formation, supernovae feedback and AGN feedback. While muddying any cosmological information that is contained in small scale cosmic shear measurements, this does mean that cosmic shear has the potential to constrain baryonic physics and galaxy formation. We perform an analysis of the Dark Energy Survey (DES) Science Verification (SV) cosmic shear measurements, now extended to smaller scales, and using the Mead et al. 2015 halo model to account for baryonic feedback. While the SV data has limited statistical power, we demonstrate using a simulated likelihood analysis that the final DES data will have the statistical power to differentiate among baryonic feedback scenarios. We also explore some of the difficulties in interpreting the small scales in cosmic shear measurements, presenting estimates of the size of several other systematic effects that make inference from small scales difficult, including uncertainty in the modelling of intrinsic alignment on nonlinear scales, `lensing bias', and shape measurement selection effects. For the latter two, we make use of novel image simulations. While future cosmic shear datasets have the statistical power to constrain baryonic feedback scenarios, there are several systematic effects that require improved treatments, in order to make robust conclusions about baryonic feedback.

We present a comprehensive method for determining stellar mass functions, and apply it to samples in the local Universe. We combine the classical 1/Vmax approach with STY, a parametric maximum likelihood method and SWML, a non-parametric maximum likelihood technique. In the parametric approach, we are assuming that the stellar mass function can be modelled by either a single or a double Schechter function and we use a likelihood ratio test to determine which model provides a better fit to the data. We discuss how the stellar mass completeness as a function of z biases the three estimators and how it can affect, especially the low mass end of the stellar mass function. We apply our method to SDSS DR7 data in the redshift range from 0.02 to 0.06. We find that the entire galaxy sample is best described by a double Schechter function with the following parameters: $\log (M^{*}/M_\odot) = 10.79 \pm 0.01$, $\log (\Phi^{*}_1/\mathrm{h^3\ Mpc^{-3}}) = -3.31 \pm 0.20$, $\alpha_1 = -1.69 \pm 0.10$, $\log (\Phi^{*}_2/\mathrm{h^3\ Mpc^{-3}}) = -2.01 \pm 0.28$ and $\alpha_2 = -0.79 \pm 0.04$. We also use morphological classifications from Galaxy Zoo and halo mass, overdensity, central/satellite, colour and sSFR measurements to split the galaxy sample into over 130 subsamples. We determine and present the stellar mass functions and the best fit Schechter function parameters for each of these subsamples.

We measure the projected density profile, shape and alignment of the stellar and dark matter mass distribution in 11 strong-lens galaxies. We find that the projected dark matter density profile - under the assumption of a Chabrier stellar initial mass function - shows significant variation from galaxy to galaxy. Those with an outermost image beyond $\sim 10$ kpc are very well fit by a projected NFW profile; those with images within 10 kpc appear to be more concentrated than NFW, as expected if their dark haloes contract due to baryonic cooling. We find that over several half-light radii, the dark matter haloes of these lenses are rounder than their stellar mass distributions. While the haloes are never more elliptical than $e_{dm} = 0.2$, their stars can extend to $e_* > 0.2$. Galaxies with high dark matter ellipticity and weak external shear show strong alignment between light and dark; those with strong shear ($\gamma \gtrsim 0.1$) can be highly misaligned. This is reassuring since isolated misaligned galaxies are expected to be unstable. Our results provide a new constraint on galaxy formation models. For a given cosmology, these must explain the origin of both very round dark matter haloes and misaligned strong-lens systems.

We present photometric redshift estimates for galaxies used in the weak lensing analysis of the Dark Energy Survey Science Verification (DES SV) data. Four model- or machine learning-based photometric redshift methods -- ANNZ2, BPZ calibrated against BCC-Ufig simulations, SkyNet, and TPZ -- are analysed. For training, calibration, and testing of these methods, we construct a catalogue of spectroscopically confirmed galaxies matched against DES SV data. The performance of the methods is evaluated against the matched spectroscopic catalogue, focusing on metrics relevant for weak lensing analyses, with additional validation against COSMOS photo-zs. From the galaxies in the DES SV shear catalogue, which have mean redshift $0.72\pm0.01$ over the range $0.3<z<1.3$, we construct three tomographic bins with means of $z=\{0.45, 0.67, 1.00\}$. These bins each have systematic uncertainties $\delta z \lesssim 0.05$ in the mean of the fiducial SkyNet photo-z $n(z)$. We propagate the errors in the redshift distributions through to their impact on cosmological parameters estimated with cosmic shear, and find that they cause shifts in the value of $\sigma_8$ of approx. 3%. This shift is within the one sigma statistical errors on $\sigma_8$ for the DES SV shear catalog. We further study the potential impact of systematic differences on the critical surface density, $\Sigma_{\mathrm{crit}}$, finding levels of bias safely less than the statistical power of DES SV data. We recommend a final Gaussian prior for the photo-z bias in the mean of $n(z)$ of width $0.05$ for each of the three tomographic bins, and show that this is a sufficient bias model for the corresponding cosmology analysis.

Spatially-varying depth and characteristics of observing conditions, such as seeing, airmass, or sky background, are major sources of systematic uncertainties in modern galaxy survey analyses, in particular in deep multi-epoch surveys. We present a framework to extract and project these sources of systematics onto the sky, and apply it to the Dark Energy Survey (DES) to map the observing conditions of the Science Verification (SV) data. The resulting distributions and maps of sources of systematics are used in several analyses of DES SV to perform detailed null tests with the data, and also to incorporate systematics in survey simulations. We illustrate the complementarity of these two approaches by comparing the SV data with the BCC-UFig, a synthetic sky catalogue generated by forward-modelling of the DES SV images. We analyse the BCC-UFig simulation to construct galaxy samples mimicking those used in SV galaxy clustering studies. We show that the spatially-varying survey depth imprinted in the observed galaxy densities and the redshift distributions of the SV data are successfully reproduced by the simulation and well-captured by the maps of observing conditions. The combined use of the maps, the SV data and the BCC-UFig simulation allows us to quantify the impact of spatial systematics on $N(z)$, the redshift distributions inferred using photometric redshifts. We conclude that spatial systematics in the SV data are mainly due to seeing fluctuations and are under control in current clustering and weak lensing analyses. The framework presented here is relevant to all multi-epoch surveys, and will be essential for exploiting future surveys such as the Large Synoptic Survey Telescope (LSST), which will require detailed null-tests and realistic end-to-end image simulations to correctly interpret the deep, high-cadence observations of the sky.

Weak lensing by large-scale structure is a powerful technique to probe the dark components of the universe. To understand the measurement process of weak lensing and the associated systematic effects, image simulations are becoming increasingly important. For this purpose we present a first implementation of the $\textit{Monte Carlo Control Loops}$ ($\textit{MCCL}$; Refregier & Amara 2014), a coherent framework for studying systematic effects in weak lensing. It allows us to model and calibrate the shear measurement process using image simulations from the Ultra Fast Image Generator (UFig; Berge et al. 2013). We apply this framework to a subset of the data taken during the Science Verification period (SV) of the Dark Energy Survey (DES). We calibrate the UFig simulations to be statistically consistent with DES images. We then perform tolerance analyses by perturbing the simulation parameters and study their impact on the shear measurement at the one-point level. This allows us to determine the relative importance of different input parameters to the simulations. For spatially constant systematic errors and six simulation parameters, the calibration of the simulation reaches the weak lensing precision needed for the DES SV survey area. Furthermore, we find a sensitivity of the shear measurement to the intrinsic ellipticity distribution, and an interplay between the magnitude-size and the pixel value diagnostics in constraining the noise model. This work is the first application of the $\textit{MCCL}$ framework to data and shows how it can be used to methodically study the impact of systematics on the cosmic shear measurement.

We present a forward-modelling simulation framework designed to model the data products from the Dark Energy Survey (DES). This forward-model process can be thought of as a transfer function -- a mapping from cosmological and astronomical signals to the final data products used by the scientists. Using output from the cosmological simulations (the Blind Cosmology Challenge), we generate simulated images (the Ultra Fast Image Simulator, Berge et al. 2013) and catalogs representative of the DES data. In this work we simulate the 244 sq. deg coadd images and catalogs in 5 bands for the DES Science Verification (SV) data. The simulation output is compared with the corresponding data to show that major characteristics of the images and catalogs can be captured. We also point out several directions of future improvements. Two practical examples, star/galaxy classification and proximity effects on object detection, are then used to demonstrate how one can use the simulations to address systematics issues in data analysis. With clear understanding of the simplifications in our model, we show that one can use the simulations side-by-side with data products to interpret the measurements. This forward modelling approach is generally applicable for other upcoming and future surveys. It provides a powerful tool for systematics studies which is sufficiently realistic and highly controllable.