The bright-star masks for the HSC-SSP survey

TL;DR

This paper details the development and validation of bright-star masks for the HSC-SSP survey, utilizing Gaia and Tycho-2 star catalogs, with automated mask creation and thorough validation procedures.

Contribution

It introduces a new automated method for creating and validating bright-star masks for the HSC-SSP survey, improving upon previous versions and providing publicly available tools and masks.

Findings

Star catalog purity is 99.2% down to G_Gaia < 18.

Masks effectively identify saturated stars in the survey footprint.

Validation confirms masks' accuracy through visual and statistical comparisons.

Abstract

We present the procedure to build and validate the bright-star masks for the Hyper-Suprime-Cam Strategic Subaru Proposal (HSC-SSP) survey. To identify and mask the saturated stars in the full HSC-SSP footprint, we rely on the Gaia and Tycho-2 star catalogues. We first assemble a pure star catalogue down to $G_{\rm Gaia} < 18$ after removing $\sim1.5\%$ of sources that appear extended in the Sloan Digital Sky Survey (SDSS). We perform visual inspection on the early data from the S16A internal release of HSC-SSP, finding that our star catalogue is $99.2\%$ pure down to $G_{\rm Gaia} < 18$. Second, we build the mask regions in an automated way using stacked detected source measurements around bright stars binned per $G_{\rm Gaia}$ magnitude. Finally, we validate those masks from visual inspection and comparison with the literature of galaxy number counts and angular two-point correlation functions. This version (Arcturus) supersedes the previous version (Sirius) used in the S16A internal and DR1 public releases. We publicly release the full masks and tools to flag objects in the entire footprint of the planned HSC-SSP observations at this address: ftp://obsftp.unige.ch/pub/coupon/brightStarMasks/HSC-SSP/.

Figures

Equations

$$G_{\mathrm{Gaia}}=$$

