データ・レスポンスの読み込み

先ほどの``XSELECT''で作ったスペクトルファイルを読み込みます。

   XSPEC> data CasA_gis2.pi
    Net count rate (cts/cm^2/s) for file   1   45.29    +/-  5.0208E-02

間違ったファイルを指定するとエラーが出ます。

次に、レスポンス を読み込ます。 RMF 名を仮に response.rmf、 ARF を eff.arf とします。

   XSPEC> resp response.rmf
   XSPEC> arf eff.arf

XSELECTによって生成されるGISのスペクトルファイルは通常1024channels ですので、 レスポンスもこれに対応したものでなければなりません。

読み込んだデータの内容を確認するには``show''を使います。

  XSPEC> show all
   00:04:48 25-Sep-94
   Fit statistic in use is Chi-Squared
   Minimization technique is Lev-Marq
   Convergence criterion =   1.0000000000000000E-02
   Querying enabled
   Prefit-renorming enabled
  
   Information for file   1
    belonging to plot group   1, data group   1, det id = ASCA GIS2
    Current data file: CasA_gis2.pi
    No current background
    No current correction
    Response (RMF) file    : gis2v3_1.rmf
    Auxiliary (ARF) file    : CasA_gis2.arf
    Noticed channels   1 to  256
    File observed count rate     45.29    +/-5.02082E-02 cts/cm^2/s
                                 45.29    +/-5.02082E-02 cts/s
    After correction of   0.0000E+00;  Model predicted rate:  0.0000E+00

とりあえず、読み込んだスペクトルを見たいのならば、

  XSPEC> plot    または   XSPEC> iplot

とします。 ``plot''では、ただプロットするだけですが、 ``iplot''では、PLT$>$のプロンプトがあらわれ、qdpのコマンドが使えます。 これだと横軸が channel になっています。 横軸をエネルギーにしたい時は、

  XSPEC> setplot energy

とします。channel に戻したいときは、``setplot channel''です。

実は``plot''コマンドにはいくつもの種類があります。以下のようなものです。

  XSPEC> plot ?
  Choose from the following `plot' sub-commands:
 data      counts    ldata     residuals chisq     delchi    ratio       
 summary   model     emodel    eemodel   contour   efficien  ufspec      
 eufspec   eeufspec  dem       insensitv sensitvty genetic   foldmodel   
Insert selection: (data)

例えば、縦軸を log スケールにしたいときは、

  XSPEC> plot ldata

とします。一度、``plot ***''を実行すると、次回以降``plot''とするだけで、 ``plot ***''と同じことになります。

このままだとエネルギーのかなり低いところまで表示されていると思います。 いらないデータ点( bin という)を無視するために以下のようにします。

  XSPEC> ignore 1-50 881-**

1ch から 50ch と 881ch 以上を除外したことになります。 いったん ignore した bin をもとに戻すには、

  XSPEC> notice 40-50  
   A total of    11 channels will be noticed
   Net count rate (cts/cm^2/s) for file   1   44.88    +/-  4.9982E-02
  XSPEC> plot

とすれば、40 から 50 bin が再び表示されます。 また、エネルギーで指定したい時には、実数で指定します。例えば、1keV以下のエネルギーを無視するには、

  XSPEC> ignore **-1.0

と打ちます。