10月21日 月
円環ではなく普通の円で領域をとった時のパイルアップを調べてみた
1500 eV以上になるとパイルアップの影響が大きくなってしまうことがわかる
スペクトル解析(ログ)
xspecを起動して
範囲を指定しただけ
fittingもやってみようとするがよくわからない
いろいろ試してみる
森下さんがphabs*diskbbでやっていたのでおなじようにやってみる
XSPEC12>model phabs*diskbb
phabs:nH>10^22
diskbb:Tin>keV
diskbb:norm>
XSPEC12>fit
========================================================================
Model phabs<1>*diskbb<2> Source No.: 1 Active/On
Model Model Component Parameter Unit Value
par comp
1 1 phabs nH 10^22 7.90007E-02 +/- 2.45314E-03
2 2 diskbb Tin keV 0.807479 +/- 6.56873E-03
3 2 diskbb norm 50.8641 +/- 1.73095
________________________________________________________________________
Fit statistic : Chi-Squared = 3045.45 using 285 PHA bins.
***Warning: Chi-square may not be valid due to bins with zero variance
in spectrum number(s): 1
Test statistic : Chi-Squared = 3045.45 using 285 PHA bins.
Reduced chi-squared = 10.7995 for 282 degrees of freedom
Null hypothesis probability = 0.000000e+00
***Warning: Chi-square may not be valid due to bins with zero variance
in spectrum number(s): 1
もっといろいろ指定したほうがいいのか？
もう一回おなじことやったはずなのに結果が違う
========================================================================
Model phabs<1>*diskbb<2> Source No.: 1 Active/On
Model Model Component Parameter Unit Value
par comp
1 1 phabs nH 10^22 1.79623E-02 +/- 1.58749E-03
2 2 diskbb Tin keV 0.882787 +/- 6.73058E-03
3 2 diskbb norm 34.1792 +/- 1.03775
________________________________________________________________________
Fit statistic : Chi-Squared = 470.74 using 284 PHA bins.
***Warning: Chi-square may not be valid due to bins with zero variance
in spectrum number(s): 1
Test statistic : Chi-Squared = 470.74 using 284 PHA bins.
Reduced chi-squared = 1.6752 for 281 degrees of freedom
Null hypothesis probability = 9.289287e-12
***Warning: Chi-square may not be valid due to bins with zero variance
in spectrum number(s): 1
わけがわからん？？？
その時の
10月25日 金
前回のログを見返してバックグラウンドファイルが読み込まれていないことに気づいた
グルーピングをやり直す
変わっているようには見えないが…
別のモデルを使ってみる
model
|
|
diskbb
|
phabs*diskbb
|
wabs*(powerlaw+diskbb)
|
|
|
|
|
|
phabs
|
nH(10^22)
|
|
7.81769E-02 +/- 2.54565E-03
|
|
wabs
|
nH(10^22)
|
|
|
0.231297 +/-
2.10344E-02
|
powerlaw
|
PhoIndex
|
|
|
3.51043
+/- 0.192467
|
|
norm
|
|
|
8.87709E-02 +/- 7.61189E-03
|
diskbb
|
Tin(keV)
|
1.41548
+/- 1.30833E-02
|
0.799616 +/-
6.61784E-03
|
0.960157 +/-
1.40238E-02
|
|
norm
|
3.70838
+/- 0.115225
|
52.4223
+/- 1.85291
|
18.9902
+/- 1.65145
|
|
|
|
|
|
Chi-Squared
|
|
15732.19 using 285 PHA bins.
|
2966.24 using 285 PHA bins.
|
2122.44 using 285 PHA bins.
|
Reduced chi-squared
|
|
55.59076 for 283 degrees of
freedom
|
10.5186 for 282 degrees of
freedom
|
7.58015 for 280 degrees of
freedom
|
初期値はなにもいれていない
χ自乗が大きすぎる
log
10月28日 月
各モデルの形がどういうものか見てみた→まとめ
初期値に先行研究の結果の値を入れてみるがほとんど変化なし
χ自乗が大きすぎるのは
低エネルギーでのイベントを除去できていなかったからだと判明
エネルギー表示にして除去
phabs*diskbb(赤)でフィッティング
========================================================================
Model phabs<1>*diskbb<2> Source No.: 1
Active/On
Model Model Component Parameter
Unit Value
par comp
1 1
phabs
nH
10^22 1.88192E-02 +/-
1.61595E-03
2 2
diskbb
Tin
keV 0.873766
+/- 6.79051E-03
3 2
diskbb
norm
35.3836 +/-
1.09560
________________________________________________________________________
Fit statistic : Chi-Squared
= 441.26 using 284
PHA bins.
***Warning: Chi-square may not be valid due to bins with zero
variance
in spectrum number(s): 1
Test statistic : Chi-Squared
= 441.26 using 284
PHA bins.
Reduced chi-squared
= 1.5703 for 281 degrees
of freedom
Null hypothesis probability = 3.110225e-09
***Warning: Chi-square may not be valid due to bins with zero
variance
in spectrum number(s): 1
phabsの初期値を0.06にした
規格化されたχ自乗が小さくなった
wabs*(powerlaw+diskbb)
========================================================================
Model wabs<1>(powerlaw<2> + diskbb<3>) Source
No.: 1 Active/On
Model Model Component Parameter
Unit Value
par comp
1 1
wabs
nH
10^22 6.81521E-02 +/-
7.46193E-03
2 2
powerlaw
PhoIndex
2.34293 +/-
0.105106
3 2
powerlaw
norm
3.26286E-02 +/- 3.64280E-03
4 3
diskbb
Tin
keV 0.865664
+/- 1.49736E-02
5 3
diskbb
norm
29.1996 +/-
1.90359
________________________________________________________________________
Fit statistic : Chi-Squared
= 329.68 using 284
PHA bins.
***Warning: Chi-square may not be valid due to bins with zero
variance
in spectrum number(s): 1
Test statistic : Chi-Squared
= 329.68 using 284
PHA bins.
Reduced chi-squared
= 1.1816 for 279 degrees
of freedom
Null hypothesis probability = 1.991149e-02
***Warning: Chi-square may not be valid due to bins with zero
variance
in spectrum number(s): 1
log
LMC X-3は標準的なブラックホールのモデルで表せることがわかった
wabsとphabsの違いがわからない