上記の例では、モデルフィットした時のパラメタはすべて共有されています。 データごとにパラメタを独立させたいときは以下のようにして読み込みます。
XSPEC> data 1:1 spec_1.pi Net count rate (cts/cm^2/s) for file 1 45.29 +/- 5.0210E-02 XSPEC> data 2:2 spec_2.pi Net count rate (cts/cm^2/s) for file 2 4.529 +/- 5.0210E-03 XSPEC> resp 1:1 resp_1.rmf XSPEC> resp 2:2 resp_2.rmf XSPEC> arf 1:1 arf_1.arf XSPEC> arf 2:2 arf_2.arfこうしてから、モデルを決めると、
XSPEC> mo bre wab
:
---------------------------------------------------------------------------
mo = wabs[2] (bremss[1])
Model Fit Model Component Parameter Value Data
par par comp group
1 1 1 bremss kT(keV) 7.00000 +/- 0.00000E+00 1
2 2 1 bremss norm 1.00000 +/- 0.00000E+00 1
3 3 2 wabs nH 10^22 1.00000 +/- 0.00000E+00 1
4 1 3 bremss kT(keV) 7.00000 = par 1 2
5 2 3 bremss norm 1.00000 = par 2 2
6 3 4 wabs nH 10^22 1.00000 = par 3 2
---------------------------------------------------------------------------
デフォルトでは、data 1 と data 2 ではパラメタがリンクされています。
data 2 の bremss の normarization をリンクからはずして独立させたい
時は、
XSPEC> new 5 0.1
:
XSPEC> show
:
---------------------------------------------------------------------------
mo = wabs[2] (bremss[1])
Model Fit Model Component Parameter Value Data
par par comp group
1 1 1 bremss kT(keV) 7.00000 +/- 0.00000E+00 1
2 2 1 bremss norm 1.00000 +/- 0.00000E+00 1
3 3 2 wabs nH 10^22 1.00000 +/- 0.00000E+00 1
4 1 3 bremss kT(keV) 7.00000 = par 1 2
5 4 3 bremss norm 0.100000 +/- 0.00000E+00 2
6 3 4 wabs nH 10^22 1.00000 = par 3 2
---------------------------------------------------------------------------
のようにします。
あとは、fit とするだけです。