CANDELS Multi-wavelength Catalogs: Source
Identification and Photometry in the CANDELS
COSMOS Survey Field

H. Nayyeri, S. Hemmati, B. Mobasher, H. C. Ferguson, A. Cooray, G. Barro, S. M. Faber, M. Dickinson, A. M. Koekemoer, M. Peth, et al.

► To cite this version:

H. Nayyeri, S. Hemmati, B. Mobasher, H. C. Ferguson, A. Cooray, et al. CANDELS Multi-wavelength Catalogs: Source Identification and Photometry in the CANDELS COSMOS Survey Field. The Astrophysical Journal Supplement Series, 2017, 228, 10.3847/1538-4365/228/1/7. insu-03666222

HAL Id: insu-03666222

https://insu.hal.science/insu-03666222v1

Submitted on 4 Sep 2025

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

CANDELS MULTI-WAVELENGTH CATALOGS: SOURCE IDENTIFICATION AND PHOTOMETRY IN THE CANDELS COSMOS SURVEY FIELD

H. NAYYERI<sup>1,2</sup>, S. HEMMATI<sup>2,3</sup>, B. MOBASHER<sup>2</sup>, H. C. FERGUSON<sup>4</sup>, A. COORAY<sup>1</sup>, G. BARRO<sup>5,6</sup>, S. M. FABER<sup>5</sup>, M. DICKINSON<sup>7</sup>, A. M. KOEKEMOER<sup>4</sup>, M. PETH<sup>8</sup>, M. SALVATO<sup>9</sup>, M. L. N. ASHBY<sup>10</sup>, B. DARVISH<sup>2,11</sup>, J. DONLEY<sup>12</sup>, M. DURBIN<sup>4</sup>, S. FINKELSTEIN<sup>13</sup>, A. FONTANA<sup>14</sup>, N. A. GROGIN<sup>4</sup>, R. GRUETZBAUCH<sup>15</sup>, K. HUANG<sup>16</sup>, A. A. KHOSTOVAN<sup>2</sup>, D. KOCEVSKI<sup>17</sup>, D. KODRA<sup>18</sup>, B. LEE<sup>19</sup>, J. NEWMAN<sup>18</sup>, C. PACIFICI<sup>20,30</sup>, J. PFORR<sup>21,22</sup>, M. STEFANON<sup>23</sup>, T. WIKLIND<sup>24</sup>, S. P. WILLNER<sup>10</sup>, S. WUYTS<sup>25</sup>, M. CASTELLANO<sup>14</sup>, C. CONSELICE<sup>26</sup>, T. DOLCH<sup>27</sup>, J. S. DUNLOP<sup>28</sup>, A. GALAMETZ<sup>9</sup>, N. P. HATHI<sup>21</sup>, R. A. LUCAS<sup>4</sup>, and H. YAN<sup>29</sup>

<sup>1</sup>Department of Physics and Astronomy, University of California Irvine, Irvine, CA 92697, USA

<sup>2</sup>Department of Physics and Astronomy, University of California Riverside, Riverside, CA 92521, USA

<sup>3</sup>Infrared Processing and Analysis Center, California Institute of Technology, MS 100-22, Pasadena, CA 91125, USA

<sup>4</sup>Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA

<sup>5</sup>UCO/Lick Observatory, Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064, USA

<sup>6</sup>Department of Astronomy, University of California at Berkeley, Berkeley, CA 94720-3411, USA

<sup>7</sup>National Optical Astronomy Observatory, Tucson, AZ 85719, USA

<sup>8</sup>Department of Physics and Astronomy, The Johns Hopkins University, 366 Bloomberg Center, Baltimore, MD 21218, USA

<sup>9</sup>Max-Planck-Institut für Extraterrestrische Physik, Giessenbachstrasse 1, D-85748 Garching bei Munchen, Germany

<sup>10</sup>Harvard Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA

<sup>11</sup>Cahill Center for Astrophysics, California Institute of Technology, 1216 East California Boulevard, Pasadena, CA 91125, USA

<sup>12</sup>Los Alamos National Laboratory, Los Alamos, NM 87544, USA

<sup>13</sup>Department of Astronomy, The University of Texas at Austin, Austin, TX 78712, USA

<sup>14</sup>INAF—Osservatorio Astronomico di Roma, via Frascati 33, I-00040 Monte Porzio Catone (RM), Italy

<sup>15</sup>Center for Astronomy and Astrophysics, Observatorio Astronomico de Lisboa, Tapada da Ajuda, 1349-018 Lisboa, Portugal

<sup>16</sup>Department of Physics, University of California, Davis, CA 95616, USA

<sup>17</sup>Colby College, 4000 Mayflower Hill, Waterville, ME 04901, USA

<sup>18</sup>Department of Physics and Astronomy and PITT PACC, University of Pittsburgh, Pittsburgh, PA 15260, USA

<sup>19</sup>Department of Astronomy, University of Massachusetts, 710 North Plesant Street, Amherst, MA 01003, USA

<sup>20</sup>Goddard Space Flight Center, Code 665, Greenbelt, MD, USA

<sup>21</sup>Aix Marseille Université, CNRS, LAM (Laboratoire d'Astrophysique de Marseille) UMR 7326, F-13388, Marseille, France

<sup>22</sup>ESA/ESTEC SCI-S, Keplerlaan 1, 2201 AZ, Noordwijk, The Netherlands

<sup>23</sup>Huygens Laboratory Niels Bohrweg 2 2333 Ca Leiden, Netherlands

<sup>24</sup>Department of Physics, Catholic University of America, 20064 Washington Dc, USA

<sup>25</sup>Department of Physics, University of Bath, Claverton Down, Bath, BA2 7AY, UK

<sup>26</sup>School of Physics and Astronomy, University of Nottingham, Nottingham, UK

<sup>27</sup>Department of Physics, Hillsdale College 33 E. College Street, Hillsdale, MI 49242

<sup>28</sup>Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ, UK

<sup>29</sup>Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211, USA

Received 2016 May 17; revised 2016 October 13; accepted 2016 November 24; published 2017 January 13

ABSTRACT

We present a multi-wavelength photometric catalog in the COSMOS field as part of the observations by the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. The catalog is based on *Hubble Space Telescope* Wide Field Camera 3 (HST/WFC3) and Advanced Camera for Surveys observations of the COSMOS field (centered at R.A.: 10h00m28s, Decl.: +02°12′21″). The final catalog has 38671 sources with photometric data in 42 bands from UV to the infrared (~0.3–8 µm). This includes broadband photometry from HST, CFHT, Subaru, the Visible and Infrared Survey Telescope for Astronomy, and *Spitzer* Space Telescope in the visible, near-infrared, and infrared bands along with intermediate- and narrowband photometry from Subaru and medium-band data from Mayall NEWFIRM. Source detection was conducted in the WFC3 F160W band (at 1.6 µm) and photometry is generated using the Template FITting algorithm. We further present a catalog of the physical properties of sources as identified in the HST F160W band and measured from the multi-band photometry by fitting the observed spectral energy distributions of sources against templates.

**Key words:** catalogs – galaxies: high-redshift – galaxies: photometry – methods: data analysis – techniques: image processing

**Supporting material:** machine-readable tables

1. INTRODUCTION

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS: PI. S. Faber and H. Ferguson; see Grogin et al. 2011 and Koekemoer et al. 2011) is the largest Multi-Cycle Treasury program ever approved on the *Hubble Space Telescope* (HST), with more than 900 orbits, and it was

designed to use deep observations by the Wide Field Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) instruments to study galaxy formation and evolution throughout cosmic time in five fields in many different bands. The observations were done by the HST/WFC3 as the main mode with ACS observations in parallel. The CANDELS images are publicly available, and multi-wavelength photometric catalogs are made available by the CANDELS team following the release of the

<footnote n="30"><sup>30</sup> NASA Postdoctoral Program Fellow.</footnote>
images. All CANDELS photometric catalogs were selected based on the WFC3 F160W band and this is the reference image for all the other HST and non-HST data. This provides a data set with consistent photometry and physical properties across all the fields targeted as part of the survey. The CANDELS catalogs for the first two observed fields of UDS and GOODS-S (Galametz et al. 2013; Guo et al. 2013) are already publicly available<fnref n="31" /> and the three fields of COSMOS (this work), EGS (M. Stefanon et al. 2017, in preparation), and GOODS-N (G. Barro et al. 2017, in preparation) are in progress. The CANDELS observations are aimed at achieving several major science goals that could only be attained with data at the depth and resolution of CANDELS. These include studying the most distant objects in the universe at the epoch of reionization in the cosmic dawn (e.g., Finkelstein et al. 2012a; Grazian et al. 2012; Yan et al. 2012; Lorenzoni et al. 2013; Oesch et al. 2013; Duncan et al. 2014; Bouwens et al. 2015; Caputi et al. 2015; Finkelstein et al. 2015; Giallongo et al. 2015; Mitchell-Wynne et al. 2015; Roberts-Borsani et al. 2016; Song et al. 2016), understanding galaxy formation and evolution during the peak epoch of star formation in the cosmic high noon (e.g., Bell et al. 2012; Bruce et al. 2012; Kocevski et al. 2012; Lee et al. 2013; Wuyts et al. 2013; Barro et al. 2014; Hemmati et al. 2014, 2015; Whitaker et al. 2014; Williams et al. 2015), and studying star formation from deep UV observations and cosmological studies from supernova observations (e.g., Jones et al. 2013; Teplitz et al. 2013; Rodney et al. 2014, 2016; Strolger et al. 2015). These main science goals are described in more detail by Grogin et al. (2011) and Koekemoer et al. (2011).

One of the major goals of modern observational cosmology is to study the formation and evolution of galaxies with cosmic time. Recent advances in this frontier have been enabled by the availability of observations in different wavelengths, targeting different populations of galaxies (e.g., York et al. 2000; Wolf et al. 2003; Bell et al. 2004; Giavalisco et al. 2004; Skrutskie et al. 2006; Faber et al. 2007; Lawrence et al. 2007; Scoville et al. 2007b; Ilbert et al. 2013; Khostovan et al. 2015; Shivaei et al. 2015; Hemmati et al. 2017; Vasei et al. 2016). The advent of HST benefited many such studies by making it possible to have the deepest observations of the sky in multiple bands. In particular the installation of the WFC3 on board HST initiated a new stage for studying galaxy evolution at new extremes (McLeod et al. 2015; Oesch et al. 2016).

Galaxy populations at different look-back times range from very blue star-forming galaxies to red dusty or very old systems. Understanding the evolution of these populations relies on the availability of multi-wavelength photometric data from the bluest to the reddest bands possible. The CANDELS multi-wavelength catalogs combine the best and deepest observations by HST with the deepest ground-based observations and Spitzer Space Telescope data (Galametz et al. 2013; Guo et al. 2013). These catalogs of tens of thousands of extragalactic sources, consistently measured across many bands from ~0.3 to 8 µm, bring a unique opportunity to study galaxy evolution.

The Cosmic Evolution Survey (COSMOS; Scoville et al. 2007b) centered at R.A.: $10^{h}00^{m}28^{s}$, Decl.: $+02^{\circ}12'21''$ is a 2 deg² field located near the celestial equator. It was initially picked to maximize the visibility from observatories from both hemispheres and was specifically chosen to avoid any bright

X-ray, UV, or radio sources (Scoville et al. 2007b) and to be large enough for studies of large-scale structure (e.g., Scoville et al. 2007a, 2013; Kovač et al. 2010; Darvish et al. 2014, 2015b).

The COSMOS field was targeted by CANDELS in a north-south strip, lying within the central ultra-deep strip of the UltraVISTA imaging (McCracken et al. 2012) and hence also the Spitzer SEDS imaging (Ashby et al. 2013) in order to ensure the best possible supporting data at longer near-infrared wavelengths. The HST observations cover an area of ≈216 arcmin² in the WFC3/IR with parallel ACS observations. The catalog was selected in the HST/WFC3 F160W band and has multi-band data for 38671 objects from ~0.3 to 8 µm. These fluxes are measured consistently across all these bands and in agreement with photometry measurement techniques adopted by all CANDELS catalogs. We use the unprecedented depth and resolution provided by the HST for measuring the flux of the faintest targets across all bands. By fitting this multi-waveband information with template libraries, we also measured photometric redshift and stellar mass for each object.

The paper is organized as follow. Section 2 presents the data used to compile the catalog. Section 3 describes our photometry of the HST, Spitzer, and ground-based bands. Section 4 is devoted to data quality checks for our TFIT photometry. Physical parameter estimation using the measured photometry is presented in Section 5. In Section 6 we investigate the applications of our deep photometry on studies of high-redshift star-forming and quiescent galaxies. We summarize our results in Section 7. Throughout this paper we assume a cosmological model with $H_0 = 70 \text{ km s}^{-1} \text{ Mpc}^{-1}$, $\Omega_m = 0.3$, and $\Omega_\Lambda = 0.7$. All magnitudes are in the AB system where $m_{AB} = 23.9 - 2.5 \times \log(f_v/1\mu Jy)$ (Oke & Gunn 1983).

The CANDELS COSMOS photometry catalog will be publicly available through the CANDELS Web site<fnref n="32" /> along with the physical properties estimates and all the documentations. These will also be available on the Mikulski Archive for Space Telescopes (MAST),<fnref n="33" /> via the online version of the catalog and through Centre de Données astronomiques de Strasbourg (CDS). We will also make these data available through the Rainbow Database<fnref n="34" /> (Pérez-González et al. 2008; Barro et al. 2011).

## 2. DATA

The COSMOS field (Scoville et al. 2007b) has been observed in many different wavelengths from the X-ray to the far-infrared. There are observations from the CFHT/MegaPrime in the $u^*$, $g^*$, $r^*$, $i^*$, and $z^*$ (Gwyn 2012), from the Subaru/Suprime-Cam in the B, $g^+$, V, $r^+$, $i^+$, and $z^+$ (Taniguchi et al. 2007), from the VLT/Visible and Infrared Survey Telescope for Astronomy (VISTA) in the Y, J, H, and $K_s$ bands (McCracken et al. 2012), from Spitzer in the four IRAC bands (Sanders et al. 2007; Ashby et al. 2013), VLA (Schinnerer et al. 2007), XMM (Cappelluti et al. 2007; Hasinger et al. 2007), Chandra (Elvis et al. 2009; Civano et al. 2016), and GALEX (Schiminovich et al. 2005). There are also numerous medium- and narrowband observations available in the COSMOS from Mayall NEWFIRM (Whitaker et al. 2011) and Subaru Suprime-Cam (Taniguchi et al. 2015). The full COSMOS field has been observed with HST in F814W (Koekemoer et al. 2007) and contains more than 2 million galaxies with multi-band data from the UV to the far-IR (Capak et al. 2007; Mobasher et al. 2007; Ilbert et al. 2009, 2013). Furthermore,

<footnote n="31">³¹ http://vizier.cfa.harvard.edu/viz-bin/VizieR?-source=J/ApJS/206/10 and http://vizier.cfa.harvard.edu/viz-bin/VizieR?-source=J/ApJS/207/24</footnote>

<footnote n="32">³² http://candels.ucolick.org/</footnote>

<footnote n="33">³³ https://archive.stsci.edu/prepds/candels/COSMOS.html</footnote>

<footnote n="34">³⁴ https://rainbowx.fis.ucm.es/Rainbow_navigator_public/</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Transmission curves for the CFHT and ACS visible (top left), *Subaru* optical broadband and narrowband (top right), WFC3 and UltraVISTA near-infrared and *Spitzer* infrared (bottom left) and the *Subaru* intermediate- and NEWFIRM medium-band filters (bottom right) used in the COSMOS CANDELS TFIT catalog. These cover observations at ~0.3–8 µm in 42 filters. The filters are adopted from http://cosmos.astro.caltech.edu/page/filterset. The effective wavelength of the filters are reported in Tables 1 and 2.</figcaption></figure>

