- 7.5.1 Processing overview for chopped observations
- 7.5.2 Setting up data for generic pattern construction
- 7.5.3 Generic pattern construction
- 7.5.4 Determination of source signal
- 7.5.5 Correction for chopped signal loss: Sky measurement
- 7.5.6 Correction for chopped signal loss: FCS measurement

7.5 Signal Processing: Chopped PHT-P and PHT-C

** Detailed description:** see also
Section 5.2.3

Analysis of many chopped observations performed with PHT-P and PHT-C (involving AOTs PHT03, PHT04, and PHT22) has shown that the `conventional' processing method had to be changed drastically. The main driver for this is that chopped measurements do not give stabilised signals which cause significant losses on the true difference signal. Most of the processing steps which were commonly shared between the staring and chopped observations became obsolete for the chopped observations. Instead, a separate SPD level processing was used.

The present signal derivation relies on the analysis of signals from pairs of consecutive readouts rather than signals per ramp. This gives better statistics of the signals per chopper plateau, since in many chopped measurements each chopper plateau covers only a few (typically 4) ramps.

To increase further the robustness in determining the difference signal, the repeated pattern of off-source and on-source chopper plateaux is converted into a `generic pattern'. The generic pattern consists of only 1 off- and 1 on-source plateau and is generated using an outlier resistant averaging of all plateaux. The shape of the generic pattern determines the correction factors with regard to stabilised staring measurements.

The FCS measurement is used to determine the responsivity of a given observation. The chopped FCS measurements are treated in the same way as the chopped sky measurements, namely for both types of measurement a generic pattern is constructed.

Using a symmetry assumption, the corrected signal level of FCS1 is determined from the measured signal levels of FCS1 and FCS2 and applying the same signal loss correction as for the sky measurements. Following the creation of a generic pattern, the DERIVE_SPD processing proceeds as for staring mode observations from reset interval correction to flux calibration, which uses the corrected FCS1 signal of the chopped FCS measurement.

After reading the ERD the following processing steps are performed for Derive_SPD:

- ramp linearisation,
- generic pattern creation,
- reset interval correction,
- dark signal subtraction,
- signal linearisation,
- determination of the difference signal from pattern,
- determination of the difference signal corrected for chopped signal losses,
- flux calibration.

Steps 2, 6, and 7 are explained in more detail in the following sections.

None

** Detailed description:** none

A * chopper unit* comprises two consecutive chopper plateaux, in
particular either plateaux on the background and source or on FCS2 and
FCS1. The number of chopper units and * chopper
dwell time* of a measurement
in a given chopper mode is determined according to:

(7.21) | |||

(7.22) |

and

(7.23) |

where:

- in s is the measurement time,
- is the number of ramps per chopper plateau,
- in s is the reset interval time,
- is the index representing the dwell time in the calibration files.

In rectangular mode one full chopper cycle corresponds to one chopper unit:

(7.24) |

in triangular mode one chopper cycle corresponds to two consecutive chopper units:

(7.25) | |||

(7.26) |

and in sawtooth mode, one chopper cycle ( ) also corresponds to two consecutive chopper units, but the last plateau of the second unit is regarded as missing and filled up with zeros:

(7.27) | |||

(7.28) |

The time interval of a * logical ramp* is defined as 1/8 of the
duration of a chopper unit and can be derived from:

(7.29) |

where:

- is the total number of readouts, which includes both destructive and non-destructive readouts;
- NDR is the number of non-destructive readouts;
- is the actual number of ramps per chopper plateau in the telemetry and is derived from where is the commanded number of ramps per chopper plateau and DAT_RED is the applied data reduction factor.

Before the signal derivation, readouts which can cause systematic deviations are `cleaned' or flagged `bad':

- destructive readouts, see Section 7.2.1;
- a fixed fraction of readouts at the beginning of a ramp, the same fraction is used as for staring measurements, see Section 7.2.2;
- saturated readouts, determined in the same way as described in Section 7.2.7;
- If for a given measurement in case of rectangular mode or in case of triangular or sawtooth modes, then the first 25% of the whole measurement is discarded;
- In case of a gap in the data stream due to e.g. a telemetry drop, the entire chopper unit in which the gap occurs is discarded. The start of the next unit is determined from the start time (in ITK) of the last chopper unit before the gap.

For all available pairs of consecutive readouts the * difference
signal* is computed:

(7.30) |

where and in V are the voltages and and are the times of consecutive readouts. Note that the time difference ( ) is always the same for a given measurement. After calculation, each is stored in an array in which each element represents a logical ramp in the measurement. In total there are elements in this array. In case of sawtooth chopper mode, the array elements belonging to the logical ramps of the second plateau of the even chopper units (2, 4, 6, ...) are set to zero.

None

7.5.3 Generic pattern construction

** Detailed description:** Section 5.2.3 for
an overview

The generic pattern is constructed by stacking the chopper units onto one single unit. However, long term transients or detector drifts which usually last longer than several consecutive chopper units introduce an unwanted noise in the generic pattern. Therefore, before stacking the chopper units, the measurement is corrected for a drift by normalising the chopper units (=1, 2, 3,...,):

