Reading Data

Reading Data Functions

1. sitar_readasm( file [; tstart=#, tstop=#, maxchi2=#, chnl=#, mjd, jd] )
   Read events from an RXTE-All Sky Monitor FITs file.
   Run as: isis> event = sitar_readasm(file;tstart=#,tstop=#,maxchi2=#,chnl,mjd);
   Variables in [] are optional qualifiers.
   Omitted qualifiers(i.e., variables after the semi-colon) take on default values.
   
   Inputs:
   • file : ASM filename, e.g. 'xa_cygx1_d1.lc', or 'xa_cygx1_d1.col'. (Be sure to use the proper file if reading colors.)
   Optional Qualifier 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 : If set, changes from default of Mission Elapsed Time (days) to MJD = JD - 2,400,000.5
   • jd : If set, changes from default of Mission Elapsed Time (days) to Julian Data. Supersedes mjd flag.
   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)