$$r\,\mathrm{[arcsec]}=\left\{\begin{array}[]{ll}708.9\,\times\exp(-G_{\mathrm{Gaia}}/8.41),\,\,\mathrm{Gaia}<9\\ 694.7\,\times\exp(-G_{\mathrm{Gaia}}/4.04),\,\,\mathrm{Gaia}\geq 9\end{array}\right.$$

Full text

1 Introduction

Extragalactic studies rely extensively on multi-wavelength imaging data to measure galaxy photometric redshifts, physical properties and shapes, used in a broad range of probes to explore the cosmology and the formation and evolution of galaxies (see e.g. [Heymans et al. (2012), Magnier et al. (2013), de Jong et al. (2015), Albareti et al. (2016), Lang et al. (2016), Dark Energy Survey Collaboration et al. (2016), Aihara et al. (2017)]). These probes, such as weak gravitational lensing ([Kuijken et al. (2015)] and references therein) and galaxy clustering ([Crocce et al. (2016)] and references therein), require high-precision data reduction and measurement algorithms, whose development and validation are driven by the increasing statistical power of the surveys and the increasing relative importance of the various systematics.

Among the systematics affecting the images are the multiple instrumental artefacts that cause spurious detections (satellite/plane trails, diffraction spikes, fast-moving objects, etc.) and real source contamination or screening (bright stars, ghosts, cross talks, etc.). All of these effects need automated processes to deal with the large amount of data generated in wide angular field surveys. But if many artefacts can be curated during data processing (Bertin & Arnouts, 1996; Ivezić et al., 2004; Magnier et al., 2004; Erben et al., 2009; Schirmer, 2013; Jurić et al., 2015), some others cannot. This is typically the case for the bright stars that saturate on the CCD images. The challenge with bright stars is that the pixels systematically saturate on each individual exposure (so one cannot use co-addition and sigma-clipping techniques to exclude them) and that the brightness information is lost beyond the saturation limit (making it difficult to assess the star properties). Hence, a strategy commonly adopted in the literature is to use external star catalogues and build photometric masks around the bright stars, at the expense of an inevitable loss in area that ranges between 10% and 20% (Coupon et al., 2012; Heymans et al., 2012). For example, the Terapix group was among the first to develop automated tools to mask the bright stars in the Canada-France-Hawaii Telescope Legacy Survey111http://www.cfht.hawaii.edu/Science/CFHTLS/ (CFHTLS), later refined by the VIPERS team (Guzzo et al., 2014; Granett et al., 2012) to improve their spectroscopic target selection. Almost all of the modern wide-field surveys such as the Sloan Digital Sky Survey (SDSS, York et al. (2000)), panSTARRS (Chambers et al., 2016), KiDS (de Jong et al., 2013), or the Dark Energy Survey (DES, Dark Energy Survey Collaboration et al. (2005)), now use similar strategies. In particular the masks produced by the THELI pipeline (Schirmer, 2013; Erben et al., 2009, 2013), used in the CFHTLenS data, were validated by a systematic visual inspection by the CFHTLenS team (Heymans et al., 2012), making sure that the stringent requirements for cosmic shear studies (Kilbinger et al., 2013) were met.

For very sensitive CCD cameras like the recently installed Hyper-Suprime-Cam (HSC, Miyazaki et al. (2012)), mounted on the 8-m Subaru telescope at Mauna Kea in Hawaii, the additional challenge when building a bright-star sample is to find a sufficiently deep and pure star catalogue to cover the entire magnitude range over which stars saturate, but without masking any extended source of interest. Indeed, typical exposures of a few minutes with HSC saturate stars as faint as the 19th (AB) magnitude, and sometimes even fainter when the seeing reaches sharp values under \timeform0.6”.

In this work we describe the procedure to automatically build and validate the bright-star masks for the broad-band filters ($grizY$) of the Hyper-Suprime-Cam Strategic Subaru Proposal (HSC-SSP, HSC collaboration et al. (in prep)) survey. HSC-SSP is an extragalactic survey whose aim is to observe $1\,400$ deg2 down to $grizY\sim 25-26$ (as well as a number of smaller and deeper fields). To do so, we use the data from the Gaia DR1 (Gaia Collaboration et al., 2016b), Tycho-2 (Høg et al., 2000; Pickles & Depagne, 2010), and SDSS DR13 (Albareti et al., 2016) catalogues, to build a complete and pure star catalogue down to the typical HSC saturation limit. We first describe the problem in Section 2, we then describe how we gather a sample of bright stars in Section 3 and how we build the masks in Section 4. In a last Section 5 we estimate the purity and completeness of the bright stars masks, quantify the improvement on a number of observables, and describe their limitations. We conclude in Section 6.

This work supersedes an early version (hereafter the Sirius version), described in Mandelbaum et al. (in prep.), of the bright-star masks that were used in the S15B and S16A internal data releases, as well as in the first public data release DR1 (Aihara et al., 2017). This early version contained about 9% of bright galaxies (making it problematic for a number of science cases), and was over conservative in the choice of mask radius for stars brighter than the 5th magnitude, causing an unnecessary loss in survey area. The version of the bright-star masks described in this work is called the Arcturus version, and will be implemented in subsequent HSC pipeline versions and data releases. The full masks and tools to flag objects in the HSC-SSP footprint can be downloaded at ftp://obsftp.unige.ch/pub/coupon/brightStarMasks/HSC-SSP/.

2 Bright stars in the HSC-SSP survey

2.1 The HSC-SSP survey

Here we briefly describe the HSC-SSP survey. For more details we refer the reader to the camera (Miyazaki et al., 2012, in prep.), the survey overview222http://hsc.mtk.nao.ac.jp/ssp/ (HSC collaboration et al., in prep), the HSC software pipeline (Bosch et al., in prep.), the DR1333https://hsc-release.mtk.nao.ac.jp/doc/ (Aihara et al., 2017), and the photometric redshifts (Tanaka et al., 2017) papers. The HSC-SSP data are calibrated with the PanSTARRS1 data. We refer the reader to Tonry et al. (2012), Schlafly et al. (2012) and Magnier et al. (2013) for additional details on the PanSTARRS1 data.

The HSC-SSP project is a 300-night deep imaging survey which started in March 2014, carried out at the Subaru telescope with the 1.5 deg-diameter HSC camera. The goal of the survey is to explore the properties of the dark universe (dark matter and dark energy) and the formation and evolution of the galaxies from the local Universe until the epoch of reionisation, using a number of different probes, including weak and strong gravitational lensing, galaxy clustering and galaxy abundances.

The survey is split into three layers:

• •

a Wide layer ($r<26$444$5\sigma$ point-source limiting magnitudes. Note that all magnitudes quoted in this work are in the AB system.) over $1\,400$ deg2, composed of three main fields (Spring, Fall, and Northern) and one pointing for calibration (Aegis),

• •

a Deep layer ($r<27$) over 28 deg2 composed of four fields (XMM-LSS, Elais-N1, E-COSMOS and DEEP2-3),

• •

and an UltraDeep layer ($r<28$) over 4 deg2, composed of two fields (SXDS/UDS and COSMOS).

The three layers are observed in the five optical and near infrared (NIR) broad-band filters $g$, $r$, $i$, $z$ and $Y$. In addition, the Deep layer is observed in the three narrow-band filters NB387, NB816, and NB921, and the UltraDeep layer is observed in the three narrow-band filters NB816, N921, and NB101.

2.2 Why masking the bright stars?

Our goal is to limit the impact of saturated stars on nearby sources, to prevent the use of spurious sources in the science analyses, and to record the loss in area due to the occultation of the bright stars for extragalactic studies. The adopted strategy is to flag all sources near a bright star and propagate the information all the way from the data processing to the database (i.e. we keep all the sources but we flag those affected), so that one can filter a sample of interest during data retrieval from the database.

It is also important to record the fraction of the survey that is occulted by the bright-star masks, so in parallel we produce a catalogue of points randomly placed in the survey area, flagged in the same way as the sources. In short, the random points are drawn with a density of 100 points per square arcmin, for each patch and each filter. Then, the survey information (field geometry, number of co-added exposures, PSF, pixel variance, etc.) is recorded at each random-point position using a python module running within the HSC pipeline environment555The source code can be accessed at https://github.com/jcoupon/maskUtils and the random-point sample is described at https://hsc-release.mtk.nao.ac.jp/doc/index.php/random-points-for-dr1/.. Finally, we add a flag indicating whether a point is inside or outside the bright-star masks (more details on the flagging procedure are given in Appendix B.7 ). The masked fraction is computed from the ratio between the random points inside the masks and the total number of points. With a density of 100 points per square arcmin and assuming a mask coverage of $\sim 10\%$, the precision on the masked area is $0.5\%$ on the scale of a degree.

2.3 The impact of saturated stars on HSC images

The main artefacts caused by the saturated bright stars can be divided into optical effects (diffraction pattern, optical reflections and ghosts) and electronics effects (electrons leaking to nearby pixels, NIR photons scattering, brighter-fatter effect). On the image, these effects result in (1) Point Spread Function (PSF)-shaped luminous patterns, (2) extended luminous haloes and ghosts, (3) linear “bleed trails”, and (4) spikes in the $Y$-band images in the row direction in addition to bleed trails in the column direction. This is the result of scattering off the channel stops at the back of the detector. The sensors are too opaque to bluer photons for this effect to appear in other bands. An example of a typical saturated bright star on a HSC image is shown in the top panel of Figure 1.

The size of the PSF diffraction pattern and the bleed trails depend on the seeing and the position on the focal plane. The radial extent of the luminous PSF pattern depends on the exposure time and the brightness of the star. The brighter-fatter effect (Antilogus et al., 2014) also mildly increases the extent of the diffraction pattern as a function of star brightness. Optical reflections in the instrument are responsible for bright images (ghosts) that spread in various places on the image, also depending on the position on the focal plane. In particular, the ghosts include a circular extended halo generated by the incident light which is reflected at the CCD surface and again reflected by the filter surface. The size of this extended halo is determined solely by the distance between the two reflecting surfaces, being in theory independent on the brightness of the ghost source star, but appear only for very bright stars.

In the bottom panel of Figure 1, we show an example of the radial luminosity profile created by the instrumental and the atmospheric PSF, as well as a number of optical reflections. The blue line is the Kolmogorov profile which describes the atmospheric turbulence with a full width half maximum (FWHM) seeing size of \timeform0.5”. The combination of $\beta=7$ and $2$ Moffat functions (plotted as blue dotted and dash-dotted lines, respectively) is used to describe the Kolmogorov profile (see Racine (1996)). The instrumental PSF, which is simplified as a Gaussian with FWHM=\timeform0.357” (Miyazaki et al., in prep.) and plotted as the green line, is smaller compared to atmospheric seeing profile and can be neglected (it only affects the image on scales on the order of \timeform0.1”).

When the contrast reaches $10^{-7}$, the circular ghosts emerge. The size and contrast are calculated for different combinations of reflecting plane surfaces (CCD, and incident surface, S1) and exit surface (S2) of the window and the filter, based on the optics model assuming that the reflectance of each surface is known. The profiles are represented as red dotted lines and the significant ones are labeled. The sum of all components is the black solid line. Note that not all of the low-contrast components – in particular those with no theoretical models (including the aureole, Racine (1996)) – are included in this figure (it must be obtained from the real images). Also, the extended halo is not exactly circular except in the central part of the camera, due to the vignetting of the wide-field corrector cutting the pupil image of the primary mirror.

2.4 Saturation limits in the HSC-SSP survey

In this work we limit ourselves to the sample of bright stars that saturate on HSC-SSP images (the information on non-saturated stars – such as brightness and PSF-size – can be easily gathered directly from the data through the HSC pipeline estimates). Here, we evaluate the typical saturation limits in the HSC-SSP survey.

The saturation limit on a given co-added image mainly depends on the single-visit exposure time, which is different for each survey layer as shown in Table 2.4,

and the PSF size (the sharper the PSF, the fainter the saturation limit). The PSF at a given position is computed by the HSC pipeline from fitting the (non-saturated) stars with a smooth spatially varying PSF model, per single exposure, and co-added using a sigma-weighted mean (see Bosch et al. (in prep.)).

To evaluate the typical PSF sizes on co-added images in HSC-SSP, we record, for many random positions in the S16A release footprint, the PSF sizes from the second moment of the best-fit PSF model, multiplied by $2\sqrt{2}$ to convert into the commonly used FWHM estimate, under the assumption that the PSF luminosity profile is gaussian666Although it is not really gaussian, this assumption is good enough for the relative comparison between two positions in the survey.. The top row of Figure 2 shows the PSF size (FWHM) distribution in each survey layer and each of the $grizY$ filters.

The bottom row of Figure 2 shows the median saturation magnitudes in each survey layer and each filter. Saturation magnitude limits are computed from selecting all bright point sources detected by the pipeline, with a successfully measured magnitude (PSF estimate), imposing that the central pixel is flagged as saturated and that the source is not deblended. The typical magnitude limits range from $17$ to $20$. The figure shows the expected trend that the saturation limit becomes fainter as the PSF gets smaller. We also confirm that the longer exposure times in the Ultra Deep layer result in fainter saturation limits as compared to the Wide layer, although the difference is not very large. The changing limit observed in the different filters is linked to the varying filter transparency and chromatic effects of the PSF.

The saturation limits in best-PSF data are summarised in Table 2.4. These represent the faintest possible magnitudes at which stars saturate in the HSC-SSP data. We observe that, depending on the filter, the faintest magnitude occurs in a different survey layer. Several effects are competing with each other: on the one hand, we expect that the longer single-exposure times in the Deep and Ultra-Deep layers will lead to fainter saturation limits (this is the case for the $z$ and $Y$ bands). On the other hand, the shorter exposure times in the Wide layer are compensated by usually better PSFs (by construction, e.g. for the $i$ band, or by chance), so that the saturation limits are equivalent or fainter than in the deeper layers (note for example the $i$-band magnitude of $20.7$ that occurs in very good PSF data, $<\timeform{0.5^{\prime\prime}}$).

Eventually, the goal is to build a bright-star sample which goes as faint as the best-PSF saturation limits, so that all saturated stars are masked in the HSC-SSP footprint. We will see that the Gaia sample is well suited for this purpose.

3 The bright-star sample

The reference bright-star777Here, other bright point sources, such as Quasi Stellar Objects (QSOs), are included in the term bright stars. sample should meet several requirements to be used for masking:

• •

it should be pure and complete over the magnitude range where stars saturate on HSC images,

• •

it must cover the entire area of the HSC-SSP survey,

• •

and it must not rely on HSC-SSP data (since the bright star masks are made before data reduction and because we do not trust the brightness estimates where pixels saturate).

We build the star sample from the Gaia Data Release 1 (DR1), completed with the Tycho-2 star sample to address a number of issues associated with the first Gaia release. We use the $G_{\rm Gaia}$ magnitude – a broad-band filter covering the range 4000–10000 Å – as the primary characterisation of star brightness, and we use the SDSS to remove the extended sources wrongly identified as stars by Gaia in the DR1. We detail below the steps we follow to build the sample. A quick summary is given in Section 3.5.

3.1 Gaia

Gaia (Gaia Collaboration et al., 2016a) is an ESA space mission launched in 2013 whose goal is to observe one billion stars (1% of the total Milky Way stars), and to provide accurate positions, radial velocities, and spectrophotometry in optical888The spectrophotometry will be limited to relatively bright stars, $G_{\rm Gaia}<15$.. The first Gaia data release (Gaia Collaboration et al., 2016b) was publicly delivered in September 2016. It includes position measurements and $G_{\rm Gaia}$-band photometry over the entire sky.

The Gaia star sample is ideally suited for our purpose. With a survey covering the entire sky at a limiting magnitude of $G_{\rm Gaia}<20.7$ and with optical spectrophotometry, it perfectly meets our requirements. However, the Gaia DR1 comes with a number of limitations:

• •

only the $G_{\rm Gaia}$ photometry is provided,

• •

the sample is incomplete at bright magnitude,

• •

the limiting magnitude is ill-defined and varies across the sky,

• •

some large areas are missing on the ecliptic due to quality filtering and not-yet fully scanned regions,

• •

and some galaxies are included in the sample towards faint magnitudes.

3.2 Tycho-2

To complete the Gaia sample at bright magnitude and in under-dense (or empty) areas, we use the Tycho-2 star sample.

The Tycho-2 star sample (Høg et al., 2000) is a complete all-sky catalogue of the 2.5 million brightest stars to a magnitude $V<11.5$, with accurate positions and proper motions, observed by the ESA Hypparcos satellite. It provides two broad-band photometry in $B_{T}$ and $V_{T}$ filters. Here we use the catalogue released by Pickles & Depagne (2010), which provides emulated SDSS-like $ugriz$ photometry for all Tycho-2 stars. To be consistent with the Gaia stars, we transform the emulated SDSS magnitudes for each Tycho-2 star into a $G_{\rm Gaia}$ magnitude, according to the formula given in Jordi et al. (2010):

[TABLE]

We check the accuracy of this expression in Appendix A.

Since a star’s brightness varies across the spectrum, it would be, in principle, required to set the size of the masks in each band, separately, using the color information when available. However, we take the simpler (and limited by the Gaia DR1 release) approach to characterise the star brightness in all HSC bands using the unique $G_{\rm Gaia}$ magnitude. We come back to this issue in Section 5.2.

3.3 Gathering a complete star sample

In practice, we first obtain the DR1 Gaia and Tycho-2 catalogues in an area enclosing the full planned HSC-SSP footprint, extended to one degree around the borders to account for bright stars outside the survey footprint that can potentially affect the sources inside. The full HSC-SSP footprint will be eventually composed of three Wide fields: Spring, Fall, and Northern, as well as one pointing in the AEGIS field used for photometric redshift calibration. All the Deep and Ultra Deep fields are within the Wide footprint, except the Elais-N1 Deep field, in which we add to our star sample the corresponding Gaia and Tycho-2 stars. We show in Figure 3 the sky density of the bright-star catalogue in the planned HSC-SSP footprint. Note the higher star density near the Galactic plane at R.A. $\sim 22$h and $\sim 9$h in the Fall and Spring fields, respectively.

Then, we combine the Tycho-2 and Gaia catalogues. We select all sources from the Gaia DR1 release, filtering out those flagged as duplicated_source. When a source is found in both Gaia and Tycho-2 catalogues, we keep the Gaia position and magnitude. For the Tycho-2 stars that are not in the Gaia sample, we emulate the $G_{\rm gaia}$ magnitude following Equation 3.2.

In Figure 4 we show the magnitude distribution of Gaia and Tycho-2 sources. The Gaia completeness is especially affected at bright magnitude due to the selection imposed in DR1, but also at all magnitudes due to some spatial incompleteness, as illustrated in Figure 3. The addition of the Tycho-2 sample ensure at least a complete star sample at $G_{\rm gaia}\lesssim 12$.

3.4 Ensuring a pure star sample

To remove the galaxies incorrectly identified as stars by Gaia, we make use of the photometric information provided by the SDSS survey (York et al., 2000; Albareti et al., 2016). To do so, we match the Gaia+Tycho-2 sample previously assembled to the SDSS photometric sources downloaded from the DR13 release database999http://www.sdss.org/dr13/. We define a source as extended in the SDSS using the photometric information provided in the three $g,r$ and $i$ bands with the highest signal-to-noise ratio as folows:

clean == 1 AND (obj.flags_g & dbo.fPhotoFlags(’SATUR_CENTER’))==0 AND (obj.flags_r & dbo.fPhotoFlags(’SATUR_CENTER’))==0 AND (obj.flags_i & dbo.fPhotoFlags(’SATUR_CENTER’))==0 AND ( (0.1 < psfMag_g-cModelMag_g AND psfMag_g-cModelMag_g < 20.0) OR (0.1 < psfMag_r-cModelMag_r AND psfMag_r-cModelMag_r < 20.0) OR (0.1 < psfMag_i-cModelMag_i AND psfMag_i-cModelMag_i < 20.0) )

The first three conditions impose that the center of the source is not flagged as saturated in any of the $gri$ bands. We do not expect galaxies – even bright – to saturate in the SDSS. The next three conditions are similar to the definition adopted by the SDSS collaboration to define extended sources, but are more conservative to ensure a 100% complete sample of extended sources, the goal being to remove all galaxies from the bight-star sample. In short, this condition means that all sources with a luminosity profile 10% brighter than the PSF luminosity profile in any of the $gri$ bands are flagged as extended.

To validate our approach we also match the star catalogue to the S16A HSC-SSP data. The HSC images are much deeper and with a better resolution than SDSS, allowing to do a more accurate extended/point-source separation. As mentioned above, our objective is to build an HSC-SSP-agnostic bright-star sample, so the comparison with HSC-SSP serves only as validating our approach.

We flag extended sources in HSC-SSP using exclusively the $i$-band photometry in the Wide, as the $i$-band observations are always done in best-seeing conditions. We define an extended source as:

parent_id == 0 AND NOT iflags_pixel_saturated_center AND 0.01 < imag_psf-icmodel_mag AND imag_psf-icmodel_mag < 20.0

This selection closely follows that of the SDSS, with the exception that it is only done in the $i$-band and that the luminosity profile should only change by 1% to be flagged as extended. The first condition ensures that the source is not deblended. Overall such a definition is also extremely conservative to ensure that all galaxies are properly identified and removed from the star sample.

The fraction and magnitude distribution of extended sources in Gaia, as a function of $G_{\rm gaia}$ magnitude, are shown in Figure 5.

Extended sources start contaminating the star sample in the regime fainter than $G_{\rm gaia}\sim 12$, but remains below 3–5%. The fraction starts becoming significantly larger beyond $G_{\rm gaia}\sim 18$, reaching at least 10% at the 19th magnitude. The fraction increases significantly beyond 19 if using the SDSS, probably due to the highly conservative cut in extendedness. We conclude that building a pure star sample (at least with the Gaia DR1 data) requires that we limit the Gaia star sample at $G_{\rm gaia}<18$, which is slightly brighter than the ideal limit for HSC-SSP (as shown in Section 2.4).

In the regime $14<G_{\rm gaia}<18$, we remove all sources ($1.5\%$) that are flagged as extended in the SDSS according to the criterium given above. At $G_{\rm gaia}<14$, visual inspection confirms that the “bump” seen in the SDSS source distribution is caused by stars with a significantly extended diffraction pattern.

The purity of our star sample is further investigated in Section 5.1.

3.5 Summary

We summarise here the construction of the bright-star reference sample to be used for masking:

• •

our starting catalogue is the Gaia DR1 all-sky survey sample,

• •

we restrict the sample to the footprint of the planned HSC-SSP survey, increased by one degree around the borders,

• •

we complete the missing bright stars with the Tycho-2 star sample ($V<11.5$) and emulate the $G_{\rm Gaia}$ magnitude from the multi-band information,

• •

the sample is cut at $G_{\rm Gaia}<18$ beyond which both SDSS and HSC-SSP images show increased extended source contamination ($>10\%$),

• •

we remove all sources ($1.5\%$) flagged as extended in the SDSS between $14<G_{\rm Gaia}<18$.

The final sample contains $1\,812\,106$ stars.

4 Building the masks

We describe in this section the adopted metric to empirically evaluate the impact of bright stars on the sources, and how we build the masks. Our goal is to account for the PSF-shaped luminous patterns, the extended luminous haloes, the $Y$-band vertical spikes and the linear bleed trails. The removal of the large-scale reflection ghosts, however, will be implemented in a future version of the HSC pipeline.

4.1 Measuring the impact on detected sources

We show in Figures 6 and 7 the 3$\sigma$-clipped stacked HSC images of the bright stars, stacked per bin of magnitude, and shown for each of the $g$, $r$, $i$, $z$ and $Y$ filters.

We use all bright stars at the brightest magnitudes ($G_{\rm Gaia}=3.5,4.5$ and 5.5), whereas for fainter stars, to ease the computation, we only select stars in a narrow magnitude range around the magnitudes $G_{\rm Gaia}=6.5,7.5,9.5,11.5,13.5,15.5$ and 17.5. Our binning scheme is summarised in Table 4.1.

We see that for the brightest stars ($G_{\rm Gaia}<9$), the main effect is due to the circular bright halos. For stars fainter than $9$, the star halos are too faint to cause significant issues, and the dominant effect is the PSF-shaped pattern.

First, we assume that the impact of bright stars on nearby sources is isotropic and, to measure the extent of such impact, we match the bright-star sample, assembled in Section 3 and binned in $G_{\rm Gaia}$ magnitudes, with the HSC-SSP sources from the S16A release. We then measure the mean density of sources as a function of distance from the star position, normalised to the density measured in an outer annulus.

The HSC pipeline proceeds in two steps to detect the sources (see Bosch et al. (in prep.)). First it detects all aggregated pixels (footprints) above the background, called the “parents”. Second, the footprints featuring several significant peaks are deblended into “children”, whereas the sources that are not deblended are called “single”. We measure the mean density of sources around the bright stars for two separate samples: (1) the parent sources, and (2) the single and children sources, referred to as “primary” sources. The drop in density of the parent sample is an indicator of the typical radius below which sources are being hidden by the PSF-shaped luminous pattern of the nearby bright star. The variation in density of the primary sample indicates where the HSC pipeline deblends artefacts (instead of real sources) caused by the bright star, and is an indicator of potential problems affecting the detection of real sources and their photometric measurements.

We show in Figure 8 the normalised detected sources density for the parent (in blue) and primary sample (in green), averaged over the bright stars, binned per $G_{\mathrm{Gaia}}$ magnitude.

For the brightest stars ($G_{\mathrm{Gaia}}<9$), the impact on sources occurs at larger radius, due to the concentric luminous halos around the star which artificially increase the density of primary sources compared to the mean density. For fainter stars ($G_{\mathrm{Gaia}}>9$), the impact occurs at lower radius, and is characterised by a rapid drop of parent sources towards the star position.

4.2 Setting up the mask sizes

To quantify the isotropic extent of the masks, we use both the parent-source and the primary-source indicators. First, we fit an error function to the parent-source density distribution, where the free parameters are the position and width of the slope. The best-fit error function is shown as a blue dashed line on Figure 8. We then measure the radius where the best-fit completeness reaches 95%, and the radius where the primary-source density is 20% above the mean source density. We finally set the mask radius to the most extended indicator, to guarantee that both the primary and parent source samples are unaffected by the bright star. The final mask radius is marked by a black arrow on Figure 8.

Figure 9 shows the mask radius as a function of star magnitude. The measured radii are shown as blue dots. The former relation (Sirius version) is also shown for comparison as the green dashed line: one can see the diverging radius at magnitudes brighter than $\sim 5$ which explains the very large masks seen in S15B, S16A and DR1 releases.

To build a circular mask for each bright star as a function of magnitude, we model the magnitude dependence of the radius by fitting two exponential functions below and above $G_{\mathrm{Gaia}}=9$, finding the following relations:

[TABLE]

The best-fit relation is shown as the blue solid line on Figure 9, and the resulting circular masks are shown as white circles in Figures 6 and 7.

4.3 Accounting for the anisotropic effects

Next, we account for the saturated bleed trails that appear as horizontal black lines on the images in all filters. We mask these with horizontal rectangles centred on the star position, and whose length and width are $1.5$ and $0.15$ times the radius of the circular mask, respectively. Similarly, to mask out the $Y$-band vertical spikes, we use vertical rectangles whose length and width are also $1.5$ and $0.15$ times the radius of the circular mask, respectively. The choice for the rectangle sizes are based on visual inspection. Although the vertical spikes are observed only on the $Y$-band images, we replicate the vertical rectangles for the other filter masks, for simplicity.

These two additional mask components are used for stars fainter than $9$ and shown as the horizontal and vertical rectangles in Figure 7. Around stars brighter than $9$, the extent of the luminous halo goes beyond the typical extent of the bleed trails and vertical spikes, so we do not add any rectangle.

4.4 The reflection ghosts

Finally, we note that the reflection ghosts for very bright stars (one can see an example of a ghost in the left part of the top $g$-band image in Figure 6) will be directly removed from the images in a future version of the HSC pipeline.

In brief, ghosts are caused by internal reflections inside the optics system and their shape strongly depends on the position of the star on the focal plane. Since the surface brightness of the circular ghost is almost uniform within the circle, it can be subtracted from the image once the boundary of the circular ghost is detected. By fitting the boundary with two ellipses, one corresponding to the pupil image and the other to the vignetting, Yagi, Utsumi & Komiyama (in press) successfully fitted the circular ghosts featuring a high signal-to-noise ratio. Then, they subtracted them from the Suprime-Cam images. The method will be implemented in a future version of the HSC pipeline as an extension of the ghost subtraction procedure which is already dealing with the cometary-shape and arc-shape ghosts (see Figure 18 of Aihara et al. (2017)). If a perfect subtraction of the circular ghosts is carried out, the atmospheric component of the PSF remains. Hence, the mask size for the brightest stars would follow the same relation as for fainter stars ($G_{\rm gaia}>9$, as seen in Figure 9), resulting in a reduction in size by a factor of $\sim$3.

4.5 Handling the masks in the HSC pipeline

In practice the masks are recorded in DS9 region (“.reg”) format. A file is produced for each individual patch (the smaller subdivision of the HSC-SSP survey). The mask file is read by the HSC pipeline, during data processing. Then the footprint of the mask is converted into pixels and the corresponding bit value is set to the general mask plane (see Bosch et al. (in prep.), for more details). When extracting the source properties, the pipeline records if the centre or any pixel of a detected source overlaps with the bright-star mask (but the masks are not used during data processing).

5 Masks validation

In this section we assess the quality of the bright-star catalogue and the masks using the data from the S16A release. We perform visual inspection on the images, we address the potential issues with seeing and color dependence of the bright-star artefacts, and we run a number of validation tests. We end this section with an estimate of the masked fraction.

5.1 The bright-star sample purity

For many science cases, the purity of the star sample is essential to avoid masking bright galaxies. To test the robustness of our method in selecting a pure star sample (as described in Section 3.4), we match the bright-star catalogue with the HSC-SSP sources from the S16A release, assuming that extended galaxies do not saturate on the images. With image quality superior to SDSS, the HSC-SSP images are more robust to identify extended sources towards faint magnitudes.

As seen from the bottom panel of Figure 5, the SDSS-based criterium tends to flag as extended a higher number of sources than the HSC-SSP-based one, but it does not necessarily means that all truly extended sources are safely removed. Therefore, we select the sources (from the parent sources samples) that are flagged as extended in HSC-SSP but not in SDSS. Out of a total of $325\,028$ sources in the S16A HSC-SSP footprint, extended sources in HSC-SSP amount to $2\,609$ ($0.8\%$). We visually inspect all of the $2\,609$ objects and find that most of them are:

• •

binary stars resolved in HSC-SSP (unlike in SDSS),

• •

stars randomly aligned with other objects (we remind that our HSC criterium takes only non-deblended sources),

• •

or stars near artefacts.

Those represent a nuisance for nearby galaxies and do not need to be removed from the star sample. We show some examples in the top row of Figure 10.

The galaxies that we do not want to mask are shown in the bottom row of Figure 10. We found a total of 13 cases ($0.004\%$). Among those, a few are most probably galaxies hosting a bright quasar. The need for masking those sources depends on the science case: the photometry (hence the photometric redshift and physical properties) is likely affected, but identifying galaxies hosting a quasar is obviously a very important goal. Although the number of such cases is small (most are safely identified as extended in the SDSS), it may be a concern for studies aiming at finding quasars in detected galaxies. One must also keep in mind that bright quasars with no detected host are most likely in the bright-star sample and therefore masked, so that quasar studies must use the masks with caution.

5.2 Filter dependence of star brightness

For each individual star, we set identical size and shape for all five HSC filters, based on its Gaia $G$-band brightness, which covers the wide optical range from $4\,000$ to $10\,000$ Å. Although this range covers the HSC bands from $g$ to $z$ (and slightly into the $Y$ band), one may ask the question whether a star significantly brighter in one band could impact the nearby sources beyond the mask limits, that are computed from the source density averaged over the five bands.

Fortunately, the HSC pipeline first detects sources in all bands independently, hence if some artefacts artificially increase the number of detected/deblended sources in one band and not the others, these fake sources will appear in the catalogue, and hence will be caught by our estimate based on the increase in source density. So we believe that our approach properly accounts for such variations in brightness from one filter to another. We illustrate this in Figure 11, which shows the bright-star masks on top of the five-band HSC images. This field is located in the currently most star-crowded region in S16A (tract 9799 and patch 6,5), located in the “GAMA09H” area, close to the Galactic plane. We see that the masks safely enclose the bright stars in all filters.

However, some stars may be significantly brighter in one filter but may not pass the Gaia magnitude cut ($G_{\mathrm{Gaia}}<18$), especially NIR-bright stars. This means that some saturated stars in the $Y$ band may not be properly masked. This caveat must be kept in mind in the studies extensively relying on the $Y$-band photometry (e.g. when selecting high redshift dropout galaxies).

5.3 PSF-size dependence

Similarly, bad seeing (which increases the size of the PSF luminous pattern) or sharp seeing (which will saturate fainter stars and enlarge the seeing pattern due to the brighter-fatter effect) may require enlarging the required size of the masks. To check the impact of seeing on the masks, we recompute our estimates (see Section 4.2) for the $i$-band $10\%$ best- and $10\%$ worst-PSF data.

Table 5.3 shows the mean PSF size values for the full star sample, the $10\%$ best- and $10\%$ worst-PSF samples and the results are shown in Figure 9.

The red triangles represent the $10\%$ best-PSF results, whereas the orange squares represent the $10\%$ worst-PSF results. For both cases, we observe no significant change over the full magnitude range.

Overall, we find that the PSF-size dependence has little impact on the measured radius. This is partly due to the fact that all five bands are used in the source detection (as discussed in the previous section) and the effect of the PSF is averaged over the five bands. We note however that these tests are made with the $i$ band which has a constraint on good seeing for shape measurement.

5.4 Number counts

We now perform a number of comparisons with the CFHTLenS dataset (Heymans et al., 2012). The CFHTLenS project is a cosmic-shear driven data reduction of the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS), a medium-deep large-scale cosmological survey covering about $150$ deg2 observed in five optical ($ugriz$) bands. Results from the CFHTLenS are now a reference in a number of extragalactic and cosmological topics (see e.g. Kilbinger et al. (2013)). One remarkable feature of the CFHTLenS is that the masking was carefully done and refined by hand by the team to meet the stringent requirements on shape measurement for cosmic shear. It is therefore the ideal dataset to assess the quality of our masks. To do a fair comparison between the two surveys, we select all sources brighter than $i<22.5$ in both the CFHTLenS and HSC-SSP surveys, in a common area covering about 20 deg2, located in the CFHTLenS “W4” field. Both $i$-band magnitudes are total magnitudes (MAG_AUTO for CFHTLenS and cmodel for HSC-SSP), and the $i$-band filter response is similar for both.

