10 aug 11 epsilon Ori /home/cohen/xspec/2008/eori vinf = 1600 beta = 1 Fe XVII @ 15.014 Emma used 14.90 to 15.124 (but it looks like maybe we should cut some out from the red wing) n = 2, norm = 1.90e-3 reproduced Emma's fit (C = 237.9) betaSob = 0 tau0star free ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 0.423252 +/- 0.707950 3 1 windprof u0 0.397733 +/- 0.226657 4 1 windprof h 0.0 frozen 5 1 windprof tau0star 4.30582 +/- 21.9464 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 0 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 15.0140 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 1600.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.09591E-04 +/- 2.21852E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.90000E-03 frozen ________________________________________________________________________ C-statistic = 234.92 using 274 PHA bins and 270 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_tau0free_bsob0.dat PLT> exit XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.122503 0.959075 (-0.294876,0.541696) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.310848 0.499751 (-0.0912126,0.0976905) XSPEC12>error 1. 5 Parameter Confidence Range (1) 5 0.494407 40.1721 (-3.92178,35.7559) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000201319 0.00022443 (-7.21447e-06,1.58966e-05) Now, fixed values of tua0star: tau0star = 1 ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 0.257803 +/- 0.198026 3 1 windprof u0 0.466650 +/- 3.03890E-02 4 1 windprof h 0.0 frozen 5 1 windprof tau0star 1.00000 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 0 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 15.0140 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 1600.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.08710E-04 +/- 2.21649E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.90000E-03 frozen ________________________________________________________________________ C-statistic = 235.47 using 274 PHA bins and 271 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_tau0star1_bsob0.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.0877686 0.46618 (-0.170034,0.208377) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.433981 0.495942 (-0.0326693,0.029291) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000201134 0.000224886 (-7.57657e-06,1.61759e-05) tau0star = 5: ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 0.443497 +/- 0.328873 3 1 windprof u0 0.390998 +/- 3.74799E-02 4 1 windprof h 0.0 frozen 5 1 windprof tau0star 5.00000 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 0 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 15.0140 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 1600.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.08459E-04 +/- 2.21793E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.90000E-03 frozen ________________________________________________________________________ C-statistic = 235.00 using 274 PHA bins and 271 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_tau0star5_bsob0.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.163162 0.785233 (-0.280335,0.341735) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.35377 0.428784 (-0.0372279,0.0377864) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000200906 0.000224731 (-7.55377e-06,1.62712e-05) tau0star = 10 ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 0.539721 +/- 0.396490 3 1 windprof u0 0.362415 +/- 4.12642E-02 4 1 windprof h 0.0 frozen 5 1 windprof tau0star 10.0000 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 0 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 15.0140 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 1600.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.08136E-04 +/- 2.21870E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.90000E-03 frozen ________________________________________________________________________ C-statistic = 235.19 using 274 PHA bins and 271 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_tau0star10_bsob0.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.20502 0.955642 (-0.334702,0.41592) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.320354 0.404231 (-0.0420606,0.0418163) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000200663 0.000224584 (-7.4723e-06,1.64484e-05) Now, betaSob = 1 and repeat free tau0star and three fixed values: tau0star free ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 0.211090 +/- 0.281809 3 1 windprof u0 0.464637 +/- 0.119436 4 1 windprof h 0.0 frozen 5 1 windprof tau0star 6.78576 +/- 22.1611 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 1.00000 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 0 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 15.0140 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 1600.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.09098E-04 +/- 2.21920E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.90000E-03 frozen ________________________________________________________________________ C-statistic = 236.90 using 274 PHA bins and 270 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_tau0free_bsob1.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Identical values of the parameter give different values of the statistic. Please check your result for the high end of the confidence range. 2 0.0500702 0.479024 (-0.158082,0.270872) XSPEC12>steppar 2 0.4 0.8 20 C-Statistic Delta taustar C-Statistic 2 237.44 0.3395 0 0.4 237.59 0.48944 1 0.42 237.76 0.65759 2 0.44 237.86 0.75922 3 0.46 237.97 0.86073 4 0.48 238.07 0.96583 5 0.5 238.18 1.0713 6 0.52 238.28 1.179 7 0.54 238.39 1.2873 8 0.56 238.5 1.3968 9 0.58 238.61 1.5068 10 0.6 238.72 1.6177 11 0.62 238.83 1.7287 12 0.64 238.94 1.8397 13 0.66 239.06 1.9505 14 0.68 239.17 2.0611 15 0.7 239.28 2.1712 16 0.72 239.39 2.2808 17 0.74 239.49 2.3898 18 0.76 239.6 2.4983 19 0.78 239.71 2.6059 20 0.8 XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.384885 0.556196 (-0.0798174,0.0914943) XSPEC12>error 1. 5 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 6.77239, 0.372393 and delta stat 0, 1.02235 but latest trial 0.442722 gives 1.64191 Suggest that you check this result using the steppar command. 5 3.57239 61.0657 (-2.89312,54.6002) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000201282 0.000223819 (-6.68896e-06,1.58483e-05) tau0star = 1: ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 0.105180 +/- 0.109813 3 1 windprof u0 0.527447 +/- 3.26775E-02 4 1 windprof h 0.0 frozen 5 1 windprof tau0star 1.00000 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 1.00000 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 0 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 15.0140 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 1600.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.08467E-04 +/- 2.21661E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.90000E-03 frozen ________________________________________________________________________ C-statistic = 237.68 using 274 PHA bins and 271 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_tau0star1_bsob1.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.0138194 0.225577 (-0.0913603,0.120398) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.492566 0.558984 (-0.034881,0.0315378) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.00020094 0.000224745 (-7.5273e-06,1.62783e-05) tau0star = 5 ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 0.190171 +/- 0.173137 3 1 windprof u0 0.475533 +/- 4.14129E-02 4 1 windprof h 0.0 frozen 5 1 windprof tau0star 5.00000 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 1.00000 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 0 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 15.0140 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 1600.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.08108E-04 +/- 2.21782E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.90000E-03 frozen ________________________________________________________________________ C-statistic = 236.97 using 274 PHA bins and 271 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_tau0star5_bsob1.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.0511165 0.371992 (-0.139054,0.181822) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.43512 0.517127 (-0.0404127,0.0415941) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000200622 0.000224539 (-7.48564e-06,1.64313e-05) tau0star = 10 ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 0.237407 +/- 0.209120 3 1 windprof u0 0.451526 +/- 4.58447E-02 4 1 windprof h 0.0 frozen 5 1 windprof tau0star 10.0000 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 1.00000 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 0 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 15.0140 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 1600.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.07732E-04 +/- 2.21839E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.90000E-03 frozen ________________________________________________________________________ C-statistic = 237.06 using 274 PHA bins and 271 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_tau0star10_bsob1.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.0712312 0.457918 (-0.166176,0.220511) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.40372 0.498475 (-0.0478057,0.0469489) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000200329 0.000224361 (-7.40244e-06,1.66295e-05)