Issue 
A&A
Volume 542, June 2012



Article Number  A5  
Number of page(s)  31  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201117625  
Published online  24 May 2012 
Galaxy clustering in the CFHTLSWide: the changing relationship between galaxies and haloes since z ~ 1.2^{⋆,}^{⋆⋆}
^{1} Astronomical Institute, Graduate School of Science, Tohoku University, Sendai 9808578, Japan
email: jean.coupon@gmail.com
^{2} Institute of Astronomy and Astrophysics, Academia Sinica, PO Box 23141, Taipei 10617, Taiwan
^{3} Excellence Cluster Universe, Technische Universität München, Boltzmannstr. 2, 85748 Garching, Germany
^{4} UniversitätsSternwarte München, Scheinerstr. 1, 81679 München, Germany
^{5} Institut d’Astrophysique de Paris, UMR 7095 CNRS, Université Pierre et Marie Curie, 98bis boulevard Arago, 75014 Paris, France
^{6} Laboratoire d’Astrophysique de Marseille (UMR 6110), CNRSUniversité de Provence, 38 rue Frédéric JoliotCurie, 13388 Marseille Cedex 13, France
^{7} CanadaFranceHawaii Telescope, 65–1238 Mamalahoa Highway, Kamuela, HI 9674, USA
^{8} INAF – Osservatorio Astronomico di Torino, 10025 Pino Torinese, Italy
^{9} SUPA, Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK
Received: 4 July 2011
Accepted: 28 February 2012
It has become increasingly apparent that studying how dark matter haloes are populated by galaxies can provide new insights into galaxy formation and evolution. In this paper, we present a detailed investigation of the changing relationship between galaxies and the dark matter haloes they inhabit from z ~ 1.2 to the present day. We do this by comparing precise galaxy clustering measurements over 133 deg^{2} of the “Wide” component of the CanadaFranceHawaii Telescope Legacy Survey (CFHTLS) with predictions of an analytic halo occupation distribution (HOD) model where the number of galaxies in each halo depends only on the halo mass. Starting from a parent catalogue of ~3 × 10^{6} galaxies at i′_{AB}<22.5 we use accurate photometric redshifts calibrated using ~10^{4} spectroscopic redshifts to create a series of typeselected volumelimited samples covering 0.2 < z < 1.2. Our principal result, based on clustering measurements in these samples, is a robust determination of the luminositytohalo mass ratio and its dependence on redshift and galaxy type. For the full sample, this reaches a peak at low redshifts of M_{h}^{peak} = 4.5×10^{11} h^{1} M_{⊙} and moves towards higher halo masses at higher redshifts. For redder galaxies the peak is at higher halo masses and does not evolve significantly over the entire redshift range of our survey. We also consider the evolution of bias, average halo mass and the fraction of satellites as a function of redshift and luminosity. Our observed growth of a factor of ~2 in satellite fraction between z ~ 1 and z ~ 0 is testament to the limited role that galaxy merging plays in galaxy evolution for ~10^{12} h^{1} M_{⊙} mass haloes at z < 1. Qualitatively, our observations are consistent with a picture in which red galaxies in massive haloes have already accumulated most of their stellar mass by z ~ 1 and subsequently undergo little evolution until the present day. The observed movement of the peak location for the full galaxy population is consistent with the bulk of starformation activity migrating from higher mass haloes at high redshifts to lower mass haloes at lower redshifts.
Key words: cosmology: observations / largescale structure of Universe / galaxies: distances and redshifts / galaxies: halos / galaxies: evolution
Based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, at the CanadaFranceHawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des 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 TERAPIX and the Canadian Astronomy Data Centre as part of the CanadaFranceHawaii Telescope Legacy Survey, a collaborative project of NRC and CNRS.
Appendices are available in electronic form at http://www.aanda.org
© ESO, 2012
1. Introduction
In our current paradigm of galaxy formation, haloes of dark matter grow from tiny imperfections in the early Universe. Against the background (at the present day) of an accelerating Universe, structures grow “hierarchically”: small haloes form first and merge to build up larger ones, resulting in a complex filamentary network where the most massive haloes lie at the nodes of the cosmic web (Springel et al. 2006). Galaxies form as baryons fall into the centres of dark matter haloes and the inflow of cold gas provides the fuel for star formation (White & Rees 1978; Fall & Efstathiou 1980). While this scenario broadly reproduces the observed galaxy distribution over a large range of scales and cosmic time, several open questions remain concerning which processes regulate star formation inside haloes and how this leads to the Hubble sequence observed at the present day.
These physical processes which drive galaxy formation and regulate star formation are expected to correlate closely to the mass of dark matter haloes which host galaxies. Therefore, a profitable avenue to pursue to better understand the physics of galaxy formation is to determine how the relationship between the stellar mass M_{star} and the dark matter halo mass M_{h} changes over the lifetime of the Universe. In the local Universe, abundancematching studies have shown that at z < 0.1 the stellartohalo mass ratio (SHMR, M_{star}/M_{h}) reaches a maximum of a few percent at M_{h} ~ 10^{12} M_{⊙}, much smaller than the universal baryon fraction (Guo et al. 2010). This indicates that the conversion efficiency of baryons to luminous galaxies is highest at these halo masses. In less massive haloes, supernovae winds are responsible for gas reheating, which reduces star formation and flattens the galaxy luminosity function at the faint end (Benson et al. 2003). In more massive haloes, feedback from active galactic nuclei (Hopkins et al. 2006; Somerville et al. 2008) can “quench” star formation, leading to formation of the “red sequence” of passively evolving galaxies (Bower et al. 2006; Croton et al. 2006). The relative importance of these different processes as well as their dependence on environment, redshift and galaxy type is still very much an open question.
Ideally we would like to make a detailed comparison between the luminous content of dark matter haloes and the halo masses themselves as a function of redshift. Several methodologies exist to estimate halo masses: for example, using satellite kinematics (More et al. 2011), Xray temperatures (Peterson & Fabian 2006; Cattaneo et al. 2009), or gravitational lensing (Heymans et al. 2006; Mandelbaum et al. 2006; Auger et al. 2010; Leauthaud et al. 2012b). However, the range in redshifts over which these techniques can be applied is limited, and amassing large samples is not always easy.
At the price of making some assumptions concerning the profiles of dark matter haloes and their evolution in number density, another way of inferring halo masses is from the observed galaxy clustering and mass functions. Usually, one makes the simplifying assumption that the number of galaxies in a given dark matter halo only depends on the halo mass, an assumption which has been verified with Nbody simulations (Berlind et al. 2003, 2005; Berlind & Weinberg 2002; Kravtsov et al. 2004). Furthermore, the galaxy population of a halo, described by the “halo occupation distribution” (HOD), is assumed to be independent of environment and assembly history of the halo. Although some authors have shown that this “halo assembly bias” can have an observable effect on galaxy clustering measurements, the magnitude of this effect is, for the moment, much smaller than the size of systematic errors (Croton et al. 2007). This will be an important effect to test in future surveys.
In this paper we will use this model to investigate the changing relationship of dark matter and luminous matter from z ~ 1.2 to the present day. Until now, the majority of papers which have used HOD modelling to interpret galaxy clustering have analysed either large, low redshift surveys such as the SDSS (see, for e.g., Zehavi et al. 2011; Blake et al. 2008), or smaller (~1 deg^{2}) deep fields such as the COMBO17 and COSMOS surveys (Phleps et al. 2006; Simon et al. 2009; Leauthaud et al. 2012b). For the moment, spectroscopic surveys at higher redshifts are also limited to small fields of view (Zheng et al. 2007; Abbas et al. 2010). Brown et al. (2008) has interpreted clustering measurements in the framework of the halo model of slightly larger redshift baseline, although their survey only covered ~ 7 deg^{2} and used fiveband photometric redshifts. None of these surveys had the required combination of depth and areal coverage to make reliable measurements from low redshift to z ~ 1 of large samples of galaxies selected by type and intrinsic luminosity. In contrast, this is an ideal task for the CanadaFranceHawaii Legacy Survey (CFHTLS) with its unique combination of depth, area, and image quality.
Fig. 1 Sky coverage of W1, W2, W3 and W4 fields. Each blue square represents one MegaCam pointing. Spectroscopic data (VVDS F02, F22, and DEEP2) are in red. CFHTLS Deep fields overlapping with the wide survey (D1 and D3) are shown in green. 
The fiveband u^{∗},g′,r′,i′,z′ CFHTLS photometry allows us to measure photometric redshifts out to z ~ 1.2 with ~3% accuracy and small systematic errors (Coupon et al. 2009). The large sky coverage, ~133 deg^{2} after masking, enables us to probe a large range of halo masses. In addition, from the four independent fields of the CFHTLS we can obtain a reliable estimate of cosmic variance from the data itself. In this paper we construct a series of volumelimited samples selected by type, luminosity, and redshift, and use the halo model combined with the HOD formalism to infer how each galaxy sample populates their hosting dark matter haloes, and how this changes with lookback time.
The paper is organized as follows. In Sect. 2, we present our data set and catalogue production, including data reduction and photometric redshift measurement. In Sect. 3 we outline the techniques we use to estimate galaxy clustering, and in Sect. 4 we present our model. Results are described in Sect. 5 and discussed in Sect. 6. Conclusions are presented in Sect. 7.
Throughout the paper we use a flat ΛCDM cosmology (Ω_{m} = 0.25, Ω_{Λ} = 0.75, H_{0} = 70 km s^{1} Mpc^{1} and σ_{8} = 0.8) with h = H_{0}/100 km s^{1} Mpc^{1}. All magnitudes are given in the AB system.
2. Observations, reductions, and catalogue production
2.1. The CanadaFranceHawaii Legacy Survey
We use the “T0006” release of the CFHTLS Wide^{1} (Goranova et al. 2009). The CFHTLS Wide has the mean limiting AB magnitudes (measured as the 50% completeness for point sources) ~25.3,25.5,24.8,24.48,23.60 in u^{∗},g′,r′,i′,z′, respectively. An important feature in T0006 compared to previous releases is that for each MegaCam pointing (“tile”), the TERAPIX group provides a magnitude offset in each filter to account for tiletotile variations in magnitude zero points. These offsets are calculated using a stellar locus regression technique (High et al. 2009). In each tile, objects are detected on a griχ^{2} image (Szalay et al. 1999). Galaxies are selected using SEXtractor “mag_auto” magnitudes (Bertin & Arnouts 1996), and colours are measured in 3″ apertures which are all well matched to the median size of galaxies at our i′ < 22.5 magnitude limit. The sky coverage of the four CFHTLS Wide fields is shown in Fig. 1.
The MegaCam iband filter broke in 2006, and subsequent observations were made with a new iband filter denoted by “y”. Data from these two filters were treated separately with separate filter curves. Using a small data set observed with both the new y and the old iband, we detected no significant difference in photometric redshift accuracy. For the rest of the paper, we use i′ to represent both filters.
2.2. Photometric redshift estimation
Our photometric redshifts were measured using LePhare (Arnouts et al. 2002) following the procedures outlined in Ilbert et al. (2006) and Coupon et al. (2009). Our primary template set is composed of the four Coleman et al. (1980) (CWW) observed spectra (Ell, Sbc, Scd, Irr) complemented with two observed starburst spectral energy distributions (SED) from Kinney et al. (1996). This primary set of templates is optimized using the VVDS Deep spectroscopic sample. We performed an automatic calibration of the zeropoints using spectroscopic redshifts in the W1, W3 and W4 fields. The calibration is obtained by comparing the observed and modelled fluxes, and done iteratively until the zeropoints converge. For spectral types later than Sbc, we introduced a reddening E(B − V) = 0 to 0.35 using the Calzetti et al. (2000) extinction law.
Without nearinfrared photometry, Lyman and Balmer breaks create degeneracies between high and low photometric redshift solutions. Following a Bayesian approach similar to Benítez (2000) and Ilbert et al. (2006), we adopt a “prior” based on the observed spectroscopic redshift distribution to reduce the number of catastrophic failures. No redshift solution is allowed which will produce a galaxy brighter than M_{g} = −24.
Our photometric redshifts were calibrated with spectroscopic samples. The VVDS Deep survey (Le Fèvre et al. 2005) is available in the W1 field. The VVDS Deep is a pure magnitudelimited sample at I_{AB} < 24 and we used 5926 secure redshifts (flags 4 and 5) for photometric calibration, and for building the redshift distribution prior. From the VVDS F22 (“VVDS WIDE”, Garilli et al. 2008) in W4 we used 4514 galaxies limited to I_{AB} < 22.5. This sample has a success rate of 92%. Finally, in the W3 field, 1267 redshifts from the DEEP2 survey (Davis et al. 2003) selected at R_{AB} < 24.1 were used for photometric calibration. No spectroscopic redshifts were available in W2.
To estimate photometric redshifts, we adopted a slightly different method than for previous studies to measure redshifts from the probability distribution function (PDF). Rather than directly using the redshift that minimizes the χ^{2} distribution function, we chose the median of the PDF. It has been found that the maximum likelihood estimate can lead to a systematic concentration of photometric redshifts around discrete values (“redshift focusing”), although we checked that the effect is very small for redshift bins of size Δz > 0.1.
The redshift accuracy (“σ”) of our samples with spectroscopic redshift z_{s} and photometric redshift z_{p} is defined as σ_{Δz}/(1 + z_{s}), where Δz = z_{p} − z_{s} represents the difference between spectroscopic and photometric redshift. We use the normalised median absolute deviation as accuracy estimator (NMAD, Hoaglin et al. 1983), expressed as σ_{Δz} = 1.48 × median(z_{p} − z_{s}). “Catastrophic” redshifts are defined as objects with z_{p} − z_{s}/(1 + z_{s}) > 0.15. The percentage of catastrophic redshifts is denoted by η. The dispersion and the fraction of catastrophic redshifts increase from bright to faint samples. The accuracy is σ_{Δz}/(1 + z_{s}) = 0.034 (W1) to 0.037 (W4) at i′ < 21.5 and reaches 0.039 (W1) to 0.051 (W3) at 21.5 < i′ < 22.5. The failure rate increases from η ~ 1.4% at i′ < 21.5 to η ~ 5.5% at 21.5 < i′ < 22.5. Galaxies with redder colours (El, Sbc) have more accurate photometric redshift estimates (σ_{Δz}/(1 + z_{s}) = 0.031, η = 1.64%) than bluer ones (Scd, Irr and starbursts, σ_{Δz}/(1 + z_{s}) = 0.038, η = 3.69%).
We constructed a merged catalogue with unique objects from individual tile catalogues provided by TERAPIX. We identified duplicate sources positionally and chose the object closest to the centre of its original tile. Although photometric redshifts were estimated to i′ < 24, we selected objects brighter than i′ = 22.5 to maintain a low outlier rate (see Table 3 of Coupon et al. 2009) and to ensure catalogue completeness. This sample covers a total effective area of 133 deg^{2}. We show the photometric redshift accuracy in the secure redshift range covered in this study, 0.2 < z < 1.2, compared to spectroscopic redshifts in Fig. 2.
Stars appear as pointlike sources on the image and can easily be identified using their halflight radius at bright magnitudes; however, at fainter magnitudes this estimator can be confused by unresolved galaxies. We refined our stargalaxy separation by combining this profilemagnitude criteria with a colourbased criteria. We computed the halflight radius limit for each tile from the median and dispersion of the brightest objects and classified as stars all objects brighter than i′ = 20 with flux radii limits smaller than this amount. In the magnitude range 20 < i′ < 22.5, we classified stars as objects smaller than the halflight radius and with χ^{2}(star) < χ^{2}(galaxy), where χ^{2}(object) is estimated from the stellar and galaxy template libraries.
We eliminated objects observed with less than three filters and with χ^{2}(galaxy) > 100. We suspect that these objects are false detections or rare objects with wrong photometric redshift estimates. Most of these objects reside in specific areas, such as around bright stars and field borders where the masking procedure has failed.
Using the photometric redshift, the associated bestfitting template and the observed apparent magnitude in a given band, we can directly measure the kcorrection and the absolute magnitude in any restframe band. Since at high redshifts the kcorrection depends strongly on the galaxy SED, it is the main source of systematic error in determining absolute magnitudes. To minimise kcorrection uncertainties, we derive the restframe luminosity at λ using the object’s apparent magnitude closer to λ × (1 + z). We use either the u∗, g′, r′, i′, or z′ observed apparent magnitudes according to the redshift of the galaxy (the procedure is described in Appendix A of Ilbert et al. 2005). For this reason the bluest absolute magnitude estimate makes full advantage of the complete observed magnitude set. However, as the uband flux has larger photometric errors, we decided to use M_{g}band magnitudes.
Fig. 2 Photometric redshift accuracy in the 0.2 < photo  z < 1.2 interval. The figure shows the number density of photometric redshifts in W1, W3 and W4, versus spectroscopic redshifts from VVDS and DEEP2 surveys. The blue line is a linear fit to photoz versus spectroz. 
Fig. 3 Top: type selection based on best fitting templates for the “red” sample (left) and the “blue” sample (right) in the CFHTLS Wide. We show the colour distribution (M_{u} − M_{r}) as function of absolute magnitude (M_{g}) and redshift (top to bottom) and number counts for “red” and “blue” objects on the right panels. Bottom: number of red and blue galaxies as function of magnitude and redshift. 
Fig. 4 Sample selection in the full (left), red (centre) and blue (right) galaxy samples. In each panel, the galaxy number density in the plane M_{g}/z is shown; red rectangles represent the luminosity threshold samples. 
3. Clustering measurements
3.1. Sample selection and galaxy number density estimation
Objects are selected with i′ < 22.5 and divided into five redshift bins: 0.2 < z < 0.4, 0.4 < z < 0.6, 0.6 < z < 0.8, 0.8 < z < 1.0, and 1.0 < z < 1.2. A large bin width (Δ_{z} = 0.2) ensures a low bintobin contamination from random errors (σ_{z} < 0.1). We investigate effects from systematic errors in Sect. 3.2. The complete sample contains 2 924 730 galaxies. To separate galaxies into “red sequence” and “blue cloud” types, we used a criterion based on bestfitting galaxy templates. We define “red” galaxies as objects with bestfitting CWW templates (see Sect. 2.2) estimated as El and Sbc (earlytype), and “blue” galaxies estimated as Sbc, Scd, Im, SB1 and SB2 (latetype). Selecting galaxies by bestfitting template is a more robust way of selecting galaxies by type, as compared to a simple colour selection, and it is less sensitive to the effects of dust extinction. Our resulting colour selection in the absolute magnituderedshift plane is shown in Fig. 3. We note that a small amount of redclassified galaxies at high redshift z > 0.8 have slightly bluer colours. Their distribution peaks near the blue galaxy distribution, suggesting that these objects could be blue galaxies erroneously identified as red galaxies due to photometric redshift errors. We also note that a simple colour cut would not exactly reproduce our selection. However, in the interests of simplicity and clarity we keep the “red”/”blue” labels for the rest of the paper.
We extract volumelimited luminosityselected samples for each of the “full” (or “all galaxies”), “red” and “blue” samples, using M_{g} absolute magnitude thresholds (hereafter denoted as “luminosity threshold samples”), from M_{g} − 5log h = −17.8 (fainter threshold in the range 0.2 < z < 0.4) to M_{g} − 5log h = −22.8 (brighter threshold in the range 1.0 < z < 1.2). The sample selection is illustrated in Fig. 4. For the rest of this paper we will refer to these samples as simply full, red and blue. Due to low numbers of pairs at small scales, luminous blue samples were discarded. We are left with 45 samples, each comprising on average ~153 000, ~70 000 and ~129 000 galaxies for a typical full, red and blue sample, respectively. The sample properties are displayed in Tables B.1–B.3.
Finally, in each redshift interval [z_{min};z_{max}] we compute the galaxy number density: (1)where Ω represents the solid angle subtended by the survey, and dV/dz the volume element. Errors are estimated from the weighted galaxy number density fieldtofield variance.
3.2. Photometric redshift uncertainties
Our modelled twopoint correlation function is projected using the measured redshift distributions. In order to take into account statistical errors on redshifts, we select galaxies in the redshift range considered and convolve the observed redshift distributions with the estimated photometric redshift errors, derived from the probability distribution functions. We construct a Gaussian error distribution for each galaxy centred on the median redshift of the PDF, with a width corresponding to the 68% confidence limits of the PDF, and normalised to unity. We then sum these Gaussians to construct the redshift distribution resampled to a redshift bin width of 0.04. Redshift distributions for each sample are illustrated in Fig. 5.
Fig. 5 Redshift distributions for the full (left), red (centre) and blue (right) galaxy samples as function of absolute magnitude threshold. 
To further assess the quality of photometric redshifts, we perform the crosscorrelation analysis introduced in Benjamin et al. (2010). The measurement of the angular correlation functions for galaxies in different photoz bins is used to constrain the fraction of galaxies that are scattered into “wrong” redshift bins due to photoz errors. We measure the bintobin cross correlation function for the full, red and blue samples, respectively. A nonzero correlation between adjacent redshift bins is present in all cases. This may be due to the presence of largescale structures extending over several redshift bins. More importantly, this is also due to photometric redshift scatter, which results in the leakage of galaxies into neighbouring bins. This error contribution to the redshifts is taken into account by the convolution of the redshift distribution with the errors, therefore, we do not consider adjacent bins further in the analysis of photometric redshifts uncertainties.
The angular crosscorrelation of galaxies in nonadjacent redshift bins is much lower, indicating a small fraction of catastrophic outliers. We use the “global pairwise analysis” method to measure the contamination between two redshift bins i and j. In this approximation, the following linear combinations of the angular crosscorrelation function w_{ij} and the two autocorrelation functions, w_{ii} and w_{jj}, respectively, are expected to cancel for all angular scales θ_{t}, (2)Here, N_{i}(N_{j}) is the observed number of galaxies in bin i(j). The contamination f_{ij} is the number of galaxies with true redshift in bin i, but misidentified into bin j, as a fraction of the true number of galaxies in bin i. For each bin pair (i,j), the leakage of the other redshift bins is neglected. This approximation is valid for contamination fractions of up to 10% (Benjamin et al. 2010). With this, the fraction of galaxies which stay in their bin is f_{ii} = 1 − f_{ij}. We fit the two parameters f_{ij} and f_{ji} in Eq. (2) by performing a χ^{2} nulltest on d_{t}. For the covariance ⟨ d_{t}d_{s} ⟩ , we take into account the correlation between angular scales for each of the three correlation functions using a Jackknife estimate (see next section). We neglect the subdominant covariance between different correlation functions. This corresponds to using the first three terms in Eq. (A4) of Benjamin et al. (2010).
Due to degeneracies between the parameters f_{ij} and f_{ji}, large values for f_{ij}(i > j) cannot be ruled out in principle. However, for the full and red galaxy samples, the contamination fractions are consistent with zero in most cases. The blue galaxy samples are slightly worse, but contaminations are consistent with values between 2% and 10%. Together with the very low outlier rate for the subsample with spectroscopic redshifts (see Fig. 2), these results further strengthen our confidence in the photometric redshift estimates in this work and to use them to measure angular correlation functions. Our measurements for the full sample are shown in Fig. 6.
Bintobin mixing caused by photometric redshift errors will reduce our measured clustering signal, given that objects in separate bins are uncorrelated. This would lead to an underestimation of halo masses, as low mass haloes tend to be less clustered than more massive ones. However, it is not trivial to predict how galaxy number density estimates would be affected by such errors, which may result in incomplete or contaminated samples. Since the halo number density decreases monotonically with halo mass, galaxy number density decreases with increasing host halo mass. Therefore, a contaminated sample will underestimate halo mass fitting estimates. Conversely, an incomplete sample will tend to overestimate halo masses.
To understand how photometric redshift errors could affect our results as function of redshift and luminosity threshold, we compared our measurements in the CFHTLS Wide with those from the “Deep” component of the CFHTLS (Goranova et al. 2009), where more precise photometric redshifts are available (see Table 3 in Coupon et al. 2009). Photometric redshifts in the CFHTLS Deep were computed using the method described in Sect. 2.2 and to avoid cosmic variance uncertainties we performed our tests in overlapping areas between the Wide and the Deep, D1/W1 and D3/W3, over a total area of 2 deg^{2}.
These tests show that faint luminosity threshold samples are slightly incomplete at low redshift, probably as a consequence of catastrophic errors moving low redshift galaxies to higher redshift. We report that the difference in galaxy density between the Wide and the Deep estimates rarely exceeds the fieldtofield variance error estimate. However, brighter and brighter samples become more and more contaminated by spurious objects, and in the worst case (the brightest sample in the range 0.8 < z < 1.0), the number of objects in the Wide is more than three times higher than in the Deep, much larger than the 40% fieldtofield variance. A higher galaxy number density estimate will enhance the effect of a reduced clustering signal, which may result in underestimated halo masses for brighter samples up to a few sigmas. In the worst case, this represents a bias of ~0.5 in log M_{min}.
3.3. Angular correlation function
We measure the twopoint angular correlation function w using the Landy & Szalay (1993) estimator (3)where ⟨ DD ⟩ is the number of galaxy pairs, ⟨ RR ⟩ the number of random pairs, ⟨ DR ⟩ the number of galaxyrandom pairs, all in a bin around the angular separation θ. N_{d} and N_{r} are the number of galaxies and random objects, respectively. A random catalogue is generated for each sample. In order to find an acceptable compromise between small Poisson errors and computational requirements we scale N_{r} depending on the number of galaxies in our input catalogues, from 50 × N_{d} in the case of low density data catalogues to 2 × N_{d} for high density catalogues. We measure w using a fast twodimensional tree code. At large separations θ, instead of counting individual galaxy pairs, we correlate boxes by multiplying the number of objects in the two boxes. For that, we define a threshold angle α_{w} = θ/d_{b}, where d_{b} is the box size. Below this threshold, the number of objects inside the box is taken into account instead of individual objects. We found that α_{w} = 0.05 gives accurate results at low computational cost.
Fig. 6 Cross correlation analysis between redshift bins for the full sample. Top: autocorrelation in the redshift bin 0.4 < z < 0.6 (straight line) and crosscorrelation between the bins 0.4 < z < 0.6 and 0.8 < z < 1.0 (dotted line) and between 0.4 < z < 0.6 and 1.0 < z < 1.2 (dashed line). Bottom: quantitative estimates of the contamination (percentage of galaxies scattered) from a pairwise analysis between redshift bins 0.4 < z < 0.6 and 0.8 < z < 1.0 (left) and between 0.4 < z < 0.6 and 1.0 < z < 1.2 (right). The contours show the 68.3 (blue), 95.5 (green) and 99.7 (red) confidence regions. 
We correct our galaxy correlation measurements for the “integral constraint” (Groth & Peebles 1977), a correction factor which arises from the finite area of our survey Ω^{2}: (4)The correction factor w_{C} can be estimated as follows: (5)assuming a simple power law fitted on the data with slope γ and amplitude A_{w}: (6)As A_{w} varies with each sample, we first compute C using Monte Carlo integration over random pairs (7)and then correct w_{mes}(θ) using (8)The integral constraint for the four Wide fields is C ~ 0.5. Assuming γ = 1.8, this leads to a correction of ~10% at θ = 0.1 deg and up to a factor of two at θ = 1.0 deg. We note that assuming a powerlaw for w(θ) to estimate C could be a source of error; it is neglected in the current analysis as the correction is smaller than our Jackknife error estimates.
We combine our galaxy samples from the four CFHTLS fields into a single catalogue for the correlation function measurement. Computing w(θ) on the four fields independently, and using a pairweighted average would not lead to the same clustering estimate at large scales. Relative photometric offsets are expected to vary more from one field to another than within the same field, as the four fields are nonoverlapping. Since our parent catalogue is selected by apparent magnitude cut (i′ < 22.5), even a small difference in photometric offsets will lead to a spurious fieldtofield variation of the galaxy number density. In addition, different photometric offsets could bias the photometric redshift selection and also increase the fieldtofield density fluctuations. A 0.02 magnitude offset variation would result in a ~2% variation in galaxy number density, which may bias the twopoint correlation at very large scales. Because we expect these effects to occur at the scale of individual field (a few degrees), we adopt a conservative cut and limit our measurements to θ < 1.5deg. We do not measure w at separations below to avoid blended objects. Twopoint correlation function measurements and errors are provided in Appendix C.
3.4. w(θ) error estimates
We estimate statistical errors on the twopoint correlation function using the Jackknife internal estimator. We divide all samples into N = 68 subsamples of about 2deg^{2}. Removing one subsample at a time for each Jackknife realisation, we compute the covariance matrix as (9)where N is the total number of subsamples, the mean correlation function, and w^{l} the estimate from the lth Jackknife sample. Since w(θ) is computed from a combined catalogue including all four fields, our Jackknife estimate leads to a noisy but fair estimate of cosmic variance. We apply the correction factor given by Hartlap et al. (2007) to the inverse covariance matrix to compensate for the bias introduced by the noise. We show the correlation coefficient of the covariance matrix, for two red and blue samples, respectively, in Fig. 7.
Fig. 7 Covariance matrix correlation coefficients of w(θ) from Jackknife estimates, for galaxy samples with 0.4 < z < 0.6 and M_{g} −5log h < −19.8. Left panel: red sample. Right panel: blue sample. 
Fig. 8 Example of a measured w(θ) (for all galaxies in the redshift range 0.4 < z < 0.6 and for M_{g} − 5log h < −19.8), as well as the bestfitting model, as described in Sect. 4. Left: w(θ) measurement and model. Right: N(M), showing the central term N_{cent} and the satellite term N_{s}. 
4. Filling haloes with galaxies: the halo occupation model
4.1. The halo occupation distribution
In order to relate galaxies to the dark matter haloes which host them we have implemented an analytic model of galaxy clustering, the halo model (for a review see Cooray & Sheth 2002), which contains at its core a prescription for how galaxies populate haloes, namely the halo occupation distribution. Our model follows closely the approaches used in recent works; further details and references can be found in Appendix A. The key assumption underlying our halo occupation distribution function is that the number of galaxies N in a given dark matter halo depends only on the halo mass M; it does not depend on environment of formation history of the haloes. Furthermore, following Zheng et al. (2005), we express N(M) as the sum of two terms, corresponding to the contribution from the central galaxy N_{c} and the satellite galaxies N_{s}: (10)where The smooth transition for central galaxies expresses the uncertainties in the galaxy formation process (Zheng et al. 2007). The factor N_{c}(M) for the satellite number in N(M) accounts for the fact that a halo cannot be populated by satellite galaxies without the presence of a central galaxy.
For a given cosmology and dark matter halo profile, our model has five adjustable parameters. M_{min} is the mass scale for which 50% of haloes host a galaxy. To reflect the scatter in the luminosityhalo mass relation, a smooth transition of width σ_{log M} is used. Zheng et al. (2007) show that if this scatter is small, the above expression takes the identical form as the distribution of central galaxies given a halo mass M, integrated over the entire luminosity range above the luminosity threshold. Thus M_{min} also represents the halo mass scale for central galaxies, whose mean luminosity ⟨ L_{c} ⟩ is equal to the luminosity threshold L_{min}.
This simple relation between M_{min} and L_{c} is based on the hypothesis that stellar mass (or luminosity) has a power law dependence on halo mass, which may not be exact over the entire mass range. Leauthaud et al. (2012a) recently proposed a model in which this relation assumes a more realistic form. The authors showed that in their model M_{min} and σ_{log M} take different values than those computed with other models. However, in the mass range over which we compare our results with theirs in Sect. 6 (M_{h} ~ 10^{12} h^{1} M_{⊙}), M_{min} values do not differ by more than 10% (see Fig. 3 in Leauthaud et al. 2012a).
The number of satellite galaxies as function of halo mass follows a power law with slope α and amplitude M_{1}. M_{1} then represents the characteristic scale for haloes hosting one satellite galaxy. At lower masses, the dependence becomes steeper and the transition mass scale occurs at M ~ M_{0}.
We show in Fig. 8 an example of measured w(θ) and its bestfitting model, together with the bestfitting HOD function N(M). The total galaxy correlation function is the sum of two terms. At distances much smaller than the virial radius, the onehalo term contains contributions from galaxy pairs within a single halo, whereas at large distances the twohalo term contains contributions from pairs in separate haloes.
4.2. Deduced parameters
From the HOD model we obtain deduced parameters describing galaxy properties. The mean galaxy bias b_{g} at redshift z is the massintegral over the halo bias b_{h} (see Sect. A) weighted by the number of galaxies, (13)The darkmatter mass function n is given in Eq. (A.6), and the halo bias b_{h} in Eq. (A.14). The total number of galaxies is (14)Similarly to the galaxy bias, the mean halo mass for a galaxy population is (15)The fraction of central galaxies per halo is (16)Consequently, the fraction of satellite galaxies is (17)
4.3. Population Monte Carlo sampling
We use the “Population Monte Carlo” (PMC) technique to sample likelihood space (Wraith et al. 2009; Kilbinger et al. 2010, 2011) and have implemented it in a publiclyavailable code, CosmoPMC^{2}. Contrary to the widelyused Monte Carlo Markov Chain method (MCMC), PMC is an adaptive importancesampling technique (Cappé et al. 2004, 2008), which has two principal advantages: first, the parallel sampling algorithm combined with our fast HOD code means that each run can be quickly computed. Secondly, PMC does not have issues with chain convergence. Instead, the perplexity (defined below) is a diagnostic for the reliability of a given sampling run.
Points are sampled from a simple, socalled importance sampling function q. Each sample point θ_{n} is attributed a weight which is the ratio of the posterior π to the importance function, (18)The initial importance sampling function q is a mixture of seven multivariate Gaussians. A PMC run consists of a number of iterations, in the course of which the importance function is adapted to better match the posterior distribution. We run PMC for typically 10 iterations using 10 000 sample points per iteration.
As a stopping criterion for the iterations we take the perplexity p which is defined as (19)The perplexity p ∈ [0;1] is a measure of the distance between the posterior and the importance function, and it approaches unity if the two distributions are identical. The perplexity is a measure of the adequacy and efficiency of the sampling. Values of p above about 0.6 indicate a good agreement between the posterior and the importance function.
The evaluation of the likelihood function, which is the most timeconsuming process for most sampling tasks, is easily parallelisable in importance sampling with little overhead. Furthermore, we have optimized several aspects of the numerical computations of the HOD model. For example, we employ the FFTLog (Hamilton 2000) algorithm to perform Fourier transforms and make use of tabulated values to improve efficiency. The computation of the angular correlation function w on a range of scales is performed in under a second on a standard desktop. This module is part of the latest public version (v1.1) of CosmoPMC. Typically, on an eightcore desktop (16 threads) a single PMC run takes around 90 min for a total of 100 000 sample points. Our PMC technique allows us to efficiently sample the parameter space for the large number of galaxy samples used in this work.
4.4. Likelihood function
For each galaxy sample we simultaneously fit both the projected angular correlation function w and the number density of galaxies n_{gal}, by summing both contributions in loglikelihood: (20)where is given by Eq. (14), at the mean redshift of the sample. The data covariance matrix C is approximated by Eq. (9). The error on the galaxy number density σ_{ngal} contains Poisson noise and cosmic variance. The latter is estimated from the fieldtofield variance between the four Wide patches.
The likelihood function is (21)The product of prior P and likelihood defines the posterior distribution, π = Pℒ. We sample a fivedimensional space where { log _{10}M_{min},log _{10}M_{1},log _{10}M_{0},α,σ_{log M} } is the parameter vector. We use flat priors for all parameters; the ranges are log _{10}M_{min} ∈ [11;15] ,log _{10}M_{1} ∈ [12;17] ,log _{10}M_{0} ∈ [8;15] ,α ∈ [0.6;2] and σ_{log M} ∈ [0.1;0.6] .
Fig. 9 w(θ) measurements (left) with their best fit models for all galaxies, and the resulting galaxy occupation function (right) with respect to halo mass. We illustrate the redshift evolution (top) in the luminosity threshold M_{g} − 5log h < − 19.8, and the luminosity dependence (bottom) in the redshift range 0.4 < z < 0.6. 
Fig. 10 w(θ) measurements (left) with their best fit models for red galaxies, and the resulting galaxy occupation function (right) with respect to halo mass. We illustrate the redshift evolution (top) in the luminosity threshold M_{g} − 5log h < − 19.8, and the luminosity dependence (bottom) in the redshift range 0.4 < z < 0.6. 
5. Results
5.1. Clustering measurements
Our measurements are presented in Figs. 9–11 for the full, red and blue samples, respectively. For clarity, we show three samples selected with the same luminosity threshold (M_{g} − 5log h < −19.8) and three samples covering the same redshift range (0.4 < z < 0.6). At the same luminosity and redshift, red samples are more clustered on small scales than the full sample, and blue galaxies are the least clustered. At a constant luminosity threshold, w(θ) increases with redshift for the full, red and blue samples.
Fig. 11 w(θ) measurements for blue galaxies. We illustrate the redshift evolution (left) at the luminosity threshold M_{g} − 5log h < − 19.8, and the luminosity dependence (right) in the redshift range 0.4 < z < 0.6. 
Fig. 12 Top panels: halo mass estimates M_{min} and M_{1} for all (left) and red (right) galaxy samples as function of luminosity threshold (M_{g}). Bottom panels: M_{1} versus M_{min} in different redshift bins. The dotted lines represent a linear relation and the dashed lines are a powerlaw fit to M_{1} versus M_{min} in each redshift bin. 
On the left panels of Figs. 9 and 10, we plot our best fit w(θ) models. Model parameters and deduced quantities are shown in Tables B.1 (full) and B.2 (red). An example of sample distributions of bestfitting HOD parameters is given in Fig. C.1. Since the halo model aims to reproduce the clustering in threshold samples as a sum of contributions from central and satellite galaxies, both types of galaxy must be present. We presume that in general, central galaxies are brighter than their satellites, so that the latter assumption is valid for full samples selected by luminosity threshold. Also, we assume that red satellites are systematically associated with a (brighter) red central galaxy. On the contrary blue satellites may belong to a red central galaxy. For this reason, we defer fitting our model to our pure blue samples. We plan to address the blue population in a future work, perhaps considering separate halo occupation functions for both red and blue populations (Simon et al. 2009; Skibba & Sheth 2009).
In most cases, we find a good agreement between the data and our model, with χ^{2}/d.o.f. from 1.77 to 10.21 for full samples and from 1.20 to 6.21 for red samples. We report a few large χ^{2} values in some full samples, mostly at high redshift and for bright luminosity thresholds, where the galaxy number density is small. We suspect systematic photometric redshift errors to be responsible for an underestimated clustering signal and systematic errors in the galaxy number density estimates (see Sect. 3.2). Perplexity (see Sect. 4.3) is larger than 0.6 for most samples, indicating that our model accurately describes the observations.
5.2. Halo model fitting
In the right panels of Figs. 9 and 10 we show our best fit N(M), the halo occupation function. First, we note that at a given redshift interval, M_{min}, the average mass for which 50% of haloes contain one galaxy increases with luminosity. Secondly, we note that the evolution in N(M) with redshift at a fixed luminosity is less pronounced than the evolution with luminosity at a fixed redshift, whereas the amplitude of w(θ) significantly decreases with redshift.
Halo masses M_{1} and M_{min} are tightly constrained, whereas M_{0}, σ_{log M} and α are poorly constrained. M_{min} and M_{1} measurements are displayed in the top panels of Fig. 12 with respect to luminosity threshold (M_{g}). For every sample, halo masses increase with luminosity. M_{min} covers the range 10^{11} to 10^{14} h^{1} M_{⊙}. M_{1} is significantly larger than M_{min}, ranging from 10^{12.5} to 10^{16} h^{1} M_{⊙}. These trends are consistent with those found in both the local (Zehavi et al. 2011, 2005) and distant Universe studies (Zheng et al. 2007; Abbas et al. 2010). At high luminosity thresholds, where the number of objects is small, error bars on our bestfitting parameters are correspondingly larger. This trend is especially pronounced for α, where the small number of satellites in high mass haloes makes it difficult to constrain this parameter. However, we note α increases with luminosity threshold, suggesting that for massive haloes, bright galaxies will be more numerous than faint ones. σ_{log M} remains relatively stable for red galaxies (~0.3), as also seen in other data sets and numerical simulations. For the full sample, values are slightly larger (up to ~ 0.5 in some cases). Our fitted σ_{log M} values are slightly larger than previous works. We also made fits where values σ_{log M} > 0.6 were allowed, and in some cases large values ( > 1.0) were found. We therefore decided to impose the restriction that σ_{log M} < 0.6. Sample incompleteness (due to errors in the photometric redshifts), leading to missing central galaxies, could explain such high σ_{log M} values.
At constant luminosity threshold, in both red and full samples, M_{1} and M_{min} decrease with increasing redshift. This effect is partially due to a selection effect caused by the dimming and reddening of stellar populations with time; for a constant luminosity cut one selects less and less massive galaxies at higher redshifts. We will return to this measurement in Sect. 6 where we attempt to separate this passive stellar evolution effect from the intrinsic stellartohalo mass evolution.
Fig. 13 M_{min} as a function of galaxy number density for all (left) and red (right) galaxies. The dashed line in each panel represents a power law fit to the data, and the dotted line represents a power law fit to the results given in Brown et al. (2008) (red galaxies) and in the SDSS by Zehavi et al. (2011) (all galaxies). We converted M_{200} halo masses to virial masses M_{vir} when necessary. 
At faint luminosities, M_{min} values are higher for red galaxies than for full samples. As we explained previously, one assumption of our model is that our red threshold samples have luminous red central galaxies. On the other hand, it is also possible, especially for fainter samples, that blue satellite galaxies inhabit haloes with red central galaxies. In such a case for the faint red galaxy sample, a red central galaxy will be seen “alone” in a more massive halo. This would also increase the halo mass estimate found for red central galaxies as compared to those in the full sample.
The bottom panels of Fig. 12 show M_{1} with respect to M_{min} for full and red samples. The parameters are strongly correlated and we fitted a power law (dashed line) in each redshift bin. For the full sample, we compare our results with a linear relation, M_{1} = 17 × M_{min} (dotted line), based on the coefficient observed in the SDSS (z ~ 0.1) by Zehavi et al. (2011). Our results agree with them in the lowest redshift bin, z ~ 0.3, but depart from this value at higher redshift.
We find that although M_{1}/M_{min} relation does not vary significantly with luminosity threshold it does show some indications of evolution with redshift. We note that M_{1}/M_{min} decreases with redshift for M_{min} ≲ 10^{13} h^{1} M_{⊙} and increases with redshift at higher mass. For red galaxies, M_{1} is no longer proportional to M_{min} but follows a power law. For M_{min} ~ 10^{12} h^{1} M_{⊙}, the ratio M_{1}/M_{min} is ~12 (as compared to ~10 for luminous red galaxies found by Zheng et al. 2009) and shows a much smaller redshift evolution over the entire mass range.
Figure 13 shows the relationship between M_{min} and the galaxy number density n_{gal} for red and full samples. As one would expect, progressively rarer objects are found in more massive haloes. We do not detect a significant evolution of the slope of our M_{min} versus n_{gal} relationship with redshift. For red galaxies, we overplot the measurements from Brown et al. (2008) with which we are in excellent agreement, even though it covers a much smaller area. For the full galaxy sample we compare with SDSS measurements (Zehavi et al. 2011). The agreement is good at high surface densities; however we note a departure for low density samples. This difference suggests that at constant density, galaxies tend to reside in more massive haloes at low redshift (z ~ 0) than at higher redshift (z > 0.6). Part of this difference could be caused by systematic photometric redshift errors. As we discussed in Sect. 3.2, bright samples may be contaminated by fainter objects, and the combined effect of higher density and lower clustering signal will result in reducing M_{min}. In the most pessimistic case outlined in Sect. 3.2, where the density estimate would be overestimated by a factor 3, log M_{min} would be biased low by a few sigmas (up to ~0.5).
The relationship between M_{min} and n_{gal} can be well approximated by a power law and we fitted galaxy samples at all redshifts simultaneously. We found: (22)for all galaxies and (23)for red galaxies.
6. Discussion
The key issue we would like to address is to understand how galaxy formation and evolution depend on the properties of the underlying dark matter haloes. Ideally, to relate halo masses to their stellar content, one would select our samples by stellar mass. However, stellar mass estimates computed from our fiveband optical data alone would suffer from large uncertainties. Keeping this limitation in mind, we adopt an intermediate approach where we determine an empirical relation to convert the observed luminosity threshold of each sample to an approximately massselected sample using a reference sample with accurate stellar masses. This method allows to transform luminosities into masses and provides an estimate of the likely errors in such a transformation.
6.1. Transforming to stellar mass threshold samples
As a consequence of passive stellar evolution, the characteristic restframe M_{g}band luminosity of a given galaxy type evolves significantly over the redshift range 0.2 < z < 1.2 (Ilbert et al. 2005; Zucca et al. 2006; Faber et al. 2007).
In order to account for these effects, we have established a simple approximation to relate stellar mass and luminosity based on COSMOS 30band photometry (Ilbert et al. 2010). We used the COSMOS Bband luminosities in order to be as close as possible to our M_{g}band selected samples. Red galaxies are selected as M(NUV) − M_{R} > 3.5 and blue galaxies as M(NUV) − M_{R} ≤ 3.5. We then fit the stellar masstoluminosity ratio with powerlaw functions as function of redshift between z = 0 and z = 1.5 for five stellar mass bins from 10^{9.0} to 10^{11.5}. These results are displayed in Fig. 14.
Due to the ageing of stellar populations, the average stellar masstolight ratio increases with time. For red galaxies, the slope depends weakly on stellar mass at intermediate masses (10^{9.5}−10^{11.5} h^{1} M_{⊙}) and ranges from −0.56 to −0.44. The slope is steepest in the lowest mass bin, but the number of objects in this bin is small. Based on these results we adopt the following relationship for the red galaxy population: (24)where the 0.1 error accounts for the scatter between mass bins. This evolution is consistent with the luminosity function evolution measured for red galaxies in the DEEP2 and COMBO17 surveys (Faber et al. 2007; Willmer et al. 2006), where the characteristic absolute magnitude decreases by about 1.2–1.3 mag per unit redshift (M_{B} ∝ − 2.5log L_{B}).
Fig. 14 Stellar masstoluminosity relations for red (left) and blue (right) galaxies, as a function of stellar mass and redshift as measured in the COSMOS survey (Ilbert et al. 2009). Each panel shows a separate slice in stellar mass. Samples are fitted with a power law (dashed line). 
For the blue sample, the evolution in redshift is more pronounced than for the red sample, as also seen in DEEP2 (Willmer et al. 2006) and VVDS (Zucca et al. 2006). Unfortunately, in this case, the slope depends on stellar mass, ranging from − 0.69 to − 1.10 in the mass range [10^{9.5},10^{11.5}] , and has a larger scatter. The full sample is dominated by blue galaxies at faint luminosities but by red galaxies at bright luminosity (see Fig. 3). Therefore a simple relation between luminosity and stellar mass valid for all luminosities cannot be obtained. For the full sample we choose to apply the same correction as for the red sample. As blue galaxies evolve more rapidly than red ones such a correction underestimates the expected evolution of the full samples at faint luminosity and therefore cannot be used to construct stellarmass selected samples. In the following sections, we discuss how our results depend on this correction.
In the following sections we use these expressions to calculate “corrected” luminosities, denoted as L′. Note that these corrections are equivalent to shifting observed galaxy luminosities to redshift z = 0. Our results are displayed as a function of relative luminosity threshold L′/L_{∗}, where L_{∗} corresponds to and −19.75 for the full and red samples, respectively, measured in the local Universe (Faber et al. 2007, and references therein). It is also worth mentioning that if galaxies would experience a pure passive luminosity evolution with cosmic time, this would result in constant HOD parameters for a given L′/L_{∗} as function of redshift. However, deduced parameters (galaxy bias, mean halo mass and satellite fraction) are expected to vary due to the evolving halo mass function and halo bias (see Eqs. (13)–(17)).
6.2. The stellartohalo mass relationship
Fig. 15 Left panels: halo mass estimates M_{min} and M_{1} for all (top) and red (bottom) galaxy samples, as function of luminosity threshold, corrected for passive redshift evolution to approximate stellar mass selected samples. The dashed lines correspond to Eq. (25). Right panels: lighttohalo mass ratios (Eq. (26) with identical parameters as those fitted with Eq. (25)) as function of halo mass. 
In our model, M_{min} represents the hosting halo mass for central galaxies whose median luminosity equals the luminosity threshold of the sample (see Sect. 4.1). The ratio L′/M_{min} is therefore the ratio of luminosity to dark matter mass for central galaxies. Zehavi et al. (2011) showed that the relation between central galaxy luminosity (corresponding to our L′) and the halo mass could be approximated by: (25)where A, α_{M} and M_{t} are free parameters and M_{h} represents M_{min} or M_{1}.
This expression encapsulates the idea that there exists a “transition halo mass” (~M_{t}) where the ratio of central galaxy luminosity to dark matter halo mass reaches a maximum. This transition mass represents the halo mass for which baryons have been most efficiently converted into stars until the time of observation, resulting from past and ongoing star formation over cosmic history.
We fit Eq. (25) for full and red samples for each redshift slice. For our two highest redshift bins where the number of points are less than the number of free parameters, we perform a simultaneous fit over both redshift bins. Results are plotted in Fig. 15 and bestfitting parameters given in Table C.6.
In our lowest redshift bin (0.2 < z < 0.4), our full sample has bestfitting parameters slightly different from in the SDSS (Zheng et al. 2007; Zehavi et al. 2011), A = 0.20 ± 0.01 (compared to 0.32) and M_{t} = 2.24 ± 0.1 × 10^{11} h^{1} M_{⊙} (compared to 3.08 × 10^{11}). We note that the difference is likely due to the different selection (M_{r} in the SDSS). At high luminosities, the difference with their relation increases. For a given luminosity, we find lower halo masses, L_{c} ∝ M^{0.5 ± 0.01}, compared to L_{c} ∝ M^{0.26 − 0.28}, which follows the trends observed in Fig. 13, and might be due to contaminated samples at low density. In the red sample, we measure L_{c} ∝ M^{0.31 ± 0.16}, whereas Zheng et al. (2009) found L_{c} ∝ M^{0.5} for luminous red galaxies in the SDSS. Both M_{min} and M_{1} increase with redshift for the full sample. The dependence with redshift is weak for red galaxies, although we detect a slight increase in the very highest redshift bins.
For the full sample we assumed a correction based on the evolution of the stellar masstoluminosity ratio of red galaxies (Eq. (24)). As the characteristic luminosity of blue galaxies evolves more rapidly than red ones and as the faint sample is dominated by blue galaxies a legitimate concern is the dependence of our results on our corrections to stellar mass. Several corrections assuming the blue galaxy luminosity evolution (using slopes given in Fig. 14) were tested. If a larger correction is applied, M_{1} and M_{min} show a similar but more pronounced trend. For this reason we conclude that the increase of halo mass parameters with redshift for the full sample is a robust result (although the exact amount of evolution observed depends on the correction we apply).
6.3. Redshift evolution of the SHMR and comparison with literature values
Fig. 16 Position of the peak corresponding to for the full (filled circles) and red (filled diamonds) samples as function of redshift. Error bars are derived from the nearest data point to the peak. We compare our results with measurements in SDSS (Zehavi et al. 2011), in COSMOS (Leauthaud et al. 2012b) and in the NEWFIRM Medium Band Survey (NMBS, Wake et al. 2011). We convert M_{200} halo masses (corresponding to a sphere with overdensity 200 times the background density) into virial masses M_{vir} (our definition) using the method described in Appendix C of Hu & Kravtsov (2003). We also display numerical predictions based on abundancematching method from Behroozi et al. (2010) and Moster et al. (2010), where the shaded region shows the 1 − σ uncertainty. 
To study the redshift evolution of the stellartohalo mass ratio, we consider the relationship between corrected galaxy luminosity and the host halo mass from Eq. (25): (26)where the maximum is located at: (27)The right panels of Fig. 15 show the relation for full and red samples as function of halo mass fitted by Eq. (25). Best fitting parameters are listed in Table C.6. For less massive haloes, the relative stellar mass content increases with halo mass and reaches a maximum around M_{min} ~ 10^{12} h^{1} M_{⊙}. In more massive haloes, the stellartohalo mass ratio sharply decreases. Our measurements are consistent with recent observations measuring the stellartohalo mass ratio as function of halo mass (Mandelbaum et al. 2006; Zheng et al. 2007; Foucaud et al. 2010; Zehavi et al. 2011; More et al. 2011; Wake et al. 2011; Leauthaud et al. 2012b), and “abundancematching” techniques based on Nbody simulations (Conroy & Wechsler 2009; Moster et al. 2010; Behroozi et al. 2010). Similar trends are found in semianalytic simulations (see, for e.g., Somerville et al. 2008), in which physical processes (such as supernovae and AGN feedback) preventing gas from cooling and quenching star formation must be taken into account to match observations.
Our measured peak values (Eq. (27), with error bars from the nearest data points) are plotted as a function of redshift in Fig. 16. Comparing to literature measurements at lower redshifts, the peak stellartohalo mass relation is , similar to the local Universe value (≈4.2 × 10^{11} h^{1} M_{⊙}, Zehavi et al. 2011). For our full sample, the peak location gradually increases with redshift.
At higher redshifts, we compare our measurements with those of Moster et al. (2010) and Behroozi et al. (2010). These authors derived the stellartohalo mass relation using an abundancematching technique, linking dark matter halo mass functions from Nbody simulations with observed galaxy stellar mass functions. Moster et al. (2010) estimate the redshift evolution of the stellartohalo mass ratio (their Table 7) and from this we derived analytically the position of the peak, as well as associated errors. Values from Behroozi et al. (2010) at redshift 0.1, 0.5 and 1.0 were supplied to us (Behroozi, priv. comm.). Both Moster et al. and Behroozi et al.’s measurements are based primarily on galaxy stellar mass measurements covering relatively small areas (a few hundred square arcminutes). Our full sample is in excellent agreement with Behroozi et al. but differs significantly from Moster et al., although in both cases the two data sets show the same general trend: shifts gradually towards higher halo masses at higher redshift. The Moster et al. points do not seem to agree very well with most observations and moreover overestimate the position of the peak at low redshift.
Measurements in the COSMOS field are also shown (Leauthaud et al. 2012b); their points are in general larger than most measurements, with the exception of Moster et al. It has already been noted that there are several rich structures in the COSMOS field (McCracken et al. 2007; Meneux et al. 2009; de la Torre et al. 2010) which lead to higher correlation functions, especially at z ~ 0.8; cosmic variance could perhaps be the origin of this discrepancy.
For red galaxies, the peak is at larger halo masses, . This is not surprising as M_{min} is larger for red galaxies at faint luminosities. The red galaxy peak position undergoes only little evolution to z ~ 1. Our intermediateredshift red galaxy points are consistent with higher redshift measurements (Wake et al. 2011); these authors used a nearinfrared survey to select galaxies by stellar mass, and we would expect their selection to be dominated by red galaxies.
The results presented here represent the first time a single data set covering a statistically significant area has been used to derive the evolution of the peak position as a function of redshift and galaxy type in a selfconsistent manner. The peak position, , can be interpreted as representing the halo mass at which the stellar mass content is most efficiently accumulated in haloes, either by star formation or merger processes. The movement of this peak towards higher halo masses at higher redshifts is consistent with a picture in which the sites of efficient stellar mass growth migrate from lowmass to highmass haloes at higher redshift. This “antihierarchical” evolution, frequently referred to as “halo downsizing”, resembles closely the scenario first sketched by Cowie et al. (1996), in which the maximum luminosity of galaxies undergoing starformation declines steadily to the current day. Conversely, the fact that we observe constant transition masses for the red population suggests that the stellartohalo mass ratio does not significantly evolve in haloes with mass ~10^{12} h^{1} M_{⊙} and supports observations that intermediatemass red galaxies have experienced no significant stellar mass growth by star formation since z ~ 1.2.
We tried several corrections to convert to stellar mass thresholds based on the slope measured for blue galaxy evolution (see Fig. 14), and we found that the peak luminosity depends only weakly on these corrections. Only the highestredshift point changes by around one sigma if we apply the extreme “blue galaxy” correction to our full sample. For this reason we conclude that the increase of with redshift for the full sample is a robust result.
6.4. Redshift evolution of halo satellite fraction
The fraction of galaxies which are satellites in a given dark matter halo is a sensitive probe of the past evolutionary history of the halo which may be modified by processes such as major mergers (which can decrease the satellite fraction with time) or environmental effects which operate primarily to reduce the number of satellites, such as, but not limited to, the restriction of gas supply or “strangulation” (Larson et al. 1980) or rampressure stripping (Gunn & Gott 1972). Satellite galaxies may also merge with their central galaxy. Observational evidence based on pair fraction measurements and theoretical studies using numerical simulations now shows that mergers seem to play a significant (although not dominant) role in the evolution of massive galaxies since z ~ 1 (White et al. 2007; Wang & Jing 2010; van Dokkum et al. 2010; Robaina et al. 2010; Zehavi et al. 2012). Our objective in this section is to see if measurements of the satellite fraction can provide new insights into the evolutionary history of haloes.
The left panels of Fig. 17 show the satellite fraction f_{s}, computed using Eq. (17). As before, each sample is converted to an approximately masslimited sample by “correcting” for passive luminosity evolution (see Sect. 6.1).
Fig. 17 Satellite fractions as function of luminosity and redshift, for all and red samples (top and bottom panels). The left panels show the measured satellite fraction as function of corrected luminosity. Dashed lines in the middle panels show satellite fractions extrapolated in redshift, based on our HOD parameters measured in the redshift range 0.4 < z < 0.6. We compare our results with the literature for all galaxies at z = 0.1 and z = 1.0 (top, Zheng et al. 2007) and with galaxygalaxy lensing estimates at z = 0.1 (bottom, Mandelbaum et al. 2006) for red galaxies. In the right panels we provide values extrapolated in redshift for several L′/L_{∗} ratios. 
At all redshifts, as the luminosity threshold increases, corresponding to higher mass thresholds, the satellite fractions become progressively smaller. The similarity between the full and red samples is simply a consequence of the dominance of red galaxies at bright luminosity thresholds.
At fainter luminosity thresholds, the satellite fraction rises and eventually reaches a plateau. This effect can be understood by considering both the halo mass function and the halo occupation function: the increase of the latter with decreasing luminosity threshold (see Fig. 8) is rapidly compensated by a sharp decrease in massive halo number density, so that faint galaxies are preferentially central galaxies in smaller haloes rather than being satellites in more massive haloes.
We compare our results with satellite fractions found in SDSS (z ~ 0.1) and DEEP2 (z ~ 1). To explore a complete range in redshift and luminosity, we extrapolate our bestfitting HOD parameters measured in the redshift range 0.4 < z < 0.6; that is, at a given luminosity, we keep the HOD parameters fixed and compute Eq. (17) at different redshifts. The choice for this intermediate redshift range is motivated by the fact that it contains the most robust photometric redshift estimates and our measured points span a large range in luminosity. Each HOD parameter set is constructed from M_{min} and M_{1} interpolated in luminosity using Eq. (25) with bestfitting parameters given in Table C.6. For red galaxies, M_{0} is set to M_{min}/100, σ_{log M} is fixed to 0.3 and α = 1.1 for L < L_{∗} or α = 1.4 otherwise. For all galaxies, M_{0} is set to M_{min}/100, σ_{log M} = 0.4 and α = 1.1 for L < L_{∗}, and σ_{log M} = 0.3 and α = 1.3 for brighter luminosity thresholds.
In the middle panels of Fig. 17 we compare our extrapolated satellite fractions for all galaxies (in dashed lines) with results derived from SDSS and DEEP2, using a HOD analysis similar to ours (Zheng et al. 2007). Our satellite fractions are consistent with these measurements at both redshifts, although we slightly overpredict satellite fractions at z ~ 0.1 for samples brighter than ~0.5 L_{∗}. An overall increase of a factor of two is measured between z ~ 1 and z ~ 0 for samples with L > L_{∗}; we note that this value is lower than the factor of three from z ~ 0.5 to z ~ 0 measured in the VVDS (Abbas et al. 2010).
We fit the HOD parameters at each redshift independently. Therefore, our model provides us with “snapshots” of the occupation function at various epochs. If any physical process had played a significant role in reducing the satellite number density between z ~ 1 and z ~ 0, our fixedHOD model calibrated at z ~ 0.5 would overpredict the satellite fraction observed in the local Universe and underpredict the satellite fraction at redshift one. The excellent match found for bright luminosity threshold samples shows that our satellite fraction measurements follow closely a redshiftindependent HOD. Since a redshift evolution of the merging rate would change the occupation function, we therefore conclude that our measurements imply a constant merging rate since z ~ 1, which confirms the constant (albeit minor) role of merging for intermediatemass galaxies (e.g., Bundy et al. 2009). As outlined by White et al. (2007), if halo merging was not followed by galaxy merging, we would observe an increase of halo masses at constant galaxy density. Instead, no redshift evolution of the M_{min}n_{gal} relationship (see Fig. 13) supports a scenario with constant galaxy merging with time.
We compare our red satellite fraction with the galaxygalaxy lensing analysis from SDSS (Mandelbaum et al. 2006). These samples were constructed from luminosity slices, rather than the threshold samples we use here. We then computed the number of satellites per luminosity bin as the difference of satellite fractions between two threshold samples: (28)As we found for the full galaxy sample, the good agreement supports a nonevolving HOD from z = 0.5 to z = 0.
Finally, in the right panels of Fig. 17, we provide extrapolated satellite fractions as function of redshift for several luminosity thresholds. It is interesting to examine our results in the context of recent results studying the role of satellite galaxies in galaxy formation and evolution (Peng et al. 2010, 2011). These works show measurements from numerical simulations and SDSS group catalogues which indicate that the fraction of satellites at a fixed overdensity is independent of halo mass and redshift. In this context, the growth of the satellite fraction we observe is simply a consequence of average growth in overdensity in rich structures between z ~ 1 to z ~ 0, reflecting the growing importance of environment at low redshift.
6.5. Galaxy bias
We compute the galaxy bias (Eq. (13)) from our bestfitting HOD parameters. We fit the bias measurements as function of luminosity and redshift, adopting a linear relationship (Norberg et al. 2001): (29)where L′ represents the corrected luminosity (computed from Eq. (24) for the full and red samples). We fit the two highest redshift bins together (0.8 < z < 1.0, 1.0 < z < 1.2), as only one point is available in the highest one. Results are plotted in Fig. 18 and bestfitting parameters are given in Table C.7.
Fig. 18 Galaxy bias as function of luminosity and redshift. Left panels: bias for full and red samples (top and bottom panels) as function of corrected luminosity threshold L′. The dashed line represents Eq. (29). Right panels: computed bias for several luminosity threshold samples as function of redshift. Error bars were derived from the parametric form shown on the left (dashed line), whereas the dotted lines show the bias evolution computed from our model, assuming a redshiftindependent HOD parameter set. 
This parameterised bias relation reproduces our results best at z < 0.6. At a fixed redshift, for L′/L_{∗} < 1, the bias is a very weak function of luminosity; for L′/L_{∗} > 1, bias depends strongly on luminosity. For the faint, lowredshift fullsample galaxies, dominated by more weaklyclustered bluer galaxies, b ~ 1; the faint red samples have higher biases, of around b ~ 1.5. These results are a consequence of the higher clustering amplitudes observed for red samples compared to the full galaxy sample for L′/L_{∗} < 1, and are consistent with similar measurements made by smaller surveys at intermediate redshifts (see for example Meneux et al. 2009).
The right panels of Fig. 18 illustrate the redshift evolution of the bias at several luminosity thresholds. We compare our results (solid lines) derived from Eq. (29) with a galaxy bias evolution computed using a constant set of HOD parameters evaluated at z ~ 0.5 (see Sect. 6.4). The slight difference observed between the measured HOD and the HOD fixed model at z = 0.5 comes from the fact that we set σ_{log M} to 0.3 and α to 1.0. For the full galaxy population, a model with fixed HOD parameters provides an excellent fit to the observations, with more luminous objects undergoing a much more rapid evolution in bias. For the red galaxy population, there is some evidence that the bias evolution is more rapid than the model with fixed HOD parameters.
6.6. Mean halo mass
From Eq. (15) we compute the average dark matter halo mass for each of our samples. Note that this quantity represents the halo mass averaged over samples containing galaxies above the luminosity threshold and differs from M_{min} which represents the halo mass for central galaxies whose luminosity corresponds to the luminosity threshold. As before, after correcting luminosities to a corresponding stellar mass threshold, we fit the mean halo mass as function of luminosity and redshift, adopting a linear relation: (30)where L′ is the corrected luminosity. Once again, we fit the two highest redshift bins simultaneously. Results are displayed in Fig. 19 and bestfitting parameters are given in Table C.7.
Fig. 19 Mean halo mass as function of luminosity and redshift. Left panels: mean halo mass for full and red samples (top and bottom panels) as function of corrected luminosity threshold L′. The dashed line represents Eq. (30). Right panels: computed mean halo mass for several luminosity threshold samples as function of redshift. Error bars were derived from the parametric form shown on the left (dashed line), whereas the dotted lines show the mean halo mass evolution computed from our model, assuming a redshiftindependent HOD parameter set. 
For fainter samples with L′ < 2 L_{∗} the mean halo mass gradually increases with luminosity and changes more rapidly for brighter samples. Faint red galaxy samples on average have higher halo masses than the full sample which is dominated by bluer galaxies at these luminosities. This is consistent with a scenario in which high numbers of red satellite galaxies in faint samples reside in more massive haloes.
In the right panels of Fig. 19 we show the redshift evolution of the mean halo mass for several luminosity thresholds. As before, we compare mean halo masses derived from our measured HOD parameters (solid lines) with a model of fixed halo model parameters set measured at redshift 0.5. At all luminosities, the mean halo mass decreases with increasing redshift. This is an indication of the fact that haloes merge through cosmic history and the effective mass of haloes grows with time. In bright samples, the mean halo mass evolves more slowly than in fainter ones. At low luminosity thresholds and high redshifts, both the full sample and red galaxy sample deviate slightly from our model with fixed HOD parameters.
7. Summary and conclusions
We have made one of the most precise measurements to date of the angular correlation function w and its dependence on luminosity and restframe colour from z ~ 0 to z ~ 1. Our measurements were constructed from a series of i′ < 22.5 volumelimited samples, containing ~3 × 10^{6} galaxies in four independent fields. These samples were created using accurate fiveband photometric redshifts in the CFHTLS Wide survey calibrated with ~10^{4} spectroscopic redshifts. A crosscorrelation analysis showed consistency with a relatively small contamination fraction arising from scatter between photometric redshift bins.
We interpreted these measurements in the framework of a model in which the number of galaxies inside each dark matter halo is parametrised by a simple analytic function (the halo occupation distribution) with five adjustable parameters. In this model, central and satellite galaxies are treated separately. For each galaxy sample, we performed a full likelihood analysis and explored model parameter space using an efficient Population Monte Carlo analysis. The very large survey volume probed by the CFHTLS Wide allowed us to place robust constraints on the redshift evolution of the halo model parameters, and make accurate estimates for halo properties of many different galaxy samples over more than two orders of magnitude in luminosity and three orders of magnitude in halo mass, from 10^{12} to 10^{15}(h^{1} M_{⊙}).
Using 30band photometric data from the COSMOS survey, we derived an empirical relation between stellar mass, luminosity and redshift and used this to convert our luminosityselected samples to stellarmass limited samples. All of our conclusions (following below) are independent of this correction.
These are our principal results:
 1.
