NAME
pcadeadspect2 - Apply PCA dead-time correction to Standard2 spectra
USAGE
pcadeadspect infile phafile outfile gtiorfile gtiandfile
timecol timemin timemax timeint
DESCRIPTION
pcadeadspect2 applies a dead-time correction to a PCA Standard2
spectrum file.
It is assumed that the user has already run pcadeadcalc2 in order
to compute the dead-time related quantities for the observation,
and it is also assumed that the user has already extracted a spectrum
for the times and detectors they are interested in.
pcadeadspect2 is used to adjust the exposure time of the extracted
spectrum so that it reflects the true live time of the instrument.
pcadeadspect2 operates on both the "ONE" and "MANY" products
of saextrct (i.e. Type I and Type II PHA files).
IMPORTANT NOTE: The user must supply to pcadeadspect2 the
*identical* time filtering GTI files which were used to make the
spectrum. If this is not done, then the calculated exposures will
be less accurate than the original values!
Other selections like which PCUs were extracted, and
the type of spectrum file, are determined automatically.
It is not necessary to supply energy selection settings because
dead-time does not depend on pulse-height.
PARAMETERS
- infile [string - file name or @filename]
-
Name of a FITS file or @filename.txt file name list, which
is identical to the infile passed to 'saextrct' to produce
the input spectrum
- phafile [string - file name]
-
Name of Standard2 spectrum produced by 'saextrct' using the
above 'infile' setting.
- outfile [string - output file name]
-
Name of revised Standard2 spectrum file, whose exposure has
been corrected for dead time.
- gtiorfile = "APPLY" [string]
-
Name of GTI "OR" time filtering, which is identical to the
gtiorfile passed to 'saextrct' to produce the input spectrum.
- gtiandfile = "-" [string]
-
Name of GTI "AND" time filtering, which is identical to the
gtiandfile passed to 'saextrct' to produce the input spectrum.
- (gticols = "START STOP") [string]
-
Name of GTI start/stop column names, which is identical to the
gticols passed to 'saextrct' to produce the input spectrum.
- timecol = "TIME" [string]
-
Name of TIME column used for extraction, which is identical to
the timecol passed to 'saextrct' to produce the input spectrum.
- timemin = INDEF [double precision real]
-
Minimum start time (MET) used for accumulation, which is
identical to the timemin passed to 'saextrct' to produce the
input spectrum.
- timemax = INDEF [double precision real]
-
Maximum start time (MET) used for accumulation, which is
identical to the timemax passed to 'saextrct' to produce the
input spectrum.
- timeint = INDEF [double precision real]
-
Comma-separated list of time intervals (MET), which is
identical to the timeint passed to 'saextrct' to produce the
input spectrum.
- (deadcorrtype = "LIVE") [string]
-
Name of live-time correction type to use. Use "LIVE" for the true
live-time correction, which results in a true dead-time corrected
spectrum. Use "ON" for the on-time column, which results in a
standard non-dead-time corrected spectrum (the basic detector
on-time exposure is calculated).
- (lcbinarray = INDEF) [string]
-
Setting passed to 'saextrct'. The default value is 400000. Set
to a larger value for better performance on larger data sets.
- (gtiarray = INDEF) [string]
-
Setting passed to 'saextrct'. The default value is 30000. Set to
a larger value if you exceed the limit.
- (maxmiss = INDEF) [string]
-
Setting passed to 'saextrct'. Not recommended to change unless needed.
- (mfracexp = INDEF) [string]
-
Setting passed to 'saextrct'. Not recommended to change unless needed.
- (cleanup = yes) [boolean]
-
Clean up scratch files?
- (chatter = 2) [int] range 0-5
-
Verbosity level of output
- (clobber = no) [boolean]
-
Overwrite output file?
- (history = yes) [boolean]
-
Write standard HEADAS parameter history into output file?
EXAMPLES
- saextrct infile=FS4a_DT outroot=spect gtiandfile=myfile.gti
-
Sample 'saextrct' call which creates spect.pha spectrum with
a user-generated GTI file (myfile.gti). Note: not all parameters
are shown!
- pcadeadspect2 infile=FS4a_DT phafile=spect.pha outfile=spect_DT.pha gtiandfile=myfile.gti
-
Result of step 1 is dead-time corrected using the identical
time filtering criteria as used in step 1. Note: not all parameters
are shown!
CAVEATS
You must use identical time filtering selections for both
your run to saextrct and pcadeadspect2.
BUGS
Please report problems to xtehelp@athena.gsfc.nasa.gov.
SEE ALSO
saextrct, pcadeadcalc2
CATEGORY
Jan95 ftools.xte