photon event fileからはじめた方で、unscreendのイベントファイルを使っている場合、以下のmake obscatと chooseを、またscreendのイベントファイルを使っている場合は、さらにsiscleanというコマンドまで実行する必要はありません。
以下のようにして``bright''モードのfileのみを選びます
![[*]](foot_motif.png)
。
さらに、今の場合は HIGH bit rate のみを選んでいます。
CasA_s0:ASCA-SIS0 > make obscat
> Enter the Event file dir >[ ] ../frf
> Default filter for the catalogue (NONE for none, DEF for default) >[DEF]
:
CasA_s0:ASCA-SIS0 > select obscat "DATAMODE=='BRIGHT'&& BIT_RATE=='HIGH'"
CasA_s0:ASCA-SIS0 > show obscat
******************** Observation Catalogue ********************
Data Directory is: /tmp_mnt/uzuki/d2/sonobe/frf/
HK Directory is: /tmp_mnt/uzuki/d2/sonobe/frf/
OBJECT DATAMODE BIT_RATE DATE-OBS TIME-OBS ONTIME NEVENTS S0CCDPOW
1 Cas A BRIGHT HIGH 01/08/93 13:55:19 .316E+04 152782 0100
2 Cas A BRIGHT HIGH 01/08/93 15:38:15 .676E+03 40193 0100
3 Cas A BRIGHT HIGH 01/08/93 15:50:51 .247E+04 143852 0100
4 Cas A BRIGHT HIGH 01/08/93 17:06:35 .147E+04 83547 0100
5 Cas A BRIGHT HIGH 01/08/93 17:33:27 .209E+04 119642 0100
6 Cas A BRIGHT HIGH 01/08/93 18:44:11 .349E+04 200812 0100
7 Cas A BRIGHT HIGH 01/08/93 20:22:35 .338E+04 193444 0100
8 Cas A BRIGHT HIGH 01/08/93 21:59:23 .267E+04 154151 0100
CasA_s0:ASCA-SIS0 > choose 1-**
最初に、ホットピクセルを取り除きましょう。
これには``sisclean''を使います。
これは、ある一定の閾値を越えたものをホットピクセル
として取り除きます。
CasA_s0:ASCA-SIS0-BRIGHT > sisclean
> Which method (1=sisclean, 2=cleansis)? >[1] 2 ←(￥)
> Clean cell size for local background estimation (odd integer > 1)>[5]
> Enter the log poission prob for clean threshold (must be < 0) >[-5.24]
> Zero background threshold >[3]
> Lower cutoff for hot pixel search >[0]
> Upper cutoff for hot pixel search >[4095]
-----------------------------------------
Making EVENTS list for cleansis
:
CLEANSIS_V1.3
allocating image arrays...
PROGRAM TO MAKE AN SIS SCIENCE FILE CLEANED OF ANOMALOUS PIXELS.
* Anomalous pixels may consist of at least two populations.
1. Persistent HOT pixels are removed by comparing to the chip mean.
2. Flickering WARM pixels are removed by comparing to the cell mean.
3. In faint areas (zero bgd), flickering pixels are removed with a
cutoff threshold.
Hints:
Choose cell size and thresholds based on the expected background and the PSF.
Use the DIRTYSIS option to examine the anomalous pixel spectra.
Try a multiple pass clean:
Choose a PHA cut to optimize the S/N of the flickering pixels,
then a broad band clean.
Be suspicious of extended source cleans.
For very bright sources you may need to turn off the iteration option.
See the help page for further info (fhelp cleansis)
:
:
updating NEVENTS keywords...
closing data file...
closing clean file...
CasA_s0:ASCA-SIS0-BRIGHT >
(￥)で``2''を選びます。以下のパラメタは多くの場合デフォルトでいいと思います。
ホットピクセルについての説明もでます。良く読んで下さい。
この方法でホットピクセルを取り除いた events file を自動的に
``CasA_s0_in_event.xsl''という名前で作り、以後これにアクセスします。
ここまで来たらあとはGISの時と同じようにできます。 すなわち、
CasA_s0:ASCA-SIS0-BRIGHT > extract image
:
CasA_s0:ASCA-SIS0-BRIGHT > saoimage
CasA_s0:ASCA-SIS0-BRIGHT > save image s0.image.fits
CasA_s0:ASCA-SIS0-BRIGHT > extract curve
:
CasA_s0:ASCA-SIS0-BRIGHT > plot curve
:
CasA_s0:ASCA-SIS0-BRIGHT > save curve s0.curve.fits
:
etc.
なお、SIS ではデフォルトではイメージは 4x4 が一つのピクセルに
なって表示されます。
それを最小単位のピクセル
にするには
CasA_s0:ASCA-SIS0-BRIGHT > set xybinsize 1とします。
スペクトルをセーブするときには、
CasA_s0:ASCA-SIS0-BRIGHT > save spec CasA_sis0.pi
The data will be rebinned up to 512 with the BRIGHT2LINEAR binning.
The resultant spectrum will have the low 13 channels flagged as bad.
> Group ( or rebin ) the spectra before outputting? >[yes] ←(∞)
:
Wrote spectrum to CasA_sis0.pi
となります。(∞)で no と答えると 2048 チャンネル変則binnig でセーブされます。
yes とこたえると RBNPHA
という FTOOLS が働いて linear な 512 チャンネルの
スペクトルがセーブされます。