FP - Reduces an 8-frame Fabry-Perot observation

Description:

This script reduces a Fabry-Perot observation with UFTI data. It takes an imaging observation comprising eight object frames and a dark frame to make a continuum-subtracted and sky-subtracted, untrimmed mosaic automatically.

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

$$\frac{[(F1 - F2)-(F4 - F3)] + [(F5 - F6)-(F8 - F7)]}{\rm Flat field}$$

where Fn is the bad-pixel masked and dark subtracted frame n. In practice, the flat field is applied to each differenced pair, such as (F4 - F3), when the pair becomes available, rather than waiting until all eight frames have been observed.

Notes:

• 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.

Output Data :

• The integrated mosaic in gf<date>_<group_number>_mos, where <date> is the UT date in yyyymmdd format.

• A mosaic for each cycle of eight in
  gf<date>_<group_number>_mos<cycle_number>, where <cycle_number>
  counts from 0.

• The individual flat-fielded frames in f<date>_<obs_number>_ff.

Configurable Steering Parameters:

NPAIRS = INTEGER

The number of frame pairs to be differenced. It must be a multiple of 2 otherwise 4 is assumed. A value of four or more is assumed to indicate sky subtraction. [4]

NUMBER = INTEGER

The number of spatial jitter positions. For each spatial position there are NPAIRS pairs of frames. A value of 1 also dictates that no jittering has occurred. To make a master mosaic combining spatial positions NUMBER should be at least 3.

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]

USEVAR = LOGICAL

Whether or not to create and propagate variance arrays. [0]

Related Recipes :

SKY_FLAT_FP, FP_JITTER, FP_JITTER_NO_SKY.

Implementation Status:

• The processing engines are from the Starlink packages: CCDPACK, KAPPA, and FIGARO.

• Uses the Starlink NDF format.

• History is recorded within the data files.

• The title of the data is propagated through intermediate files to the mosaic.

• Error propagation is controlled by the USEVAR parameter.