on d139: O VIII 19.969 A (MEG only) continuum on 18.70:18.82; 19.10:19.25 ======================================================================== Model powerlaw<1> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 powerlaw PhoIndex 2.00000 frozen 2 1 powerlaw norm 2.64153E-03 +/- 9.11961E-04 ________________________________________________________________________ C-statistic = 66.56 using 104 PHA bins and 103 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>plot XSPEC12>goodness 100 nosim 0.00% of realizations are < best fit statistic 66.56 *note: when I repeated the Ne X 12.134 fit on d139, I got MC 0.0% too, for the identical fit that on ruby had something like 89%... is there a bug in the new xspec? XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.00210841 0.00326068 (-0.000533119,0.00061915) XSPEC12>error 2.7 2 Parameter Confidence Range (2.7) 2 0.00180368 0.00370509 (-0.000837849,0.00106356) fitting the line on 18.81:19.12 ======================================================================== 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 8.60295E-15 +/- 0.170284 3 1 windprof u0 0.572291 +/- 2.66219E-02 4 1 windprof h 0.0 frozen 5 1 windprof tau0star 0.0 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 18.9690 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2350.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 8.64390E-04 +/- 9.14571E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.64000E-03 frozen ________________________________________________________________________ C-statistic = 145.74 using 122 PHA bins and 119 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0 0.0217129 (-4.17968e-15,0.0217129) XSPEC12>error 2.7 2 Parameter Confidence Range (2.7) ***Warning: Parameter pegged at hard limit: 0 2 0 0.0577702 (-1.9682e-15,0.0577702) XSPEC12>error 4. 2 Parameter Confidence Range (4) ***Warning: Parameter pegged at hard limit: 0 2 0 0.0843764 (-8.62538e-16,0.0843764) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.538385 0.600214 (-0.0341555,0.0276738) XSPEC12>error 2.7 3 Parameter Confidence Range (2.7) 3 0.516985 0.618255 (-0.0556078,0.0456623) XSPEC12>error 4. 3 Parameter Confidence Range (4) 3 0.504804 0.627912 (-0.0678385,0.0552699) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000818194 0.000917949 (-4.61906e-05,5.35637e-05) XSPEC12>error 2.7 17 Parameter Confidence Range (2.7) 17 0.000787755 0.000951386 (-7.66292e-05,8.70017e-05) XSPEC12>error 4. 17 Parameter Confidence Range (4) 17 0.000771199 0.000970386 (-9.31851e-05,0.000106003) XSPEC12>iplot PLT> wdata /home/cohen/xspec/2009/iori/oviii_18969.dat XSPEC12>log /home/cohen/xspec/2009/iori/oviii_18969.log XSPEC12>steppar 2 0.01 4.5 35 3 0.48 0.99 35 Now, redoing, focused in on the relevant region of parameter space: XSPEC12>log /home/cohen/xspec/2009/iori/oviii_18969_20by20.log Logging to file:/home/cohen/xspec/2009/iori/oviii_18969_20by20.log XSPEC12>steppar 2 0.01 0.40 20 3 0.33 0.72 20 Seeing how much difference it makes if we treat the two components of the line separately XSPEC12>newpar windprof + windprof + pow ======================================================================== Model windprof<1> + windprof<2> + powerlaw<3> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 8.39224E-15 +/- 0.164317 3 1 windprof u0 0.586783 +/- 2.57554E-02 4 1 windprof h 0.0 frozen 5 1 windprof tau0star 0.0 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 18.9671 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2350.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 5.79500E-04 +/- 6.22298E-05 18 2 windprof q 0.0 frozen 19 2 windprof taustar 8.39224E-15 = windprof[1]:taustar 20 2 windprof u0 0.586783 = windprof[1]:u0 21 2 windprof h 0.0 frozen 22 2 windprof tau0star 0.0 frozen 23 2 windprof beta 1.00000 frozen 24 2 windprof betaSob 0.0 frozen 25 2 windprof numerica 0 frozen 26 2 windprof anisotro 0 frozen 27 2 windprof rosselan 0 frozen 28 2 windprof expansio 0 frozen 29 2 windprof thick 0 frozen 30 2 windprof waveleng A 18.9725 frozen 31 2 windprof shift mA 0.0 frozen 32 2 windprof velocity km/s 2350.00 frozen 33 2 windprof verbose 0 frozen 34 2 windprof norm 2.89750E-04 = 0.5*windprof[1]:norm 35 3 powerlaw PhoIndex 2.00000 frozen 36 3 powerlaw norm 2.64000E-03 frozen ________________________________________________________________________ C-statistic = 146.03 using 122 PHA bins and 119 degrees of freedom. Note that the fit (including uo) hardly changed at all For consistency, we'll fit it as a single line, as we usually do.