We begin the comparison with the galaxy number counts in both catalogues. The star/galaxy separation is done using color and morphology for CFHTLS (Coupon et al., 2015), whereas it is done using the extendedness estimate in HSC-SSP (Bosch et al., in prep.), based on PSF vs cmodel flux comparison. Results are shown in Figure 12.

The top panel shows the old mask version (Sirius), whereas the bottom panel shows the new mask version (Arcturus) developed in this work. In the former case, many bright galaxies are missing due to a number of extended sources that were wrongly included in the bright-star sample (see Appendix C). In the latter case, the HSC-SSP number counts are in remarkable agreement with the CFHTLenS in the magnitude range $16<i<22.5$. We note however a few missing sources brighter than $i=16$, which are due to cmodel measurement failures for large radius galaxies (Huang et al., in prep.).

5.5 The two-point correlation function

Next, we compare the galaxy angular two-point correlation function (TPCF). The TPCF is an observable sensitive to the photometric artefacts, as it captures the local changes in density. The spurious detections in the luminous halos around very bright stars or the drop in source density observed near less-bright stars may cause significant biases in the TPCF.

We measure the TCPFs in six $i$-band magnitude bins in the range $17.5<i<22.5$, filtering sources using the new (Arcturus) and the old masks (Sirius), displayed in Figure 13 as blue and green dots, respectively.

