Reading Data

Reading Data Functions

1. sitar_readasm( file [, tstart, tstop, maxchi2, chnl, mjd] )
   Read events from an RXTE-All Sky Monitor FITs file.
   Run as: isis> event = sitar_readasm(file,tstart,tstop,maxchi2,chnl,mjd);
   Variables omitted or set to 0 take on default values.
   Variables in [] are optional, but are order specific.
   
   Inputs:
   • file : Root of asm filename, e.g. 'xa_cygx1_d1', leaving off the '.lc' or '.col'
   Optional Inputs:
   • tstart : Start of time-interval to be read (units of MET or MJD) (Default = file beginning)
   • tstop : Stop of time-interval to be read (units of MET or MJD) (Default = file ending)
   • maxchi2 : The maximum reduced chi2 for an "ASM solution" which will be accepted to consider a data point
   • chnl : 0 for total band (source.lc files = default)
     !=0 for three ASM colors + total (source.col files)
     Energy bands are as follows:
     ch1: ASM channel 410...1188 = 1.3 - 3.0 keV
     ch2: ASM channel 1189...1860 = 3.0 - 5.0 keV
     ch3: ASM channel 1861...4750 = 5.0 - 12.2 keV
     Note: One does *not* have to pre-process the ASM color file into separate channel files using the ftools 'asmchannel'. This will work directly on the 'source.col' files.
   • mjd : !=0 changes from default of Mission Elapsed Time (days) to MJD = JD - 2,400,000.5
   Outputs:
   • event.time : Time of each event
   • event.rate : Total ASM count-rate
   • event.err : Uncertainty in count rate
   • event.chi2 : Chi2 of ASM solution
   Optional Outputs:
   • event.[ch1,ch2,ch3]_rate : ASM count rate in each channel
   • event.[ch1,ch2,ch3]_err : ASM count rate error in each channel
   • event.[ch1,ch2,ch3]_chi2 : Chi2 for ASM solution in each channel (Overrides event.chi2, which won't be output)