Extract product files from the entered data, using the entered filters and selections. Extract is the new name for the old bin command. Bin still works for the nonce.
The extract command works on the latest set of product event files that you have produced. So if you have run GISCLEAN, SISCLEAN, or SELECT EVENTS, extract will operate on the cleaned or selected files, not the original files. Similarly, if you have done a EXTRACT EVENTS, a later extract will use the output events list.
The possible products are a light curve (EXTRACT CURVE), in either QDP or FITS format (or both), a spectrum (EXTRACT SPEC), an image (EXTRACT IMAGE), or a filtered FITS events list (EXTRACT EVENTS). Any combination of these is also possible simultaneously. Also, EXTRACT ALL is equivalent to EXTRACT ``CURVE IMAGE SPEC''
Extract assumes that there is a GTI extension (whose extname is set in the mission database file) in the same fits file which stores the event data. It accepts timing filters in three forms, ASCII or FITS GTI form, and Xronos window format. It also accepts SAOIMAGE format region descriptor files, which are concatenated into one file called xsel_region.xsl, and fed to the extractor program.
If sisclean has been run, then the cleaned events list (stored in xsel_in_events.xsl) will be used. To clear this use clear sisclean.
If select hk or select mkf have been run, the resulting GTI's will be applied automatically. To clear this use clear hksel, or clear mkfsel.
If a timing selection has been made using filter time cursor, this too will be applied automatically.
The various bin sizes, etc, are inherited from the set command. Some of them can also be changed on the fly, by entering hidden parameters, but these values only apply to this invocation of extract.
When extracting events only the events and GTI extensions are written by default. If the copyall=yes option is set then the extra extensions from the (first) input events file are appended to the output file.
Binned spectra, images and lightcurves can be plotted with the PLOT command.
Syntax:
xsel > extract <bin_what> [use_qdp] [binsize_t] [pharebin_t]
[xybinsize_t] [phalcut_t] [phahcut_t]
[xcenter_t] [ycenter_t] [xysize_t] [exposure]
[use_qdp] [wmapcalfile]
Examples :
a)
xsel > extract "CURVE EVENT"
Produces a light curve, and also a filtered event list. To use this filtered event list for the next run of extract, say
xsel > save event my_evtlist.fits yes
This is equivalent to saving the event list, and re-entering it using the read command.
b)
xsel > extract ALL
this is equivalent to
xsel > extract "CURVE IMAGE SPECTRUM"
c)
xsel > extract "CURVE SPECT" bins = 30 pharebin = 4
Extracts a light curve and a spectrum, and temporarily resets the binsize and pha- rebinning factors.
d)
xsel > extract "CURVE SPECT" exposure=0.0
The exposure parameter controls the rejection of bins from the light curve (only), based on how much of the bin lies within the GTI's. If exposure=1, then a bin is only written to the light curve if it lies wholly within a GTI. If exposure=0, then bins are written to the light curve if any portion of them lies within a GTI. Exposure can have any value between these two extremes. Of course, if a bin does not fall inside any of the GTI's, it will not be written to the light curve, no matter what the value of exposure is.
xsel > extract "CURVE SPECT" bins = 30 pharebin = 4