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