Data Analysis Tutorial for the Fermi Gamma-ray Burst Monitor (GBM) Hoi-Fung Yu (MPE) for the Fermi GBM team 4 th FAN Workshop, University of Hong Kong, Hong Kong (12 July 2013)
Contents • Basic facts of the Fermi GBM • GBM data types • Analysis tool: RMFIT • Retrieving GBM data • Spectral analysis using RMFIT
What is Fermi GBM? • Fermi GBM is a joint project of MPE and UAH • Joint analysis of spectra and time histories of Gamma-ray Bursts (GRBs) observed by both the GBM and LAT (Meegan et al. 2009) • Provide near-real time burst locations onboard to permit repointing of the spacecraft to obtain LAT observations of delayed emission from bursts, and to disseminate burst locations rapidly to the community of ground-based observers • 12 thallium activated sodium iodide, NaI(Tl), scintillation detectors and 2 bismuth germanate, BGO, scintillation detectors • NaIs cover lower energy spectrum (~8 keV - ~1 MeV), BGOs cover higher energy spectrum (~200 keV - ~40 MeV) • The axes of the NaI detectors are oriented such that the positions of GRBs can be derived from the measured relative counting rates, a technique previously employed by CGRO BATSE
These detectors and the data read-out systems were actually built at MPE! NaI detector (Meegan et al. 2009) BGO detector (Meegan et al. 2009) Locations and orientations of the GBM detectors (Meegan et al. 2009)
GBM data types • TTE (Time-tagged events): Event data. Individual photons (or count rates) with time tag and energy. 128 Energy channels (adjustable boundaries). Each row in event extension in the FITS file is a time and an energy channel. • CSPEC: 4.096 s (nominal) and 1.024 s (during burst) temporal resolution, 128 energy channels (adjustable boundaries, same as TTE). Each row in spectrum extension in the FITS file is a 128 energy channel spectrum. • CTIME: 0.256 s (nominal) and 64 ms (adjustable from 64 ms to 1.024 s during burst) temporal resolution, 8 energy channels (adjustable boundaries). Not suitable for spectroscopy.
Which data type should I use? • Look at the quicklook lightcurve to determine the approximate duration of the burst. • Spectral analysis: • For short (~< few seconds) or highly structured events, choose TTE data. (TTE data has only time range from ~T0-25 s to ~T0+300 s, and is slower to process. • For longer or less structured events, choose CSPEC data. • PHA files from other experiments, e.g. Swift BAT, can also be analyzed using RMFIT, providing they have a single associated response file, and conform to standard FITS format. BAT files are single PHA files, i.e. one spectrum per file. CSPEC/CTIME files are PHA2 files, i.e. multiple spectra per file, each spectrum belonging to a time bin.
• Naming convention of GRB: GRB yymmddX • If there is IPN/Swift detection, X = A, B, C... • For GBM detection only, X is a three digit number • For each GRB, at the FSSC you can find... • 14 TTE files (per each version, e.g. glg_tte_b0_bnyymmddxxx_v01.fit) [FIT] • 14 CSPEC files (per each version, e.g. glg_cspec_b0_bnyymmddxxx_v01.pha) [PHA] • 14 CTIME files (per each version, e.g. glg_ctime_b0_bnyymmddxxx_v01.pha) [PHA] • 1 Trigdat file (e.g. glg_trigdat_all_bn130601547_v02.fit) for burst localization • Various GIF and PDF files for quicklook to the lightcurves • and TCAT files for internal use of burst information
Analysis tool: RMFIT version 4.3.2 (should have installed before this tutorial!) • Operating system: Mac OS X (10.6 or newer), Linux • Supporting softwares: IDL 8.1 (virtual machine) [Good news: this is a self-contained distribution, no external library is needed!] • Download the RMFIT package (version 4.3.2) and the installation instructions at the FSSC user contribution page: http://fermi.gsfc.nasa.gov/ssc/data/analysis/rmfit/ • Choose The Mac OS X Package if you are using OS X • Choose The Linux Distribution if you are using Linux • Extract the tarball and follow the installation instructions • To run RMFIT do one of the following options: • Double click the rmfit432.command file in /Applications/rmfit_v432/ • At a terminal prompt cd into the rmfit_v432 directory and type ./rmfit432.command • Create an alias: alias rmfit=$directory_to_file/rmfit432.command. Then at the command prompt type rmfit. Put this inside your .bashrc to use from all terminals. • You cannot use idl -vm to call RMFIT. You must use one of the two options above. The IDL virtual machine included in this install can only be accessed by RMFIT through the rmfit432.command file. If RMFIT will not run make sure you have X11 installed.
Retrieving GBM data From the spectroscopy tutorial: http://fermi.gsfc.nasa.gov/ssc/data/analysis/user/vc_rmfit_tutorial.pdf (RMFIT team, April 2010) • Directly by FTP at: http://heasarc.gsfc.nasa.gov/FTP/fermi/data/gbm/triggers if you know the full name of the burst e.g. 080916009. Look for the files you want under the “current” directory for the burst. • Using the FSSC interface: http://fermi.gsfc.nasa.gov/ssc/data/access/ under “GBM Burst Catalog”. You will see a page with explanations about GBM data. Click the “Browse this Table” at the top left-hand corner. Search for your burst by entering “GRB080916*” under the “name” field in the table. Select your burst from the returned triggers using your knowledge of the trigger time. “Preview and Retrieve” will allow you to generate the directory listing the available data files for this burst, or you can just retrieve the whole available data package in a tar file. • CSPEC (CTIME) files are “glg_cspec_nx_bn080916009_v0y.pha” (“glg_ctime_nx_bn080916009_v0y.pha”) where x is the NaI detector number 0-9, a, b or “glg_cspec_bx_bn080916009_v0y.pha” (“glg_ctime_bx_bn080916009_v0y.pha”) with x =0 or 1. Choose b0 if you have low-numbered NaI, b1 for NaI with x>5. If you have a mixture, select both BGOs and choose better one after looking at data. • TTE files are “glg_tte_nx_bn080916009_v0y.fit” or “glg_ctime_bx_bn080916009_v0y.pha”)... • In either case, you will need a detector response file for each detector you are analyzing “glg_cspec_nx_bn080916009_v0y.rsp” or “glg_cspec_bx_bn080916009_v0y.rsp”. CSPEC response files apply to both CSPEC and TTE data files: the number of energy channels is 128 for both. • Quicklook lightcurves are called glg_lc_all_bn080916009.gif • The version number v0y is unimportant and can be different for each file type.
Spectral analysis I will use a recent burst GRB 130305A (130305486) as an example (a Fermi and Swift detected burst): http://fermi.gsfc.nasa.gov/ssc/data/access/
GBM Products
1. Choose the brightest 3-4 NaIs, and also the suitable BGO. Usually the triggered NaIs are bright. Use the blocked BGO.
2. Open RMFIT, then load all the desired CSPEC files [PHA].
3. Zoom to around the trigger time. (You can do this just once and then load the lookup files to all other detectors!) You can toggle to switch between the lightcurve and energy channel display windows. You can rebin according to temporal resolution (must be integer multiple of the finest bin) or S/N ratio, or click “Rebin: Full Resolution” to return to full resolution.
3. Toggle to the energy channels display, click “Select Source: Source Interactive” in order to select usable energy ranges. 4. Click on the right-most region to delete the current selection or click on “EXIT” at the bottom-left corner of the display window to confirm selection. For NaIs, usually we use ~8-900 keV. Avoid the overflow channels.
5. In the lightcurve display, click “Fit Background” and then select the background intervals before and after the burst. 6. Click on the right-most region to delete the current selection or click on “EXIT” at the bottom-left corner of the display window to confirm selection. 7. Choose the order of polynomial used to fit the background intervals (RMFIT will then extrapolate the fit into the burst region). 8. Look at the reduced Chi 2 values for each energy channel, they should be all around the value 1. Click “Dismiss” to return to the lightcurve display.
9. Zoom closer to the burst to check whether the fitted background (blue dotted line) is reasonable. If not, repeat the background fitting process by selecting another pair of intervals. 10. Click “Select Source: Source Interactive” to select the source interval. Click “Exit” to confirm or the right-most region to cancel selection. In this tutorial I just choose a ~15 s interval containing the brightest part of the burst. Important: The source interval interval must not intercept with the background intervals!
RMFIT will display the selected energy range at the upper-right corner, and the bottom, and the source time interval at the bottom. You can also toggle to the energy channel display to see the fitted background rates for each energy channel (remember the energy channel display reveals mainly the effective area effect).
11. Click “File” and then “Lookup -> Save Lookup” to save the current binning, zooming, background and source selections. The lookup file will be saved with the same name as your current CSPEC/ TTE file with extension “.lu”. You can also choose your desired name for it.
12. Read the lookup file to every other detectors. For consistency you should always use the same lookup file for all detectors!
Recommend
More recommend
