D. Cohen 19sep08 Testing the sensitivity of the profile fits to the terminal velocity Fit quality as a function of terminal velocity; can the data themselves help us address the issue of the sensitivity of taustar to vinf? Note that we might consider fitting on a different (wider) range of wavelengths, especially for bigger vinf values. But at first, we'll just use the same ranges. 15.014 - MEG and HEG simultaneously fit with the canonical non-porous models reproduced the fit on the project website: ======================================================================== 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.94220 +/- 0.368279 3 1 windprof u0 0.645934 +/- 6.04798E-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" 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 5.24469E-04 +/- 3.90823E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.71000E-03 frozen ________________________________________________________________________ C-statistic = 280.83 using 308 PHA bins and 305 degrees of freedom. First, try thawing vinf - just treat it like a free parameter ======================================================================== 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.62835 +/- 0.426271 3 1 windprof u0 0.650178 +/- 4.71260E-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" 2351.61 +/- 103.183 16 1 windprof verbose 0 frozen 17 1 windprof norm 5.27640E-04 +/- 3.91732E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.71000E-03 frozen ________________________________________________________________________ C-statistic = 279.63 using 308 PHA bins and 304 degrees of freedom. XSPEC12>goodness 1000 nosim 23.50% of realizations are < best fit statistic 279.63 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/fexvii_15014_smooth_vinf_free.dat XSPEC12>error 1. 15 Parameter Confidence Range (1.0000000) 15 2258.895185 2453.247928 (-92.679873,101.672869) XSPEC12>error 2.7 15 Parameter Confidence Range (2.700000) 15 2203.121209 2530.949260 (-148.464227,179.363823) XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) 2 1.275298 2.023774 (-0.353256,0.395221) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 1.089345 2.302623 (-0.539205,0.674074) XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) 3 0.608937 0.691728 (-0.041239,0.041553) XSPEC12>error 2.7 3 Parameter Confidence Range (2.700000) 3 0.584584 0.732412 (-0.065592,0.082237) XSPEC12>error 1. 17 Parameter Confidence Range (1.000000) 17 0.000509 0.000550 (-0.000019,0.000022) XSPEC12>error 2.7 17 Parameter Confidence Range (2.700000) 17 0.000496 0.000564 (-0.000032,0.000036) Now, let's see how the taustar-uo confidence limits/regions change with vinf free XSPEC12>log /home/cohen/xspec/2007/zpup/fexvii_15014_smooth_vinf_free_35by35.log Logging to file:/home/cohen/xspec/2007/zpup/fexvii_15014_smooth_vinf_free_35by35.log XSPEC12>steppar 2 0.01 4.5 35 3 0.48 0.99 35 Note: what if vinf=3000? ======================================================================== 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.719850 +/- 0.146791 3 1 windprof u0 0.709112 +/- 2.08152E-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" 3000.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 5.39999E-04 +/- 3.94356E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.71000E-03 frozen ________________________________________________________________________ C-statistic = 295.86 using 308 PHA bins and 305 degrees of freedom. XSPEC12>goodness 1000 nosim 26.50% of realizations are < best fit statistic 295.86 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/fexvii_15014_smooth_vinf_3000.dat and if vinf=2000? ======================================================================== 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.09594 +/- 0.581824 3 1 windprof u0 0.655139 +/- 0.219323 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" 2000.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 5.14449E-04 +/- 3.89927E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.71000E-03 frozen ________________________________________________________________________ C-statistic = 296.97 using 308 PHA bins and 305 degrees of freedom. XSPEC12> goodness 1000 nosim 74.40% of realizations are < best fit statistic 296.97 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/fexvii_15014_smooth_vinf_2000.dat Next line, O VIII 18.969 ======================================================================== 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.93910 +/- 1.34251 3 1 windprof u0 0.876124 +/- 3.74638 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.9690 frozen 14 1 windprof shift "mA" 0.0 frozen 15 1 windprof velocity "km/s" 2096.17 +/- 123.752 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.68864E-04 +/- 5.20017E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 149.75 using 130 PHA bins and 126 degrees of freedom. XSPEC12>goodness 1000 nosim 86.20% of realizations are < best fit statistic 149.75 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/oviii_18969_smooth_vinf_free.dat XSPEC12>error 1. 15 Parameter Confidence Range (1.0000000) 15 1994.581020 2240.226408 (-101.811121,143.834267) XSPEC12>steppar 15 1700. 2400. 35 C-Statistic Delta velocity C-Statistic 15 160.64324 10.79429 0 1700.00000 159.33760 9.48865 1 1720.00000 158.17190 8.32295 2 1740.00000 157.13467 7.28572 3 1760.00000 156.50740 6.65845 4 1780.00000 156.22848 6.37953 5 1800.00000 155.86114 6.01219 6 1820.00000 155.33968 5.49073 7 1840.00000 154.77994 4.93099 8 1860.00000 154.21733 4.36838 9 1880.00000 153.42294 3.57399 10 1900.00000 152.82295 2.97400 11 1920.00000 152.22100 2.37205 12 1940.00000 151.60994 1.76099 13 1960.00000 151.06221 1.21326 14 1980.00000 150.70474 0.85579 15 2000.00000 150.64005 0.79110 16 2020.00000 150.10682 0.25787 17 2040.00000 149.88433 0.03538 18 2060.00000 149.83226 -0.01669 19 2080.00000 149.79089 -0.05806 20 2100.00000 149.96281 0.11386 21 2120.00000 149.85472 0.00577 22 2140.00000 149.87754 0.02859 23 2160.00000 150.05018 0.20123 24 2180.00000 150.25990 0.41095 25 2200.00000 150.46552 0.61657 26 2220.00000 150.74259 0.89364 27 2240.00000 151.21931 1.37036 28 2260.00000 151.37858 1.52963 29 2280.00000 151.73964 1.89069 30 2300.00000 152.12627 2.27732 31 2320.00000 152.69598 2.84703 32 2340.00000 153.05250 3.20355 33 2360.00000 153.50769 3.65873 34 2380.00000 153.98300 4.13405 35 2400.00000 XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) 2 2.826168 5.353310 (-1.096110,1.431032) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 2.262671 7.350014 (-1.628823,3.458520) XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) ***Warning: Parameter pegged at hard limit: 0.99 3 0.591087 0.000000 (-0.314566,-0.905653) XSPEC12>error 2.7 3 Parameter Confidence Range (2.700000) 3 0.409275 0.000000 (-0.471450,-0.880724) XSPEC12>log /home/cohen/xspec/2007/zpup/oviii_18969_smooth_vinf_free_35by35.log Logging to file:/home/cohen/xspec/2007/zpup/oviii_18969_smooth_vinf_free_35by35.log XSPEC12>steppar 2 0.01 4.5 35 3 0.48 0.99 35 and, with vinf=3000 frozen ======================================================================== 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.884169 +/- 0.301588 3 1 windprof u0 0.703619 +/- 4.43016E-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.9690 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 3.64429E-04 +/- 5.36235E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 165.40 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/oviii_18969_smooth_vinf_3000.dat and with vinf=1500. ======================================================================== 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 99.9995 +/- 481.061 3 1 windprof u0 0.100021 +/- 2706.57 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.9690 frozen 14 1 windprof shift "mA" 0.0 frozen 15 1 windprof velocity "km/s" 1500.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.39502E-04 +/- 5.01209E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 187.68 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/oviii_18969_smooth_vinf_1500.dat XSPEC12>goodness 1000 nosim 100.00% of realizations are < best fit statistic 187.68 Next Line, Mg XII at 8.421 A ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 2.89669 +/- 1.15894 3 1 windprof u0 0.990000 +/- 5.65832 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" 0.0 frozen 15 1 windprof velocity "km/s" 1725.50 +/- 200.736 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.82631E-05 +/- 4.64028E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 184.80 using 188 PHA bins and 184 degrees of freedom. XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/mgxii_8421_smooth_vinf_free.dat XSPEC12>error 1. 15 Parameter Confidence Range (1.0000000) Fit Parameter: 3 has pegged Fit Parameter: 3 has pegged Fit Parameter: 3 has pegged 15 1543.783453 1992.161520 (-180.668878,267.709189) Fit Parameter: 3 has pegged Fit Parameter: 3 has pegged Fit Parameter: 3 has pegged XSPEC12>error 2.7 15 Fit Parameter: 3 has pegged ***Warning: Parameter pegged at hard limit: 5000 15 1437.966347 1992.161520 (-286.490438,267.704735) Fit Parameter: 3 has pegged let's try a steppar 193.66339 7.83041 0 1200.00000 192.42146 6.58848 1 1220.00000 192.57999 6.74701 2 1240.00000 191.54450 5.71152 3 1260.00000 191.03969 5.20672 4 1280.00000 190.64077 4.80779 5 1300.00000 190.25450 4.42153 6 1320.00000 189.72625 3.89328 7 1340.00000 189.50769 3.67471 8 1360.00000 189.18424 3.35126 9 1380.00000 188.69093 2.85796 10 1400.00000 188.19064 2.35766 11 1420.00000 187.70636 1.87338 12 1440.00000 187.22739 1.39441 13 1460.00000 186.73665 0.90368 14 1480.00000 186.40517 0.57219 15 1500.00000 186.10428 0.27131 16 1520.00000 185.83949 0.00652 17 1540.00000 185.61224 -0.22074 18 1560.00000 185.48950 -0.34348 19 1580.00000 185.25564 -0.57733 20 1600.00000 185.16559 -0.66739 21 1620.00000 185.08133 -0.75164 22 1640.00000 184.97867 -0.85431 23 1660.00000 184.91077 -0.92221 24 1680.00000 184.85044 -0.98253 25 1700.00000 184.81138 -1.02160 26 1720.00000 184.79447 -1.03851 27 1740.00000 184.80087 -1.03211 28 1760.00000 184.84698 -0.98600 29 1780.00000 184.89531 -0.93767 30 1800.00000 184.97610 -0.85688 31 1820.00000 185.06159 -0.77139 32 1840.00000 185.16717 -0.66581 33 1860.00000 185.26697 -0.56601 34 1880.00000 185.36488 -0.46810 35 1900.00000 185.45787 -0.37511 36 1920.00000 185.54657 -0.28641 37 1940.00000 185.64503 -0.18795 38 1960.00000 185.74653 -0.08644 39 1980.00000 185.84928 0.01630 40 2000.00000 185.94711 0.11413 41 2020.00000 186.03860 0.20563 42 2040.00000 186.12458 0.29160 43 2060.00000 186.20585 0.37287 44 2080.00000 186.28043 0.44745 45 2100.00000 186.34456 0.51158 46 2120.00000 186.42248 0.58950 47 2140.00000 186.46832 0.63534 48 2160.00000 186.50665 0.67367 49 2180.00000 186.53396 0.70099 50 2200.00000 XSPEC12>steppar 15 2200 2600 20 C-Statistic Delta velocity C-Statistic 15 188.36056 2.52759 0 2200.00000 186.54985 0.71687 1 2220.00000 186.55701 0.72403 2 2240.00000 186.55576 0.72278 3 2260.00000 186.54616 0.71318 4 2280.00000 186.52702 0.69404 5 2300.00000 186.49926 0.66629 6 2320.00000 186.48722 0.65424 7 2340.00000 186.44451 0.61154 8 2360.00000 186.39578 0.56280 9 2380.00000 186.34150 0.50852 10 2400.00000 186.28249 0.44952 11 2420.00000 186.21952 0.38654 12 2440.00000 186.15408 0.32111 13 2460.00000 186.08858 0.25560 14 2480.00000 186.04632 0.21334 15 2500.00000 185.98066 0.14768 16 2520.00000 185.91623 0.08326 17 2540.00000 185.85468 0.02170 18 2560.00000 185.79620 -0.03678 19 2580.00000 185.74065 -0.09232 20 2600.00000 OK; it does seem unconstrained at the high end. Other parameters: XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) Fit Parameter: 3 has pegged 2 1.632746 4.134742 (-1.278568,1.223428) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) Apparent non-monotonicity in statistic space detected. Current bracket values 1.134812, 0.000000 and delta stat 1.784537, 10.837802 but latest trial 0.625928 gives 0.496223 Suggest that you check this result using the steppar command. Fit Parameter: 3 has pegged 2 0.567406 6.702712 (-2.343908,3.791398) So, basically unconstrained problems with uo; can try again if we need these values. arbitrarily choosing vinf=2600 frozen, the best-fit model has Delta C ~ 1 ======================================================================== 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.847171 +/- 0.372801 3 1 windprof u0 0.758226 +/- 6.04637E-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" 0.0 frozen 15 1 windprof velocity "km/s" 2600.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.06809E-05 +/- 4.70443E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 185.74 using 188 PHA bins and 185 degrees of freedom. XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/mgxii_8421_smooth_vinf_2600.dat and with vinf=3000 frozen, the fit is a bit better, even: ======================================================================== 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.625810 +/- 0.278733 3 1 windprof u0 0.782553 +/- 4.38057E-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" 0.0 frozen 15 1 windprof velocity "km/s" 3000.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.12171E-05 +/- 4.81255E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 185.28 using 188 PHA bins and 185 degrees of freedom. XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/mgxii_8421_smooth_vinf_3000.dat and at the other extreme, vinf=1500 frozen, the fit is similarly not too bad: ======================================================================== 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.95760 +/- 1.64378 3 1 windprof u0 0.989940 +/- 577.691 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" 0.0 frozen 15 1 windprof velocity "km/s" 1500.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.83198E-05 +/- 4.44563E-06 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.65000E-03 frozen ________________________________________________________________________ C-statistic = 186.49 using 188 PHA bins and 185 degrees of freedom. XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/mgxii_8421_smooth_vinf_1500.dat Fourth line: Ne X Lyman alpha at 12.134 reproduced documented best fit: C=191.35 ======================================================================== 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.68142 +/- 0.375494 3 1 windprof u0 0.683263 +/- 4.77307E-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" 2374.71 +/- 99.4246 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.71793E-04 +/- 1.77022E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.64000E-03 frozen ________________________________________________________________________ C-statistic = 190.91 using 176 PHA bins and 172 degrees of freedom. XSPEC12>error 1. 15 Parameter Confidence Range (1) 15 2207.9 2582.02 (-169.988,204.123) XSPEC12>error 2.7 15 Parameter Confidence Range (2.7) 15 2124.44 2699.02 (-255.758,318.819) XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/nex_12134_smooth_vinf_free.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 1.2001 2.193 (-0.465243,0.527656) XSPEC12>error 2.7 2 Parameter Confidence Range (2.7) 2 1.035 2.58425 (-0.627676,0.921567) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.638661 0.718264 (-0.045137,0.0344656) XSPEC12>error 2.7 3 Parameter Confidence Range (2.7) 3 0.61208 0.774712 (-0.0717887,0.0908436) XSPEC12>log /home/cohen/xspec/2007/zpup/nex_12134_smooth_vinf_free_35by35.log Logging to file:/home/cohen/xspec/2007/zpup/nex_12134_smooth_vinf_free_35by35.log XSPEC12>steppar 2 0.01 4.5 35 3 0.48 0.99 35 best-fit models for a couple of fixed vinf values: vinf=1800: ======================================================================== 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.70953 +/- 0.517103 3 1 windprof u0 0.990000 +/- 108.443 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" 1800.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.64115E-04 +/- 1.69517E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.64000E-03 frozen ________________________________________________________________________ C-statistic = 219.50 using 176 PHA bins and 173 degrees of freedom. XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/nex_12134_smooth_vinf_1800.dat vinf=3000: ======================================================================== 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.868394 +/- 0.137552 3 1 windprof u0 0.730569 +/- 2.20521E-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" 3000.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.80849E-04 +/- 1.78951E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.64000E-03 frozen ________________________________________________________________________ C-statistic = 198.84 using 176 PHA bins and 173 degrees of freedom. XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/nex_12134_smooth_vinf_3000.dat Final line to test - Fe XVII at 16.780A reproduced 2250 fit with C = 102.88 Now, with vinf free ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 2.30832 +/- 0.941620 3 1 windprof u0 0.741533 +/- 0.204506 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" 0.0 frozen 15 1 windprof velocity "km/s" 2372.27 +/- 180.338 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.31390E-04 +/- 3.06108E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.20000E-03 frozen ________________________________________________________________________ C-statistic = 102.67 using 104 PHA bins and 100 degrees of freedom. XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/fexvii_1678_smooth_vinf_free.dat XSPEC12>error 1. 15 Parameter Confidence Range (1) 15 2148.13 2656.79 (-224.348,284.306) XSPEC12>error 2.7 15 Parameter Confidence Range (2.7) ***Warning: Number of trials exceeded before convergence. Current trial values 2068.09, 2068.08 and delta statistic 2.32075, 2.74006 ***Warning: Identical values of the parameter give different values of the statistic. Please check your result for the low end of the confidence range. ***Warning: Identical values of the parameter give different values of the statistic. Please check your result for the low end of the confidence range. 15 2068.08 2981.44 (-304.435,608.915) now, using steppar: XSPEC12>steppar 15 2000. 3000. 50 C-Statistic Delta velocity C-Statistic 15 106.8 4.1094 0 2000 106.17 3.4866 1 2020 105.63 2.9418 2 2040 105.14 2.4523 3 2060 104.71 2.0237 4 2080 104.33 1.6464 5 2100 104.01 1.3204 6 2120 103.73 1.0448 7 2140 103.5 0.81277 8 2160 103.3 0.61312 9 2180 103.13 0.4469 10 2200 103.01 0.32677 11 2220 102.92 0.23022 12 2240 102.85 0.16082 13 2260 102.8 0.11482 14 2280 102.74 0.056354 15 2300 102.7 0.010618 16 2320 102.66 -0.025556 17 2340 102.64 -0.047604 18 2360 102.63 -0.055109 19 2380 102.64 -0.049199 20 2400 102.66 -0.026546 21 2420 102.7 0.012263 22 2440 102.75 0.065957 23 2460 102.82 0.13287 24 2480 102.9 0.20958 25 2500 102.98 0.29875 26 2520 103.08 0.39699 27 2540 103.19 0.50272 28 2560 103.3 0.61463 29 2580 103.42 0.73147 30 2600 103.54 0.8512 31 2620 103.66 0.97273 32 2640 103.78 1.0956 33 2660 103.91 1.219 34 2680 104.03 1.342 35 2700 104.15 1.4634 36 2720 104.27 1.583 37 2740 104.39 1.7004 38 2760 104.5 1.8153 39 2780 104.61 1.9271 40 2800 104.72 2.0358 41 2820 104.83 2.1412 42 2840 104.93 2.2434 43 2860 105.03 2.3421 44 2880 105.12 2.4374 45 2900 105.22 2.5294 46 2920 105.3 2.6179 47 2940 105.39 2.7032 48 2960 105.47 2.7852 49 2980 105.55 2.864 50 3000 so, 2050 to 2940 is roughly the 90% confidence range confidence limits on other parameters: XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 1.51006 3.30734 (-0.797545,0.999742) XSPEC12>error 2.7 2 Parameter Confidence Range (2.7) 2 1.15246 3.88516 (-1.15509,1.57761) XSPEC12>error 1. 3 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0.99 3 0.620976 0 (-0.120514,-0.74149) XSPEC12>error 2.7 3 Parameter Confidence Range (2.7) ***Warning: Parameter pegged at hard limit: 0.99 3 0.56401 0 (-0.17748,-0.74149) and now, saving best-fit models with specific, fixed vinf values: vinf=1800: ======================================================================== 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 6.23048 +/- 1.50415 3 1 windprof u0 0.989917 +/- 5565.63 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" 0.0 frozen 15 1 windprof velocity "km/s" 1800.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.21653E-04 +/- 2.96265E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.20000E-03 frozen ________________________________________________________________________ C-statistic = 115.62 using 104 PHA bins and 101 degrees of freedom. XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/fexvii_1678_smooth_vinf_1800.dat ======================================================================== 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.21949 +/- 0.354316 3 1 windprof u0 0.746090 +/- 6.05201E-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" 0.0 frozen 15 1 windprof velocity "km/s" 3000.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 2.37641E-04 +/- 3.13814E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.20000E-03 frozen ________________________________________________________________________ C-statistic = 105.61 using 104 PHA bins and 101 degrees of freedom. XSPEC12>iplot PLT> wdata /home/cohen/xspec/2007/zpup/fexvii_1678_smooth_vinf_3000.dat XSPEC12>log /home/cohen/xspec/2007/zpup/fexvii_1678_smooth_vinf_free_35by35.log Logging to file:/home/cohen/xspec/2007/zpup/fexvii_1678_smooth_vinf_free_35by35.log XSPEC12>steppar 2 0.01 4.5 35 3 0.48 0.99 35