The JWST Resolved Stellar Populations Early Release Science Program V.
DOLPHOT Stellar Photometry for NIRCam and NIRISS

Daniel R. Weisz Department of Astronomy, University of California, Berkeley, CA 94720, USA Andrew E. Dolphin Raytheon, 1151 E. Hermans Rd., Tucson, AZ 85756 Steward Observatory, University of Arizona, 933 N. Cherry Avenue, Tucson, AZ 85719, USA Alessandro Savino Department of Astronomy, University of California, Berkeley, CA 94720, USA Kristen B. W. McQuinn Department of Physics and Astronomy, Rutgers, the State University of New Jersey, 136 Frelinghuysen Road, Piscataway, NJ 08854, USA Max J. B. Newman Department of Physics and Astronomy, Rutgers, the State University of New Jersey, 136 Frelinghuysen Road, Piscataway, NJ 08854, USA Benjamin F. Williams Department of Astronomy, University of Washington, Box 351580, U.W., Seattle, WA 98195-1580, USA Nitya Kallivayalil Department of Astronomy, University of Virginia, 530 McCormick Road, Charlottesville, VA 22904, USA Jay Anderson Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA Martha L. Boyer Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA Matteo Correnti INAF Osservatorio Astronomico di Roma, Via Frascati 33, 00078, Monteporzio Catone, Rome, Italy ASI-Space Science Data Center, Via del Politecnico, I-00133, Rome, Italy Marla C. Geha Department of Astronomy, Yale University, New Haven, CT 06520, USA Karin M. Sandstrom Department of Astronomy & Astrophysics, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, USA Andrew A. Cole School of Natural Sciences, University of Tasmania, Private Bag 37, Hobart, Tasmania 7001, Australia Jack T. Warfield Department of Astronomy, University of Virginia, 530 McCormick Road, Charlottesville, VA 22904, USA Evan D. Skillman University of Minnesota, Minnesota Institute for Astrophysics, School of Physics and Astronomy, 116 Church Street, S.E., Minneapolis, MN 55455, USA Roger E. Cohen Department of Physics and Astronomy, Rutgers, the State University of New Jersey, 136 Frelinghuysen Road, Piscataway, NJ 08854, USA Rachael Beaton Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA Department of Astrophysical Sciences, Princeton University, 4 Ivy Lane, Princeton, NJ 08544, USA The Observatories of the Carnegie Institution for Science, 813 Santa Barbara St., Pasadena, CA 91101, USA Alessandro Bressan SISSA, Via Bonomea 265, 34136 Trieste, Italy Alberto Bolatto Department of Astronomy, University of Maryland, College Park, MD 20742, USA Joint Space-Science Institute, University of Maryland, College Park, MD 20742, USA Michael Boylan-Kolchin Department of Astronomy, The University of Texas at Austin, 2515 Speedway, Stop C1400, Austin, TX 78712-1205, USA Alyson M. Brooks Department of Physics and Astronomy, Rutgers, the State University of New Jersey, 136 Frelinghuysen Road, Piscataway, NJ 08854, USA Center for Computational Astrophysics, Flatiron Institute, 162 Fifth Avenue, New York, NY 10010, USA James S. Bullock Department of Physics and Astronomy, University of California, Irvine, CA 92697 USA Charlie Conroy Center for Astrophysics — Harvard & Smithsonian, Cambridge, MA, 02138, USA Michael C. Cooper Department of Physics and Astronomy, University of California, Irvine, CA 92697 USA Julianne J. Dalcanton Department of Astronomy, University of Washington, Box 351580, U.W., Seattle, WA 98195-1580, USA Center for Computational Astrophysics, Flatiron Institute, 162 Fifth Avenue, New York, NY 10010, USA Aaron L. Dotter Department of Physics and Astronomy, Dartmouth College, 6127 Wilder Laboratory, Hanover, NH 03755, USA Tobias K. Fritz Department of Astronomy, University of Virginia, Charlottesville, 530 McCormick Road, VA 22904-4325, USA Christopher T. Garling Department of Astronomy, University of Virginia, 530 McCormick Road, Charlottesville, VA 22904, USA Mario Gennaro Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA The William H. Miller III Department of Physics & Astronomy, Bloomberg Center for Physics and Astronomy, Johns Hopkins University, 3400 N. Charles Street, Baltimore, MD 21218, USA Karoline M. Gilbert Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA The William H. Miller III Department of Physics & Astronomy, Bloomberg Center for Physics and Astronomy, Johns Hopkins University, 3400 N. Charles Street, Baltimore, MD 21218, USA Leo Girardi Padova Astronomical Observatory, Vicolo dell’Osservatorio 5, Padova, Italy Benjamin D. Johnson Center for Astrophysics — Harvard & Smithsonian, Cambridge, MA, 02138, USA Cliff Johnson Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA) and Department of Physics and Astronomy, Northwestern University, 1800 Sherman Avenue, Evanston, IL 60201, USA Jason Kalirai John Hopkins Applied Physics Laboratory, 11100 Johns Hopkins Road, Laurel, MD 20723, USA Evan N. Kirby Department of Physics, University of Notre Dame, Notre Dame, IN 46556, USA Dustin Lang Perimeter Institute for Theoretical Physics, Waterloo, ON N2L 2Y5, Canada Paola Marigo Department of Physics and Astronomy G. Galilei, University of Padova, Vicolo dell’Osservatorio 3, I-35122, Padova, Italy Hannah Richstein Department of Astronomy, University of Virginia, 530 McCormick Road, Charlottesville, VA 22904, USA Edward F. Schlafly Space Telescope Science Institute, 3700 San Martin Dr., Baltimore, MD 21218, USA Erik J. Tollerud Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA Andrew Wetzel Department of Physics and Astronomy, University of California, Davis, CA 95616, USA

(Accepted February 2, 2024)

