Combining the CLAUDS&HSC-SSP surveys: U+grizy(+YJHKs) photometry and photometric redshifts for 18M galaxies in the 20 deg2 of the HSC-SSP Deep and ultraDeep fields

We present the combination of the Canada-France-Hawaii Telescope (CHFT) Large Area $U$-bands Deep Survey (CLAUDS) and the Hyper-Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP) data over their four deep fields. We provide photometric catalogs for $u$, $u^*$ (CFHT--MegaCam), $g$, $r$, $i$, $z$, and $y$ (Subaru--HSC) bands over $\sim 20~{\rm deg}^2$, complemented in two fields by data from the Visible and Infrared Survey Telescope for Astronomy (VISTA) Deep Extragalactic Observations (VIDEO) survey and the UltraVISTA survey, thus extending the wavelength coverage toward near-infrared with VIRCAM $Y$, $J$, $H$, and $K_s$ observations over $5.5~{\rm deg}^2$. The extraction of the photometry was performed with two different softwares: the HSC pipeline hscPipe and the standard and robust SExtractor software. Photometric redshifts were computed with template-fitting methods using the new Phosphoros code for the hscPipe photometry and the well-known Le Phare code for the SExtractor photometry. The products of these methods were compared with each other in detail. We assessed their quality using the large spectroscopic sample available in those regions, together with photometry and photometric redshifts from COSMOS2020, the latest version of the Cosmic Evolution Survey catalogs. We find that both photometric data sets are in good agreement in $Ugrizy$ down to magnitude$\sim26$, and to magnitude$\sim24.5$ in the $YJHK_s$ bands. We achieve good performance for the photometric redshifts, reaching precisions of $\sigma_{NMAD} \lesssim 0.04$ down to ${m}_i\sim25$, even using only the CLAUDS and HSC bands. At the same magnitude limit, we measured an outlier fraction of $\eta \lesssim 10\%$ when using the $Ugrizy$ bands, and down to $\eta \lesssim 6\%$ when considering near-infrared data. [abridged]


