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

heasarc.nasa.gov
from heasarc.nasa.gov More from this publisher
30.01.2013 Views

5.5. MASK WEIGHTING (RAY TRACING) 37 Checking for old software versions The SDC has used several versions of the software throughout the mission life-time, and they are not always the best versions. Versions of the software before HEASOFT 6.0.3 contained an error which affects the flux determination of far off-axis bursts. HEASOFT 6.1.2 introduced some new filtering techniques which should slightly improve sensitivity. You can check which version of the software was used with the following command: fkeyprint sw00145675000bevshsp_uf.evt.gz BATCREAT You should see the following result: # EXTENSION: 1 BATCREAT= ’batmaskwtevt 1.11’ / BAT Program that modified this FITS file If you see a version earlier than 1.11, then you should definitely re-run the mask-weighting task. If you see a version earlier than 1.16, then you might consider re-running the mask weighting task, which could gain a slight improvement in sensitivity. What is an auxiliary raytracing file? An auxiliary raytracing file contains information about the position of the GRB and the mask geometry as a function of time within the BAT field of view. This information is needed to make an accurate spectral response before, during and after the spacecraft is slewing to the burst. As of late 2005, the SDC has started to make an auxiliary raytracing file for each observation. Before that time, the file will not be present, and you will need to re-run the mask-weighting tool in order to produce one. 5.5.3 Prerequisites We will practice the mask weighting on GRB 050713A (obsid 00145675 000). Here is what is needed: • Known position of the GRB (RA and Dec) • Attitude file (obsid/aux/swNNNNNNNNNNNsat.fits.gz) • BAT event data (in obsid/bat/event/swNNNNNNNNNNNbevshsp uf.evt.gz); • BAT quality map (in obsid/bat/hk/swNNNNNNNNNNNbcbdq.hk.gz) In this case, let us enter the revised position as reported by Castro-Tirado (GCN 3584): RA=320.5397083, Dec=+77.074861, and apply the mask weighting operation to the original event list. 5.5.4 Initial Steps Change directories to the location of the event data.

38 CHAPTER 5. BAT ANALYSIS PROCEDURES cd 00145675000/bat/event If the event file is compressed, unzip it. The mask-weighting task modifies the event file in place, and gzipped files cannot be modified. gunzip sw*.evt.gz 5.5.5 Compute Mask Weights for BAT Events The mask-weighting task for events is called ‘batmaskwtevt’. Here is an example run, batmaskwtevt infile=sw00145675000bevshsp_uf.evt \ attitude=../../auxil/sw00145675000sat.fits.gz \ ra=320.5397083 dec=+77.074861 \ detmask=../hk/sw00145675000bcbdq.hk.gz \ rebalance=YES corrections=default \ auxfile=/local/data/gcn5b/craigm/sw00145675000bevtr.fits clobber=YES • infile gives the event file name • attitude gives the attitude file name • ra and dec give the revised sky position of the GRB (see above) • detmask is a quality filter which excludes disabled and noisy detectors; this is important so than an accurate accounting of the number of illuminated detectors can be computed • rebalance=YES means that the resulting light curves will be background-subtracted • corrections=default applies the default corrections to the mask weights • auxfile gives the output auxiliary raytracing file, to be used for spectral analysis This analysis may take some time (several minutes), as each event needs to be ray-traced. The speed of the task depends on the number of slews that occur during the observation (i.e. more slews require more intensive calculations). 5.6 Light Curve from Event Data This section describes how to generate a light curve from BAT event data. Later sub-sections describe how to make a light curve with custom time binning. Figure 5.2 shows schematically how light curves are created using the batbinevt task. The same task is used to create spectra and detector plane images and detector plane histograms. 5.6.1 Prerequisites • BAT event data (in obsid/bat/event/swNNNNNNNNNNNbevshsp uf.evt.gz); • BAT quality map (in obsid/bat/hk/swNNNNNNNNNNNbcbdq.hk.gz).

5.5. MASK WEIGHTING (RAY TRACING) 37<br />

Checking for old software versions<br />

<strong>The</strong> SDC has used several versions of the software throughout the mission life-time, and they are<br />

not always the best versions. <strong>Version</strong>s of the software before HEASOFT 6.0.3 contained an error<br />

which affects the flux determination of far off-axis bursts. HEASOFT 6.1.2 introduced some new<br />

filtering techniques which should slightly improve sensitivity. You can check which version of the<br />

software was used with the following command:<br />

fkeyprint sw00145675000bevshsp_uf.evt.gz <strong>BAT</strong>CREAT<br />

You should see the following result:<br />

# EXTENSION: 1<br />

<strong>BAT</strong>CREAT= ’batmaskwtevt 1.11’ / <strong>BAT</strong> Program that modified this FITS file<br />

If you see a version earlier than 1.11, then you should definitely re-run the mask-weighting task.<br />

If you see a version earlier than 1.16, then you might consider re-running the mask weighting task,<br />

which could gain a slight improvement in sensitivity.<br />

What is an auxiliary raytracing file?<br />

An auxiliary raytracing file contains information about the position of the GRB and the mask<br />

geometry as a function of time within the <strong>BAT</strong> field of view. This information is needed to make<br />

an accurate spectral response before, during and after the spacecraft is slewing to the burst. As<br />

of late 2005, the SDC has started to make an auxiliary raytracing file for each observation. Before<br />

that time, the file will not be present, and you will need to re-run the mask-weighting tool in order<br />

to produce one.<br />

5.5.3 Prerequisites<br />

We will practice the mask weighting on GRB 050713A (obsid 00145675 000).<br />

Here is what is needed:<br />

• Known position of the GRB (RA and Dec)<br />

• Attitude file (obsid/aux/swNNNNNNNNNNNsat.fits.gz)<br />

• <strong>BAT</strong> event data (in obsid/bat/event/swNNNNNNNNNNNbevshsp uf.evt.gz);<br />

• <strong>BAT</strong> quality map (in obsid/bat/hk/swNNNNNNNNNNNbcbdq.hk.gz)<br />

In this case, let us enter the revised position as reported by Castro-Tirado (GCN 3584):<br />

RA=320.5397083, Dec=+77.074861, and apply the mask weighting operation to the original event<br />

list.<br />

5.5.4 Initial Steps<br />

Change directories to the location of the event data.

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

Saved successfully!

Ooh no, something went wrong!