18oct11 dhc Re-analyzing the xi Per spectrum, looking into the effects of the = 70km/s radial velocity of the star Fe XVII @ 15.014A Emma found a continuum level of 8.75e-4 (n=2) and fit the line on 14.85 to 15.05 I'm going to use less of the blue continuum and more of the red wing 14.88 to 15.08 ======================================================================== 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.156076 +/- 0.126474 3 1 windprof u0 0.652852 +/- 3.29305E-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 15.0140 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 1 frozen 17 1 windprof norm 9.55654E-05 +/- 1.42144E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 8.75000E-04 frozen ________________________________________________________________________ C-statistic = 96.17 using 78 PHA bins and 75 degrees of freedom. Yes; very consistent with Emma's result XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.0495155 0.304287 (-0.10655,0.148222) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.60572 0.690159 (-0.0471402,0.0372986) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 8.89505e-05 0.000103409 (-6.61492e-06,7.84316e-06) XSPEC12>iplot PLT> wdata 15.014.dat (Note: I only am fitting the MEG data here; Emma fit HEG too) Next, look at the effect of a shift (associated with the RV of the star, of 70 km/s -- not corrected to the solar system barycenter) ======================================================================== 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.289693 +/- 0.169727 3 1 windprof u0 0.646576 +/- 3.55252E-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 15.0140 frozen 14 1 windprof shift mA 3.50000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 9.61596E-05 +/- 1.43478E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 8.75000E-04 frozen ________________________________________________________________________ C-statistic = 96.38 using 78 PHA bins and 75 degrees of freedom. So, taustar doubles with the 3.5 mA shift associated with 70 km/s! XSPEC12>iplot PLT> wdata 15.014_shift3.5mA.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.144201 0.481939 (-0.145547,0.192191) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.598269 0.688362 (-0.048264,0.0418293) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 8.94884e-05 0.000104092 (-6.67107e-06,7.93267e-06) OK, before going on, let's read in the HEG data too, and fit it for this line, and any others where it may contribute significant information. We can use the above fits to make a direct comparison. (my confidence limits are very, very similar to Emma's). [cohen@max xper]$ mv 15.014.dat 15.014_meg.dat [cohen@max xper]$ mv 15.014_shift3.5mA.dat 15.014_shift3.5mA_meg.dat OK, loading in all 4 spectra [14.88-15.08] ======================================================================== 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.122041 +/- 0.101043 3 1 windprof u0 0.645530 +/- 2.87465E-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 15.0140 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 9.87310E-05 +/- 1.30029E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 8.75000E-04 frozen ________________________________________________________________________ C-statistic = 229.55 using 236 PHA bins and 233 degrees of freedom. Hmmm, taustar somewhat lower than Emma's... XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.0331668 0.241649 (-0.088823,0.119659) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.602122 0.681869 (-0.0434588,0.0362878) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 9.26731e-05 0.000105931 (-6.0496e-06,7.20826e-06) XSPEC12>iplot PLT> wdata 15.014.dat Before moving on, do a quick check of the wavelength range effect - XSPEC12>ignore 1-4:15.05-** ======================================================================== 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.68696E-02 +/- 0.201585 3 1 windprof u0 0.639189 +/- 3.68565E-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 15.0140 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 1.02330E-04 +/- 1.86078E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 8.75000E-04 frozen ________________________________________________________________________ C-statistic = 195.88 using 200 PHA bins and 197 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 2 0 0.189953 (-0.0479125,0.142041) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.563214 0.683322 (-0.0766052,0.0435026) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 9.35594e-05 0.000115621 (-8.61463e-06,1.34465e-05) Back to the regular wavelength range: XSPEC12>notice 1-4:15.-** XSPEC12>ignore 1-4:15.08-** Now, applying the 3.5 mA shift: ======================================================================== 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.242847 +/- 0.137938 3 1 windprof u0 0.638771 +/- 3.06682E-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 15.0140 frozen 14 1 windprof shift mA 3.50000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 9.98807E-05 +/- 1.31261E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 8.75000E-04 frozen ________________________________________________________________________ C-statistic = 229.58 using 236 PHA bins and 233 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.120427 0.398971 (-0.123055,0.15549) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.59878 0.679493 (-0.0394764,0.0412371) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 9.31989e-05 0.000106571 (-6.16224e-06,7.2095e-06) XSPEC12>iplot PLT> wdata 15.014_shift3.5mA.dat NEXT, NE X 12.134 Emma excluded a lot of the red wing; I'm going to experiment with the noticed wavelength range I'll use both MEG and HEG and use Emma's continuum model - OK, reproduced her results. - short wavelength limit should be 12.01 A - testing long limit of 12.20 --> results are very similar to 12.0-12.15 fit Looking at the data...I think, let's cut things off at 12.18 A So, (12.01:12.18) ======================================================================== 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.538379 +/- 0.199509 3 1 windprof u0 0.648332 +/- 3.50937E-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 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.99124E-05 +/- 5.01594E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 3.58000E-04 frozen ________________________________________________________________________ C-statistic = 184.47 using 200 PHA bins and 197 degrees of freedom. So, that's equivalent to Emma's fit (taustar is 20% higher) Just for comparison: (12.01:12.15) is even higher, hmmm... I see best-fit taustars jumping around (0.4...0.7) well within the 68% confidence limits OK, now fit with a shift. ======================================================================== 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.776362 +/- 0.256188 3 1 windprof u0 0.644207 +/- 4.13858E-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 2.80000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 4.02330E-05 +/- 5.07857E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 3.58000E-04 frozen ________________________________________________________________________ C-statistic = 185.49 using 200 PHA bins and 197 degrees of freedom. XSPEC12>iplot PLT> wdata 12.134_shift2.8mA.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.534518 1.05025 (-0.241844,0.273891) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.59728 0.682839 (-0.0469272,0.0386316) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 3.80767e-05 4.32189e-05 (-2.15632e-06,2.9859e-06) ********************************** OK, going to use 57 km/s as our standard for xi Per; so redo the above line fits (just the ones with the shifts) Ne X 12.1339 A: on (12.01:12.18) ======================================================================== 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.730293 +/- 0.245223 3 1 windprof u0 0.644669 +/- 4.00649E-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 2.30000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 4.01773E-05 +/- 5.06722E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 3.58000E-04 frozen ________________________________________________________________________ C-statistic = 185.28 using 200 PHA bins and 197 degrees of freedom. XSPEC12>iplot PLT> wdata 12.134_shift2.3mA.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.499233 0.993741 (-0.23106,0.263449) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.599034 0.681561 (-0.0456352,0.0368917) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 3.80142e-05 4.31418e-05 (-2.16306e-06,2.96447e-06) Fe XVII @ 15.014 A on (14.88:15.08) ======================================================================== 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.216411 +/- 0.130908 3 1 windprof u0 0.639617 +/- 3.02729E-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 15.0140 frozen 14 1 windprof shift mA 2.80000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 9.92222E-05 +/- 1.31049E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 8.75000E-04 frozen ________________________________________________________________________ C-statistic = 229.59 using 236 PHA bins and 233 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_shift2.8mA.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.0996812 0.364331 (-0.116729,0.14792) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.599279 0.680012 (-0.0403387,0.0403941) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 9.30669e-05 0.000106441 (-6.15534e-06,7.21872e-06) Mg XII @ 8.421 (8.34:8.50) -- somewhat narrower than Emma'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 1.56407E-05 +/- 0.448943 3 1 windprof u0 0.647955 +/- 8.55252E-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 8.42100 frozen 14 1 windprof shift mA 1.60000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.61324E-06 +/- 9.24637E-07 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.97000E-04 frozen ________________________________________________________________________ C-statistic = 145.88 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.198054 (-5.14745e-06,0.198048) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 2.17867e-06 3.09657e-06 (-4.345e-07,4.83404e-07) XSPEC12>error 1. 3 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 0.647992, 0.562467 and delta stat 0, 1.02108 but latest trial 0.563355 gives 1.07718 Suggest that you check this result using the steppar command. 3 0.60523 0.716513 (-0.0427782,0.0685048) upper OK, but lower should be 0.52 XSPEC12>iplot PLT> wdata 8.421_shift1.9mA.dat extra: XSPEC12>error 2.7 2 Parameter Confidence Range (2.7) ***Warning: Parameter pegged at hard limit: 0 2 0 0.507802 (-5.42413e-12,0.507802) XSPEC12>error 4. 2 Parameter Confidence Range (4) ***Warning: Parameter pegged at hard limit: 0 2 0 0.763166 (-4.26904e-13,0.763166) Ne IX @ 11.544 A (11.44:11.65) -- narrower than Emma'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.150804 +/- 0.231703 3 1 windprof u0 0.667285 +/- 5.88138E-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 11.5440 frozen 14 1 windprof shift mA 2.20000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 9.73759E-06 +/- 2.47674E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 4.63000E-04 frozen ________________________________________________________________________ C-statistic = 178.86 using 244 PHA bins and 241 degrees of freedom. XSPEC12>iplot PLT> wdata 11.544_shift2.2mA.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 2 0 0.501645 (-0.150807,0.350838) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.57218 0.733169 (-0.0951056,0.065883) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 8.24311e-06 1.08904e-05 (-1.49446e-06,1.15279e-06) extra: XSPEC12>error 4. 17 Parameter Confidence Range (4) 17 7.09083e-06 1.22718e-05 (-2.64674e-06,2.53423e-06) XSPEC12>error 4. 2 Parameter Confidence Range (4) ***Warning: Parameter pegged at hard limit: 0 Apparent non-monotonicity in statistic space detected. Current bracket values 0.150807, 1.71938 and delta stat 0, 4.53585 but latest trial 1.56612 gives 4.60596 Suggest that you check this result using the steppar command. 2 0 0.935093 (-0.150807,0.784286) better upper limit is 1.30 Si XIII @ 6.7 A (6.58:6.81) ======================================================================== Model hewind<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hewind q 0.0 frozen 2 1 hewind taustar 2.39653E-27 +/- 0.266405 3 1 hewind u0 0.679648 +/- 4.18497E-02 4 1 hewind h 0.0 frozen 5 1 hewind tau0star 0.0 frozen 6 1 hewind beta 1.00000 frozen 7 1 hewind betaSob 0.0 frozen 8 1 hewind G 1.25631 +/- 0.222822 9 1 hewind numerica 0 frozen 10 1 hewind anisotro 0 frozen 11 1 hewind rosselan 0 frozen 12 1 hewind expansio 0 frozen 13 1 hewind thick 0 frozen 14 1 hewind Z 14 frozen 15 1 hewind shift mA 1.30000 frozen 16 1 hewind velocity km/s 2450.00 frozen 17 1 hewind phiratio 7.83000 frozen 18 1 hewind verbose 0 frozen 19 1 hewind norm 9.67377E-06 +/- 1.36610E-06 20 2 powerlaw PhoIndex 2.00000 frozen 21 2 powerlaw norm 1.29000E-04 frozen ________________________________________________________________________ C-statistic = 297.19 using 272 PHA bins and 268 degrees of freedom. XSPEC12>iplot PLT> wdata 6.7_shift1.3mA.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 2 0 0.0290608 (-4.69782e-28,0.0290608) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.639793 0.715039 (-0.0396884,0.0355574) XSPEC12>error 1. 8 Parameter Confidence Range (1) 8 1.07215 1.49983 (-0.184084,0.243596) XSPEC12>error 1. 19 Parameter Confidence Range (1) 19 8.95436e-06 1.02962e-05 (-7.19372e-07,6.2248e-07) extra: XSPEC12>error 2.7 2 Parameter Confidence Range (2.7) ***Warning: Parameter pegged at hard limit: 0 2 0 0.0877309 (-1.82112e-29,0.0877309) XSPEC12>error 4. 2 Parameter Confidence Range (4) ***Warning: Parameter pegged at hard limit: 0 2 0 0.13806 (-3.1591e-30,0.13806) Mg XI @ 9.2 A (9.08:9.40) ======================================================================== Model hewind<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hewind q 0.0 frozen 2 1 hewind taustar 0.400238 +/- 0.226088 3 1 hewind u0 0.629022 +/- 3.17608E-02 4 1 hewind h 0.0 frozen 5 1 hewind tau0star 0.0 frozen 6 1 hewind beta 1.00000 frozen 7 1 hewind betaSob 0.0 frozen 8 1 hewind G 0.956976 +/- 0.129295 9 1 hewind numerica 0 frozen 10 1 hewind anisotro 0 frozen 11 1 hewind rosselan 0 frozen 12 1 hewind expansio 0 frozen 13 1 hewind thick 0 frozen 14 1 hewind Z 12 frozen 15 1 hewind shift mA 1.70000 frozen 16 1 hewind velocity km/s 2450.00 frozen 17 1 hewind phiratio 84.5000 frozen 18 1 hewind verbose 0 frozen 19 1 hewind norm 2.42126E-05 +/- 2.41743E-06 20 2 powerlaw PhoIndex 2.00000 frozen 21 2 powerlaw norm 2.60000E-04 frozen ________________________________________________________________________ C-statistic = 367.75 using 380 PHA bins and 376 degrees of freedom. XSPEC12>iplot PLT> wdata 9.2_shift1.7mA.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.218528 0.624001 (-0.18126,0.224213) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.597345 0.658517 (-0.0316521,0.0295198) XSPEC12>error 1. 8 Parameter Confidence Range (1) 8 0.834733 1.09062 (-0.122172,0.133714) XSPEC12>error 1. 19 Parameter Confidence Range (1) 19 2.28632e-05 2.53107e-05 (-1.34934e-06,1.09815e-06) Ne IX @ 13.5 (13.32:13.55) -- note excluding about half of the complex due to potential blending problems ======================================================================== Model hewind<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hewind q 0.0 frozen 2 1 hewind taustar 4.86847E-07 +/- 0.206411 3 1 hewind u0 0.614444 +/- 2.74145E-02 4 1 hewind h 0.0 frozen 5 1 hewind tau0star 0.0 frozen 6 1 hewind beta 1.00000 frozen 7 1 hewind betaSob 0.0 frozen 8 1 hewind G 0.965856 +/- 0.153634 9 1 hewind numerica 0 frozen 10 1 hewind anisotro 0 frozen 11 1 hewind rosselan 0 frozen 12 1 hewind expansio 0 frozen 13 1 hewind thick 0 frozen 14 1 hewind Z 10 frozen 15 1 hewind shift mA 2.60000 frozen 16 1 hewind velocity km/s 2450.00 frozen 17 1 hewind phiratio 733.000 frozen 18 1 hewind verbose 0 frozen 19 1 hewind norm 1.16772E-04 +/- 1.50415E-05 20 2 powerlaw PhoIndex 2.00000 frozen 21 2 powerlaw norm 6.19000E-04 frozen ________________________________________________________________________ C-statistic = 260.81 using 272 PHA bins and 268 degrees of freedom. iXSPEC12>iplot PLT> wdata 13.5_shift2.6mA.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 2 0 0.0354898 (-1.99014e-07,0.0354896) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.587473 0.635698 (-0.0266503,0.0215744) XSPEC12>error 1. 8 Parameter Confidence Range (1) 8 0.828726 1.13995 (-0.136363,0.174858) XSPEC12>error 1. 19 Parameter Confidence Range (1) 19 0.000110339 0.000125038 (-6.41054e-06,8.28864e-06) extra: XSPEC12>error 2.7 2 Parameter Confidence Range (2.7) ***Warning: Parameter pegged at hard limit: 0 2 0 0.0959856 (-1.59716e-10,0.0959856) XSPEC12>error 4. 2 Parameter Confidence Range (4) ***Warning: Parameter pegged at hard limit: 0 2 0 0.142648 (-1.9515e-11,0.142648) NOW HAVE JUST MEG DATA, NO HEG (for all lines longer than 15.014 A) Emma fit 16.78, 17.05, 17.10 all together, I'll do 16.78 separately Fe XVII @ 16.780 (16.63:16.91) ======================================================================== 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.249613 +/- 0.237760 3 1 windprof u0 0.744674 +/- 4.89696E-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 16.7800 frozen 14 1 windprof shift mA 3.20000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.94939E-05 +/- 1.00818E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 5.65000E-04 frozen ________________________________________________________________________ C-statistic = 109.23 using 110 PHA bins and 107 degrees of freedom. XSPEC12>iplot PLT> wdata 16.780_shift3.2mA.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.0589727 0.563115 (-0.19064,0.313502) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.681649 0.795363 (-0.0630246,0.0506893) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 3.35069e-05 4.40548e-05 (-5.98711e-06,4.56081e-06) O VIII @ 18.9689 (18.80:19.13) ======================================================================== 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.404445 +/- 0.224827 3 1 windprof u0 0.600131 +/- 3.93902E-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 18.9689 frozen 14 1 windprof shift mA 3.60000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 1.91124E-04 +/- 2.81932E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 4.70000E-04 frozen ________________________________________________________________________ C-statistic = 107.39 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata 18.969_shift3.6mA.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.213476 0.655405 (-0.190964,0.250965) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.558571 0.63261 (-0.0415582,0.0324807) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000176389 0.000205273 (-1.47355e-05,1.41487e-05) Fe XVII @ 17.051,17.096 (16.91:17.25) with .10/.05 ratio set to 0.8: ======================================================================== 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.02684 +/- 0.485924 3 1 windprof u0 0.574698 +/- 7.42933E-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 17.0510 frozen 14 1 windprof shift mA 3.20000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 6.45656E-05 +/- 9.78583E-06 18 2 windprof q 0.0 frozen 19 2 windprof taustar 1.02684 = 2 20 2 windprof u0 0.574698 = 3 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 17.0960 frozen 31 2 windprof shift mA 3.20000 frozen 32 2 windprof velocity km/s 2450.00 frozen 33 2 windprof verbose 0 frozen 34 2 windprof norm 5.16525E-05 = 0.8*17 35 3 powerlaw PhoIndex 2.00000 frozen 36 3 powerlaw norm 5.65000E-04 frozen ________________________________________________________________________ C-statistic = 128.83 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata 17.05_17.10_ratio0.8_shift3.2mA.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.608464 2.31142 (-0.41838,1.28457) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.500534 0.667674 (-0.0741635,0.0929761) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 5.98633e-05 6.96972e-05 (-4.70232e-06,5.13156e-06) now with a ratio of 0.9: ======================================================================== 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.16861 +/- 0.536469 3 1 windprof u0 0.578142 +/- 8.14742E-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 17.0510 frozen 14 1 windprof shift mA 3.20000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 6.12698E-05 +/- 9.27065E-06 18 2 windprof q 0.0 frozen 19 2 windprof taustar 1.16861 = 2 20 2 windprof u0 0.578142 = 3 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 17.0960 frozen 31 2 windprof shift mA 3.20000 frozen 32 2 windprof velocity km/s 2450.00 frozen 33 2 windprof verbose 0 frozen 34 2 windprof norm 5.51428E-05 = 0.9*17 35 3 powerlaw PhoIndex 2.00000 frozen 36 3 powerlaw norm 5.65000E-04 frozen ________________________________________________________________________ C-statistic = 128.73 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata 17.05_17.10_ratio0.9_shift3.2mA.dat XSPEC12> error 1. 2 Parameter Confidence Range (1) 2 0.705055 2.51654 (-0.463562,1.34793) XSPEC12> error 1. 3 Parameter Confidence Range (1) 3 0.49871 0.699731 (-0.0794314,0.121589) XSPEC12> error 1. 17 Parameter Confidence Range (1) 17 5.67689e-05 6.6086e-05 (-4.5009e-06,4.81626e-06) Looking at Mauche and Liedahl 2001, looks like the ratio could be as low as 0.6; so -- ratio = 0.6: ======================================================================== 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 0.730849 +/- 0.385972 3 1 windprof u0 0.571560 +/- 6.22526E-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 17.0510 frozen 14 1 windprof shift mA 3.20000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 7.24878E-05 +/- 1.09591E-05 18 2 windprof q 0.0 frozen 19 2 windprof taustar 0.730849 = 2 20 2 windprof u0 0.571560 = 3 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 17.0960 frozen 31 2 windprof shift mA 3.20000 frozen 32 2 windprof velocity km/s 2450.00 frozen 33 2 windprof verbose 0 frozen 34 2 windprof norm 4.34927E-05 = 0.6*17 35 3 powerlaw PhoIndex 2.00000 frozen 36 3 powerlaw norm 5.65000E-04 frozen ________________________________________________________________________ C-statistic = 131.21 using 134 PHA bins and 131 degrees of freedom. XSPEC12>iplot PLT> wdata 17.05_17.10_ratio0.6_shift3.2mA.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.405318 1.19672 (-0.324426,0.466981) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.508645 0.637621 (-0.0632015,0.0657743) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 6.70791e-05 7.81085e-05 (-5.40635e-06,5.62298e-06) Ne X Ly-beta @ 10.239 (fit with two windprofs, just because that's what was defined; also, fit continuum simultaneously) -- Emma did not do this line (10.12:10.34) ======================================================================== 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 0.418726 +/- 0.650489 3 1 windprof u0 0.677850 +/- 0.127723 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 10.2390 frozen 14 1 windprof shift mA 1.90000 frozen 15 1 windprof velocity km/s 2450.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.51813E-06 +/- 1.03476E-06 18 2 windprof q 0.0 frozen 19 2 windprof taustar 0.418726 = 2 20 2 windprof u0 0.677850 = 3 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 10.2400 frozen 31 2 windprof shift mA 1.90000 = 14 32 2 windprof velocity km/s 2450.00 frozen 33 2 windprof verbose 0 frozen 34 2 windprof norm 1.51088E-06 = 0.6*17 35 3 powerlaw PhoIndex 2.00000 frozen 36 3 powerlaw norm 4.38945E-04 frozen ________________________________________________________________________ C-statistic = 111.72 using 86 PHA bins and 83 degrees of freedom. XSPEC12>iplot PLT> wdata 10.239_shift1.9mA.dat steppar on taustar gives 0.03 to 1.26 XSPEC12> error 1. 3 Parameter Confidence Range (1) 3 0.586667 0.770475 (-0.0913557,0.0924522) XSPEC12> error 1. 17 Parameter Confidence Range (1) 17 1.9391e-06 2.97855e-06 (-5.79029e-07,4.60419e-07)