We compare these with the CFHTLenS TPCF, shown as a blue solid line. The agreement is remarkably good in the case where the new masks are applied (in blue). This result shows the robustness of the masks, both in its ability to remove problematic data around bright stars, but without wrongly removing bright galaxies.

In the case of the Sirius masks (shown in green), we observe a drop in density at small scales, especially at bright magnitude, due to the removal of bright galaxies. We also show the TPCF measured inside the bright-star masks, to illustrate the effects of the nearby bright stars, that lead to an estimate with significant systematics.

5.6 CAMIRA clusters

Finally, we check the cluster member estimates in the cluster catalogue generated by the CAMIRA cluster finder (Oguri et al., 2017). Figure 14 shows the relative fraction of cluster members after and before applying the masks. The top and bottom panels show the Arcturus and Sirius versions, respectively.

After applying the Sirius version mask, low-redshift clusters feature a significant lower number of cluster members at low redshift, due to the masking of bright galaxies. This trend disappears with the (new) Arcturus version, showing that the latter is well suited for cluster detection and cluster member identification.

5.7 The masked fraction

Finally, we compute the masked fraction as a function of position on the sky. We generate a sample of random points and compute the fraction of points that fall inside the bright-star masks, compared to the total number of random points in a given area. The result for the full planned HSC-SSP footprint is shown in Figure 15.

