Suzaku RIKEN Home
Question ID = 0028

Assigned Member (担当者)

2006-11-10, P. Ranalli
2006-11-14 A. Bamba

Name and Institute (質問者/機関)

Paolo Esposito (INAF)

Question (質問)

Favorite Target; if possible:
Obs ID: 401092010
Obs date: 2006-09-09T22:13:43
Target Name SGR 1806-20

Question [Free format]:

I provide to the task " xissimarfgen" the input information requested as it follows:

xissimarfgen
xissimarfgen version 2006-08-02
Written by Y.ISHISAKI (TMU)
Instrument Name (XIS0,XIS1,XIS2,XIS3)[XIS0] 
source mode (SKYFITS,DETFITS,J2000,SKYXY,DETXY,UNIFORM)[J2000] 
source position R.A. (deg)[272.163833]
source position DEC. (deg)[-20.410972]
number of ARF regions[1] 
region mode (SKYFITS,DETFITS,SKYREG,DETREG,etc)[SKYREG] 
region file #1[xis0_src_arf.reg]xis0_src_arf.reg
limit mode (ACCURACY,NUM_PHOTON)[MIXED]
number of photons for each energy bin[400000]
calculation accuracy for each energy bin[0.005]
input pha file[xis0_src.pha]  xis0_src.pha
XIS det-coordinates mask image file[none] 
input GTI file[ae401092010xi0_0_3x3n000a_cl.evt]
input attitude file[/auxil/ae401092010.att]
input rmf file[xis0.rmf]   
energy step file[estep-default.list] 

The region "xis0_src_arf.reg" is in the coordinate system "Physical"(is this correct?) and the rmf "xis0.rmf" was created with the XIS RMF generator ftool "xisrmfgen" (beta-release). I access the HEASARC CALDB remotely.
The result follows.

     **************************
     ***** Analysis chain *****
     **************************

      ANL HEADAS converter 1.71 for ANL version 1.72

     [ 1] SimASTE_Root            version 2.1
     [ 2] SimASTE_XISarfPhotonGen version 2.1
     [ 3] SimASTE_XRTsim          version 2.0
     [ 4] SimASTE_XRTOUTtoDET     version 1.3
     [ 5] SimASTE_XISarfBuild     version 2.1

HBOOK_BUF_SIZE = 2200000
BNK_BUF_SIZE   = 120000
VERBOSE_LEVEL  = -1
PROFILE_MODULE = YES
CHATTER        = 2

SimASTE_Root:  *** show parameter ***

    TELESCOP            'SUZAKU'
    INSTRUME            'XIS0'
    RAND_SEED           7
    RAND_SKIP           0
    SIMULATION_MODE     1:WEIGHTING
    TELDEF              
'ftp://legacy.gsfc.nasa.gov/caldb/data/suzaku/xis/bcf/ae_xi0_teldef_20060125.fits' 
(CALDB)
    LEAPFILE            
'ftp://legacy.gsfc.nasa.gov/caldb/data/gen/bcf/leapsec_010905.fits' 
(CALDB)

atMissionTime: reading leapsec file 
'ftp://legacy.gsfc.nasa.gov/caldb/data/gen/bcf/leapsec_010905.fits' ...
01: 1972-01-01 00:00:00.000   +0.0 sec
02: 1972-07-01 00:00:00.000   +1.0 sec
03: 1973-01-01 00:00:00.000   +1.0 sec
04: 1974-01-01 00:00:00.000   +1.0 sec
05: 1975-01-01 00:00:00.000   +1.0 sec
06: 1976-01-01 00:00:00.000   +1.0 sec
07: 1977-01-01 00:00:00.000   +1.0 sec
08: 1978-01-01 00:00:00.000   +1.0 sec
09: 1979-01-01 00:00:00.000   +1.0 sec
10: 1980-01-01 00:00:00.000   +1.0 sec
11: 1981-07-01 00:00:00.000   +1.0 sec
12: 1982-07-01 00:00:00.000   +1.0 sec
13: 1983-07-01 00:00:00.000   +1.0 sec
14: 1985-07-01 00:00:00.000   +1.0 sec
15: 1988-01-01 00:00:00.000   +1.0 sec
16: 1990-01-01 00:00:00.000   +1.0 sec
17: 1991-01-01 00:00:00.000   +1.0 sec
18: 1992-07-01 00:00:00.000   +1.0 sec
19: 1993-07-01 00:00:00.000   +1.0 sec
20: 1994-07-01 00:00:00.000   +1.0 sec
21: 1996-01-01 00:00:00.000   +1.0 sec
22: 1997-07-01 00:00:00.000   +1.0 sec
23: 1999-01-01 00:00:00.000   +1.0 sec
24: 2006-01-01 00:00:00.000   +1.0 sec
atMissionTime: 24 lines were read

SimASTE_Root: reading teldef file 
'ftp://legacy.gsfc.nasa.gov/caldb/data/suzaku/xis/bcf/ae_xi0_teldef_20060125.fits'
xissimarfgen: reading gtifile 'ae401092010xi0_0_3x3n000a_cl.evt'
   tstart=211158625.0, tstop=211262385.3, ontime=49104.9, ngti=28, 
nt=10000
xissimarfgen: reading attitude '/auxil/ae401092010.att'
   EA=( 272.1652 110.4732 180.0557 ) at t=211158625.0 
[2006-09-09T23:10:24]
   EA=( 272.1626 110.4680 180.0504 ) at t=211262385.3 
[2006-09-11T03:59:44]
xissimarfgen: WARNING: reading RA_NOM, DEC_NOM failed, calculate fromEULER
    RA_NOM=272.1646, DEC_NOM=-20.4718
xis_contami: reading 
'ftp://legacy.gsfc.nasa.gov/caldb/data/suzaku/xis/bcf/ae_xi0_contami_20061016.fits' 
..

   ngp=522, t0=177242000.0, t1=492342800.0, nen=7900, e0=0.201, e1=15.999
xissimarfgen: reading estepfile 'estep-default.list'
   num_ebin=7, emin_min=0.200, emax_max=16.000
xissimarfgen: making region maps
xissimarfgen: fits_find_rows() failed
Dumping CFITSIO error stack:
--------------------------------------------------
ffgcnn could not find column: QPARAM
ffgtcs could not find all the celestial coordinate keywords
Error reading Region file
--------------------------------------------------
CFITSIO error stack dump complete.
Task xissimarfgen 2006-08-02 terminating with status -1

Am I doing something wrong?

Annotation: using a region in the coordinate system "Image", the simulation starts but, after a very long time, the task creates an arf that doesn't work.

Answer (答え)

Answer from: Paolo Esposito

I found the problem with xissimarfgen: when the task reads the region xis0_src_arf.reg it finds 'physical' instead of 'circle(...)'. During the analysis of the Suzaku data I changed the ds9 version from 3.0 to 4.0, and the latter writes the region files in a slightly different way (that's why the simulation started using the region in the 'image' region: it was generated using the DS9 3.0). 2006-11-10 Accept 
2006-11-10 Forward to SOFT team
2006-11-10 Answer from Ishisaki-san sent to the user
2006-11-14 Done


Suzaku Official Page Suzaku Help Desk Japan Suzaku Help Desk GOF Suzaku RIKEN Page RIKEN Makishima Cosmic Radiation Laboratory