30.01.2013 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!