27jan11 Maurice was skeptical about the low significance of the S XV flux; I don't seem to have records of that fit... so: continuum fit on 4.6:4.9; 5.2:5.4 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.14305E-04 +/- 6.55333E-05 ________________________________________________________________________ C-statistic = 140.68 using 196 PHA bins and 195 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.000182694 0.000249541 (-3.16116e-05,3.52357e-05) ======================================================================== 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.318027 +/- 1.29489 3 1 hewind u0 0.859127 +/- 0.235416 4 1 hewind h 0.0 frozen 5 1 hewind tau0star 0.0 frozen 6 1 hewind beta 0.700000 frozen 7 1 hewind betaSob 0.0 frozen 8 1 hewind G 0.700000 frozen 9 1 hewind numerica 1 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 16 frozen 15 1 hewind shift mA 0.0 frozen 16 1 hewind velocity km/s 3200.00 frozen 17 1 hewind phiratio 4.00000 frozen 18 1 hewind verbose 0 frozen 19 1 hewind norm 4.11096E-06 +/- 2.28010E-06 20 2 powerlaw PhoIndex 2.00000 frozen 21 2 powerlaw norm 2.14300E-04 frozen ________________________________________________________________________ C-statistic = 105.76 using 118 PHA bins and 115 degrees of freedom. XSPEC12>error 1. 19 Parameter Confidence Range (1) 19 2.99573e-06 5.24677e-06 (-1.08967e-06,1.16137e-06) XSPEC12>error 2.7 19 Parameter Confidence Range (2.7) 19 2.38755e-06 6.08879e-06 (-1.69785e-06,2.00338e-06) XSPEC12>error 4. 19 Parameter Confidence Range (4) 19 2.07477e-06 6.58202e-06 (-2.01063e-06,2.49662e-06) XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 Apparent non-monotonicity in statistic space detected. Current bracket values 1.61497, 1.74129 and delta stat 0.907805, 1.08337 but latest trial 1.67704 gives 0.90774 Suggest that you check this result using the steppar command. 2 0 1.67813 (-0.318027,1.3601) XSPEC12>error 1. 3 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0.99 3 0.623916 0 (-0.235211,-0.859127) XSPEC12>iplot PLT> wdata SXV_hewind_b0.7_v3200_megonly.dat Next, same model, but allowing G to be a free parameter Pretty significant change: ======================================================================== 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.537676 +/- 2.98403 3 1 hewind u0 0.844984 +/- 0.386083 4 1 hewind h 0.0 frozen 5 1 hewind tau0star 0.0 frozen 6 1 hewind beta 0.700000 frozen 7 1 hewind betaSob 0.0 frozen 8 1 hewind G 1.84671 +/- 2.05576 9 1 hewind numerica 1 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 16 frozen 15 1 hewind shift mA 0.0 frozen 16 1 hewind velocity km/s 3200.00 frozen 17 1 hewind phiratio 4.00000 frozen 18 1 hewind verbose 0 frozen 19 1 hewind norm 4.26244E-06 +/- 2.28320E-06 20 2 powerlaw PhoIndex 2.00000 frozen 21 2 powerlaw norm 2.14300E-04 frozen ________________________________________________________________________ C-statistic = 104.14 using 118 PHA bins and 114 degrees of freedom. XSPEC12>error 1. 19 Parameter Confidence Range (1) 19 3.19664e-06 5.46518e-06 (-1.07012e-06,1.19842e-06) XSPEC12>error 4. 19 Parameter Confidence Range (4) 19 2.24912e-06 6.81879e-06 (-2.01061e-06,2.55906e-06) XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 Apparent non-monotonicity in statistic space detected. Current bracket values 0.535612, 3.51964 and delta stat 0, 1.39797 but latest trial 3.05941 gives 1.47972 Suggest that you check this result using the steppar command. 2 0 2.02763 (-0.535612,1.49201) XSPEC12>error 1. 3 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0.99 3 0.589501 0 (-0.258698,-0.848199) XSPEC12>error 1. 8 Parameter Confidence Range (1) 8 0.853377 5.07051 (-0.93802,3.27911) XSPEC12>iplot PLT> wdata SXV_hewind_b0.7_v3200_Gfree_megonly.dat Next, we can test Gaussian profile fits: hegauss ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 20.0000 +/- 329.272 2 1 hegauss G 1.39215 +/- 1.08081 3 1 hegauss sigma_v km/s 493.399 +/- 358.712 4 1 hegauss delta_v km/s -692.787 +/- 239.592 5 1 hegauss Z 16.0000 frozen 6 1 hegauss norm 4.14964E-06 +/- 2.20014E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 2.14000E-04 frozen ________________________________________________________________________ C-statistic = 99.28 using 118 PHA bins and 113 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>error 1. 1 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 20 1 2.86162 0 (-17.1384,-20) So... low f/i is rule out (or at least not indicated) XSPEC12>error 2.7 1 Parameter Confidence Range (2.7) ***Warning: Parameter pegged at hard limit: 20 1 1.41883 0 (-18.5812,-20) XSPEC12>error 4. 1 Parameter Confidence Range (4) ***Warning: Parameter pegged at hard limit: 20 1 1.08396 0 (-18.916,-20) XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.831252 2.74963 (-0.561588,1.35679) Parameter Confidence Range (1) 3 191.646 772.663 (-301.619,279.398) XSPEC12>error 1. 4 Parameter Confidence Range (1) 4 -942.16 -452.263 (-249.042,240.855) XSPEC12>error 1. 6 Parameter Confidence Range (1) 6 3.12475e-06 5.33247e-06 (-1.02486e-06,1.18287e-06) XSPEC12>error 2.7 6 Parameter Confidence Range (2.7) 6 2.52745e-06 6.15076e-06 (-1.62215e-06,2.00116e-06) XSPEC12>error 4. 6 Parameter Confidence Range (4) 6 2.22024e-06 6.65385e-06 (-1.92937e-06,2.50425e-06) XSPEC12>iplot PLT> wdata SXV_hegauss_v1.dat Testing 95% confidence extreme model with R = f/i = 1.08: ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 1.08000 frozen 2 1 hegauss G 1.89950 +/- 1.36891 3 1 hegauss sigma_v km/s 472.441 +/- 344.704 4 1 hegauss delta_v km/s -688.047 +/- 369.807 5 1 hegauss Z 16.0000 frozen 6 1 hegauss norm 4.02383E-06 +/- 2.19842E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 2.14000E-04 frozen ________________________________________________________________________ C-statistic = 103.02 using 118 PHA bins and 114 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.977462 4.56408 (-0.940462,2.64616) XSPEC12>error 4. 2 Parameter Confidence Range (4) ***Warning: Parameter pegged at hard limit: 10 2 0.494429 0 (-1.40659,-1.90102) XSPEC12>error 1. 3 Parameter Confidence Range (1) ***Warning: Zero alpha-matrix diagonal element for parameter 4 ... Suggest that you check this result using the steppar command. 3 68.5486 1046.09 (-393.344,584.194) XSPEC12>error 1. 4 Parameter Confidence Range (1) 4 -1068.87 -237.265 (-381.846,449.758) XSPEC12>error 1. 6 Parameter Confidence Range (1) 6 3.03456e-06 5.22087e-06 (-9.84581e-07,1.20172e-06) XSPEC12>error 4. 6 Parameter Confidence Range (4) 6 2.12959e-06 6.53309e-06 (-1.89462e-06,2.50888e-06) XSPEC12>iplot PLT> wdata SXV_hegauss_v2.dat