12月17日　火 (log)
現状報告

光度を出すときはphabsを抜いた場合のも載せる
ライトカーブはビンニングしないと時間変動がわからない

他の観測データも解析しよう
とりあえずXMM-Newtonで
澤口さんの論文より
Obs.ID　　　　優先順位(PNの観測時間の長い順)
0093641501　　10
0093641601　　PNは無し
0093641701　　PNは無し
0200670101　　8
0200670201　　9
0200670301　　7
0200670401　　5
0401360101　　6
0401360201　　3
0401360301　　3
0500730201　　1
0500730101　　2

とりあえず0500730101をダウンロード
0500730201長谷しゅんのとこからコピー


12月19日　木 (log)

emchain,epchainがなんか成功した


mos1
 mos2
 PN
0500730101




こいつだけsmoothingした
0500730201







12月20日　 (log)
X線変動の発表をした


12月23日　 (log)

領域指定をミスっていた
無駄になった1時間くらいのログ


先に0500730101(2番目に長い方)

範囲はこれまでと同じ(12″〜17″)


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.639814     +/-  5.19921E-02  
   2    2   diskbb     Tin        keV      1.74769      +/-  0.155909     
   3    2   diskbb     norm                2.40202E-03  +/-  9.04937E-04  
   4    3   powerlaw   PhoIndex            5.63016      +/-  0.349985     
   5    3   powerlaw   norm                9.57727E-04  +/-  1.46863E-04  
                           Data group: 2
   6    1   phabs      nH         10^22    0.639814     = 1
   7    2   diskbb     Tin        keV      1.74769      = 2
   8    2   diskbb     norm                2.40202E-03  = 3
   9    3   powerlaw   PhoIndex            5.63016      = 4
  10    3   powerlaw   norm                9.57727E-04  = 5
                           Data group: 3
  11    1   phabs      nH         10^22    0.639814     = 1
  12    2   diskbb     Tin        keV      1.74769      = 2
  13    2   diskbb     norm                2.40202E-03  = 3
  14    3   powerlaw   PhoIndex            5.63016      = 4
  15    3   powerlaw   norm                9.57727E-04  = 5
________________________________________________________________________


Fit statistic : Chi-Squared =         338.05 using 277 PHA bins.

Test statistic : Chi-Squared =         338.05 using 277 PHA bins.
 Reduced chi-squared =         1.2428 for    272 degrees of freedom 
 Null hypothesis probability =   3.925321e-03

ノルムが小さい
PhoIndexがでかい

光度は
 
XSPEC12>lumi 0.2 12.0 0.000133
 Spectrum Number: 1
 Data Group Number: 1
Model Luminosity 2.6119e+37 ergs/s  (0.20000 - 12.000 keV rest frame)
 Spectrum Number: 2
 Data Group Number: 2
Model Luminosity 2.6119e+37 ergs/s  (0.20000 - 12.000 keV rest frame)
 Spectrum Number: 3
 Data Group Number: 3
Model Luminosity 2.6116e+37 ergs/s  (0.20000 - 12.000 keV rest frame)
     (z = 0.000133 H0 =  70.0 q0 =  0.00 Lambda0 = 0.730)

光度
内縁半径は
 fornax:~/bin% perl radius 2.40202E-03 1.74769 50 600
 R=3.66780139480806 [km]
小さい



パワースペクトルの作り方→　す ざくファーストステップ

試しにやってみる(log)
12秒の円領域のライトカーブをパワースペクトルにした

設定は適当に決めた


12月24日　(log)


次に0500730201(長い方)

範囲は



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.551814     +/-  4.69209E-02 
   2    2   diskbb     Tin        keV      1.60678      +/-  0.132070    
   3    2   diskbb     norm                3.04449E-03  +/-  1.13561E-03 
   4    3   powerlaw   PhoIndex            5.01733      +/-  0.327690    
   5    3   powerlaw   norm                7.50436E-04  +/-  1.03136E-04 
________________________________________________________________________

光度は

XSPEC12>lumi 0.2 12.0 0.000133
 Spectrum Number: 1
 Data Group Number: 1
Model Luminosity 2.5717e+37 ergs/s  (0.20000 - 12.000 keV rest frame)
 Spectrum Number: 2
 Data Group Number: 2
Model Luminosity 2.5717e+37 ergs/s  (0.20000 - 12.000 keV rest frame)
 Spectrum Number: 3
 Data Group Number: 3
Model Luminosity 2.5711e+37 ergs/s  (0.20000 - 12.000 keV rest frame)
     (z = 0.000133 H0 =  70.0 q0 =  0.00 Lambda0 = 0.730)