Introduction
Deep, wide-area multiband imaging surveys play a pivotal role in our studies of the Universe and of the formation and evolution of its content (e.g., CANDELS, Grogin et al. 2011;COSMOS, Scoville et al. 2007; CFHTLS, Hudelot et al. 2012;DES, A&A 670, A82 (2023) modern, large-area mosaic imagers, and is also significantly less affected by selection biases.
Much can be learned from photometry, particularly when samples are large and contain many faint objects.In particular, photometry can be used to estimate redshifts of vast samples of distant galaxies (e.g., Beck et al. 2016;Tanaka et al. 2018;Moutard et al. 2020), to characterize key physical properties such as their stellar masses and star formation rates (SFRs; e.g., Sawicki & Yee 1998;Papovich et al. 2001;Laigle et al. 2016), and even constrain star formation histories (e.g., Pacifici et al. 2016;Iyer et al. 2019).These measurements are key ingredients in studies of how galaxies form and evolve over cosmic time and motivate many of the modern photometric surveys.
The Hyper Suprime-Cam Strategic Survey Program (HSC-SSP, Aihara et al. 2018a), which recently completed observations on the Subaru telescope, represents the current state of the art in deep, wide-area imaging surveys.Of particular interest here are HSC-SSP's Deep and ultraDeep layers, which cover 26 and 3.5 deg 2 , respectively, in the grizy filters to depths of ∼25-28 AB (signal to noise ratio S /N = 5σ in 2 ′′ apertures).This is an unprecedented combination of area and depth that allows studies of galaxy evolution as a function of redshift and environment that are largely insensitive to cosmic variance and that will be unsurpassed until several years into LSST's science operations.
However, the HSC-SSP grizy photometry spans only the optical and lacks both the shorter and the longer wavelengths that are critical for a number of science applications.In particular, complementing grizy imaging with U-band photometry provides better measurements of SFRs at all redshifts and significantly improves photometric redshift performance at z < 0.7 and at z ∼ 2-3 because it allows the combined filter set to span the Balmer and Lyman breaks, respectively, at these redshifts (see Fig. 9 and Sect.4.2 of Sawicki et al. 2019 for more discussion on the U-band benefits for photo-zs).The SFR and photometric redshift measurements in the HSC-SSP Deep and ultraDeep layers were indeed key motivators for the Canada-France-Hawaii Telescope (CFHT) Large Area U-band Deep Survey (CLAUDS; Sawicki et al. 2019), a major (68 dedicated nights plus archival data) imaging survey that complements, with U-band, the grizy imaging in the Deep and ultraDeep layers of HSC-SSP with an overlap of ∼20 deg 2 between both surveys.Meanwhile, at longer wavelengths, complementing the CLAUDS+HSC-SSP U+grizy photometry with near-infrared data would make measurements of galaxy stellar masses beyond z ∼ 1 possible and further improve photometric redshift performance by adding coverage of the Balmer break at higher redshifts.
These ingredients -uniform photometry, photometric redshifts, and galaxy physical parameters -are key elements needed for many studies of distant galaxy populations and their evolution.For this reason, in this paper, we present our latest multiband catalogs that combine CLAUDS, HSC-SSP, andwhere available -VIDEO (VISTA Deep Extragalactic Observations;Jarvis et al. 2013) and UltraVISTA (McCracken et al. 2012) data in the HSC-SSP Deep and ultraDeep fields.In previous papers (Sawicki et al. 2019;Moutard et al. 2020), we presented CLAUDS+HSC-SSP catalogs for these fields based on shallower HSC-SSP photometry and less sophisticated technical treatment.These earlier catalogs have already been used for several science applications: for example, Halevi et al. (2019) used this earlier CLAUDS+HSC-SSP photometry to characterize the properties of an active galaxy nucleus (AGN) containing lowmass galaxy; Moutard et al. (2020) and Harikane et al. (2021) constrained ultra-violet luminosity functions; and Iwata et al. (2022) used them to constrain the Lyman continuum escape fraction from AGNs.Additionally, the data were used for star-galaxy classification (Golob et al. 2021), selection of lowluminosity galaxies (Cheng et al. 2021), photometric redshift estimation (Sawicki et al. 2019;Moutard et al. 2020;Huang et al. 2020), and constraining galaxy-galaxy merger fraction evolution (Thibert et al. 2021).
These earlier CLAUDS+HSC-SSP catalogs already demonstrated the usefulness of a deep, U-band-enhanced wide-area photometric dataset in the HSC-SSP Deep and ultraDeep fields.However, because of their developmental nature, we did not make these earlier catalogs publicly available to the community.Moreover, the depth of grizy imaging used in their creation was significantly shallower than what has become available as the HSC-SSP survey accumulated data over time.Finally, for simplicity, our earlier catalogs did not include any NIR information even though deep NIR imaging.
The goal of the present paper is then to provide wellcharacterized catalogs that contain U + grizy (+NIR, where available) photometry, photometric redshifts in the Deep and ultraDeep areas of the CLAUDS+HSC-SSP surveys.To build these catalogs, we use two sets of codes.We use recently developed photometry and photo-z methods, namely a modified version of the HSC-SSP photometric pipeline (hscPipe, Bosch et al. 2018) for photometry extraction and Phosphoros (Paltani et al., in prep.;Euclid Collaboration 2020) for photo-z computation.We also provide photometry and photometric redshifts measured with widely used and long-established codes, namely SExtractor (Bertin & Arnouts 1996) and Le Phare (Arnouts et al. 2002;Ilbert et al. 2006).The comparison of the results is allowing the validation of the different methods.With the present paper, we are making both catalogs publicly available for the community.This work shows the overall agreement between both catalogs in terms of photometry and photometric redshifts.The availability of the two catalogs allows the users to select the one most suitable for their own interests or to use them in tandem to assess the reliability of results as advocated by Weaver et al. (2022) for COSMOS2020.
Throughout the paper, magnitudes are provided in the AB system.CLAUDS u and u * bands are referred to as U bands when no distinction is made between them (see Sect. 2.1 of this paper or Sect.2.1 of Sawicki et al. 2019).We also assume the Planck Collaboration XIII (2016) cosmology with H 0 = 67.74km s −1 Mpc −1 , Ω M = 0.3089, and Ω Λ = 0.6911.

Overview of photometric data
The data consist of four fields (E-COSMOS, DEEP2-3, ELAIS-N1, and XMM-LSS; see Tables 1 and 2) and come from four different surveys.All the fields are covered by CLAUDS (Sawicki et al. 2019) and the HSC-SSP (Aihara et al. 2018a), but E-COSMOS and XMM-LSS are also partly covered by two nearinfrared (NIR) surveys, UltraVISTA (McCracken et al. 2012) and VIDEO (Jarvis et al. 2013), respectively.The transmissions of the photometric bands of these data are shown in Fig. 1, while their coverage and depth are presented in Fig. 2 and Table 2.
CLAUDS is a deep survey using the Canada-France-Hawaii-Telescope (CFHT) MegaCam imager (Boulade et al. 2003) in the u and u * bands 1 .This program combines observing time from Canada, France, and China and consists of a total  The HSC-SSP survey is a recently completed survey conducted with the Hyper Suprime-Cam camera (Miyazaki et al. 2018) on the Subaru telescope in Hawaii.In this work, we consider the public data release 2 (PDR2) Deep and ultraDeep layers of the survey presented in Aihara et al. (2019).The four fields have been observed in the g, r, i, z, and y bands (Kawanomoto et al. 2018) with median depths ranging from m g = 27.3 to m y = 25.3 for point sources detected at 5 σ (all the depths are presented in Aihara et al. 2019, Table 2).
The UltraVISTA survey is a NIR survey covering 1.5 deg 2 in the central region of the E-COSMOS field.We use the publicly available data release three images in the Y, J, H, and K s bands 2 .The images were acquired by the VISTA Telescope (Emerson et al. 2004) with the VIRCAM instrument (Dalton et al. 2006).The data reach depths of m Y ∼ 25 and m J,H,K s ∼ 24 for 2 ′′ -diameter apertures at 5 σ (McCracken et al. 2012).
2 https://ultravista.orgFor the XMM-LSS field, we use the public VIDEO survey data release 4 images in the NIR Y, J, H, and K s bands 3 .As for UltraVISTA, the images were produced using the VIRCAM instrument on the VISTA Telescope.The VIDEO data cover around 4.5 deg 2 , down to depths ranging from m Y = 24.6 to m K s = 23.8 for 5 σ in 2 ′′ apertures (Jarvis et al. 2013).
A summary of the data available in the different fields is presented in Table 2.The depths shown in the table are the medians of the 5σ scatter measured by us in 2 ′′ apertures thrown in the background of the images in all the patches (see Sect. 2.3 for patches definition).We note that our reported depths can differ from the ones announced by other teams, depending on the adopted definition.Indeed, the depths computed from point sources are deeper than those computed from the standard deviations of background fluxes in 2 ′′ apertures.
We also include private spectroscopic observations from the spectroscopic follow-up of faint low redshift galaxies (r ≤ 22 and z ≤ 0.5) in three of the HSC fields with the Hectospec spectrograph on the MMT telescope (Cheng et al. 2021); and from the spectroscopic redshift catalog from the COSMOS team (M.Salvato, private communication), which consists of several optical and NIR spectroscopic follow-ups of X-ray to far-IR/radio sources, high-redshift star-forming and passive galaxies, as well as poorly represented galaxies in multidimensional color space (C3R2, Masters et al. 2019).We also include the low-resolution spectroscopic redshifts from slitless spectroscopy with the near-infrared grism from the 3DHST survey data release v4.1.5using only those classified as secure grism redshift measurements (Momcheva et al. 2016;Skelton et al. 2014).
For all the above redshift surveys, we compare our photometric redshifts with the most secure spectroscopic redshifts only, identified with high signal-to-noise and with several spectral features (i.e., equivalent to VVDS or VIPERS redshift flags 3 and 4).

Data preparation
The PDR2 of HSC-SSP has been processed with version 6 of the HSC/LSST pipeline (Bosch et al. 2018;Aihara et al. 2019).We process the U and NIR bands using the same version of the pipeline, which has been modified in order to process non-HSC data.
The first step is to project the images onto the HSC reference pixel grid.As presented in Bosch et al. (2018), the hscPipe divides the sky into tracts, themselves divided in patches.A patch has a resolution of 4200 × 4100 pixels, with a pixel scale 3 http://www.eso.org/sci/observing/phase3/data_releases.htmlA82, page 3 of 25 A&A 670, A82 (2023)  of 0.168 ′′ pixel −1 .A tract is composed of 9 × 9 patches, all having a 100-pixel overlap with the adjacent patches and sharing the same tract coordinate system.The tracts form a mosaic on the sky, with some overlaps varying depending on the declination.The U-band data reduction uses the same HSC data format (i.e., tract patch and pixel grid) for the final products.The full details of U-band image processing are given in Sawicki et al. (2019), but we review the key steps below.The detrended Mega-Cam images were processed using the MEGAPIPE data pipeline (Gwyn 2008) for astrometric and photometric calibration and for the stacking of the images.Two versions of the CLAUDS images are produced, one using the Pan-STARRS astrometry and another one using the more recent Gaia astrometry (Gaia Collaboration 2016a,b).The processing done with the hscPipe described in Sect.3.1 uses the Pan-STARRS calibration for the U band and the processing done with SExtractor described in Sect.3.2 uses the images calibrated with Gaia.Individual images underwent photometric calibration before the production of stacks, generated directly on the HSC pixel grid.The images were resampled using Swarp (Bertin et al. 2002), the background subtracted using a 128 × 128 pixel mesh grid done with Swarp, and stacked to produce the images and weight maps matching the HSC patches.In the case of the NIR data, the fully calibrated public mosaics (image and weight maps) are projected onto the HSC pixel grid with Swarp for all the patches covered by the NIR data.For the VIDEO data, some tiles are duplicated due to the layout and overlap of the different observations.The duplicated tiles are combined by averaging them over the overlaps to obtain a single tile per patch.The U and NIR band processed data are formatted into coadd files that can be handled by hscPipe (see Bosch et al. 2018).
A82, page 4 of 25 Fig. 3. Flow chart and illustration of the process of mask creation.Images are smoothed and converted into binary masks with a threshold to identify bright extended objects.To avoid the masking of bright extragalactic objects, the masks are eroded (the pixels at the edges are set to 0), which removes small-scale objects.The masks are then dilated to recover the full size of the original masks for bright stars.Contours are created around masked areas (plotted in green) and we compare our new masks with the Coupon et al. (2018) ones (displayed as white circles and rectangles).This is performed with a module developed for hscPipe4 .For each patch and band, the corresponding coadd file is gathering the images, and the variance map is computed from the weight map.A mask of the bright stars found in the patch is created using the prescription of Coupon et al. (2018).The point spread function (PSF) of the patch image is estimated.Finally, the coadd file is provided to the hscPipe for source detection and photometry extraction.

Masking
2.4.1.Bright-star masks Aihara et al. (2019) identified some issues in the bright-star masks in the HSC PDR2, due to the change in the sky subtraction algorithm from HSC PDR1 (Aihara et al. 2018b).An outcome of this change is that the bright-star masks based on Gaia and Tycho-2 catalogs presented in Coupon et al. (2018) are not conservative enough for PDR2.Moreover, stars in the NIR bands suffer more diffusion, implying that their halos are not completely covered by the masks.In addition to this problem, some bright stars are missing from the masks used by hscPipe5 , requiring some improvements in the masking procedure.To solve these two issues, we combine the original masks with new masks derived using a simple image-based approach that follows the flowchart described in Fig. 3. Swarp is used to smooth the images suffering the most scattering (H-band if available, HSC y-band otherwise) using the LANCZOS3 method down to ∼5 ′′ pixel scale (Bertin et al. 2002).Then a threshold is applied to the pixel flux values to identify very bright objects and the images are converted onto binary masks (1 for masked and 0 otherwise) matching approximately the HSC ones.Finally, we apply to the binned masks a 5 × 5 erosion-dilation process.This process removes small masked areas and smooths the edges of the larger ones.Contours are converted into DS9 polygon regions that can be handled by the VENICE code6 (Coupon et al. 2018) and applied to the source catalogs.
As shown in the comparison with the masks of Coupon et al. (2018; bottom left plot of Fig. 3), our image-based approach does not mask all stars, but manages to mask the missing bright (m g < 16.5) stars and increases the size of the mask that were undersized compared to the observations.The final bright-star masks are built by combining the Coupon et al. (2018) masks with those we generated.

Satellite trails
Despite the update of HSC-Pipe that identifies and clips transient artifacts before coaddition, PDR2 images (as well as CLAUDS images) still exhibit an important number of satellite trails (≫10 deg −2 ).These trails could be easily handled via median-type stacking, but this would lead to a depth loss of 0.3 mag.An alternative is to use artifact rejection, but it misses a significant number of trails, which affects both source detection and photometry estimation.Indeed, SExtractor often fails at detecting sources close to bright satellite trails or generates spurious detections in those cases.This creates a bias in the number of counts per surface element.In addition, even when objects are detected around satellite trails, the photometry in the band with the satellite trail can be biased.This bias will propagate in the redshift estimation accuracy by generating catastrophic failures.
We decided to flag the regions possibly affected by these trails.We used the convolutional neural network (CNN) Maxi-Mask (Paillassa et al. 2020) to detect satellite trails in astronomical images.We used 0.6% prior for satellite detections and thresholded the probability mapping to the minimum to get most of the trails.On 4200 × 4200 images, 0.6% corresponds to approximately two satellite trails of 10 pixels width crossing the image diagonally.Maxi-Mask can mismatch edge-on galaxies with track and generate false satellite detections that are significantly smaller than the satellite trails.We get rid of these false detections by removing any detection shorter than 1 ′ .
We only applied the code to the detection images, as most bright satellite trails can be detected from them.To avoid galaxies with photometry contaminated by the trails, we perform a 2 ′′ binary dilation on the satellite mask images and use VENICE to convert these binary images into the ST_TRAILS flag in the released catalogs.

Source detection and photometry
Source extraction and photometric analyses are two crucial steps to optimize the detection of faint sources in deep and crowded regions and to estimate the colors required for the photometric redshifts.As mentioned by Weaver et al. (2022), photometric redshift estimates are more sensitive to the input flux catalogs than to the photo-z codes themselves.Hereafter we perform these two steps with two different codes: the one implemented in hscPipe (Bosch et al. 2018), which has been adapted to handle images from other facilities (e.g., U and NIR images), and the well tested and widely used SExtractor software (Bertin & Arnouts 1996).In the following part we describe the main steps followed by each software to generate the multiband catalogs and we compare the photometry obtained by the two catalogs as well as with the external catalog from COSMOS2020 (Weaver et al. 2022).

HSC pipeline
The source detection process using the hscPipe is described in detail in Bosch et al. (2018).It consists first of source detection in all individual bands.The configuration to run the detection on both U-band images is the same as that used in Aihara et al. (2019) on the HSC PDR2.However, for the NIR-data, a custom detection threshold was required to perform efficiently7 , as the default configuration led to an important number of false detections of background fluctuations in these bands.Then, the lists of detected objects are merged and the photometry of all sources is extracted in all bands.Two measurements of the photometry are produced for the detected sources: one independent in all the bands, and one forced photometry across all the bands; we use the latter in the following sections.The full raw hscPipe data products, divided in tracts, patches (see Sect. 2.3), and bands, are available through the HSC data access portal (see Sect. 5).However, those data contain redundant sources due to the overlap between tracts and between patches that need to be cleaned to construct the final catalogs.In the E-COSMOS field, both U-bands were treated as two different sets of images, with detection and flux measurements made in both bands.
For each field, we include several photometric flux measurements: 2 ′′ and 3 ′′ aperture photometry, PSF photometry, Kron photometry (Kron 1980), and cmodel photometry (Abazajian et al. 2004).Among these measurements both cmodel and PSF photometries take into account the PSF variation across patches and bands.Quality flags from the output data are also included for each band to ensure the quality of the photometry.Sources identified by the hscPipe as duplicated are removed.The masks for bright stars used by the pipeline are the same as those used in Aihara et al. (2019) and thus suffer from the same issues as discussed in Sect.2.4.1.We thus apply the updated star masks presented in Sect.2.4 as well as the satellite trail masks using the VENICE code (Coupon et al. 2018).
We check the consistency of our photometry with that obtained in Aihara et al. (2019) for the HSC PDR2 release.This comparison is done in the center of the E-COSMOS field where all the bands are available and where pipeline issues are most likely to happen, in particular in the processing of non-HSC data.The photometry of the hscPipe is in excellent agreement with the HSC PDR2 one.No bias is observed with magnitudes, and the scatter remains below ∼0.011 mag in all the HSC bands.This is also true for the cmodel magnitudes, where the scatter never exceeds 0.04 mag.
A compactness flag is computed in each band by comparing the PSF photometry to the cmodel one.A source that meets the following criterion: is considered as compact in band X.If this criterion applies to all HSC bands, the source is considered as compact and flagged accordingly.

SExtractor
We also produce a version of the HSC-CLAUDS catalogs with the detection and photometry based on the SExtractor software (Bertin & Arnouts 1996).
The source detection is performed on a detection image constructed by a nonlinear combination of the U + grizy bands.When near-infrared observations are available, the detection image also includes the K s band.To compute the detection image we adopt the χ mean combination proposed by Drlica-Wagner et al. (2018) as an alternative to the standard χ 2 image (Szalay et al. 1999), because it minimizes the discontinuities between regions with a different number of input images.
The detection image, χ mean , is defined as : and f i is the image in the band i, σ 2 i its variance and n the number of input images used in the combination.The χ mean detection images are shown in Fig. 2.
Before constructing the χ mean images, we first checked the homogenization of the variance maps, as the procedures used to generate the images and their associated weight maps are different for the different datasets (CLAUDS, HSC, VIDEO or UltraVISTA).To ensure that the weight map properly scales as a variance map, we measured the standard deviation, σ(Im), of 5000 randomly thrown 2 ′′ apertures in the background of the image, Im, and compared with the median value, median(Var), A82, page 6 of 25 G. Desprez et al.: Combining the CLAUDS and HSC-SSP surveys obtained in its associated variance map, Var.For each patch and passband, we rescaled the variance map by the factor: σ 2 (Im) /median(Var).
For CLAUDS and HSC bands, the scaling factors are around ∼2.5-3, which means that without the rescaling we would underestimate the photometric errors by a factor ∼1.5.For the NIR data, the variance maps are rescaled by a factor ∼4 in the deepest areas and up to ∼10 in the shallowest areas, implying an underestimation of the photometric errors by a factor ∼2-3.
The source detection is performed on the χ mean images convolved with a Gaussian with a 4 × 4 pixel kernel.A source is detected if a minimum of 10 contiguous pixels is above a signal-to-noise ratio of 0.8σ per pixel.The source deblending is controlled by the DEBLEND_MINCONT parameter, setting the fraction of the total intensity to be considered as a separated source, while the number of deblending thresholds is set by the DEBLEND_NTHRESH parameter.We found that the values of 0.0003 and 64 for the two parameters provide efficient deblending.
The detections and flux measurements are performed using the SExtractor dual mode.The χ mean image is used to detect and measure the shape parameters, while flux measurements are performed in individual images.For each object, we provide two fixed aperture magnitudes (MAG_APER) with 2 ′′ and 3 ′′ diameters and the Kron magnitude ( MAG_AUTO Kron 1980).The Kron magnitude automatically adapts to the first-order moment of the light distribution, and its elliptical aperture is scaled to provide a pseudo-total magnitude (Bertin & Arnouts 1996).However, fixed-aperture magnitudes are expected to be less noisy for faint sources and less impacted by blended sources, offering less noisy colors, a key input ingredient for the photometric redshift estimates (Sawicki et al. 1997;Hildebrandt et al. 2012;Moutard et al. 2016).As proposed by Moutard et al. (2016), for each object, we compute a single offset (the same for all the bands) that allows us to convert from aperture to total magnitudes, that are required for physical parameters measurements: m i = m APER,i + δ m , with the offset defined as: where i runs through the Ugrizy and J, K s filters and w i , the weight, defined as 1 where f AUTO and f APER are the fluxes associated to the Kron and aperture magnitudes respectively and σ AUTO and σ APER are their respective uncertainties.Being defined as a weighted average over all bands, the offset is more robust than if it were estimated for each band individually.However, some objects can have an unrealistic error in one band, which then accounts for more than 99% of the total weight.To prevent this situation, we iteratively set to 0 the individual weights representing more than 95% of the total weight.The offsets for the 2 and 3 ′′ aperture magnitudes are given in the final catalog.All the magnitudes provided in the catalog are corrected for galactic extinction based on the Schlegel extinction maps (Schlegel et al. 1998).

Comparison and quality assessment
The HSC-CLAUDS photometry is assessed by comparing the measurements of our hscPipe and SExtractor catalogs.We also compare to a third and distinct method, delivered in the COSMOS2020 catalog (Weaver et al. 2022), based on Farmer software (Weaver et al., in prep.) which uses The Tractor (Lang et al. 2016) to measure the photometry.Farmer performs model photometry of individual sources with a simultaneous optimization of a group of nearby sources with additional constraints from the multi-channels.The COSMOS2020 catalog uses, among others, the same CLAUDS, HSC and UltraVISTA data as we did.The comparison is thus performed in the central E-COSMOS field, which is also the only region where all filters are available.
In Figs. 4, 5 and 6, we compare the magnitudes.The cmodel from the hscPipe with the SExtractor aperture corrected magnitudes (Fig. 4); the cmodel with Farmer magnitudes A82, page 7 of 25 (Fig. 5) and the SExtractor corrected aperture with the Farmer magnitudes (Fig. 6).The panels show for each filter the magnitude differences as a function of magnitude for the matched sources and the median of the distributions (solid lines) and their ±34% intervals around the median (dashed lines) down to the 5σ depth measured in the ultraDeep region of E-COSMOS.Several points are noteworthy.First, in the optical bands (bluer than y band), we observe small offsets (≤−0.1 mag) in the bright end, with SExtractor magnitudes being systematically brighter than hscPipe and Farmer photometry, while almost no offset is observed between hscPipe and Farmer.The two latter magnitudes are based on modeled magnitudes and are not sensitive to seeing variations between bands, in contrast to the SExtractor aperture corrected magnitudes (see Eq. ( 3)), where no PSF homogenization has been applied, and the mean offset can be slightly over-estimated due to the worse seeing images.
In the optical bands also, hscPipe magnitudes are getting gradually brighter than SExtractor and Farmer magnitudes for the faint objects (m ≳ 25.5).The effect is much less pronounced between SExtractor and Farmer photometry.This faint-end trend could arise if the background correction in hscPipe is underestimated.This can be seen as well, although less pronounced, when comparing 2 ′′ aperture magnitudes between SExtractor and hscPipe optical magnitudes (Fig. C.1).In contrast to the optical bands, in the near-infrared bands (redder than the Y band), the three different sets of magnitudes are consistent with only small systematic biases (≤0.05) with no trend with magnitude.
Finally, for all the bands, the distribution of the difference between the cmodel and SExtractor magnitudes shows a skewness with a large negative tail for the brighter SExtractor magnitudes.As a modeled photometry, the hscPipe may better handle the contribution of neighbors than the SExtractor aperturecorrected magnitudes.Consistently with this interpretation, the distributions for the two modeled magnitudes, Farmer and hscPipe (Fig. 5 while the distributions between Farmer and SExtractor also show a slight asymmetry in the same direction. Since the photo-z computation is made using 2 ′′ aperture photometry with both catalogs (see Sects.4.1 and 4.2), we make a comparison of the color apertures in Fig. 7.We show the color differences between SExtractor and hscPipe photometry for six colors, namely u − g, u * − g, g − r, i − y, z − K s , and J − K s .We note the good agreement between the colors provided by both catalogs up to magnitude ∼24 in all the tested colors with no bias and a scatter ≲ 0.1.For fainter objects, a bias appears in colors involving optical and NIR bands, accompanied by a strong increase of the scatter.This can be due to the difference in background treatment, but also to the differences in depths.For i − y and z − K s colors, this difference is ≳1 mag.The large scatter observed would be due to strong noise in the NIR band magnitudes of the faint objects compared to the optical ones.
To check the consistency of the photometric calibration between the four different regions of HSC-CLAUDS, we selected the stellar population and compare their stellar loci in color diagrams.We perform this test only with the hscPipe catalog, where a PSF magnitude is provided to get optimal stellar colors.Figure 8 displays the contours encompassing 68% of the stars in four color-color diagrams for all the fields (color coded) available in a given color combination.The stellar loci for the different fields show an excellent overlap between each other in the optical bands.Combined with NIR color (bottom right panel, z − K s vs. i − z), the stellar locus in the XMM-LSS field shows a larger scatter and a shift of ∼0.05) compared to the E-COSMOS field.The figure shows that the observed loci are not properly aligned with the main-sequence star colors computed from the Pickles template library (Pickles 1998).This indicates a need for zero-point calibration to correct colors.From the diagrams, it appears that the U bands and the NIR bands require the largest corrections.Although, for the bluest bands, the shift can also be due to Milky-Way extinction which is not accounted for because of the difficulty of computing it for objects (stars) inside the Galaxy.We conclude that the photometry in the optical bands is consistent across the four fields but requires further correction in the form of zero-points calibration.The NIR band photometry might present some differences between E-COSMOS and XMM-LSS fields, which can be taken into account in the analysis by adjusting the zero-points on a field-to-field basis.
Figure 9 shows the number counts in all bands for the SExtractor and hscPipe catalogs.The number counts in the two ultraDeep regions (COSMOS-UD and XMM-LSS-UD) are presented separately, while those from the deep regions are the averages of the four fields.The comparison of the ultraDeep region of E-COSMOS and COSMOS2020 shows a good agreement, especially in the HSC bands, which is expected as the data used to build the catalogs are the same in the three cases.More differences arise in the NIR bands, first due to the difference in data, as COSMOS2020 uses the DR4 of UltraVISTA and our catalogs are built with the DR3, which is shallower.The number counts in Fig. 9 show the counts after applying a cut at S /N > 3 in the concerned bands, which peak at brighter magnitudes than the COSMOS2020 counts.This could be explained by the rescaling of the errors, which affects strongly the NIR A82, page 10 of 25 Notes.The corrections are given in magnitudes, to add to the measured one.We chose that by convention the i-band correction would be null.
For the Phosphoros computation, the corrections in the E-COSMOS and XMM-LSS fields were measured in both Ugrizy and Ugrizy+NIR configurations.
bands, leading to larger errors and hence pushing some objects below the S/N cut.The ultraDeep uncut counts match rather well with the COSMOS2020 ones, which is further evidence that the discrepancy is due to the cut in S/N.For the bluest bands, there is a shift between the COSMOS2020 counts and ours.This shift could indicate some uncertainties in the estimation of the area that would produce a vertical shift, or systematic overestimation of the flux that would cause a horizontal shift, but it could also be explained by the differences in detection schemes used in the different pipelines.The hscPipe performs source detection in all the bands, while SExtractor uses a χ 2 image constructed with bands from U to z, with the addition of the K s band when available.The detection in COSMOS2020 uses only the reddest bands from i to K s .This difference makes our catalog more sensitive to blue objects, which explains why the discrepancy between the curves in Fig. 9 is stronger in the blue bands at faint magnitudes.

Photometric redshifts
The photometric redshifts (photo-zs) are estimated with two template-fitting codes for the two photometric catalogs based on the 2 ′′ aperture flux in all bands.The hscPipe catalog is processed with the Phosphoros code, created to be part of the Euclid photo-z pipeline (Paltani et al., in prep.; Euclid Collaboration 2020), while the SExtractor catalog is processed with the well-tested code, Le Phare (Arnouts et al. 2002;Ilbert et al. 2006).As this is the first application of Phosphoros on deep imaging surveys destined to be released, we also run Phosphoros on SExtractor photometry to distinguish the differences related to photometry issues from template-fitting code issues; this is detailed in Appendix D.
For both codes, several configurations were considered and tested.We explored using different templates, prior information, photometry, and errors.The configuration presented in this section is the result of this exploration.We make the choice to use similar configurations, most of the differences being due to differences in the parameters both codes can adjust.
Both code use the templates library of Ilbert et al. ( 2013) for galaxies, consisting of a total of 33 spectral energy distributions (SED): seven elliptical galaxies, twelve spiral galaxies (Polletta et al. 2007;Ilbert et al. 2009), twelve starburst galaxies, and two elliptical SEDs generated with the Bruzual & Charlot (2003) stellar population synthesis models.Intrinsic extinction is also configured in a similar manner in both codes.Extinction was added as a free parameter for a reddening excess E(B-V) = 0, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.4, 0.5 and by considering four attenuation laws: starburst-like Calzetti et al. (2000) and two dusty versions of the Calzetti law (with a bump at 2175 Å) for templates bluer than SB3, SMC-like Prevot et al. (1984) for templates redder than SB3 and no extinction for templates redder than Sb.
Emission lines are added to the templates by assuming an empirical relation between the UV light and Hα line flux (Kennicutt 1998) and ratios between the different hydrogen lines and oxygen lines as measured from high redshift spectroscopic surveys (Ilbert et al. 2009) or SDSS survey (Paltani et al., in prep.).
As proposed by Ilbert et al. (2006), to account for photometric variations in the four HSC-CLAUDS fields, the spectroscopic redshifts of bright sources (19 ≤ i ≤ 21.5) are used to apply a band-per-band zero-point correction in each field separately.These corrections are given in Table 3 for each field and combination of bands used.A small difference between the two configurations is that Phosphoros computed the zero points for both the Ugrizy and Ugrizy+NIR configurations in the E-COSMOS and XMM-LSS fields.Also, for both codes, the photo-zs computed in the central part of the E-COSMOS field use both u and u * data separately when available.

Phosphoros
Phosphoros is a template-fitting code that implements all the main functionalities of other template-fitting implementations.The focus on Phosphoros is the fully Bayesian treatment of the multi-dimensional likelihoods, with the possibility to apply complex priors and to marginalize over any parameters.Phosphoros uses the configuration provided in Sect. 4.
Phosphoros uses fluxes uncorrected for the Milky Way extinction, as the Galactic extinction is handled internally by Phosphoros by applying it directly on the fitted SED, following the prescription of Galametz et al. (2017).Photometric errors in the different bands x, σ f,x , are modified using the A82, page 11 of 25 A&A 670, A82 (2023) following equation: where f is the flux, x the band and α is the uncertainty rescaling factor and is 2 for Ugrizy bands and 4 for NIR bands, and β corresponds to a systematic uncertainty (which could come from the photometry or from the models) and is set to 0.03 for all the bands.We also consider several prior information.When computing photo-zs using the NIR bands, a top-hat prior in absolute magnitude is applied, allowing only results with −24 < M g < −8, as well as a volume prior to take into account that the probability of a source to be in a slice of redshift increases with the volume of the slice.When the number of available bands is restricted to the optical wavelengths, we apply a more informative prior based on luminosity functions from Zucca et al. (2006).The rest-frame templates are sorted into four types according to their B − I colors.For each type, we construct a prior assuming a Schechter function (Schechter 1976), with the α, the index of the power law at low luminosity, and ϕ * , the source density, parameters given by Zucca et al. (2006) in the rest-frame B-band for galaxies between redshift z = 0.4-0.9.However, the M * values are modified, setting them 3 magnitudes brighter, to avoid an overly strong suppression of luminous galaxies at high redshift.We introduce this modification because the strong evolution of M * makes the values determined by Zucca et al. (2006), which are essentially obtained at low redshift, not applicable even at moderately high redshifts.Any small variation of this shift value (±0.5 mag) only has a small impact on the results.We prefer this prior to that on the magnitude-redshift distribution from Benítez (2000) because its shape is physically better justified, which is especially important at high redshifts.
For each object, the probability distribution function of the redshift (PDZ) is obtained by marginalizing the posterior on the redshift dimension.The PDZs are defined with a fixed redshift grid (z = 0-6) with a step δz = 0.01.The photo-zs point estimates are defined as the median values of the PDZs.To deliver a consistent dataset across the entire HSC-CLAUDS survey, we provide a complete set of photo-zs with their PDZs based on the Ugrizy photometry alone, and an additional set based on UgrizyY JHK s in the regions where near-infrared bands are available.
We also provide a star/galaxy separation by fitting the sources with stellar templates from the Pickles (1998) library.In the catalog, all object with χ 2 star < χ 2 gal are flagged as potential stars.Objects that are both potential stars and compact according to Eq. ( 1) are flagged as stars.

Le Phare
Le Phare (Arnouts et al. 2002;Ilbert et al. 2006) is a well known and standard template-fitting method.With Le Phare, we adopted the same configuration as with Phosphoros, in terms of the SED library, attenuation laws and reddening excess, emission lines and zero-point calibrations (see above, Sect.4).
The SED fitting is run on SExtractor 2 ′′ aperture photometry, corrected for the Galactic extinction.Photometric errors are amplified by a factor 1.5, on top of which systematic magnitude errors of 0.02 mag and 0.05 mag are added in quadrature for optical and NIR bands respectively.Depending on whether NIR photometry is used, a different prior set is adopted.In the case of Ugrizy+NIR, a g-band absolute magnitude is applied to reject solutions with M g < −24.When only the optical bands are used, a fine-tuned prior similar to the Benítez (2000) prior is used (Ilbert et al. 2006).
The redshift output of Le Phare is the median of the PDZ as point estimate, or, when the PDZ is not defined, the redshift is associated to the lowest χ 2 value in the grid (the redshift grid is the same as the Phosphoros one, see Sect.4.1).The photozs errors are computed from the 68% confidence interval of the PDZ.
In addition to the galaxy template library, we also run the SED fitting with stellar and quasar template libraries.The stellar library includes normal spectral types at solar abundance, metalweak and metal-rich F-K dwarfs and G-K giants (Pickles 1998), low mass stars with different effective temperatures, gravity, and types (Chabrier & Baraffe 2000) and white dwarfs (Bohlin et al. 1995).The quasar library is a compilation of observed quasar templates (Polletta et al. 2007) and templates with different contribution of galaxy and AGN spectra (Salvato et al. 2009(Salvato et al. , 2011)).A first star/QSO/galaxy classification is provided, where compact objects with best fit obtained for a star or QSO template are flagged as such (see Appendix E).

Comparisons
In this section, we compare our photometric redshift estimates with spectroscopic redshifts (described in Sect.2.2) and the latest 30-band photometric redshifts in the COSMOS field (COS-MOS2020, Weaver et al. 2022).As mentioned in Weaver et al. (2022), more than the differences in the photometric redshift codes, the main differences in the photo-z vs. spec-z comparisons are due to the input photometric catalogs.This is especially true in our work where we use two SED fitting codes with similar approaches but two quite different photometry codes.For this reason in the following we perform the comparisons using the photometric redshifts derived with the combinations hscPipe+Phosphoros (hereafter HPH) and SExtractor+Le Phare (hereafter SLP).In Appendix D, we show that the results from the two photo-z codes are similar when run on the same SExtractor catalog (see Fig. D.1).
The photometric redshift point estimates are delivered with an uncertainty based on the 68% confidence interval.To assess the relevance of the errors, for each galaxy i, we compute the absolute scaled residuals: where δ 68% i is the 68% confidence interval for the galaxy i.The cumulative distribution, D z , should reach 0.68 for D z = 1.
We also assess the quality of the PDZs with the Probability Integral Transform (PIT) statistic (PIT plot, Dawid 1984;D'Isanto & Polsterer 2018).It is based on the shape of the distribution of the cumulative probabilities (CDF) at the true A82, page 12 of 25 Fig. 10.Comparisons of the photo-zs based on the Ugrizy photometry with spectroscopic redshifts down to i ≤ 22.5 in the four HSC-CLAUDS deep fields.Top panels show the photo-zs computed by Phosphoros code and the hscPipe catalogs, while the bottom panels show the photo-zs computed with Le Phare and SExtractor catalogs.The solid gray lines show the 1:1 line and the dashed lines show the z photo = z spec ± 0.15(1 + z spec ) corresponding to the limits for catastrophic failure.value: for galaxy i with spectroscopic redshift z i and probability distribution function PDZ i .If the spec-z's can be randomly drawn from the PDZs then a flat PIT distribution is expected, suggesting that the PDZs are statistically correct, whereas convex or concave PIT distributions point to under-or over-dispersed PDZs, respectively (Polsterer et al. 2016).

Comparisons with spectroscopic redshifts
Figure 10 presents density plots comparing the photo-zs from the HPH (top panels) and SLP (bottom panels) runs to spec-z's for sources with m i < 22.5 in the four fields separately and based only on the Ugrizy photometry.A zero-point calibration has been applied to each field separately.The performance appears similar in the two configurations, with a precision varying between σ = 0.02-0.04 and outlier fractions η = 2-3% across the four fields.However, we note a systematically higher scatter for the HPH results compared to the SLP ones by ∼0.01.In Fig. 11, the four HSC-CLAUDS fields are combined together and the Ugrizy photo-zs are compared with spectroscopic redshifts in faint magnitude bins.However, we note that the ELAIS-N1 field has almost no spectroscopic redshift with m i ≥ 22.5 and the DEEP-23 field only down to m i ≤ 24.The vast majority of the faint spectroscopic sources come from the XMM-LSS and E-COSMOS fields.
The precision of the photometric redshifts gradually deteriorates with magnitude, passing from σ ∼ 0.03 at m i ≤ 22.5 to σ ∼ 0.09 at m i ≤ 26, with slightly worse results for the HPH implementation.Similarly, the fraction of catastrophic failures gradually increases from η ∼ 4% to η ∼ 29%.
Although the fraction is similar between the two implementations, there are differences in the location of the outliers.Due to the Balmer and Lyman break confusion, two approximately symmetric clouds of outliers are generally produced, one for low spec-z sources with a z phot ≥ 2, and the other one for high spec-z galaxies with a low photo-zs (z phot ≤ 0.3).The SLP implementation presents a visible cloud for sources with high spec-z's and low photo-zs, while the HPH one leads to fewer sources in this cloud, but presents a higher population at low spec-z's and high photo-zs than the SLP implementation.This difference is due to the application of a volume prior in the Phosphoros code since this volume prior favors the high redshift values for poorly constrained photometric sources.
A striking feature is present with spectroscopic sources around z spec ∼ 1, to which higher photometric redshifts z phot ∼ 1.5-2 are assigned.The vast majority of those sources show a red z − y color that can be interpreted either as caused by the 4000 Å break at z ∼ 1.5, or by the contribution of emission lines (e.g., H β and [O III] lines) at z ∼ 0.9 (putting the 4000 Å break in the r − i color).Without the near-infrared, this double solution cannot easily be broken by the N(z) prior and occasionally leads to favor the wrong solution.
In Fig. 12, we show the comparisons of photometric redshifts measured by taking into account the NIR data.All the performance metrics are improved and we still observe slightly better results for the SLP implementation compared to the HPH one.For the faint galaxy population, the scatter and the catastrophic rate are improved by almost a factor of two.Including NIR data also alleviates the degeneracy observed around z spec ∼ 0.9 with z phot ∼ 1.5-2.This illustrates the benefit of the NIR observations to provide robust photometric redshifts over the entire redshift range.
The point estimate photometric redshifts are given with an uncertainty based on the 68% confidence level interval derived

Spectroscopic redshifts
Fig. 12. Same as Fig. 11, but for the XMM-LSS and COSMOS fields only and including the near-infrared photometry.
from the individual redshift probability distribution function (PDZ).To assess the relevance of the errors, for each source we compute the absolute scaled residual, D z (Eq.( 6), Sect.4.3.1).
Figure 13 presents the cumulative distribution of D z for the two implementations (HPH and SLP) and the two imaging datasets (optical, optical+NIR).If the photo-z uncertainties are correct we expect the cumulative curves to reach the 68% level at D(z) = 1.Indeed, all the cumulative distributions cross the 68% level at higher D(z), suggesting that the estimated photo-z uncertainties are underestimated.In both configurations, the effect seems even worse when the near-infrared photometry is included, probably because of the additional constraints it provides on the PDZs.Despite the fact that scaling factors have been applied to enlarge the photometric errors (see Sects.4.1 and 4.2), they seem insufficient to provide realistic photo-z uncertainties.We note that SLP configuration seems to provide more reliable error estimates than HPH.As it is still the case when Phosphoros is run on SExtractor photometry, this difference is probably due to the code configurations, and in particular the rescaling of the errors, and not only to the photometry.We, therefore, caution  that the current photometric-redshift uncertainties are on average underestimated.
Alternatively, instead of the 68% confidence interval, we can directly estimate the quality of the PDZs with a PIT plot (see Sect. 4.3.1).In Fig. 14, we show the histogram of the HPH C i values.The concave shape of the histograms shows that the PDZs are too narrow, leading to an underestimation of the errors.Such behavior has already been reported with several template-fitting methods on similar data sets (Tanaka et al. 2018;Euclid Collaboration 2020).We also see that using the NIR data has the effect of constraining the PDZs even more.These observations are in agreement with the conclusions drawn from Fig. 13.

Comparison with 3DHST and COSMOS2020 surveys
At faint magnitudes, the spectroscopic samples are often selected with specific photometric criteria and are therefore not a representative sample of the overall photometric population (Masters et al. 2015).To address this issue, we compare our results to the 3DHST (Skelton et al. 2014) and COSMOS2020 Fig. 15.Photo-z comparison with 3DHST redshifts for galaxies down to m H ∼ 24 for HPH (top panels) and SLP (bottom panels) catalogs, with the optical+NIR (left panels) and optical only (right panels) bands.The metrics are reported in each panel.(Weaver et al. 2022) surveys.The 3DHST is a slitless lowresolution grism survey that provides spectroscopic redshifts of near-infrared sources down to m H ≤ 24.However, 3DHST is still impacted by a selection bias due to the success rate of reliable redshift identifications, which is strongly dependent on the galaxy type, and in particular on the strength of the emission lines.COSMOS2020 provides robust photometric redshifts over the entire color-space thanks to its 30 photometric bands, but these redshifts are less precise than the spectroscopic ones.These two samples are used to assess the quality of our photometric redshifts for the faint population.
In Fig. 15, we compare the photometric redshifts from the HPH and SLP catalogs with the 3DHST redshifts, down to m H ∼24. Two 3DHST fields overlap with the HSC-CLAUDS observations in the E-COSMOS and XMM-LSS fields and the sources are matched with a 1 ′′ tolerance.A small fraction of spectroscopic redshifts, especially when determined using slitless grisms, can be wrongly assigned; however, the impact should be negligible on the estimated scatter, as the median absolute deviation (MAD) is robust against outliers.The results are consistent with the trends observed in Sect.4.3.2 with optically selected spectroscopic sources.HPH photo-zs have slightly higher scatter and outlier fractions in the two adopted sets of filters.The catastrophic-failure fraction is reduced by a factor of two when including the NIR bands and the mismatch in the redshift domain 1 ≤ z ≤ 2 is noticeably reduced.On the code side, the implementation of priors used in Phosphoros for the optical configuration reduces the number of low photo-z outliers compared to those from the Le Phare code.
In Fig. 16, we show the comparisons with the COSMOS2020 photometric redshifts.As described in Weaver et al. (2022), four versions of photometric redshifts are available based on two photometric extraction softwares (Farmer and SExtractor) and two photometric redshift codes (Le Phare and EAzY; Brammer et al. 2008).In an attempt to be as independent as possible from  source extraction and photometric redshift methods, in Fig. 16 (top 2 rows) we first compare our measurements based on the optical+NIR bands with the COSMOS2020 Farmer/EAzY catalog.However we point out that we use the same deep images (HSC, CLAUDS, and UltraVISTA) as COSMOS2020; these images provide the strongest constraints on the photo-zs at faint magnitudes.But COSMOS2020 benefits from the ultradeep IRAC images and the narrow and intermediate Suprime-Cam bands (Taniguchi et al. 2007(Taniguchi et al. , 2015)).The former will reduce the number of catastrophic failures, while the latter will improve the scatter at bright magnitudes.In Fig. 16, the sources are restricted to cleaned flags (FLAG_COMBINED=0) and sorted according to the Farmer i-band magnitude (Weaver et al. 2022) to get the same number of sources in each magnitude bin.The comparisons show similar trends in the metrics to those in Figs. 12 and 15, with slightly better results for the SLP implementation.No systematic bias is observed up to redshift z ∼ 6 and down to m i ∼ 26 with a catastrophic-failure fraction below 15%.In the faintest magnitude bin, 26 ≤ m i ≤ 27, the scatter almost doubles and the catastrophic-failure rate becomes significant, reaching 30%; but those numbers are comparable to the level of consistency between the different photometric redshifts in the COSMOS2020 catalogs (Weaver et al. 2022, see their Fig.14).The metrics are in good agreement with the values reported for the spectroscopic redshifts up to m i ∼ 26 (see Fig. 12) and confirm that we do not expect major photometric-redshift issues for the whole galaxy population up to this depth.However two noticeable features are observed: first, a larger catastrophicfailure fraction at bright magnitudes with COSMOS2020 with respect to the comparison with spectroscopic redshifts; second, a mismatched photo-z population at low redshift (z FARMER/EAzY < 0.5) for which our two HSC-CLAUDS configurations estimate higher redshifts (above the 1:1 line).These two points are due to the fact that we are comparing our photo-zs with other photo-zs from COSMOS2020 that have their own inaccuracies.
To mitigate this effect, we make another comparison with the COSMOS2020 results by considering the four different versions of the photometric redshifts provided by Weaver et al. (2022).For all the sources with four photo-zs, we compute the mean redshift (⟨z⟩) and its standard deviation, σ(z).We then select  8) and ( 9)).C 20 F/E are the COSMOS2020 FARMER/EAzY photo-zs.
the COSMOS2020 self-consistent photo-zs by imposing a cut in the standard deviation σ(z) < 0.1(1 + ⟨z⟩).Strictly speaking, this does not constitute a guarantee to select secure redshifts, but rather it reduces issues related to either photometric extraction in problematic sources or sources with degenerate solutions, which may not be apprehended the same way by the two codes used in Weaver et al. (2022).
The comparisons with the mean COSMOS2020 photo-zs are shown in Fig. 16 (two bottom rows).With respect to the whole COSMOS2020 catalog, more than 85% and 70% of COSMOS2020 photo-zs are consistent with each other up to m i ∼ 25 and 26 respectively.Our own photo-zs also appear consistent with the mean COSMOS2020 photo-zs with a scatter lower than σ = 0.03 and a small fraction of catastrophic failures (η ≤ 5%).In the faintest magnitude bin, 26 ≤ m i ≤ 27, only 40% of the COSMOS2020 redshifts are consistent with each other (Weaver et al. 2022).For this population, our photozs are still in good agreement with a scatter σ ∼ 0.05 and a catastrophic-failure fraction not exceeding η = 10%.As already noted, we observed with the COSMOS2020 photo-zs that the distribution of catastrophic failures for the SLP and HPH catalogs are different.For the former, most of the catastrophic sources are high-z COSMOS2020 photo-zs which are assigned to low redshifts, while the latter mainly assign high redshifts to low-z COSMOS2020 photo-zs.
In Table 4, we report the scatter and the catastrophic-failure fractions for the different configurations for different i-band cuts tested against the spectroscopic sample and COSMOS2020.In the former case, the scatter and outlier fraction values are misleading, because the spec-z sample is strongly biased toward brighter magnitudes compared to the whole photometric catalog.Therefore we provide scatter σ and outlier fraction η values weighted by the number of photometric objects in each magnitude bin, which provides a much more realistic view of the performance on the full photometric sample (Euclid Collaboration 2020; Hartley et al. 2022).There is no such problem when comparing with COSMOS2020 since it contains essentially all photometric sources.The weighted σ and η for the spec-z sample are thus computed as follows: and: where j represents the magnitude bins (m i < 22.5, 22.5 ≤ m i < 24, 24 ≤ m i < 25, etc.), n j is the number of sources in the photometric catalog in the magnitude bin j, and σ j and η j are the metrics for bin j.The weighting scheme doubles the outlier fractions measured in the spectroscopic sample with a m i < 25 cut and triples it at m i < 26.We note that the scatter values with the weighted scheme are consistent with those of COSMOS2020 FARMER/EAzY.For the outlier fractions, we note some differences especially for the bright cuts, as by construction the results in the first magnitude bin do not change.It is with a cut at m i < 25 that we have the most consistent results.With this cut, we have a scatter σ ≲ 0.04, even when using only Ugrizy data and after weighting the spectroscopic sample results or using the COSMOS2020 data as a reference.With the same cut, we obtain an outlier fraction η < 0.05 for the spectroscopic sample, regardless of the photometry used.At m i < 25, the catastrophic-failure rate remains below 10%, both when compared with the spec-z sample and with COSMOS2020.