COSMOS has been followed up spectroscopically by the larger and most sensitive telescopes/instruments like VLT/VIMOS (e.g., Lilly et al. 2009; Le Fèvre et al. 2015) and Keck/DEIMOS and MOSFIRE among others.<fnref n="35" /> This makes it possible to study different populations of galaxies from blue and young systems at shorter wavelengths to red dusty or old objects at longer wavelengths. These ancillary data are accompanied by high-resolution observations from CANDELS using the *HST* in both visible and near-infrared. Figure 1 shows the transmission curves for all bands included in the CANDELS COSMOS catalog.

## 2.1. CANDELS HST Observations

The CANDELS COSMOS field was observed by the WFC3 in F125W and F160W ($J_{125}$ and $H_{160}$) and in parallel by the ACS in the F606W and F814W filters ($V_{606}$ and $i_{814}$). The WFC3 observations covered a rectangular grid of 4 × 11 tiles (~8.6 × 23.8) running north to south, allowing for maximum contiguous coverage in the near-infrared. The field was observed over two epochs with each tile observed for one orbit in each epoch. The one-orbit observations were divided into two exposures in F125W (~1/3 orbit depth) and F160W (~2/3 orbit depth) along with parallel ACS observations in the F606W and F814W (Grogin et al. 2011; Koekemoer et al. 2011). The exposures in each orbit were dithered using a small-scale pattern providing half-pixel subsampling of the PSF and also ensuring that the hot pixels and persistences were moved

around. The output is a calibrated and astrometry-corrected mosaics of all the exposures in the four individual *HST* bands. The astrometry is based on the CFHT/MegaCam $i^*$ imaging supplemented by deep *Subaru*/Suprime-Cam $i^+$ imaging with absolute astrometry registered to the VLA 20 cm survey of the COSMOS field (Schinnerer et al. 2007; see also Koekemoer et al. 2007). This is also the adopted reference grid by the COSMOS team (Capak et al. 2007). All the ground-based and *Spitzer* data (described in the next Section) were aligned to the *HST* data astrometry using SWARP. For the present work we used the V0.5 release of the *HST*/ACS and WFC3 data available from the CANDELS website.<fnref n="36" /> The observation depth, effective wavelength, and the PSF information for each of the *HST* filters are summarized in Table 1.

## 2.2. Ground-based Observations

The full COSMOS field was targeted by the Canada-France-Hawaii Telescope (CFHT) 3.6 m telescope MegaPrime instrument in the $u^*, g^*, r^*, i^*,$ and $z^*$ optical bands as part of the CFHT Legacy Survey<fnref n="37" /> with COSMOS being in the second Deep field. The MegaCam/MegaPrime camera used for the observations has a pixel scale of 0.187 arcsec pixel⁻¹ (Boulade et al. 2003). The final images are from MegaPipe<fnref n="38" /> (Gwyn 2008) with zero-point adjustments. The image processing and stacking

<footnote n="35">³⁵ For a complete list of ancillary data on the entire field visit: http://astro.caltech.edu/~cosmos.</footnote>

<footnote n="36">³⁶ http://candels.ucolick.org/data_access/Latest_Release.html</footnote>

<footnote n="37">³⁷ http://www.cfht.hawaii.edu/Science/CFHTLS/</footnote>

<footnote n="38">³⁸ http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/en/megapipe/cfhtls/index.html</footnote>
<table><caption>**Table 1**
Summary of the CANDELS COSMOS Broadband Data</caption><thead><tr><th>Instrument</th><th>Filter<sup>a</sup></th><th>Effective<br>Wavelength<sup>b</sup><br>(Å)</th><th>PSF<br>FWHM<br>(arcsec)</th><th>5σ limiting depth<sup>c</sup><br>(AB Magnitude)</th><th>Version</th><th>References</th></tr></thead><tbody><tr><td rowspan="5">CFHT/MegaPrime</td><td>u<sup>a</sup></td><td>3817</td><td>0.93</td><td>27.31</td><td rowspan="5">July 2009</td><td rowspan="5">Gwyn (2012)</td></tr><tr><td>g<sup>a</sup></td><td>4860</td><td>1.08</td><td>27.69</td></tr><tr><td>r<sup>a</sup></td><td>6220</td><td>0.84</td><td>27.18</td></tr><tr><td>i<sup>a</sup></td><td>7606</td><td>0.85</td><td>27.23</td></tr><tr><td>z<sup>a</sup></td><td>8816</td><td>0.84</td><td>26.18</td></tr><tr><td rowspan="6">Subaru/Suprime-Cam</td><td>B</td><td>4448</td><td>0.95</td><td>27.98</td><td rowspan="6">V2</td><td rowspan="6">Taniguchi et al. (2007)</td></tr><tr><td>g<sup>+</sup></td><td>4761</td><td>1.58</td><td>26.78</td></tr><tr><td>V</td><td>5470</td><td>1.33</td><td>26.86</td></tr><tr><td>r<sup>+</sup></td><td>6276</td><td>1.05</td><td>27.18</td></tr><tr><td>i<sup>+</sup></td><td>7671</td><td>0.95</td><td>26.95</td></tr><tr><td>z<sup>+</sup></td><td>9096</td><td>1.15</td><td>25.55</td></tr><tr><td rowspan="2">HST/ACS</td><td>F606W</td><td>5919</td><td>0.10</td><td>28.34</td><td rowspan="2">V0.5</td><td rowspan="2">Koekemoer et al. (2011)</td></tr><tr><td>F814W</td><td>8060</td><td>0.10</td><td>27.72</td></tr><tr><td rowspan="2">HST/WFC3</td><td>F125W</td><td>12486</td><td>0.14</td><td>27.72</td><td rowspan="2">...</td><td rowspan="2">...</td></tr><tr><td>F160W</td><td>15369</td><td>0.17</td><td>27.56</td></tr><tr><td rowspan="4">VISTA/VIRCAM</td><td>Y</td><td>10210</td><td>1.17</td><td>25.47</td><td rowspan="4">DR1</td><td rowspan="4">McCracken et al. (2012)</td></tr><tr><td>J</td><td>12524</td><td>1.07</td><td>25.26</td></tr><tr><td>H</td><td>16431</td><td>1.00</td><td>24.87</td></tr><tr><td>K<sub>s</sub></td><td>21521</td><td>0.98</td><td>24.83</td></tr><tr><td rowspan="4">Spitzer/IRAC</td><td>3.6 µm</td><td>35569</td><td>1.80</td><td>24.41</td><td rowspan="4">V1.2</td><td rowspan="4">Ashby et al. (2013)</td></tr><tr><td>4.5 µm</td><td>45020</td><td>1.86</td><td>24.40</td></tr><tr><td>5.8 µm</td><td>57450</td><td>2.13</td><td>21.28</td></tr><tr><td>8.0 µm</td><td>79158</td><td>2.29</td><td>21.20</td></tr></tbody></table>

**Notes.**

<sup>a</sup> Filters adopted from http://cosmos.astro.caltech.edu/page/filterset.

<sup>b</sup> Calculated as $\lambda_{eff} = \sqrt{\int(S(\lambda)\lambda d\lambda) / \int(S(\lambda)\lambda^{-1}d\lambda)}$ with $S(\lambda)$ the filter response function (Tokunaga & Vacca 2005).

<sup>c</sup> The 5σ limiting magnitude calculated within a circular aperture with a radius $r_{ap} = \text{FWHM}$ of the PSF in each filter.

of the data is further described by Gwyn (2012). In this work we used the CFHTLS D2 mosaics 2009.<sup>39</sup>

The COSMOS field was also observed by the Subaru/Suprime-Cam in the B, g+, V, r+, i+ and z+ broadband filters. The Suprime-Cam has a field of view of $34' \times 27'$ with a pixel scale of 0.202 arcsec pixel$^{-1}$. The data were processed using the IMCAT package.<sup>40</sup> The individual frames were combined, flat-fielded, and photometry and astrometry calibrated. We refer the reader to Capak et al. (2007) and Taniguchi et al. (2007, 2015) for further description of the observations and data processing. There is also further observations of the COSMOS field by Subaru/Suprime-Cam in twelve intermediate bands along with narrowband data across two filters (Capak et al. 2007; Taniguchi et al. 2007, 2015) covering the wavelength range of ~4000–8500 Å. These observations were processed similarly to the broadband optical data discussed above (Taniguchi et al. 2015). Although slightly shallower than the optical broadband observations, these filters have higher resolving power than the former (with $R = \lambda/\Delta\lambda \sim 23$; Taniguchi et al. 2015), equivalent to low-resolution spectroscopy in the optical. The resolving power is even higher for the two narrowband filters ($R \sim 50–100$; Taniguchi et al. 2015) although with smaller wavelength coverage. This provides a unique data set for studying emission-line galaxies and high-redshift systems such as Ly$\alpha$ emitters (Shimasaku et al. 2006; Iwata et al. 2009; Koyama et al. 2014). Table 2 summarizes the

Subaru intermediate and narrowband observations. We used the Subaru V2 mosaics for broadband and NB816 observations and V1 mosaics for intermediate and NB711 data available from the IRSA<sup>41</sup> in this work.

The ground-based near-infrared observations are from the VISTA (Emerson & Sutherland 2010) 4.1 m telescope VIRCAM large-format array camera (Dalton et al. 2006) in the Y, J, H, and K<sub>s</sub> bands with mean pixel scale of 0.34 arcsec pixel$^{-1}$. The complete contiguous ~1.5 deg$^2$ of UltraVISTA observations were done using a stripes pattern with ~0.7 deg$^2$ of the field observed with longer exposure (in four stripes) separating the observations into deep and ultra-deep regions (McCracken et al. 2012) with the CANDELS HST observations inside one of the ultra-deep stripes. The data were pre-processed at CASU,<sup>42</sup> which includes dark subtraction, flat fielding, gain normalization, and initial sky subtraction (Irwin et al. 2004; McCracken et al. 2012). The data were further processed at TERAPIX using an iterative sky-background removal technique and resampled to a pixel scale of 0.15 arcsec pixel$^{-1}$. McCracken et al. (2012) give more details on the data processing. We used the final DR1 mosaics of the UltraVISTA data for the CANDELS multi-wavelength catalog.<sup>43</sup>

The COSMOS field also has been observed by the NOAO Extremely Wide-Field Infrared Imager (NEWFIRM) on the Mayall 4 m telescope as part of the NEWFIRM Medium Band

<footnote n="41"><sup>41</sup> http://irsa.ipac.caltech.edu/data/COSMOS/images/subaru/mosaics/</footnote>

<footnote n="42"><sup>42</sup> http://casu.ast.cam.ac.uk/surveys-projects/vista/technical/data-processing</footnote>

<footnote n="43"><sup>43</sup> http://irsa.ipac.caltech.edu/data/COSMOS/images/Ultra-Vista/</footnote>

<footnote n="39"><sup>39</sup> http://www.cadc-ccda.hia-ihna.cnrncrc.gc.ca/en/search/?collection=CFHTMEGAPIPE&noexec=true</footnote>

<footnote n="40"><sup>40</sup> http://www.ifa.hawaii.edu/~kaiser/imcat/</footnote>
<table><caption>**Table 2**
Summary of the CANDELS COSMOS Medium-band and Narrowband Data</caption><thead><tr><th>Instrument</th><th>Filter</th><th>Effective<br>Wavelength<br>(Å)</th><th>PSF<br>FWHM<br>(arcsec)</th><th>5σ Limiting Depth<br>(AB Magnitude)</th><th>Version</th><th>References</th></tr></thead><tbody><tr><td rowspan="14">Subaru/Suprime-Cam</td><td>IA484</td><td>4849</td><td>1.14</td><td>26.34</td><td rowspan="14">V1</td><td rowspan="14">Taniguchi et al. (2007, 2015)</td></tr><tr><td>IA527</td><td>5261</td><td>1.60</td><td>26.04</td></tr><tr><td>IA624</td><td>6232</td><td>1.05</td><td>26.21</td></tr><tr><td>IA679</td><td>6780</td><td>1.58</td><td>25.65</td></tr><tr><td>IA738</td><td>7361</td><td>1.08</td><td>25.94</td></tr><tr><td>IA767</td><td>7684</td><td>1.65</td><td>25.33</td></tr><tr><td>IB427</td><td>4263</td><td>1.64</td><td>25.91</td></tr><tr><td>IB464</td><td>4635</td><td>1.89</td><td>25.66</td></tr><tr><td>IB505</td><td>5062</td><td>1.44</td><td>25.82</td></tr><tr><td>IB574</td><td>5764</td><td>1.71</td><td>25.66</td></tr><tr><td>IB709</td><td>7073</td><td>1.58</td><td>25.79</td></tr><tr><td>IB827</td><td>8244</td><td>1.74</td><td>25.44</td></tr><tr><td>NB711</td><td>7120</td><td>0.79</td><td>25.56</td></tr><tr><td>NB816</td><td>8149</td><td>1.00</td><td>26.10</td></tr><tr><td rowspan="6">Mayall/NEWFIRM</td><td>J<sub>1</sub></td><td>10460</td><td>1.19</td><td>24.60</td><td rowspan="6">DR1</td><td rowspan="6">Whitaker et al. (2011)</td></tr><tr><td>J<sub>2</sub></td><td>11946</td><td>1.17</td><td>24.32</td></tr><tr><td>J<sub>3</sub></td><td>12778</td><td>1.12</td><td>24.26</td></tr><tr><td>H<sub>1</sub></td><td>15601</td><td>1.03</td><td>23.86</td></tr><tr><td>H<sub>2</sub></td><td>17064</td><td>1.24</td><td>23.45</td></tr><tr><td>K</td><td>21700</td><td>1.08</td><td>23.80</td></tr></tbody></table>

Survey<sup>44</sup> (NMBS; Whitaker et al. 2011). The NEWFRIM observations in the COSMOS cover an area of 27.6 × 27.6 encompassing the CANDELS HST observations. The observations are over five medium-band filters of $J_1$, $J_2$, $J_3$, $H_1$, and $H_2$ covering the wavelength of 1–1.8 µm and the K filter centered at 2.2 µm. The three medium-band $J_1$, $J_2$, and $J_3$ filters are a single broadband J filter split into three and the two medium-band $H_1$ and $H_2$ filters combine into a single broadband H (van Dokkum et al. 2009; Whitaker et al. 2011). The final mosaic has a pixel scale of 0.3 arcsec pixel<sup>-1</sup>. Whitaker et al. (2011) further discuss the data processing. Here we used the first data release of the NMBS data (DR1) of the COSMOS field available from the NOAO science archive.<sup>45</sup>

## 2.3. Spitzer Infrared Observations

The COSMOS field was observed by the Spitzer Space Telescope IRAC instrument (Fazio et al. 2004) at 3.6 µm, 4.5 µm, 5.8 µm, and 8.0 µm as part of the S-COSMOS (Sanders et al. 2007). The 3.6 µm and 4.5 µm bands have much deeper data from the Spitzer Extended Deep Survey<sup>46</sup> (SEDS; Ashby et al. 2013). The SEDS observations cover a strip of 10' × 1° oriented north–south coinciding with the deep VISTA data mentioned above and incorporate previous 3.6 µm and 4.5 µm data from the S-COSMOS providing a uniform depth of 26 mag (3σ) for all observations (Ashby et al. 2013). The 5σ limiting magnitude and FWHM size of the PSF in each IRAC band are reported in Table 1. In this work we used the SEDS first data release (V1.2) (Ashby et al. 2013) available from https://www.cfa.harvard.edu/SEDS/data.html for photometry measurements in the 3.6 µm and 4.5 µm bands and the S-COSMOS Spitzer/IRAC 5.8 µm and 8.0 µm GO2 (V2) data (Sanders et al. 2007) available from

http://irsa.ipac.caltech.edu/data/SPITZER/S-COSMOS/. Figure 2 shows the sky coverage of the different ground-based and space data in the COSMOS CANDELS.

## 3. CATALOG PHOTOMETRY

In generating the multi-wavelength catalog, the high-resolution data (HST/ACS and WFC3) were treated differently from the low-resolution data (ground-based and Spitzer/IRAC).

### 3.1. HST Photometry

