7May12 DOCUMENTING POROUS MODEL FITS TO THE NE X 12.134 LINE - reproduced documented non-porous fit that's on the website: taustar = 1.81 uo = 0.623 norm = 3.10e-4 cont norm = 1.45e-3 chisq = 24.04 for N = 19 RGS2 only being used ANISOTROPIC POROSITY (8, 9, 10 = 1) - hinf = 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 1.78019 +/- 2.33933 3 1 windprof u0 0.619498 +/- 0.266016 4 1 windprof h 2.84860E-07 +/- 1.33530 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 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 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 3.09848E-04 +/- 1.75174E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.44552E-03 +/- 6.44261E-04 ________________________________________________________________________ Chi-Squared = 24.06 using 19 PHA bins. Reduced chi-squared = 1.718 for 14 degrees of freedom Null hypothesis probability = 4.509565e-02 XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 1.59108 3.3688 (-0.190831,1.58689) XSPEC12>error 1. 3 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0.99 3 0.565906 0 (-0.0539869,-0.619893) XSPEC12>error 1. 4 Parameter Confidence Range (1) 4 0 0.283233 (-4.13918e-10,0.283233) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000301775 0.000318698 (-8.11497e-06,8.80816e-06) XSPEC12>iplot PLT> wdata nex_12.134_rgs1_aniso_hfree.dat - hinf = 0.5 (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 3.92218 +/- 1.64630 3 1 windprof u0 0.743153 +/- 0.273072 4 1 windprof h 0.500000 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 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 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 3.17742E-04 +/- 1.61921E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.15220E-03 +/- 5.98554E-04 ________________________________________________________________________ Chi-Squared = 27.27 using 19 PHA bins. Reduced chi-squared = 1.818 for 15 degrees of freedom Null hypothesis probability = 2.661977e-02 XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 3.24863 4.74736 (-0.669308,0.829419) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.641609 0.865415 (-0.10146,0.122346) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000309693 0.000325869 (-8.03103e-06,8.14516e-06) XSPEC12>iplot PLT> wdata nex_12.134_rgs1_aniso_h0.5.dat - h_inf = 1 (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 5.99254 +/- 2.75273 3 1 windprof u0 0.706140 +/- 0.190870 4 1 windprof h 1.00000 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 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 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 3.23278E-04 +/- 1.69431E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 9.69536E-04 +/- 6.41047E-04 ________________________________________________________________________ Chi-Squared = 34.22 using 19 PHA bins. Reduced chi-squared = 2.282 for 15 degrees of freedom Null hypothesis probability = 3.167593e-03 XSPEC12>iplot PLT> wdata nex_12.134_rgs1_aniso_h1.dat XSPEC12>error 1. 2 Cannot do error calc: Reduced Chi^2 (= 2.28157) > maximum (2) - h_inf = 2 (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 12.6144 +/- 5.88192 3 1 windprof u0 0.686585 +/- 0.139696 4 1 windprof h 2.00000 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 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 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 3.30706E-04 +/- 1.76983E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 7.13631E-04 +/- 6.87138E-04 ________________________________________________________________________ Chi-Squared = 47.21 using 19 PHA bins. Reduced chi-squared = 3.147 for 15 degrees of freedom Null hypothesis probability = 3.410039e-05 XSPEC12>iplot PLT> wdata nex_12.134_rgs1_aniso_h2.dat - h_inf = 5 (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 64.8152 +/- 21.0016 3 1 windprof u0 0.743409 +/- 9.94079E-02 4 1 windprof h 5.00000 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 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 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 3.43868E-04 +/- 1.77705E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.07889E-04 +/- 7.15626E-04 ________________________________________________________________________ Chi-Squared = 66.88 using 19 PHA bins. Reduced chi-squared = 4.459 for 15 degrees of freedom Null hypothesis probability = 1.597345e-08 XSPEC12>iplot PLT> wdata nex_12.134_rgs1_aniso_h5.dat - ISOTROPIC POROSITY - h_inf 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 1.97331 +/- 3.43696 3 1 windprof u0 0.621274 +/- 6.79586E-02 4 1 windprof h 0.327596 +/- 6.22583 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 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 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 3.09989E-04 +/- 1.65889E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.43463E-03 +/- 6.20294E-04 ________________________________________________________________________ Chi-Squared = 24.04 using 19 PHA bins. Reduced chi-squared = 1.717 for 14 degrees of freedom Null hypothesis probability = 4.532718e-02 XSPEC12>error 1. 2 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 5.46792, 5.51593 and delta stat 0.585968, 3.57835 but latest trial 5.47822 gives 0.518331 Suggest that you check this result using the steppar command. 2 1.61887 5.49193 (-0.354521,3.51853) XSPEC12>steppar 2 1.0 6.0 40 Chi-Squared Delta taustar Chi-Squared 2 47.926 23.887 0 1 37.83 13.791 1 1.125 32.325 8.2862 2 1.25 28.627 4.5878 3 1.375 26.17 2.1309 4 1.5 24.74 0.70065 5 1.625 24.132 0.093099 6 1.75 24.201 0.16157 7 1.875 24.04 0.00047709 8 2 24.042 0.0033998 9 2.125 24.049 0.0096047 10 2.25 24.058 0.01873 11 2.375 24.069 0.030017 12 2.5 24.082 0.043042 13 2.625 24.097 0.057516 14 2.75 24.112 0.07322 15 2.875 24.129 0.089975 16 3 24.147 0.10763 17 3.125 24.165 0.12608 18 3.25 24.184 0.14519 19 3.375 24.204 0.16491 20 3.5 24.224 0.18514 21 3.625 24.245 0.20582 22 3.75 24.266 0.2269 23 3.875 24.287 0.24834 24 4 24.309 0.27009 25 4.125 24.331 0.29211 26 4.25 24.353 0.31438 27 4.375 24.376 0.33688 28 4.5 24.399 0.35956 29 4.625 24.421 0.38242 30 4.75 24.444 0.40543 31 4.875 24.468 0.42858 32 5 24.491 0.45185 33 5.125 24.514 0.47524 34 5.25 24.538 0.49871 35 5.375 24.561 0.52228 36 5.5 24.585 0.54592 37 5.625 24.609 0.56963 38 5.75 24.632 0.5934 39 5.875 24.656 0.61723 40 6 24.688 0.64454 1 6.15 24.715 0.67163 2 6.3 24.743 0.69961 3 6.45 24.772 0.72801 4 6.6 24.8 0.75663 5 6.75 24.829 0.78537 6 6.9 24.858 0.81418 7 7.05 24.887 0.84302 8 7.2 24.916 0.8719 9 7.35 24.945 0.90078 10 7.5 24.973 0.92968 11 7.65 25.002 0.95857 12 7.8 25.031 0.98747 13 7.95 25.06 1.0164 14 8.1 25.089 1.0452 15 8.25 25.118 1.0741 16 8.4 25.147 1.103 17 8.55 25.176 1.1318 18 8.7 25.204 1.1606 19 8.85 25.233 1.1895 20 9 XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.563593 0.707378 (-0.0576801,0.0861054) XSPEC12>error 1. 4 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 4 0 3.83175 (-0.327714,3.50404) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000301702 0.000318436 (-8.27048e-06,8.46279e-06) XSPEC12>iplot PLT> wdata nex_12.134_rgs1_iso_hfree.dat **JUST REALIZED I'VE BEEN SAVING FILES WITH 'RGS1' IN THEIR NAMES; WILL CHANGE THEM ALL TO RGS2 now will run a steppar to make a taustar-hinf contour plot XSPEC12>log nex_12.134_rgs2_hinf_vs_taustar_35by35.log Logging to file:nex_12.134_rgs2_hinf_vs_taustar_35by35.log XSPEC12>steppar 4 0. 5. 35 2 0. 10. 35 NOW BACK TO FIXED H_INF MODELS - h_inf = 0.5 (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 2.07239 +/- 0.299528 3 1 windprof u0 0.620934 +/- 6.31493E-02 4 1 windprof h 0.500000 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 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 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 3.10048E-04 +/- 1.62455E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.42916E-03 +/- 5.35087E-04 ________________________________________________________________________ Chi-Squared = 24.04 using 19 PHA bins. Reduced chi-squared = 1.603 for 15 degrees of freedom Null hypothesis probability = 6.441126e-02 XSPEC12>iplot PLT> wdata nex_12.134_rgs2_iso_h0.5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 1.80356 2.40594 (-0.268836,0.333545) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.568028 0.689458 (-0.0529048,0.0685252) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000301796 0.000318457 (-8.19236e-06,8.46858e-06) - h_inf = 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 2.41085 +/- 0.395580 3 1 windprof u0 0.620970 +/- 6.01925E-02 4 1 windprof h 1.00000 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 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 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 3.10079E-04 +/- 1.61712E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.41779E-03 +/- 5.34173E-04 ________________________________________________________________________ Chi-Squared = 24.06 using 19 PHA bins. Reduced chi-squared = 1.604 for 15 degrees of freedom Null hypothesis probability = 6.407281e-02 XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 2.06761 2.86081 (-0.343233,0.449968) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.56942 0.683616 (-0.0515521,0.0626447) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000301977 0.000318571 (-8.2174e-06,8.37703e-06) XSPEC12>iplot PLT> wdata nex_12.134_rgs2_iso_h1.dat - h_inf = 2 (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 3.41630 +/- 0.754905 3 1 windprof u0 0.624094 +/- 5.77789E-02 4 1 windprof h 2.00000 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 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 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 3.10770E-04 +/- 1.59838E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.37141E-03 +/- 5.31756E-04 ________________________________________________________________________ Chi-Squared = 24.19 using 19 PHA bins. Reduced chi-squared = 1.613 for 15 degrees of freedom Null hypothesis probability = 6.188545e-02 XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 2.81059 4.28968 (-0.605713,0.873379) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.573929 0.682122 (-0.0501634,0.05803) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000302604 0.000318997 (-8.06583e-06,8.32711e-06) XSPEC12>iplot PLT> wdata nex_12.134_rgs2_iso_h2.dat - h_inf = 5 (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 12.9894 +/- 6.50068 3 1 windprof u0 0.633055 +/- 6.32182E-02 4 1 windprof h 5.00000 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 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 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 3.16185E-04 +/- 1.54224E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.12011E-03 +/- 5.56196E-04 ________________________________________________________________________ Chi-Squared = 26.31 using 19 PHA bins. Reduced chi-squared = 1.754 for 15 degrees of freedom Null hypothesis probability = 3.489595e-02 XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 8.8685 19.6101 (-4.05066,6.69096) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.5836 0.686875 (-0.0491557,0.0541192) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000307811 0.000323646 (-7.67837e-06,8.15672e-06) XSPEC12>iplot PLT> wdata nex_12.134_rgs2_iso_h5.dat OK, ONE LAST THING - GOING BACK TO THE ANISO CASE AND RUNNING A STEPPAR GRID TO MAKE A CONTOUR PLOT FOR THE CONFIDENCE LIMITS FOR THAT CASE, TOO XSPEC12>log nex_12.134_rgs2_aniso_hinf_vs_taustar_35by35.log Logging to file:nex_12.134_rgs2_aniso_hinf_vs_taustar_35by35.log XSPEC12> steppar 4 0. 5. 35 2 0. 10. 35