Redshift distributions
In Fig. 17  The different estimates lead to reasonably consistent redshift distributions, and are in agreement with the COSMOS2020 redshift distributions.The selection band impacts the overall shape of the distributions, with a bounded N(z) below z ∼ 2.5-3 with the u band selection and an increasing population at high redshift when adopting a redder selection band.
We also observe some specific differences.With the NIR selection, a bump at z ∼ 3.5 is visible for the point estimates in the HPH configuration.This bump might be due to poorly informative PDZs (nearly flat) skewed toward high redshift because of the priors, putting the median point estimate of the PDZ in this redshift range.This effect vanishes when adopting the sum of the PDZs.
With the U-band selection, we observe a significant difference in the shape of the N(z) between the photo-zs based on the optical bands and those using all bands, the former showing a higher fraction of sources between 1.3 ≤ z ≤ 2. This reflects the biased population seen in Fig. 11 (Ugrizy plots) and discussed in Sect.4.3.2.It means that these features are due to the lack of constraints in the absence of NIR data.However, the use of the stacked PDZ helps to reduce this effect.
Finally, we also observe differences between the SLP and HPH redshift distributions in the case of the all-band based photozs.At low redshift, the SLP distributions increase more rapidly while at high redshift they tend to decrease faster.While the Le Phare code does not apply the Benítez (2000) prior on the redshift likelihood when the NIR bands are included in the computation, Phosphoros always uses the volume prior which penalizes low redshift solutions and favors higher redshift ones.This is the same effect we observe in Fig. 16 for the faintest sources.

