Suzaku WAM 7 Step Manual ------------------------- ver 0.1 2008 June 16, Kazu Yamaoka (Aoyama Gakuin Univ.) ver 0.2 2009 Jan. 22, Kazu Yamaoka (Aoyama Gakuin Univ.) ver 0.3 2009 Feb. 4, Kazu Yamaoka (Aoyama Gakuin Univ.) 1. Obtain the data and information Please see the HXD 7 step guide. http://www.astro.isas.jaxa.jp/suzaku/analysis/7step_HXD_20080501.txt 2. Reprocessing with the latest CALDBs. We do not need to reprocess the WAM TRN data in current versions, but need to reprocess the BST data as shown in the Appendix A. 3. Make light curves for TRN and BST data 3.1) TRN data NOTE: a)energy channels (min_channel and max_channel) should be specified. The approximate energy channels are calculated by Energy (keV) = channel * 511/15.5 b)deadtime should be corrected. --> dt_cor=yes e.x.) hxdmkwamlc input_name=ae101005040hxd_0_wam_uf.evt.gz outroot=ae101005040hxd_0 (as you like) leapfile=CALDB min_channel=2 max_channel=54 dt_cor=yes 3.2) BST data NOTE: a)Only TH data can be processed. PH data are not supported yet. b)Deadtime correction is not supported for TH data. e.x.) hxdmkbstlc input_name=ae101005040hxd_0_bst01_uf.evt.gz outroot=ae101005040hxd_0_bst01(as you like) leapfile=CALDB tpu_board=-1 th_mode=1 3.3) After making the light curves. The FITS file of the light curves will be created, then you can enjoy the timing analysis using the XRONOS package. The TIME column means ASTETIME, i.e. the time (seconds) since 2000-01-01 00:00:00 (UTC). You can convert the ASTETIME into UTC as follows e.x.) sec2time datezero=2000-01-01 timezero=00:00:00 offset=1.8e8 leapfile=FTOOLS or by the date converter at the HEASARC (http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/xTime/xTime.pl). 4. Make energy spectra for TRN and BST data. 4.1) TRN data NOTE: a)Deadtime should be corrected. (dt_cor=yes) b)The integration time interval should be given in ASTETIME. (in time_min and time_max) e.x.) hxdmkwamspec input_name="ae20040824_1853_1959_hxd_wam.uff" \ outroot="ae20040824_1853_1959" tpu_board=-1 time_min=146689200\ time_max=146691000 dt_cor=yes 4.2) BST data NOTE: a)Only PH data can be processed. TH 4 energy band spectra are not supported. (th_mode=0) b)Deadtime should be corrected.(dt_cor=yes) c)The time interval for spectral integration should be given in ASTETIME. (in time_min and time_max) e.x.) hxdmkbstspec input_name="ae20040824_1853_1959_hxd_wam.uff" \ outroot="ae20040824_1853_1959" tpu_board=-1 th_mode=0 \ time_min=146689200 time_max=146691000 dt_cor=yes 4.3) After making the spectral files. You should also make the background spectrum in the same way. 5. Response matrices. The WAM response generator is not available yet. The energy response for some GRBs will be available at the web site http://www.astro.isas.jaxa.jp/suzaku/HXD-WAM/WAM-GRB/ , but please request the WAM team for your interested sources. 6. Spectral fitting. The spectra fitting will be done by XSPEC using the source plus background spectrum, background spectrum and response file. The current response has large uncertainties in low energy channels below ~120 keV, so ignore **-3 High energy channels which we should choose depend on statistics. Please enjoy it. 7. Other useful WAM sites. 6.1) The WAM GRB(Gamma-Ray Burst) table Trigger events: http://www.astro.isas.jaxa.jp/suzaku/HXD-WAM/WAM-GRB/grb/trig/grb_table.html Un-trigger events: http://www.astro.isas.jaxa.jp/suzaku/HXD-WAM/WAM-GRB/grb/untrig/grb_table.html 6.2) The WAM SGR (Soft Gamma Repeater) table Trigger events: http://www.astro.isas.jaxa.jp/suzaku/HXD-WAM/WAM-GRB/sgr/trig/sgr_table.html Un-trigger events: http://www.astro.isas.jaxa.jp/suzaku/HXD-WAM/WAM-GRB/sgr/untrig/sgr_table.html 6.3) The WAM solar flare table Trigger events: http://www.astro.isas.jaxa.jp/suzaku/HXD-WAM/WAM-GRB/solar/trig/solar_table.html Un-trigger events: http://www.astro.isas.jaxa.jp/suzaku/HXD-WAM/WAM-GRB/solar/untrig/solar_table.html 6.4) The WAM event candidate list http://www.astro.isas.jaxa.jp/suzaku/HXD-WAM/WAM-GRB/list/wam_eventlist.html Appendix A. How to reprocess the WAM data A.1) TRN data hxdwamtime read_iomode="overwrite" input_name="ae20041215_1055_1341_hxd_wam.uff" \ leapfile="leapsec.fits" hklist_name="@hk.list" \ tim_filename="ae20041215_1055_1341.tim" A.2) BST data hxdbsttime read_iomode="overwrite" input_name="ae20040824_1853_1959_hxd_bst_0.uff" \ hklist_name="@hk.list" leapfile="refdata/leapsec.fits" \ tim_filename="ae20041215_1055_1341.tim" bstidt_fname=CALDB