The fraction ranges from a few percents to nearly 18% in the most star-crowded regions.

6 Conclusions

We gathered a pure star sample used to build photometric masks around saturated stars in the full planned HSC-SSP survey footprint. The masks were validated on the S16A internal release data. We performed visual inspection on the HSC images and we measured a number of observables that we compared with the literature.

Given the high sensitivity of the HSC camera, the typical good-seeing conditions at the telescope location (e.g., a median value of \timeform0.6” for the $i$ band in the Wide layer) and the exposure times adopted in the HSC-SSP survey, we found that the bright-star saturation limit magnitudes, in best-seeing conditions, are as faint as $19.2,20.0,20.7,18.5,18.2$, for the $g,r,i,z,Y$ bands, respectively, and depend on the PSF. The effect of saturated stars translates into a PSF-shaped luminous pattern (whose extent increases as a function of star brightness due to the increasing contrast), extended luminous haloes and ghosts, linear horizontal bleed trails, and vertical spikes (in the $Y$ band).

We built a sample of bright stars from the Gaia DR1 release, completed at bright magnitude ($<12$) with the Tycho-2 star catalogue. The $G_{\mathrm{Gaia}}$ magnitude, measured in a broad-band filter that extends from the HSC $g$ band to $z/Y$ band, is our primary star brightness estimate. For the Tycho-2 stars with no Gaia observations, we computed the $G_{\mathrm{Gaia}}$ using an analytical colour transformation applied to SDSS filter-emulated magnitudes from the literature. We matched the sample with the SDSS source catalogue and measured the fraction of sources that appear extended in the SDSS. For Gaia sources brighter than $G_{\mathrm{Gaia}}=18$, we found that about $1.5\%$ appear as extended in the SDSS, which we removed from the star sample. Using the HSC-SSP S16A data we confirmed that our cleaning strategy leads to a star sample which is $99.2\%$ pure and, from visual inspection, the remaining sources that appear extended in HSC are blended binary stars, stars near artefacts, and overlapped with extended sources. Beyond $G_{\mathrm{Gaia}}>18$ where we do not use Gaia stars, the fraction of extended sources lies between $10$ and $20\%$, depending on the HSC-SSP or SDSS criterium. Due to our conservative approach to identify and withdraw from the star catalogue all sources that appear as extended in the SDSS, and the incomplete Gaia observations between $14<G_{\mathrm{Gaia}}<18$ in a few small areas in the HSC-SSP footprint, our bright-star sample suffer from some incompleteness towards faint magnitudes.

