Question ID = 0045
Question ID = 0045
Assigned Member (担当者)2007-01-31 A. BambaName and Institute (質問者/機関)寺田幸功 (理研)Question (質問)HEADAS 6.1.2 に入っている xissimarfgen version 2006-11-26 に関する質問です。 下記のような感じで、DET 座標 (1-bin サイズ) の region file を用いて XIS の Arf を作成すると、添付図 黒のように、ほぼ期待通りの arf が得られます。 $ cat reg/xis0_2005src_det.reg circle(283,500,250) -box(513,771,1024,366,0) -box(514,265,1026,386,0) $ xissimarfgen instrume=XIS${XISID} \ source_mode=DETXY source_x=283 source_y=500 \ num_region=2 region_mode=DETREG \ regfile1=${REG_NAME_SRC} regfile2=${REG_NAME_BGD} \ arffile1=${ARF_NAME_SRC} arffile2=${ARF_NAME_BGD} \ limit_mode=MIXED num_photon=${SIM_MUM_PHOTON} accuracy=0.005 \ phafile=${PHA_NAME} detmask=${CALDM_CALMASK} gtifile=none \ date_obs=${DATE_OBS} \ attitude=${ATTFILE} rmffile=${RMF_NAME} estepfile="sparse" 一方で、下記のように SKY 座標 (1-bin サイズ)で arf を作成すると、添付図赤のように、 エネルギー依存だけでなく、絶対値自体が3桁ちかくおかしな arf が生成されます。 $ cat reg/xis0_2005src.reg circle(796,996,250) -box(630,774,200,1500,355) -box(956,768,200,1500,355) $ xissimarfgen instrume=XIS${XISID} \ source_mode=SKYXY source_x=796 source_y=996 \ num_region=2 region_mode=SKYREG \ regfile1=${REG_NAME_SRC} regfile2=${REG_NAME_BGD} \ arffile1=${ARF_NAME_SRC} arffile2=${ARF_NAME_BGD} \ limit_mode=MIXED num_photon=${SIM_MUM_PHOTON} accuracy=0.005 \ phafile=${PHA_NAME} detmask=${CALDM_CALMASK} gtifile=none \ date_obs=${DATE_OBS} \ attitude=${ATTFILE} rmffile=${RMF_NAME} estepfile="sparse" これは、使いかたがなにか間違えているでしょうか。(単なる座標のミスならすいません) ログの最初のほうも添付しておきます。 1. DET_X,Y を使用した場合 ----------------------- xissimarfgen version 2006-11-26 Written by Y.ISHISAKI (TMU) .. SimASTE_Root: reading teldef file '/caldb//data/suzaku/xis/bcf/ae_xi0_teldef_20060125.fits' xissimarfgen: reading detmask '/caldb//data/suzaku/xis/bcf/ae_xi0_calmask_20060731.fits' TELESCOP= 'SUZAKU ' / Telescope (mission) name INSTRUME= 'XIS0 ' / Instrument name CTYPE1 = 'ACTX ' / Name of the X axis CTYPE2 = 'ACTY ' / Name of the Y axis CDELT1 = 1 / X-axis increment CDELT2 = 1 / Y-axis increment CRPIX1 = 1 / Position of the X reference pixel CRPIX2 = 1 / Position of the Y reference pixel CRVAL1 = 0 / Value of the X reference pixel CRVAL2 = 0 / Value of the Y reference pixel 10 keywords detected xissimarfgen: WARNING: image type is ACT, converting to DET 1024x1024 image read mask_ratio_ccd = 942793.0 / 1048576 = 0.899117 xissimarfgen: reading attitude 'att/aeaqr2005.att.gz' EA=( 310.0432 91.0615 185.1804 ) at t=184023556.0 [2005-10-30T21:39:16] EA=( 310.0432 91.0615 185.1804 ) at t=184023556.0 [2005-10-30T21:39:16] RA_NOM=310.0437, DEC_NOM=-0.9297 xis_contami: reading '/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 num_calc=55, ne_max=1000 xissimarfgen: making region maps regexpr=regfilter("reg_detxy/xis0_2005src_det.reg"), n_good_rows=63477 regexpr=regfilter("reg_detxy/xis0_2005bgd_det.reg"), n_good_rows=67406 xissimarfgen: *** show parameter *** POINTING AUTO REF_ALPHA (J2000) 310.0437 (20h40m10.488s) REF_DELTA (J2000) -0.9297 (-00d55m46.92s) REF_ROLL (deg) 0.0000 ABERRATION YES APERTURE_COSINE YES MINANGLE 0.000 (deg) MAXANGLE 360.000 (deg) MINRADIUS 57.96755 (mm) MAXRADIUS 200.10644 (mm) GEOMAREA 1152.4100 (cm2) SOURCE_MODE DETXY SOURCE_DETX,Y ( 283.0 , 500.0 ) REGION_MODE DETREG 1: REGFILE reg_detxy/xis0_2005src_det.reg ARFFILE arf/aeaqr2005_xis0src.arf 2: REGFILE reg_detxy/xis0_2005bgd_det.reg ARFFILE arf/aeaqr2005_xis0bgd.arf LIMIT_MODE MIXED NUM_PHOTON 500000 (photons) ACCURACY 0.005000 PHAFILE pha/aeaqr2005_xis0src.pha DETMASK /caldb//data/suzaku/xis/bcf/ae_xi0_calmask_20060731.fits GTIFILE none DATE_OBS '2005-10-30T21:39:16' (t=184023556.0) ATTITUDE att/aeaqr2005.att.gz RMFFILE rmf/aeaqr2005_xis0.rmf CONTAMIFILE /caldb//data/suzaku/xis/bcf/ae_xi0_contami_20061016.fits (CALDB) ESTEPFILE sparse 1: 0.200 - 1.000 keV : Fixed Energy in Each 0.200 keV 2: 1.000 - 7.000 keV : Fixed Energy in Each 0.500 keV 3: 7.000 - 8.000 keV : Fixed Energy in Each 1.000 keV 4: 8.000 - 16.000 keV : Fixed Energy in Each 2.000 keV 5: 0.250 - 0.300 keV : Fixed Energy in Each 0.050 keV 6: 2.100 - 2.400 keV : Fixed Energy in Each 0.050 keV 7: 2.200 - 2.250 keV : Fixed Energy in Each 0.010 keV 8: 2.300 - 2.350 keV : Fixed Energy in Each 0.010 keV 9: 2.600 - 2.800 keV : Fixed Energy in Each 0.020 keV 10: 3.100 - 3.300 keV : Fixed Energy in Each 0.100 keV 11: 11.800 - 12.000 keV : Fixed Energy in Each 0.100 keV CLOBBER YES SimASTE_XRTsim: *** show parameter *** SIMULATION_MODE Follow SimASTE_Root XRT_ID 'XRT-I0' MIRRORFILE '/caldb//data/suzaku/xrt/bcf/ae_xrt0_mirror_20060710.fits' (CALDB) MIRROR 'mirror' OBSTRUCT 'obstruct' QUADRANT 'quadrant' PCOL 'pcol' FOILTHICKNESS 0.17800000 MISSALIGNMODE 'gauss' MISSALIGNMENT 0.00800000 PCOLTHICKNESS 0.12000000 PCOLMISSALIGNMODE 'gauss' PCOLMISSALIGNMENT 0.00500000 REFLECTFILE '/caldb//data/suzaku/xrt/bcf/ae_xrta_reflect_20060710.fits' (CALDB) BACKPROFFILE '/caldb//data/suzaku/xrt/bcf/ae_xrta_backprof_20060719.fits' (CALDB) SCATTER ASTROE NvfSw 1.00000000 GauSigma 0.11700000 Ngau 1.02710000 ExpSigma 0.56970000 Nexp 0.73690000 LorGamma 6.46200000 Nlor 0.00000000 ScatSw 0.00000000 S3 0.00000000 C3 0.00000000 SHIELDFILE '/caldb//data/suzaku/xrt/bcf/ae_xrta_shield_20051129.fits' (CALDB) SimASTE_XRTsim: reading mirror & reflection file ... Reading '/caldb//data/suzaku/xrt/bcf/ae_xrt0_mirror_20060710.fits' ... Reading '/caldb//data/suzaku/xrt/bcf/ae_xrta_reflect_20060710.fits' ... NAXIS2=885, binAngleNum=4000, binAngleZero=0.000000, binAngleDelta=1.745320e-05 NAXIS2=107, binAngleNum=4000, binAngleZero=0.000000, binAngleDelta=1.745320e-05 NAXIS2=36, binAngleNum=4000, binAngleZero=0.000000, binAngleDelta=1.745320e-05 Reading '/caldb//data/suzaku/xrt/bcf/ae_xrta_backprof_20060719.fits' ... Reading '/caldb//data/suzaku/xrt/bcf/ae_xrta_shield_20051129.fits' ... finished xissimarfgen: creating ARF files ... Creating arf/aeaqr2005_xis0src.arf Creating arf/aeaqr2005_xis0bgd.arf 2. SYY_X,Y を使用した場合 ----------------------- xissimarfgen version 2006-11-26 Written by Y.ISHISAKI (TMU) .. SimASTE_Root: reading teldef file '/caldb//data/suzaku/xis/bcf/ae_xi0_teldef_20060125.fits' xissimarfgen: reading detmask '/caldb//data/suzaku/xis/bcf/ae_xi0_calmask_20060731.fits' TELESCOP= 'SUZAKU ' / Telescope (mission) name INSTRUME= 'XIS0 ' / Instrument name CTYPE1 = 'ACTX ' / Name of the X axis CTYPE2 = 'ACTY ' / Name of the Y axis CDELT1 = 1 / X-axis increment CDELT2 = 1 / Y-axis increment CRPIX1 = 1 / Position of the X reference pixel CRPIX2 = 1 / Position of the Y reference pixel CRVAL1 = 0 / Value of the X reference pixel CRVAL2 = 0 / Value of the Y reference pixel 10 keywords detected xissimarfgen: WARNING: image type is ACT, converting to DET 1024x1024 image read mask_ratio_ccd = 942793.0 / 1048576 = 0.899117 xissimarfgen: reading attitude 'att/aeaqr2005.att.gz' EA=( 310.0432 91.0615 185.1804 ) at t=184023556.0 [2005-10-30T21:39:16] EA=( 310.0432 91.0615 185.1804 ) at t=184023556.0 [2005-10-30T21:39:16] RA_NOM=310.0437, DEC_NOM=-0.9297 xis_contami: reading '/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 num_calc=55, ne_max=1000 xissimarfgen: making region maps regexpr=regfilter("reg/xis0_2005src.reg"), n_good_rows=62049 regexpr=regfilter("reg/xis0_2005bgd.reg"), n_good_rows=63756 xissimarfgen: *** show parameter *** POINTING AUTO REF_ALPHA (J2000) 310.0437 (20h40m10.488s) REF_DELTA (J2000) -0.9297 (-00d55m46.92s) REF_ROLL (deg) 0.0000 ABERRATION YES APERTURE_COSINE YES MINANGLE 0.000 (deg) MAXANGLE 360.000 (deg) MINRADIUS 57.96755 (mm) MAXRADIUS 200.10644 (mm) GEOMAREA 1152.4100 (cm2) SOURCE_MODE SKYXY SOURCE_X,Y ( 796.0 , 996.0 ) REGION_MODE SKYREG 1: REGFILE reg/xis0_2005src.reg ARFFILE arf/aeaqr2005_xis0src.arf 2: REGFILE reg/xis0_2005bgd.reg ARFFILE arf/aeaqr2005_xis0bgd.arf LIMIT_MODE MIXED NUM_PHOTON 50000 (photons) ACCURACY 0.005000 PHAFILE pha/aeaqr2005_xis0src.pha DETMASK /caldb//data/suzaku/xis/bcf/ae_xi0_calmask_20060731.fits GTIFILE none DATE_OBS '2005-10-30T21:39:16' (t=184023556.0) ATTITUDE att/aeaqr2005.att.gz RMFFILE rmf/aeaqr2005_xis0.rmf CONTAMIFILE /caldb//data/suzaku/xis/bcf/ae_xi0_contami_20061016.fits (CALDB) ESTEPFILE sparse 1: 0.200 - 1.000 keV : Fixed Energy in Each 0.200 keV 2: 1.000 - 7.000 keV : Fixed Energy in Each 0.500 keV 3: 7.000 - 8.000 keV : Fixed Energy in Each 1.000 keV 4: 8.000 - 16.000 keV : Fixed Energy in Each 2.000 keV 5: 0.250 - 0.300 keV : Fixed Energy in Each 0.050 keV 6: 2.100 - 2.400 keV : Fixed Energy in Each 0.050 keV 7: 2.200 - 2.250 keV : Fixed Energy in Each 0.010 keV 8: 2.300 - 2.350 keV : Fixed Energy in Each 0.010 keV 9: 2.600 - 2.800 keV : Fixed Energy in Each 0.020 keV 10: 3.100 - 3.300 keV : Fixed Energy in Each 0.100 keV 11: 11.800 - 12.000 keV : Fixed Energy in Each 0.100 keV CLOBBER YES SimASTE_XRTsim: *** show parameter *** SIMULATION_MODE Follow SimASTE_Root XRT_ID 'XRT-I0' MIRRORFILE '/caldb//data/suzaku/xrt/bcf/ae_xrt0_mirror_20060710.fits' (CALDB) MIRROR 'mirror' OBSTRUCT 'obstruct' QUADRANT 'quadrant' PCOL 'pcol' FOILTHICKNESS 0.17800000 MISSALIGNMODE 'gauss' MISSALIGNMENT 0.00800000 PCOLTHICKNESS 0.12000000 PCOLMISSALIGNMODE 'gauss' PCOLMISSALIGNMENT 0.00500000 REFLECTFILE '/caldb//data/suzaku/xrt/bcf/ae_xrta_reflect_20060710.fits' (CALDB) BACKPROFFILE '/caldb//data/suzaku/xrt/bcf/ae_xrta_backprof_20060719.fits' (CALDB) SCATTER ASTROE NvfSw 1.00000000 GauSigma 0.11700000 Ngau 1.02710000 ExpSigma 0.56970000 Nexp 0.73690000 LorGamma 6.46200000 Nlor 0.00000000 ScatSw 0.00000000 S3 0.00000000 C3 0.00000000 SHIELDFILE '/caldb//data/suzaku/xrt/bcf/ae_xrta_shield_20051129.fits' (CALDB) SimASTE_XRTsim: reading mirror & reflection file ... Reading '/caldb//data/suzaku/xrt/bcf/ae_xrt0_mirror_20060710.fits' ... Reading '/caldb//data/suzaku/xrt/bcf/ae_xrta_reflect_20060710.fits' ... NAXIS2=885, binAngleNum=4000, binAngleZero=0.000000, binAngleDelta=1.745320e-05 NAXIS2=107, binAngleNum=4000, binAngleZero=0.000000, binAngleDelta=1.745320e-05 NAXIS2=36, binAngleNum=4000, binAngleZero=0.000000, binAngleDelta=1.745320e-05 Reading '/caldb//data/suzaku/xrt/bcf/ae_xrta_backprof_20060719.fits' ... Reading '/caldb//data/suzaku/xrt/bcf/ae_xrta_shield_20051129.fits' ... finished xissimarfgen: creating ARF files ... Creating arf/aeaqr2005_xis0src.arf Creating arf/aeaqr2005_xis0bgd.arf 以上、どうぞよろしくおねがいします。 
Answer from(回答の文責)ソフトチーム 石崎Answer (答え)見てみたところ、姿勢ファイルのヘッダから自動取得される skyref が> REF_ALPHA (J2000) 310.0437 (20h40m10.488s) > REF_DELTA (J2000) -0.9297 (-00d55m46.92s) > REF_ROLL (deg) 0.0000 姿勢ファイルの、指定時刻 2005-10-30T21:39:16 におけるオイラー角が > xissimarfgen: reading attitude 'att/aeaqr2005.att.gz' > EA=( 310.0432 91.0615 185.1804 ) at t=184023556.0 [2005-10-30T21:39:16] となっていて、90 - EA2 = -1.0615 と REF_DELTA の値と 0.13 deg ほど食い違っています。 指定した時刻においては、姿勢が安定してなくて、region ファイルで 設定した領域内にソースが入ってないのだと思います。 作成した ARF の primary extension に入っている画像を ds9 で見れば、設定 region 内に正しくソースが入っているか どうかは簡単にチェックできますので、確認してみて下さい。 Status （詳細なステータス）2007-01-31 Accept 2007-02-01 Done
|