We performed photometry on the high-resolution (ACS + WFC3) data using SExtractor software version 2.8.6 (Bertin & Arnouts 1996) in dual mode with the WFC3 F160W as the detection band consistently with the multi-wavelength catalogs in the other four CANDELS fields (Galametz et al. 2013; Guo et al. 2013). SExtractor software was modified in several ways to enhance the sky measurement, add a new cleaning procedure, and fix isophotal-corrected magnitude calculations as discussed by Galametz et al. (2013). In order to measure the photometry in the visible bands we PSF-matched the ACS and WFC3 images and extracted the photometry from the matched images in dual mode.

As shown in the photometric studies of the CANDELS UDS and CANDELS GOODS-S fields (Galametz et al. 2013; Guo et al. 2013), it is impossible to identify all galaxies using a single set of parameters (signifying the area, signal to noise, background, etc.) for the extraction. The challenge is to detect the brightest targets while avoiding blending and also detect the faintest objects without introducing spurious sources into the catalog. To this end we used two sets of SExtractor input parameters. One set of parameters is aimed at bright source detection with a focus on deblending extended sources (cold mode), and a second set on faint galaxies (hot mode). The two catalogs generated by the hot and cold parameters were then combined following a routine adopted from GALAPAGOS<sup>47</sup>

<footnote n="44"><sup>44</sup> http://www.astro.yale.edu/nmbs/Overview.html</footnote>

<footnote n="45"><sup>45</sup> http://r2.sdm.noao.edu/nsa/nsa_form.html</footnote>

<footnote n="46"><sup>46</sup> The S-CANDELS data (Ashby et al. 2015) were not used here because they were not available at the time of catalog compilation.</footnote>

<footnote n="47"><sup>47</sup> http://astro-staff.uibk.ac.at/~m.barden/galapagos/</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Sky coverage of data in the COSMOS field. The WFC3 F160W mosaics are shown as the grey shaded region. The entire WFC3 footprint is covered by the CFHT, *Subaru*, and UltraVISTA observations (which are much larger than the scales of this figure).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Left: the $1\sigma$ limiting magnitude distribution (per bin of 0.01) in the WFC3 F160W detection band in each pixel normalized to an area of $1 \text{ arcsec}^2$. Middle: cumulative distribution of the area with a sensitivity greater than a given $1\sigma$ limiting magnitude. Right: distribution of the exposure time.</figcaption></figure>

(Barden et al. 2012). The combined catalog includes all the sources from the cold-mode catalog plus sources in the hot-mode catalog that do not exist in the cold mode as identified by the Kron ellipse of a cold mode detected source as discussed by Galametz et al. (2013). Figure 3 shows the $1\sigma$ detection limit of the combined catalog computed over a circular aperture with a radius of $1 \text{ arcsec}$ along with the cumulative distribution of the detection area and the exposure time distribution in the F160W band. Figure 4 shows the magnitude distribution of the sources in the hot, cold, and combined catalogs along with the comparison of the F160W combined counts with the CANDELS UDS (Galametz et al. 2013) and 3D-HST (Skelton et al. 2014). Table 3 gives the number counts in magnitude bins for the combined catalog along with the associated Poissonian uncertainties.

At this stage we also assigned a photometry flag to every object in the catalog. The flagging system is the same as that adopted by

the other CANDELS fields (Galametz et al. 2013; Guo et al. 2013) and discussed in detail by Galametz et al. (2013). We use a zero for a good photometry in the flagging and assigned a value of one for bright stars and spikes associated with those stars as photometry for objects contaminated by this would be unreliable. A photometric flag of two is associated with the edges of the image as measured from the F160W rms maps.

## 3.2. Ground-based and Spitzer Photometry

In order to measure the photometry in the ground-based and Spitzer bands, we used the Template FITting method (TFIT; Laidler et al. 2007) similarly to the other CANDELS multi-wavelength catalogs (Galametz et al. 2013; Guo et al. 2013, M. Stefanon et al. 2017, in preparation). TFIT is a robust algorithm for measuring photometry in mixed-resolution data sets. Sources that are well separated in the high-resolution image (HST) could be blended in the low-resolution image
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Left: the number of galaxies in the CANDELS COSMOS catalog in bins of F160W magnitude (black filled circles). The number counts of the cold-mode selected galaxies (bright sample) and hot-mode-only selected galaxies (faint sample) are shown in blue and red, respectively. The uncertainties associated with the total counts are Poisson errors. The counts and the associated uncertainties are reported in Table 3. Right: number counts of the combined CANDELS COSMOS catalog compared to the CANDELS UDS (Galametz et al. 2013) and the 3D-HST COSMOS (Skelton et al. 2014).</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Table 3**
HST/WFC3 F160W Number Counts from the Combined Hot+Cold Catalog</figcaption></figure>

<table><thead><tr><th>WFC3 F160W<sup>a</sup></th><th>N(deg<sup>-2</sup> mag<sup>-1</sup>)</th><th>Poisson Uncertainty</th></tr></thead><tbody><tr><td>15.25</td><td>70</td><td>50</td></tr><tr><td>15.75</td><td>387</td><td>117</td></tr><tr><td>16.25</td><td>563</td><td>141</td></tr><tr><td>16.75</td><td>633</td><td>149</td></tr><tr><td>17.25</td><td>1196</td><td>205</td></tr><tr><td>17.75</td><td>1653</td><td>241</td></tr><tr><td>18.25</td><td>2638</td><td>305</td></tr><tr><td>18.75</td><td>3166</td><td>334</td></tr><tr><td>19.25</td><td>4960</td><td>418</td></tr><tr><td>19.75</td><td>6965</td><td>495</td></tr><tr><td>20.25</td><td>9990</td><td>593</td></tr><tr><td>20.75</td><td>13789</td><td>696</td></tr><tr><td>21.25</td><td>17377</td><td>782</td></tr><tr><td>21.75</td><td>27086</td><td>976</td></tr><tr><td>22.25</td><td>33735</td><td>1089</td></tr><tr><td>22.75</td><td>44956</td><td>1257</td></tr><tr><td>23.25</td><td>62720</td><td>1485</td></tr><tr><td>23.75</td><td>82841</td><td>1707</td></tr><tr><td>24.25</td><td>106200</td><td>1933</td></tr><tr><td>24.75</td><td>133426</td><td>2166</td></tr><tr><td>25.25</td><td>156045</td><td>2343</td></tr><tr><td>25.75</td><td>178769</td><td>2508</td></tr><tr><td>26.25</td><td>201810</td><td>2664</td></tr></tbody></table>

**Note.**

(ground-based or *Spitzer*). TFIT uses position and light profiles from the high-resolution image to calculate templates that are used to measure the photometry in the low-resolution image. It does that by smoothing the high-resolution image to match the PSF of the low-resolution image using a convolution kernel (Galametz et al. 2013; Guo et al. 2013). Fluxes in the low-resolution image are then measured using these templates while fitting the sources simultaneously.

TFIT requires some pre-processing of low-resolution image in terms of orientation and pixel scale. The individual steps taken are as follows.

**Background Subtraction:** The low-resolution images must be background subtracted before running TFIT. We used a background subtraction routine with several iterations that included a first estimate through smoothing the image on large scales followed by PSF smoothing and source masking that led to a noise map that was interpolated to determine the background (Galametz et al. 2013).

**Image Scale:** TFIT requires the low-resolution image pixel scale to be an integer multiple of the high-resolution detection image (0.06 arcsec pixel⁻¹ for the F160W) and that both images have the same orientation. Furthermore, the astrometry of the low-resolution images must be consistent with the high-resolution observations. We used SWARP to resample low-resolution data sets to the next larger pixel scale that is an integer multiple of the WFC3 mosaic and also used it for astrometry and image alignment.

**Point-spread Function and Kernel:** The point-spread function of both the high-resolution and low-resolution images are needed for the TFIT pre-processing. We constructed the PSF by stacking isolated and unsaturated stars in each band using custom IDL routines. We also constructed a kernel to convolve the high-resolution templates to the low-resolution ones. The kernel was constructed using a Fourier space analysis technique similar to Galametz et al. (2013), which takes the ratio of the Fourier transform of each PSF. This gives the Fourier transform of the kernel, which is then transformed back into normal space generating the kernel. As discussed by Galametz et al. (2013) a low passband filter is applied in the Fourier space to cancel the high-frequency fluctuations and remove the effect of noise. For *Spitzer*/*IRAC*, which has the largest difference in resolution from *HST*, one could use the PSFs directly as the convolution kernel (Galametz et al. 2013). We generated a model PSF by averaging a set of oversampled PSFs that measure PSF variations across the detector. The final PSF is a boxcar kernel

<footnote><sup>a</sup> Bin center magnitude.</footnote>
smoothed and flux normalized model, which also incorporates all the PAs associated with different Astronomical Observation Requests (AORs). We refer the reader to Galametz et al. (2013) and Guo et al. (2013) for more detail.

**Dilation Correction for High Resolution Images:** TFIT uses the area of the Galaxy identified from the high-resolution HST bands to measure the photometry. These are pixels defined by the high-resolution segmentation maps and are fed into TFIT in the form of the isophotal areas of the high-resolution image. However, as demonstrated by Galametz et al. (2013) and previously by De Santis et al. (2007), SExtractor usually underestimates the isophotal area of faint or small galaxies, and this leads to an underestimate of the flux for such systems. Galametz et al. (2013) performed extensive tests to quantify and correct for this effect, the so-called “dilation correction,” by refining and applying the public DILATE code (De Santis et al. 2007). These simulations showed that the correction factor is negligible for objects with large isophotal areas and that it is largest for objects with area <60 pixels. The original isophotal area size from SExtractor hence defines the dilation factor applied. We used the same criteria to correct our high-resolution SExtractor segmentation maps as outlined by Galametz et al. (2013) and refer the reader to this work for further details. Even after this correction, the total flux measurement is always bound by uncertainties incorporated into the above assumptions and this constitutes one of the limitations of photometry estimation.

We ran TFIT in two stages. During the first step TFIT measured any remaining misalignment in the form of distortion or mis-registration between the high-resolution and low-resolution images in the form of shifted kernels (Laidler et al. 2007; Galametz et al. 2013; Guo et al. 2013). In the second step, TFIT used the kernels measured from the first step to correct the misalignment and construct a difference residual map of the low-resolution bands. Figure 5 shows the TFIT residual maps in the low-resolution visible, near-infrared, and infrared bands along with the high-resolution F160W detection band. The residual maps in the visible and near-infrared are close to zero and only show residuals in the center of very bright objects. However, as argued by Guo et al. (2013), the residual maps are qualitative representations of the TFIT photometry measurements, and we later show in our data quality checks that the photometry is properly measured for these bright objects. Tables 1 and 2 summarize the PSF FWHMs used for the high- and low-resolution images.

The final F160W SExtractor catalog has 38671 sources over the 216 arcmin² area of the CANDELS COSMOS field with WFC3 F160W observations. TFIT keeps the original F160W SExtractor ID and coordinates for each object in the combined hot+cold catalog and therefore we only need to combine corresponding entries from the high-resolution and low-resolution catalogs. Figure 6 shows the 5σ limiting depth for all the filters in the catalog as computed and tabulated in Tables 1 and 2. We further derive and report a weight for each target in the catalog calculated from the F160W rms maps at the SExtractor positions for each object as described in Guo et al. (2013).

# 4. DATA QUALITY CHECKS

We checked the TFIT measured photometry by comparing it with other independently measured photometry in the field.

Additionally we checked the colors of point sources in the catalog against model predictions and color-color plots.

## 4.1. Stars Color Checks

The color of stars changes as a function of their spectral type, which in turn depends on the mass (and hence temperature) among other parameters (Kurucz 1979; Vandenberg 1985; Baraffe et al. 1998). Using this, we could compare the measured color of the point-like objects in our catalog against predictions of the colors of stars derived from stellar physics. The predicted colors of stars were computed from the stellar library of BaSeL (Lejeune et al. 1997; Westera et al. 2002) to measure the model stellar colors. We present the comparisons on color-color plots using several observed filters. To measure the predicted photometry from the templates, we integrated the model stellar SED over the wavelength range of each filter taking into account the filter response functions. This provides the predicted colors of point sources. Figure 7 shows our TFIT measured colors for the point-like objects compared to the colors from the BaSeL stellar library. The color trend of our point-like objects agrees with the general distribution of colors predicted by the stellar models. This further confirms our measured photometry, specifically for the brighter sources, and shows no systematic bias in the photometry. The scatter at the redder color is mostly associated with the fainter sources in the catalog and also due to the intrinsic scatter of colors inherent to the library because of the degeneracies among the different populations of stars.

## 4.2. Infrared Color Validation Check

A validation check of the infrared colors of objects in the catalog was done by using the Spitzer/IRAC TFIT measured photometry of galaxies to identify luminous active galactic nuclei (AGNs). Mid-infrared observations of galaxies have been used extensively to identify and study AGN host galaxies (e.g., Laurent et al. 2000; Farrah et al. 2007; Petric et al. 2011; Yan et al. 2013; Lacy et al. 2015). Several recent studies have used wide and deep observations in the mid-infrared by Spitzer/IRAC to successfully identify large samples of AGNs using flux ratios (Lacy et al. 2004, 2007; Stern et al. 2005, 2012; Donley et al. 2012; Messias et al. 2012). These color selections are based on the power-law behavior of the mid-infrared continuum of luminous AGNs caused by heated dust, producing a thermal continuum (Neugebauer et al. 1979; Ivezić et al. 2002; Donley et al. 2012; Messias et al. 2012).

Figure 8 shows the TFIT-measured IRAC color distributions of galaxies in our catalog that have a S/N > 5 in all four channels. According to these selections the red IRAC colors are expected to be dominated by emission from AGN-heated dust, as also predicted by previous studies of SDSS and radio-selected quasars (Lacy et al. 2004), whereas any stellar components usually shifts the $S_{5.8}/S_{3.6}$ ratio to bluer colors. The Donley et al. (2012) criteria are more conservative in selecting AGNs by removing star-forming and quiescent galaxies identified through other selection methods from optical and near-infrared observations (such as the BzK and LBG selections; Madau et al. 1996; Giavalisco 2002). Figure 8 also shows the IRAC color distribution of the X-ray detected sources in COSMOS (Cappelluti et al. 2007) along with the IRAC color distribution of point sources. The majority of the X-ray detected AGNs have IRAC color distributions consistent
<figure><img src="image_8.png" /><figcaption>**Figure 5.** WFC3 F160W source selection band (top left) along with TFIT residual maps in the optical (top right), near-infrared (bottom left) and IRAC infrared (bottom right). The maps all show the same area as the WFC3 F160W band, are in µJy/pixel units, and are all scaled linearly as shown by the respective color bars. The background noise for the *Subaru V*, *UltraVISTA Y*, and IRAC 3.6 µm are at the levels of 0.001 µJy, 0.004 µJy, and 0.006 µJy, respectively (images are background subtracted as discussed in Section 3). The arrow marked on the WFC3 maps shows a reference object at a flux density of 74 µJy.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** 5σ limiting magnitude of the different observations in the CANDELS COSMOS (Tables 1 and 2) as a function of the wavelength. The symbol sizes are proportional to the filter response function widths.</figcaption></figure>

with the selections by Lacy et al. (2004) and Donley et al. (2012). As discussed by Donley et al. (2012), not all the X-ray luminous and IRAC detected sources are identified by AGN color criteria. In fact, Donley et al. (2012) argue that the X-ray-detected QSOs that fall outside the color criteria seem to be more heavily obscured with lower luminosity AGNs such that the host galaxy contributes more to the optical-near-IR flux. The stellar sources and the general population, however, follow different IRAC color distributions as demonstrated by Lacy et al. (2004) and Donley et al. (2012). We also present a check of stellar colors using IRAC 3.6 µm and 4.5 µm bands in Section 6.2 and Figure 22, where our photometry is found to be consistent with synthetic colors.

### 4.3. Validation Checks with Public Photometry

We also checked the TFIT photometry by comparing it with the public photometry available from the 3D-HST survey<fnref n="48" />

<footnote n="48">⁴⁸ http://3dhst.research.yale.edu/Home.html</footnote>
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Color–color diagrams showing the TFIT color of stars (determined from *SExtractor* as objects with CLASS_STAR > 0.9 and *H*<sub>160</sub> < 22 mag) in CANDELS COSMOS (blue) and model stars from the BaSeL library (black) (Lejeune et al. 1997; Westera et al. 2002). The model colors of stars were computed in each filter by integrating the model SED of stars from the library over the filter transmission curves while the observed colors of the stars are directly from the TFIT catalog with no SED inferred zero-point corrections (Table 5) applied. The filters used are from http://cosmos.astro.caltech.edu/page/filterset.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 8.** IRAC color–color diagram and the AGN selection criteria from Lacy et al. (2004; dashed magenta line) and Donley et al. (2012; solid magenta line). Black circles are galaxies from the catalog with S/N > 5 in all four IRAC bands. Sources with CLASS_STAR > 0.95 and *H*<sub>160</sub> < 21 are shown by yellow circles. X-ray point-like detected sources from *XMM-Newton* wide-field observations of COSMOS (Brusa et al. 2010) are shown by large blue circles.</figcaption></figure>

(Skelton et al. 2014). The 3D-*HST* photometry was measured by *SExtractor* on PSF-matched combined *HST*/WFC3 images in three bands (F125W, F140W, and F160W) as the detection (Skelton et al. 2014). The *HST* images were reduced similarly to CANDELS using the same pixel scale and tangent

