Seven Steps to the Suzaku HXD analysis
Authors: | The HXD team |
---|---|
Date: | Apr., 2012 |
Version: | Processing version 2.2.x.x, HEADAS version 6.10 or later |
Contanct: | suzaku_help(at)crab.riken.jp |
Table of contents
1 Purpose
The purpose of the document is to guide HXD users to appropriate sections in the ABC guide of Suzaku or other web pages.
This document is valid for the version 2.x products. If you dare want to use v1.x products, please check Version 1.X processed data.
2 Setup
2.1 HEADAS
Before starting the Suzaku analysis, you need to install the latest version of the HEADAS package from lheasoft.
3 Correct Data and Information
3.1 Make sure the last information
Users should check important events and processing version.
Go to 4.1 Important Events and Processing Version.
Also, users should check any updates.
Go to 4.2 Checking for Updates.
3.2 Obtain the desired data set of your observation
Go to 3.1.1 Retrieving the data
For example, when your OBS_ID is 101005040, you can download whole directory tree with the wget command as follows:
unix% wget -rL -nH --cut-dirs=4 ftp://darts.isas.jaxa.jp/pub/suzaku/ver2.1/101005040
5 HXD-GSO
5.1 Download NXB data
You should check New GSO Gain Annouce Page to get appropriate background files.
5.2 Check if reprocess is needed
You should check New GSO Gain Annouce Page to see if you need reprocess your data.
When you need to reprocess the data, then
6 Frequantly Asked Question
6.1 HXD Perl scprits do not work on Mac, why?
HXD Perl scripts refer to 7.5.6 hxdpinxbpi and 7.8.1 hxdpinxblc or 7.7.2 hxdgsoxbpi. and 7.9.1 hxdgsoxblc and also you need 5.4.1 Reprocessing and Screening : Aepipeline for reprocessing your data. All of them are written in Perl scripts, so you should set the following enironmental variable
export LHEAPERL=/usr/bin/perl
for bash users:
setenv LHEAPERL /usr/bin/perl
for tcsh users.
If you do so, then these scripts should work.
6.2 Cross Normalization
Go to 5.7.2 Flux Cross Calibration: XIS and HXD.
Generally, crossnormarization among the three detectors should be as follows.
Nominal | XIS | HXD-PIN | HXD-GSO |
---|---|---|---|
XISnom | 1.0 | 1.16 | 1.16 |
HXDnom | 1.0 | 1.18 | 1.18 |
As for crossnormalization amond the XISs,
6.3 Systematic Uncertainties on NXB
If you would like to examine the systematic error of the background model for your object, the best way is to extract the NXB spectrum by use of the earth-occulted period and compare it with the model. Please note that there is not always earth-occultation.
You can estimate the duration with the EHK file:
unix% fplot auxil/ae101005040.ehk.gz TIME ELV - /XW pltcmd=" "
- (In this example, there is no earth-occultation because this is a very
- short observation, actually, the step-point for the large maneuver.)
When the elevation of the source (ELV) is less than -5 degree, the FOV of the detector is occulted by the earth and hence the spectrum during this period can be treated as same as the NXB component. Thus, what you should do is to change ELV>5 into ELV<-5 when extracting data. The procedure can be easily done as follows.
- Go to 7.3 HXD Reprocessing and Screening: General Remarks and download pin_mkf.sel ( Note that gso_mkf.sel is the same as pin_mkf.sel).
- Chenge ELV>5 into ELV<-5 in pin_mkf.sel.
- Run aepipeline with specifying hxd_pin_expr= @pin_mkf.sel hxd_gso_expr= @pin_mkf.sel hxd_pse_expr= @pin_mkf.sel.
- To compare the earth with the NXB data, go to 7.5.6 hxdpinxbpi and 7.8.1 hxdpinxblc or 7.7.2 hxdgsoxbpi. and 7.9.1 hxdgsoxblc.
7 History
- August, 2006 by Terada, Kokubun (v0.7)
- Dec, 2006 by Fukazawa, hTakahashi, Terada, Kokubun (v1.2)
- Mar, 2007 by Fukazawa
- Oct, 2007 by Terada (v2.0)
- May, 2008 by Terada, hTakahashi (v2.x) / rev.4
- Oct, 2008 by Terada (v2.x)/ epoch 5
- Jun, 2011 by Nishino / epoch 6-11
- Feb, 2012 by Yamada / new gsogain and pin noise updated
- Apr, 2012 by Yamada / link the background page, and addd a note on LHEAPERL