The SWIFT BAT Software Guide Version 6.3 30 ... - HEASARC - Nasa
The SWIFT BAT Software Guide Version 6.3 30 ... - HEASARC - Nasa
The SWIFT BAT Software Guide Version 6.3 30 ... - HEASARC - Nasa
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
44 CHAPTER 5. <strong>BAT</strong> ANALYSIS PROCEDURES<br />
5.7.3 Change into Event Directory<br />
Change directories to the location of the event data.<br />
cd 00145675000/bat/event<br />
5.7.4 Make the spectrum<br />
<strong>BAT</strong> spectra are made with the ‘batbinevt’ task. Here is an example for the flare discussed above:<br />
batbinevt infile=sw00145675000bevshsp_uf.evt.gz outfile=flare.pha outtype=PHA \<br />
timedel=0.0 timebinalg=u tstart=142921843 tstop=142921868 \<br />
energybins=CALDB:80 outunits=RATE \<br />
detmask=../hk/sw00145675000bcbdq.hk.gz clobber=YES<br />
• the outtype=PHA means make a “spectrum” (PHA = pulse height analyzer)<br />
• the timedel=0.00 means accumulate all the data in the specified time range into one spectrum<br />
• the timebinalg=u means make them uniform bins (really “don’t care” in this case)<br />
• the tstart and tstop give the start and stop time of the spectrum as discussed above<br />
• the energybins=CALDB:80 means make a default 80-channel spectrum<br />
• the outunits=RATE makes a rate spectrum, which is usually required since <strong>BAT</strong> spectra<br />
have fractional counts, which XSPEC does not accept<br />
• the detmask is required to filter out noisy detectors<br />
5.7.5 Apply Corrections to the Spectrum<br />
<strong>The</strong>re are several corrections that need to be applied to a spectrum before it can be fit by XSPEC.<br />
<strong>The</strong> first correction involves inserting the proper ray-tracing keywords into the file.<br />
Several geometric parameters are required to make an accurate response matrix, such as the<br />
position of the burst in instrument coordinates. Because the spacecraft is often slewing during<br />
the burst, it is not possible to know for sure what those geometric parameters should be until<br />
the desired spectrum has been created. <strong>The</strong> ‘batupdatephakw’ task does this for you, using the<br />
auxiliary ray-tracing file described above. If the auxiliary raytracing file does not exist in your<br />
observation, then you will have to create it using the Mask-Weighting thread. If you do not apply<br />
this correction, then fluxes may be in error by up to 50%, and power law indices may be off by up<br />
to 0.5.<br />
batupdatephakw flare.pha sw00145675000bevtr.fits.gz<br />
Next, you should apply the <strong>BAT</strong>-recommended systematic error to your spectrum. <strong>The</strong> ‘batphasyserr’<br />
task is used to do this, and the error vector is stored in the calibration database.<br />
batphasyserr flare.pha CALDB