Data access/catalogs
Two catalogs are presented in this work.They correspond to two different processing methods applied to the same data, providing overall similar data products of equivalent quality as demonstrated in the previous sections.However, the two catalogs also present differences in the data products (e.g., types of photometry, PDZs, or classification) which can be suited for different purposes.We discuss here these differences that can lead to the selection of a particular catalog for different science cases, though, we recommend the use of both catalogs.
The SLP catalog is extremely similar to the Classic/ Le Phare version of the COSMOS2020 catalogs in its production.Therefore, for applications where a close comparison with COSMOS2020 results is desired, this catalog is the most suited.This catalog will also benefit from source physical parameters computed in a similar way as for the COSMOS2020 catalogs (Picouet et al., in prep.).
The HPH catalog provides different types of photometry, not present in the SLP one.For instance, it includes cmodel photometry and PSF photometry, that take into account PSF variation across the fields.On the photo-z side, it includes individual PDZ for each object in the catalog.Another interesting feature it the similarities between the HSC pipeline and the LSST pipeline.With its depth and wavelength coverage, the CLAUDS+HSC combination is close to the expected LSST ten-year dataset (Ivezić et al. 2019).Thus, the HPH catalog is most similar to the future LSST data.
The catalogs presented in this work will be publicly released and available for download on several sites8 .For the SExtractor photometry, several catalogs are provided with photometric redshifts estimated with Le Phare run with the 6 optical bands and all-bands configurations.
For the hscPipe photometry, all the catalogs are provided with the photometric redshifts estimated with the Phosphoros code.When run with the Phosphoros code, all the catalogs have their full PDZ available in separate files.
The full description of the columns for the SLP and HPH configurations are given in Appendix E together with some practical information not detailed in the core of the paper.This release is expected to evolve over time and other catalogs will be available.Additional information such as physical parameters derived with the Le Phare code (Picouet et al., in prep.) and photometric redshifts based on convolutional neural network (Pasquet et al. 2019;Ait-Ouahmed et al., in prep.) will also be delivered.

