Reduces an 8-frame Fabry-Perot observation FP
The sequence of frames expected in the observations are tabulated below.
Frame | Position | Wavelength
|
1 | On source | On line |
2 | Off source | On line |
3 | Off source | Off line, positive offset |
4 | On source | Off line, positive offset |
5 | On source | On line |
6 | Off source | On line |
7 | Off source | Off line, negative offset |
8 | On source | Off line, negative offset |
It performs a null debiassing, bad-pixel masking, dark subtraction, pairwise frame differencing, flat-field division, integer shifts of origin to register, and mosaicking. The desired result is given by
where is the bad-pixel masked and dark subtracted frame . In practice, the flat field is applied to each differenced pair, such as (), when the pair becomes available, rather than waiting until all eight frames have been observed.
A World Co-ordinate System (WCS) using the AIPS convention is created in the headers should no WCS already exist.
The bad-pixel mask applied is $ORAC_DATA_CAL/bpm.
Each dark-subtracted frame has thresholds applied beyond which pixels are flagged as bad. The lower limit is 5 standard deviations below the mode, but constrained to the range 100 to 1. The upper limit is 1000 above the saturation limit for the detector in the mode used.
You should use SKY_FLAT_FP to make the flat field.
Registration is performed using the telescope offsets transformed to pixels.
There is no resampling, merely integer shifts of origin.
The recipe makes the mosaics by applying offsets in intensity to give the most consistent result amongst the overlapping regions. The mosaic is not trimmed to the dimensions of a single frame, thus the noise will be greater in the peripheral areas having received less exposure time. The mosaic is not normalised by its exposure time (that being the exposure time of a single frame).
For each cycle of eight, the recipe creates a mosaic, which is then added into a master mosaic of improving signal to noise. The exposure time is also summed and stored in the mosaic’s corresponding header. Likewise the end airmass and end UT headers are updated to match that of the last-observed frame contributing to the mosaic.
Intermediate frames are deleted except for the flat-fielded (_ff suffix) frames.
The integrated mosaic in gfdate_group_number_mos, where date is the UT date in yyyymmdd format.
A mosaic for each cycle of eight in
gfdate_group_number_moscycle_number,
where cycle_number
counts from 0.
The individual flat-fielded frames in fdate_obs_number_ff.
If NUMBER is absent, the number of offsets, as given by internal header NOFFSETS, minus one is used. An error state arises if the resulting number of jittered frames is fewer than 3, and a default of 3 is assumed.
If neither NUMBER nor NOFFSETS is defined, 1 is used. [1]