point (Skelton et al. 2014). Photometry on the low-resolution bands was performed using the MOPHONGO code (Labbé et al. 2005, 2006; Wuyts et al. 2007; Labbé et al. 2013), which takes into account the variations in the PSF size across different filters and in particular the source confusion problem in the low-resolution images by using a combined PSF-matched WFC3 images as a high-resolution image prior for photometry estimation (Labbé et al. 2005; Skelton et al. 2014). The 3D-*HST* adjusted the AUTO fluxes by an aperture correction derived from growth curves and furthermore performed Galactic extinction corrections. The Galactic extinction correction was measured at the center of each filter and was based on Finkbeiner et al. (1999). These corrections are relatively small (≤0.07; Table 5 in Skelton et al. 2014). We took this into account when comparing our fluxes with that of the 3D-*HST*. All fluxes in the 3D-*HST* catalogs were converted to AB magnitudes using a zero-point of 25 (Skelton et al. 2014).

Figure 9 shows the comparison between TFIT measured photometry and the public photometry from the 3D-*HST*. For the source matching we used a radius equal to the FWHM size of the PSF in the F160W (~0.17 arcsec). We find good agreement between the measured fluxes in our catalog and that of the 3D-*HST*. The offset is generally ~0.1 mag. There is, however, a magnitude-dependent trend when comparing the CANDELS and 3D-*HST* photometry. This is related to the difference in the photometry extraction between the CANDELS (TFIT) and the 3D-*HST* (aperture photometry with fixed
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Photometry comparison between CANDELS and 3D-HST (Skelton et al. 2014). The grayscale density map shows all sources and the magenta shows point sources (identified with SExtractor CLASS_STAR > 0.95 and $H_{160} < 25$ mag). The thick and thin cyan lines show the median of the distribution and the corresponding $1\sigma$ confidence intervals. The number reported in each panel represents the median of the offset for the bright end of the distribution between the CANDELS and the 3D-HST photometry (arbitrarily chosen to be $H_{160} < 22$ mag).</figcaption></figure>

apertures for each band; Tables 4–8 of Skelton et al. 2014). Figure 10 shows the comparison between CANDELS COS-MOS TFIT measured (Band − F160W) color and the corresponding colors of sources measured from the 3D-HST catalog

(Skelton et al. 2014) as a function of the F160W magnitude. For both catalogs the F160W band is taken as the reference band for measuring the color. While comparing the *Spitzer* − F160W color for the different fields as a function
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Color comparison of the CANDELS COSMOS and the 3D-HST (Skelton et al. 2014) as a function of the F160W magnitude. Each panel shows the difference of the (Band-F160W) color between CANDELS and 3D-HST (with F160W being the reference band for calculating the colors in both). The cyan lines show the median and $1\sigma$ variations in the color difference. Point sources in the plot are shown as magenta data points and the median of the offset is reported at the bottom left of each panel.</figcaption></figure>

of photometric redshift, we noticed a deviation of ~0.6 mag between the colors in the COSMOS and GOODS-S fields for the red objects in the *Spitzer* 8.0 µm band. There is an

offset (though smaller) in similar color between the 3D-HST COSMOS and GOODS-S fields. Looking at the variation of this color difference between the CANDELS COSMOS and
<figure><img src="image_14.png" /><figcaption>**Figure 11.** Comparison between photometric redshift and spectroscopic redshift for the CANDELS, COSMOS (Ilbert et al. 2013), and the 3D-HST (Skelton et al. 2014). The CANDELS photometric redshift correspond to the median value measured by the different methods.</figcaption></figure>

3D-HST as a function of the F160W, we notice that most of the difference is associated with objected fainter than ~22 mag, which is similar to our 5σ detection limit.

## 5. PHOTOMETRIC REDSHIFT AND STELLAR MASS ESTIMATES

The methods we used to measure the photometric redshifts and stellar masses are presented in Dahlen et al. (2013) and Mobasher et al. (2015), respectively. Measurements of these parameters are carried out for other CANDELS fields, including CANDELS UDS, CANDELS GOODS-S (Santini et al. 2015), CANDELS EGS (M. Stefanon et al. 2017, in preparation), and CANDELS GOODS-N (G. Barro et al. 2017, in preparation).

### 5.1. Photometric Redshifts

We provided the CANDELS COSMOS photometric catalog to individual teams in the collaboration. The teams were asked to estimate photometric redshifts to galaxies using a calibrating sample containing spectroscopic redshifts. The spectroscopic redshifts were taken from the zCOSMOS compilation (Lilly et al. 2007). Only redshifts for sources in the CANDELS COSMOS area with clear emission-line features were used with uncertain redshifts left out. The spectroscopic redshifts used here are all in public domain. However, for training purposes, a larger sample of unpublished spectroscopic redshifts were used from zCOSMOS (M. Salvato 2016, private communication). Measurements from different teams are in good agreement and also agreed with an independent sample of 448 high-quality spectroscopic redshifts not used to calibrate the photometric redshift methods. A total of six individuals participated. Methods included various fitting codes using minimum χ² and MCMC along with varying Star Formation Histories (SFH). The details of these methods are outlined in Table 7 in the Appendix. Using simulations, we showed that one could determine redshifts to an accuracy of 0.025 in $\Delta z = (z_{phot} - z_{spec})/(1 + z_{spec})$. For each galaxy, we derived the median redshift from different methods and considered that as the redshift estimate for that galaxy. The confidence intervals were measured by combining the intervals from different methods, following the procedure used by Dahlen

<figure><img src="image_15.png" /><figcaption>**Figure 12.** Spectroscopic redshift distributions of the training set used to calibrate the photometric redshifts (in blue) and the independent comparison set (in red).</figcaption></figure>

<table><caption>**Table 4**
Spectroscopic Redshift Comparison Table</caption><thead><tr><th>Survey</th><th>OLF<sup>a</sup></th><th>σ<sub>F</sub><sup>b</sup></th><th>σ<sub>NMAD</sub><sup>c</sup></th><th>σ<sub>O</sub><sup>d</sup></th><th>Number of Galaxies<sup>e</sup></th></tr></thead><tbody><tr><td>CANDELS COSMOS</td><td>0.008</td><td>0.035</td><td>0.011</td><td>0.016</td><td>506</td></tr><tr><td>COSMOS Team</td><td>0.05</td><td>0.071</td><td>0.008</td><td>0.017</td><td>504</td></tr><tr><td>3D-HST</td><td>0.012</td><td>0.045</td><td>0.008</td><td>0.015</td><td>499</td></tr></tbody></table>

**Notes.**

<sup>a</sup> Defined as fraction of objects with $|\Delta z|/(1 + z_{spec}) > 0.15$ where $\Delta z = (z_{phot} - z_{spec})/(1 + z_{spec})$ (Dahlen et al. 2013).

<sup>b</sup> $\sigma_F \equiv \text{rms}(\Delta z/(1 + z_{CANDELS}))$.

<sup>c</sup> $\sigma_{NMAD} \equiv 1.48 \times \text{median}(|\Delta z|/(1 + z_{CANDELS}))$.

<sup>d</sup> $\sigma_O \equiv \text{rms}(\Delta z/(1 + z_{CANDELS}))$ after removing the outliers.

<sup>e</sup> with reliable spectroscopic redshift used for the comparison. This is for objects with $0 < z < 1$ and hence the reported numbers are valid for this redshift range where high quality spectroscopic redshifts are available.

et al. (2013). Taking the median of the photometric redshifts does not necessarily produce a better measurement as this depends heavily on the codes and templates used in calculating the redshift and the corresponding scatter when compared to the spectroscopic redshifts. Comparing with an independent sample of spectroscopic redshifts not used for photometric redshift training, we find that while taking the median reduces the outlier fractions compared to some of the individual codes, marginally lower outlier fractions are obtained using only a subset of the codes (those of Wuyts, Gruetzbach, and Salvato for this field). Given the relatively small number of spectroscopic redshifts available for this test (262), we do not consider the difference significant and choose to report the median of all the codes as our recommended best photometric redshift. Figure 11 shows a direct comparison between the photometric and spectroscopic redshifts for galaxies in the CANDELS COSMOS field, only using high quality spectroscopic redshifts. The training was done using a larger sample of spectroscopic redshifts in CANDELS COSMOS area from the zCOSMOS survey. These do not overlap with the comparison sample here, which is a smaller subsample that is in public domain. Figure 12 shows the spectroscopic redshift distributions of the training and comparison samples. The high-quality spectroscopic redshifts used for the training of the photometric redshifts are available over the range $0 < z < 1$ and
<figure><img src="image_16.png" /><figcaption>**Figure 13.** Photometric redshift comparison plots of CANDELS COSMOS with 3D-HST (left) and COSMOS (right) catalogs. The comparison is shown in the form of a 2D histogram with logarithmic bins to help see the small outlier fraction. The sub-panels in both plots show $\Delta z/(1 + z_{\text{CANDELS}})$ as a function of CANDELS photometric redshifts, where $\Delta z = (z_{\text{CANDELS}} - z_{\text{other}})$ for photometric redshift measurements.</figcaption></figure>

<table><caption>**Table 5**
SED Fitting Measured Photometric Offsets</caption><thead><tr><th>Filter</th><th>Median Offset<sup>a</sup><br>(AB Mag)</th></tr></thead><tbody><tr><td>CFHT-u<sup>*</sup></td><td>0.041</td></tr><tr><td>CFHT-g<sup>*</sup></td><td>-0.044</td></tr><tr><td>CFHT-r<sup>*</sup></td><td>-0.005</td></tr><tr><td>CFHT-i<sup>*</sup></td><td>-0.025</td></tr><tr><td>CFHT-z<sup>*</sup></td><td>0.000</td></tr><tr><td>Subaru-B</td><td>-0.004</td></tr><tr><td>Subaru-r<sup>+</sup></td><td>-0.071</td></tr><tr><td>Subaru-i<sup>+</sup></td><td>-0.074</td></tr><tr><td>Subaru-z<sup>+</sup></td><td>-0.150</td></tr><tr><td>ACS-F606W</td><td>0.072</td></tr><tr><td>ACS-F814W</td><td>0.019</td></tr><tr><td>WFC3 F125W</td><td>0.104</td></tr><tr><td>WFC3 F160W</td><td>0.091</td></tr><tr><td>UVISTA-Y</td><td>-0.027</td></tr><tr><td>UVISTA-J</td><td>0.046</td></tr><tr><td>UVISTA-H</td><td>0.067</td></tr><tr><td>UVISTA-K<sub>s</sub></td><td>-0.031</td></tr><tr><td>IRAC-3.6 µm</td><td>-0.061</td></tr><tr><td>IRAC-4.5 µm</td><td>-0.176</td></tr><tr><td>IRAC-5.8 µm</td><td>-0.069</td></tr><tr><td>IRAC-8.0 µm</td><td>-0.697</td></tr><tr><td>NEWFIRM-J<sub>1</sub></td><td>-0.061</td></tr><tr><td>NEWFIRM-J<sub>2</sub></td><td>-0.027</td></tr><tr><td>NEWFIRM-J<sub>3</sub></td><td>0.000</td></tr><tr><td>NEWFIRM-H<sub>1</sub></td><td>-0.025</td></tr><tr><td>NEWFIRM-H<sub>2</sub></td><td>-0.060</td></tr><tr><td>NEWFIRM-K</td><td>-0.073</td></tr></tbody></table>

**Note.**

hence photometric redshifts reported are most reliable for that redshift range. The values of the redshift difference, as defined in Dahlen et al. (2013), along with the outlier fractions are listed in Table 4. As mentioned above, reliable spectroscopic redshifts are available out to $z \sim 1$ and the numbers listed in Table 4 are only representative for galaxies in this range.

We cross-compared the CANDELS photometric redshift catalog with those from COSMOS (Ilbert et al. 2013) and 3D-HST (Skelton et al. 2014). Figure 11 shows the comparison between the photometric and spectroscopic redshifts in the COSMOS and 3D-HST catalogs. The rms values are listed in Table 4 and show similar trends as CANDELS COSMOS. We compare photometric redshifts between the CANDELS COSMOS and those from the COSMOS (Ilbert et al. 2013) and 3D-HST (Skelton et al. 2014) catalogs, in Figure 13. The two photometric redshift comparison plots present the full scatter ($\sigma_F \equiv \text{rms}(\Delta z/(1 + z_{\text{CANDELS}}))$), the normalized median absolute deviation ($\sigma_{\text{NMAD}} \equiv 1.48 \times \text{median}(|\Delta z|/(1 + z_{\text{CANDELS}}))$) and the scatter after excluding outliers ($\sigma_O \equiv \text{rms}(\Delta z/(1 + z_{\text{CANDELS}}))$) as well as the outlier fraction (defined as fraction of objects with $|\Delta z|/(1 + z_{\text{CANDELS}}) > 0.15$). There is consistency in the photometric redshift measurements for the majority of galaxies. The spread seen in Figure 13, is due to log-binning of the histograms and does not indicate any large inconsistency. Table 5 shows the photometric offsets measured from SED fitting for different bands. The offsets are measured through the SED fitting, simultaneously with photometric redshifts. Different independent SED fitting methods estimated the offsets and they agreed fairly well. The magnitude offsets were then used in the photometry to estimate the final photometric redshifts. The photometric redshifts have the correction for the photometric offsets, but the photometry presented here is not corrected for the offset.

<footnote><sup>a</sup> Positive offset: measured flux fainter than expected from template.</footnote>
<figure><img src="image_17.png" /><figcaption>**Figure 14.** Top: distribution of the photometric redshift difference for the random pair (black) and the close pairs, defined as objects with separation less than 15 arcsec (blue). Bottom: distribution of the overdensity of the photometric redshift difference for close pairs. The distribution is observed in excess of the photometric redshift difference distribution of random galaxies that are subtracted (top panel; Dahlen et al. 2013). A Gaussian fit to the distribution gives the uncertainties associated with the redshifts. This test only characterizes the core of the photometric redshift errors and not the outlier rate.</figcaption></figure>

One method for estimating the photometric redshift uncertainties that is specially useful for fainter flux limits where not enough spectroscopic redshift information is available is the pair statistics estimates (Quadri & Williams 2010; Dahlen et al. 2013; Huang et al. 2013; Hsu et al. 2014). The method, as outlined in Dahlen et al. (2013), relies on the fact that spatially close pairs (defined as objects with separation less than 15 arcsec) have a high probability of being associated with each other and therefore being at similar redshifts (Dahlen et al. 2013). This close-pair association would show up as excess power at small separations when compared to the distribution based on random galaxies (Dahlen et al. 2013). Figure 14 shows the random and close-pair photometric redshift difference distribution along with the excess distribution of the close pairs in photometric redshift after subtracting the photometric redshift distribution of the random galaxies. By fitting a Gaussian function to this excess distribution we measured an uncertainty of 0.017 ($\sqrt{2} \times \sigma_{\text{Gaussian}}$) in the photometric redshift distribution.

## 5.2. Stellar Masses