Summary
In this work, we presented the first public data release of the combined data over the 20 deg 2 of the CLAUDS and HSC surveys.For the four deep and ultraDeep fields of the HSC survey, we processed the CLAUDS data along with the HSC data.
For the E-COSMOS and XMM-LSS fields, we also added NIR data parts of the UltraVISTA and VIDEO surveys to extend the wavelength coverage of our catalogs, especially in the ultraDeep regions.
We produced catalogs containing the photometry extracted in two ways.One uses the well-known SExtractor code.The other uses a modified version of the hscPipe that allows the handling of the non-HSC data.CLAUDS and NIR data are preprocessed to match the HSC data format, then are treated by both codes to detect sources and extract their photometry across all bands.We then test the obtained photometry against one another and compare it to the COSMOS2020 photometry.Overall, we find that the photometry from the different methods is in good agreement down to mag ∼ 26 for Ugrizy bands, and down to mag ∼ 24.5 for the NIR ones.For fainter sources, the agreement declines, with increased bias and scatter for the photometry of both codes.For the hscPipe photometry, we found that the obtained photometry in the grizy bands is extremely consistent with that of HSC PDR2.However, the comparisons with our SExtractor photometry and the COSMOS2020 Farmer photometry show some differences probably related to the background determination and subtraction.hscPipe seems to systematically over-subtract the background, leading to flux loss for all the sources.We also found good agreement when comparing our galaxy number counts with those from Weaver et al. (2022) Farmer catalog.The strongest differences are found in the U and NIR bands, which are mainly due to the difference in detection strategy and error computation.
We computed photo-zs from both our photometry catalogs.Two template-fitting codes were used: Phosphoros was used to produce photo-zs from both hscPipe and SExtractor photometry, and Le Phare was used to process SExtractor photometry.The comparison of the different results shows that the major factor for the differences found between the photo-z catalogs is the photometry, which is similar to the observation made by Weaver et al. (2022) with their four different versions of COSMOS2020 photo-zs.Here, we found that the photo-zs from hscPipe and SExtractor photometry are in good agreement.Whatever the reference redshift they are compared to, we find a scatter σ ≲ 0.05 down to m i ∼ 25 or m K s ∼ 24.However, we find that photo-zs derived from hscPipe photometry present a systematically higher scatter, of around 0.01, probably due to the difference in the photometry, as both codes perform very similarly on the SExtractor photometry.Another difference in the results, this time attributed to the photo-z code configurations, is the systematic lack of low photo-z solutions from Phosphoros for faint galaxies, due to the priors used in its configuration, which for badly constrained PDZ is favoring higher redshift solutions.All the photo-z comparisons made show that Phosphoros is as efficient as Le Phare and ready to be used as a stand-alone photo-z code.
The catalogs presented in this work have been publicly released.For the hscPipe photometry, all the data from the images to the fluxes can be fully retrieved in the HSC data access portal.Otherwise, all the CLAUDS images and data product catalogs are accessible on the CLAUDS and DEEPDIP websites and will be updated with upcoming works on photometric redshift and physical parameter estimates.asymmetric scatter is present.This asymmetry seems to indicate that the hscPipe fluxes are lower than the SExtractor ones, and this effect is stronger for fainter sources.As these sources are more sensitive to background noise, we interpret the increase in asymmetry toward fainter magnitudes as a sign of stronger subtraction of the background by hscPipe compared to SExtractor.The same conclusion was reached when analyzing Fig. 5, 4 and 7 in Sect.3.3.

