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と打ちます。