Abstract

We present NIRCam and NIRISS modules for DOLPHOT, a widely-used crowded field stellar photometry package. We describe details of the modules including pixel masking, astrometric alignment, star finding, photometry, catalog creation, and artificial star tests (ASTs). We tested these modules using NIRCam and NIRISS images of M92 (a Milky Way globular cluster), Draco II (an ultra-faint dwarf galaxy), and WLM (a star-forming dwarf galaxy). DOLPHOT’s photometry is highly precise and the color-magnitude diagrams are deeper and have better definition than anticipated during original program design in 2017. The primary systematic uncertainties in DOLPHOT’s photometry arise from mismatches in the model and observed point spread functions (PSFs) and aperture corrections, each contributing $\lesssim 0.01$ mag to the photometric error budget. Version 1.2 of WebbPSF models, which include charge diffusion and interpixel capacitance effects, significantly reduced PSF-related uncertainties. We also observed minor ( $\lesssim 0.05$ mag) chip-to-chip variations in NIRCam’s zero points, which will be addressed by the JWST flux calibration program. Globular cluster observations are crucial for photometric calibration. Temporal variations in the photometry are generally $\lesssim 0.01$ mag, although rare large misalignment events can introduce errors up to 0.08 mag. We provide recommended DOLPHOT parameters, guidelines for photometric reduction, and advice for improved observing strategies. Our ERS DOLPHOT data products are available on MAST, complemented by comprehensive online documentation and tutorials for using DOLPHOT with JWST imaging data.

James Webb Space Telescope (2291); Stellar photometry (1620); Hertzsprung Russell diagram (725)

^†^†journal: ApJS^†^†facilities: JWST(NIRCAM), JWST(NIRISS)^†^†software: This research made use of routines and modules from the following software packages: Astropy (Astropy Collaboration et al., 2013, 2018, 2022), DOLPHOT (Dolphin, 2016), IPython (Perez & Granger, 2007), Matplotlib (Hunter, 2007), NumPy (van der Walt et al., 2011) and SciPy (Virtanen et al., 2020)

1 Introduction

JWST has the potential to resolve millions of stars in thousands of galaxies out to large distances (e.g., $D\sim 100$ Mpc). Such data will enable new foundational science in a broad range of areas such as the cosmic distance ladder and local $H_{0}$ measurements, reionization, globular cluster formation, dark matter, the stellar initial mass function, galaxy assembly, the effects of rare red stars that can effect the spectral energy distributions (SEDs) of galaxies at all cosmic epochs, and much more (e.g., see discussion in Weisz et al. 2023).

Much of this science comes from observations of resolved stars in crowded fields. In crowded fields, neighboring stars have overlapping point spread functions (PSFs), which can lead to confusion over the number of stars and their relative contributions to the observed flux in a given pixel. Recovering accurate and precise photometry for large numbers of stars in the limit of modest-to-severe crowding is technically daunting and requires highly optimized observations and sophisticated analysis tools (e.g., Dalcanton et al., 2012a; Williams et al., 2014).

Fortunately, crowded field stellar photometry is a mature field based on a rich history of development dating back nearly $\sim 50$ years. Early crowded field photometry routines combined pioneering work on photoelectric detectors with innovative approaches to simultaneously modeling the stellar light profiles of adjacent stars, resulting in a number of codes in the 1980s that could photometer thousands of stars in a field (e.g., Buonanno et al., 1979; Tody, 1980; Stryker, 1983; Lupton & Gunn, 1986; Penny & Dickens, 1986; Schechter et al., 1993). A major achievement of this era was the creation of the legacy software package DAOPHOT (Stetson, 1987).¹¹1As discussed by Stetson (1987), DAOPHOT was a version of the photometric routine POORMAN written by Mould & Shortridge that was improved to handle a higher density of stars. Unfortunately, there is no bibliographic record for POORMAN.

Continued improvements in crowded field photometry were catalyzed by the launch of the Hubble Space Telescope (HST). Via the Hubble Key Project aimed at measuring $H_{0}$ , several independent photometric routines were developed to gauge systematics in the photometry (e.g., Stetson 1994; and see discussion in Freedman et al. 2001). Similarly, the ground-breaking sensitivity and precision of HST/WFPC2, along with its notoriously undersampled PSF, motivated the development of specialized photometric and astrometric routines aimed at dense stellar fields (e.g., Holtzman et al., 1995; Lauer, 1999; Anderson & King, 2000; Dolphin, 2000). More recently, stellar surveys of crowded fields such as the Galactic plane and M31 have provided important gains in the speed and flexibility of crowded field codes (e.g., Dalcanton et al., 2012a; Schlafly et al., 2018). In the context of nearby galaxies, the Panchromatic Hubble Andromeda Treasury (PHAT) survey provided substantial new additions (e.g., simultaneous multi-camera, multi-wavelength crowded field photometry) to DOLPHOT (Dolphin, 2000, 2016), a crowded field photometric package that has produced photometry for millions of stars in hundreds of galaxies in the Local Group (LG) and Local Volume (e.g., Holtzman et al., 2006; Rizzi et al., 2007; Weisz et al., 2008; Dalcanton et al., 2009; McQuinn et al., 2010; Radburn-Smith et al., 2011; Dalcanton et al., 2012b, a; Williams et al., 2014; Jang & Lee, 2017; McQuinn et al., 2017; Skillman et al., 2017; Sabbi et al., 2018; Anand et al., 2021; Jang et al., 2021; Williams et al., 2021; Lee et al., 2022; Savino et al., 2022; Riess et al., 2023; Williams et al., 2023).