Appendix D: Phosphoros results on SExtractor photometry
Throughout the paper, the results that are presented are those of Phosphoros/hscPipe and Le Phare/SExtractor.This choice makes it difficult to attribute the differences observed in the photo-zs to the codes or to the photometry.To address this problem, we apply here Phosphoros to the SExtractor photometry to create a version of the photo-zs for which only the photo-z codes differ.The configuration used to compute these photo-zs is the same as that presented in Sect.4.1.We use SExtractor 2 ′′ aperture photometry uncorrected for galactic extinction, that we fit with the templates, extinction law, and priors described in Sect.4.1.As a result, PDZs are produced, from which the point estimates (median) are computed.Figure D.1 presents a comparison of the results of SExtractor/Phosphoros (middle panels) with both hscPipe/Phosphoros (top panels) and SExtractor/Le Phare (bottom panels).In all three cases, Ugrizy photo-zs are compared with the spectroscopic sample.In these plots, we can see that the two sets of photo-zs computed from the SExtractor photometry have the same quality when comparing the scatter and the outlier fractions.The scatter from the Phosphoros results is still larger in the bright bins, whereas for the fainter sources the Phosphoros results present a better scatter and outlier fraction than Le Phare photo-zs.This is probably due to the choices in the priors applied by both codes.
These results show that the majority of the differences observed between the photo-zs presented in this work are due to the differences in photometry.Outside of the effect of the priors, there are a few differences between Phosphoros and Le Phare photo-zs, which is expected due to the strong similarities between the two template-fitting codes.All these comparisons thus validate that Phosphoros is able to provide sensible results and is ready to be used by the community.

