10月4日金
アメリカ政府の機能停止騒ぎのせいでNASAのサーバーがダウン
今日もHEASARCにはつながらない いつ復旧するのか
NASAで
はなくESAからデータを
落とすことにした (手順)
落としたデータは0116900501(森下さんと同じもの)
各.tarファイルを解凍
% tar xvf 0030_0116900501.TAR
SASの環境設
定をやり直した
emchain,epchainを やってみた
10月6日 日
イメージ作成
coma:~/analysis/0116900501/PPS% evselect table=P0116900501M1S007MIEVLI0000.FTZ:EVENTS
withimageset=yes imageset=m1_image.fits
xcolumn=X ycolumn=Y imagebinning=imageSize ximagesize=600
yimagesize=600
を実行すると
m1_image.fitsというファイルができる
ファイル名に「EV」と入っているやつじゃないとダメ
できたファイルをds9で開く
coma:~/analysis/0116900501/PPS% ds9
m1_image.fits
coma:~/analysis/0116900501/PPS%
ds9 m2_image.fits
先輩にならって色を「b」にして「解析→平滑化設定」からカーネルを10でスムージング
他の設定はそのまま
(これで合ってるのか不安)
ds9でいろいろいじってみた結果
ライトカーブ作成
SAS起動して
coma:~/analysis/0116900501/PPS% evselect table=P0116900501M1S007MIEVLI0000.FTZ:EVENTS
withrateset=yes rateset=m1_Lc.fits
maketimecolumn=yes timecolumn=TIME timebinsize=50
makeratecolumn=yes
centaurus:~/analysis/0116900501/PPS% dsplot
table=m1_Lc_01.fits x=TIME y=RATE.ERROR
xmgrace: Symbol `_XmStrings' has different size in shared
object, consider re-linking
実は失敗した（グラフが真っ白）
→次の日解決
10月7日 月
イメージ作成のつづき
pnのデータでもやってみる（昨日と同様）
logスケール 平滑化なし 色ｂ
試行錯誤
ライトカーブ作成のつづき
失敗の原因
mos1
mos2
pn
10月8日 火
ds9とライトカーブをいじってた
フィルターのかけ方を調べる
10月9日 水
フィルター
centaurus:~/analysis/0116900501/PPS% evselect
table=P0116900501M1S007MIEVLI0000.FTZ:EVENTS
withfilteredset=yes expression='(PATTERN
<= 12)&&(PI in [200:12000])&&#XMMEA_EM'
filteredset=m1_filtered_200-12000.fits
filtertype=expression keepfilteroutput=yes updateexposure=yes
filterexposure=yes
とりあえずmos1だけお試し
線形
0.2〜12keVの範囲
(10/17 追記)イメージ作成してねえよ！だから変なんだよ！
image
gtiファイルを作成
hercules:~/analysis/0116900501/PPS% tabgtigen table=m1_Lc_100.fits:RATE gtiset=m1_gti_time.fits timecolumn=TIME
expression='(TIME in [66354826.9 :
66360826.9])'
gtiファイルでフィルター
hercules:~/analysis/0116900501/PPS% evselect table=m1_filtered_200-12000.fits:EVENTS
withfilteredset=yes expression='GTI(m1_gti_time.fits,TIME)'
filteredset=m1_filt_time.fits filtertype=expression
keepfilteroutput=yes updateexposure=yes filterexposure=yes
ライトカーブ,イメージを作って一応確認
(追記ここまで)
ライトカーブから不要な部分を取り除こうとしてみるがうまくいかない…
ログ