12月10日 火 (log)
現状報告
先週のフィルターがけについて
ブラックホールの時間変動までカットしてしまう可能性が
対象領域のライトカーブからバックグラウンドのライトカーブを引くのが良い
光度について
ULXにしては光度が低い？
XSPEC
のluminosity,fluxコマンドで出す
luminosityはJじゃなくergで書こう
とりあえず前のmos1のデータでluminosityとfluxを出してみるモデルはph*(di+po)
XSPEC12>flux 0.2 10
Model Flux 0.00029708 photons (6.8298e-13 ergs/cm^2/s) range
(0.20000 - 10.000 keV)
XSPEC12>lumi 0.2 10.0 0.000133
Model Luminosity 2.6460e+37 ergs/s
(0.20000 - 10.000 keV rest frame)
(z = 0.000133 H0 = 70.0 q0 =
0.00 Lambda0 = 0.730)
計算は間違ってなかったがやっぱり低い
先行研究では10^40近くある
12月13日 金 (log)
X線の時間変動について調べました→.pdf
来週ちゃんと発表します
フィルターかける前のライトカーブをとってみる
中心領域(〜12″)
|
バックグラウンド領域(12~17″)
|
全領域
|
|
|
|
フィルターかける必要なかった？
mos1,mos2,PNを一緒にフィット
GTIフィルターかける前のデータを使ってみる
log
phabs*(diskbb+powerlaw)
========================================================================
Model phabs<1>(diskbb<2> + powerlaw<3>) Source No.: 1 Active/On
Model Model Component Parameter Unit Value
par comp
Data group: 1
1 1 phabs nH 10^22 0.395808 +/- 2.04845E-02
2 2 diskbb Tin keV 0.168427 +/- 7.57009E-03
3 2 diskbb norm 228.531 +/- 87.2970
4 3 powerlaw PhoIndex 2.41768 +/- 4.33297E-02
5 3 powerlaw norm 2.79330E-04 +/- 1.52021E-05
Data group: 2
6 1 phabs nH 10^22 0.395808 = 1
7 2 diskbb Tin keV 0.168427 = 2
8 2 diskbb norm 228.531 = 3
9 3 powerlaw PhoIndex 2.41768 = 4
10 3 powerlaw norm 2.79330E-04 = 5
Data group: 3
11 1 phabs nH 10^22 0.395808 = 1
12 2 diskbb Tin keV 0.168427 = 2
13 2 diskbb norm 228.531 = 3
14 3 powerlaw PhoIndex 2.41768 = 4
15 3 powerlaw norm 2.79330E-04 = 5
________________________________________________________________________
Fit statistic : Chi-Squared = 864.41 using 842 PHA bins.
Test statistic : Chi-Squared = 864.41 using 842 PHA bins.
Reduced chi-squared = 1.0327 for 837 degrees of freedom
Null hypothesis probability = 2.485570e-01
内縁半径は
centaurus:~/bin% perl radius
228.531 0.168427 50 600
R=1131.33197481263
[km]
L=1.33060998754365e+31 [J/s] =1.33060998754365e+38
[erg/s]
XSPEC12>flux 0.2 10.0
Spectrum Number: 1
Data Group Number: 1
Model Flux 0.00031849 photons (7.8068e-13 ergs/cm^2/s) range
(0.20000 - 10.000 keV)
Spectrum Number: 2
Data Group Number: 2
Model Flux 0.00031849 photons (7.8068e-13 ergs/cm^2/s) range
(0.20000 - 10.000 keV)
Spectrum Number: 3
Data Group Number: 3
Model Flux 0.00031822 photons (7.8042e-13 ergs/cm^2/s) range
(0.20000 - 10.000 keV)
XSPEC12>lumi 0.2 10.0 0.000133
Spectrum Number: 1
Data Group Number: 1
Model Luminosity 3.0245e+37
ergs/s (0.20000 - 10.000 keV rest frame)
Spectrum Number: 2
Data Group Number: 2
Model Luminosity 3.0245e+37
ergs/s (0.20000 - 10.000 keV rest frame)
Spectrum Number: 3
Data Group Number: 3
Model Luminosity 3.0235e+37
ergs/s (0.20000 - 10.000 keV rest frame)
(z = 0.000133 H0 = 70.0 q0 =
0.00 Lambda0 = 0.730)
ULX(>10^39)にしてはまだ小さい
phabsを抜くと
XSPEC12>flux 0.2 12
Spectrum Number: 1
Data Group Number: 1
Model Flux 0.00038245 photons (8.3181e-13 ergs/cm^2/s) range
(0.20000 - 12.000 keV)
Spectrum Number: 2
Data Group Number: 2
Model Flux 0.00038245 photons (8.3181e-13 ergs/cm^2/s) range
(0.20000 - 12.000 keV)
Spectrum Number: 3
Data Group Number: 3
Model Flux 0.00038245 photons (8.3181e-13 ergs/cm^2/s) range
(0.20000 - 12.000 keV)
XSPEC12>lumi 0.2 12.0 0.000133
Spectrum Number: 1
Data Group Number: 1
Model Luminosity 3.2225e+37
ergs/s (0.20000 - 12.000 keV rest frame)
Spectrum Number: 2
Data Group Number: 2
Model Luminosity 3.2225e+37
ergs/s (0.20000 - 12.000 keV rest frame)
Spectrum Number: 3
Data Group Number: 3
Model Luminosity 3.2225e+37
ergs/s (0.20000 - 12.000 keV rest frame)
(z = 0.000133 H0 = 70.0 q0 =
0.00 Lambda0 = 0.730)
12月17日 火 (log)
LMC X-3でも同じように解析しなおしてみる
mos1とmos2のみ統合できた
========================================================================
Model phabs<1>(diskbb<2> + powerlaw<3>) Source No.: 1 Active/On
Model Model Component Parameter Unit Value
par comp
Data group: 1
1 1 phabs nH 10^22 6.73269E-02 +/- 5.40744E-03
2 2 diskbb Tin keV 0.835958 +/- 1.08804E-02
3 2 diskbb norm 31.8850 +/- 1.69147
4 3 powerlaw PhoIndex 2.20841 +/- 6.06434E-02
5 3 powerlaw norm 3.64316E-02 +/- 3.05694E-03
Data group: 2
6 1 phabs nH 10^22 6.73269E-02 = 1
7 2 diskbb Tin keV 0.835958 = 2
8 2 diskbb norm 31.8850 = 3
9 3 powerlaw PhoIndex 2.20841 = 4
10 3 powerlaw norm 3.64316E-02 = 5
Reduced chi-squared = 1.2061 for 555 degrees of freedom
内縁半径は
R=45.1673081427641 [km]
L=1.28709336933463e+31 [J/s] =1.28709336933463e+38 [erg/s]
こっちは大き杉ないか
(追記 2014/1/6)ログ
レポートのために見直し
PNの解析ができなかった原因はコマンドミスかと思ったが
それを直してもスペクトルがおかしい