6 December 2005
unix% mkdir n2110
Also download the housekeeping files (click on EHK) and save them to your working directory.
unix% tar xvf 100012010.tar
unix% gunzip */*
unix> pset xselect imagedisp="ds9 -unix_only"Also, if you would like to define the default mission to Suzaku:
unix> pset xselect default_mission="SUZAKU"
unix% xselect ** XSELECT V2.3 ** > Enter session name >[xsel] n2110 n2110:ASCA > set mission suzaku n2110:SUZAKU > set instr xis0 n2110:SUZAKU-XIS0 > set datadir xis Setting data directory to /home/lgallo/SUZAKU/ngc2110/xis/ Setting mkf directory to /home/lgallo/SUZAKU/ngc2110/auxil/
n2110:SUZAKU-XIS0 > read events ae20050916_0850_1956_xis0_3x3_cl2.evt Notes: XSELECT set up for SUZAKU Time keyword is TIME in units of s Default timing binsize = 16.000 Setting... Image keywords = X Y with binning = 8 WMAP keywords = X Y with binning = 8 Energy keyword = PI with binning = 1 Getting Min and Max for Energy Column... Got min and max for PI: 0 4095 Got the minimum time resolution of the read data: 8.0000 MJDREF = 5.1544000742870E+04 with TIMESYS = TT Number of files read in: 1 ******************** Observation Catalogue ******************** Data Directory is: /home/lgallo/SUZAKU/ngc2110/xis/ HK Directory is: /home/lgallo/SUZAKU/ngc2110/xis/ OBJECT DATE-OBS TIME-OBS DATE-END TIME-END EDITMODE 1 NGC_2110 2005-09-16 08:49:42 2005-09-18 19:55:50 3x3In this example I am only considering one event file for simplicity. In principle, since mode-specific responses are not yet available, the 5x5, 3x3 and 2x2 data for each telescope can be combined, in which case one would group the event files with quotations, i.e. :
read events "file1.evt file2.evt".Notice that the event file *_cl2.evt is the file which has been cleaned with "cleansis". Also, keep in mind that XSELECT assumes the housekeeping files are in the directory ../auxil, relative to the working directory (usually xis/).
n2110:SUZAKU-XIS0-STANDARD > extract all extractor v4.56 7 Sep 2005 Getting FITS WCS Keywords Doing file: /home/lgallo/SUZAKU/ngc2110/xis/ae20050916_0850_1956_xis0_3x3_cl2.evt 100% completed Total Good Bad: Region Time Phase Grade Cut 235840 234116 0 0 0 0 1724 Grand Total Good Bad: Region Time Phase Grade Cut 235840 234116 0 0 0 0 1724 in 61064. seconds Fits light curve has 234116 counts for 3.834 counts/sec Thresholding removed significant counts Try exposure=0.0 on the extract command in xselect or lcthresh=0.0 if running extractor stand-alone Thresholding removed more than half the bins Try exposure=0.0 on the extract command in xselect or lcthresh=0.0 if running extractor stand-alone Spectrum has 234116 counts for 3.834 counts/sec ... written the PHA data Extension Image has 234116 counts for 3.834 counts/sec
To view the light curve, spectrum, image:
xsel:SUZAKU-XIS-1-STANDARD > plot curve xsel:SUZAKU-XIS-1-STANDARD > plot spectrum n2110:SUZAKU-XIS0-STANDARD > plot imageThe "standard" (still evolving) filter criteria have already been applied to this data; they correspond to:
xsel:SUZAKU-XIS-1-STANDARD > select mkf "SAA==0 && ELV > 5 && DYE_ELV > 20 && T_SAA > 256"You can also set any other time, space, or energy filters as will be demonstrated below.
The goal here is to produce a pha file which can be analysed in XSPEC.
n2110:SUZAKU-XIS0-STANDARD > plot image
Region --> file format --> SAOimage
n2110:SUZAKU-XIS0-STANDARD > filter region src.regTo confirm that you have done this, you can type "show stat" at anytime:
n2110:SUZAKU-XIS0-STANDARD > show stat *** Status of XSELECT *** Plot device is /XW *** MISSION *** SUZAKU NONE XIS0 STANDARD Time keyword is TIME in units of s Default timing binsize = 16.000 Minimum timing resolution: 8.0000 Energy keyword is PI with binning 1 Max and min for PI : 0 4095 Image keywords = X Y with binning = 8 WMAP keywords = X Y with binning = 8 *** DATA *** The data directory is /home/lgallo/SUZAKU/ngc2110/xis/ The mkf directory is /home/lgallo/SUZAKU/ngc2110/auxil/ The obscat n2110_read_cat.xsl has been made. One data file has been read in. GTI files are in use. *** PRODUCTS *** A spectrum has been accumulated. A curve has been accumulated. An image has been accumulated. *** SELECTIONS *** NONE *** FILTERS *** Spatial selections have been read in from: src.reg
n2110:SUZAKU-XIS0-STANDARD > extract spectrum extractor v4.56 7 Sep 2005 Getting FITS WCS Keywords Doing file: /home/lgallo/SUZAKU/ngc2110/xis/ae20050916_0850_1956_xis0_3x3_cl2.evt 100% completed Total Good Bad: Region Time Phase Grade Cut 235840 179064 56637 0 0 0 139 Grand Total Good Bad: Region Time Phase Grade Cut 235840 179064 56637 0 0 0 139 in 61064. seconds Spectrum has 179064 counts for 2.932 counts/sec ... written the PHA data Extension
n2110:SUZAKU-XIS0-STANDARD > save spectrum n2110_xis0_3x3_srcThe .pha extension is added.
n2110:SUZAKU-XIS0-STANDARD > filter region -src.regNow apply the background extraction region
n2110:SUZAKU-XIS0-STANDARD > filter region bkg.regIn principle, multiple files can be handled at one time thus we could remove one region and add another in one step. However, care must be taken in placing the files in the correct order. Alternatively, one can simply type:
n2110:SUZAKU-XIS0-STANDARD > clear regionto remove all existing regions, before adding a new region.
n2110:SUZAKU-XIS0-STANDARD > show stat *** Status of XSELECT *** Plot device is /XW *** MISSION *** SUZAKU NONE XIS0 STANDARD Time keyword is TIME in units of s Default timing binsize = 16.000 Minimum timing resolution: 8.0000 Energy keyword is PI with binning 1 Max and min for PI : 0 4095 Image keywords = X Y with binning = 8 WMAP keywords = X Y with binning = 8 *** DATA *** The data directory is /home/lgallo/SUZAKU/ngc2110/xis/ The mkf directory is /home/lgallo/SUZAKU/ngc2110/auxil/ The obscat n2110_read_cat.xsl has been made. One data file has been read in. GTI files are in use. *** PRODUCTS *** A spectrum has been accumulated. A curve has been accumulated. An image has been accumulated. *** SELECTIONS *** NONE *** FILTERS *** Spatial selections have been read in from: bkg.reg
n2110:SUZAKU-XIS0-STANDARD > extract spectrum extractor v4.56 7 Sep 2005 Getting FITS WCS Keywords Doing file: /home/lgallo/SUZAKU/ngc2110/xis/ae20050916_0850_1956_xis0_3x3_cl2.evt 100% completed Total Good Bad: Region Time Phase Grade Cut 235840 1240 234522 0 0 0 78 Grand Total Good Bad: Region Time Phase Grade Cut 235840 1240 234522 0 0 0 78 in 61064. seconds Spectrum has 1240 counts for 2.0307E-02 counts/sec ... written the PHA data Extension
n2110:SUZAKU-XIS0-STANDARD > save spec n2110_xis0_3x3_bkg
n2110:SUZAKU-XIS0-STANDARD > plot image
Region --> file format --> SAOimage
n2110:SUZAKU-XIS0-STANDARD > filter region src.regTo confirm that you have done this, you can type "show stat" at anytime:
n2110:SUZAKU-XIS0-STANDARD > show stat *** Status of XSELECT *** Plot device is /XW *** MISSION *** SUZAKU NONE XIS0 STANDARD Time keyword is TIME in units of s Default timing binsize = 16.000 Minimum timing resolution: 8.0000 Energy keyword is PI with binning 1 Max and min for PI : 0 4095 Image keywords = X Y with binning = 8 WMAP keywords = X Y with binning = 8 *** DATA *** The data directory is /home/lgallo/SUZAKU/ngc2110/xis/ The mkf directory is /home/lgallo/SUZAKU/ngc2110/auxil/ The obscat n2110_read_cat.xsl has been made. One data file has been read in. GTI files are in use. *** PRODUCTS *** A spectrum has been accumulated. A curve has been accumulated. An image has been accumulated. *** SELECTIONS *** NONE *** FILTERS *** Spatial selections have been read in from: src.regNotice that the default bin size is 16 s.
n2110:SUZAKU-XIS0-STANDARD > extract curve extractor v4.56 7 Sep 2005 Getting FITS WCS Keywords Doing file: /home/lgallo/SUZAKU/ngc2110/xis/ae20050916_0850_1956_xis0_3x3_cl2.evt 100% completed Total Good Bad: Region Time Phase Grade Cut 235840 179064 56637 0 0 0 139 Grand Total Good Bad: Region Time Phase Grade Cut 235840 179064 56637 0 0 0 139 in 61064. seconds Fits light curve has 179064 counts for 2.932 counts/sec Thresholding removed significant counts Try exposure=0.0 on the extract command in xselect or lcthresh=0.0 if running extractor stand-alone Thresholding removed more than half the bins Try exposure=0.0 on the extract command in xselect or lcthresh=0.0 if running extractor stand-alone
n2110:SUZAKU-XIS0-STANDARD > save curve n2110_xis0_3x3_srcb_16sThe extension .lc is added automatically.
n2110:SUZAKU-XIS0-STANDARD > filter region "-src.reg bkg.reg"or,
n2110:SUZAKU-XIS0-STANDARD > clear region n2110:SUZAKU-XIS0-STANDARD > filter region bkg.regAnd check to see that this was done correctly:
n2110:SUZAKU-XIS0-STANDARD > show stat *** Status of XSELECT *** Plot device is /XW *** MISSION *** SUZAKU NONE XIS0 STANDARD Time keyword is TIME in units of s Default timing binsize = 16.000 Minimum timing resolution: 8.0000 Energy keyword is PI with binning 1 Max and min for PI : 0 4095 Image keywords = X Y with binning = 8 WMAP keywords = X Y with binning = 8 *** DATA *** The data directory is /home/lgallo/SUZAKU/ngc2110/xis/ The mkf directory is /home/lgallo/SUZAKU/ngc2110/auxil/ The obscat n2110_read_cat.xsl has been made. One data file has been read in. GTI files are in use. *** PRODUCTS *** A spectrum has been accumulated. A curve has been accumulated. An image has been accumulated. *** SELECTIONS *** NONE *** FILTERS *** Spatial selections have been read in from: bkg.reg
n2110:SUZAKU-XIS0-STANDARD > save curve n2110_xis0_3x3_bkg_16sThe background light curve can now be subtracted from the source (+background) light curve in the standard way using the XRONOS package.
The process is essential the same as extracting an XIS spectrum. If available, cleaned event files for the PIN (*_hxd_wel_pin_cl.evt.gz) and/or GSO (*_hxd_wel_gso_cl.evt.gz) are found in the source directory.
After unzipping the files, each event file can be read into xselect as usual:
n2110:SUZAKU-HXD-STANDARD > read events *_hxd_wel_pin_cl.evt
Without need to define an extraction region, a source plus background spectrum can be extracted and saved:
n2110:SUZAKU-HXD-STANDARD > extract spectrum n2110:SUZAKU-HXD-STANDARD > save spectrum *_hxd_wel_pin_clThe created pha file can now be read by XSPEC. The background and response matrices are available at http://www.astro.isas.jaxa.jp/suzaku/pub/hxd/
This document is work in progress a will be improved over the following days.