27mar11 Using James's final reduction of the full (combined GTO + GO) datasets. Total exposure time: 172,200 seconds Default loaddata_both uses "back none" XSPEC12>wdata cen1a_dataonly.dat Looks like data really only on 2.5 - 9.5 A. ======================================================================== Model bapec<1>*TBabs<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 bapec kT keV 1.84384 +/- 6.87595E-02 2 1 bapec Abundanc 1.00000 frozen 3 1 bapec Redshift 0.0 frozen 4 1 bapec Velocity km/s 1597.97 +/- 156.321 5 1 bapec norm 1.79499E-03 +/- 1.66947E-04 6 2 TBabs nH 10^22 2.25330 +/- 0.124918 ________________________________________________________________________ C-statistic = 6145.82 using 8396 PHA bins and 8392 degrees of freedom. XSPEC12>error 1. 1 Parameter Confidence Range (1) 1 1.78188 1.90782 (-0.0631192,0.0628203) XSPEC12>error 1. 4 Parameter Confidence Range (1) 4 1452.94 1812.03 (-141.009,218.079) XSPEC12>error 1. 5 Parameter Confidence Range (1) 5 0.00170397 0.00189096 (-8.93313e-05,9.76608e-05) XSPEC12>error 1. 6 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 2.13685, 2.13137 and delta stat 0.959614, 1.06352 but latest trial 2.1352 gives 0.926546 Suggest that you check this result using the steppar command. 6 2.13411 2.37436 (-0.117311,0.122938) (these look accurate, based on the steppar) XSPEC12>iplot PLT> wdata apec1T_solar.dat Now allowing the abundance to be free: ======================================================================== Model bapec<1>*TBabs<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 bapec kT keV 1.49179 +/- 7.04878E-02 2 1 bapec Abundanc 0.500365 +/- 8.74696E-02 3 1 bapec Redshift 0.0 frozen 4 1 bapec Velocity km/s 1404.67 +/- 154.935 5 1 bapec norm 3.02420E-03 +/- 5.84672E-04 6 2 TBabs nH 10^22 2.74101 +/- 0.154429 ________________________________________________________________________ C-statistic = 6111.78 using 8396 PHA bins and 8391 degrees of freedom. XSPEC12>error 1. 1 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 1.54142, 1.54143 and delta stat 0.978991, 1.01226 but latest trial 1.54143 gives 1.01556 Suggest that you check this result using the steppar command. 1 1.44835 1.54143 (-0.0434477,0.049631) fine steppar grid showed good agreement with this XSPEC12>error 1. 2 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 0.500368, 0.587838 and delta stat 0, 1.38835 but latest trial 0.583524 gives 1.79769 Suggest that you check this result using the steppar command. 2 0.438816 0.544103 (-0.0615521,0.0437348) lower limit is fine, but steppar shows upper limit is 0.58 XSPEC12>error 1. 4 Parameter Confidence Range (1) 4 1258.41 1573.36 (-146.078,168.875) XSPEC12>error 1. 5 Parameter Confidence Range (1) 5 0.0027738 0.00327722 (-0.000253264,0.000250154) Suggest that you check this result using the steppar command. 6 2.60592 2.81821 (-0.135077,0.0772147) lower limit OK, but upper is closer to 2.87 XSPEC12>iplot PLT> wdata cen1a_apec1T_nonsolar.dat OK... adding a secont temperature component: (The plots of the data and 1-T apec models certainly look like another temperature component would be warranted) I wasn't really intending to do a solar abundance model, but did do one along the way to a global best-fit 2-T model. So, here's the best-fit solar abundance 2-T bapec model (but no confidence limits): ======================================================================== Model (bapec<1> + bapec<2>)TBabs<3> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 bapec kT keV 0.469892 +/- 3.27973E-02 2 1 bapec Abundanc 1.00000 frozen 3 1 bapec Redshift 0.0 frozen 4 1 bapec Velocity km/s 1549.50 +/- 138.678 5 1 bapec norm 3.01087E-03 +/- 1.18053E-03 6 2 bapec kT keV 2.43199 +/- 0.161518 7 2 bapec Abundanc 1.00000 = 2 8 2 bapec Redshift 0.0 frozen 9 2 bapec Velocity km/s 1549.50 = 4 10 2 bapec norm 1.46610E-03 +/- 1.87117E-04 11 3 TBabs nH 10^22 2.93398 +/- 0.181134 ________________________________________________________________________ C-statistic = 5931.06 using 8396 PHA bins and 8390 degrees of freedom. And now with the free abundance parameter: ======================================================================== Model (bapec<1> + bapec<2>)TBabs<3> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 bapec kT keV 0.470049 +/- 3.28875E-02 2 1 bapec Abundanc 0.900062 +/- 0.188663 3 1 bapec Redshift 0.0 frozen 4 1 bapec Velocity km/s 1544.23 +/- 147.111 5 1 bapec norm 3.38155E-03 +/- 1.74952E-03 6 2 bapec kT keV 2.37960 +/- 0.195861 7 2 bapec Abundanc 0.900062 = 2 8 2 bapec Redshift 0.0 frozen 9 2 bapec Velocity km/s 1544.23 = 4 10 2 bapec norm 1.53242E-03 +/- 3.18448E-04 11 3 TBabs nH 10^22 2.97060 +/- 0.198367 ________________________________________________________________________ C-statistic = 5930.86 using 8396 PHA bins and 8389 degrees of freedom. XSPEC12>error 1. 1 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 0.505869, 0.507502 and delta stat 0.920228, 1.17117 but latest trial 0.506631 gives 0.902745 Suggest that you check this result using the steppar command. 1 0.442734 0.506685 (-0.0273111,0.0366404) yes, these limits look good, based on steppar XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.752219 1.11183 (-0.147908,0.211706) XSPEC12>error 1. 4 Parameter Confidence Range (1) 4 1391.86 1675.06 (-151.061,132.145) XSPEC12>error 1. 5 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 0.00456983, 0.00457328 and delta stat 0.982695, 1.01935 but latest trial 0.00457186 gives 0.978792 Suggest that you check this result using the steppar command. 5 0.00236516 0.00457155 (-0.0010161,0.00119028) OK, the lower limit looks fine...and so does the high one XSPEC12>error 1. 6 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 2.47763, 2.57556 and delta stat 0.323362, 1.06809 but latest trial 2.56847 gives 1.31091 Suggest that you check this result using the steppar command. 6 2.21573 2.5266 (-0.163966,0.146895) The lower limit is good, the upper one should be 2.57 XSPEC12>error 1. 10 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 0.00169384, 0.00169433 and delta stat 0.977638, 1.01639 but latest trial 0.00169408 gives 0.968198 Suggest that you check this result using the steppar command. 10 0.00136484 0.00169409 (-0.000167536,0.000161713) These look OK XSPEC12>error 1. 11 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 3.15609, 3.16007 and delta stat 0.98285, 1.02257 but latest trial 3.15785 gives 1.02337 Suggest that you check this result using the steppar command. 11 2.7716 3.15808 (-0.199028,0.18746) Looks fine XSPEC12>iplot PLT> wdata cen1a_apec2T.dat and the fluxes: with NH = 2.97: XSPEC12>flux 0.5 10. Spectrum Number: 1 Data Group Number: 1 Model Flux 0.00016948 photons (7.7077e-13 ergs/cm^2/s) range (0.50000 - 10.000 keV) Lower range bound 0.5 reset by matrix bound to 0.577208 Spectrum Number: 3 Data Group Number: 1 Model Flux 0.0001695 photons (7.709e-13 ergs/cm^2/s) range (0.57721 - 10.000 keV) XSPEC12>iplot ***Warning: Fit is not current. PLT> wdata cen1a_apec2T_noabsorb.dat PLT> exit XSPEC12>flux 0.5 10. Spectrum Number: 1 Data Group Number: 1 Model Flux 0.0042557 photons (6.4162e-12 ergs/cm^2/s) range (0.50000 - 10.000 keV) Lower range bound 0.5 reset by matrix bound to 0.577208 Spectrum Number: 3 Data Group Number: 1 Model Flux 0.0039456 photons (6.1483e-12 ergs/cm^2/s) range (0.57721 - 10.000 keV) Note: flux is almost 10 X higher OK, NOW LINE FITTING looks like Si XIV at 6.18 A is the best unblended line. Si XIII complex is the strongest, but blended. try continuum on 5.9-6.05; 6.3-6.45 ======================================================================== 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.74836E-04 +/- 3.58295E-05 ________________________________________________________________________ C-statistic = 237.38 using 352 PHA bins and 351 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.000157265 0.000193641 (-1.75715e-05,1.88043e-05) Model wgauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 wgauss sigma_l mA 30.6013 +/- 5.65257 2 1 wgauss delta_l mA -1.84277 +/- 5.70798 3 1 wgauss waveleng A 6.18200 frozen 4 1 wgauss norm 3.05540E-06 +/- 9.05836E-07 5 2 powerlaw PhoIndex 2.00000 frozen 6 2 powerlaw norm 1.75000E-04 frozen ________________________________________________________________________ C-statistic = 265.35 using 296 PHA bins and 293 degrees of freedom. XSPEC12>goodness 100 0.00% of realizations are < best fit statistic 265.35 (nosim) XSPEC12>error 1. 1 Parameter Confidence Range (1) 1 25.2674 37.9124 (-5.33397,7.31109) XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 -7.67683 3.75158 (-5.83406,5.59434) XSPEC12>error 1. 4 Parameter Confidence Range (1) 4 2.52684e-06 3.47446e-06 (-5.2856e-07,4.19055e-07) XSPEC12>iplot PLT> wdata cen1a_sixiv_6.18_gauss.dat WINDPROFILE WITH VINF = 3000 fit on pretty broad region - 6.05 to 6.30 ======================================================================== 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.83540E-02 +/- 0.445765 3 1 windprof u0 0.489087 +/- 0.104519 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 6.18200 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 3000.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.83127E-06 +/- 8.17196E-07 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.75000E-04 frozen ________________________________________________________________________ C-statistic = 263.92 using 296 PHA bins and 293 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 2 0 0.500882 (-0.0383643,0.462518) XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 2 0 0.500882 (-0.0383643,0.462518) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.395461 0.578667 (-0.0936248,0.0895805) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.40726e-06 3.22929e-06 (-4.24012e-07,3.98021e-07) XSPEC12>iplot PLT> wdata cen1a_sixiv_6.18_windprof3000.dat (no, overwrote below) refitting with 6.10 to 6.26 A range: ======================================================================== 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 4.06209E-02 +/- 0.440717 3 1 windprof u0 0.489852 +/- 0.104665 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 6.18200 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 3000.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.87636E-06 +/- 8.06014E-07 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.75000E-04 frozen ________________________________________________________________________ C-statistic = 175.13 using 188 PHA bins and 185 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 2 0 0.496225 (-0.040625,0.4556) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.393793 0.572503 (-0.0960592,0.0826513) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.39612e-06 3.22848e-06 (-4.80243e-07,3.52118e-07) NEXT, SAME MODEL BUT WITH VINF = 2500 ======================================================================== 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.145858 +/- 0.654710 3 1 windprof u0 0.317896 +/- 0.142844 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 6.18200 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2500.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.85020E-06 +/- 7.92908E-07 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.75000E-04 frozen ________________________________________________________________________ C-statistic = 174.69 using 188 PHA bins and 185 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 2 0 1.35897 (-0.145861,1.21311) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.166325 0.446071 (-0.151571,0.128175) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.37045e-06 3.19709e-06 (-4.79754e-07,3.46878e-07) XSPEC12>iplot PLT> wdata cen1a_sixiv_6.18_windprof2500.dat 94 MEG counts and 40 HEG counts Wait, just realized that I fit the wgauss model over a wider wavelength range; so now refitting wgauss + pow on the same range as I fit the windprofile model: 6.10 to 6.26 A. ======================================================================== Model wgauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 wgauss sigma_l mA 28.7121 +/- 5.40971 2 1 wgauss delta_l mA -1.81715 +/- 5.24756 3 1 wgauss waveleng A 6.18200 frozen 4 1 wgauss norm 2.87671E-06 +/- 9.02395E-07 5 2 powerlaw PhoIndex 2.00000 frozen 6 2 powerlaw norm 1.75000E-04 frozen ________________________________________________________________________ C-statistic = 177.20 using 188 PHA bins and 185 degrees of freedom. XSPEC12>iplot PLT> wdata cen1a_sixiv_6.18_gauss.dat XSPEC12>error 1. 1 Parameter Confidence Range (1) 1 24.3639 34.4534 (-4.3482,5.7413) XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 -7.09946 3.52223 (-5.28231,5.33938) XSPEC12>error 1. 4 Parameter Confidence Range (1) 4 2.47116e-06 3.36203e-06 (-4.05547e-07,4.8532e-07) NEXT, Si XIII continuum fit on 6.3-6.55; 6.80-7.0 ======================================================================== 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.84585E-04 +/- 2.53167E-05 ________________________________________________________________________ C-statistic = 405.85 using 532 PHA bins and 531 degrees of freedom. eliminating 0.1 A on either end and fitting a smaller and narrower continuum: ======================================================================== 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.83383E-04 +/- 3.41851E-05 ________________________________________________________________________ C-statistic = 210.07 using 292 PHA bins and 291 degrees of freedom. nicely similar... looks like a good, conservative range is 6.55 to 6.83 (based on v = 1500 km/s hegauss model) ... note overlap with long side of continuum fit range. maybe should refit that on 6.3-6.5, 6.85-7.0 ======================================================================== 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.79142E-04 +/- 2.84558E-05 ________________________________________________________________________ C-statistic = 318.83 using 412 PHA bins and 411 degrees of freedom. just a small change. now, fit line complex on 6.54 to 6.84 ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 1.20851 +/- 0.798496 2 1 hegauss G 1.26255 +/- 0.751853 3 1 hegauss sigma_v km/s 1473.56 +/- 265.676 4 1 hegauss delta_v km/s -93.2211 +/- 461.667 5 1 hegauss Z 14.0000 frozen 6 1 hegauss norm 7.14500E-06 +/- 1.26143E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 1.79000E-04 frozen ________________________________________________________________________ C-statistic = 336.37 using 356 PHA bins and 351 degrees of freedom. XSPEC12>error 1. 1 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 0.525964, 0.45843 and delta stat 0.982646, 1.10201 but latest trial 0.518577 gives 1.17695 Suggest that you check this result using the steppar command. 1 0.492197 2.77964 (-0.716405,1.57104) no, more like 0.46 to 2.76 XSPEC12> error 1. 2 Parameter Confidence Range (1) 2 0.635335 2.56678 (-0.627048,1.3044) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 1197.22 1891.52 (-276.321,417.977) XSPEC12>error 1. 4 Parameter Confidence Range (1) 4 -565.737 489.625 (-472.439,582.922) XSPEC12>error 1. 6 Parameter Confidence Range (1) 6 6.55891e-06 7.83205e-06 (-5.85942e-07,6.87191e-07) XSPEC12>iplot PLT> wdata cen1a_sixiii_hegauss.dat ** 20 JAN 2012 ** back to the lines again; try a bit wider wavelength range on 6.18 A (to match Sierra's for CEN1B) [6.07:6.29] use the same continuum (on 5.9:6.05; 6.30:6.45) giving norm = 1.75e-4 for the n=2 continuum ======================================================================== Model wgauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 wgauss sigma_l mA 30.3732 +/- 5.61801 2 1 wgauss delta_l mA -1.83686 +/- 5.65921 3 1 wgauss waveleng A 6.18200 frozen 4 1 wgauss norm 3.04514E-06 +/- 9.02952E-07 5 2 powerlaw PhoIndex 2.00000 frozen 6 2 powerlaw norm 1.75000E-04 frozen ________________________________________________________________________ C-statistic = 242.21 using 260 PHA bins and 257 degrees of freedom. XSPEC12>error 1. 1 Parameter Confidence Range (1) 1 25.1812 37.2651 (-5.19201,6.89191) XSPEC12>error 2.7 1 Parameter Confidence Range (2.7) 1 22.5233 42.6411 (-7.84993,12.2679) XSPEC12>error 4. 1 Parameter Confidence Range (4) 1 21.1919 46.1114 (-9.18137,15.7381) XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 -7.57477 3.71137 (-5.73791,5.54823) XSPEC12>error 2.7 2 Parameter Confidence Range (2.7) 2 -11.6204 7.31645 (-9.78354,9.15331) XSPEC12>error 4. 2 Parameter Confidence Range (4) 2 -14.5715 9.38128 (-12.7347,11.2181) XSPEC12>error 1. 4 Parameter Confidence Range (1) 4 2.52123e-06 3.45975e-06 (-5.23915e-07,4.14611e-07) XSPEC12>error 2.7 4 Parameter Confidence Range (2.7) 4 2.2576e-06 3.78896e-06 (-7.87544e-07,7.43814e-07) XSPEC12>error 4. 4 Parameter Confidence Range (4) 4 2.11654e-06 3.98079e-06 (-9.28601e-07,9.35645e-07) XSPEC12>iplot PLT> wdata cen1a_sixiv_6.18_gauss_wider.dat Re-do windprofile fits v_inf = 3000 km/s first ======================================================================== 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 4.17810E-02 +/- 0.425498 3 1 windprof u0 0.487406 +/- 9.91151E-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 6.18200 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 3000.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.51605E-06 +/- 8.80462E-07 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.75000E-04 frozen ________________________________________________________________________ C-statistic = 241.28 using 260 PHA bins and 257 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 2 0 0.669081 (-0.0413113,0.62777) XSPEC12>error 2.7 2 Parameter Confidence Range (2.7) ***Warning: Parameter pegged at hard limit: 0 2 0 1.25571 (-0.0415872,1.21412) XSPEC12>error 4. 2 Parameter Confidence Range (4) ***Warning: Parameter pegged at hard limit: 0 2 0 1.75832 (-0.0414245,1.71689) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.366648 0.595754 (-0.120785,0.108321) XSPEC12>error 2.7 3 Parameter Confidence Range (2.7) 3 0.29542 0.631769 (-0.192018,0.14433) XSPEC12>error 4. 3 Parameter Confidence Range (4) 3 0.248576 0.651558 (-0.238859,0.164123) XSPEC12>error 1. 17 Parameter Confidence Range (1) XSPEC12>steppar 17 2.3e-6 2.8e-6 25 C-Statistic Delta norm C-Statistic 17 242.37 1.0798 0 2.3e-06 242.23 0.94608 1 2.32e-06 242.11 0.81864 2 2.34e-06 241.98 0.69737 3 2.36e-06 241.87 0.5822 4 2.38e-06 241.76 0.47309 5 2.4e-06 241.66 0.36993 6 2.42e-06 241.56 0.27252 7 2.44e-06 241.47 0.18084 8 2.46e-06 241.38 0.094818 9 2.48e-06 241.3 0.014332 10 2.5e-06 241.23 -0.060604 11 2.52e-06 241.16 -0.13006 12 2.54e-06 241.09 -0.19409 13 2.56e-06 241.03 -0.25281 14 2.58e-06 240.98 -0.3062 15 2.6e-06 240.93 -0.35434 16 2.62e-06 240.89 -0.39727 17 2.64e-06 240.85 -0.43506 18 2.66e-06 240.82 -0.46777 19 2.68e-06 240.79 -0.49543 20 2.7e-06 240.77 -0.51809 21 2.72e-06 240.75 -0.5358 22 2.74e-06 240.74 -0.54861 23 2.76e-06 240.73 -0.55656 24 2.78e-06 240.73 -0.55971 25 2.8e-06 XSPEC12>steppar 17 2.8e-6 3.0e-6 20 C-Statistic Delta norm C-Statistic 17 240.73 -0.5543 0 2.8e-06 240.73 -0.55459 1 2.81e-06 240.73 -0.55367 2 2.82e-06 240.73 -0.55154 3 2.83e-06 240.74 -0.54821 4 2.84e-06 240.74 -0.54368 5 2.85e-06 240.75 -0.53795 6 2.86e-06 240.75 -0.53104 7 2.87e-06 240.76 -0.52294 8 2.88e-06 240.77 -0.51367 9 2.89e-06 240.78 -0.50322 10 2.9e-06 240.79 -0.49161 11 2.91e-06 240.81 -0.47884 12 2.92e-06 240.82 -0.46492 13 2.93e-06 240.84 -0.44984 14 2.94e-06 240.85 -0.43362 15 2.95e-06 240.87 -0.41627 16 2.96e-06 240.89 -0.39778 17 2.97e-06 240.91 -0.37816 18 2.98e-06 240.93 -0.35742 19 2.99e-06 240.95 -0.33562 20 3e-06 So, the best fit norm is ~2.80e-6 (but for some reason -- poor data quality -- "fit" always goes to a worse fit with norm ~ 2.5e-6) So for now, freeze norm = 2.80e-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 3.60880E-02 +/- 0.451678 3 1 windprof u0 0.489805 +/- 0.101421 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 6.18200 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 3000.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.80000E-06 frozen 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.75000E-04 frozen ________________________________________________________________________ C-statistic = 240.73 using 260 PHA bins and 258 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 2 0 0.499906 (-0.0361679,0.463738) XSPEC12>error 2.7 2 Parameter Confidence Range (2.7) ***Warning: Parameter pegged at hard limit: 0 2 0 1.0752 (-0.0362466,1.03896) XSPEC12>error 4. 2 Parameter Confidence Range (4) ***Warning: Parameter pegged at hard limit: 0 2 0 1.54804 (-0.0362544,1.51178) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.394222 0.577806 (-0.0955701,0.0880134) XSPEC12>error 2.7 3 Parameter Confidence Range (2.7) 3 0.315938 0.619964 (-0.173855,0.130172) XSPEC12>error 4. 3 Parameter Confidence Range (4) ***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.27002 0.641201 (-0.219773,0.151408) XSPEC12>steppar 17 2.3e-6 3.3e-6 20 C-Statistic Delta norm C-Statistic 17 242.37 1.0798 0 2.3e-06 242.04 0.75723 1 2.35e-06 241.76 0.47309 2 2.4e-06 241.51 0.226 3 2.45e-06 241.3 0.014499 4 2.5e-06 241.12 -0.16236 5 2.55e-06 240.98 -0.30462 6 2.6e-06 240.87 -0.41625 7 2.65e-06 240.79 -0.49512 8 2.7e-06 240.74 -0.54253 9 2.75e-06 240.73 -0.5593 10 2.8e-06 240.74 -0.5457 11 2.85e-06 240.78 -0.50446 12 2.9e-06 240.85 -0.43425 13 2.95e-06 240.95 -0.33609 14 3e-06 241.08 -0.21021 15 3.05e-06 241.23 -0.058971 16 3.1e-06 241.41 0.11869 17 3.15e-06 241.61 0.32194 18 3.2e-06 241.84 0.55031 19 3.25e-06 242.09 0.80227 20 3.3e-06 XSPEC12>iplot PLT> wdata cen1a_sixiv_6.18_windprof3000_wider.dat ** next, with vinf = 2500 km/s ======================================================================== 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.139368 +/- 0.669815 3 1 windprof u0 0.316716 +/- 0.142798 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 6.18200 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2500.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.79651E-06 +/- 8.05371E-07 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.75000E-04 frozen ________________________________________________________________________ C-statistic = 240.32 using 260 PHA bins and 257 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 2 0 1.32428 (-0.135607,1.18867) XSPEC12>error 2.7 2 ***Warning: Identical values of the parameter give different values of the statistic. Please check your result for the high end of the confidence range. 2 0 5.0822 (-0.136775,4.94542) XSPEC12>steppar 2 0.01 2.0 50 C-Statistic Delta taustar C-Statistic 2 240.38 -0.62053 0 0.01 240.35 -0.65289 1 0.0498 240.33 -0.67203 2 0.0896 240.32 -0.6806 3 0.1294 240.32 -0.68073 4 0.1692 240.33 -0.67403 5 0.209 240.34 -0.66171 6 0.2488 240.36 -0.64474 7 0.2886 240.38 -0.62388 8 0.3284 240.4 -0.59974 9 0.3682 240.43 -0.57283 10 0.408 240.46 -0.54358 11 0.4478 240.49 -0.51233 12 0.4876 240.52 -0.47938 13 0.5274 240.56 -0.445 14 0.5672 240.59 -0.40938 15 0.607 240.63 -0.3727 16 0.6468 240.67 -0.33512 17 0.6866 240.71 -0.29693 18 0.7264 240.74 -0.25834 19 0.7662 240.78 -0.21946 20 0.806 240.82 -0.1804 21 0.8458 240.86 -0.14125 22 0.8856 240.9 -0.1021 23 0.9254 240.94 -0.063 24 0.9652 bigger range: XSPEC12>steppar 2 0.01 6.0 50 C-Statistic Delta taustar C-Statistic 2 240.38 -0.62204 0 0.01 240.32 -0.68216 1 0.1298 240.34 -0.66305 2 0.2496 240.4 -0.60078 3 0.3694 240.49 -0.51306 4 0.4892 240.59 -0.40981 5 0.609 240.71 -0.29711 6 0.7288 240.82 -0.18018 7 0.8486 240.94 -0.062444 8 0.9684 241.06 0.054142 9 1.0882 241.17 0.16825 10 1.208 241.28 0.27908 11 1.3278 241.39 0.38605 12 1.4476 241.49 0.48448 13 1.5674 241.59 0.58248 14 1.6872 241.68 0.67642 15 1.807 241.77 0.76968 16 1.9268 241.86 0.85411 17 2.0466 241.94 0.93246 18 2.1664 242.01 1.0051 19 2.2862 242.71 1.7048 20 2.406 242.76 1.7601 21 2.5258 242.82 1.812 22 2.6456 242.86 1.8596 23 2.7654 242.91 1.9031 24 2.8852 242.95 1.9427 25 3.005 242.43 1.4264 26 3.1248 242.47 1.4663 27 3.2446 242.51 1.5041 28 3.3644 242.54 1.5402 29 3.4842 242.58 1.5748 30 3.604 242.64 1.6334 31 3.7238 242.67 1.6662 32 3.8436 242.7 1.6987 33 3.9634 242.78 1.7721 34 4.0832 242.72 1.7146 35 4.203 242.75 1.7443 36 4.3228 242.78 1.7786 37 4.4426 242.85 1.8442 38 4.5624 242.86 1.8584 39 4.6822 242.9 1.891 40 4.802 242.93 1.9291 41 4.9218 242.97 1.9676 42 5.0416 243.01 2.0065 43 5.1614 243.05 2.0457 44 5.2812 243.08 2.0717 45 5.401 243.1 2.0917 46 5.5208 243.14 2.1326 47 5.6406 243.18 2.1735 48 5.7604 243.22 2.2145 49 5.8802 243.26 2.2555 50 6 So, DeltaC = 2.7 comes at 5.19 XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.168737 0.444811 (-0.147365,0.128709) XSPEC12>error 2.7 3 Parameter Confidence Range (2.7) ***Warning: Parameter pegged at hard limit: 0.1 3 0 0.510488 (-0.316109,0.194379) XSPEC12>error 4. 3 Parameter Confidence Range (4) ***Warning: Parameter pegged at hard limit: 0.1 3 0 0.542387 (-0.316112,0.226275) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.38098e-06 3.1927e-06 (-4.16032e-07,3.95689e-07) XSPEC12>error 2.7 17 Parameter Confidence Range (2.7) 17 2.14014e-06 3.4784e-06 (-6.56872e-07,6.81396e-07) XSPEC12>error 4. 17 Parameter Confidence Range (4) 17 2.01196e-06 3.64076e-06 (-7.8505e-07,8.43757e-07) XSPEC12>iplot PLT> wdata cen1a_sixiv_6.18_windprof2500_wider.dat