'sxsanticolc' is a script that takes an SXS antico event file and generates a light curve after optional cleaning. 'sxsanticolc' assumes that the input file has the PI column already calculated using the 'sxsanticopi' task.
The PSP can be selected using the pspid
parameter.
If pspid
is not set to A or B, then all events are used.
Before calculating the light curve the event file is
filtered for all GTI. Any input gtifile
is merged with a GTI attached to the input event file as well as events between
parameters lcstart
and lcstop
.
The event file is then screened for events based on the expr
parameter.
Once the event file has been filtered, the 'extractor' tool is run to create
a light curve for each picut
.
The light curve will contain bins between
lcstart
and lcstop
. If lcstart
or
lcstop
are set to -1 (default) then TSTART or TSTOP are
used for the light curve start and stop times, respectively. Additional
control parameters are lcthresh
and lctzero
.
The columns RATE and ERROR are appended to the output lightcurve
for each PI range with the range in the column name. For example,
if the picut
="0-12200,0-1200,1201-12200" then the
RATE column names in the output file will be:
RATE0_12200, RATE0_1200 and RATE1201_12200.
Optionally, the extract
parameter will create a spectrum using
'extractor' and the full PI range for antico events (PI=0:12200).
1. Extract a lightcurve (sxsanticolc.lc) and spectrum (sxsanticolc.pha) using PSP A, specified GTI file and expression. Use a 20 second time bin with three energy ranges.
sxsanticolc "infile=ah100050050sxs_a0ac_uf.evt" "outroot=sxsanticolc" "bintime=20" "antpsp=A" "expr=FLG_BASELINE==b0&&PI>=0&&PI<=12200&&DURATION>2&&DURATION<19" "gtifile=input/ah100050050sxs_el.gti[4]" "numlc=3" "picut=0-12200,0-1200,1201-12200"
2. Extract a lightcurve using PSP A with a 1 second bin time and constricting the lightcurve between lcstart and lcstop. The default energy range is used (0-12200)
sxsanticolc "infile=input/ah100050050sxs_a0ac_uf.evt" "outroot=output/sxsanticolc" "bintime=1" "antpsp=A" "lcstart=70203712" "lcstop=70205402" "expr=NONE"