Ne X 9.7082 (Ly gamma) 9.5:9.6 and 9.78:9.93 ======================================================================== 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.16305E-03 +/- 2.02995E-04 ________________________________________________________________________ C-statistic = 253.92 using 292 PHA bins and 291 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>plot XSPEC12>goodness 1000 nosim 11.30% of realizations are < best fit statistic 253.92 XSPEC12>plot XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 0.002042 0.002289 (-0.000121,0.000126) [9.62:9.80] lower 68% conf on norm = 7.4e-6 ======================================================================== 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.05465 +/- 1.16191 3 1 windprof u0 0.990000 +/- 1.69921 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" 9.70820 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 9.39735E-06 +/- 3.82490E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.16000E-03 frozen ________________________________________________________________________ Using energies from responses. C-statistic = 197.14 using 212 PHA bins and 209 degrees of freedom. XSPEC12> goodness 1000 nosim 12.90% of realizations are < best fit statistic 197.14 ======================================================================== 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 0.614312 +/- 0.745567 3 1 windprof u0 0.677314 +/- 0.127609 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" 9.70820 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 9.44968E-06 +/- 3.77774E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.16000E-03 frozen ________________________________________________________________________ C-statistic = 195.79 using 212 PHA bins and 209 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. ***Warning: New best fit found, fit parameters will be set to new values. 2 0.102120 1.672433 (-0.481542,1.088771) XSPEC12>goodness 1000 nosim 10.30% of realizations are < best fit statistic 195.81 PLT> wdata /home/cohen/xspec/2007/zpup/nex_97082_smooth.dat XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 0.095300 1.689813 (-0.517561,1.076952) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) ***Warning: Parameter pegged at hard limit: 0 2 0.095300 2.973798 (-0.485395,2.393102) XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) 3 0.572299 0.777004 (-0.104660,0.100045) XSPEC12>error 2.7 3 Parameter Confidence Range (2.700000) ***Warning: Parameter pegged at hard limit: 0.99 3 0.478108 0.777004 (-0.192728,0.106169) XSPEC12>steppar 2 0.01 0.15 30 C-Statistic Delta taustar C-Statistic 2 197.55887 1.11363 0 0.01000 197.69499 1.24975 1 0.01467 197.47043 1.02518 2 0.01933 197.41929 0.97405 3 0.02400 197.36795 0.92271 4 0.02867 197.31095 0.86571 5 0.03333 197.26614 0.82090 6 0.03800 197.22142 0.77618 7 0.04267 197.17852 0.73327 8 0.04733 197.13731 0.69207 9 0.05200 197.09768 0.65244 10 0.05667 197.05951 0.61427 11 0.06133 197.02271 0.57747 12 0.06600 196.97990 0.53466 13 0.07067 196.95060 0.50536 14 0.07533 196.91215 0.46691 15 0.08000 196.87816 0.43292 16 0.08467 196.84696 0.40172 17 0.08933 196.81671 0.37147 18 0.09400 196.78739 0.34215 19 0.09867 196.75896 0.31372 20 0.10333 196.73138 0.28613 21 0.10800 196.70461 0.25936 22 0.11267 196.67861 0.23337 23 0.11733 196.65337 0.20812 24 0.12200 196.62884 0.18359 25 0.12667 196.60499 0.15975 26 0.13133 196.58182 0.13657 27 0.13600 196.55927 0.11403 28 0.14067 196.53735 0.09211 29 0.14533 196.51602 0.07078 30 0.15000 XSPEC12>steppar 2 1.5 3.5 20 C-Statistic Delta taustar C-Statistic 2 196.57815 0.15848 0 1.50000 196.70876 0.28909 1 1.60000 196.83360 0.41394 2 1.70000 196.95974 0.54007 3 1.80000 197.08982 0.67015 4 1.90000 197.22391 0.80424 5 2.00000 197.36142 0.94176 6 2.10000 197.50418 1.08451 7 2.20000 197.41654 0.99687 8 2.30000 197.55728 1.13761 9 2.40000 197.70854 1.28888 10 2.50000 197.86697 1.44730 11 2.60000 198.03028 1.61061 12 2.70000 198.19672 1.77705 13 2.80000 198.36494 1.94527 14 2.90000 198.53386 2.11419 15 3.00000 198.70262 2.28295 16 3.10000 Fit Parameter: 3 has pegged 198.87052 2.45085 17 3.20000 Fit Parameter: 3 has pegged Fit Parameter: 3 has pegged 199.03699 2.61732 18 3.30000 Fit Parameter: 3 has pegged Fit Parameter: 3 has pegged 199.20160 2.78193 19 3.40000 Fit Parameter: 3 has pegged Fit Parameter: 3 has pegged 199.36398 2.94431 20 3.50000 XSPEC12>steppar 3 0.75 0.9 30 C-Statistic Delta u0 C-Statistic 3 196.44505 -0.00018 0 0.75000 196.52759 0.08236 1 0.75500 196.60316 0.15793 2 0.76000 196.67497 0.22974 3 0.76500 196.74306 0.29782 4 0.77000 196.80754 0.36231 5 0.77500 196.86874 0.42350 6 0.78000 196.92701 0.48177 7 0.78500 196.96128 0.51605 8 0.79000 197.01249 0.56726 9 0.79500 197.06068 0.61544 10 0.80000 197.10615 0.66092 11 0.80500 197.14922 0.70399 12 0.81000 197.19016 0.74493 13 0.81500 197.22918 0.78394 14 0.82000 197.26492 0.81969 15 0.82500 197.29690 0.85167 16 0.83000 197.32526 0.88002 17 0.83500 197.35009 0.90485 18 0.84000 197.37143 0.92619 19 0.84500 197.38927 0.94404 20 0.85000 197.40366 0.95842 21 0.85500 197.41551 0.97028 22 0.86000 197.42502 0.97978 23 0.86500 197.43201 0.98678 24 0.87000 197.43630 0.99106 25 0.87500 197.43757 0.99233 26 0.88000 197.43544 0.99021 27 0.88500 197.42988 0.98465 28 0.89000 197.42326 0.97803 29 0.89500 197.41574 0.97051 30 0.90000 XSPEC12>steppar 3 0.9 0.99 9 C-Statistic Delta u0 C-Statistic 3 197.47735 1.05769 0 0.90000 197.39650 0.97683 1 0.91000 197.36763 0.94796 2 0.92000 197.32923 0.90956 3 0.93000 197.29243 0.87276 4 0.94000 197.25465 0.83498 5 0.95000 197.20892 0.78925 6 0.96000 197.16348 0.74381 7 0.97000 197.12228 0.70261 8 0.98000 197.08742 0.66775 9 0.99000 XSPEC12>error 1. 17 Parameter Confidence Range (1.000000) 17 0.000008 0.000011 (-0.000002,0.000002) XSPEC12>error 2.7 17 Parameter Confidence Range (2.700000) 17 0.000006 0.000013 (-0.000003,0.000003) XSPEC12>log /home/cohen/xspec/2007/zpup/nex_97082_smooth_35by35.log Logging to file:/home/cohen/xspec/2007/zpup/nex_97082_smooth_35by35.log XSPEC12>steppar 2 0.01 4.5 35 3 0.48 0.99 35 OK... this looks reasonable, except the fit looks high in the wings. Let's try lowering the continuum by a factor of two, and refitting: ======================================================================== 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.91550 +/- 1.18836 3 1 windprof u0 0.440471 +/- 0.109834 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" 9.70820 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.02844E-05 +/- 4.12506E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.08000E-03 frozen ________________________________________________________________________ C-statistic = 209.45 using 212 PHA bins and 209 degrees of freedom. XSPEC12>goodness 1000 nosim 56.30% of realizations are < best fit statistic 209.45 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/nex_97082_smooth_lowcont.dat XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 0.982246 3.660245 (-0.933391,1.744607) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 0.608298 6.401649 (-1.307365,4.485986) XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) 3 0.327004 0.545863 (-0.113440,0.105419) XSPEC12>error 2.7 3 Parameter Confidence Range (2.700000) ***Warning: Parameter pegged at hard limit: 0.99 3 0.200635 0.545863 (-0.239808,0.105419) XSPEC12>error 1. 17 Parameter Confidence Range (1.000000) 17 0.000018 0.000023 (-0.000002,0.000002) XSPEC12>error 2.7 17 Parameter Confidence Range (2.700000) 17 0.000017 0.000024 (-0.000003,0.000004) XSPEC12>log /home/cohen/xspec/2007/zpup/nex_97082_smooth_lowcont_35by35.log Logging to file:/home/cohen/xspec/2007/zpup/nex_97082_smooth_lowcont_35by35.log XSPEC12>steppar 2 0.01 4.5 35 3 0.48 0.99 35