big plot on 11.9:12.4 Looks like a definite blend near 12.2; Here are the ranges we will try to use: continuum: 11.94-12.02; 12.31-12.34 line: 12.03-12.18 ======================================================================== 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.87704E-03 +/- 4.01420E-04 ________________________________________________________________________ C-statistic = 109.92 using 124 PHA bins and 123 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>plot XSPEC12>goodness 100 nosim 83.00% of realizations are < best fit statistic 109.92 XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 0.001649 0.002125 (-0.000228,0.000248) Now, eliminating the longer wavelength section: ======================================================================== 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.38571E-03 +/- 4.11441E-04 ________________________________________________________________________ C-statistic = 67.33 using 92 PHA bins and 91 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>goodness 100 nosim 44.00% of realizations are < best fit statistic 67.33 XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 0.001165 0.001633 (-0.000220,0.000248) Well, there's a tiny bit of overlap (almost) in the 68% confidence regions; so since it's not clear whether it's legit to include the longer wavelength region, let's go with the value that's consistent with both: 1.64e-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 2.02984 +/- 0.306142 3 1 windprof u0 0.680813 +/- 6.48780E-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" 12.1339 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.69344E-04 +/- 1.74818E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.64000E-03 frozen ________________________________________________________________________ C-statistic = 191.35 using 176 PHA bins and 173 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>plot XSPEC12>goodness 1000 nosim 63.10% of realizations are < best fit statistic 191.35 XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 1.752831 2.267835 (-0.276946,0.238058) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 1.608393 2.450939 (-0.421423,0.421123) XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) 3 0.631357 0.730591 (-0.049443,0.049792) XSPEC12>error 2.7 3 Parameter Confidence Range (2.700000) 3 0.608263 0.783731 (-0.072541,0.102927) XSPEC12>error 1. 17 Parameter Confidence Range (1.000000) 17 0.000262 0.000280 (-0.000008,0.000011) XSPEC12>error 2.7 17 Parameter Confidence Range (2.700000) 17 0.000256 0.000286 (-0.000013,0.000016) --- *As an extra test, we'll try fitting the non-porous model with q free: ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.190810 +/- 0.156177 2 1 windprof taustar 2.19353 +/- 0.317019 3 1 windprof u0 0.649658 +/- 6.47279E-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" 12.1339 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.64701E-04 +/- 1.81979E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.64000E-03 frozen ________________________________________________________________________ C-statistic = 189.47 using 176 PHA bins and 172 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>plot XSPEC12>goodness 1000 nosim 64.90% of realizations are < best fit statistic 189.47 Modest improvement in goodness of fit, but not highly significant (Delta C ~ 2) And q not that different from zero; testing the confidence limits (on all parameters of interest; to see if a free q strongly affects the constraints on other parameters): XSPEC12>error 1. 1 Parameter Confidence Range (1.0000000) 1 0.054447 0.356785 (-0.136337,0.166000) XSPEC12>error 2.7 1 Parameter Confidence Range (2.700000) 1 -0.034891 0.464240 (-0.225683,0.273449) XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) 2 1.890062 2.498826 (-0.303504,0.305259) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 1.720252 2.719722 (-0.473313,0.526157) XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) 3 0.601458 0.703180 (-0.048219,0.053503) XSPEC12>error 2.7 3 Parameter Confidence Range (2.700000) 3 0.574924 0.757611 (-0.074753,0.107934) XSPEC12>error 1. 17 Parameter Confidence Range (1.000000) 17 0.000257 0.000276 (-0.000008,0.000011) XSPEC12>error 2.7 17 Parameter Confidence Range (2.700000) 17 0.000251 0.000282 (-0.000013,0.000017) taustar vs. uo log file... --- compare: taustar vs. uo log file for q=0 frozen Source File: acis640_pha2.fits{9} Net count rate (cts/s) for Spectrum:1 5.751e-03 +/- 3.762e-04 Assigned to Data Group 1 and Plot Group 1 Noticed Channels: 5958-5986 Telescope: CHANDRA Instrument: ACIS Channel Type: PI Exposure Time: 67636.0078 sec Using Response (RMF) File meg_m1.rmf for Source 1 Using Auxiliary Response (ARF) File acis640MEG_-1_garf.fits Source spectrum counts: 389.000000 Model predicted rate: 6.06914E-03 Source File: acis640_pha2.fits{10} Net count rate (cts/s) for Spectrum:2 6.151e-03 +/- 3.860e-04 Assigned to Data Group 1 and Plot Group 2 Noticed Channels: 5958-5986 Telescope: CHANDRA Instrument: ACIS Channel Type: PI Exposure Time: 67636.0078 sec Using Response (RMF) File meg_p1.rmf for Source 1 Using Auxiliary Response (ARF) File acis640MEG_1_garf.fits Source spectrum counts: 416.000000 Model predicted rate: 6.06123E-03 Source File: acis640_pha2.fits{3} Net count rate (cts/s) for Spectrum:3 2.528e-03 +/- 3.267e-04 Assigned to Data Group 1 and Plot Group 3 Noticed Channels: 3722-3780 Telescope: CHANDRA Instrument: ACIS Channel Type: PI Exposure Time: 67636.0078 sec Using Response (RMF) File heg_m1.rmf for Source 1 Using Auxiliary Response (ARF) File acis640HEG_-1_garf.fits Source spectrum counts: 171.000000 Model predicted rate: 2.20590E-03 Source File: acis640_pha2.fits{4} Net count rate (cts/s) for Spectrum:4 1.479e-03 +/- 2.863e-04 Assigned to Data Group 1 and Plot Group 4 Noticed Channels: 3722-3780 Telescope: CHANDRA Instrument: ACIS Channel Type: PI Exposure Time: 67636.0078 sec Using Response (RMF) File heg_p1.rmf for Source 1 Using Auxiliary Response (ARF) File acis640HEG_1_garf.fits Source spectrum counts: 100.000000 Model predicted rate: 1.49232E-03 --- Now repeating the fitting but excluding bins > 12.16 (cutting out more of the wing) Here's the best-fit model from above, just superimposed (compared to, statistically; no fitting yet) on the newly truncated data: ======================================================================== 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.02981 +/- 0.306146 3 1 windprof u0 0.680803 +/- 6.48844E-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" 12.1339 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.69344E-04 +/- 1.74817E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.64000E-03 frozen ________________________________________________________________________ Using energies from responses. C-statistic = 169.42 using 152 PHA bins and 149 degrees of freedom. Now, fitting: ======================================================================== 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.96048 +/- 0.315589 3 1 windprof u0 0.660180 +/- 6.29936E-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" 12.1339 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.73813E-04 +/- 1.85884E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.64000E-03 frozen ________________________________________________________________________ C-statistic = 169.01 using 152 PHA bins and 149 degrees of freedom. Only a very slight change XSPEC12>goodness 1000 nosim 72.70% of realizations are < best fit statistic 169.01 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/nex_12134_smooth_1216.dat XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 1.689714 2.210592 (-0.273589,0.247288) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 1.530822 2.396726 (-0.440724,0.425180) XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) 3 0.614504 0.711653 (-0.048259,0.048891) XSPEC12>error 2.7 3 Parameter Confidence Range (2.700000) 3 0.589617 0.762448 (-0.073151,0.099680) XSPEC12>error 1. 17 Parameter Confidence Range (1.000000) 17 0.000264 0.000284 (-0.000008,0.000012) XSPEC12>error 2.7 17 Parameter Confidence Range (2.700000) 17 0.000258 0.000290 (-0.000014,0.000018)