6 aug 11 Mg XII 8.421 A Ly-alpha using range: 8.34 to 8.50 n = 2, norm = 1.65e-3 ======================================================================== 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.17846 +/- 0.520742 3 1 windprof u0 0.736062 +/- 9.20900E-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 1 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 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.00061E-05 +/- 4.59592E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 186.56 using 188 PHA bins and 185 degrees of freedom. *SEE BELOW (NEAR CONTOUR CALC) A SLIGHTLY (.06) BETTER NEW BEST FIT MODEL* XSPEC12>iplot PLT> wdata 8.421_nonporous.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.773355 1.74401 (-0.405107,0.565546) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.660188 0.892164 (-0.0758747,0.156102) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.71283e-05 3.1912e-05 (-2.87781e-06,1.90589e-06) ANISOTROPIC POROSITY hinf free floats to zero aniso, hinf = 0.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 1.72325 +/- 0.975522 3 1 windprof u0 0.730852 +/- 8.41435E-02 4 1 windprof h 0.500000 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 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 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 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.81296E-05 +/- 4.83828E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 189.80 using 188 PHA bins and 185 degrees of freedom. XSPEC12>iplot PLT> wdata 8.421_aniso_h0.5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.913535 2.49035 (-0.650674,0.926146) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.648656 0.815604 (-0.0727576,0.0941911) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.73018e-05 3.21336e-05 (-2.9602e-06,1.87163e-06) 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.46393 +/- 1.68583 3 1 windprof u0 0.687475 +/- 8.70301E-02 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 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 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 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.04834E-05 +/- 4.60409E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 191.52 using 188 PHA bins and 185 degrees of freedom. XSPEC12>iplot PLT> wdata 8.421_aniso_h1.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.859889 3.04461 (-0.609072,1.57565) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.629351 0.77371 (-0.058474,0.0858849) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.72326e-05 3.22943e-05 (-3.25688e-06,1.80481e-06) aniso, 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 2.37411 +/- 6.39371 3 1 windprof u0 0.648946 +/- 5.85312E-02 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 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 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 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.97369E-05 +/- 4.71881E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 198.61 using 188 PHA bins and 185 degrees of freedom. XSPEC12>iplot PLT> wdata 8.421_aniso_h5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 1.21556, 0.572257 and delta stat 0.84401, 1.33738 but latest trial 1.00744 gives 0.53767 Suggest that you check this result using the steppar command. 2 0.893906 8.55544 (-0.820289,6.84125) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.592845 0.703387 (-0.056122,0.0544199) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.7421e-05 3.2092e-05 (-2.36606e-06,2.30487e-06) ISOTROPIC POROSITY hinf free: hinf floats to zero making the 2-D steppar grid first: XSPEC12>log 8.421_iso_hinf_vs_taustar_35by35.log Logging to file:8.421_iso_hinf_vs_taustar_35by35.log XSPEC12>steppar 4 0. 5. 35 2 0. 11. 35 well, new non-porous best-fit: ======================================================================== 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.22271 +/- 2.08211 3 1 windprof u0 0.747156 +/- 0.104896 4 1 windprof h 0.0 +/- 10.9613 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 1 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 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.94421E-05 +/- 4.68821E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 186.50 using 188 PHA bins and 184 degrees of freedom. XSPEC12>iplot PLT> wdata 8.421_nonporous_better.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.776881 2.00183 (-0.426081,0.798867) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.660774 0.883956 (-0.0792218,0.143961) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.71859e-05 3.18355e-05 (-2.33394e-06,2.3157e-06) OK, now iso, hinf = 0.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 1.31179 +/- 0.646960 3 1 windprof u0 0.727089 +/- 7.90840E-02 4 1 windprof h 0.500000 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 1 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 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.00133E-05 +/- 4.59707E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 186.83 using 188 PHA bins and 185 degrees of freedom. XSPEC12>iplot PLT> wdata 8.421_iso_h0.5dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.829801 2.04649 (-0.481993,0.734695) XSPEC12>error 1. 3 Parameter Confidence Range (1) ***Warning: Number of trials exceeded before convergence. Current trial values 0.656567, 0.656557 and delta statistic 0.98786, 1.67747 ***Warning: Identical values of the parameter give different values of the statistic. Please check your result for the low end of the confidence range. 3 0.656562 0.829943 (-0.0705266,0.102854) XSPEC12>steppar 3 0.61 0.69 16 C-Statistic Delta u0 C-Statistic 3 189.67 2.4138 0 0.61 189.47 2.2093 1 0.615 189.22 1.967 2 0.62 189 1.7388 3 0.625 188.78 1.5247 4 0.63 188.55 1.2985 5 0.635 188.37 1.1094 6 0.64 188.19 0.93269 7 0.645 188.02 0.76818 8 0.65 187.87 0.6156 9 0.655 187.73 0.47472 10 0.66 187.6 0.34534 11 0.665 187.48 0.2274 12 0.67 187.38 0.1215 13 0.675 187.28 0.027169 14 0.68 187.2 -0.056237 15 0.685 187.13 -0.12933 16 0.69 XSPEC12> error 1. 17 Parameter Confidence Range (1) 17 2.71416e-05 3.19241e-05 (-2.87172e-06,1.91083e-06) iso, 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.45784 +/- 0.802925 3 1 windprof u0 0.719400 +/- 7.09794E-02 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 1 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 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.00226E-05 +/- 4.59862E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 187.12 using 188 PHA bins and 185 degrees of freedom. XSPEC12>iplot PLT> wdata 8.421_iso_h1.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.881646 2.39009 (-0.576194,0.932245) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.652539 0.803325 (-0.0668609,0.0839252) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.71555e-05 3.19403e-05 (-2.86703e-06,1.91775e-06) 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 3.19080 +/- 4.00043 3 1 windprof u0 0.681713 +/- 5.21045E-02 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 1 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 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.00990E-05 +/- 4.61308E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 189.74 using 188 PHA bins and 185 degrees of freedom. XSPEC12>iplot PLT> wdata 8.421_iso_h5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 1.34782 7.51832 (-1.8442,4.3263) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.628611 0.728586 (-0.0531107,0.0468639) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.72767e-05 3.2061e-05 (-2.82308e-06,1.96118e-06) 13may12 ADDING THE H_INF = 2 CASE - ANISO ======================================================================== 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.82520 +/- 2.60371 3 1 windprof u0 0.671389 +/- 7.00926E-02 4 1 windprof h 2.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 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 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 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.97601E-05 +/- 4.70491E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 194.13 using 188 PHA bins and 185 degrees of freedom. XSPEC12>iplot PLT> wdata 8.421_aniso_h2.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.875335 3.90563 (-0.948588,2.0817) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.617655 0.728359 (-0.0536932,0.057011) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.74462e-05 3.21151e-05 (-2.31335e-06,2.3555e-06) - ISO CASE ======================================================================== 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.79491 +/- 1.23267 3 1 windprof u0 0.706854 +/- 6.18022E-02 4 1 windprof h 2.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 1 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 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.00443E-05 +/- 4.60211E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 187.73 using 188 PHA bins and 185 degrees of freedom. XSPEC12>iplot PLT> wdata 8.421_iso_h2.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.992485 3.20768 (-0.802427,1.41277) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.645617 0.772787 (-0.0612362,0.0659333) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.71884e-05 3.19739e-05 (-2.85587e-06,1.92956e-06)