To set up the mask sizes as a function of star brightness, we first matched the bright-star positions with all the sources detected by the HSC pipeline. Second, we computed the mean source density as a function of distance from the star position, and measured two kinds of estimates: the parent source density and the primary source density. The former indicator is used to capture the drop in density caused by the PSF-shaped luminous pattern that screens the foreground sources. The latter indicator reveals the typical extent of the luminous circular halo. We set the circular mask radius as the largest value between the two. We found that for stars brighter than $G_{\mathrm{Gaia}}<9$, nearby sources are primarily affected by the luminous halo, whereas for fainter stars they are primarily affected by the screening effect from the PSF-shaped luminous pattern. To automate the process of building individual masks, we parametrised the required mask radius as a function of Gaia magnitude, using two exponential functions below and above $G_{\mathrm{Gaia}}=9$. To mask the additional anisotropic effects caused by the bleed trails and the $Y$-band vertical spikes, we added two horizontal and vertical rectangles scaled to the size of the circular mask.

After applying the masks to the latest data release from HSC-SSP (S16A), we performed a number of validation checks to make sure of (1) the purity of the bright-star sample and (2) the robustness of the masks. We found no indication of a significant loss of bright galaxies due to the masks and we measured good agreements, within the error bars, between HSC-SSP (after masking) and CFHTLenS, for the galaxy number counts and the two-point correlation functions.

