N VII Ly alpha 24.781 *note: blend with N VI He-beta (24.89) ignore this for now... continuum on: 24.40:24.55, 25.0:25.2 ======================================================================== 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.13366E-03 +/- 8.06839E-04 ________________________________________________________________________ C-statistic = 74.58 using 136 PHA bins and 135 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>plot XSPEC12>goodness 1000 nosim 61.50% of realizations are < best fit statistic 74.58 XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 0.001681 0.002659 (-0.000453,0.000525) fitting line on: 24.58:24.98 ======================================================================== 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.973802 +/- 0.516089 3 1 windprof u0 0.328008 +/- 4.71463E-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" 24.7810 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 7.59753E-04 +/- 9.35080E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.13000E-03 frozen ________________________________________________________________________ C-statistic = 172.86 using 158 PHA bins and 155 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>plot XSPEC12>goodness 1000 nosim 73.50% of realizations are < best fit statistic 172.86 XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 0.518641 1.606418 (-0.455162,0.632615) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 0.312186 2.192212 (-0.661617,1.218409) XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) 3 0.281080 0.376870 (-0.046928,0.048862) XSPEC12>error 2.7 3 Parameter Confidence Range (2.700000) 3 0.245910 0.404663 (-0.082098,0.076655) XSPEC12>error 1. 17 Parameter Confidence Range (1.000000) 17 0.000712 0.000813 (-0.000047,0.000053) XSPEC12>error 2.7 17 Parameter Confidence Range (2.700000) 17 0.000682 0.000847 (-0.000078,0.000087) made a uo-taustar confidence grid; but using standard range of grid params the best C value on the grid is 184.8 because this line favors much lower u_o values saved the files associated with this exercise with the "stdrange" in their names; and remade the grid on: XSPEC12>log /home/cohen/xspec/2007/zpup/nvii_24781_smooth_35by35.log Logging to file:/home/cohen/xspec/2007/zpup/nvii_24781_smooth_35by35.log XSPEC12>steppar 2 0.01 4.5 35 3 0.08 0.59 35 Source File: acis640_pha2.fits{9} Net count rate (cts/s) for Spectrum:1 3.238e-03 +/- 3.737e-04 Assigned to Data Group 1 and Plot Group 1 Noticed Channels: 3398-3476 Telescope: CHANDRA Instrument: ACIS Channel Type: PI Exposure Time: 6.764e+04 sec Using Response (RMF) File meg_m1.rmf for Source 1 Using Auxiliary Response (ARF) File acis640MEG_-1_garf.fits Source spectrum counts: 219 Model predicted rate: 0.0 Source File: acis640_pha2.fits{10} Net count rate (cts/s) for Spectrum:2 7.836e-04 +/- 2.849e-04 Assigned to Data Group 1 and Plot Group 2 Noticed Channels: 3398-3476 Telescope: CHANDRA Instrument: ACIS Channel Type: PI Exposure Time: 6.764e+04 sec Using Response (RMF) File meg_p1.rmf for Source 1 Using Auxiliary Response (ARF) File acis640MEG_1_garf.fits Source spectrum counts: 53 Model predicted rate: 0.0 ---- Now, we're fitting two windprofile models, including one to the N VI beta line at 24.89; First we tie their uo and taustars together, and also tie the norms together but scaling the 2nd one to be 1/5 the first one. ======================================================================== Model windprof<1> + windprof<2> + powerlaw<3> 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.37183 +/- 1.03492 3 1 windprof u0 0.352047 +/- 7.23797E-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" 24.7810 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 6.48417E-04 +/- 8.03214E-05 18 2 windprof q 0.0 frozen 19 2 windprof taustar 2.37183 = windprof[1]:taustar 20 2 windprof u0 0.352047 = windprof[1]:u0 21 2 windprof h 0.0 frozen 22 2 windprof tau0star 0.0 frozen 23 2 windprof beta 1.00000 frozen 24 2 windprof betaSob 0.0 frozen 25 2 windprof numerica 0 frozen 26 2 windprof anisotro 0 frozen 27 2 windprof rosselan 0 frozen 28 2 windprof expansio 0 frozen 29 2 windprof thick 0 frozen 30 2 windprof waveleng "A" 24.8900 frozen 31 2 windprof shift "mA" 0.0 frozen 32 2 windprof velocity "km/s" 2250.00 frozen 33 2 windprof verbose 0 frozen 34 2 windprof norm 1.29683E-04 = 0.2*windprof[1]:norm 35 3 powerlaw PhoIndex 2.00000 frozen 36 3 powerlaw norm 2.13000E-03 frozen ________________________________________________________________________ C-statistic = 176.10 using 158 PHA bins and 155 degrees of freedom. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>goodness 1000 nosim 74.20% of realizations are < best fit statistic 176.10 So, the results hardly changed... XSPEC12>goodness 1000 nosim 74.20% of realizations are < best fit statistic 176.10 XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 1.473133 3.710390 (-0.898695,1.338562) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 1.043974 4.956129 (-1.327851,2.584304) XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) 3 0.290944 0.419905 (-0.061104,0.067858) XSPEC12>error 2.7 3 Parameter Confidence Range (2.700000) 3 0.250809 0.506673 (-0.101239,0.154625) XSPEC12>error 1. 17 Parameter Confidence Range (1.000000) 17 0.000607 0.000692 (-0.000042,0.000044) XSPEC12>error 2.7 17 Parameter Confidence Range (2.700000) 17 0.000581 0.000721 (-0.000068,0.000073) PLT> wdata /home/cohen/xspec/2007/zpup/nvii_24781_smooth_twolines02.dat XSPEC12>log /home/cohen/xspec/2007/zpup/nvii_24781_smooth_twolines02_35by35.log Logging to file:/home/cohen/xspec/2007/zpup/nvii_24781_smooth_twolines02_35by35.log XSPEC12>steppar 2 0.01 4.5 35 3 0.08 0.59 35 Re-running the exercise with a normalization ratio of 0.4 (instead of 0.2): ======================================================================== Model windprof<1> + windprof<2> + powerlaw<3> 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.37996 +/- 1.88941 3 1 windprof u0 0.472287 +/- 0.266104 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" 24.7810 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 5.61943E-04 +/- 7.00585E-05 18 2 windprof q 0.0 frozen 19 2 windprof taustar 4.37996 = windprof[1]:taustar 20 2 windprof u0 0.472287 = windprof[1]:u0 21 2 windprof h 0.0 frozen 22 2 windprof tau0star 0.0 frozen 23 2 windprof beta 1.00000 frozen 24 2 windprof betaSob 0.0 frozen 25 2 windprof numerica 0 frozen 26 2 windprof anisotro 0 frozen 27 2 windprof rosselan 0 frozen 28 2 windprof expansio 0 frozen 29 2 windprof thick 0 frozen 30 2 windprof waveleng "A" 24.8900 frozen 31 2 windprof shift "mA" 0.0 frozen 32 2 windprof velocity "km/s" 2250.00 frozen 33 2 windprof verbose 0 frozen 34 2 windprof norm 2.24777E-04 = 0.4*windprof[1]:norm 35 3 powerlaw PhoIndex 2.00000 frozen 36 3 powerlaw norm 2.13000E-03 frozen ________________________________________________________________________ C-statistic = 178.65 using 158 PHA bins and 155 degrees of freedom. XSPEC12>goodness 1000 nosim 77.50% of realizations are < best fit statistic 178.65 PLT> wdata /home/cohen/xspec/2007/zpup/nvii_24781_smooth_twolines04.dat XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 2.980393 6.055502 (-1.407823,1.667286) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 2.277181 7.540812 (-2.115105,3.148526) XSPEC12>plot XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) ***Warning: Parameter pegged at hard limit: 0.99 3 0.339865 0.506673 (-0.133642,0.033166) XSPEC12>error 2.7 3 Parameter Confidence Range (2.700000) ***Warning: Parameter pegged at hard limit: 0.99 3 0.288157 0.506673 (-0.185635,0.032882) XSPEC12>error 1. 17 Parameter Confidence Range (1.000000) 17 0.000526 0.000600 (-0.000037,0.000037) XSPEC12>error 2.7 17 Parameter Confidence Range (2.700000) 17 0.000503 0.000625 (-0.000059,0.000062) XSPEC12>log /home/cohen/xspec/2007/zpup/nvii_24781_smooth_twolines04_35by35.log Logging to file:/home/cohen/xspec/2007/zpup/nvii_24781_smooth_twolines04_35by35.log XSPEC12>steppar 2 0.01 4.5 35 3 0.08 0.59 35 ---- Now, let's do a weaker N VI line: 0.1 normalization ratio ======================================================================== Model windprof<1> + windprof<2> + powerlaw<3> 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.54026 +/- 0.726460 3 1 windprof u0 0.333809 +/- 5.59526E-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" 24.7810 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 7.01191E-04 +/- 8.64976E-05 18 2 windprof q 0.0 frozen 19 2 windprof taustar 1.54026 = windprof[1]:taustar 20 2 windprof u0 0.333809 = windprof[1]:u0 21 2 windprof h 0.0 frozen 22 2 windprof tau0star 0.0 frozen 23 2 windprof beta 1.00000 frozen 24 2 windprof betaSob 0.0 frozen 25 2 windprof numerica 0 frozen 26 2 windprof anisotro 0 frozen 27 2 windprof rosselan 0 frozen 28 2 windprof expansio 0 frozen 29 2 windprof thick 0 frozen 30 2 windprof waveleng "A" 24.8900 frozen 31 2 windprof shift "mA" 0.0 frozen 32 2 windprof velocity "km/s" 2250.00 frozen 33 2 windprof verbose 0 frozen 34 2 windprof norm 7.01191E-05 = 0.1*windprof[1]:norm 35 3 powerlaw PhoIndex 2.00000 frozen 36 3 powerlaw norm 2.13000E-03 frozen ________________________________________________________________________ C-statistic = 174.42 using 158 PHA bins and 155 degrees of freedom. XSPEC12>goodness 1000 nosim 72.40% of realizations are < best fit statistic 174.42 PLT> wdata /home/cohen/xspec/2007/zpup/nvii_24781_smooth_twolines01.dat XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 0.904686 2.457330 (-0.635190,0.917453) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 0.605089 3.293975 (-0.934603,1.754282) XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) 3 0.282092 0.385956 (-0.051656,0.052208) XSPEC12>error 2.7 3 Parameter Confidence Range (2.700000) 3 0.245042 0.419441 (-0.088705,0.085693) XSPEC12>error 1. 17 Parameter Confidence Range (1.000000) 17 0.000656 0.000749 (-0.000045,0.000048) XSPEC12>error 2.7 17 Parameter Confidence Range (2.700000) 17 0.000628 0.000780 (-0.000073,0.000079) Logging to file:/home/cohen/xspec/2007/zpup/nvii_24781_smooth_twolines01_35by35.log XSPEC12>steppar 2 0.01 4.5 35 3 0.08 0.59 35