Mobasher et al. (2015) studied stellar mass measurement using CANDELS data. With extensive simulations, they explored different sources of uncertainty in stellar mass measurements and estimated the error budget associated

with them. The stellar masses were measured from different methods independently and were compared with their expected (input) mass. All methods produced stellar masses in good agreement. We use the median mass (among all the methods) as the reported CANDELS estimate. As we discussed earlier, taking the median does not necessarily produce a more robust mass estimate because the results of individual codes used are not instances of the same stochastic process.

The TFIT multi-waveband photometric catalog for the CANDELS COSMOS field was provided to the CANDELS teams outlined in Table 8 and they measured the stellar masses through separate SED fitting techniques. For all the independent measurements, redshifts were fixed to their median values (as described in the previous section). Given that the photometric redshifts are calibrated with a spectroscopic sample with $0 < z < 1$, the stellar mass estimates are also well calibrated and most robust within this redshift range, as discussed in the previous section, although we report stellar mass estimates out to $z \sim 5$ in this catalog. Some of the methods included nebular emission when fitting the SEDs as outlined in Table 8. A total of eight entries were received from different teams. We measured the median stellar mass between different methods. We used the Hodges-Lehmann<fnref n="49" /> method to estimate the median stellar mass which accounts for the small number of entries when measuring the median value.

Figure 15 presents a direct comparison between different methods used to measure stellar masses. Here, we plot estimates from each method against the median from all the rest of the methods. There is good internal consistency between different methods and the tails seen in Figure 15 help identify where the outliers are in each method compared to the rest. These outliers will not bias our measurements as we report the median of all methods for our final stellar mass measurements as demonstrated in Mobasher et al. (2015).

The median stellar mass estimates from the CANDELS are compared to those from the COSMOS (Ilbert et al. 2013) and 3D-HST catalogs (Skelton et al. 2014) in Figure 16. In these comparisons, we measure the scatter similarly to the photometric redshift where $\sigma_F \equiv \text{rms}(\log(M_{\text{CANDELS}}) - \log(M_{\text{other}}))$ and $\sigma_O \equiv \text{rms}(\log(M_{\text{CANDELS}}) - \log(M_{\text{other}}))$ after removing the outliers where outlier fraction is defined as the fraction of objects with $\Delta \log(M) \equiv (\log(M_{\text{CANDELS}}) - \log(M_{\text{other}})) > 0.5$ (Mobasher et al. 2015), as shown on Figure 16. The CANDELS stellar mass measurements are consistent with both COSMOS and 3D-HST stellar masses. The stellar mass offsets between CANDELS median measurement and 3D-HST and COSMOS are also plotted as a function of F160W magnitude in Figure 17. As expected, the larger discrepancies between the stellar mass estimates occur at fainter magnitudes. The rms values of the stellar mass estimate offsets as defined above are presented in Table 6 as a function of H-band magnitude, confirming a good agreement at brighter magnitudes and an overall consistency for the majority of galaxies.

When measuring the stellar mass, we fit the SEDs by fixing redshifts to the median of the photometric redshifts from the CANDELS COSMOS team, as discussed in Section 5.1. Therefore, if the redshifts for the same galaxies are different in those from 3D-HST and COSMOS teams, the effect would propagate to the estimated stellar masses. As a result, the observed offsets between the stellar mass values in Figures 16 and 17 could partly be explained by the discrepancy between the

<footnote n="49">⁴⁹ http://www.jstor.org/stable/2238406</footnote>
<figure><img src="image_18.png" /><figcaption>**Figure 15.** Comparisons of the reported stellar masses using eight different methods outlined in Table 8 vs. the median stellar mass of all the other methods. The plots are shown as 2D histograms with logarithmic bins. We report the variations in the mass difference and the outliers as defined in Section 5.1 in each panel.</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>**Figure 16.** Stellar mass comparison plots of CANDELS COSMOS with 3D-HST (left) and COSMOS (right) measurements. The comparison is shown in the form of 2D-histogram with logarithmic bins to help see the small outlier fraction. In each panel we report the variations in the mass difference and the outlier fractions as defined in Section 5.1. The sub-panels show $\Delta \log(M)$ as a function of CANDELS $\log(M_*/M_{\odot})$ where $\Delta \log(M) = \log(M_{\text{CANDELS}}) - \log(M_{\text{other}})$ for stellar mass measurement.</figcaption></figure>

redshifts. To explore this, we studied the residual diagrams between redshifts and stellar masses for the three measurements. For this we look at the ratio of the stellar masses as measured by different methods ($\log(M_{\text{CANDELS}}) - \log(M_{\text{other}})$) as a function of the redshift difference $(z_{\text{CANDELS}} - z_{\text{other}})/(1 + z_{\text{CANDELS}})$ in Figure 18 showing a 0.25 dex scatter in stellar mass for

galaxies with similar redshifts. This is in agreement with results of Mobasher et al. (2015), who measured the combined error budget in stellar mass values due to different parameters. Therefore, the distribution in residual mass here is consistent with the expected uncertainties in the stellar mass measurements (i.e., the vertical scatter). The galaxies with deviant redshifts
<figure><img src="image_20.png" /><figcaption>**Figure 17.** Stellar mass offsets between CANDELS and 3D-HST (left) and COSMOS (right) as a function of the F160W magnitude. Both plots are in forms of 2D histograms with logarithmic bins, showing the fraction of inconsistent measurements increase at fainter magnitudes.</figcaption></figure>

<table><caption>**Table 6**
Variations and Outlier Fractions in Mass Measurements</caption><thead><tr><th></th><th>Magnitude Cut (AB)</th><th>OLF<sup>a</sup></th><th>σ<sub>F</sub></th><th>σ<sub>NMAD</sub></th><th>σ<sub>O</sub></th></tr></thead><tbody><tr><td rowspan="4">COSMOS</td><td>21 &lt; H &lt; 22</td><td>0.032</td><td>0.307</td><td>0.124</td><td>0.149</td></tr><tr><td>22 &lt; H &lt; 23</td><td>0.048</td><td>0.309</td><td>0.138</td><td>0.153</td></tr><tr><td>23 &lt; H &lt; 24</td><td>0.077</td><td>0.384</td><td>0.179</td><td>0.175</td></tr><tr><td>24 &lt; H &lt; 25</td><td>0.172</td><td>0.549</td><td>0.251</td><td>0.204</td></tr><tr><td rowspan="4">3D-HST</td><td>21 &lt; H &lt; 22</td><td>0.049</td><td>0.597</td><td>0.133</td><td>0.131</td></tr><tr><td>22 &lt; H &lt; 23</td><td>0.054</td><td>0.463</td><td>0.138</td><td>0.144</td></tr><tr><td>23 &lt; H &lt; 24</td><td>0.089</td><td>0.563</td><td>0.148</td><td>0.148</td></tr><tr><td>24 &lt; H &lt; 25</td><td>0.138</td><td>0.687</td><td>0.201</td><td>0.179</td></tr></tbody></table>

**Note.**

also have deviant stellar mass estimates, partly explaining the observed scatter between the CANDELS and 3D-HST and CANDELS and COSMOS stellar masses.

Figure 19 compares the difference between stellar mass measurements with and without correction for nebular emission lines. This shows a small scatter (0.25 dex) in the stellar mass, consistent with Mobasher et al. (2015), but no significant offset over the whole population over the large redshift range.

### 5.3. Mass Completeness

We use the method introduced by Pozzetti et al. (2010) to estimate the stellar mass completeness limit of the general population of galaxies (see also Ilbert et al. 2013 and Darvish et al. 2015a). Given the magnitude limit of the sample ($H_{lim}$ = 27.56; $5\sigma$ limiting magnitude in the F160W detection band), we assigned a limiting stellar mass ($M_{lim}$) to each galaxy. $M_{lim}$ is the stellar mass that a galaxy would have at its estimated redshift, if its apparent magnitude was the same as the magnitude limit of our sample ($H_{lim}$ = 27.56). This was evaluated by $\log(M_{lim}/M_{\odot}) = \log(M/M_{\odot}) + 0.4(H - H_{lim})$, where M is the estimated stellar mass of the Galaxy with its apparent magnitude H. This results in a distribution of $M_{lim}$ values at any given redshift. The 90% (70%) stellar mass completeness limit at each redshift is therefore equivalent to the mass with 90% (70%) of the galaxies having their $M_{lim}$ value below the stellar mass completeness limit. In general, the stellar mass completeness limit depends on the M/L ratio and is

higher for quiescent and dusty galaxies. The stellar mass estimate is also sensitive to the presence of the Balmer break in the SED of galaxies. At $z \gtrsim 3$ the WFC3 H-band will no longer be probing the Balmer break and a redder filter (such as $K_s$-band or 3.6 $\mu$m) is more suited at estimating the stellar mass at these redshifts (Ilbert et al. 2013). Using the UltraVISTA $K_s$ band limiting magnitude reported in Table 1, we found that while the completeness estimates change at $z \gtrsim 3$, the deviations are at the level of 0.1–0.2 dex, which is within the stellar mass uncertainties. Figure 20 shows the distribution of stellar mass as a function of redshift, along with the estimated 90% and 70% completeness limits.

## 6. APPLICATION TO $2 < Z < 5$ GALAXIES

### 6.1. Star-forming Galaxies at $z \sim 4-5$

One of the main methods of identifying high-redshift galaxies is by targeting the pronounced Lyman break at rest-frame 912 Å that exist in the SED of these galaxies (e.g., Pettini et al. 2002; Shapley et al. 2003; Bolton & Haehnelt 2013; Faucher-Giguère et al. 2015; Williams et al. 2015). The break is caused by the absorption of the UV light from hot and young stars by the neutral Hydrogen (Madau 1995; Madau et al. 1996). Because of this break in the SED, these galaxies would be undetected in the bluer bands and appear in the redder filters (the so-called dropout technique or LBG selection; Madau et al. 1996; Steidel et al. 1999; Giavalisco 2002; Stark et al. 2009). This technique has been used extensively over the past few years in conjunction with deep multi-waveband data and spectroscopic observations to identify and study star-forming galaxies all the way to the cosmic dawn and epoch of reionization (e.g., Yan & Windhorst 2004; Capak et al. 2011; Stark et al. 2011; Steidel et al. 2011; Oesch et al. 2013, 2015; Treu et al. 2013; Bouwens et al. 2014).

Here we use the following LBG color selection from Bouwens et al. (2015) to identify candidates at $\langle z \rangle \sim 4$ and $\langle z \rangle \sim 5$ respectively:

$$ (Subaru(B) - Subaru(V)) > 1, \qquad (1a) $$

$$ (\mathrm{CFHT}(i^*) - \mathrm{WFC3(F125W)}) < 1, \qquad (1b) $$

$$ (Subaru(B) - Subaru(V)) > 1.6 \times (\mathrm{CFHT}(i^*) - \mathrm{WFC3(F125W)}) + 1. \qquad (1c) $$

for the B-dropout and

$$ (Subaru(V) - \mathrm{CFHT}(i^*)) > 1.2, \qquad (2a) $$

<footnote><sup>a</sup> Defined as $|\Delta \log(M)| > 0.5$ (Mobasher et al. 2015).</footnote>
<figure><img src="image_21.png" /><figcaption>Figure 18. Stellar mass offsets vs. photometric redshifts offsets between CANDELS measurements and those from COSMOS (left) and 3D-HST team (right). The solid red lines show the 1:1 relations.</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>Figure 19. Difference of the median of the stellar mass measured with and without nebular emission as a function of the redshift. The dashed red line shows the 1:1 relation. The median and 1σ variations are shown with the light and dark blue respectively. The median is consistent with no evolution as a function of redshift for the two mass estimates.</figcaption></figure>

$$ (\text{CFHT}(z^*) - \text{WFC3(F160W)}) < 1.3, \qquad (2b) $$

$$ (\text{Subaru}(V) - \text{CFHT}(i^*)) > 0.8 \times (\text{CFHT}(z^*) - \text{WFC3(F160W)}) + 1.2 \qquad (2c) $$

for the V-dropout. Figure 21 shows the color–color diagrams for the B-dropout and V-dropout selections and the corresponding photometry used along with the photometric redshift distribution of the selected candidates, as measured in the previous section. The LBG selected sources have photometric redshifts consistent with the selections. There are five spectroscopically confirmed sources (three in the B-dropout and two in the V-dropout) among the selected candidates with colors consistent with the corresponding criteria.

## 6.2. Passively Evolving Galaxies at $z > 1$

The presence of passively evolving galaxies at $z \sim 2$ has been established for some time now (Daddi et al. 2005; Kriek et al. 2006, 2008, 2009; Papovich et al. 2006, 2015; Ilbert et al. 2013; Muzzin et al. 2013). These galaxies are argued to be the progenitors of the most massive systems that form the red sequence at low redshifts (Cassata et al. 2013; Barro et al. 2014; Williams et al. 2014). In the past few years there have

<figure><img src="image_23.png" /><figcaption>Figure 20. Stellar mass as a function of redshift. Points show measured stellar masses for the CANDELS COSMOS sample. Red and blue lines show the 90% and 70% stellar mass completeness limits for the general population of galaxies. The 90% (70%) completeness limit is defined in a sense that only <10% (30%) of the galaxies could be missed in the low-mass end of the Galaxy distribution.</figcaption></figure>

been several studies that predict the presence of these systems at $z \sim 3$ (Nayyeri et al. 2014; Wang et al. 2016) and some predictions as high as $z \sim 5-6$ (Mobasher et al. 2005; Wiklind et al. 2008). Identification of these objects requires deep observations in infrared wavelengths where they have the peak of their emission due to colder and redder stars (Wuyts et al. 2007; Nayyeri et al. 2014). The CANDELS COSMOS deep HST/WFC3 near infrared observations and the Spitzer/IRAC TFIT photometry measurements are crucial for separating red galaxies from blends of fainter sources and lower redshift dusty star-forming galaxies. The multi-band data in the CANDELS COSMOS, especially the wealth of data available on both sides of the Balmer/4000 Å break help constrain the SEDs of old systems and separate them from the red dusty star-forming galaxies.

One well-known color–color diagram is the BzK color plot that is used to identify star-forming and quiescent galaxies at $z \sim 1-2$ (Daddi et al. 2004). The BzK method uses $(z - K)$ versus $(B - z)$ colors of galaxies to separate the two populations at high redshift (Daddi et al. 2004, 2007). This is mostly based on the presence of the Balmer/4000 Å break in
<figure><img src="image_24.png" /><figcaption>**Figure 21.** LBGs selected from CANDELS COSMOS catalog using the selection criteria from Bouwens et al. (2015) as outlined in the text. B-dropout and V-dropout galaxies are shown by blue points in the shaded selection areas with other sources shown by orange points. Sources that fall within the color selection but are not detected do not satisfy the non-detection S/N limit on the blue bands and mostly sit at the boundary of the selection. The subplots in each panel shows the distribution of photometric redshifts of selected galaxies.</figcaption></figure>

<figure><img src="image_25.png" /></figure>

the SED of older galaxies that is observed in the (z − K) color at z ~ 1–2. The color-selection diagram is such that objects with a BzK ≡ (z − K) − (B − z) > −0.2 are identified as actively star-forming galaxies at z ~ 1–2 while objects with a BzK < −0.2 and (z − K) > 2.5 are identified as passively evolving systems at similar redshifts. Variations of the BzK diagram can be used to extend this method to higher redshifts (Daddi et al. 2004; Guo et al. 2012). The so-called VJL and iHM diagrams use the (J − 3.6 µm) versus (V − J) and (H − 4.5 µm) versus (i − H) colors to identify star-forming and quiescent galaxies at z ~ 2.5 and 3.5 respectively. It is based on the same principle as the BzK diagram with the (J − 3.6 µm) and (H − 4.5 µm) colors probing the 4000 Å break at higher redshifts.