These procedures are well suited for wide-field surveys in which building (and validating) all individual masks by hand is too time consuming given the large amount of imaging data. In this work we have shown that the great quality of the Gaia star sample and the remarkable data homogeneity of HSC allowed to build robust star masks to be used for a large number of science cases, which is promising in the context of the future LSST (Ivezić et al., 2008) and Euclid (Laureijs et al., 2011) surveys.

{ack}

The Hyper Suprime-Cam (HSC) collaboration includes the astronomical communities of Japan and Taiwan, and Princeton University. The HSC instrumentation and software were developed by the National Astronomical Observatory of Japan (NAOJ), the Kavli Institute for the Physics and Mathematics of the Universe (Kavli IPMU), the University of Tokyo, the High Energy Accelerator Research Organization (KEK), the Academia Sinica Institute for Astronomy and Astrophysics in Taiwan (ASIAA), and Princeton University. Funding was contributed by the FIRST program from Japanese Cabinet Office, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), the Japan Society for the Promotion of Science (JSPS), Japan Science and Technology Agency (JST), the Toray Science Foundation, NAOJ, Kavli IPMU, KEK, ASIAA, and Princeton University.

This paper makes use of software developed for the Large Synoptic Survey Telescope. We thank the LSST Project for making their code available as free software at http://dm.lsst.org.

The Pan-STARRS1 Surveys (PS1) have been made possible through contributions of the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, Queen’s University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation under Grant No. AST-1238877, the University of Maryland, and Eotvos Lorand University (ELTE) and the Los Alamos National Laboratory.

Based in part on data collected at the Subaru Telescope and retrieved from the HSC data archive system, which is operated by Subaru Telescope and Astronomy Data Center at National Astronomical Observatory of Japan.

This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

We thank Laurent Eyer for useful discussions about the Gaia DR1 and the referee for her/his useful comments.

Appendix A Emulated Gaia magnitude in Tycho-2

We show in Figure 16 the emulated magnitude for Tycho-2 stars into $G_{\rm gaia}$ magnitude using Equation 3.2, versus the true $G_{\rm gaia}$ magnitude for the stars in common. The transformation gives unbiased $G_{\rm gaia}$-like magnitudes.

Appendix B Content of the released bright-star catalogue and masks

The bright-star catalogue, the masks, and the tools to flag sources are distributed as a unique archive, available at this location101010Also available at http://jeancoupon.com/brightStarMasks.: ftp://obsftp.unige.ch/pub/coupon/brightStarMasks/HSC-SSP/. The repository contains a permanent link to the latest version (HSC-SSP_brightStarMask_latest.tgz) and a README file.

B.1 Version history

The current mask version described in this work is Arcturus, which differs from previous versions as follows:

• •

