|
|
|
WIRCam Direct Imaging Exposure Time Calculator (DIET)
|
DIET for WIRCam is a direct port from the "CFHT12K/MegaCam DIET calculator."
The basics:
Table of contents:
The calculator

Click on this figure to launch the calculator.
Quick WIRCam photometric performance table
| Filter (click for details) | Y | J | H | Ks | Low OH- 1 | Low OH- 2 | CH4 Off | CH4 On | H2 v=1-0 S(1) | K continuum |
| Point source in median sky brightness - MagAB - Optimal ap. | 22.8 | 22.9 | 22.6 | 22.5 | 20.9 | 21.0 | 22.0 | 22.0 | 21.1 | 20.9 |
| Field galaxy in median sky brightness - MagAB - 2.2" ap. | 22.1 | 22.3 | 21.9 | 21.8 | 20.2 | 20.3 | 21.3 | 21.3 | 20.4 | 20.2 |
| Conversion from AB to Vega magnitude system (mag) | -0.66 | -0.96 | -1.40 | -1.99 | -0.69 | -0.87 | -1.35 | -1.47 | -1.97 | -2.08 |
This table summarizes the expected camera performance for a 10 sigma
detection in a 1 hour exposure under 0.7 arcsecond seeing with 1.5 airmass.
Those are AB instrumental magnitudes. These numbers are now based on standard
star measurements done during the commissionning. Read the important note below
if you observe an extended object or a field that requires a sky-target-sky
observing strategy: these estimates then need to be offset to significantly
brighter magnitudes.
Important note on sky subtraction and limiting
magnitudes with WIRCam. If you need to use the nodding (sky-target-sky)
observing strategy (i.e. if your target is larger than 10 arcmin) then your
request should be for twice the output of the calculator, plus overheads. We
urge you to read further. The quality of construction and
subtraction of the sky frames strongly depends on your chosen observing
strategy and target. The DIET estimates DO NOT take the sky subtraction
into account. In the worst case scenario, if the constructed sky frames have
the same SNR as the individual science frames, the limiting magnitude drops by
a factor sqrt(2) or ~0.4 magnitudes compared to the DIET estimates. And
remember that DIET's estimates is based on ON TARGET exposure times. So if you
plan to use a sky-target-sky strategy, in addition to the DIET exposure time,
your time proposal needs to typically double to account for the sky time.
Sky-Target-Sky strategies spend typically half their time on the target such
that another factor sqrt(2) is lost (for a given telescope allocation). So, if
a sky-target-sky strategy has to be used and if only one image is used in the
construction of the sky frame, then the limiting magnitude is ~0.8 magnitudes
brighter that DIET's estimates. In a case where sky frames can be built from a
larger set of science frames and the subtraction be applied to each individual
science frame, then the DIET estimates do not need to be offset.
Caution! it was found by three different teams that the Vega
to AB conversions were off by up to 0.2 mag in Ks. The numbers given here will
need to be revised. Meanwhile, here are the numbers other teams are
obtaining:
Emmanuel Bertin at Terapix:
Y Wircam: 0.583
J Wircam: 0.899
H Wircam: 1.335
Ks Wircam: 1.802
Stephane Arnouts at CFHT:
# FILTER NAME IDENT Lbda_mean Lbeff(Vega) FWHM AB-cor
Y engineering cfh8001.pb 1 0.9982 0.9945 0.1016 0.577
Y cfh8002.pb 2 1.0259 1.0221 0.1106 0.609
J (in use) cfh8101.pb 3 1.2545 1.2481 0.1588 0.924
Low OH2 cfh8102.pb 4 1.1892 1.1892 0.0113 0.836
J (spare) cfh8103.pb 5 1.2540 1.2478 0.1572 0.923
Low OH1 cfh8104.pb 6 1.0629 1.0628 0.0101 0.655
H (in use) cfh8201.pb 7 1.6310 1.6161 0.2910 1.352
H (spare) cfh8202.pb 8 1.6369 1.6222 0.2895 1.358
CH4 On cfh8203.pb 9 1.6919 1.6897 0.1062 1.418
CH4 Off cfh8204.pb 10 1.5892 1.5875 0.0950 1.331
Ks (spare) cfh8301.pb 11 2.1519 2.1358 0.3280 1.826
Ks (in use) cfh8302.pb 12 2.1497 2.1338 0.3270 1.824
K Cont cfh8303.pb 13 2.2246 2.2244 0.0327 1.898
H2 cfh8304.pb 14 2.1300 2.1298 0.0294 1.816
Br Gamma cfh8305.pb 15 2.1720 2.1720 0.0295 1.894
WIRCam Filter Page (click for details)
Mark Dickinson - COSMOS group:
J 0.937
H 1.365
K 1.846
|
1 How DIET works
|
DIET (Version 2.0) is a calculator allowing the observer to compute the
exposure time (in seconds) required to reach a given signal-to-noise ratio in
various observing conditions (source type, magnitude, filter, seeing, sky
background, airmass, and atmospheric transmission). Encouraged by the accurate
results of DIET for Megacam, DIET for WIRCam is simply a copy.
DIET's current computations are now based on actual on-sky measurements. The
zero points have been determined and the sky median background has been
measured over the full 2005B semester for all filters except LowOH2 and H
continuum.
The interactive graphical interface allows the user to experiment with
some custom parameters. This iterative process can be time consuming
for a single set of parameters: it is recommended to use the keyword
"Range" for both parameters "Seeing" and "Sky" to obtain small tables
exploring the domain of input conditions.
The magnitude system in DIET for WIRCam is AB. The following link
(courtesy of D. Patton) provides information on the various magnitude
systems and ways to go from one to another:
http://www.astro.utoronto.ca/~patton/astro/mags.html.
|
1.1 Measured zero-points and sky brightnesses
|
The following table gives the measured zero-points. It also gives the median
measured sky flux (adu/sec/pixel) through each of the 10 WIRCam filters during
semester 2005B and the corresponding sky brigthness (ABmag/arcsec2). All along,
a detector gain of 2.5e-/adu is assumed (it was the gain from Nov 2005 to Aug
2006 - since then the gain is 3.7+-0.1e-/adu, but that does not affect the
following measurements). The conversion from the AB to the Vega magnitude
system is also given (for example, J_Vega = J_AB - 0.96).
WIRCam Zero-Point and Sky Brightness in AB Magnitudes
| Filter | Y | J | H | Ks | Low OH- 1 | Low OH- 2 | CH4 On | CH4 Off | H2 v=1-0 S(1) | K Continuum |
| Zero Point (AB mag) | 25.31 | 25.98 | 26.58 | 26.42 | 22.40 | 22.48 | 25.31 | 25.30 | 23.75 | 23.71 |
| Sky Brightness (AB mag/arcsec2) | 17.3 | 16.7 | 15.5 | 16.2 | 18.8 | 18.5 | 15.6 | 15.3 | 15.9 | 15.6 |
| Detected Sky Flux (adu/sec/pix) | 60 | 200 | 1000 | 465 | 1.0 | ~1.5 | 275 | 380 | 50 | 90 |
| AB to Vega Magnitudes | -0.66 | -0.96 | -1.40 | -1.99 | -0.69 | -0.87 | -1.47 | -1.35 | -1.97 | -2.08 |
The sky flux and brightness are the the medians for semester 2005B. They
correspond to GREY conditions in DIET.
The sky variability in the near infrared is important. As a canonical rule, the
level varies by 10% over a 10 minute timescale. The H band is notoriously
variable by +/-0.5 mag. In J and Ks, it +/-0.3 mag.
|
2. Information on how DIET computes exposure times vs. magnitudes and SNRs
|
The following paragraphs propose a tutorial on the magnitude
and signal to noise calculation schemes used in DIET.
2.1 Object types
There are three classes of objects considered in DIET:
- Point sources: stars & QSOs
- Galaxies: distant compact to sub-compact sources - seeing dominated profile
- Nearby galaxy: profile dominated by the object itself
- Extended source: uniform illumination over square arcsecond scales
2.2 Objects profiles
As provided in Iraf's Imexamine function, a Moffat function is used to fit
all the profiles:
I = Ic (1 + (r / alpha)^2)^(-beta) [1]
where Ic is the peak value, r is the radius, and the parameters
are alpha, and beta. The alpha value is equal to half the FWHM
of the object (FWHM = seeing for seeing dominated profiles, and
FWHM = object width at half maximum for large nearby galaxies),
and beta defines the profile type, say how much energy is distributed
in the wings of the profile. Based on a large sample, the median
value for beta on WIRCam images, using Iraf's Imexamine, is:
- Point sources: beta = 3.8
- Galaxies: beta = 2.5
- Nearby galaxies: beta = 1.8
- Extended source: non applicable
2.3 Exposure time and SNR computation
Equation [1] can be integrated over the profile and the flux within
an aperture of r=R is: (Mathematica integration)
F(R) = (Pi Ic alpha^2)/(1 - beta) ( (1 + (R / alpha)^2)^(-beta + 1) - 1) [2]
If R is infinite, the flux represents the total flux received from
the object:
F(total) = (Pi Ic alpha^2)/(beta - 1) [3]
Equation [2] becomes:
F(R) = F(total) ( 1 - (1 + (R / alpha)^2)^(-beta + 1)) [4]
And F(total) can be computed simply from the magnitude equation,
considering flux in electrons per second (Fe):
MagAB = Zero - k(a - 1) - 2.5log10 (Fe) [5]
With Zero the zero photometric point in electron per second, k
the airmass term, a the airmass and Fe the total flux in electrons
generated per second by the object of magnitude MagAB.
Fe = 10^((Zero - MagAB - k(a - 1)) / 2.5) [6]
Integrated over the exposure time Texp, we have:
F(total) = Fe Texp [7]
Now let us compute the signal to noise ratio over the aperture
of radius R (flux considered in electrons):
SNR = F(R) / ( sqrt( F(R) + n Pi R^2 (Se Texp + noiseccd^2)) ) [8]
Where Se is the flux in electrons per second and per pixel coming
from the sky background and noiseccd is the CCD readout noise.
Pi R^2 represents the circular area of integration of the profile.
Since the sky background to be subtracted from underneath the
object is more often than not totally noise-free due
to some modelling problem or crowding, the factor n is introduced.
n can range from 1 for a perfect sky subtraction to 2 for a 1-1
pixel type subtraction. n could be higher than 2 actually in
some crowded fields when some earby object flux is subtracted to
the object itself. In DIET, a value of n=1.5 is used.
Now solving the quadratic equation from [8] for Texp, one gets:
Texp = ( -b + sqrt (b^2 - 4 a c) ) / 2 a [9]
a = Fe^2, b = -SNR^2 (Fe + n Pi R^2 Se), c = -SNR^2 n Pi R^2 noiseccd^2
2.4 Optimal radius of flux integration
Maximizing the SNR over R for point spread functions (beta = 3.8),
the relationship
Ropt = 1.45 * alpha is obtained. This results in measuring 96%
of the flux within the aperture. To integrate 96% of the flux for
the galaxies (beta = 2.5) which carry more weight in their wings, the relation
becomes Ropt = 2.80 * alpha (which corresponds to a 2.2 arcsec
aperture for a 0.8 arcsec seeing). For nearby galaxies (beta = 1.8) which cover
around 10 arcsec in diameter (and for which seeing plays little effect
for SNR computation over the profile), the relation to get 96% of
the flux within the diameter of integration leads to a
generous Ropt = 13 * alpha (which corresponds to a 11 arcsec aperture
for a 0.8 arcsec seeing).
By picking a radius for all three cases that lead to the integration of
96% of the light, the difference in magnitude between the total magnitude
and the magnitude within the aperture is 0.044 mag. This will be of
importance when comparing DIET results to real measurement by SExtractor
which will have to be corrected for the 0.044 mag (subtraction).
For extended sources, the SNR is computed over a single pixel though
the magnitudes have to be provided per square arcsecond.
2.5 Direct relation between the SNR and the error on the magnitude
Let us consider a simplified expression of the magnitude vs. flux
measurement: m = -2.5 log(f). The error Em on m can be written Em = |dm/df| Ef
where Ef is the error on f. Since |dm/df| = (1 / f)(2.5 / ln(10)), one
gets Em = 0.92 (1 / (f / Ef)). Since SNR = f / Ef, the final relation is:
SNR = 0.92 / Em [10]
Where SNR is the signal to noise ratio on the object and Em the error on
the magnitude on that object. For example, a error in magnitude of 0.1
mag gives a SNR of ~9, say a ~10% error.
When using a source extraction software like
SExtractor (the most commonly used software today for large images) the error on
the magnitude is estimated within the defined aperture. SExtractor derives the variance of
the background (should it be sky background and/or detector read noise) regardless of its level.
Hence it assumes the image is not affected by a convolution and/or resampling. In that regard,
all the testing of DIET was done directly on raw data, knowing that further pre-processing and
processing should increase the detection performance and the quality of the photometry.
2.6 Relation between SNR and photometry quality
- SNR = 3 : Detection - Photometry error = 33%
- SNR = 7 : Fair detection - Photometry error = 15%
- SNR = 15 : Good detection - Photometry error = 7%
- SNR = 25 : Quality photometry - Photometry error = 4%
- SNR = 100 : High quality photometry - Photometry error = 1%
|
3. Examples of WIRCam photometric performance and image samples
3.1 Camera and site characteristics used in DIET
| Magnitude system | AB |
| Instrumental magnitude equation | MagAB = Zp[e-/sec] - 2.5log(Flux[ADU]) -2.5log(Gain[e-/ADU]) + 2.5log(ExpTime) - k(airmass - 1) |
| Gain (e- / ADU) | 2.5 |
| Filters | Y, J, H, Ks |
| Zero points (e-/sec) at 1 airmass (Zp) | 25.31, 25.98, 26.58, 26.42 |
| Airmass term | 0.02, 0.05, 0.03, 0.05 (expected) |
| Median sky brightness @ zenith (Mag/Arcsec^2)) | 17.3, 16.7, 15.5, 16.2 |
| Sky brightness airmass dependency offset (e-/s/pixel/airmass) | Estimates are: 20, 40, 180, 170 |
| Transition from AB to Vega magnitude system (mag) | -0.66, -0.96, -1.40, -1.99 |
| Pixel scale (arcsec / pixel) | 0.3 |
| Read noise (e- / pixel) | 30 |
3.2 Example of performance on point sources
Here is the output from DIET with the given input parameters.
J filter - MagAB=23.5 - SNR=5
Filter=J Sky=GREY Bin=1pix
Type=point source SNR=5.0 Trans=100%
Mag=23.5 SNR ap=Optimal
______________________________________________________________________
Seeing | Airmass
| 1.0 1.5
| |------------------------------------------------------------
0.4 | 761 823
0.5 | 1265 1367
0.6 | 1769 1912
0.7 | 2273 2457
0.8 | 3028 3274
0.9 | 3784 4092
1.0 | 4540 4909
1.1 | 5548 5999
1.2 | 6556 7088
Ks filter - MagAB=22.5 - SNR=7
Filter=KS Sky=GREY Bin=1pix
Type=point source SNR=7.0 Trans=100%
Mag=22.5 SNR ap=Optimal
______________________________________________________________________
Seeing | Airmass
| 1.0 1.5
| |------------------------------------------------------------
0.4 | 505 544
0.5 | 839 905
0.6 | 1174 1267
0.7 | 1509 1628
0.8 | 2011 2170
0.9 | 2513 2711
1.0 | 3015 3253
1.1 | 3684 3975
1.2 | 4354 4698
3.3 Example of performance on field galaxies
Here is the output from DIET with the given input parameters.
Ks filter - MagAB=22.0 - SNR=5
Filter=KS Sky=GREY Bin=1pix
Type=galaxy SNR=5.0 Trans=100%
Mag=22.0 SNR ap=Optimal
______________________________________________________________________
Seeing | Airmass
| 1.0 1.5
| |------------------------------------------------------------
0.4 | 374 404
0.5 | 578 623
0.6 | 815 879
0.7 | 1120 1208
0.8 | 1459 1574
0.9 | 1831 1976
1.0 | 2238 2415
1.1 | 2712 2927
1.2 | 3221 3475
H filter - MagAB=20.0 - SNR=20
Filter=H Sky=GREY Bin=1pix
Type=galaxy SNR=20.0 Trans=100%
Mag=20.0 SNR ap=Optimal
______________________________________________________________________
Seeing | Airmass
| 1.0 1.5
| |------------------------------------------------------------
0.4 | 138 146
0.5 | 212 225
0.6 | 299 316
0.7 | 410 434
0.8 | 534 565
0.9 | 670 709
1.0 | 818 866
1.1 | 991 1050
1.2 | 1177 1246
3.4 Example for units conversion from Jansky to AB mag for time proposal purposes
By definition, the zero point is the magnitude of a source which yields 1
e-/sec (detected). The AB system uses for mag=0 the constant flux of 3720
Jansky (i.e. 1 Jansky = 10^-26 J/s/m2/Hz). So, if you know your source flux in
Jansky, you can compare it to the zero point magnitude and deduce the flux of
detected photons that you should expect for your source. For example:
In H2, a zero point of 23.75 means that a flux of 3720 Jy / 2.512^(23.75) =
1.2x10-6 Jy produces a flux of 1 detected electron/sec. Say your source flux at
the wavelength of H2 is 10^-3 Jy, then your source has a magnitude of magAB =
2.5*log10(3720/10^-3) = 16.43. And a magAB=16.43 in H2 produces a flux of
2.512^(23.75-16.43) = 847 detected photons/sec.
If you integrate for 100 sec, you should expect 85000 photons on your source.
|
4 How to optimally fragment the total exposure time
|
To efficiently remove the cosmic rays and cosmetic defaults of the mosaic
(gaps between the detectors, bad columns), a minimum of 4 dithered exposures
per field is required. To obtain a reasonably uniform SNR across the gap
areas, we recommand using at least 5 positions whenever a contiguous coverage
is desired. The number of positions and their placement, however, will
often be driven higher by sky-subtraction considerations: the sky background
in the IR is both much higher and more variable (specially in the first few
hours after sunset) than in the optical. The sky signal is usually estimated
from the neigbouring positions in the dither pattern, which then need to
be spaced by more than the size of the largest structure of interest in
the image.
Each exposure needs to be in the sky photon noise regime, so that their
coaddition will produces the expected signal to noise ratio for the total
exposure time. WIRCam has a low readout noise (by IR camera standards), which
in the broad band filters is very quickly dominated by the sky photon noise:
except (marginally) at Y, the minimum exposure time of 5s actually prevents
obtaining a sky exposure that has significant readout noise. This is, on
the other hand, a significant consideration for the narrow band filters,
and in particular for the low OH filters in the J band.
The saturation level is usually a more important consideration: exposing
too long will indeed save a couple of minutes by skipping some
readouts, but will result in high sky background and too many
objects (and even possibly the sky) reaching saturation. For the broad band
filters, the longest reasonable exposure time (e.g. 20s at Ks) would often
result in very large overheads if the telescope would move between every
exposures (we estimate that each offset will initially take 60s, to be
later optimised down to ~20s). Multiple exposures are then taken at a
given telescope position, before moving to the next position in the pattern.
The maximum number of
exposures at each position, which must be a multiple of 4 if microstepping is
used, is then set by the timescale of the background fluctuations.
We estimate that the following exposure times and times per position
will be good compromises to achieve low overheads while keeping
the signal in a reasonable range:
- Y band: 2 mn and 8 mn
- J band: 45s and 3 mn
- H band: 10s and 1.5 mn
- Ks band: 20s and 1.5 mn
- CH4on band: 30s and 2 mn
- CH4off band: 30s and 2 mn
- LowOH1 band: 10 mn and 10 mn
- H2 band: 3 mn and 12 mn
|
4.1 Maximum Exposure Times
|
The following table gives the maximum exposure times allowed for WIRCam
designed so that the sky background uses no more than half the detector
potential well. The threshold adopted is based on the 75% percentile of the
first year of WIRCam observations.
Maximum WIRCam Exposure Times (per single slice)
| Filter | Max. etime |
| Y | 150s |
| J | 60s |
| H | 15s |
| Ks | 25s |
| Low OH1 | 3000s |
| Low OH2 | 3000s |
| CH4 On | 50s |
| CH4 Off | 50s |
| H2 | 200s |
| K Cont. | 200s |
| Br. Gamma | 200s |
|
4.2 Overheads
|
The following table gives the expected overheads for WIRCam.
Expected WIRCam Overheads
| Action | Overhead | Charged to CFHT or PI? | Comment |
| Readout Overhead | 10 sec | PI | With final 32-amp controller |
| Micro-Dithering | ~0.1s | CFHT | |
| Small Telescope Offsets | 10s | CFHT | |
| Large Sky Offset | 60 sec | PI | Charged twice: going to Sky and returning from Sky |
| Filter Wheel Change | 20 sec | CFHT | Recommend completing one filter before going to next |
| Focusing | minutes | CFHT | As part of QSO |
| Standard stars | minutes | CFHT | As part of QSO unless needed photometry is better than 3-5% |
|
|
|