Appendix E: Description of the released catalogs
The current release (Sect.5) consists of a set of catalogs for the four HSC-CLAUDS regions.Several versions are available depending on the adopted photometry (based on the HSC pipeline (Sect.3.1) or SExtractor run (Sect.3.2)), the photometric code used ( Phosphoros (Sect.4.1) or Le Phare (Sect.4.2)) and the number of photometric bands considered for the photometric redshift runs (Ugrizy and all-bands).Only subregions in the XMM-LSS and E-COSMOS fields are covered with the near-infrared VISTA imaging.
The description of the catalog outputs is given in Tab E.1 for the SExtractor -Le Phare catalogs and in Tab E.2 for the hscPipe -Phosphoros catalogs 9 .
As the footprints of the different surveys are different, the keyword FLAG_FIELD_BINARY allows the user to restrict the data to a specific dataset by combining its 7 boolean elements (described in Tab E.1).A version defined with a single integer is also available (FLAG_FIELD for Le Phare catalogs).To guarantee the selection of clean objects, the flags MASK and ST_TRAIL must be set to 0 (for SExtractor catalogs) or ISOUTSIDEMASK to "True" (for hscPipe catalogs).
All the catalogs include a stellar classification based on χ 2 estimates between the Pickles stellar templates (Pickles 1998) and the galaxy templates.If the sources are compact and better fitted by a stellar template they are flagged as stars (OBJ_TYPE=2 or ISSTAR=True).In addition the Le Phare code also performs a classification using a QSO library (Salvato et al. 2011).If the sources are compact and better fitted by a QSO template, they are flagged as QSO (OBJ_TYPE=1).
For the catalogs based on the Phosphoros code, the PDZs of all sources are provided in separate catalogs for each tract of the 9 the catalogs do not contain all the information produced by the hscPipe output, but the raw photometry file are available through the HSC collaboration distribution portal.https://hsc-release.mtk.A82, page 25 of 25