Arcturus (April 21th, 2017): identical to Canopus, with tract- and patch-based region files included.

• •

Canopus (early April 2017): Gaia DR1, Tycho-2 and SDSS, pure star sample, a few areas with lower Gaia-star density due to low-scanned areas.

• •

Sirius (March 2016, used in S15B, S16A and DR1 releases): “old” masks, 9% of bright galaxies are masked, the size of a dozen masks is over conservative for magnitude $\sim 5$ stars.

B.2 Uncompressing the archive

To uncompress the archive, run:

$ tar xzvf HSC-SSP_brightStarMask_VERSION.tgz

(where VERSION is the latest version).

B.3 Installing venice

venice is a mask utility program that reads a mask file (DS9 or fits type) and a catalogue of objects to:

• •

create a pixelized mask,

• •

find objects inside/outside a mask,

• •

or generate a random catalogue of objects inside/outside a mask.

The code is optimised to deal with a large number of objects and input masks. Currently, venice allows to flag 150M objects with 5M-region mask file on a desktop machine in less than an hour.

The code source is in

HSC-SSP_brightStarMask_VERSION/venice-V.V.V/

(where V.V.V is the latest stable version of venice111111For most recent versions of venice, see https://github.com/jcoupon/venice.). To compile it, you first need to install the gsl and cfitsio libraries (http://www.gnu.org/software/gsl/, http://heasarc.gsfc.nasa.gov/fitsio/). Then, go to the venice directory:

$ cd HSC-SSP_brightStarMask_VERSION/venice-V.V.V/

and run:

$ make

or, if gsl and cfisio libraries are installed in a different directory than /usr/local, run:

$ make PREFIX_GSL=DIRECTORY_NAME
PREFIX_CFITSIO=DIRECTORY_NAME

gcc is the default compile. If you wish to use a different compiler than gcc, run:

$ make CC=COMPILER_NAME

The compiled program is automatically installed in HSC-SSP_brightStarMask_VERSION/venice-V.V.V/bin/.

B.4 The bright-star catalogue

The bright-star catalogue is located in HSC-SSP_brightStarMask_VERSION/star. It contains the following columns:

1: source_id(Long) 2: ra(Double)/Angle[deg] 3: dec(Double)/Angle[deg] 4: G_Gaia(Double)/Magnitude[mag] 5: origin(String)/Gaia or Tycho-2 6: G_Gaia_SDSS(Double)/Magnitude[mag]

The source identification for the Gaia sources is the released source_id column, however, for Tycho-2 stars, it is a new running number independent from the released catalogue and the original Tycho-2 sample. So, to get a unique identification number from this catalogue, it is necessary to select both the source_id and origin identifiers.

B.5 flagging a catalogue using venice

To flag a catalogue using venice, run:

$ venice-V.V.V/bin/venice
-m reg/masks_all.reg -f all
-cat MY_INPUT_CAT
-xcol RA_COLUMN_NAME -ycol DEC_COLUMN_NAME
-o MY_OUTPUT_CAT

where reg/masks_all.reg is the master region mask file.

Note: venice can read both fits files (default) and ascii files. For input and output ascii files, set -ifmt ascii and -ofmt ascii, respectively.

B.6 Tract and patch region files

The reg/patches.tgz and reg/tracts.tgz archives contain the region mask files split in tracts and patches for the five filters. To untar it, run:

$cd HSC-SSP_brightStarMask_VERSION/reg$ tar xzvf tracts.tgz $ tar xzvf patches.tgz

The tracts archive expands into $5\,715$ files and the patches archive expands into $356\,550$ files.

The paths and names for the patches and tracts files are:

tract: HSC-SSP_brightStarMask_VERSION/reg/tracts/ BrightStarMask-TRACT-FILTER.reg patch: HSC-SSP_brightStarMask_VERSION/reg/tracts/ TRACT/BrightStarMask-TRACT-PATCH-FILTER.reg

Currently, the mask is identical for each filter, so $g$, $r$, $z$, and $Y$ filter masks are symbolic links pointing to the $i$ filter mask.

B.7 A concrete example: flagging objects in tract 9376

First, create a catalogue of random points to be used for this example:

$ venice-V.V.V/bin/venice
-r -xmin 221.476 -xmax 222.967
-ymin -1.5014 -ymax 0.0442215
-coord spher -o tract_9376.fits

Then, run the catalogue through the mask file:

$ venice-V.V.V/bin/venice
-m tracts/9376/BrightStarMask-9376-HSC-I.reg
-cat tract_9376.fits -xcol ra -ycol dec
-f all -flagName isOutsideMask
-o tract_9376_flagged.fits

Options:

• •

-m tracts/9376/BrightStarMask-9376-HSC-I.reg: masks in region format

• •

-cat tract_9376.fits: catalogue to flag

• •

-xcol ra -ycol dec: names of the input coordinates columns

• •

-f all: keep all objects, 1: outside the mask, 0: inside the mask

• •

-flagName isOutsideMask: name of the flag column

• •

-o tract_9376_flagged.fits: the output file

Appendix C Galaxy contamination in the (previously released) Sirius bright-object catalogue

An earlier version of the bright-star mask (Sirius version) was used in the S15B, S16A and DR1 releases. The description of the star catalogue and the construction of the masks are detailed in Mandelbaum et al. (in prep.). The code source used for the Sirius version can be retrieved at https://bitbucket.org/czakon/hsc-data-analysis.

To evaluate the fraction of galaxies in the Sirius catalogue we match it to the SDSS catalogue, in a single HSC-SSP tract (9376). First, we compute the emulated $G_{\mathrm{Gaia}}$ magnitudes using the SDSS PSF magnitudes, and we identify the sources flagged as extended in the SDSS. In the left panel of Figure 17, we show the magnitude distribution of point-like and extended sources, as the blue and green curves, respectively. We also show the distribution of the Arcturus sample as the black solid curve.

One can see that for extended sources, the emulated $G_{\rm Gaia}$ magnitudes extend to the fainter part, as it is computed from the SDSS PSF magnitudes (hence in the central, fainter, part of the galaxies). Visual inspection reveals that most of the extended sources are bright galaxies with a prominent bulge, as illustrated in the right panel of Figure 17.

Out of $3\,255$ objects (in the HSC-SSP tract 9376) from the Sirius catalogue, we find that $20\%$ are not matched with none of the $G_{\rm Gaia}<18$ or the bright SDSS sample. This is mostly due to some sources from the Sirius catalogue that scatter to fainter magnitude. $2\,609$ sources are matched with at least Gaia or SDSS, among which $211$ (8.9%) are flagged as extended. Note that this number is an underestimate since we did not measure the fraction of galaxies beyond the limiting magnitude of the bright-star Gaia and SDSS catalogues.