1月14日　火(log)

現状報告
エネルギー範囲を論文に合わせること
澤口さんの論文によると
0.2keV付近は検出器由来？
荷電粒子が検出器周りの物質と相互作用することで起こるフレアーが10-12keVで起こり易いノイズがあるらしい



他の観測データも解析
パワースペクトルの作成


1月16日　木(log)


エネルギー範囲を0.3keV〜10keVにしてフィット


mos1(黒)の低エネルギー(0.3keVギリギリ？)にすごいズレがある？

パラメーター
========================================================================
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.401934     +/-  2.05925E-02  
   2    2   diskbb     Tin        keV      0.167425     +/-  7.25632E-03  
   3    2   diskbb     norm                258.336      +/-  96.4468      
   4    3   powerlaw   PhoIndex            2.42732      +/-  4.33613E-02  
   5    3   powerlaw   norm                2.85801E-04  +/-  1.55150E-05  
                           Data group: 2
   6    1   phabs      nH         10^22    0.401934     = 1
   7    2   diskbb     Tin        keV      0.167425     = 2
   8    2   diskbb     norm                258.336      = 3
   9    3   powerlaw   PhoIndex            2.42732      = 4
  10    3   powerlaw   norm                2.85801E-04  = 5
                           Data group: 3
  11    1   phabs      nH         10^22    0.401934     = 1
  12    2   diskbb     Tin        keV      0.167425     = 2
  13    2   diskbb     norm                258.336      = 3
  14    3   powerlaw   PhoIndex            2.42732      = 4
  15    3   powerlaw   norm                2.85801E-04  = 5
________________________________________________________________________


Fit statistic : Chi-Squared =         852.29 using 817 PHA bins.

Test statistic : Chi-Squared =         852.29 using 817 PHA bins.
 Reduced chi-squared =         1.0496 for    812 degrees of freedom 
 Null hypothesis probability =   1.585940e-01

内縁半径：Rin=13643.6 km
光度：L=3.43e+39 [erg/s]
質量：M=1515.6 太陽質量


前のエネルギー範囲との比較
0.2-12.0keV
 0.3-10.0keV



低エネルギーの一点はやはりノイズ？





1月20日　月(log)

ノイズらしきものを取り除いてみる



========================================================================
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.399536     +/-  2.05147E-02 
   2    2   diskbb     Tin        keV      0.167953     +/-  7.30716E-03 
   3    2   diskbb     norm                250.153      +/-  93.4246     
   4    3   powerlaw   PhoIndex            2.42557      +/-  4.34084E-02 
   5    3   powerlaw   norm                2.85031E-04  +/-  1.54974E-05 
                           Data group: 2
   6    1   phabs      nH         10^22    0.399536     = 1
   7    2   diskbb     Tin        keV      0.167953     = 2
   8    2   diskbb     norm                250.153      = 3
   9    3   powerlaw   PhoIndex            2.42557      = 4
  10    3   powerlaw   norm                2.85031E-04  = 5
                           Data group: 3
  11    1   phabs      nH         10^22    0.399536     = 1
  12    2   diskbb     Tin        keV      0.167953     = 2
  13    2   diskbb     norm                250.153      = 3
  14    3   powerlaw   PhoIndex            2.42557      = 4
  15    3   powerlaw   norm                2.85031E-04  = 5
________________________________________________________________________


Fit statistic : Chi-Squared =         844.54 using 816 PHA bins.

Test statistic : Chi-Squared =         844.54 using 816 PHA bins.
 Reduced chi-squared =         1.0414 for    811 degrees of freedom
 Null hypothesis probability =   2.010058e-01


まあほんのちょっとは良くなったかな


内縁半径R=13420.5 [km]
光度　　L=3.43354e+39 [erg/s]
質量　　M=1491.17 [M]


パワースペクトル関係のログ


1月21日　火(log)
他の観測データについてもエネルギー範囲を変更した
0500730101
 各パラメータ

===========================================
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.381419     +/-  4.45857E-02 
   2    2   diskbb     Tin        keV      0.165346     +/-  1.44165E-02 
   3    2   diskbb     norm                239.744      +/-  182.463     
   4    3   powerlaw   PhoIndex            2.18378      +/-  9.41099E-02 
   5    3   powerlaw   norm                1.87948E-04  +/-  2.23799E-05 
                           Data group: 2
   6    1   phabs      nH         10^22    0.381419     = 1
   7    2   diskbb     Tin        keV      0.165346     = 2
   8    2   diskbb     norm                239.744      = 3
   9    3   powerlaw   PhoIndex            2.18378      = 4
  10    3   powerlaw   norm                1.87948E-04  = 5
                           Data group: 3
  11    1   phabs      nH         10^22    0.381419     = 1
  12    2   diskbb     Tin        keV      0.165346     = 2
  13    2   diskbb     norm                239.744      = 3
  14    3   powerlaw   PhoIndex            2.18378      = 4
  15    3   powerlaw   norm                1.87948E-04  = 5
________________________________________________________________________

Fit statistic : Chi-Squared =         291.27 using 272 PHA bins.

Test statistic : Chi-Squared =         291.27 using 272 PHA bins.
 Reduced chi-squared =         1.0909 for    267 degrees of freedom
 Null hypothesis probability =   1.471586e-01
Reduced chi-squared　
1.1147　→　1.0909　(ちょっと改善された)

内縁半径
 R=13138.3 [km]
光度
 L=3.00138 e+39 [erg/s]
質量
 M=1459.81 [M太陽]


0500730201
 各パラメータ

===========================================
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.456700     +/-  5.25525E-02 
   2    2   diskbb     Tin        keV      0.138630     +/-  1.14051E-02 
   3    2   diskbb     norm                1019.53      +/-  873.414     
   4    3   powerlaw   PhoIndex            2.51910      +/-  8.43995E-02 
   5    3   powerlaw   norm                2.81937E-04  +/-  2.84741E-05 
                           Data group: 2
   6    1   phabs      nH         10^22    0.456700     = 1
   7    2   diskbb     Tin        keV      0.138630     = 2
   8    2   diskbb     norm                1019.53      = 3
   9    3   powerlaw   PhoIndex            2.51910      = 4
  10    3   powerlaw   norm                2.81937E-04  = 5
                           Data group: 3
  11    1   phabs      nH         10^22    0.456700     = 1
  12    2   diskbb     Tin        keV      0.138630     = 2
  13    2   diskbb     norm                1019.53      = 3
  14    3   powerlaw   PhoIndex            2.51910      = 4
  15    3   powerlaw   norm                2.81937E-04  = 5
________________________________________________________________________


Fit statistic : Chi-Squared =         277.54 using 316 PHA bins.

Test statistic : Chi-Squared =         277.54 using 316 PHA bins.
 Reduced chi-squared =        0.89242 for    311 degrees of freedom
 Null hypothesis probability =   9.140103e-01
Reduced chi-squared　 0.87485　→　0.89242　(ちょっと改善された)

内縁半径
 R=27093.6 [km]
光度
 L=2.91894 e+39 [erg/s]
質量
 M=3010.40 [M太陽]


まとめ
対象天体
LMC X-3

NGC6946 X-1
観測データ
0116900501
0691570101
 0500730101
 0500730201
内縁半径 Rin[km]
39.9
13420.5 13138.3 27093.6
光度 L[erg/s]
1.28 e+38
3.43354 e+39
 3.00138 e+39
 2.91894  e+39
質量 M[太陽質量]
4.43
1491.17 1459.81 3010.40
Reduced chi-squared
1.2061
1.0414 1.0909 0.89242