We use the *BzK*, VJL, and *iHM* selections (Daddi et al. 2004; Guo et al. 2012) to identify passively evolving systems in the CANDELS COSMOS at *z* > 1. Figure 22 shows the *BzK* color–color plot along with the higher redshift *VJL* and *iHM* plots. The star-forming and quiescent galaxies in each plot are identified in the blue and red regions, respectively. The inset in each plot shows the photometric redshift distribution. As expected, these sources are very red in both the (*B* − *z*) and (*z* − *K*) colors and the corresponding ones in the *VJL* and *iHM* diagrams. This is indicative of the pronounced Balmer/4000 Å break in the SED of these galaxies and the general old stellar population. Figure 22 inset shows the redshift distribution of the quiescent galaxies in red. The *BzK*-identified passive systems have a mean photometric redshift of *z* = 1.61 and a distribution that is consistent with the population being at *z* > 1. Furthermore this population has a mean stellar mass of 10<sup>10.97</sup> *M*<sub>☉</sub>. The photometric redshift and stellar mass distributions and mean values along with the colors of the passive *BzK* galaxies are consistent with the quiescent galaxy selection at high redshift. One of the passive galaxies identified in the CANDELS COSMOS area has been spectroscopically confirmed at *z* = 1.265. This galaxy has a stellar mass of 10<sup>10.98</sup> *M*<sub>☉</sub> and photometric redshift of *z*<sub>phot</sub> = 1.29 consistent with the spectroscopic redshift and with the Galaxy being massive and old at *z* > 1. The *VJL*-selected passive systems have mean photometric redshift of *z* = 2.33 and mean stellar

mass of 10<sup>11.03</sup> M<sub>☉</sub>. The photometric redshift and mass distributions are consistent with passive old galaxy selection at z ~ 2–3. The iHM redshift distribution shows a bi-modality with a large fraction of lower redshift galaxies in the selection. Guo et al. (2012) discussed the 50% contaminant fraction in the iHM color selection. We have 16 spectroscopically confirmed star-forming galaxies in the BzK plot at z ~ 1.5 and two at z ~ 2.5 in the VJL. All these sources have colors consistent with the expected values from the selection criteria. There is only one source in the spectroscopic sample with z<sub>spec</sub> = 1.564 with inconsistent (B − z) and (z − K) colors compared to what is expected for the sBzK sample. This source has z<sub>phot</sub> = 0.682 and is one of the outliers in Figure 11.

The BzK color diagram is well calibrated using large and highly complete spectroscopic samples of galaxies (Daddi et al. 2004) and furthermore it predicts the location of stars on the (z − K) versus (B − z) plane (through the criterion (z − K) < 0.3(B − z) − 0.5; Daddi et al. 2004). Guo et al. (2012) recently extended this analysis to higher redshifts using redder filters and in particular the *Spitzer* 3.6 µm and 4.5 µm observations. This allows us to examine the colors of the stars measured in our catalog against predictions made by the BzK, VJL, and iHM color–color diagrams. Figure 22 shows the BzK colors of point sources from our catalog, identified from the *SExtractor* CLASS_STAR parameter, along with the predicted locus from Daddi et al. (2004). We further show the (J − 3.6 µm) versus (V − J) and (H − 4.5 µm) versus (i − H) colors of stars in our catalog along with the corresponding synthetic colors of stars from BaSeL stellar library (Lejeune et al. 1997; Westera et al. 2002) similar to Guo et al. (2012). We see from Figure 22 that the measured near-infrared and specifically infrared *Spitzer* colors of the stars in our catalog are consistent with the predictions from stellar models.

7. SUMMARY

We used CANDELS HST/WFC3 and ACS observations over a 216 arcmin<sup>2</sup> area of the COSMOS field to construct a multi-wavelength catalog of galaxies that is selected in the WFC3 F160W band. The catalog contains photometry for 38671 sources from ~0.3to8 µm along with physical properties.
<figure><img src="image_26.png" /><figcaption>**Figure 22.** BzK (left), VJL (middle), and iHM (left) color–color plots (Daddi et al. 2004; Guo et al. 2012, 2013) showing the star-forming and passive galaxies positions at $z \sim 1.5$, $z \sim 2.5$ and $z \sim 3.5$, respectively. The star-forming and quiescent population are plotted in blue and red in both diagrams. The corresponding colors of stars in our catalog (identified from SExtractor CLASS_STAR parameter) are shown in each plot with gray symbols. The BzK plot shows the predicted colors of stars as reported by Daddi et al. (2004). We further show the colors of model stars from the BaSeL library (Lejeune et al. 1997; Westera et al. 2002) as magenta diamonds in the VJL and iHM plots. The inset in the plots show the photometric redshift distribution (Section 5) of the star-forming and quiescent populations identified by the corresponding color selections. The redshift distributions of the identified sources are consistent with expectations from the well calibrated color selections further verifying our photometry.</figcaption></figure>

1. We used SExtractor to measure the photometry of objects in the high-resolution (HST/WFC3 and ACS) bands. This also becomes the reference catalog to measure the photometry in the low-resolution (ground-based and Spitzer/IRAC) bands. SExtractor was run using two sets of parameters (the so-called hot and cold modes) that were adjusted to detect both faint/small and bright/blended objects with the combined catalog getting contributions from both.

2. We used TFIT to measure the photometry in the low-resolution bands. This involves using prior information from the high-resolution SExtractor runs (the position and light distribution profiles) to construct templates that were fit to objects in the low-resolution image and from which the photometry was measured.

3. Our combined final catalog contains photometry in all four HST bands observed by CANDELS (two ACS and two WFC3), the optical broadband and narrowband observations by Subaru, optical data by CFHT, near-infrared medium-band and broadband observations by NEWFIRM and UVISTA, respectively, and infrared observations by Spitzer/IRAC in all four channels.

4. We measured the photometric redshift and physical properties of all the objects in our catalog through SED fitting using various codes and reporting the median values with corresponding uncertainties.

5. We selected star-forming galaxies from LBG selections using our measured photometry at $3 < z < 5$ and quiescent systems from BzK and VJL selections at $z > 1$. The photometric redshift distribution of the candidates along with spectroscopic confirmations further verifies the measured photometry.

We wish to thank the anonymous referee for carefully reading the original manuscript and providing very useful suggestions. We also thank S. Fleming for his help with the MAST archive. This work is based on observations taken by the CANDELS Multi-Cycle Treasury Program with the NASA/ESA HST, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. This work is based in part on observations made with the Spitzer Space

Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. Financial support for this work was provided by NSF through AST-1313319 for H.N. and A.C. H.N. further acknowledges support from NASA (grant No. NNX16AF39G). This work is based in part on data products from observations made with ESO Telescopes at the La Silla Paranal Observatories under ESO programme ID 179.A-2005 and on data products produced by TER-APIX and the Cambridge Astronomy survey Unit on behalf of the UltraVISTA consortium. This study was based in part on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii Telescope (CFHT), which is operated by the National Research Council (NRC) of Canada, the Institut National des Sciences de l'Université de la Centre Nationale 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 Canada-France-Hawaii Telescope Legacy Survey, a collaborative project of NRC and CNRS.

## APPENDIX A

### THE HOT/COLD MODE SExTRACTOR PARAMETERS

The cold-mode SExtractor parameters were optimized to detect the bright galaxies while avoiding deblended sources. For this reason we use a Tophat filter with a small deblending threshold and minimum count. The hot-mode was adjusted to detect low surface brightness and small galaxies. Therefore we chose a Gaussian filter which is more suitable for smaller targets and larger deblending parameters. The main SExtractor parameters are listed below with the differences in the cold/hot mode marked in bold.

#—Catalog—

CATALOG_TYPE ASCII_HEAD

#—Extraction—

DETECT_TYPE CCD

FLAG_TYPE OR

DETECT_MINAREA 5.0/10.0
(Continued)

DETECT_THRESHOLD 0.75/0.70

ANALYSIS_THRESHOLD 5.0/0.70

FILTER Y

FILTER_NAME tophat_9.0_9x9 conv/Gauss_4.0_7x7 conv

DEBLEND_NTHRESH 16/64

DEBLEND_MINCONT 0.0001/0.001

CLEAN Y

CLEAN_PARAM 1.0

MASK_TYPE CORRECT

#——Photometry ——

PHOT_FLUSHFRAC 0.2, 0.5, 0.8

PHOT_APERTURES 1.47, 2.08, 2.94, 4.17, 5.88, 8.34,
11.79, 16.66, 23.57, 33.34, 47.13

SATUR_LEVEL 120.0/3900.0

PIXEL_SCALE 0.060

MAG_GAMMA 4.0

#——Star/Galaxy Separation ——

SEEING_FWHM 0.18/0.19

STARNW_NAME default.nnw

#——Background ——

BACK_SIZE 256/128

BACK_FILTERSIZE 9/5

BACKPHOTO_TYPE LOCAL

BACKPHOTO_THICK 100/48

#——Check Image ——

(Continued)

CHECKIMAGE_TYPE SEGMENTATION

#——Memory ——

MEMORY ObJSTACK 4000

MEMORY_PIXSTACK 400000

MEMORY_BUFSIZE 5000

#——Miscellaneous ——

VERBOSE_TYPE NORMAL

#——New Stuff ——

WEIGHT_TYPE MAP_RMS, MAP_RMS

WEIGHT_THRESHOLD 10000.0, 10000.0

PHOT_AUTOPARAMS 2.5, 3.5

GAIN 3070.790

MAG_ZEROPOINT 25.960

## APPENDIX B

### SUMMARY OF THE ASSUMPTIONS USED IN MEASURING THE PHOTOMETRIC REDSHIFT AND STELLAR MASS

This appendix features an outline of the codes and templates used for calculating the photometric redshifts and stellar masses along with the relevant references.

<table><caption>**Table 7**
Summary of the SED Fitting Codes Used in Estimating the Photometric Redshift</caption><thead><tr><td>PI</td><td>Code</td><td>Fitting Method</td><td>Template Set</td><td>Emission Line</td><td>References</td></tr></thead><tbody><tr><td>Finkelstein</td><td>EAZY</td><td>min χ<sup>2</sup></td><td>EAZY+BX418</td><td>Yes</td><td>Brammer et al. (2008), Erb et al. (2010)</td></tr><tr><td>Gruetzbauch</td><td>EAZY</td><td>min χ<sup>2</sup></td><td>EAZY</td><td>Yes</td><td>Brammer et al. (2008)</td></tr><tr><td>Pforr</td><td>HyperZ</td><td>min χ<sup>2</sup></td><td>Maraston05</td><td>No</td><td>Bolzonella et al. (2000), Maraston (2005)</td></tr><tr><td>Salvato</td><td>LePhare</td><td>min χ<sup>2</sup></td><td>BC03+Polletta</td><td>Yes</td><td>Arnouts & Ilbert (2011), Bruzual & Charlot (2003), Polletta et al. (2007)</td></tr><tr><td>Wiklind</td><td>WikZ</td><td>min χ<sup>2</sup></td><td>BC03</td><td>No</td><td>Wiklind et al. (2008)</td></tr><tr><td>Wuyts</td><td>EAZY</td><td>min χ<sup>2</sup></td><td>EAZY</td><td>Yes</td><td>Brammer et al. (2008)</td></tr></tbody></table>

**Note.** See Dahlen et al. (2013) for more detail.

<table><caption>**Table 8**
Summary of the SED Fitting Codes Used in Estimating the Stellar Mass</caption><thead><tr><td>ID</td><td>PI</td><td>Code</td><td>Fitting Method</td><td>Template Set</td><td>Emission Line</td><td>IMF</td><td>References</td></tr></thead><tbody><tr><td>M2</td><td>Barro</td><td>FAST</td><td>min χ<sup>2</sup></td><td>BC03</td><td>No</td><td>Chabrier</td><td>Kriek et al. (2009), Barro et al. (2013)</td></tr><tr><td>M4</td><td>Finkelstein</td><td>own code</td><td>min χ<sup>2</sup></td><td>CB07</td><td>Yes</td><td>Salpeter</td><td>Bruzual (2007), Finkelstein et al. (2012b)</td></tr><tr><td>M6</td><td>Fontana</td><td>zphot</td><td>min χ<sup>2</sup></td><td>BC03</td><td>Yes<sup>a</sup></td><td>Chabrier</td><td>Giallongo et al. (1998), Fontana et al. (2000), Fontana et al. (2006)</td></tr><tr><td>M14</td><td>Lee</td><td>SpeedyMC</td><td>MCMC</td><td>BC03</td><td>Yes</td><td>Chabrier</td><td>Acquaviva et al. (2012)</td></tr><tr><td>M10</td><td>Pforr</td><td>HyperZ</td><td>min χ<sup>2</sup></td><td>M05</td><td>No</td><td>Chabrier</td><td>Bolzonella et al. (2000), Maraston (2005), Maraston et al. (2006), Daddi et al. (2005), Pforr et al. (2012), Pforr et al. (2013)</td></tr><tr><td>M11</td><td>Salvato</td><td>LePhare</td><td>median of PDFs</td><td>BC03</td><td>Yes</td><td>Chabrier</td><td>Arnouts & Ilbert (2011)</td></tr><tr><td>M12</td><td>Wiklind</td><td>WikZ</td><td>min χ<sup>2</sup></td><td>BC03</td><td>No</td><td>Chabrier</td><td>Wiklind et al. (2008)</td></tr><tr><td>M13</td><td>Wuyts</td><td>FAST</td><td>min χ<sup>2</sup></td><td>BC03</td><td>No</td><td>Chabrier</td><td>Kriek et al. (2009), Wuyts et al. (2011)</td></tr></tbody></table>

**Note.** See Mobasher et al. (2015) and Santini et al. (2015) for more detail.

<footnote><sup>a</sup> also without nebular emission included.</footnote>
APPENDIX C
CATALOG ENTRIES

This appendix features CANDELS COSMOS multi-band photometric catalog entries with optical observations from HST/ACS, CFHT and Subaru, near infrared data from HST/WFC3, NEWFIRM and UltraVISTA and infrared observations from Spitzer/IRAC.

<table><caption>**Table 9**
Photometry Catalog Entries</caption><thead><tr><th>Column Number</th><th>Column Name</th></tr></thead><tbody><tr><td># 1 ID</td><td></td></tr><tr><td># 2</td><td>IAU designation</td></tr><tr><td># 3</td><td>R.A.</td></tr><tr><td># 4</td><td>Decl.</td></tr><tr><td># 5</td><td>APCOR</td></tr><tr><td># 6</td><td>CFHT_uS_FLUSH</td></tr><tr><td># 7</td><td>CFHT_uS_FLUSHERR</td></tr><tr><td># 8</td><td>CFHT_gS_FLUSH</td></tr><tr><td># 9</td><td>CFHT_gS_FLUSHERR</td></tr><tr><td># 10</td><td>CFHT_rS_FLUSH</td></tr><tr><td># 11</td><td>CFHT_rS_FLUSHERR</td></tr><tr><td># 12</td><td>CFHT_iS_FLUSH</td></tr><tr><td># 13</td><td>CFHT_iS_FLUSHERR</td></tr><tr><td># 14</td><td>CFHT_zS_FLUSH</td></tr><tr><td># 15</td><td>CFHT_zS_FLUSHERR</td></tr><tr><td># 16</td><td>Subaru_B_FLUSH</td></tr><tr><td># 17</td><td>Subaru_B_FLUSHERR</td></tr><tr><td># 18</td><td>Subaru GP_FLUSH</td></tr><tr><td># 19</td><td>Subaru GP_FLUSHERR</td></tr><tr><td># 20</td><td>Subaru V_FLUSH</td></tr><tr><td># 21</td><td>Subaru V_FLUSHERR</td></tr><tr><td># 22</td><td>Subaru RP_FLUSH</td></tr><tr><td># 23</td><td>Subaru RP_FLUSHERR</td></tr><tr><td># 24</td><td>Subaru IP_FLUSH</td></tr><tr><td># 25</td><td>Subaru IP_FLUSHERR</td></tr><tr><td># 26</td><td>Subaru ZP_FLUSH</td></tr><tr><td># 27</td><td>Subaru ZP_FLUSHERR</td></tr><tr><td># 28</td><td>ACS_F606W_FLUSH</td></tr><tr><td># 29</td><td>ACS_F606W_FLUSHERR</td></tr><tr><td># 30</td><td>ACS_F814W_FLUSH</td></tr><tr><td># 31</td><td>ACS_F814W_FLUSHERR</td></tr><tr><td># 32</td><td>WFC3_F125W_FLUSH</td></tr><tr><td># 33</td><td>WFC3_F125W_FLUSHERR</td></tr><tr><td># 34</td><td>WFC3_F160W_FLUSH</td></tr><tr><td># 35</td><td>WFC3_F160W_FLUSHERR</td></tr><tr><td># 36</td><td>Ultravista_Y_FLUSH</td></tr><tr><td># 37</td><td>UltraVISTA_Y_FLUSHERR</td></tr><tr><td># 38</td><td>UltraVISTA_J_FLUSH</td></tr><tr><td># 39</td><td>UltraVISTA_J_FLUSHERR</td></tr><tr><td># 40</td><td>UltraVISTA_H_FLUSH</td></tr><tr><td># 41</td><td>UltraVISTA_H_FLUSHERR</td></tr><tr><td># 42</td><td>UltraVISTA_Ks_FLUSH</td></tr><tr><td># 43</td><td>UltraVISTA_Ks_FLUSHERR</td></tr><tr><td># 44</td><td>IRAC_Ch1_FLUSH</td></tr><tr><td># 45</td><td>IRAC_Ch1_FLUSHERR</td></tr><tr><td># 46</td><td>IRAC_Ch2_FLUSH</td></tr><tr><td># 47</td><td>IRAC_Ch2_FLUSHERR</td></tr><tr><td># 48</td><td>IRAC_Ch3_FLUSH</td></tr><tr><td># 49</td><td>IRAC_Ch3_FLUSHERR</td></tr></tbody></table>