A main goal of our JWST Resolved Stellar Populations Early Release Science is to provide the astronomy community with an easy-to-use and efficient means for performing crowded field photometry on JWST imaging that will ultimately help to realize JWST’s full potential for the resolved Universe. Specifically, we have developed modules for DOLPHOT that are tailored to the characteristics of NIRCam and NIRISS, which are important imaging instruments for studies of resolved stellar populations with JWST. DOLPHOT is a well-tested, widely-used, and publicly available package that already supports modules specific to several HST cameras (WFPC2, ACS, WFC3/UVIS and IR), has been a testing ground for the Roman Space Telescope, and includes general purpose routines that can be used on virtually any images of resolved stars. The addition of JWST modules will enable a wide array of JWST-specific and cross-facility (e.g., JWST and HST) science, some of which has already been demonstrated using early versions of our JWST DOLPHOT modules (e.g., Chen et al., 2023; Lee et al., 2023a, b; McQuinn et al., 2023; Riess et al., 2023; Van Dyk et al., 2023; Warfield et al., 2023; Peltonen et al., 2024; Li et al., 2024).

Table 1: A summary of our JWST Early Release Science (ERS) observations taken in 2022. Although we acquired 4 exposures for M92, the 3rd exposure produces poor photometry due to larger-than-normal jitter in the telescope stability that occurred only during this exposure. We examine this issue in Appendix A. The first entry for M92 in this table reflects all observations taken, while the second entry is without the 3rd exposure. More details on the exact observations (e.g., dither pattern, readout mode) are given in Weisz et al. (2023) and are available in our public Phase II file in APT.

\topruleTarget	Date	Camera	Filter	$t_{exp}$ [s]	Groups	Integrations	Dithers
\topruleM92	June 20–21	NIRCam	F090W/F277W	1245.465	6	1	4
		NIRCam	F150W/F444W	1245.465	6	1	4
		NIRISS	F090W	1245.465	7	1	4
		NIRISS	F150W	1245.465	7	1	4
M92 (no 3rd exp)	June 20–21	NIRCam	F090W/F277W	934.099	6	1	3
		NIRCam	F150W/F444W	934.099	6	1	3
		NIRISS	F090W	934.099	7	1	3
		NIRISS	F150W	934.099	7	1	3
Draco II	July 3	NIRCam	F090W/F480M	11810.447	7	4	4
		NIRCam	F150W/F360M	5883.75	7	2	4
		NIRISS	F090W	11123.294	9	7	4
		NIRISS	F150W	5883.75	10	3	4
WLM	July 23–24	NIRCam	F090W/F430M	30492.427	8	9	4
		NIRCam	F150W/F250M	23706.788	8	7	4
		NIRISS	F090W	26670.137	17	9	4
		NIRISS	F150W	19841.551	19	6	4
\toprule

In this paper, we describe the NIRCam and NIRISS stellar photometry modules for DOLPHOT. DOLPHOT’s underlying algorithms are already well-documented in the literature, along with rigorous tests of their accuracy in a variety of crowded and uncrowded fields (e.g., Radburn-Smith et al., 2011; Williams et al., 2014, 2023). Accordingly, our focus is on describing the details specific to the DOLPHOT NIRCam and NIRISS modules and providing examples of its application to the three Early Release Science (ERS) targets in the Local Group: M92, Draco II, and WLM. As described in Weisz et al. (2023), these targets, and the associated observing strategies, were carefully selected to benchmark the development of DOLPHOT in a variety of regimes that we anticipate will be common for resolved star science, and thus need to be vetted for study with JWST. This paper is designed to describe the modules and provide examples of their application to the ERS data. As part of the ERS program, we have created an extensive set of deliverables including online documentation and data products that allow interested readers to reduce ERS data identically to what is done in this paper, as well as explore the various aspects of the data for their own purposes (e.g., to customize catalog culling criteria). Essential DOLPHOT input and output data associated with the photometric reductions in this paper are hosted as high-level science products on MAST²²2https://archive.stsci.edu/hlsp/jwststars/, while step-by-step guides for our DOLPHOT reductions can be found on our DOLPHOT documentation page³³3https://dolphot-jwst.readthedocs.io.

This paper is organized as follows. We summarize the ERS observations in §2. In §3, we describe the DOLPHOT NIRCam and NIRISS modules and provide a general outline of how to apply these new modules to JWST imaging in order to produce stellar catalogs. We illustrate the application of these modules to ERS data in §4. In §5, we examine the time variability of the PSF and compare the DOLPHOT SNR estimates with expectations from the JWST ETC. Finally, we summarize the paper and highlight future areas for improvement in §6.

2 Observations

Extensive details of our ERS survey and observations are provided in Weisz et al. (2023). Here, we briefly summarize the observations and list their basic characteristics in Table 1.

In June and July, 2022 our program acquired NIRCam and NIRSS imaging of three LG targets: globular cluster M92, ultra-faint dwarf galaxy Draco II, and LG star-forming dwarf galaxy WLM. These targets were selected to satisfy a number of science and technical goals including the development and testing of DOLPHOT in a variety of scenes (e.g., crowded and uncrowded fields, varying surface brightness, varying degrees of saturation, a representative set of wide and medium filters). In all cases, the NIRCam fields were placed centrally on each target with locations and orientations set to maximize overlap with archival HST imaging and schedulability early in the ERS window. The NIRISS fields were acquired in parallel. Table 1 summarizes basic characteristics of our NIRCam and NIRISS observations.

In the process of analyzing our M92 data, we found that the third exposure of M92 appears to be “corrupt” in the sense that although the 3rd exposure of M92 visually looks fine, it results in remarkably poor photometry for both NIRCam and NIRISS, despite extensive efforts to fix it. We therefore have excluded it from the DOLPHOT reductions in this paper. A later analysis of the fine guidance sensor data revealed instabilities in the telescope only during this exposure. We discuss details of the third exposure in Appendix A.

