HCI ETC Instructions
The JWST Exposure Time Calculator (ETC) develops and evaluates the complex, multi-source astronomical scenes that are characteristic of JWST high-contrast imaging (HCI).
On this page
ETC functionalities for coronagraphy
The standard coronagraphic sequence combines two complementary PSF subtraction strategies:
- Referenced differential imaging (RDI), which involves subtracting a coronagraphic image of a nearby PSF reference star.
- Angular differential imaging (ADI), which involves differencing two coronagraphic images of the science target taken at different spacecraft orientations (e.g., seperated by a 10° roll).
Thus, when fully implemented, the standard coronagraphic sequence involves a minimum of three observations:
- A science observation with the host centered in the coronagraphic mask.
- A second science observation following a telescope roll, with the host centered in the coronagraphic mask.
- A PSF reference observation with the PSF reference star centered in the coronagraphic mask.
Adding small grid dithers (SGDs) of the PSF reference star is a future variant of the standard coronagraphic sequence. SGDs can increase the confidence that misalignments between the position of the science PSF and reference PSF relative to the coronagraph masks have a minimal impact on the final contrast.
The example ETC computations below capture the current scope of functionalities for coronagraphy and are limited to only the standard coronagraphic sequence elements 1 and 3, listed above, because the ADI and SGD processing is not currently supported by the ETC (V. 4.0).
The ETC does not support either ADI or SGD—only RDI.
In its current implementation (V. 4.0), the ETC is mainly useful for 2 tasks: (1) investigating detector saturation and (2) computing the signal-to-noise ratio (SNR) of a faint companion source under the ideal contrast assumption. Ideal contrast is the most optimistic assumption possible because it assumes one type of noise is dominant: the counting statistics of collected photons (shot noise). See HCI Contrast Considerations for more information about the "ideal" assumption.
The ETC treats residual flat field noise. The flat field error is a division by ~1 (the flat field is normalized), with a variance of 1/f_electrons. Note that the value of the flat field response is constant for multiple exposures and multiple integrations, so Nexposures > 1 does not decrease the residual flat field noise. To reduce it, a user either has to improve the flat field or dither with >1 pixel offsets. The most apparent effect for everyday ETC use is that residual flat field noise sets an upper limit on the achievable SNR.
The SNR source must lie within a square centered on the coronagraphic mask, aligned with detector rows and columns, with sides of 101 pixels for NIRCam or 81 pixels for MIRI.
This article is about ETC functionalities specific to coronagraphy. The reader is encouraged to become acquainted with the general ETC documentation, which covers the underlying algorithms, synthetic astrophysical scenes, simulated JWST exposures, step-by-step ETC operations, and best practices. The reader should take particular cognizance of this article: JWST ETC Coronagraphy Strategy.
Avoiding saturation
Words in bold are GUI menus/
panels or data software packages;
bold italics are buttons in GUI
tools or package parameters.
The ETC separately flags pixels that are expected to saturate at some point within a ramp. This feature lets users make subtle distinctions to deal with saturation. For example, if a ramp goes into saturation at, say, Ngroups = 10, the slope may still be accurately recoverable by discarding from the analysis only the individual frames where saturation has occurred. See JWST ETC Saturation for more information and specific saturation limits.
The user checks for saturation using the Saturation tab in the Images pane on the Calculations page in the ETC.
In some cases, warnings in the ETC Reports panel will indicate whether some ramps may still be useful. The user must proceed by trial and error, varying the readout pattern and Ngroups until all pixels at the expected position of the companion are not saturated. Note that if at least one pixel in the scene is saturated, the ETC will produce a warning.
PSF subtraction strategy
The ETC engine is designed to support predictions of coronagraphic performance, including estimates of the PSF using one or multiple reference images. The Strategy tab of the GUI calls for choices of reference target and PSF calibration method. Since March 2018 (patch release ETC 1.2.2), users can download the "Unsubtracted Science Scene" or the "PSF Subtraction Source only" as well as perform sub-optimal subtractions as described in the JWST ETC Coronagraphy Strategy page. The image registration, treatment of the noise and background have not changed and are somewhat optimistic.
Currently, the default PSF subtraction strategy is called "optimal." Under this option, the ETC assumes that thermal and dynamical changes do not occur in the optical system between observations of the science target and the reference PSF star exposures. Residual noise, in this case, is driven solely by the shot noise in the wings of the host star and reference PSF source profiles, and SNR for the detection of the faint source depends on the exposure time. Ideally, the reference PSF star would be brighter and will be re-scaled for PSF subtraction. However, as currently implemented in the ETC, the user must set the magnitude and stellar type of the reference star to closely match that of the target star in order to yield a realistic result. Choosing different magnitudes and/or stellar types for the target and reference star will not trigger a warning in the ETC, but the user should be aware that such results might not be accurate.
Example of ETC calculations for high-contrast imaging
The next section provides an example of ETC calculations under the "ideal" assumption.
Tables A–D show the input values and computational results for an ETC study of an analog of the β Pictoris system, comprising a circumstellar disk and a self-luminous planetary companion. These observations use the MIRI four-quadrant phase mask (4QPM) coronagraph that is optimized at λ0 = 11.3 μ.
For PSF subtraction in this example, only referenced differential imaging (RDI) is used, with a PSF reference star at an apparent separation significantly larger than that of the host and companion sources.
In Tables A–C, the column headers identify the suite of ETC inputs for coronagraphic observations, presented here in approximately the same order as the user encounters them on the ETC interface.
The column footers, in italics, assign an ordinal label to each input, to facilitate the descriptions and comments.
This section assumes some user familiarity with the general ETC documentation.
Web users may click on Tables A–D for a larger view.
A1–A2: source identifiers
A3–A4: the spectrum or color of a source, expressed as a spectral type in the Phoenix system or as the effective temperature (Teff) of a blackbody (BB)
A5–A6: the apparent Vega magnitude of the source in K band. Note that, for now, the value of the brightness of the reference star must be identical to the value for the host star. Otherwise the ETC will give unphysical negative SNRs
A7–A8: the shape of each source, point or extended. If extended—referring now to the disk—A8 gives the standard deviations (A+, A–) of an equivalent dual-Gaussian distribution
A9–A11: the X and Y offsets of a source from (0,0), and the rotation-in-place of the source (not meaningful for point sources)
The information on the host and reference stars comes from a catalog or outside research.
The X-Y offsets and PA value are arbitrary and purely notional.
B1: calculation identifier
B2: sources included in the scene
B3–B5: zodiacal light foreground (ignored)
B6–B7: instrument setup (selected coronagraph type and filter)
B7: detector subarray
B8–B10: detector readout pattern and numbers of groups, integrations, and exposures
rotates the position of a point source by theta
- rotates the position and orientation of an extended sources (although generally they are centered at 0,0)
C2: select the reference source, selecting a source identifier from A1. Here, with only the RDI PSF strategy available, we must choose source #3
C3: The PSF subtraction strategy
C4: select which companion source is being observed (disk or planet)
C5: "contrast azimuth" is
- The azimuthal direction along which the contrast vs separation figure is produced
- The azimuthal direction used for the scalar calculation of contrast in the text form ETC report
C6: radius of the virtual photometric aperture, which is centered on the position of the SNR source
C7: "contrast separation" is the radial separation used for the scalar calculation of contrast in the text form ETC report
C8–C9: specify the annular, virtual photometric aperture, which collects stray light from around the SNR source, preparing for its subtraction in post-processing
Discussion
The ETC output plots show no evidence of saturation for these computations.
The results of computations 1 and 2—planet and disk—show a reasonable job of detecting both the planet and disk, with SNR = 38 and SNR = 3, respectively, in 1,422 s exposure time. See Figures 1 and 2 for the two-dimensional images on the detector.
\