For a given luminosity threshold, galaxies with redderrestframe colours are more clustered than bluer ones. Clusteringstrength increases with increasing luminosity, reflecting thatbright galaxies reside in more strongly clustered massive haloes.
 2.
We consider the redshift and luminosity evolution of halo parameters M_{1}, representing the characteristic halo mass required to host at least one satellite galaxy, and M_{min}, representing the halo mass where on average 50% of haloes contain one galaxy. For our full sample, M_{1} is closely approximated by ~17 × M_{min}. This ratio becomes smaller at higher redshifts for lowermass haloes. For red galaxies, M_{1}/M_{min} is ~12 at intermediate mass scale, and remains constant for all halo masses over 0.2 < z < 0.8. This lower ratio for redder galaxies reflects the higher abundance of satellites in these haloes; in general, redder populations reside in more massive haloes.
 3.
By fitting a simple analytic relation between central galaxy luminosity and halo mass to our observations, we find a maximum in the stellartohalo mass ratio. In the lowest redshift bin z ~ 0.3, this peak mass is for the full samples, and for the red samples, respectively. This transition represents the halo masses where baryons were most efficiently converted into stars, and is in good agreement with measurements from other studies.
 4.
For the full sample, shifts to higher halo masses at higher redshifts. For the red galaxy sample, the peak position evolves less rapidly with redshift. These results can be understood qualitatively from the lack of ongoing starformation in the red galaxy population which means that the stellar mass content in these massive haloes changes very slowly. For the full galaxy sample, which is expected to contain galaxies still undergoing active starformation, the shift of the transition mass to progressively more massive haloes at higher redshifts is a manifestation of “antihierarchical” galaxy formation. These results also indicate that the massive, passive galaxies in our survey are already fully formed by z ~ 1.2.
 5.
