3apr11 HD 93129A - fits to Mg XII with h = 0, 1, 5 beta = 1 non-porous; hinf = 0: ======================================================================== 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 1.36532 +/- 1.01139 3 1 windprof u0 0.557386 +/- 0.119959 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 8.42100 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 3200.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 4.17620E-06 +/- 1.04739E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.12000E-04 frozen ________________________________________________________________________ C-statistic = 216.86 using 260 PHA bins and 257 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.591322 2.82853 (-0.773402,1.4638) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.45414 0.744142 (-0.102009,0.187994) OK, now doing porous models (iso and aniso; hinf = 1 and 5): iso, hinf=1: OK, taustar goes up a bit and fit gets slightly better: ======================================================================== 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 1.70459 +/- 1.45437 3 1 windprof u0 0.559115 +/- 0.116104 4 1 windprof h 1.00000 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 8.42100 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 3200.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 4.05200E-06 frozen 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.12000E-04 frozen ________________________________________________________________________ C-statistic = 216.74 using 260 PHA bins and 258 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.708003 3.63923 (-0.996583,1.93464) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.45997 0.682685 (-0.0991457,0.123569) XSPEC12>iplot PLT> wdata mgxii_beta1_iso_h1.dat Next, aniso, hinf=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 1.71758 +/- 2.07107 3 1 windprof u0 0.580759 +/- 0.201708 4 1 windprof h 1.00000 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 1 frozen 10 1 windprof rosselan 0 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 8.42100 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 3200.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 4.05200E-06 frozen 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.12000E-04 frozen ________________________________________________________________________ C-statistic = 216.70 using 260 PHA bins and 258 degrees of freedom. (note: norm was frozen I realized later, but when I refit with it thawed, it didn't change the fit) XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.731823 3.23895 (-0.986344,1.52078) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.469764 0.723525 (-0.111051,0.14271) XSPEC12>iplot PLT> wdata mgxii_beta1_aniso_h1.dat Next, iso, hinf=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 5.36270 +/- 8.81454 3 1 windprof u0 0.574697 +/- 0.116602 4 1 windprof h 5.00000 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 8.42100 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 3200.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 4.16895E-06 +/- 1.04866E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.12000E-04 frozen ________________________________________________________________________ C-statistic = 216.75 using 260 PHA bins and 257 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 4.97807, 22.5794 and delta stat 0.0294321, 1.75067 but latest trial 5.94343 gives 0.0123229 Suggest that you check this result using the steppar command. 2 1.35568 13.7787 (-3.59461,8.82843) lower limit OK, upper limit more like 16.7 XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.477704 0.663766 (-0.0940699,0.0919925) XSPEC12>iplot PLT> wdata mgxii_beta1_iso_h5.dat Next, aniso, hing=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 2.45069 +/- 7.21385 3 1 windprof u0 0.561773 +/- 0.126849 4 1 windprof h 5.00000 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 1 frozen 10 1 windprof rosselan 0 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 8.42100 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 3200.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 4.05660E-06 +/- 1.05377E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.12000E-04 frozen ________________________________________________________________________ Using energies from responses. C-statistic = 219.17 using 260 PHA bins and 257 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.442952 12.7642 (-2.18433,10.1369) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.463933 0.66177 (-0.103819,0.0940178) XSPEC12>iplot PLT> wdata mgxii_beta1_aniso_h5.dat