Suzaku RIKEN Home
Question ID = 0052

Assigned Member (担当者)

2007-04-19 Y.Terada, continue from Q=0049.

Name and Institute (質問者/機関)

尾崎 正伸 / ISAS

Question (質問)

ありがとうございます。このように変更し、その後色々と試行錯誤して以下の状態まで辿り着きました。
  • HK フィルタリング条件は、HK ファイル毎に分けて記述しなければならない。
    その際には、今まで読んでいたイベントおよび HK ファイルを一旦破棄して、改めて直前までフィルタしたイベントおよび新しい HK ファイルを読み込んで処理を継続する必要がある。(でないと、xselect が古い HK と新しい HK をマージしようとして、両者のカラム構成が一緒でない限り失敗する。)
    つまり、
    (前略)
    xsel:SUZAKU-XIS0-STANDARD > filter column "GRADE=0:0 2:4 6:6 STATUS=0:262143"
    xsel:SUZAKU-XIS0-STANDARD > clear hk
    xsel:SUZAKU-XIS0-STANDARD > read hk
    > HK file directory >[./] ./
    > Enter HK file list >[ae100020010.hk.gz] ae100020010.ehk.gz
    > Are the HK files in expanded form? >[yes] yes
    xsel:SUZAKU-XIS0-STANDARD > select hk "SAA_HXD == 0 && T_SAA_HXD > 436 && ELV > -5 && NTE_ELV>100 && ANG_DIST < 1.5"
    xsel:SUZAKU-XIS0-STANDARD > bin event
    xsel:SUZAKU-XIS0-STANDARD > save event tmp1.evt clobber=yes
    > Use filtered events as input data file ? >[no] no
    xsel:SUZAKU-XIS0-STANDARD > clear data
    xsel:SUZAKU-XIS0-STANDARD > read event tmp1.evt
    xsel:SUZAKU-XIS0-STANDARD >
    xsel:SUZAKU-XIS0-STANDARD > clear hk
    xsel:SUZAKU-XIS0-STANDARD > read hk
    > HK file directory >[./] ./
    > Enter HK file list >[ae100020010.ehk.gz] ae100020010xi0_0.hk.gz
    > Are the HK files in expanded form? >[yes] yes
    xsel:SUZAKU-XIS0-STANDARD > select hk "S0_ADHST_A == 0 && S0_ADHST_B == 0 && S 0_ADHST_C == 0 && S0_ADHST_D == 0"
    (以下略)
    のようにする。
しかしここで、別の問題にぶつかりました。フィルタ条件で使うカラムがHK ファイルの 2nd extension 以降にある場合、xselect (から呼ばれるmaketime) がそのカラムを認識してくれないようです。どのようにすれば、任意の extension のカラムを認識させられるでしょうか。 末尾に
実行例のログと、 処理に失敗している HK ファイルの fstruct結果を添付します。

※ 質問を言い替えるとこんな感じ(En)

Answer from(回答の文責)

Keith Arnaud / FTOOLS HELP

Answer (答え)

I've looked into this and I don't see an easy way of changing xselect to make it work. The problem is that xselect has always assumed that the HK data is in the first extension of the file. The easiest workaround would be to make a new HK file by copying the fifth extension eg
ftcopy ae400001020xi0_0.hk+5 ae400001020xi0_0.hk5 copyall=NO
Then read ae400001020xi0_0.hk5 as the hk file.

この不具合は、Keith さんが対処されたようですので、 次の HEADAS release (6.2.x)に取り込まれることでしょう。

Status (詳細なステータス)

2007-04-11 expanded from Q=0049
2007-04-19 Accept
2007-04-19 Circulate to FTOOLS help, Dr.Keith Arnaud.
2007-04-21 Got Answer from Keith.(1)
2007-04-23 Done.
2007-04-27 Got Answer from Keith.(2)
2007-04-27 Done(2).

Suzaku Official Page Suzaku Help Desk Japan Suzaku Help Desk GOF Suzaku RIKEN Page RIKEN Makishima Cosmic Radiation Laboratory