Our observations of WLM were designed to sample RR Lyrae light curves. However, the default JWST reduction pipeline is currently not capable of producing the time series images for short period observations needed to extract flux as a function of time. The default pipeline currently only provides time series data when the time series observation (TSO) mode is used. We were unable to use TSO mode because it prohibits dithering. While it is possible to modify the JWST pipeline to produce the time series images necessary for short period variable analysis, it is a topic beyond the scope of this paper. Instead, a Cycle 2 archival proposal undertaken by some members of our team (AR-03248; PI Skillman) is developing the documentation and tools needed to recover variable star light curves from NIRCam and NIRISS imaging with DOLPHOT. These will be made publicly available upon completion.

Finally, the NIRISS observations of Draco II, were located at several half-light radii from the galaxy. As far as we can tell, the field is consistent with being blank, i.e., no obvious galaxy member stars, and we do not analyze or discuss this field in the paper.

3 NIRCam and NIRISS Photometry with DOLPHOT

In this section, we provide an overview of the new DOLPHOT NIRCam and NIRISS modules with application to imaging from our ERS program. The core workings of DOLPHOT, along with extensive tests of its functionality and reliability are well-documented in the literature (e.g., Dolphin, 2000; Dalcanton et al., 2012b, a; Williams et al., 2014; Dolphin, 2016; Williams et al., 2021). Here, we will not re-visit these details. Instead, we focus on modifications made to DOLPHOT for incorporating NIRCam and NIRISS imaging into its existing framework.

3.1 Overview

As input, DOLPHOT takes a reference image, a list of science images, and several dozen input parameters with user defined values. DOLPHOT astrometically aligns all the science images to a reference image. It then performs simultaneously, multi-wavelength photometry on all science images by fitting PSF models to the signal-to-noise ratio (SNR) peaks (and any neighboring SNR peaks in crowded fields) it detects in each of the science images. The result of this process is a set of photometric measurements for all detected objects in all science images. Stellar catalogs are created by culling the main source catalog using criteria such as SNR, how compact/extended sources are, etc., which we detail in §3.5.

For JWST, the science images we use are the stage 2 cal files, which are calibrated single exposure science images produced by the JWST pipeline. The references images are I2D files, which are resampled, stacked science images, akin to drizzled images with HST. It is possible to use cal as reference images as well. The JWST pipeline also makes available crf images, which are like cal science images, but with cosmic ray flags applied. Generally, we have found astrometric alignment of the crf images in DOLPHOT to be worse than cal images.

The reference image is used only to align each of the science images. Our general recommendation for a reference image is to select the deepest image available. All images used in this analysis were created by the standard STScI pipeline and downloaded from MAST⁴⁴4The specific observations used in his paper can be accessed via https://doi.org/10.17909/71kb-ga31 (catalog DOI: 10.17909/71kb-ga31)..

In this paper, we use FITS images with the following JWST pipeline versioning information CAL_VER $=$ 1.11.4, CRDS_VER $=$ 11.17.2, and CRDS_CTX $=$ jwst_p1147.pmap. This version includes updates to the chip-to-chip zero points, the switch from Vega to Sirus as a reference star, and updated flat fields released in late 2023.

DOLPHOT requires only a few pre-processing steps for all images. These steps include masking bad pixels (e.g., cosmic rays, hot pixels) based on the data quality (DQ) flags, multiplying by the camera specific pixel area masks, and making initial estimates of the sky and the positions of bright stars for alignment.

Following pre-processing, DOLPHOT aligns each science image to the reference image. The quality of the alignment is determined by several factors including the number of bright stars available, depth of the reference image, the fidelity of the provided WCS information, SNRs, relative orientations (e.g., images with large rotations may be harder to align), stars in common between the reference and science image (e.g., images taken in very different filters, such as ultra-violet and IR, may be hard to align as they may have few sources in common).

With all images aligned, DOLPHOT searches for objects to photometer by iteratively identifying signal-to-noise ratio peaks in the stack of science images. DOLPHOT measures the fluxes of each object by simultaneously fitting a PSF model, and a local background model, to the target object plus all neighboring objects within a user specified radius.

Upon completion of photometry, DOLPHOT provides extensive output including its position on the reference image and the flux and a number of quality assessment metrics (e.g., $\chi^{2}$ , shape of the star’s light profile relative to the PSF) for each star in each image. It also provides combined fluxes and magnitudes for each object from which stellar catalogs are usually constructed.

Characterizing uncertainties for crowded field stellar photometry requires artificial star tests (ASTs). ASTs are synthetic stars with known positions and magnitudes that are inserted into real JWST science images and then recovered by DOLPHOT. It is well-established that the difference in input and recovered flux for ASTs provides a more realistic accounting of photometric uncertainties than the Poisson noise that is reported by the crowded field photometric process alone (e.g., Stetson & Harris, 1988).

3.2 Pre-processing steps

Prior to running DOLPHOT, pre-processing is required in order to convert the data to a format suitable for PSF-fitting photometry. For the case of NIRCAM and NIRISS data, steps in this process are as follows (using the nircammask and nirissmask utilities, respectively):

•

Mask out bad or saturated pixels. At the time of this writing, bad pixels on cal and crf images are identified by having an SCI array value of NaN; previous versions of the pipeline have used SCI array values of exactly 0. The mask utilities will correctly interpret either approach. Additionally, saturated pixels in cal and crf images are identified by having a DQ array flag with a value of 2. Bad pixels on I2D images are identified by having a WHT array value of exactly 0.
•

Convert from the default calibration of MJy/sr to DN (data number). This is performed by dividing all pixel values by the FITS keyword PHOTMJSR, and subsequently multiplying by the exposure time (FITS keyword EFFEXPTM)⁵⁵5We note that as of this writing, the actual time the telescope spends collecting data slightly differs from the exposure time in the FITS keywords, such as EFFEXPTM. In practice, the ramp fitting procedure begins at the end of the first group. But currently, the FITS exposure time keywords are based on when the first group starts. This effect is generally subtle, i.e., the impact on SNR is typically $\lesssim 1$ %, but it is now factored into DOLPHOT for accuracy and completeness.. An additional step for cal and crf images is to multiply pixel values by the pixel area map (AREA array).
•

