20131203

12月3日　火 (log)

quantities について調べた

現状報告では特になし
引き続きmos2とPNの解析を進める
来週は勉強報告

栗山くんが作ったプログラムの使い方を教えてもらった

解析ログは増えたが何も進んでいない

12月6日　金 (log)

フィルターのとりかたをちょっと変えた
目算で範囲を指定するのをやめて
栗山くんのプログラムで出た範囲のオーバーしてる部分のみをカットした

mos1もやり直した

とった領域は同じ

mos1	mos2	PN

modelはphabs*(diskbb+powerlaw)でフィット

mos1

Model phabs<1>(diskbb<2> + powerlaw<3>) Source No.: 1   Active/On
Model Model Component  Parameter  Unit     Value
 par  comp
   1    1   phabs      nH         10^22    0.562913     +/-  5.30181E-02  
   2    2   diskbb     Tin        keV      1.88280      +/-  0.214919     
   3    2   diskbb     norm                1.77400E-03  +/-  8.73352E-04  
   4    3   powerlaw   PhoIndex            4.85733      +/-  0.345296     
   5    3   powerlaw   norm                9.36987E-04  +/-  1.45859E-04

log

Reduced chi-squared = 1.2524 for 154 degrees of freedom

***Warning: Zero alpha-matrix diagonal element for parameter 2
***Warning: Zero alpha-matrix diagonal element for parameter 3
という表示が多い
前の解析でもこの表示は多かったのでフィルターがけの範囲を変えたこととは関係ない

mos2

Model phabs<1>(diskbb<2> + powerlaw<3>) Source No.: 1   Active/On
Model Model Component  Parameter  Unit     Value
 par  comp
   1    1   phabs      nH         10^22    0.425996     +/-  5.58246E-02  
   2    2   diskbb     Tin        keV      0.165010     +/-  1.72365E-02  
   3    2   diskbb     norm                311.547      +/-  289.119      
   4    3   powerlaw   PhoIndex            2.38215      +/-  0.111224     
   5    3   powerlaw   norm                2.78715E-04  +/-  3.82348E-05

log

Reduced chi-squared = 1.0155 for 152 degrees of freedom

ノルムが異常な値

Due to zero model norms, the following fit parameters are temporarily frozen:2
が多いがこれは問題ない？

Model phabs<1>(diskbb<2> + powerlaw<3>) Source No.: 1   Active/On
Model Model Component  Parameter  Unit     Value
 par  comp
   1    1   phabs      nH         10^22    0.620772     +/-  3.30207E-02  
   2    2   diskbb     Tin        keV      1.57396      +/-  7.80637E-02  
   3    2   diskbb     norm                3.90907E-03  +/-  8.98451E-04  
   4    3   powerlaw   PhoIndex            5.32243      +/-  0.225627     
   5    3   powerlaw   norm                1.08393E-03  +/-  1.05657E-04

log

Reduced chi-squared = 1.2294 for 448 degrees of freedom

***Warning: Zero alpha-matrix diagonal element for parameter 2
***Warning: Zero alpha-matrix diagonal element for parameter 3
が多い

内縁半径も出したい
NGC6946　X-1は
赤方偏移　z=0.000133
0.6 Mpc = 600 kpc 　　←これが間違いだったのか(1月10日)
降着円盤の傾き：わからないのでとりあえず50°で←ちなみに澤口さんは60でやってる(1 月10日)

	mos1	mos2	PN
内縁半径 [km]	3.15	1320.9	4.68
光度 [J/s]	1.61e+30	1.67e+31	1.74e+30

どうみても異常

今日はここまで
~~データの統合はeモザイクを使う…のか？~~

12月10日　火 (log)

やっぱりmos2がおかしいんだと思う

Tinも小さい(10分の1くらい？)

Tinが10分の1
Kが10^5倍
になっているのでflux自体はそんなにおかしくないはず
解析作業でなにかミスがある?

mos1の結果をmos2の初期値として入れてみる

mos2

Model phabs<1>(diskbb<2> + powerlaw<3>) Source No.: 1   Active/On
Model Model Component  Parameter  Unit     Value
 par  comp
   1    1   phabs      nH         10^22    0.676856     +/-  6.28518E-02  
   2    2   diskbb     Tin        keV      1.49498      +/-  0.141169     
   3    2   diskbb     norm                5.09014E-03  +/-  2.13351E-03  
   4    3   powerlaw   PhoIndex            5.59535      +/-  0.406463     
   5    3   powerlaw   norm                1.25184E-03  +/-  2.28199E-04

Reduced chi-squared = 1.1821 for 152 degrees of freedom

まだノルムが大きいかと思ったけど

~~fornax:~/bin% perl radius 5.09014E-03 1.49498 50 600~~
~~R=5.33927592977933 [km]~~
~~L=1.83962822608644e+30 [J/s] =1.83962822608644e+37 [erg/s]~~

かなり良くなってる

	mos1	mos2	PN
内縁半径 [km]	3.15	5.34	4.68
光度 [J/s]	1.61e+30	1.84e+30	1.74e+30