Mg XI He-like complex using hewind 10Oct08 atomdb says: 9.1687 9.2282 9.2312 9.3143 and I'm going to use (for now) Maurice's phiratio value of 131 fitting the continuum 8.80:9.00; 9.50:9.60 ======================================================================== 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 1.54537E-03 +/- 1.50595E-04 ________________________________________________________________________ C-statistic = 312.11 using 352 PHA bins and 351 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>plot XSPEC12>goodness 1000 nosim 5.10% of realizations are < best fit statistic 312.11 XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) 2 0.001460 0.001634 (-0.000085,0.000089) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 0.001407 0.001692 (-0.000138,0.000147) let's fit the complex itself over a relatively wide range (at first) 9.07:9.40 ======================================================================== Model hewind<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hewind q 0.0 frozen 2 1 hewind taustar 0.893159 +/- 0.164803 3 1 hewind u0 0.641809 +/- 2.27599E-02 4 1 hewind h 0.0 frozen 5 1 hewind tau0star 0.0 frozen 6 1 hewind beta 1.00000 frozen 7 1 hewind betaSob 0.0 frozen 8 1 hewind G 0.700000 frozen 9 1 hewind numerica 0 frozen 10 1 hewind anisotro 0 frozen 11 1 hewind rosselan 0 frozen 12 1 hewind expansio 0 frozen 13 1 hewind thick 0 frozen 14 1 hewind Z 12 frozen 15 1 hewind shift "mA" 0.0 frozen 16 1 hewind velocity "km/s" 2250.00 frozen 17 1 hewind phiratio 131.000 frozen 18 1 hewind verbose 0 frozen 19 1 hewind norm 1.76525E-04 +/- 9.95263E-06 20 2 powerlaw PhoIndex 2.00000 frozen 21 2 powerlaw norm 1.55000E-03 frozen ________________________________________________________________________ C-statistic = 410.96 using 392 PHA bins and 389 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>plot XSPEC12>goodness 1000 nosim 41.00% of realizations are < best fit statistic 410.96 whoops; G was fixed at 0.7 for the fit above (which is the canonical value) Now, redoing the fit with G free: ======================================================================== Model hewind<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hewind q 0.0 frozen 2 1 hewind taustar 0.917765 +/- 0.173245 3 1 hewind u0 0.644817 +/- 2.41303E-02 4 1 hewind h 0.0 frozen 5 1 hewind tau0star 0.0 frozen 6 1 hewind beta 1.00000 frozen 7 1 hewind betaSob 0.0 frozen 8 1 hewind G 0.723272 +/- 5.07934E-02 9 1 hewind numerica 0 frozen 10 1 hewind anisotro 0 frozen 11 1 hewind rosselan 0 frozen 12 1 hewind expansio 0 frozen 13 1 hewind thick 0 frozen 14 1 hewind Z 12 frozen 15 1 hewind shift "mA" 0.0 frozen 16 1 hewind velocity "km/s" 2250.00 frozen 17 1 hewind phiratio 131.000 frozen 18 1 hewind verbose 0 frozen 19 1 hewind norm 1.78239E-04 +/- 9.91005E-06 20 2 powerlaw PhoIndex 2.00000 frozen 21 2 powerlaw norm 1.55000E-03 frozen ________________________________________________________________________ C-statistic = 410.31 using 392 PHA bins and 388 degrees of freedom. XSPEC12>goodness 1000 nosim 37.00% of realizations are < best fit statistic 410.31 XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 0.755779 1.111577 (-0.161988,0.193810) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 0.656486 1.238630 (-0.261281,0.320863) XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) 3 0.620795 0.670015 (-0.024022,0.025198) XSPEC12>error 2.7 3 Parameter Confidence Range (2.700000) 3 0.606456 0.687123 (-0.038361,0.042306) XSPEC12>error 1. 8 Parameter Confidence Range (1.000000) 8 0.674543 0.779233 (-0.048730,0.055960) XSPEC12>error 2.7 8 Parameter Confidence Range (2.700000) 8 0.643282 0.814351 (-0.079991,0.091079) XSPEC12>error 1. 19 Parameter Confidence Range (1.000000) infinite loop, so XSPEC12>steppar 19 1.5e-4 2.1e-4 30 C-Statistic Delta norm C-Statistic 19 452.76900 42.45830 0 0.00015 447.15019 36.83948 1 0.00015 441.96427 31.65357 2 0.00015 437.20188 26.89117 3 0.00016 432.85391 22.54321 4 0.00016 428.91155 18.60085 5 0.00016 425.36628 15.05558 6 0.00016 422.20984 11.89913 7 0.00016 419.43421 9.12350 8 0.00017 417.03164 6.72094 9 0.00017 414.99462 4.68391 10 0.00017 413.31585 3.00514 11 0.00017 411.98825 1.67755 12 0.00017 411.00498 0.69428 13 0.00018 410.35936 0.04866 14 0.00018 410.04493 -0.26577 15 0.00018 410.05540 -0.25530 16 0.00018 410.38468 0.07398 17 0.00018 411.02682 0.71612 18 0.00019 411.97607 1.66537 19 0.00019 413.22681 2.91611 20 0.00019 414.77360 4.46290 21 0.00019 416.61113 6.30043 22 0.00019 418.73423 8.42353 23 0.00020 421.13789 10.82719 24 0.00020 423.81721 13.50651 25 0.00020 426.76742 16.45672 26 0.00020 429.98389 19.67319 27 0.00020 433.46209 23.15139 28 0.00021 437.19762 26.88692 29 0.00021 441.18619 30.87548 30 0.00021 hmmmm, better fit at norm=1.81e-4; but when I refit, none of the other parameters changed by more than 0.5%; so stick with the fit above XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/mgxi_917_hewind_smooth.dat confidence limit plot (taustar vs. uo): XSPEC12>log /home/cohen/xspec/2007/zpup/mgxi_917_hewind_smooth_35by35.log Logging to file:/home/cohen/xspec/2007/zpup/mgxi_917_hewind_smooth_35by35.log XSPEC12>steppar 2 0.01 4.5 35 3 0.48 0.99 35 Note that we're letting G be a free parameter too