Readout noise and gain values are also saved into the FITS file. Details of the pedigree of that data are available in the versions.txt files that are included with the NIRCAM and NIRISS DOLPHOT modules.

The result of the pre-processing step is an image in units of DN, along with FITS keywords for gain and readout noise, allowing PSF-fitting photometry to run.

As of this writing, for the purpose of this ERS program, there was no need to incorporate information from the Advanced Scientific Data Format (ASDF) metadata provided for JWST images (Greenfield et al., 2015). ASDF has the ability to host more detailed metadata than the standard FITS format, such as improved WCS information. However, at this time all necessary information (astrometry information, photometric calibrations, pixel areas, etc.) for DOLPHOT are available via FITS header keywords. If that situation changes in the future, we will explore writing an ASDF reader for DOLPHOT.

3.3 Alignment

The first step in DOLPHOT’s reduction process is to align all of the original-sampling (cal or crf) images to a common reference frame. Normally, the common reference is an I2D file. In principle, if no dithering was used in the observations, one of the cal or crf images could be used as a reference image. In practice, typically the deepest I2D file makes for the best reference image as it allows for the most star matches with the science images leading to better astrometric alignment. We note that DOLPHOT does not re-sample or re-bin any images, which can lead to issues in conservation of intensity, for example. Instead, DOLPHOT only performs photometry on the original, non-drizzled, science images provided by the JWST pipeline.

A key difference between the JWST modules and the previously released DOLPHOT HST modules is that DOLPHOT does not apply distortion corrections based on DOLPHOT’s internal model. With the HST modules, common practice is to use the astrometric data in the header but not the distortion. Using HST with DOLPHOT’s internal distortion model requires the parameter setting UseWCS=1, for which DOLPHOT estimates only the shift, scale, and rotation of each science image relative to the reference. For JWST the astrometric data included in the FITS header is used for both alignment and application of distortion corrections, requiring a setting of UseWCS=2. In this case, DOLPHOT estimates a full distortion solution, which is beyond the shift, scale, and rotation typically used by DOLPHOT on HST images. As of this writing, NIRCAM data are provided with 3rd order SIP polynomials, while NIRISS data are provided with 4th order polynomials. DOLPHOT currently processes up to 5th order polynomials, so it can handle additional fidelity in the astrometry data, should it become available in the pipeline or added by offline astrometry tools (e.g., astrometry.net; Lang et al. 2010).

3.4 Star Detection and Photometry

As with previously available DOLPHOT modules, photometry proceeds once the images are aligned. An initial pass (and usually multiple passes) detects peaks in the SNR map across the reference image. Photometry is performed at each peak to attempt to identify a point source. As with all previous versions of DOLPHOT, the user can adjust parameters to alter the noise models, photometry modes (e.g., aperture vs. PSF-fitting), sky fitting method, etc. Our internal testing on the ERS data resulted in a set of recommended DOLPHOT parameters, which are listed in Table 4. We discuss the process by which we determined these parameters in more detail in §3.8.

A standard feature of DOLPHOT is to make adjustments to the model PSFs to improve the fit quality and photometry. As discussed in Dolphin (2000), mismatches between the shape of the model and the true PSF contribute to the photometric error budget, which can grow large for faint sources. As part of its normal operation, DOLPHOT measures a PSF residual image relative to the pre-calculated PSF model library. It then makes adjustments to the model PSFs based on comparisons to bright stars in the field to improve agreement between the model and observed PSFs. This step is performed automatically (i.e., PSFRes=1) by DOLPHOT unless the use of residual PSF images is turned off (i.e., PSFRes=0). These adjustments are made independently in each DOLPHOT science image. The amplitude of the PSF adjustments provides a means to quantify the systematic uncertainty floor on the photometry for a given set of PSF models. In §3.6, we provide the typical PSF adjustments made by DOLPHOT on ERS data relative to the default WebbPSF models.

Likewise, aperture corrections are normally computed automatically (i.e., ApCor=1), unless they are turned off in DOLPHOT (i.e., ApCor=0). This calculation will estimate the magnitude difference between instrumental PSF-fitted magnitudes and magnitudes within a standard 10-pixel radius, accounting for the WebbPSF-predicted encircled energy within that standard radius.

Finally, zero points are applied to convert from instrumental magnitudes (in DN/sec) to the VEGAMAG system. These zero points are in units of Jy for Vega (though Sirus is now the standard reference star) so they require conversion back from DN/sec to Jy using the FITS header keywords PHOTMJSR and PIXAR_SR before the calibration is applied. Zero points for NIRISS were provided relative to 1 DN/sec. Thus, they are applied directly without additional conversion. Finally, we note that DOLPHOT parameters NIRCAMvega and NIRISSvega can be set to zero to report in ABmag instead of VEGAMAG.

3.5 Post-Processing & Catalog Creation

Once complete, DOLPHOT saves all photometry data to a single ASCII file containing overall fit metrics (positions of objects in the coordinate system of the reference image, $\chi^{2}$ , S/N, sharpness, roundness, crowding, and the object type: single pixel, point source, or extended source) for all objects identified. Photometry and the same quality assessment information is also provided for all combined exposures in each filter (e.g., all F090W images), as well as for all individual exposures (e.g., each F090W images, which can be used for time domain studies, for example). The photometric data provided by filter and image include counts (DN), background, calibrated magnitude, and calibrated count rates, which can be useful if an upper limit is informative, such as in multi-wavelength SED fitting and time-domain studies (e.g., Gordon et al., 2016).