<table><caption>**Table 9**
(Continued)</caption><thead><tr><th>Column Number</th><th>Column Name</th></tr></thead><tbody><tr><td># 50</td><td>IRAC_Ch4_FLUSH</td></tr><tr><td># 51</td><td>IRAC_Ch4_FLUSHERR</td></tr><tr><td># 52</td><td>NEWFIRM_J1_FLUSH</td></tr><tr><td># 53</td><td>NEWFIRM_J1_FLUSHERR</td></tr><tr><td># 54</td><td>NEWFIRM_J2_FLUSH</td></tr><tr><td># 55</td><td>NEWFIRM_J2_FLUSHERR</td></tr><tr><td># 56</td><td>NEWFIRM_J3_FLUSH</td></tr><tr><td># 57</td><td>NEWFIRM_J3_FLUSHERR</td></tr><tr><td># 58</td><td>NEWFIRM_H1_FLUSH</td></tr><tr><td># 59</td><td>NEWFIRM_H1_FLUSHERR</td></tr><tr><td># 60</td><td>NEWFIRM_H2_FLUSH</td></tr><tr><td># 61</td><td>NEWFIRM_H2_FLUSHERR</td></tr><tr><td># 62</td><td>NEWFIRM_K_FLUSH</td></tr><tr><td># 63</td><td>NEWFIRM_K_FLUSHERR</td></tr><tr><td># 64</td><td>Subaru IB_427_FLUSH</td></tr><tr><td># 65</td><td>Subaru IB_427_FLUSHERR</td></tr><tr><td># 66</td><td>Subaru IB_464_FLUSH</td></tr><tr><td># 67</td><td>Subaru IB_464_FLUSHERR</td></tr><tr><td># 68</td><td>Subaru IB_484_FLUSH</td></tr><tr><td># 69</td><td>Subaru IB_484_FLUSHERR</td></tr><tr><td># 70</td><td>Subaru IB_505_FLUSH</td></tr><tr><td># 71</td><td>Subaru IB_505_FLUSHERR</td></tr><tr><td># 72</td><td>Subaru IA_527_FLUSH</td></tr><tr><td># 73</td><td>Subaru IA_527_FLUSHERR</td></tr><tr><td># 74</td><td>Subaru IB_574_FLUSH</td></tr><tr><td># 75</td><td>Subaru IB_574_FLUSHERR</td></tr><tr><td># 76</td><td>Subaru IA_624_FLUSH</td></tr><tr><td># 77</td><td>Subaru IA_624_FLUSHERR</td></tr><tr><td># 78</td><td>Subaru IA_679_FLUSH</td></tr><tr><td># 79</td><td>Subaru IA_679_FLUSHERR</td></tr><tr><td># 80</td><td>Subaru IB_709_FLUSH</td></tr><tr><td># 81</td><td>Subaru IB_709_FLUSHERR</td></tr><tr><td># 82</td><td>Subaru NB_711_FLUSH</td></tr><tr><td># 83</td><td>Subaru NB_711_FLUSHERR</td></tr><tr><td># 84</td><td>Subaru IA_738_FLUSH</td></tr><tr><td># 85</td><td>Subaru IA_738_FLUSHERR</td></tr><tr><td># 86</td><td>Subaru IA_767_FLUSH</td></tr><tr><td># 87</td><td>Subaru IA_767_FLUSHERR</td></tr><tr><td># 88</td><td>Subaru NB_816_FLUSH</td></tr><tr><td># 89</td><td>Subaru NB_816_FLUSHERR</td></tr><tr><td># 90</td><td>Subaru IB_827_FLUSH</td></tr><tr><td># 91</td><td>Subaru IB_827_FLUSHERR</td></tr><tr><td># 92</td><td>FWHM_IMAGE</td></tr><tr><td># 93</td><td>FLAGS</td></tr><tr><td># 94</td><td>CLASS_STAR</td></tr></tbody></table>

**Notes.** Col. (1): F160W SExtractor ID. Col. (2): IAU designation. Col. (3) & (4): Target coordinates (in degrees). Col. (5): F160W FLUX_AUTO/FLUX_ISO, applied to ACS and WFC3 bands. Col. (6)–(15): CFHT fluxes and errors from TFIT (microJansky). Col. (16)–(27): Subaru broadband fluxes and errors from TFIT (microJansky). Col. (28)–(35): HST/ACS and WFC3 fluxes and errors from TFIT (microJansky). Col. (36)–(43): UltraVISTA near-infrared fluxes and errors from TFIT (microJansky). Col. (44)–(51): Spitzer/IRAC infrared fluxes and errors from TFIT (microJansky). Col. (52)–(63): NEWFIRM medium-band fluxes and errors from TFIT (microJansky). Col. (64)–(91): Subaru medium and narrowband fluxes and errors from TFIT (microJansky). Col. (92): SExtractor F160W FWHM (pixel). Col. (93): Photometry flags. Col. (94): SExtractor stellar classification (1 = Star).

(This table is available in its entirety in machine-readable form.)
<table><caption>**Table 10**
Photometric Redshift Catalog Entries</caption><thead><tr><th>Column Number</th><th>Column Name</th></tr></thead><tbody><tr><td># 1</td><td>ID</td></tr><tr><td># 2</td><td>Spec_z</td></tr><tr><td># 3</td><td>Spec_z_dq</td></tr><tr><td># 4</td><td>Photo_z_Wuyts</td></tr><tr><td># 5</td><td>zinf68_Wuyts</td></tr><tr><td># 6</td><td>zsup68_Wuyts</td></tr><tr><td># 7</td><td>zinf95_Wuyts</td></tr><tr><td># 8</td><td>zsup95_Wuyts</td></tr><tr><td># 9</td><td>Photo_z_Pforr</td></tr><tr><td># 10</td><td>zinf68_Pforr</td></tr><tr><td># 11</td><td>zsup68_Pforr</td></tr><tr><td># 12</td><td>zinf95_Pforr</td></tr><tr><td># 13</td><td>zsup95_Pforr</td></tr><tr><td># 14</td><td>Photo_z_Wiklind</td></tr><tr><td># 15</td><td>zinf68_Wiklind</td></tr><tr><td># 16</td><td>zsup68_Wiklind</td></tr><tr><td># 17</td><td>zinf95_Wiklind</td></tr><tr><td># 18</td><td>zsup95_Wiklind</td></tr><tr><td># 19</td><td>Photo_z_Finkelstein</td></tr><tr><td># 20</td><td>zinf68_Finkelstein</td></tr><tr><td># 21</td><td>zsup68_Finkelstein</td></tr><tr><td># 22</td><td>zinf95_Finkelstein</td></tr><tr><td># 23</td><td>zsup95_Finkelstein</td></tr><tr><td># 24</td><td>Photo_z_Gruetzbauch</td></tr><tr><td># 25</td><td>zinf68_Gruetzbauch</td></tr><tr><td># 26</td><td>zsup68_Gruetzbauch</td></tr><tr><td># 27</td><td>zinf95_Gruetzbauch</td></tr><tr><td># 28</td><td>zsup95_Gruetzbauch</td></tr><tr><td># 29</td><td>Photo_z_Salvato</td></tr><tr><td># 30</td><td>zinf68_Salvato</td></tr><tr><td># 31</td><td>zsup68_Salvato</td></tr><tr><td># 32</td><td>zinf95_Salvato</td></tr><tr><td># 33</td><td>zsup95_Salvato</td></tr></tbody></table>

**Notes.** Col. (1): F160W SExtractor ID. Col. (2): Spectroscopic redshift. Col. (3): Spectroscopic redshift data quality: 1 = secure, 2 = intermediate, 3 = uncertain. Col. (4)–(8): Wuyts photometric redshift estimates and uncertainties. Col. (9)–(13): Pforr photometric redshift estimates and uncertainties. Col. (14)–(18): Wiklind photometric redshift estimates and uncertainties. Col. (19)–(23): Finkelstein photometric redshift estimates and uncertainties. Col. (24)–(28): Gruetzbauch photometric redshift estimates and uncertainties. Col. (29)–(33): Salvato photometric redshift estimates and uncertainties.

(This table is available in its entirety in machine-readable form.)

<table><caption>**Table 11**
Stellar Mass Catalog Entries</caption><thead><tr><th>Column Number</th><th>Column Name</th></tr></thead><tbody><tr><td># 1</td><td>ID</td></tr><tr><td># 2</td><td>Hmag</td></tr><tr><td># 3</td><td>PhotFlag</td></tr><tr><td># 4</td><td>CLASS_STAR</td></tr><tr><td># 5</td><td>AGNFlag</td></tr><tr><td># 6</td><td>zbest</td></tr><tr><td># 7</td><td>zspec</td></tr><tr><td># 8</td><td>q-zspec</td></tr><tr><td># 9</td><td>zphot</td></tr><tr><td># 10</td><td>zphot_l68</td></tr><tr><td># 11</td><td>zphot_u68</td></tr><tr><td># 12</td><td>zphot_l95</td></tr><tr><td># 13</td><td>zphot_u95</td></tr><tr><td># 14</td><td>zphot Ilbert</td></tr><tr><td># 15</td><td>M_med</td></tr></tbody></table>

<table><caption>**Table 11**
(Continued)</caption><thead><tr><th>Column Number</th><th>Column Name</th></tr></thead><tbody><tr><td># 16</td><td>s_med</td></tr><tr><td># 17</td><td>M_neb_med</td></tr><tr><td># 18</td><td>s_neb_med</td></tr><tr><td># 19</td><td>M_14_cons</td></tr><tr><td># 20</td><td>M_11_tau</td></tr><tr><td># 21</td><td>M_6_tau_NEB</td></tr><tr><td># 22</td><td>M_13_tau</td></tr><tr><td># 23</td><td>M_12</td></tr><tr><td># 24</td><td>M_6_tau</td></tr><tr><td># 25</td><td>M_2_tau</td></tr><tr><td># 26</td><td>M_6_deltau</td></tr><tr><td># 27</td><td>M_6_invtau</td></tr><tr><td># 28</td><td>M_10</td></tr><tr><td># 29</td><td>M_4b</td></tr><tr><td># 30</td><td>M_14_lin</td></tr><tr><td># 31</td><td>M_14_deltau</td></tr><tr><td># 32</td><td>M_14_tau</td></tr><tr><td># 33</td><td>M_14_inctau</td></tr><tr><td># 34</td><td>M_14</td></tr><tr><td># 35</td><td>M_neb_med_lin</td></tr><tr><td># 36</td><td>s_neb_med_lin</td></tr><tr><td># 37</td><td>M_med_lin</td></tr><tr><td># 38</td><td>s_med_lin</td></tr></tbody></table>

**Notes.** Col. (1): F160W SExtractor ID. Col. (2): F160W SExtractor MAG_AUTO for convenience (mag). Col. (3): 0 = good nonzero = use with caution or bad. Col. (4): SExtractor stellar classification (1 = Star). Col. (5): AGN. Col. (6): Best of photo-z or spec-z. This is the redshift used for SED fitting. Col. (7): Spectroscopic redshift. Col. (8): Quality of spectroscopic redshift (good = 1). Col. (9): photometric redshift. Col. (10)–(13): 68% and 95% confidence intervals on the photo-z. Col. (14): Photometric redshift from COSMOS catalog. Col. (15): CANDELS reference median stellar mass (dex.solMass)). Col. (16): Standard deviation on M_med (dex.solMass)). Col. (17): Median stellar mass including nebular component (dex.solMass)). Col. (18): Standard deviation on M_neb_med (dex.solMass)). Col. (19): Stellar mass from Method 14_cons (dex.solMass)). Col. (20): Stellar mass from Method 11_tau (dex.solMass)). Col. (21): Stellar mass from Method 6_tau_NEB (dex.solMass)). Col. (22): Stellar mass from Method 13_tau (dex.solMass)). Col. (23): Stellar mass from Method 12 (dex.solMass)). Col. (24): Stellar mass from Method 6_tau (dex.solMass)). Col. (25): Stellar mass from Method 2_tau (dex.solMass)). Col. (26): Stellar mass from Method 6_deltau (dex.solMass)). Col. (27): Stellar mass from Method 6_invtau (dex.solMass)). Col. (28): Stellar mass from Method 10 (dex.solMass)). Col. (29): Stellar mass from Method 4 (dex.solMass)). Col. (30): Stellar mass from Method 14_lin (dex.solMass)). Col. (31): Stellar mass from Method 14_deltau (dex.solMass)). Col. (32): Stellar mass from Method 14_tau (dex.solMass)). Col. (33): Stellar mass from Method 14_inctau (dex.solMass)). Col. (34): Stellar mass best fitted from Method 14 (dex.solMass)). Col. (35)–(36): Median stellar mass including nebular component calculated by the Hodges–Lehmann estimator in the linear space and standard deviation (dex.solMass)). Col. (37)–(38): Median stellar mass with no nebular component calculated by the Hodges–Lehmann estimator in the linear space and standard deviation (dex.solMass)).

(This table is available in its entirety in machine-readable form.)

## REFERENCES

Acquaviva, V., Gawiser, E., & Guaita, L. 2012, in IAU Symp. 284, The Spectral Energy Distribution of Galaxies—SED 2011, ed. R. J. Tuffs & C. C. Popescu (Cambridge: Cambridge Univ. Press), 42

Arnouts, S., & Ilbert, O. 2011, LePHARE: Photometric Analysis for Redshift Estimate, Astrophysics Source Code Library, ascl:1108.009

Ashby, M. L. N., Willner, S. P., Fazio, G. G., et al. 2013, ApJ, 769, 80

Ashby, M. L. N., Willner, S. P., Fazio, G. G., et al. 2015, ApJS, 218, 33

Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, A&A, 337, 403

Barden, M., Häußler, B., Peng, C. Y., McIntosh, D. H., & Guo, Y. 2012, MNRAS, 422, 449

Barro, G., Faber, S. M., Pérez-González, P. G., et al. 2013, ApJ, 765, 104
Barro, G., Faber, S. M., Pérez-González, P. G., et al. 2014, ApJ, 791, 52

Barro, G., Pérez-González, P. G., Gallego, J., et al. 2011, ApJS, 193, 13

Bell, E. F., van der Wel, A., Papovich, C., et al. 2012, ApJ, 753, 167

Bell, E. F., Wolf, C., Meisenheimer, K., et al. 2004, ApJ, 608, 752

Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393

Bolton, J. S., & Haehnelt, M. G. 2013, MNRAS, 429, 1695

Bolzonella, M., Miralles, J.-M., & Pelló, R. 2000, A&A, 363, 476

Boulade, O., Charlot, X., Abbon, P., et al. 2003, Proc. SPIE, 4841, 72

