データの読み込み

まず 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 >

のようになります。
Enter the Event file dir $>$[ ]
の[ ]の中はdefaultのディレクトリで、 それでよければそのままリターンをします。 違うときだけ正しいディレクトリを指定します 。
Default filter for the catalogue (NONE for none, DEF for default) $>$[DEF]
``DEF''では一つの science file について、 ONTIME(exposure)が100秒以上、 イベント数が 1 以上、HV_RED が off, HVH の値が 3, HVL の値が 4 のファイルだけ 選びます。

``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) にその一覧を示します。

表 7.1: XSELECT で有効な演算記号

;;	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 filter

XSELECT全体のステイタスを見るには

  CasA_g2:ASCA-GIS2-PH > show status

で見ることができます。