As the DOLPHOT output includes all sources identified on the images, it is necessary to establish criteria for good detections (i.e., stars). As part of this ERS program, Warfield et al. (2023) developed a set of criteria to identify stars using DOLPHOT reported quality parameters. We adopt this scheme for this paper and classify good stars as those that satisfy all of the following criteria:

•

$SNR_{F090W}\geq 5$
•

$SNR_{F150W}\geq 5$
•

$sharp_{F090W}^{2}\leq 0.01$
•

$sharp_{F150W}^{2}\leq 0.01$
•

$crowd_{F090W}\leq 0.5$
•

$crowd_{F150W}\leq 0.5$
•

$flag_{F090W}\leq 3$
•

$flag_{F150W}\leq 3$
•

Object Type $\leq 2$

The sharpness parameter is zero for a perfectly-fit star, positive for a star that is too sharp (i.e., the flux is concentrated in a small number of pixels, e.g., a cosmic ray), and negative for a star that is too broad (perhaps a blend, cluster, or galaxy). Our choice of SNR $\geq 5$ is lower than the SNR $\geq 10$ threshold in Warfield et al. (2023). Warfield et al. (2023) focused on optimizing the other parameters for star-galaxy separation, and therefore adopted a more conservative SNR threshold.

The crowding parameter is in magnitudes. It reports how much brighter the star would have been measured had nearby stars not been fit simultaneously. For an isolated star, the value is zero. High crowding values are generally a sign of poorly-measured stars.

Error flags are defined as: 0 is a star that is recovered extremely well; 1 is that the photometry aperture extends off chip; 2 is that there are too many bad or saturated pixels; 4 is the center of the star is saturated; 8 is an extreme case of one of the above. The DOLPHOT manual suggests using values of 3 or less in general or 2 or less for precision photometry.

Object types are as follows: 1 is a good star; 2 is a star too faint for PSF determination; 3 is an elongated object; 4 is an object that is too sharp; 5 is an extended object. As recommended by the DOLPHOT manual, we only keep object types 1 and 2 in our stellar catalogs.

As the deepest and highest angular resolution images, we found that applying these criteria to F090W and F150W photometry had the largest impact on the catalog culling (e.g., Warfield et al., 2023). For targeted science (e.g., luminous red stars at longer wavelengths) other criteria and/or application to other filters may produce more desirable results. Similarly, cuts on single bands may be useful for particular science cases beyond CMDs (e.g., stellar SED fitting; Gordon et al. 2016). Finally, as discussed by Warfield et al. (2023), these criteria were focused on purity rather than completeness. This is motivated by the large number of background galaxies present in our ERS imaging. Less stringent cuts, particularly in sharpness, can produce deeper CMDs with less conservative completeness limits, albeit with a larger degree of non-stellar contamination. We illustrate the effects of our fiducial culling criteria on our ERS targets in §4. Readers who wish to explore alternative culling criteria can download our catalogs from the MAST high levels science products page.

\topruleFilter	$N_{exp}$	Central Pixel	Central Pixel	Central Pixel	Photometric Error
		Model Mean	Model $\sigma$	Mean $\delta$
		(% flux)	(% flux)	(% flux)	(mag)
(1)	(2)	(3)	(4)	(5)	(6)
\topruleNIRCAM F090W	$88$	$24.18$	$3.98$	$-0.69$	$0.005$
NIRCAM F150W	$88$	$16.63$	$1.41$	$-0.25$	$0.001$
NIRCAM F250M	$8$	$20.19$	$2.07$	$-0.24$	$0.001$
NIRCAM F277W	$6$	$18.41$	$1.68$	$-1.0$	$0.005$
NIRCAM F360M*	$8$	$12.98$	$0.76$	$-1.57$	$0.008$
NIRCAM F430M	$8$	$10.22$	$0.44$	$0.07$	$0.000$
NIRCAM F444W	$6$	$10.05$	$0.43$	$-0.34$	$0.002$
NIRCAM F480M*	$8$	$8.48$	$0.30$	$-0.41$	$0.002$
NIRISS F090W	$7$	$31.47$	$5.57$	$-2.31$	$0.014$
NIRISS F150W	$7$	$28.53$	$4.48$	$-1.69$	$0.010$
\toprule

Table 2: Central pixel PSF data for the input WebbPSF PSF models (mean and standard deviation) and mean adjustments as measured by DOLPHOT from application to all ERS NIRCam and NIRISS images, except the 3rd exposure for M92. Values in columns 3-5 are fractions of the total stellar flux. Column 6 shows approximate photometric error in magnitudes caused by the application of the same PSF residual to all stars (see Equation 1). Asterisks for NIRCAM F360M and F480M denote the filters only observed in the very sparse Draco II field; results from those filters have correspondingly higher uncertainties.

3.6 Point Spread Function Models

The currently available NIRCAM and NIRISS modules incorporate model PSF libraries calculated using WebbPSF. For this paper, the PSF models were generated using WebbPSF version 1.2.1, which adopts “in-flight” optical performance data (as opposed to pre-launch data). The alignment and stability of JWST, including the wavefront and PSF, are known to vary over time (e.g., McElwain et al., 2023), which has the potential to affect photometry. WebbPSF incorporates time-dependent optical path delay (OPD) maps to capture changes to the wavefront and PSF over time, enabling corrections for temporal changes to JWST.

