まず Observation Catalogue を作ります。
CasA_g2:ASCA-GIS2 > make obscat > Enter the Event file dir >[ ] ../frf > Default filter for the catalogue (NONE for none, DEF for default) >[DEF] Using default selection expression: ONTIME>100&&NEVENTS>0&&HV_RED=='OFF' &&HVH_LVL==3&&HVL_LVL==4 ******************** 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 PHA_BINS 1 Cas A PH MEDIUM 01/08/93 13:50:53 .272E+03 196 1024 2 Cas A PH HIGH 01/08/93 13:55:31 .316E+04 162710 1024 3 Cas A PH MEDIUM 01/08/93 14:48:45 .304E+03 4864 1024 4 Cas A PH MEDIUM 01/08/93 15:31:09 .432E+03 6697 1024 5 Cas A PH HIGH 01/08/93 15:38:27 .672E+03 35768 1024 6 Cas A PH HIGH 01/08/93 15:52:13 .239E+04 124755 1024 7 Cas A PH HIGH 01/08/93 17:09:35 .130E+04 68990 1024 8 Cas A PH HIGH 01/08/93 17:33:33 .210E+04 107963 1024 9 Cas A PH MEDIUM 01/08/93 18:09:49 .320E+03 1445 1024 10 Cas A PH HIGH 01/08/93 18:47:43 .330E+04 172338 1024 11 Cas A PH MEDIUM 01/08/93 19:45:49 .672E+03 1897 1024 12 Cas A PH HIGH 01/08/93 20:25:51 .320E+04 166390 1024 13 Cas A PH MEDIUM 01/08/93 21:21:49 .102E+04 2102 1024 14 Cas A PH HIGH 01/08/93 22:05:03 .234E+04 125024 1024 15 Cas A PH MEDIUM 01/08/93 22:47:09 .656E+03 8715 1024 16 Cas A PH LOW 01/08/93 22:58:05 .128E+03 370 1024 17 Cas A PH LOW 01/08/93 23:00:13 .576E+03 1378 1024 CasA_g2:ASCA-GIS2 >のようになります。
``LOW'' bit rate のデータがいらない場合は、
(``LOW'' bit rate のデータは、データがサチっているので、使うことはほとんどありません。取り除いておきましょう)
CasA_g2:ASCA-GIS2 > select obscat > Enter catalogue filter >[ ] BIT_RATE.ne.'LOW' CasA_g2:ASCA-GIS2 > 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 PHA_BINS 1 Cas A PH MEDIUM 01/08/93 13:50:53 .272E+03 196 1024 2 Cas A PH HIGH 01/08/93 13:55:31 .316E+04 162710 1024 3 Cas A PH MEDIUM 01/08/93 14:48:45 .304E+03 4864 1024 4 Cas A PH MEDIUM 01/08/93 15:31:09 .432E+03 6697 1024 5 Cas A PH HIGH 01/08/93 15:38:27 .672E+03 35768 1024 6 Cas A PH HIGH 01/08/93 15:52:13 .239E+04 124755 1024 7 Cas A PH HIGH 01/08/93 17:09:35 .130E+04 68990 1024 8 Cas A PH HIGH 01/08/93 17:33:33 .210E+04 107963 1024 9 Cas A PH MEDIUM 01/08/93 18:09:49 .320E+03 1445 1024 10 Cas A PH HIGH 01/08/93 18:47:43 .330E+04 172338 1024 11 Cas A PH MEDIUM 01/08/93 19:45:49 .672E+03 1897 1024 12 Cas A PH HIGH 01/08/93 20:25:51 .320E+04 166390 1024 13 Cas A PH MEDIUM 01/08/93 21:21:49 .102E+04 2102 1024 14 Cas A PH HIGH 01/08/93 22:05:03 .234E+04 125024 1024 15 Cas A PH MEDIUM 01/08/93 22:47:09 .656E+03 8715 1024このようにしてObservation Catalogueにフィルターをかけます。
xselect ではいろいろな条件でセレクションをするのにブール代数の表記を使います。 Fortran 型と C 型の両方が使えます。表(7.1) にその一覧を示します。
;; | Fortran type | C type | ;; | Fortran type | C type |
equal | .eq. | == | not equal | .ne. | != |
less than | .lt. | < | less than or equal | .le. | ≦ |
greater than | .gt. | > | greater than or equal | .ge. | ≧ |
or | .or. | || | and | .and. | && |
negation | .not. | ! | ;; | ;; | ;; |
この中から、必要なデータを選びます。
CasA_g2:ASCA-GIS2 > choose 2,5-9すると、選ばれたファイルが表示されます。 なお、連続した番号であったら5-9のように記述できます。 すべてを選びたかったら
CasA_g2:ASCA-GIS2 > choose 1-** Setting datamode to PH Got the minimum time resolution of the chosen data: 0.50000E+00, Getting Min and Max for Energy Column... Got min and max for PI: 0 1023 Number of files read in: 15 Files currently in use: 1 ft930801_1340_2309G200370M.fits 2 ft930801_1340_2309G200770H.fits 3 ft930801_1340_2309G200970M.fits 4 ft930801_1340_2309G201370M.fits 5 ft930801_1340_2309G201770H.fits 6 ft930801_1340_2309G202570H.fits 7 ft930801_1340_2309G203570H.fits 8 ft930801_1340_2309G203770H.fits 9 ft930801_1340_2309G203970M.fits 10 ft930801_1340_2309G204870H.fits 11 ft930801_1340_2309G205070M.fits 12 ft930801_1340_2309G205970H.fits 13 ft930801_1340_2309G206170M.fits 14 ft930801_1340_2309G207070H.fits 15 ft930801_1340_2309G207270M.fits CasA_g2:ASCA-GIS2-PH >とします。choose をすると、データのモード、時間分解能、PHA や PI のコラム数も表示します。 そして、プロンプトにデータモードを表示します。
``make obscat'' を実行するまえに、``set datamode PH''を実行しておくと、 ``make obscat'' をしたときに、PH モードのデータのみ選ばれます。
一連の作業の中でObservation Catalogueが見たくなったら、
CasA_g2:ASCA-GIS2-PH > show obscat現在選び出されているファイルを知りたかったら、
CasA_g2:ASCA-GIS2-PH > show data現在どのようなフィルターがかけられ、セレクションが行なわれているかをみるには
CasA_g2:ASCA-GIS2-PH > show filterXSELECT全体のステイタスを見るには
CasA_g2:ASCA-GIS2-PH > show statusで見ることができます。