Fig. 1 .
Fig. 1.Transmission curves of the photometric bands accounting for instrument throughputs and typical atmospheric absorption.All curves are scaled to 1.The horizontal lines indicate the 5σ depths in 2 ′′ aperture measured for each band in the ultradeep regions for the E-COSMOS field.

Fig. 2 .
Fig. 2. Detection images in the four CLAUDS-HSC fields.The footprints of the different observations are over-plotted in different colors: The CLAUDS Deep layer u and u ⋆ bands in blue; the HSC-SSP grizy data in green; the VIDEO Near Infrared data in red (except in the XMM-LSS field where only two J band pointings are available, colored in pink).The footprints of the ultraDeep areas are shown with dashed lines.
Fig. A.1 shows the comparison for 2 arcsec aperture photometry.

Fig. 4 .
Fig. 4. Comparison of hscPipe cmodel magnitudes with SExtractor total magnitudes.The color shades show the density of sources.The solid lines are the median of the distributions and the dashed lines are the ±34% intervals around the medians.These values are computed down to the 5σ depth in each band in the E-COSMOS ultraDeep region (vertical dotted line).
Fig. 7. Comparison of hscPipe and SExtractor colors computed from 2 ′′ aperture photometry.In all the plots, the x-axis shows the SExtractor 2 ′′ aperture magnitudes of the second term in the colors.The solid line is the median of the distribution, and the dashed lines are the ±34% confidence interval around the median.The vertical dotted lines mark the 5σ depth in the reference band (x-axis one).

Fig. 8 .Fig. 9 .
Fig. 8. Color-color density contours of stars in the different fields computed from the hscPipe PSF photometry.The contours enclose 68% of the stellar loci.The crosses show the positions of main sequence stars from the Pickles library in these color-color diagrams.The offset between templates and observed loci indicates the need of proper extinction and zero-point corrections, which are applied in Sect. 4. A82, page 9 of 25

Fig. 13 .Fig. 14 .
Fig. 13.Absolute scaled residuals (D z ) cumulative distributions for the two implementations (HPH: blue lines and SLP: green lines) with the optical (light lines) and optical+NIR (dark lines) catalogs.The horizontal dashed gray lines show the 68% level that all the distributions should cross at D z = 1 (vertical dashed line).

Fig. 17 .
Fig.17.Galaxy redshift distributions (N(z)) in the CLAUDS-HSC regions for three selections, as indicated in the top of each panel.The three panels in the top row show the results based on estimates with all the bands and are thus limited to E-COSMOS and XMM-LSS fields.The two panels in the bottom row show those computed with Ugrizy photometry for the four deep fields.For HPH, the N(z) is either based on the redshift point estimates (median of each PDZ; light red lines) or based on the stack of the individual PDZs (dark red line).For SLP, the N(z) is obtained with the redshift point estimate (blue lines).A comparison is made with COSMOS2020 Farmer/EAzY photo-zs (black lines) in the top panels.The COSMOS2020 curves are weighted to match the CLAUDS-HSC ones.The distributions are constructed in bins of δz = 0.05.

Fig. D. 1 .
Fig. D.1.Comparison of Phosphoros results when computing photo-zs on hscPipe (Top) and SExtractor (Middle) photometry.The results are also compared to those of Le Phare and SExtractor (Bottom).The comparison is made using photo-zs computed using only the Ugrizy bands in all three cases.

Table 1 .
Right ascension and declination coordinates (J2000) of the four deep field centers.

Table 2 .
Summary of the characteristics of the survey in the different fields: filter coverage, 5σ, 2 ′′ depth and area.
Field Depth Complementary bands Area U Depth U Area HSC Depth HSC Area NIR Depth NIR (deg 2 ) u-u *

Table 3 .
Compilation of the zero-point corrections computed by Phosphoros and Le Phare.
Same as Fig.10but for all the HSC-CLAUDS fields combined and split by i-band magnitudes.The comparison is extended to faint magnitude bins, as indicated in each panel, with faint spectroscopic sources coming from the DEEP-23 (m i ≤ 24), XMM-LSS, and E-COSMOS (m i ≤ 26) fields.

Table 4 .
A82, page 16 of 25 G. Desprez et al.: Combining the CLAUDS and HSC-SSP surveys Summary of the results in the different configurations tested in this work depending on i-band cuts.The weighted column indicates if the results on the spectroscopic sample are weighted by the number of photometric sources (see Eqs. ( Notes. Table E.1.Summary of the SExtractor/Le Phare catalog.Mass of best fit with X in (BEST, MED, INF, SUP) SFR_X M ⊙ year −1 SFR of best fit with X in (BEST, MED, INF, SUP) SSFR_X year −1 sSFR of best fit with X in (BEST, MED, INF, SUP) Desprez et al.: Combining the CLAUDS and HSC-SSP surveys Table E.2.Summary of the hscPipe/Phosphoros catalog.Source outside of updated bright stars mask and satellite trail footprint FLAG_FIELD_BINARY Bool[x7] [HSC, u, u * , J, VIRCAM, u * deep , HSC deep ] Only for E-COSMOS and XMM-LSS fields as they both benefit from additional NIR data. a