3 aug 11 back to 15.014 A line fits using the same proceedure as we did for the RGS data (non-porous, then both kinds of porosity with hinf = 0.5, 1, 5) fit line on 14.87 to 15.13 and use old continuum determination: n=2 norm = 1.71e-3 @loadall.xcm @setup_12.6.xcm non-porous ======================================================================== 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.94217 +/- 0.368274 3 1 windprof u0 0.645932 +/- 6.04784E-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. Warning: Cstat statistic is only valid for Poisson data. XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 1.6117 2.25982 (-0.330472,0.317645) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.597128 0.700121 (-0.0488044,0.0541883) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000506852 0.000548221 (-1.76173e-05,2.37517e-05) XSPEC12>iplot PLT> wdata 15.014_nonporous.dat anisoporous cases aniso, hinf free: hinf = 0 is favored ======================================================================== 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.87986 +/- 0.786075 3 1 windprof u0 0.640037 +/- 9.61391E-02 4 1 windprof h 4.93082E-05 +/- 0.379279 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 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.27144E-04 +/- 3.89578E-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 304 degrees of freedom. aniso, hinf = 0.5 ======================================================================== 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.50563 +/- 0.760933 3 1 windprof u0 0.622560 +/- 6.63523E-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 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 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.22639E-04 +/- 3.89669E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.71000E-03 frozen ________________________________________________________________________ C-statistic = 293.94 using 308 PHA bins and 305 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_aniso_h0.5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 2.10965 3.18653 (-0.648897,0.427976) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.591524 0.681808 (-0.0313756,0.0589092) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000505294 0.000546459 (-1.8099e-05,2.30661e-05) aniso, hinf = 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 3.01979 +/- 1.03882 3 1 windprof u0 0.604510 +/- 5.10852E-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 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 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.22397E-04 +/- 3.88640E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.71000E-03 frozen ________________________________________________________________________ C-statistic = 303.80 using 308 PHA bins and 305 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 2.3722 3.86465 (-0.638393,0.854062) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.571808 0.646824 (-0.0326859,0.0423294) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000504649 0.000545609 (-1.77277e-05,2.32321e-05) XSPEC12>iplot PLT> wdata 15.014_aniso_h1.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 6.94714 +/- 4.38872 3 1 windprof u0 0.543211 +/- 3.36984E-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 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.20579E-04 +/- 3.87098E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.71000E-03 frozen ________________________________________________________________________ C-statistic = 338.35 using 308 PHA bins and 305 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_aniso_h5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 3.94133 12.2664 (-3.00582,5.3193) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.513681 0.579076 (-0.0295304,0.0358648) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000502691 0.000543555 (-1.78886e-05,2.29752e-05) isotropic porosity hinf free goes to zero iso, hinf = 0.5 ======================================================================== 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.18348 +/- 0.462550 3 1 windprof u0 0.632123 +/- 4.89787E-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 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.24264E-04 +/- 3.90780E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.71000E-03 frozen ________________________________________________________________________ C-statistic = 282.01 using 308 PHA bins and 305 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_iso_h0.5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 1.79317 2.60444 (-0.390305,0.420964) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.591942 0.675869 (-0.0401803,0.043746) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000506716 0.000548067 (-1.75478e-05,2.38033e-05) iso, hinf = 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.45522 +/- 0.577492 3 1 windprof u0 0.622027 +/- 4.29247E-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 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.24058E-04 +/- 3.90712E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.71000E-03 frozen ________________________________________________________________________ C-statistic = 283.35 using 308 PHA bins and 305 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_iso_h1.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 1.99265 3.00164 (-0.462425,0.546564) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.587121 0.661021 (-0.0349027,0.0389965) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000506553 0.000547901 (-1.75032e-05,2.38448e-05) iso, hinf = 5 ======================================================================== 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.69573 +/- 2.51430 3 1 windprof u0 0.575602 +/- 2.95293E-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 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.22518E-04 +/- 3.89786E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.71000E-03 frozen ________________________________________________________________________ C-statistic = 297.17 using 308 PHA bins and 305 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_iso_h5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 4.03798 8.45024 (-1.66117,2.75108) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.551821 0.606619 (-0.0237914,0.0310072) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000505042 0.000546349 (-1.74809e-05,2.38262e-05) confidence limits plot XSPEC12>log 15.014_iso_hinf_vs_taustar_35by35.log Logging to file:15.014_iso_hinf_vs_taustar_35by35.log XSPEC12>steppar 4 0. 5. 35 2 0. 11. 35 10may12 ADDING THE H_INF = 2 CASE ANISO: ======================================================================== 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.76003 +/- 1.65786 3 1 windprof u0 0.573359 +/- 4.13054E-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 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 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.22001E-04 +/- 3.87939E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.71000E-03 frozen ________________________________________________________________________ C-statistic = 317.39 using 308 PHA bins and 305 degrees of freedom. XSPEC12>iplot PLT> wdata 15.014_aniso_h2.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 2.75984 5.3211 (-1.00777,1.55349) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.545497 0.613327 (-0.0278996,0.0399309) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000503716 0.000544657 (-1.78437e-05,2.30973e-05) ISO: ======================================================================== 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.09013 +/- 0.881980 3 1 windprof u0 0.606795 +/- 3.65655E-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 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.23633E-04 +/- 3.90528E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.71000E-03 frozen ________________________________________________________________________ C-statistic = 286.44 using 308 PHA bins and 305 degrees of freedom. XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 2.43616 3.98325 (-0.653967,0.89312) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.577673 0.642133 (-0.0291215,0.0353383) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000506171 0.000547504 (-1.74622e-05,2.38709e-05) XSPEC12>iplot PLT> wdata 15.014_iso_h2.dat 15may12 JUST TO BE ON THE SAFE SIDE; REMAKING THE STEPPAR LOG FILE XSPEC12>plot XSPEC12>log fexvii_15.014_both_iso_hnuvstaustar_35by35.log Logging to file:fexvii_15.014_both_iso_hnuvstaustar_35by35.log XSPEC12>steppar 4 0. 5. 35 2 0. 10. 35