hxtlcurv - produce background-corrected RXTE HEXTE light curves
hxtlcurv [-i "hexte science data file name"] [-k "hexte dead-time coefficients file name"] [-r "hexte 16-s hk file name (contains Rates)"] [-b "bin size in seconds"] [-g "good-time-interval file name"] [-t "timeint file name, produced by timetrans"] [-l "lower energy channel number for light curve"] [-u "upper energy channel number for light curve"] [-d "data column list (archive files only)"] [-e "detector list for dead time correction"] [-c "chatter level for hxtdead"] [-p "make pha files too"] [-v (version number written to standard output)] [-h (help text written to standard output)]
hxtlcurv produces a background subtracted light curve and optionally an average spectrum from HEXTE science or archive data using the off-source rocking measurements. Standard interval selection options are supported. In addition to the data file, required inputs include the contemporaneous HEXTE rates file and deadtime coefficients files. The hxtlcurv script queries the user for each of these as needed. If no timeint file is given, the TSTART and TSTOP keyword values from the science file are used. The script then checks for source and background files having the same name as the sci- ence file with the addition of "_src" or "_bkg" suffixes, respectively. If these are not found, hxtback is run to produce them. Then either saextrct or seextrct is run to produce the light curve files, which are passed to hxtdead for deadtime correction. Finally, a background light curve is interpolated and passed to lcmath for background subtrac- tion from the source-looking light curve.
Limited by the rocking and the extract tools to time bins of 16 seconds or less, and 16 seconds only for archive spectral data. Calculated errors assume incorrectly that background estimates are independent, not quite true. The RMS is under- estimated by 20 percent in the worst case of weak source and 16-sec time bins. Some point-to-point correlation up to 32 seconds is also present. Rocking data only. No list-of-files option for input.
The script has not been made operational under VMS. Channel selection with archive multiscalar data produces zero rates. This appears to be a bug in saextrct. Use event-list data if you want a light curve with time bins shorter than 16 sec and a limited channel range.
Please report problems to xtehelp@athena.gsfc.nasa.gov.
hxtback, hxtdead, saextrct, seextrct, finterp and lcmath for more information.