Filter file をつかう

前に``mkfilter2''をつかって衛星の環境パラメタを収めたファイルをつくりました 。 そのファイル名を仮に``ft930801_1340_2309.mkf''としましょう。 ``fdump''コマンドでft930801_1340_2309.mkfを見ると ``keyword''がわかります。その``keyword''はfilterの中でつかえます。 例えばelevationが5度以上で SAA 中でないデータがほしい時はつぎのようにします 。

   CasA_g2:ASCA-GIS2-PH > select mkf
   > Boolean expression for filter file selection >[ ] ELV_MIN>5.0 && SAA==0 
   > Enter the filter file directory >[ ] .

セレクションの条件はブール代数で表現します。 いくつか条件を並べるときは.and.(または &&)で並べます(表(7.1))。 ``Enter the filter file directory''で mkf のある directory を指定します。 XSELECTはこのときft*.mkfというファイルを探して読み込みます 。 このあと実行される``extract curve''etc.にはこの条件が満たされているもののみ使われます。 セレクションをクリアしたい時はつぎのようにします。

   CasA_g2:ASCA-GIS2-PH > clear mkfsel

ゴダードからscreenedのイベントファイル(rev2)を持ってきた場合、すでに以下のfilterがかけられています。

   GIS
     ACS==0 &&ANG_DIST>0 &&ANG_DIST<0.01&&ELV>5 &&
     SAA==0&&COR>4 &&G2_H0+G2_H2+G3_H0+G3_H2<45 &&
     G2_H0+G2_H2+G3_H0+G3_H2<0.45*COR**2-13*COR+125 &&RBM_CONT<100 &&
     G[2,3]_L1>0 

   SIS
    SAA==0 &&BR_EARTH>20 &&ELV>10 &&COR>6 &&ACS==0 &&ANG_DIST>0 &&
    ((T_SAA  <0)||(T_SAA  >*)) && ANG_DIST<0.01 &&FOV==0 && 
    ((T_DY_NT<0)||(T_DY_NT>*)) && 
    ((S[0,1]_PIXL[0,1,2,3]>*)   && (S[0,1]_PIXL[0,1,2,3]<*)) &&
    ((S[0,1]_SATF[0,1,2,3]<1)||(S[0,1]_SATF[0,1,2,3]>1))

SISの*と書かれた部分は、S0、S1、又は4CCDモードか1CCDモードかなどにより、値が異なっています。詳しくは、データを取ってきたdirectoryの''aux/ad********_***_joblog.html.''というファイルに、イベントファイルを処理した時のログが残っているので、これで確認してみて下さい。

S[0,1]_SATF[0,1,2,3]は、saturation flagで、mkfilterの古いバージョンで作ったmkfファイルには、このキーワードは含まれていないようです。

mkfファイル中の役に立ちそうな項目と標準的なセレクション示しておきます。詳細は、英語版マニュアル''ASCA ABC Guide''を参照して下さい。

      KEYWORD               DESCRIPTION                   STANDARD SETTING (and Units)    
      --------------------------------------------------------------------------------    
      SAA       Passage through South Atlantic Anomaly  0 (for no)    
      T_SAA     Time after SAA passage                  > 16 (for 1CCD mode)
                                                        > 64 (for 4CCD mode)
                                                        (4 readout cycles)
      FOV       Field of View                           0 (for Sky)
      COR_MIN   Minimum value of Cut Off Rigidity       >4 for GIS, 6 for SIS (GeV/c)
      ELV_MIN   Miminum elevation angle from Earth      >5-10 (deg)
      BR_EARTH  Elevation angle from bright Earth       >15-40 for SIS0 4CCD (deg)
                                                        >15-20 for SIS1 4CCD (deg)
                                                        >10-15 for SIS[0,1] [1,2]CCD (deg)
      NSAS      Sun angle                               >0 means satellite day
      G2_L1     GIS2 events within discriminators       >0 means good data
      G3_L1     GIS3 events within discriminators       >0 means good data
                                                        <0 means satellite night
      T_DY_NT   Time after day/night transition         > 16 (for 1CCD mode)
                                                        > 64 (for 4CCD mode)
                                                        (4 readout cycles)
      HV_RED    GIS high voltage reduction              'OFF'
      Z_ALPHA   RA of satellite z-axis                  (exclude fluctuating values)
      Z_DEC     dec of satellite z-axis                 (exclude fluctuating values)
      S0_ID     List of active chips for SIS0            0123 (for 4-CCD mode)
      S1_ID     List of active chips for SIS1            0123 (for 4-CCD mode)
      Sn_PIXLm  Number of SIS Events above threshold    >0 always
                (n = sensor number, m = chip number)    <400 in 1-CCD mode
                                                        <600 in 2-CCD mode
                                                        <800 in 4-CCD mode
      Sn_EVNTm  SIS_n chip_m HK Event number
      Sn_TELMm  SIS_n chip_m telemely event number
      Gn_LDHIT  GIS_n Lower Discri Hit counting rate
      SIS_RBMF  Radiation belt monitor  0:good/1:bad    0 (for good)