ZTF¶
Zwicky Transient Facility (ZTF) is a ground-based telescope run by Caltech at Palomar Observatory. It was designed for wide-field time-domain surveys and often performs Target of Opportunity (ToO) follow-up of gravitational wave events
- Field of View:
7°.50 N-S x 7°.32 E-W; 47.7 deg \(^2\) light sensitive area
- Filters:
r-, g-, i-
- Location:
Palomar Observatory
- Readout Overhead:
8.2 s
- Maximum Slew Velocity:
2.5° s-1
- Maximum Slew Acceleration:
0.4° s-2
Todo
Other important values folks think I should add?
Objective¶
ZTF receives a gravitational wave alert, which has an event time and a 3-D probability density skymap \(w\) in right ascension, declination and luminosity distance. The skymap, along with a a model lightcurve \(L\) that describes the expected time dependence of the source’s flux and a time interval \(t\) to \(t + \Delta T\) to observe during are used to produce a schedule.
This schedule is produced by choosing a subset \(C\) of the set of 1778 fixed reference images for ZTF, denoted \(P\) and arranging them into a schedule \(S\). Restrictions \(S\) must fulfil to be a valid schedule are laid out in the Constraints section.
The objective when creating this schedule is to maximize the sum of the probability density observed since it correlates with detection probability. For a pixel to be marked as observed, an exposure long enough to have sufficient signal-to-noise ratio (≥ 5) needs to be made. The length of the exposure needed primarily depends on \(L\) and how long after the event the exposure is taken.
Todo
Use the same CCD S/N equation for PSF photometry as Dorado? It seems like this is a larger question for the group as a whole
Decision Variables¶
The decision variables consist of:
The number of observations
For each observation:
Its pointing
Its start time
Its exposure time
For each HEALPix pixel on the sky:
Whether it is observed
Constraints¶
ZTF observations are subject to the following constraints:
The telescope must
observe fields with an airmass below 2.5
observe fields greater than 20 degrees from the center of the moon
observe fields not in ZTF pointing limits
|HA| < 5.95 hours
Not HA < -17.6 deg and Dec < -22 deg
Not west of HA -17.6 deg, Dec < -45 deg
Not |HA| > 3 deg and Dec < -46
Dec \(\leq\) 87.5
The telescope can only make observations after the sun is 18° below the horizon
Observations can only be made during the provided time interval
All fields in \(C\) must be observed twice, once in r- and once in g- bands with at least 30 minutes of time between.
For a HEALPix pixel to be marked as observed, an exposure long enough to have SNR ≥ 5 of it must be taken.
Pseudocode¶
from astropy import units as u
from ligo.skymap.io import read_sky_map
import m4opt
observer = m4opt.Observer.on_ground('path/to/ztf.tle')
skymap, data = read_sky_map('path/to/skymap.fits')
start_time = data['gpstime']
m4opt.objective.ToOObjective(skymap).add(observer)
m4opt.constraints.SunElevationConstraint(-18 * u.deg).add(observer)
m4opt.constraints.MoonSeparationConstraint(20 * u.deg).add(observer)
m4opt.constraints.AirmassConstraint(2.5).add(observer)
m4opt.constraints.TimeConstraint(
start=start_time, end=start_time + 24 * u.hour).add(observer)
m4opt.constraints.CadenceConstraint(
count=2, time=30 * u.minute, filters=['r', 'g']).add(observer)
m4opt.constraints.SNRConstraint(
lightcurve='path/to/lightcurve.dat', SNR=5).add(observer)
'''
default pointing constraints automatically stored in ztf file, and
no additional ones are to be added
'''
...