* Rectangular chop mode*:

compute the median of per chopper unit :

compute

store in an array with elements

determine a scaling factor
from the average
over elements to .

* Triangular chop mode*:

compute the median of in two consecutive chopper
units : , ( is odd)

compute

store in an array with elements

determine a scaling factor
from the average
over elements to .

* Sawtooth chop mode*:

compute the median of of the first of every
two consecutive units : , ( is odd)

compute

store in an array with elements

determine a scaling factor
from the average
over elements to .

Before the generic pattern is created two intermediate patterns are
obtained by stacking of the odd (pattern 1) and even (pattern 2)
chopper units. The 8 average signals plus their associated uncertainties
which correspond to the 8 logical ramps in each pattern are obtained by
computing an outlier resistant mean of the values of for each
logical ramp.
For triangular mode the two patterns correspond to the different background
positions. For the sawtooth mode, of the logical ramps
5 to 8 of the even chopper units (pattern 2) is always zero.

The * generic pattern* for rectangular and triangular chop mode
is determined from:

(7.31) | |||

(7.32) | |||

(7.33) | |||

(7.34) |

similarly, for sawtooth mode:

(7.35) | |||

(7.36) |

where:

- =1,,8 is the logical ramp index,
- , (in V/s) are the outlier resistant means of pattern 1 and pattern 2,
- , (in V/s) are the associated uncertainties,
- and (in V/s) are intermediate uncertainties,
- (in V/s) is the final signal uncertainty for each logical ramp in the generic pattern.

The 8 signals in the generic pattern are subsequently subject to reset interval correction (Section 7.3.1), dark current subtraction (Section 7.3.2), and signal linearisation (Section 7.3.4).

None

7.5.4 Determination of source signal

** Detailed description:** none

Based on the generic pattern, the signal of the source is determined. The calibration analysis has shown that the appearance of the pattern depends on several parameters such as the detector used, the signal difference, the mean signal level, as well as the chopper dwell time. The determination of the signal depends on the detector used:

(7.37) | |||

(7.38) | |||

(7.39) | |||

(7.40) | |||

(7.41) |

To determine the median, the middle two values of the 4 elements are averaged. Depending on the operation, the signal uncertainties are determined via:

I. In case of average value:

(7.42) | |||

(7.43) |

II. For the median value, where and are the indices of
the two values averaged to obtain the median:

(7.44) | |||

(7.45) |

III. In case of minimum-maximum value, where and are the
indices of the maximum and minimum, respectively:

(7.46) | |||

(7.47) |

Finally, the source signal uncertainty is derived from

(7.48) |

None

7.5.5 Correction for chopped signal loss: Sky measurement

** Detailed description:** Section 5.2.3

The source signal is corrected for signal loss by the chopper modulation by means of a transformation:

(7.49) |

where in V/s is the corrected value. The correction function depends besides also on the detector pixel and chopper dwell time . The function is implemented via look-up tables. Each table gives the correction values for a given detector. Within each table the correction values are ordered by dwell time and detector pixel.

The zero point of takes care of the vignetting or chopper offset correction (see Section 4.5.3). As a consequence, separate vignetting correction tables are not required.

Assuming that the signal loss is symmetric, i.e. the loss from the on-source signal is gained from the off-source signal, it is possible to correct the on- and off-source signals:

(7.50) | |||

(7.51) |

then the corrected signals can be determined from:

(7.52) | |||

(7.53) |

In practice, the signal loss is asymmetric for detectors P3, C100, and C200. For these detectors the loss on the on-source signal is not equivalent to the gain from the off-source signal. Based on the previous results, the asymmetry is included by using an empirical asymmetry factor :

(7.54) | |||

(7.55) |

where are empirical constants with
0.3558,
, 2.579,
0.6472,
,
. Finally, the on- and off-signals corrected
for the signal asymmetry are computed for
and:

If
:

(7.56) | |||

(7.57) |

If
:

(7.58) | |||

(7.59) |

In case and for detectors P1 and P2, and .

Cal-G files PP1CHOPSIG, PP2CHOPSIG, PP3CHOPSIG, PC1CHOPSIG, and PC2CHOPSIG contain the look-up tables for . The description of these files is given in Section 14.10.

7.5.6 Correction for chopped signal loss: FCS measurement

** Detailed description:** none

For the FCS measurements only rectangular chop mode is used. The generic pattern is constructed in the same way as for the chopped sky measurements for which on-source corresponds to beam deflection to FCS1 and off-source corresponds to FCS2. Adopting the same notation as before but with `on'=FCS1, `off'=FCS2, the same signal loss corrections as presented in Section 7.5.5 apply.

Cal-G files PP1CHOPSIG, PP2CHOPSIG, PP3CHOPSIG, PC1CHOPSIG, and PC2CHOPSIG contain the look-up tables for . The description of these files is given in Section 14.10.