The WebbPSF model PSF library we implement in DOLPHOT consists of distorted PSF models for all the available NIRCam/NIRISS filters, oversampled by a factor of 5, calculated over a $51\times 51$ physical pixel region. The models are calculated on a $5\times 5$ spatial grid for each of the detector chips. The models were generated using OPD maps from July 24th 2022 (O2022072401-NRCA3_FP1-1.fits). We used a G5V source spectrum from the Phoenix stellar library (e.g., Husser et al., 2013) to generate the PSF, which was sampled at 21, 9, and 5 wavelengths for wide, medium, and narrow bands, respectively. v1.2.1 of WebbPSF includes the effects of charge diffusion and interpixel capacitance, which were not incorporated into previous WebbPSF models. We found that the inclusion of these effects dramatically improves the quality of DOLPHOT photometry, including reducing photometric systematics by nearly an order-of-magnitude relative to WebbPSF models without these effects. We discuss the the total photometric error budget in §5.3.

As described in §3.4, DOLPHOT makes adjustments to the model PSFs to provide improved matches to the data. We summarize the effect of these PSF adjustments on the photometry in Table 2. This table provides the mean fractional central pixel brightness in each filter (i.e., the fraction of total PSF light in the central pixel averaged across all PSF models), the scatter in the central pixel fraction flux, and the mean PSF adjustment in the central pixel measured by DOLPHOT. We computed these quantities across NIRCam and NIRISS images from our ERS program, except for the 3rd exposure of M92.

As DOLPHOT provides an average PSF correction (i.e., by computing the PSF adjustments on a set of bright, high SNR stars in each science image), we can quantify the uncertainty in this correction. Table 2 provides a simple estimate of the 1- $\sigma$ photometric error created by application of the same PSF residual image to all stars in a given science image. The PSF adjustments are computed separately for each science image for each target. That is, DOLPHOT has improved the model PSFs by adjusting the central pixel on average. This is an improvement over using the native WebbPSF models without any adjustment, but it does still yield an uncertainty floor. Considering the central pixel only, as DOLPHOT does for its PSF correction (see Dolphin 2000), the magnitude error can be estimated as

\sigma=\frac{2.5\,\sigma_{\rm model}\,\delta_{\rm mean}}{\log(10)\,\mu_{\rm model% }^{2}}

(1)

where $\mu_{\rm model}$ is the mean model flux in the central pixel across all model PSFs in a given filter, $\sigma_{\rm model}$ is the scatter in the central pixel fluxes across all model PSFs in a given filter, and $\delta_{\rm mean}$ is the mean flux adjustment made by DOLPHOT. While this estimate is reasonably accurate for highly concentrated PSFs, such as the two NIRISS filters in our data sample, it becomes increasingly conservative for broader PSFs.

As Table 2 shows, the typical PSF adjustments, and the corresponding photometric uncertainties are small. For all NIRCam filters, a modest amount of light is concentrated in the central pixel and the PSF adjustments are only typically a small fraction of the total flux. The resulting photometric uncertainties introduced by the PSF models range from 0.001 to 0.008 mag. For NIRISS, the PSFs are slightly sharper, i.e., more of the total light is centrally concentrated in the model PSFs, which leads to larger PSF adjustments. The corresponding photometric uncertainties are $\sim 0.01$ mag in the NIRISS F090W and F150W filters. Overall, the small adjustments needed by DOLPHOT to improve the WebbPSF models indicate that the models themselves are already quite good. This was not the case with PSF models from previous versions of WebbPSF, all of which were systematically too sharp and the corresponding photometric errors, i.e., systematics from the PSF models, were at times larger than the photon counting noise.

We note that the values listed in Table 2 were computed for SW and LW images independently. We found all PSF adjustments to be marginally larger when running SW and LW images simultaneously, though the qualitative finding that the PSF adjustments are quite small (i.e., $\lesssim 0.01$ mag) still holds.

3.7 Artificial Star Tests

While DOLPHOT provides an estimate of photometric uncertainties based on the goodness of fit and the noise characteristics of the data, a much better characterization of the photometric measurement uncertainties and selection function are accomplished through artificial star tests (ASTs). This long-established approach (e.g., Stetson, 1987; Stetson & Harris, 1988) represents the “gold standard” in the field of resolved stellar population photometry. It relies on the injection of mock stellar sources into the raw images, which are then recovered using the identical photomeric procedure used to construct the raw and stellar DOLPHOT catalogs. The output of such simulated data tests can be used to quantify a number of aspects of data quality. A common example is that the comparison between the input and output magnitudes of the mock catalog, as well as the fraction of mock stars that are successfully detected, provides a self-consistent characterization of photometric errors, systematic uncertainties, and photometric completeness as a function of spatial position in the images and location on the CMD. Throughout this paper, we focus on ASTs run only on the SW data, which illustrate the main points of ASTs. The same procedures we describe can be used to run ASTs in an arbitrary number of bands, though the computational time can become quite expensive for large numbers of photometric bands.

The first step in running ASTs is to create a suitable input star catalog. For each target and camera, we created a list of $\gtrsim 4\times 10^{5}$ mock stars (see Table 4.1.1), with positions drawn from uniform spatial distributions on the NIRCam/NIRISS footprints, aside from gaps between the chips and modules. For each star, we assign input magnitudes such that they are uniformly distributed over the F090W vs. (F090W-F150W) CMD with $17\leq$ F090W $\leq 31$ and $-0.5\leq$ F090W $-$ F150W $\leq 2$ .

The artificial stars are then injected into all science (i.e., cal) images using the best PSF model and realistic noise obtained from the original reduction run. The stars are injected one at the time, to avoid altering the crowding properties of the original images. The star magnitude and position is then measured by DOLPHOT, as if it were a real source. Performing this operation for all the input stars, we obtain a catalog of output magnitudes, positions, and goodness-of-fit parameters. This list is then culled using the same quality criteria applied to the original photometric catalog (cf. § 3.5).

3.8 Optimizing Photometric Parameters

