Description:
This routine analyses the jiggle map made by a bolometer
during a PHOTOM observation. If status is good on entry the routine checks that the data array holds
data for the specified bolometer. If it does not an error message will be output and the bad
status returned. All being well the routine will analyse the data by the method specified in
PHOTOM_ANALYSIS:-
AVERAGE - the ‘result’ will be the mean of the jiggle points. The variance is the sum of the variances
of the averaged data divided by the number of them squared. The result will be flat-fielded and
reported to the user. If no valid data were available for the calculation of the mean then a warning
message will be output but good status returned.
PARABOLA - the routine will call SCULIB_FIT_2D_PARABOLA to fit a 2-d parabola to the data and the
‘result’ will be the peak of the fitted curve. If the fit is good the results will be flat-fielded and reported
to the user. Otherwise the result quality is set bad but good status will be returned.
Invocation
CALL
SCULIB_ANALYSE_PHOTOM_JIGGLE (PHOTOM_ANALYSIS, BOL, N_BOLS, J_COUNT, JDATA, VARIANCE,
QUALITY, RESULT_D, RESULT_V, RESULT_Q, BADBIT, STATUS)
Arguments
PHOTOM_ANALYSIS =
CHARACTER()
(Given)
the method of analysis to be applied to the data
BOL = INTEGER (Given)
the index of the
bolometer whose data is to be analysed
N_BOLS = INTEGER (Given)
the number of bolometers
being measured
J_COUNT = INTEGER (Given)
the number of jiggles in the pattern
JIGGLE_X
= REAL (Given)
x jiggle offsets
JIGGLE_Y = REAL (Given)
y jiggle offsets
JDATA
(N_BOLS, J_COUNT) = REAL (Given)
the measured data
VARIANCE (N_BOLS, J_COUNT)
= REAL (Given)
the variance
QUALITY (N_BOLS, J_COUNT) = BYTE (Given)
the
quality
RESULT_D = REAL (Returned)
the sum of the input data
RESULT_V = REAL
(Returned)
the sum of the input variances
RESULT_Q = BYTE (Returned)
the quality on
the sum
A0 = REAL (Returned)
fitted parabola parameter
A1 = REAL (Returned)
fitted parabola parameter
X0 = REAL (Returned)
x offset of peak of fitted parabola
Y0 = REAL (Returned)
y offset of peak of fitted parabola
BADBIT = BYTE (Given)
bad bit mask
STATUS = INTEGER (Given and returned)
global status
Copyright
Copyright ©1993-1999 Particle Physics and Astronomy Research Council. All Rights Reserved.