Analysis of Early Release Data
2005 December 2
The Suzaku project at ISAS/JAXA have released selected calibration observations from the early phase of the mission. Here, we provide a quick guide on how to analyze these data.
These notes should be considered preliminary and experimental. They cover cleaning XIS hot and flickering pixels, making XIS and HXD products (especially spectra), how to use xselect to catalog all the event files for a source, how to filter by grade, how to select events based on housekeeping information, how to make summed images of all 4 XIS.
Some options in xselect can be set before running using the pset command. To change the default mission to Suzaku and to ensure thatall extensions are copied when extracting events files give the following commands from the unix command line
pset xselect default_mission="SUZAKU" pset xselect copyall=yes
The latest version of HEAsoft (6.0.4) should be installed.
cleansis.f and replacing the copy in ftools/asca/src/cleansis in your
To run cleansis on Suzaku XIS event files type from the command line cleansis chipcol=SEGMENT , give the input and output filenames and use the default values of the remaining parameters.
cleansis can be run from within xselect by giving the command sisclean and choosing option number 2.
(If you have a version of HEAsoft earlier than 6.0.4 you will need to download a new xselect.mdb file. To use this you can do one of the following
- replace the copy in your HEAsoft source distribution in ftools/xselect/common and install it using hmake install in this directory
- if you only have the HEAsoft binary distribution then install the new xselect.mdb in ftools/operating-system/bin where operating-system is the appropriate directory for you system (it is the end of the directory string written by eg printenv HEADAS).
- if you don't have write-access to your HEADAS distribution you can setenv XSELECT_MDB directory-tree/xselect.mdb where directory-tree is the location you have placed the file.
There was also a bug in earlier versions that caused events to be lost when reading both 3x3 and 5x5 event files if the 5x5 file was read first.
xsel:ASCA> set mission suzaku xsel:SUZAKU> set instr xis0 xsel:SUZAKU-XIS0 > set datadir xis Setting data directory to /Volumes/Documents/DATA/Suzaku/Fornax/xis/evt/ Setting mkf directory to /Volumes/Documents/DATA/Suzaku/Fornax/auxil/ xsel:SUZAKU-XIS0 > make cat > Default filter for the catalogue (NONE for none, DEF for default) [DEF] Using default selection expression: INSTRUME=='XIS0' ******************** Observation Catalogue ******************** Data Directory is: /Volumes/Documents/DATA/Suzaku/Fornax/xis/ HK Directory is: /Volumes/Documents/DATA/Suzaku/Fornax/xis/ OBJECT DATE-OBS TIME-OBS DATE-END TIME-END EDITMODE 1 Fornax_Clus 2005-09-13 01:29:50 2005-09-14 11:46:46 3x3 2 Fornax_Clus 2005-09-13 15:53:34 2005-09-13 22:59:10 5x5 xsel:SUZAKU-XIS0 > choose 1-2 Setting datamode to STANDARD Getting Min and Max for Energy Column... Got min and max for PI: 0 4095 Number of files read in: 2 Files currently in use: 1 ae20050913_0130_1200_xis0_3x3_cl.evt 2 ae20050913_0130_1200_xis0_5x5_cl.evt(In pre-6.0.4 versions doing this for the HXD actually caused one of the event files to be deleted !)
xsel> filter grade "0,2-6" xsel> filter column "GRADE=0:0 2:6"The latter was required for earlier Suzaku processing versions which did not have the TLMIN/MAX keywords for the GRADE column. new xselect.mdb file, or issue the command
xsel:SUZAKU-XIS0 > set mkfdir ../../auxilin xselect.
After the set-up above we can get a list of available housekeeping parameters by show param mkf. To plot lightcurves of any of these parameters use mkfbin followed by plot mkf.
To select events based on housekeeping information use select mkf and then give either a C-style or Fortran-style selection string. This will create a GTI file which will be used for any subsequent products (using extract or bin). This GTI file can be saved by save mkf_sel.
read/rebin=4 xis0.evt save read/rebin=4 xis1.evt sum save read/rebin=4 xis2.evt sum save read/rebin=4 xis3.evt sum save smoothThe rebinning and final smooth produce an image well-matched to the psf. To produce an image from a restricted range of PI bins replace the read commands by eg
If you have any questions concerning Suzaku, visit the Feedback form.