DOLPHOT is a flexible code in that it provides the user extensive control over details of the photometric reduction (e.g., PSF adjustments, sky fitting, image alignment methods, noise models). This ensures that DOLPHOT can produce high-quality photometry for a wide variety of images (e.g., crowded vs. uncrowded, presence/absence of surface brightness gradients, images from multiple filters with widely different characteristics). At the same time, this flexibility is characterized by a large number of parameters that can require tuning in order to produce “optimal” photometry. DOLPHOT parameters have been refined in the context of major HST programs over the past decade, culminating in a set of parameters recommended by the PHAT program that encompass a wide range of image properties (e.g., Williams et al., 2014). While these parameters produce excellent HST photometry, it is important that we investigate parameters that may provide better photometry for JWST.

Accordingly, we performed a large set of DOLPHOT runs to explore the effect of changing a select set of DOLPHOT parameters. We performed these runs on all three ERS targets, to experiment with different stellar density regimes: high crowding for WLM, low crowding for M92, and an almost empty field for Draco II. As the SW and LW channels have different detector characteristics (e.g., plate scale), we executed the experiments on the two channels independently. For each target and each channel, we explored different values of FitSky (which sets the method for local sky measurement), RAper (which sets the size of the aperture in which photometry is performed), and Rchi (which sets the size of the region over which the fit is evaluated). We set the value of FitSky to either 2 (fit the sky inside the PSF region but outside the photometry aperture) or 3 (fit the sky within the photometry aperture as a 2-parameter PSF fit). For a given FitSky, we let RAper range between 1 and 5 pixels, in discrete increments, and Rchi range between 0.5 and RAper, in increments of 0.5 pixels. In the FitSky=2 case, the values of Rsky2 (which set the inner and outer radius for sky computation) are also adjusted and set to {RAper+1; 2.5(RAper+1)}. For FitSky=3, we explored an additional grid, defined by RAper values of 7, 10, and 13 and Rchi values of 1.5, 2, and 3. This exploration results in 69 parameter permutations per channel, per target, totaling 414 DOLPHOT runs, which consumed nearly five years of CPU time.

For each of these runs, we used thousands of ASTs (see § 3.7) as one metric for evaluating photometric performance. ASTs were injected at various locations on the CMD (from very bright, high-SNR parts down to the detection limit). We then used these stars to quantify the mean scatter in color and magnitude, and the completeness fraction at each CMD location. Inspection of these metrics revealed that in large regions of the parameter space, DOLPHOT performance was poor. For example, most runs with $\texttt{Rchi}>3.0$ and/or FitSky=3 produced obviously poor photometry (e.g., low number of stars, poor completeness, large scatter). Adopting Fitsky=2, we identified a small region of the RAper-Rchi parameter space ( $2\leq\texttt{RAper}\leq 4$ and $1.5\leq\texttt{Rchi}\leq 2.5$ ) where DOLPHOT provided the “best” photometry. Within this parameter region, the photometric performance was fairly comparable. Different permutations of these parameters produced slight trade-offs in completeness versus photometric precision, though the differences were generally at the few percent level or less. We also found a slight trend with crowding, in that the optimal RAper value would increase as the field became less crowded, but again, this effect was small.

Given the similarity of the photometry over this parameter space, we decided to adopt a single set of parameters for each instrument and channel. Part of the motivation behind this choice is to provide the community with easy-to-use guidance for DOLPHOT reductions that also produces high-quality photometry. The recommended PHAT parameters all live within the “optimal” DOLPHOT parameters we have identified for JWST. Therefore, we decided to adopt the PHAT set up as our JWST DOLPHOT parameters. Specifically, we recommend the PHAT WFC3/IR parameters for the NIRCam SW channel and the ACS/WFC parameters for the NIRCam LW channel and NIRISS. The full parameter set is listed in Table 4.

Refer to caption — Figure 1: The spatial distribution of objects and stars detected by DOLPHOT in our NIRCam imaging of M92, excluding the 3rd exposure. Left: A density map ( $\sim 50\times 50$ pixel bins) of the $\sim 9.8\times 10^{5}$ objects reported in the raw DOLPHOT catalog. There is a density gradient toward the center of the image and cluster center. Additional features include a high density of sources that trace saturated stars and a lower density of sources in the SW chip gaps. Right: The density map of stars that passed the catalog culling criteria listed in § 3.5. These criteria removed the vast majority of obvious artifacts (e.g., corresponding to saturated stars, diffraction spikes) and reveal a clear stellar density gradient as is expected for a GC. Note that the individual SW chips are labeled.

\topruleTarget	Camera	Object Density	Stellar Density
		( $N$ / arcsec ${}^{2}$ )	( $N$ / arcsec ${}^{2}$ )
(1)	(2)	(3)	(4)
\topruleM92	NIRCam	28.1	2.8
	NIRISS	7.4	0.2
WLM	NIRCAM	48.8	13.2
	NIRISS	8.9	0.7
Draco II	NIRCAM	14.7	0.03
\toprule

\topruleDetector	Parameter	Value
\topruleNIRCam/SW	RAper	2
NIRCam/SW	Rchi	1.5
NIRCam/SW	Rsky2	“3 10”
NIRCam/LW	RAper	3
NIRCam/LW	Rchi	2.0
NIRCam/LW	Rsky2	“4 10”
NIRISS	RAper	3
NIRISS	Rchi	2.0
NIRISS	Rsky2	“4 10”
All	FitSky	2
All	PSFPhotIt	2
All	PSFPhot	1
All	SkipSky	1
All	SkySig	2.25
All	SecondPass	5
All	SigFindMult	0.85
All	MaxIT	25
All	NoiseMult	0.1
All	FSat	0.999
All	FlagMask	4
All	ApCor	1
All	Force1	0
All	PosStep	0.25
All	RCombine	1.5
All	SigPSF	5.0
All	PSFres	1
All	InterPSFlib	1
All	UseWCS	2
All	CombineChi	0
\toprule

The JWST Resolved Stellar Populations Early Release Science Program V. DOLPHOT Stellar Photometry for NIRCam and NIRISS