... 番目のX線天文衛星[*]
宇宙研の第 15 号科学衛星でもある
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 日[*]
当初の予定は 2/12 であったが 自然にやさしいロケットの都合により 8 日延期された
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... mです[*]
打ち上げ時は約 2.8 m である
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 分角(HPD[*]
半分のX線が集められる円の直径
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... SIS[*]
Solidstate Imaging Spectrometer
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... を達成している[*]
打ちあげ当初の性能であり、2000年現在では 5%くらい。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... パイルアップ、CTE[*]
Charge Transfer Efficiency:電荷転送率
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... これらについては、ここでは省略するが理解しておくことが望ましい[*]
これらのことは、小谷修論 ``ASTRO-D 搭載 SIS の性能評価''1993 東京大学、山下修論''「あすか」搭載X線CCDカメラの軌道上での較正''1995 東京大学 にくわしい
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... GIS[*]
Gas Imaging Spectrometer;撮像型蛍光比例計数管
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 観測モードには、PH[*]
Pulse Height
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... PCAL[*]
Position Calibration
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... MPC[*]
Multi-channel Pulse Count
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... そして、そのファイル[*]
FRF file、またはphoton event file
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... fileからデータ解析を行ないたい方は、frfread[*]
現在は、ちょっとした理由により、宇宙研からしか取ってこれません。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... (f2b、f2b2)[*]
f2b、f2b2は、ゴダードの向井さん製のshell scriptで、2〜3のftoolsを自動的に実行してくれます。このプログラムがなくても、ftoolsさえあれば、解析は行なえますが、あると便利です。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... caldbというパッケージ[*]
このパッケージには、あすか以外の衛星のレスポンス等も含まれています。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... caldbで、基本的な解析で使う最低限のファイルは、上に書いたようなファイルですが、ascaのdirectoryにあるデータは、全て持ってきておいた方がいいかもしれません。ただ、他の衛星のデータ(ROSATなど、、、)も含めた全てのファイルを持ってくるには、かなりのdisk容量が必要となります[*]
directoryごとftpでファイルを持ってくるには、ftp$>$get asca.tarなどと打つと、tarでまとめられたdirectoryを得ることができて便利です。ただし、ftpのバージョンによっては、できないものもあります。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
...[*][*]
HEAsoft5.0以降、GSFCのcaldbをリモートで参照できるようになりましたが、DEC/Alphaでは、まだ使えない?ようです。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... あすかのデータの、標準的な解析toolは、FTOOLS+XANADUです。最近、この2つが合体し、HEAsoftなるものができました。管理者の方にインストールさえしていただければ、.cshrcに以下の2行を書き加えるだけでOKです[*]
宇宙研では、アカウントをもらった時点で、自分の.cshrcの中に、source /usr/local/template/CSHRC という行があり、これらの解析ソフトの初期設定が自動的に行なわれます。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ftpでデータを得ることができます。たどり着くまでのdirectory名が少し違いますが、両者は基本的に同じものです[*]
PLAINセンターのデータは、NASAゴダードからAnonymous ftpで持ってきているので、リリースは少し遅れます。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... FRFからはじめる方は、fa******\_****.****、ft******\_****.****、frf.orbitというファイルが必要となります。以下で、()内は、ゴダードから持ってくる場合です[*]
ここでも、directoryごとftpでファイルを持ってくるには、ftp$>$get aux.tarなどと打つとauxというdirectoryごと、持ってくることができます。ftp$>$get aux.tar.gzも可能ですが、すでに中のデータは圧縮されいるので、あまりメリットはありません。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... FRFREADは終了までに数十分ほどかかります[*]
状況により異なります
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
...program-file節参照)。数カ月に1度、更新されるので、できれば最新のものを使いましょう[*]
解析するデータの観測された日より古いgis_temp2gain.fitsを使うと、ゲインが正しく補正されません。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... と、``.''が``\_''になります[*]
VAX/VMS では、.*** を拡張子と解釈してしまうため
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... を読みとって[*]
bright mode のファイルがない場合は 40 の値がつかわれます
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... というファイルを新しくつくります[*]
古いバージョンではft*S*11H.fitsというファイルを作ってました。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... の補正も行ないます[*]
ascatimeから計算して、補正される
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... にあります[*]
'00/06/22 現在
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ``fplot''[*]
``便利な FTOOLS'' の章参照
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... をとおしましょう[*]
ascalin がとおしてなかったことからくるエラーがよくあります
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
...。``fdump'' [*]
``便利な FTOOLS'' の章参照
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... frfread[*]
y2k問題にも対応している、frfread4.00 以上を使うこと
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... に正しい座標を書き込みます[*]
telemetry から出てくる science file には生座標(RAWX, RAWY)がついています。 それをもとに検出器座標(DETX, DETY)と sky 座標(X, Y)を書き込みます
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... のゲインの補正も行ないます[*]
telemetry から出てくる PHA の情報とゲインヒストリーをもとにゲイン補正された PI コラムを書きだします
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... PI[*]
Pulse Invariant
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... file[*]
temp2gain でつくったファイルです
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... で自動的にしかるべきファイルを読みにいきます[*]
パラメタファイルを変更する必要がある場合もあります。次節参照
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... での平均の方向が視野の方向となります[*]
いくつもの science file に対して ATT を使うと各ファイルごとに 姿勢がゆらいでいるために、わずかに視野方向がかわることがありますので KEY または USER の方がいいと思います
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... に書かれている値を視野方向とします[*]
このキーワードの値は frfread or frfreadng のパラメタファイル中の ``ranom''と``decnom''に指定しておけばよい
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... edit[*]
edit するには、直接パラメタファイルをエディッタに読み込むか、 もしくは、pset (便利なFTOOLSの章参照)というコマンドを使う
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... file[*]
前章参照
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
...``defATTpath''は姿勢ファイルのあるディレクトリを指定します [*]
ascalin はこのディレクトリをきちんと指定しておくと、読み込んだ science file に 対応する姿勢ファイルを読みにいきます
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... XSELECTは、かんたんなコマンドの入力によってFTOOLS[*]
FITSファイルを扱うutilityプログラム。 基本的にFITSファイルを入力し、 何らかの処理をしてFITSファイルに出力する。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... に[*]
後にascaarfなどで、スペクトルファイルに書き込まれているイメージを用いてレスポンスを作る場合、イメージは検出器座標でないと、まずい?
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ファイル名を仮に``saisho\_gis2.xco''[*]
.xcoを必ずつけて下さい
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ここでは、GISの標準的な観測モード、``PH''モード[*]
MPC モードもサポートされていますがここでは特に この扱いについては述べません
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 違うときだけ正しいディレクトリを指定します[*]
相対パスでも絶対パスでもよい
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... RTI[*]
Rise Time Invariant
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 作られていますので可能であるならばとおしましょう[*]
運用初期のデータには、かけられないデータもある。 また、これで削られるイベントは数% です
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ``extract''[*]
``bin''でも同じことが実行できますが、``extract''の方が正式なコマンドのようです。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... binsizeがあまりに小さい[*]
0.5などでもよい
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ライトカーブは見ることができます[*]
正しくないものが表示される時がある
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... と打ちます。''exposure''というパラメターを0にすると(デフォルト)、とにかく1秒でもbinにデータが入っていれば、そのbinは、以後のデータ処理で使われます[*]
0から1の間の値を指定することもできます。0.5であれば、50%以上データが存在すれば、そのbinは含まれることになります。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... を実行する[*]
quit でも q でも同じ
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
...とタイムリージョンを設定できる状態になります [*]
PGPLOT Window上でのポインターが太陽型から十字にかわる
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ``CasA\_g2\_cursor\_time001.xsl''というファイル[*]
最初のCasA_g2はsession名
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ASCAタイム[*]
1993/01/01を00:00:00として計る時間
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 読まれます[*]
show statusを実行すると``*** FILTERS ***''のところに タイムリージョンのファイルが読まれていることがわかります
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... とすると、GTI[*]
Good Time Interval
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... time\_region.dat[*]
FITS形式でも ASCIIでも構いません
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 前に``mkfilter2''をつかって衛星の環境パラメタを収めたファイルをつくりました[*]
または、ftpを用いて、取ってきたファイル(3.1参照)
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 中でないデータがほしい時はつぎのようにします[*]
このコマンドを実行したあとshow statusをすると``*** FILTERS ***''のところに 選ばれた条件がかかれる
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... XSELECTはこのときft*.mkfというファイルを探して読み込みます[*]
複数読み込むことが出来る
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
...[*][*]
show statusの``*** SELECTION ***''に状況が現れている
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... SISの*と書かれた部分は、S0、S1、又は4CCDモードか1CCDモードかなどにより、値が異なっています。詳しくは、データを取ってきたdirectoryの''aux/ad********\_***\_joblog.html.''というファイルに、イベントファイルを処理した時のログが残っているので、これで確認してみて下さい[*]
かなり大きなテキストファイルです。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... S[0,1]\_SATF[0,1,2,3]は[*]
[0,1,2,3]はチップの番号を表します
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 最終的に我々はスペクトルを得たいわけ[*]
そうでない人もいるらしい.....
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... その際イメージ上でどのような領域を選ぶかが重要になります[*]
sourceからのスペクトルなのかバックグラウンドなのかということ
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... それは自動的に使われます[*]
clearしてしまった場合は別である
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... また、filterを設定してある場合はそれも使われます[*]
show statusの``*** FILTERS ***''に設定されている条件全部
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ``e''[*]
excludeの意味
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... この二つを組み合わせればリング状にリージョンがとれます[*]
saoimage.regの中身はリング状であったら次のようになる
# image.xsl
# Sun Jun 27 16:39:06 1993
# shape x, y, [x dimension, y dimension], [angle]
CIRCLE(122,129,30.92)
-CIRCLE(122,129,19.91)
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... でのイメージができます[*]
dmprmf というFTOOLSのコマンドで、レスポンスファイル から、channelとエネルギーの対応する値を知ることができます。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 後で、xselectを抜けた後、かけることもできるので[*]
XSPEC の章でくわしく述べます
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ただし、ライトカーブ、イメージ、スペクトルともに同じフィルターの条件[*]
show filterで見ることができます
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... を読み込むには次のようにします[*]
あらかじめ INSTRUMENT の設定をしておいたほうがよい。 していない場合や違っている場合は聞かれるので適当に答える
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... モードのデータは別々に扱わなければなりません[*]
フォーマットを変換して一緒に扱う方法も存在する
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... なお、ここでは``Fast''モードのデータの扱いは述べません[*]
一言いっておくと、fasttimeとかいう FTOOLS をとおす必要があります。 これによって、時刻付けがされます。ほとんど bright モードと同じように扱えます が、当然ながらイメージをみることはできません
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 以下のようにして``bright''モードのfileのみを選びます[*]
あらかじめ set datamode BRIGHT としておいてから make obscat を 実行すると bright モードのデータのみが選ばれる
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... これは、ある一定の閾値[*]
Heat factor という
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... なって表示[*]
320x320 の大きさイメージ
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... それを最小単位のピクセル[*]
1280x1280 の大きさのイメージ
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... RBNPHA[*]
XSPEC の章でくわしく述べます
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... CCDID.eq.3''[*]
次のように書いても同じ。``CCDID==2 $\vert\vert$ CCDID==3''
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ``Bright''[*]
PHAのチャンネルが2048
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... または``Bright2''[*]
PHAのチャンネルが4096。 これはテレメトリのBrightモードとはフォーマットが違います。 よって、FaintモードからXSELECTでBright2モードに変換したものと もともとのBrightモードのデータを一緒に扱うことはできません
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ※のところで``n''とすると今までの過程で一時的に生成されたファイルは消去されます[*]
もちろん、saveしたファイルは消去されません。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... また、打ちあげ直後のデータ[*]
1993年5月28日のHydra Aまで、ただし、sequence numberで40007000、60038000、20000000のデータは除く。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... BGDファイルの中身の詳細については、BGDファイルの存在するdirectory内の、READMEファイルを読んで下さい[*]
Spread Discriminatorがoffの場合も、、、
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... Guide''に紹介されています[*]
または、''fhelp ascascreen''
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ``BRIGHT2LINEAR''[*]
小文字でもよい
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... XSPECでスペクトル解析を行なった場合、χ^2統計が用いられます。この統計を用いる場合、モデルとデータを比較する時の1つ1つのbinの中のイベントの数は、20以上必要となります[*]
詳しくは、統計の本を参照して下さい。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... スペクトルファイルに書き込むだけなので[*]
新しいcolumnを作り、1と-1のフラグを立てています。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... をたしますよ、という式をいれます[*]
この表現の仕方からわかるように、足し算ばかりでなく、スペクトルの四則演算ができる。 ただし、割算のときは、integer では問題があるので、% mathpha units=R とした方がよい
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... (3)では、CALC[*]
NULL とすると 1sec となり、任意に数字(40000.0など、4e5はだめ)を書くとその時間となります
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... (行)を指定[*]
コメントがいらない場合は 0 とする
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ができます[*]
mathphaを使ってたしても、正しくたされないということが昔ありました(たまに、1+0=0などとなってました)。正しく足されているかどうかを、確認した方がいいです。''fstatistic sum_spec.pi COUNTS -''で、totalカウントを調べることができます。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... で、バックグラウンドの引き算の時に使われます[*]
例えばGISで20x20pixelの領域を切った時、20x20/256x256=0.0061に近い値が書き込まれるようです。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 引かれてしまいます[*]
著者はこのような苦情を何度か聞きました
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... のレスポンスをつくり出すツールは``ascaarf''もしくは``jbldarf''です[*]
sisrmg, ascaarf は FTOOLS であり、純粋な表街道です
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... のレスポンスは各人がつくる必要はなく、すでに出来ているそれをつかいます[*]
2.2節で、とってきたファイル。宇宙研では、/usr/local/asca/caldb/data/asca/gis/cpf/95mar06/gis[2,3]v4_0.rmf
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... sis\_bldrsp2を使いたい方は、初めての表街道ver5.20を参照して下さい[*]
sis_bldrsp2は、RDDの扱い等、responseの経年変化に対応し切れていない面があるため、SISチームとしては、一般ユーザーにはsisrmgを使って欲しいとのことでした。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... (第1.2版)」という日本語マニュアルがついているので、そちらを参照して下さい。ただし、jbldarfでarfを作る時、xrt\_ea\_V2.0.fits、xrt\_psf\_V2.0.fitsファイルを使用した場合は、arffilterをかける必要があります[*]
ascaarfの場合、少なくともftools4.2以降では、defaultでarffilterがかかっているようです。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... (2)では、レスポンスを指定します[*]
SIS は、さきほどのsisrmgで作ったファイル、GISはgis[2,3]v4_0.rmf。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 存在するイメージ[*]
積分領域等の情報を得るため。 実際、スペクトルファイルの primary extension はWMAP(weighted map)のイメージになっており、 saoimage、ximage等で見ることができます
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 単位で指定するようです[*]
間違っていたらゴメンなさい
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... とします。[*]
addarfは以前は表街道ではなかったですが、今はHEAsoft(またはFTOOLS)に組み込まれているようです。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... を一つにすることができます[*]
一つにしたものは XSPEC で resp コマンドでよめます
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 足す必要があります[*]
s0c1.rmfとs1c3.rmf のゲインは当然あっていなければならない。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... リストファイルをつくってから、``addrmf''[*]
FTOOLS です
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... レスポンスのファイル名と重みを書いておきます[*]
重みの factor は自分で normalize しておく必要があります
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 先ず、コマンドラインから``xspec''[*]
小文字で入力
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... とします[*]
''cpd''は、``setplot device''でもOKです。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... setp)を打たなければ、認識してもらえません[*]
以前は、''set''という省略ができたのですが、最近は''setp''まで打たないとだめになりました。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 次に、レスポンス[*]
レスポンスについては``レスポンスを作る''の節参照
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... XSELECTによって生成されるGISのスペクトルファイルは通常1024channels[*]
タイミングに bit をとった場合などは 256channels
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... と打ちます[*]
端を**で指定すると、hung upすることがあります。こういう場合は、XSPEC$>$ ignore 1-0.7 とか、10.0-256 のように、端をchannelで指定すると、hung upを避けることができるようです。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... とやれば結果が出てきます[*]
local minimumになっていないか確認するために、何回かfitを行なってみる、又はerrorコマンドでパラメーターをふってみることをおすすめします。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... のようにします[*]
残差を見たい場合には、XSPEC$>$ plot ldata residual と打ちます。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... で、Δχ^2が2.3、4.61、9.21[*]
自由度が2の場合の、68%、90%、99% confidence level に相当しているようです。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... で終了します[*]
``q''だけでは、終了しなくなりました。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... などとします[*]
xspec ver9.00以前では、''$ls''と打ちます。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... script.xcm[*]
.xcm という拡張子をつけてください
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 以下のようにすると、マシにみえます[*]
くわしくは qdp のマニュアルを御覧ください
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... コマンドラインから``ximage''[*]
小文字で入力
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... プロットデバイスは次のようにして設定、変更をします[*]
PostScript にしたい場合は cpd /ps とする
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... であった[*]
XSELECT のデフォルトでは SIS の binning は 4 になっています。GIS は 1 です
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... いいでしょう[*]
Mushotzky談
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... contour[*]
smoothをしてからの方がよい
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... いくつか便利なコマンドを挙げておきます。各自研究して下さい[*]
ここに挙げた 以外にもたくさんある
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... coordinatesを描きます[*]
イメージのFITSファイルのヘッダに正しく書かれていなければなりません
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 年分点に変更します[*]
デフォルトでは 2000 年分点です
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 実行します[*]
ただし、pbmplusなどのソフトがインストールされている必要あり。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... これでFITS形式でファイルになります。あとでSAOIMAGE等で見ることができます[*]
枠などはセーブされません
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... IRAF[*]
Image Reduction and Analysis Facility
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... NOAO[*]
National Optical Astronomy Observatories
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... のディレクトリで``cl[*]
立ち上げる前に login.cl というファイルの``set home''と``set imdir''をエディット してください
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... PROS[*]
The Post-Reduction Off-line Software
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... というファイルができました[*]
実際には sis0_image.iraf.imh と sis0_image.iraf.pix の二つができる
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... まず、最初にもどります[*]
別に戻る必要はなく da$>$ plot としてもよい
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... ``eparm[*]
Edit PARaMeter
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... も表示されます[*]
イメージファイルに座標がきちんと書かれている場合に限る
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... を終了するには``logout''とします[*]
Ctrl-C などでは終了することはできない
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... IDL[*]
Interactive Data Language
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... RSI[*]
Research Systems, Inc
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... の能力をすべてこのマニュアルに載せるのは不可能ですので[*]
筆者はすべてをしりません
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... の方法を説明します[*]
宇宙研では IDL は 2000/07/26 現在 、使える状態には なっていないようです。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... が望ましい[*]
FFT を用いているから
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... のイメージ[*]
この Deconvolution では現在のところ PSF の位置依存性およびエネルギー依存性は取り込めません。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 使用例を挙げて説明いたします[*]
このマニュアル中ですべての FTOOLS の紹介はしてありません
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.