Bouwens, R. J., Bradley, L., Zitrin, A., et al. 2014, ApJ, 795, 126

Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2015, ApJ, 803, 34

Brammer, G. B., van Dokkum, P. G., & Coppi, P. 2008, ApJ, 686, 1503

Bruce, V. A., Dunlop, J. S., Cirasuolo, C., et al. 2013, MNRAS, 427, 1666

Brusa, M., Civano, F., Comastri, A., et al. 2010, ApJ, 716, 348

Bruzual, G. 2007, in ASP Conf. Ser. 374, From Stars to Galaxies: Building the Pieces to Build Up the Universe, ed. A. Vallenari et al. (San Francisco, CA: ASP), 303

Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000

Capak, P., Aussel, H., Ajiki, M., et al. 2007, ApJS, 172, 99

Capak, P., Mobasher, B., Scoville, N. Z., et al. 2011, ApJ, 730, 68

Cappelluti, N., Hasinger, G., Brusa, M., et al. 2007, ApJS, 172, 341

Caputi, K. I., Ilbert, O., Laigle, C., et al. 2015, ApJ, 810, 73

Cassata, P., Giavalisco, M., Williams, C. C., et al. 2013, ApJ, 775, 106

Civano, F., Marchesi, S., Comastri, A., et al. 2016, ApJ, 819, 62

Daddi, E., Cimatti, A., Renzini, A., et al. 2004, ApJ, 617, 746

Daddi, E., Dickinson, M., Morrison, G., et al. 2007, ApJ, 670, 156

Daddi, E., Renzini, A., Pirzkal, N., et al. 2005, ApJ, 626, 680

Dahlen, T., Mobasher, B., Faber, S. M., et al. 2013, ApJ, 775, 93

Dalton, G. B., Caldwell, M., Ward, A. K., et al. 2006, Proc. SPIE, 6269, 0

Darvish, B., Mobasher, B., Sobral, D., et al. 2015a, ApJ, 814, 84

Darvish, B., Mobasher, B., Sobral, D., Scoville, N., & Aragon-Calvo, M. 2015b, ApJ, 805, 121

Darvish, B., Sobral, D., Mobasher, B., et al. 2014, ApJ, 796, 51

De Santis, C., Grazian, A., Fontana, A., & Santini, P. 2007, NewA, 12, 271

Donley, J. L., Koekemoer, A. M., Brusa, M., et al. 2012, ApJ, 748, 142

Duncan, K., Conselice, C. J., Mortlock, A., et al. 2014, MNRAS, 444, 2960

Elvis, M., Civano, F., Vignali, C., et al. 2009, ApJS, 184, 158

Emerson, J., & Sutherland, W. 2010, Msngr, 139, 2

Erb, D. K., Pettini, M., Shapley, A. E., et al. 2010, ApJ, 719, 1168

Faber, S. M., Willmer, C. N. A., Wolf, C., et al. 2007, ApJ, 665, 265

Farrah, D., Bernard-Salas, J., Spoon, H. W. W., et al. 2007, ApJ, 667, 149

Faucher-Giguère, C.-A., Hopkins, P. F., Kereš, D., et al. 2015, MNRAS, 449, 987

Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10

Finkbeiner, D. P., Davis, M., & Schlegel, D. J. 1999, ApJ, 524, 867

Finkelstein, S. L., Papovich, C., Ryan, R. E., et al. 2012a, ApJ, 758, 93

Finkelstein, S. L., Papovich, C., Salmon, B., et al. 2012b, ApJ, 756, 164

Finkelstein, S. L., Ryan, R. E., Jr., Papovich, C., et al. 2015, ApJ, 810, 71

Fontana, A., D'Odorico, S., Poli, F., et al. 2000, AJ, 120, 2206

Fontana, A., Salimbeni, S., Grazian, A., et al. 2006, A&A, 459, 745

Galametz, A., Grazian, A., Fontana, A., et al. 2013, ApJS, 206, 10

Giallongo, E., D'Odorico, S., Fontana, A., et al. 1998, AJ, 115, 2169

Giallongo, E., Grazian, A., Fiore, F., et al. 2015, A&A, 578, A83

Giavalisco, M. 2002, ARA&A, 40, 579

Giavalisco, M., Ferguson, H. C., Koekemoer, A. M., et al. 2004, ApJL, 600, L93

Grazian, A., Castellano, M., Fontana, A., et al. 2012, A&A, 547, A51

Grogin, N. A., Kocevski, D. D., Faber, S. M., et al. 2011, ApJS, 197, 35

Guo, Y., Ferguson, H. C., Giavalisco, M., et al. 2013, ApJS, 207, 24

Guo, Y., Giavalisco, M., Cassata, P., et al. 2012, ApJ, 749, 149

Gwyn, S. D. J. 2008, PASP, 120, 212

Gwyn, S. D. J. 2012, AJ, 143, 38

Hasinger, G., Cappelluti, N., Brunner, H., et al. 2007, ApJS, 172, 29

Hemmati, S., Miller, S. H., Mobasher, B., et al. 2014, ApJ, 797, 108

Hemmati, S., Mobasher, B., Darvish, B., et al. 2015, ApJ, 814, 46

Hemmati, S., Yan, L., Diaz-Santos, T., et al. 2017, ApJ, 834, 36

Hsu, L.-T., Salvato, M., Nandra, K., et al. 2014, ApJ, 796, 60

Huang, J.-S., Faber, S. M., Willmer, C. N. A., et al. 2013, ApJ, 766, 21

Ilbert, O., Capak, P., Salvato, M., et al. 2009, ApJ, 690, 1236

Ilbert, O., McCracken, H. J., Le Fèvre, O., et al. 2013, A&A, 556, A55

Irwin, M. J., Lewis, J., Hodgkin, S., et al. 2004, Proc. SPIE, 5493, 411

Ivezić, Ž., Menou, K., Knapp, G. R., et al. 2002, AJ, 124, 2364

Iwata, I., Inoue, A. K., Matsuda, Y., et al. 2009, ApJ, 692, 1287

Jones, D. O., Rodney, S. A., Riess, A. G., et al. 2013, ApJ, 768, 166

Khostovan, A. A., Sobral, D., Mobasher, B., et al. 2015, MNRAS, 452, 3948

Kocevski, D. D., Faber, S. M., Mozena, M., et al. 2012, AJ, 744, 148

Koekemoer, A. M., Aussel, H., Calzetti, D., et al. 2007, ApJS, 172, 196

Koekemoer, A. M., Faber, S. M., Ferguson, H. C., et al. 2011, ApJS, 197, 36

Kovač, K., Lilly, S. J., Cucciati, O., et al. 2010, ApJ, 708, 505

Koyama, Y., Kodama, T., Tadaki, K.-i., et al. 2014, ApJ, 789, 18

Kriek, M., van der Wel, A., van Dokkum, P. G., Franx, M., & Illingworth, G. D. 2008, ApJ, 682, 896

Kriek, M., van Dokkum, P. G., Franx, M., et al. 2006, ApJL, 649, L71

Kriek, M., van Dokkum, P. G., Labbé, I., et al. 2009, ApJ, 700, 221

Kurucz, R. L. 1979, ApJS, 40, 1

Labbé, I., Bouwens, R., Illingworth, G. D., & Franx, M. 2006, ApJL, 649, L67

Labbé, I., Huang, J., Franx, M., et al. 2005, ApJL, 624, L81

Labbé, I., Oesch, P. A., Bouwens, R. J., et al. 2013, ApJL, 777, L19

Lacy, M., Petric, A. O., Sajina, A., et al. 2007, AJ, 133, 186

Lacy, M., Ridgway, S. E., Sajina, A., et al. 2015, ApJ, 802, 102

Lacy, M., Storrie-Lombardi, L. J., Sajina, A., et al. 2004, ApJS, 154, 166

Laidler, V. G., Papovich, C., Grogin, N. A., et al. 2007, PASP, 119, 1325

Laurent, O., Mirabel, I. F., Charmandaris, V., et al. 2000, A&A, 359, 887

Lawrence, A., Warren, S. J., Almaini, O., et al. 2007, MNRAS, 379, 1599

Le Fèvre, O., Tasca, L. A. M., Cassata, P., et al. 2015, A&A, 576, A79

Lee, B., Giavalisco, M., Williams, C. C., et al. 2013, ApJ, 774, 47

Lejeune, T., Cuisinier, F., & Buser, R. 1997, A&AS, 125, 229

Lilly, S. J., Le Brun, V., Maier, C., et al. 2009, ApJS, 184, 218

Lilly, S. J., Le Fèvre, O., Renzini, A., et al. 2007, ApJS, 172, 70

Lorenzoni, S., Bunker, A. J., Wilkins, S. M., et al. 2013, MNRAS, 429, 150

Madau, P. 1995, ApJS, 441, 18

Madau, P., Ferguson, H. C., Dickinson, M. E., et al. 1996, MNRAS, 283, 1388

Maraston, C. 2005, MNRAS, 362, 799

Maraston, C., Daddi, E., Renzini, A., et al. 2006, ApJ, 652, 85

McCracken, H. J., Milvang-Jensen, B., Dunlop, J., et al. 2012, A&A, 544, A156

McLeod, D. J., McLure, R. J., Dunlop, J. S., et al. 2015, MNRAS, 450, 3032

Messias, H., Afonso, J., Salvato, M., Mobasher, B., & Hopkins, A. M. 2012, ApJ, 754, 120

Mitchell-Wynne, K., Cooray, A., Gong, Y., et al. 2015, NatCo, 6, 7945

Mobasher, B., Capak, P., Scoville, N. Z., et al. 2007, ApJS, 172, 117

Mobasher, B., Dahlen, T., Ferguson, H. C., et al. 2015, ApJ, 808, 101

Mobasher, B., Dickinson, M., Ferguson, H. C., et al. 2005, ApJ, 635, 832

Muzzin, A., Marchesini, D., Stefanon, M., et al. 2013, ApJ, 777, 18

Nayyeri, H., Mobasher, B., Hemmati, S., et al. 2014, ApJ, 794, 68

Neugebauer, G., Oke, J. B., Becklin, E. E., & Matthews, K. 1979, ApJ, 230, 79

Oesch, P. A., Bouwens, R. J., Illingworth, G. D., et al. 2013, ApJ, 773, 75

Oesch, P. A., Brammer, G., van Dokkum, P. G., et al. 2016, ApJ, 819, 129

Oesch, P. A., van Dokkum, P. G., Illingworth, G. D., et al. 2015, ApJL, 804, L30

Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713

Papovich, C., Labbé, I., Quadri, R., et al. 2015, ApJ, 803, 26

Papovich, C., Moustakas, L. A., Dickinson, M., et al. 2006, ApJ, 640, 92

Pérez-González, P. G., Rieke, G. H., Villar, V., et al. 2008, ApJ, 675, 234

Petric, A. O., Armus, L., Howell, J., et al. 2011, ApJ, 730, 28

Pettini, M., Rix, S. A., Steidel, C. C., et al. 2002, ApJ, 569, 742

Pforr, J., Maraston, C., & Tonini, C. 2012, MNRAS, 422, 3285

Pforr, J., Maraston, C., & Tonini, C. 2013, MNRAS, 435, 1389

Pолletta, M., Tajer, M., Maraschi, L., et al. 2007, ApJ, 663, 81

Pozzetti, L., Bolzonella, M., Zucca, E., et al. 2010, A&A, 523, A13

Quadri, R. F., & Williams, R. J. 2010, ApJ, 725, 94

Roberts-Borsani, G. W., Bouwens, R. J., Oesch, P. A., et al. 2016, ApJ, 823, 143

Rodney, S. A., Riess, A. G., Scolnic, D. M., et al. 2016, AJ, 151, 47

Rodney, S. A., Riess, A. G., Strolger, L.-G., et al. 2014, AJ, 148, 13

Sanders, D. B., Salvato, M., Aussel, H., et al. 2007, ApJS, 172, 86

Santini, P., Ferguson, H. C., Fontana, A., et al. 2011, ApJ, 801, 97

Schiminovich, D., Ilbert, O., Arnouts, S., et al. 2005, ApJL, 619, L47

Schinnerer, E., Smolčić, V., Carilli, C. L., et al. 2007, ApJS, 172, 46

Scoville, N., Arnouts, S., Aussel, H., et al. 2013, ApJS, 206, 3

Scoville, N., Aussel, H., Benson, A., et al. 2007a, ApJS, 172, 150

Scoville, N., Aussel, H., Brusa, M., et al. 2007b, ApJS, 172, 1

Shapley, A. E., Steidel, C. C., Pettini, M., & Adelberger, K. L. 2003, ApJ, 588, 65

Shimasaku, K., Kashikawa, N., Doi, M., et al. 2006, PASJ, 58, 313

Shivaei, I., Reddy, N. A., Steidel, C. C., & Shapley, A. E. 2015, ApJ, 804, 149

Skelton, R. E., Whitaker, K. E., Momcheva, I. G., et al. 2014, arXiv:1403.3689

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

Song, M., Finkelstein, S. L., Ashby, M. L. N., et al. 2016, ApJS, 825, 5

Stark, D. P., Ellis, R. S., Bunker, A., et al. 2009, ApJ, 697, 1493

Stark, D. P., Ellis, R. S., & Ouchi, M. 2011, ApJL, 728, L2

Steidel, C. C., Adelberger, K. L., Giavalisco, M., Dickinson, M., & Pettini, M. 1999, ApJ, 519, 1
Steidel, C. C., Bogosavljević, M., Shapley, A. E., et al. 2011, ApJ, 736, 160

Stern, D., Assef, R. J., Benford, D. J., et al. 2012, ApJ, 753, 30

Stern, D., Eisenhardt, P., Gorjian, V., et al. 2005, ApJ, 631, 163

Strolger, L.-G., Dahlen, T., Rodney, S. A., et al. 2015, ApJ, 813, 93

Taniguchi, Y., Kajisawa, M., Kobayashi, M. A. R., et al. 2015, PASJ, 67, 104

Taniguchi, Y., Scoville, N., Murayama, T., et al. 2007, ApJS, 172, 9

Teplitz, H. I., Rafelski, M., Kurczynski, P., et al. 2013, AJ, 146, 159

Tokunaga, A. T., & Vacca, W. D. 2005, PASP, 117, 421

Treu, T., Schmidt, K. B., Trenti, M., Bradley, L. D., & Stiavelli, M. 2013, ApJL, 775, L29

van Dokkum, P. G., Labbé, I., Marchesini, D., et al. 2009, PASP, 121, 2

Vandenberg, D. A. 1985, ApJS, 58, 711

Vasei, K., Siana, B., Shapley, A. E., et al. 2016, ApJ, 831, 38

Wang, T., Elbaz, D., Schreiber, C., et al. 2016, ApJ, 816, 84

Westera, P., Lejeune, T., Buser, R., Cuisinier, F., & Bruzual, G. 2002, A&A, 381, 524

Whitaker, K. E., Franx, M., Leja, J., et al. 2014, ApJ, 795, 104

Whitaker, K. E., Labbé, I., van Dokkum, P. G., et al. 2011, ApJ, 735, 86

Wiklind, T., Dickinson, M., Ferguson, H. C., et al. 2008, ApJ, 676, 781

Williams, C. C., Giavalisco, M., Cassata, P., et al. 2014, ApJ, 780, 1

Williams, C. C., Giavalisco, M., Lee, B., et al. 2015, ApJ, 800, 21

Wolf, C., Meisenheimer, K., Rix, H.-W., et al. 2003, A&A, 401, 73

Wuyts, S., Förster Schreiber, N. M., Nelson, E. J., et al. 2013, ApJ, 779, 135

Wuyts, S., Förster Schreiber, N. M., van der Wel, A., et al. 2011, ApJ, 742, 96

Wuyts, S., Labbé, I., Franx, M., et al. 2007, ApJ, 655, 51

Yan, H., Finkelstein, S. L., Huang, K.-H., et al. 2012, ApJ, 761, 177

Yan, H., & Windhorst, R. A. 2004, ApJL, 612, L93

Yan, L., Donoso, E., Tsai, C.-W., et al. 2013, AJ, 145, 55

York, D. G., Adelman, J., Anderson, J. E., Jr., et al. 2000, AJ, 120, 1579