At increasing luminosities, the fraction of satellite galaxies rapidly drops to zero. For less luminous galaxies, the satellite fraction is higher in our red sample than in our full sample, reflecting the high number of faint red satellites, consistent with observations in the local Universe (Mandelbaum et al. 2006; Zehavi et al. 2011). For our full sample, the number of satellites increases by a factor of two from z ~ 1 to z ~ 0, consistent with a combined study of galaxies in DEEP2 and SDSS (Zheng et al. 2007). This result suggests that the occupation function for satellite galaxies inside intermediatemass haloes (~10^{12} h^{1} M_{⊙}) remains constant with cosmic time, and implies that the galaxy merging rate does not change since z = 1.
 6.
Above the characteristic sample luminosity, the mean galaxy bias b_{g} increases rapidly. Faint red samples have b_{g} ~ 1.5, compared to b_{g} ~ 1 for the full sample. Redder samples also have higher average halo masses, which decrease towards higher redshifts.
 7.
We compare the measured evolution of our mean halo masses, biases and satellite fractions with a model assuming constant HOD parameters over our redshift range. This model provides a remarkably good fit to our measurements for massive haloes. We note that lower mean halo masses (particularly in the full sample) show a slower decrease compared to a redshiftindependent HOD model.
Systematic photometric redshift errors have little impact on our faint luminositythreshold samples, for which we find no evidence for incompleteness larger than the fieldtofield variance estimates. Our very brightest samples however may be affected by incompleteness or contamination from fainter objects. We have shown that these errors could reduce our measured clustering amplitudes and consequently could lead to an underestimation of halo masses. The effect will be more important in contaminated samples, where the higher galaxy number density will lead to lower halo masses. In the worst case, log M_{min} might be underestimated by ~0.5.
Because of the lack of deep nearinfrared data in the CFHTLS, we cannot calculate reliable stellar masses and cleanly separate the active and passive populations in massselected samples. Such a data set, perhaps combined with approximate local density estimators calibrated using spectroscopic redshifts, and a refined halo model would allow us to understand in greater detail the physical origin of the shift of the transition masses to higher halo masses observed here. Furthermore, given the fact that the peak of the starformation and mass assembly takes place in the redshift range 1 < z < 2, it is important to push these studies to higher redshifts by constructing larger, widearea surveys with deep nearinfrared data.
In our model, we have assumed a halo mass function, halo density profile, and a cosmological model, although it is well known that the evolution of galaxy clustering depends on these parameters. Combining galaxygalaxy lensing, galaxy clustering and massfunction measurements from an expanded CFHTLS including nearinfrared data should provide new insights on both cosmology and galaxy evolution.
Acknowledgments
We acknowledge CFHT, TERAPIX, CADC and the CFHTLS Steering Group for their assistance in planning, executing and reducing the CFHTLS survey. We also acknowledge the COSMOS, VVDS and DEEP2 collaborations for making their data publicly available. We thank T. Hamana, S. Colombi, M. Brown, M. White, A. Ross, A. Leauthaud and P. Behroozi for useful comments. We thank the referee for his careful and detailed report. HJMCC acknowledges support from ANR grant “ANR07BLAN0228”. J.C. is supported by the Japanese Society for the Promotion of Science.
References
 Abbas, U., de la Torre, S., Le Fèvre, O., et al. 2010, MNRAS, 406, 1306 [NASA ADS] [Google Scholar]
 Arnouts, S., Moscardini, L., Vanzella, E., et al. 2002, MNRAS, 329, 355 [NASA ADS] [CrossRef] [Google Scholar]
 Auger, M. W., Treu, T., Bolton, A. S., et al. 2010, ApJ, 724, 511 [NASA ADS] [CrossRef] [Google Scholar]
 Behroozi, P. S., Conroy, C., & Wechsler, R. H. 2010, ApJ, 717, 379 [NASA ADS] [CrossRef] [Google Scholar]
 Benítez, N. 2000, ApJ, 536, 571 [NASA ADS] [CrossRef] [Google Scholar]
 Benjamin, J., van Waerbeke, L., Ménard, B., & Kilbinger, M. 2010, MNRAS, 408, 1168 [NASA ADS] [CrossRef] [Google Scholar]
 Benson, A. J., Bower, R. G., Frenk, C. S., et al. 2003, ApJ, 599, 38 [NASA ADS] [CrossRef] [Google Scholar]
 Berlind, A. A., & Weinberg, D. H. 2002, ApJ, 575, 587 [NASA ADS] [CrossRef] [Google Scholar]
 Berlind, A. A., Weinberg, D. H., Benson, A. J., et al. 2003, ApJ, 593, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Berlind, A. A., Blanton, M. R., Hogg, D. W., et al. 2005, ApJ, 629, 625 [Google Scholar]
 Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Blake, C., Collister, A., & Lahav, O. 2008, MNRAS, 385, 1257 [NASA ADS] [CrossRef] [Google Scholar]
 Bower, R. G., Benson, A. J., Malbon, R., et al. 2006, MNRAS, 370, 645 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, M. J. I., Zheng, Z., White, M., et al. 2008, ApJ, 682, 937 [NASA ADS] [CrossRef] [Google Scholar]
 Bundy, K., Fukugita, M., Ellis, R. S., et al. 2009, ApJ, 697, 1369 [NASA ADS] [CrossRef] [Google Scholar]
 Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682 [NASA ADS] [CrossRef] [Google Scholar]
 Cappé, O., Guillin, A., Marin, J.M., & Robert, C. 2004, J. Comput. Graph. Statist., 13, 907 [CrossRef] [MathSciNet] [Google Scholar]
 Cappé, O., Douc, R., Guillin, A., Marin, J.M., & Robert, C. 2008, Statist. Comput., 18, 447 [CrossRef] [Google Scholar]
 Cattaneo, A., Faber, S. M., Binney, J., et al. 2009, Nature, 460, 213 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Coleman, G. D., Wu, C.C., & Weedman, D. W. 1980, ApJS, 43, 393 [NASA ADS] [CrossRef] [Google Scholar]
 Conroy, C., & Wechsler, R. H. 2009, ApJ, 696, 620 [NASA ADS] [CrossRef] [Google Scholar]
 Cooray, A., & Sheth, R. 2002, Phys. Rep., 372, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Coupon, J., Ilbert, O., Kilbinger, M., et al. 2009, A&A, 500, 981 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 112, 839 [Google Scholar]
 Croton, D. J., Springel, V., White, S. D. M., et al. 2006, MNRAS, 365, 11 [NASA ADS] [CrossRef] [Google Scholar]
 Croton, D. J., Gao, L., & White, S. D. M. 2007, MNRAS, 374, 1303 [NASA ADS] [CrossRef] [Google Scholar]
 Davis, M., Faber, S. M., Newman, J., et al. 2003, in Discoveries and Research Prospects from 6 to 10MeterClass Telescopes II, ed. P. Guhathakurta, Puragra, Proc. SPIE, 4834, 161 [Google Scholar]
 de la Torre, S., Guzzo, L., Kovač, K., et al. 2010, MNRAS, 409, 867 [NASA ADS] [CrossRef] [Google Scholar]
 Faber, S. M., Willmer, C. N. A., Wolf, C., et al. 2007, ApJ, 665, 265 [NASA ADS] [CrossRef] [Google Scholar]
 Fall, S. M., & Efstathiou, G. 1980, MNRAS, 193, 189 [NASA ADS] [CrossRef] [Google Scholar]
 Foucaud, S., Conselice, C. J., Hartley, W. G., et al. 2010, MNRAS, 406, 147 [NASA ADS] [CrossRef] [Google Scholar]
 Garilli, B., Le Fèvre, O., Guzzo, L., et al. 2008, A&A, 486, 683 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Goranova, Y., Hudelot, P., Magnard, F., et al. 2009, The CFHTLS T0006 Release [Google Scholar]
 Groth, E. J., & Peebles, P. J. E. 1977, ApJ, 217, 385 [Google Scholar]
 Gunn, J. E., & Gott, III, J. R. 1972, ApJ, 176, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Guo, Q., White, S., Li, C., & BoylanKolchin, M. 2010, MNRAS, 404, 1111 [NASA ADS] [Google Scholar]
 Hamilton, A. J. S. 2000, MNRAS, 312, 257 [NASA ADS] [CrossRef] [Google Scholar]
 Hartlap, J., Simon, P., & Schneider, P. 2007, A&A, 464, 399 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Henry, J. P. 2000, ApJ, 534, 565 [Google Scholar]
 Heymans, C., Bell, E. F., Rix, H.W., et al. 2006, MNRAS, 371, L60 [Google Scholar]
 High, F. W., Stubbs, C. W., Rest, A., Stalder, B., & Challis, P. 2009, AJ, 138, 110 [NASA ADS] [CrossRef] [Google Scholar]
 Hoaglin, D. C., Mosteller, F., & Tukey, J. W. 1983, Understanding robust and exploratory data anlysis, ed. D. C. Hoaglin, F. Mosteller, & J. W. Tukey [Google Scholar]
 Hopkins, P. F., Somerville, R. S., Hernquist, L., et al. 2006, ApJ, 652, 864 [Google Scholar]
 Hu, W., & Kravtsov, A. V. 2003, ApJ, 584, 702 [NASA ADS] [CrossRef] [Google Scholar]
 Ilbert, O., Tresse, L., Zucca, E., et al. 2005, A&A, 439, 863 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ilbert, O., Arnouts, S., McCracken, H. J., et al. 2006, A&A, 457, 841 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ilbert, O., Capak, P., Salvato, M., et al. 2009, ApJ, 690, 1236 [NASA ADS] [CrossRef] [Google Scholar]
 Ilbert, O., Salvato, M., Le Floc’h, E., et al. 2010, ApJ, 709, 644 [NASA ADS] [CrossRef] [Google Scholar]
 Kilbinger, M., Wraith, D., Robert, C. P., et al. 2010, MNRAS, 405, 2381 [NASA ADS] [Google Scholar]
 Kilbinger, M., Benabed, K., Cappe, O., et al. 2011 [arXiv:1101.0950] [Google Scholar]
 Kinney, A. L., Calzetti, D., Bohlin, R. C., et al. 1996, ApJ, 467, 38 [NASA ADS] [CrossRef] [Google Scholar]
 Kitayama, T., & Suto, Y. 1996, ApJ, 469, 480 [Google Scholar]
 Kravtsov, A. V., Berlind, A. A., Wechsler, R. H., et al. 2004, ApJ, 609, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Landy, S. D., & Szalay, A. S. 1993, ApJ, 412, 64 [NASA ADS] [CrossRef] [Google Scholar]
 Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 692 [NASA ADS] [CrossRef] [Google Scholar]
 Le Fèvre, O., Vettolani, G., Garilli, B., et al. 2005, A&A, 439, 845 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Leauthaud, A., Tinker, J., Behroozi, P. S., Busha, M. T., & Wechsler, R. 2012a, ApJ, 738, 45 [Google Scholar]
 Leauthaud, A., Tinker, J., Bundy, K., et al. 2012b, ApJ, 744, 159 [NASA ADS] [CrossRef] [Google Scholar]
 Limber, D. N. 1954, ApJ, 119, 655 [NASA ADS] [CrossRef] [Google Scholar]
 Mandelbaum, R., Seljak, U., Kauffmann, G., Hirata, C. M., & Brinkmann, J. 2006, MNRAS, 368, 715 [Google Scholar]
 McCracken, H. J., Peacock, J. A., Guzzo, L., et al. 2007, ApJS, 172, 314 [NASA ADS] [CrossRef] [Google Scholar]
 Meneux, B., Guzzo, L., de la Torre, S., et al. 2009, A&A, 505, 463 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 More, S., van den Bosch, F. C., Cacciato, M., et al. 2011, MNRAS, 410, 210 [NASA ADS] [CrossRef] [Google Scholar]
 Moster, B. P., Somerville, R. S., Maulbetsch, C., et al. 2010, ApJ, 710, 903 [NASA ADS] [CrossRef] [Google Scholar]
 Nakamura, T. T., & Suto, Y. 1997, Progr. Theor. Phys., 97, 49 [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [NASA ADS] [CrossRef] [Google Scholar]
 Norberg, P., Baugh, C. M., Hawkins, E., et al. 2001, MNRAS, 328, 64 [Google Scholar]
 Peng, Y.j., Lilly, S. J., Kovač, K., et al. 2010, ApJ, 721, 193 [NASA ADS] [CrossRef] [Google Scholar]
 Peng, Y., Lilly, S. J., Renzini, A., & Carollo, M. 2011, ApJ, submitted [arXiv:1106.2546] [Google Scholar]
 Peterson, J. R., & Fabian, A. C. 2006, Phys. Rep., 427, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Phleps, S., Peacock, J. A., Meisenheimer, K., & Wolf, C. 2006, A&A, 457, 145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 [NASA ADS] [CrossRef] [Google Scholar]
 Robaina, A. R., Bell, E. F., van der Wel, A., et al. 2010, ApJ, 719, 844 [NASA ADS] [CrossRef] [Google Scholar]
 Sheth, R. K., & Tormen, G. 1999, MNRAS, 308, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Sheth, R. K., Mo, H. J., & Tormen, G. 2001, MNRAS, 323, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Simon, P., Hetterscheidt, M., Wolf, C., et al. 2009, MNRAS, 398, 807 [NASA ADS] [CrossRef] [Google Scholar]
 Skibba, R. A., & Sheth, R. K. 2009, MNRAS, 392, 1080 [NASA ADS] [CrossRef] [Google Scholar]
 Somerville, R. S., Hopkins, P. F., Cox, T. J., Robertson, B. E., & Hernquist, L. 2008, MNRAS, 391, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Springel, V., Frenk, C. S., & White, S. D. M. 2006, Nature, 440, 1137 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Szalay, A. S., Connolly, A. J., & Szokoly, G. P. 1999, AJ, 117, 68 [Google Scholar]
 Takada, M., & Jain, B. 2003, MNRAS, 344, 857 [NASA ADS] [CrossRef] [Google Scholar]
 Tinker, J. L., Weinberg, D. H., Zheng, Z., & Zehavi, I. 2005, ApJ, 631, 41 [NASA ADS] [CrossRef] [Google Scholar]
 van Dokkum, P. G., Whitaker, K. E., Brammer, G., et al. 2010, ApJ, 709, 1018 [NASA ADS] [CrossRef] [Google Scholar]
 Wake, D. A., Whitaker, K. E., Labbé, I., et al. 2011, ApJ, 728, 46 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, L., & Jing, Y. P. 2010, MNRAS, 402, 1796 [NASA ADS] [CrossRef] [Google Scholar]
 Weinberg, N. N., & Kamionkowski, M. 2003, MNRAS, 341, 251 [NASA ADS] [CrossRef] [Google Scholar]
 White, S. D. M., & Rees, M. J. 1978, MNRAS, 183, 341 [NASA ADS] [CrossRef] [Google Scholar]
 White, M., Zheng, Z., Brown, M. J. I., Dey, A., & Jannuzi, B. T. 2007, ApJ, 655, L69 [NASA ADS] [CrossRef] [Google Scholar]
 Willmer, C. N. A., Faber, S. M., Koo, D. C., et al. 2006, ApJ, 647, 853 [NASA ADS] [CrossRef] [Google Scholar]
 Wraith, D., Kilbinger, M., Benabed, K., et al. 2009, Phys. Rev. D, 80, 023507 [NASA ADS] [CrossRef] [Google Scholar]
 Zehavi, I., Zheng, Z., Weinberg, D. H., et al. 2005, ApJ, 630, 1 [Google Scholar]
 Zehavi, I., Zheng, Z., Weinberg, D. H., et al. 2011, ApJ, 736, 59 [NASA ADS] [CrossRef] [Google Scholar]
 Zehavi, I., Patiri, S., & Zheng, Z. 2012, ApJ, 746, 145 [NASA ADS] [CrossRef] [Google Scholar]
 Zheng, Z. 2004, ApJ, 610, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Zheng, Z., Berlind, A. A., Weinberg, D. H., et al. 2005, ApJ, 633, 791 [NASA ADS] [CrossRef] [Google Scholar]
 Zheng, Z., Coil, A. L., & Zehavi, I. 2007, ApJ, 667, 760 [NASA ADS] [CrossRef] [Google Scholar]
 Zheng, Z., Zehavi, I., Eisenstein, D. J., Weinberg, D. H., & Jing, Y. P. 2009, ApJ, 707, 554 [NASA ADS] [CrossRef] [Google Scholar]
 Zucca, E., Ilbert, O., Bardelli, S., et al. 2006, A&A, 455, 879 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Online material
Appendix A: The halo model
The halo occupation distribution (HOD) model (Berlind & Weinberg 2002; Kravtsov et al. 2004; Zheng et al. 2005) allows one to derive physical information about galaxy populations and the dark matter haloes which host them. The HOD prescription is based on the halo model, which describes how dark matter is distributed in space. In this framework, all of the matter is assumed to reside in virialised haloes. The HOD parametrisation specifies how many galaxies populate haloes, on average, as function of halo mass. Accordingly, the number of galaxies per halo N only depends on the mass M of the halo.
A.1. Darkmatter halo model
The three ingredients to the halo model of dark matter are the halo mass function, the halo profile and the halo bias. The dark matter halo abundance can be inferred using the Press & Schechter (1974) approach where dark matter collapses into overdense regions above the critical density δ_{c}, linearly evolved to z = 0. The mass function, which is the halo number density per unit mass, can be parametrized as (A.1)where is the mean density of matter at the present day. The new mass variable ν writes (A.2)and characterizes the peak heights of the density field as function of mass and redshift. The linear critical density δ_{c} depends on the adopted cosmology and redshift; we use the fitting formula from Weinberg & Kamionkowski (2003); see also Kitayama & Suto (1996), D(z) is the linear growth factor at redshift z, and σ(M) is the rms of density fluctuations in a tophat filter of width , computed from linear theory, (A.5)where W(x) = (3/x^{3})[sinx − xcosx] . For the mass function f(ν), we choose the parameterisation by Sheth & Tormen (1999), calibrated on simulations: (A.6)where the normalisation A is fixed by imposing: (A.7)we adopt the values p = 0.3 and . If not indicated otherwise, all integrals over the mass function are performed between M_{low} = 10^{3} h^{1} M_{⊙} and M_{high} = 10^{16} h^{1} M_{⊙}.
We describe the halo density profile by the following form Navarro et al. (1997), (A.8)The total halo mass is then written as (see Takada & Jain 2003): (A.9)where c = r_{vir}/r_{s} is the “concentration parameter”, for which we assume the following expression, (A.10)We take c_{0} = 11 and β = 0.13, and M_{ ⋆ } is defined such that ν(z = 0) = 1, i.e. δ_{c}(0) = σ(M_{ ⋆ }). The virial radius r_{vir} is given by the following relation (A.11)with Δ_{vir}(z) being the critical overdensity for virialisation at redshift z (Kitayama & Suto 1996; Nakamura & Suto 1997; Henry 2000). We take the fitting formula from Weinberg & Kamionkowski (2003)(A.12)Following Tinker et al. (2005), we use the scaledependent halo bias (A.13)where ξ_{m} is the matter correlation function. The largescale halo bias b_{h}(M,z) is given by Sheth et al. (2001) as (A.14)As in Tinker et al. (2005), we adopt the revised parameters , b = 0.35 and c = 0.8.
A.2. The galaxy correlation function
We write the correlation function as a sum of two components: The onehalo term, to express the galaxy correlation inside a halo, and the twohalo term, to account for halotohalo correlation, (A.15)the onehalo term depends on the number of galaxy pairs ⟨ N(N − 1) ⟩ per halo. This is comprised of the centralsatellite contribution ⟨ N_{c}N_{s} ⟩ and the satellitesatellite term ⟨ N_{s}(N_{s} − 1) ⟩ . Assuming a Poisson distribution, we write (A.16)The onehalo term correlation function for centralsatellite pairs is then (A.17)The lower integration limit M_{vir}(r) is the virial mass contained in a halo of radius r, computed with Eq. (A.11). This accounts for the fact that lessmassive haloes are too small to contribute to the correlation at separation r.
The onehalo satellitesatellite contribution ξ_{ss} involves the halo profile autoconvolution and is therefore easier to compute in Fourier space. The corresponding power spectrum is written as (A.18)where u_{h}(kM) is the Fourier transform of the darkmatter halo profile ρ_{h}(rM). The correlation function ξ_{ss} is then obtained via a Fourier transform. The onehalo correlation function is the sum of the two contributions, (A.19)The twohalo term is derived from the darkmatter power spectrum and the halo twopoint correlation function: (A.20)where (A.21)The upper integration limit M_{lim}(r) takes into account the halo exclusion (Zheng 2004), i.e. the fact that haloes are
nonoverlapping. We follow Tinker et al. (2005) to compute M_{lim}(r) by matching with the following expression, (A.22)where P(r,M_{1},M_{2}) is the probability that two ellipsoidal haloes of mass M_{1} and M_{2}, respectively, do not overlap. Defining x = r/[r_{vir}(M_{1}) + r_{vir}(M_{2})] as the ratio of the halo separation and the sum of the virial radii, and y = (x − 0.8)/0.29, Tinker et al. (2005) found the probability of nonoverlapping haloes to be (A.23)We Fouriertransform Eq. (A.20) for a range of tabulated values of r, to compute the twohalo term ξ_{2} of the galaxy autocorrelation function. Finally, we renormalise it to the total number of galaxy pairs: (A.24)The angular twopoint correlation function w(θ) is computed from the observed photometric redshift distribution and ξ(r) using Limber’s equation (Limber 1954): (A.25)with (A.26)and x(z), the radial comoving coordinate.
Appendix B: Bestfitting HOD parameters and deduced quantities
Description of all galaxy samples and bestfitting HOD parameters. Halo masses are given in h^{1} M_{⊙} and galaxy number densities in h^{3} Mpc^{3}.
Description of red galaxy samples and bestfitting HOD parameters. Halo masses are given in h^{1}M_{⊙} and galaxy number densities in h^{3} Mpc^{3}.
Description of blue galaxy samples. Galaxy number densities are given in h^{3} Mpc^{3}.
Appendix C: Twopoint correlation function measurements
Twopoint correlation function measurements in the range 0.2 < z < 0.4.
Twopoint correlation function measurements in the range 0.4 < z < 0.6.
Twopoint correlation function measurements in the range 0.6 < z < 0.8.
Twopoint correlation function measurements in the range 0.8 < z < 1.0.
Twopoint correlation function measurements in the range 1.0 < z < 1.2.
Fig. C.1 1D (diagonal) and 2D likelihood distributions of bestfitting HOD parameters for the full sample, in the range 0.4 < z < 0.6, and M_{g} − 5log h < −19.8. 
Bestfitting parameters (from Eq. (25)) of M_{min} and M_{1} as function of luminosity threshold corrected for passive redshift evolution to approximate stellar mass selected samples. Results are given for all and red samples, as function of redshift bins. Halo masses are given in h^{1} M_{⊙}.
All Tables
Description of all galaxy samples and bestfitting HOD parameters. Halo masses are given in h^{1} M_{⊙} and galaxy number densities in h^{3} Mpc^{3}.
Description of red galaxy samples and bestfitting HOD parameters. Halo masses are given in h^{1}M_{⊙} and galaxy number densities in h^{3} Mpc^{3}.
Description of blue galaxy samples. Galaxy number densities are given in h^{3} Mpc^{3}.
Bestfitting parameters (from Eq. (25)) of M_{min} and M_{1} as function of luminosity threshold corrected for passive redshift evolution to approximate stellar mass selected samples. Results are given for all and red samples, as function of redshift bins. Halo masses are given in h^{1} M_{⊙}.
All Figures
Fig. 1 Sky coverage of W1, W2, W3 and W4 fields. Each blue square represents one MegaCam pointing. Spectroscopic data (VVDS F02, F22, and DEEP2) are in red. CFHTLS Deep fields overlapping with the wide survey (D1 and D3) are shown in green. 

In the text 
Fig. 2 Photometric redshift accuracy in the 0.2 < photo  z < 1.2 interval. The figure shows the number density of photometric redshifts in W1, W3 and W4, versus spectroscopic redshifts from VVDS and DEEP2 surveys. The blue line is a linear fit to photoz versus spectroz. 

In the text 
Fig. 3 Top: type selection based on best fitting templates for the “red” sample (left) and the “blue” sample (right) in the CFHTLS Wide. We show the colour distribution (M_{u} − M_{r}) as function of absolute magnitude (M_{g}) and redshift (top to bottom) and number counts for “red” and “blue” objects on the right panels. Bottom: number of red and blue galaxies as function of magnitude and redshift. 

In the text 
Fig. 4 Sample selection in the full (left), red (centre) and blue (right) galaxy samples. In each panel, the galaxy number density in the plane M_{g}/z is shown; red rectangles represent the luminosity threshold samples. 

In the text 
Fig. 5 Redshift distributions for the full (left), red (centre) and blue (right) galaxy samples as function of absolute magnitude threshold. 

In the text 
Fig. 6 Cross correlation analysis between redshift bins for the full sample. Top: autocorrelation in the redshift bin 0.4 < z < 0.6 (straight line) and crosscorrelation between the bins 0.4 < z < 0.6 and 0.8 < z < 1.0 (dotted line) and between 0.4 < z < 0.6 and 1.0 < z < 1.2 (dashed line). Bottom: quantitative estimates of the contamination (percentage of galaxies scattered) from a pairwise analysis between redshift bins 0.4 < z < 0.6 and 0.8 < z < 1.0 (left) and between 0.4 < z < 0.6 and 1.0 < z < 1.2 (right). The contours show the 68.3 (blue), 95.5 (green) and 99.7 (red) confidence regions. 

In the text 
Fig. 7 Covariance matrix correlation coefficients of w(θ) from Jackknife estimates, for galaxy samples with 0.4 < z < 0.6 and M_{g} −5log h < −19.8. Left panel: red sample. Right panel: blue sample. 

In the text 
Fig. 8 Example of a measured w(θ) (for all galaxies in the redshift range 0.4 < z < 0.6 and for M_{g} − 5log h < −19.8), as well as the bestfitting model, as described in Sect. 4. Left: w(θ) measurement and model. Right: N(M), showing the central term N_{cent} and the satellite term N_{s}. 

In the text 
Fig. 9 w(θ) measurements (left) with their best fit models for all galaxies, and the resulting galaxy occupation function (right) with respect to halo mass. We illustrate the redshift evolution (top) in the luminosity threshold M_{g} − 5log h < − 19.8, and the luminosity dependence (bottom) in the redshift range 0.4 < z < 0.6. 

In the text 
Fig. 10 w(θ) measurements (left) with their best fit models for red galaxies, and the resulting galaxy occupation function (right) with respect to halo mass. We illustrate the redshift evolution (top) in the luminosity threshold M_{g} − 5log h < − 19.8, and the luminosity dependence (bottom) in the redshift range 0.4 < z < 0.6. 

In the text 
Fig. 11 w(θ) measurements for blue galaxies. We illustrate the redshift evolution (left) at the luminosity threshold M_{g} − 5log h < − 19.8, and the luminosity dependence (right) in the redshift range 0.4 < z < 0.6. 

In the text 
Fig. 12 Top panels: halo mass estimates M_{min} and M_{1} for all (left) and red (right) galaxy samples as function of luminosity threshold (M_{g}). Bottom panels: M_{1} versus M_{min} in different redshift bins. The dotted lines represent a linear relation and the dashed lines are a powerlaw fit to M_{1} versus M_{min} in each redshift bin. 

In the text 
Fig. 13 M_{min} as a function of galaxy number density for all (left) and red (right) galaxies. The dashed line in each panel represents a power law fit to the data, and the dotted line represents a power law fit to the results given in Brown et al. (2008) (red galaxies) and in the SDSS by Zehavi et al. (2011) (all galaxies). We converted M_{200} halo masses to virial masses M_{vir} when necessary. 

In the text 
Fig. 14 Stellar masstoluminosity relations for red (left) and blue (right) galaxies, as a function of stellar mass and redshift as measured in the COSMOS survey (Ilbert et al. 2009). Each panel shows a separate slice in stellar mass. Samples are fitted with a power law (dashed line). 

In the text 
Fig. 15 Left panels: halo mass estimates M_{min} and M_{1} for all (top) and red (bottom) galaxy samples, as function of luminosity threshold, corrected for passive redshift evolution to approximate stellar mass selected samples. The dashed lines correspond to Eq. (25). Right panels: lighttohalo mass ratios (Eq. (26) with identical parameters as those fitted with Eq. (25)) as function of halo mass. 

In the text 
Fig. 16 Position of the peak corresponding to for the full (filled circles) and red (filled diamonds) samples as function of redshift. Error bars are derived from the nearest data point to the peak. We compare our results with measurements in SDSS (Zehavi et al. 2011), in COSMOS (Leauthaud et al. 2012b) and in the NEWFIRM Medium Band Survey (NMBS, Wake et al. 2011). We convert M_{200} halo masses (corresponding to a sphere with overdensity 200 times the background density) into virial masses M_{vir} (our definition) using the method described in Appendix C of Hu & Kravtsov (2003). We also display numerical predictions based on abundancematching method from Behroozi et al. (2010) and Moster et al. (2010), where the shaded region shows the 1 − σ uncertainty. 

In the text 
Fig. 17 Satellite fractions as function of luminosity and redshift, for all and red samples (top and bottom panels). The left panels show the measured satellite fraction as function of corrected luminosity. Dashed lines in the middle panels show satellite fractions extrapolated in redshift, based on our HOD parameters measured in the redshift range 0.4 < z < 0.6. We compare our results with the literature for all galaxies at z = 0.1 and z = 1.0 (top, Zheng et al. 2007) and with galaxygalaxy lensing estimates at z = 0.1 (bottom, Mandelbaum et al. 2006) for red galaxies. In the right panels we provide values extrapolated in redshift for several L′/L_{∗} ratios. 

In the text 
Fig. 18 Galaxy bias as function of luminosity and redshift. Left panels: bias for full and red samples (top and bottom panels) as function of corrected luminosity threshold L′. The dashed line represents Eq. (29). Right panels: computed bias for several luminosity threshold samples as function of redshift. Error bars were derived from the parametric form shown on the left (dashed line), whereas the dotted lines show the bias evolution computed from our model, assuming a redshiftindependent HOD parameter set. 

In the text 
Fig. 19 Mean halo mass as function of luminosity and redshift. Left panels: mean halo mass for full and red samples (top and bottom panels) as function of corrected luminosity threshold L′. The dashed line represents Eq. (30). Right panels: computed mean halo mass for several luminosity threshold samples as function of redshift. Error bars were derived from the parametric form shown on the left (dashed line), whereas the dotted lines show the mean halo mass evolution computed from our model, assuming a redshiftindependent HOD parameter set. 

In the text 
Fig. C.1 1D (diagonal) and 2D likelihood distributions of bestfitting HOD parameters for the full sample, in the range 0.4 < z < 0.6, and M_{g} − 5log h < −19.8. 

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