!XSPEC12>data 1:1 xrt_wt_nopu_src_grp.pha;
#
#1 spectrum in use
#
#Spectral Data File: xrt_wt_nopu_src_grp.pha Spectrum 1
#Net count rate (cts/s) for Spectrum:1 2.034e+01 +/- 2.907e-01
# Assigned to Data Group 1 and Plot Group 1
# Noticed Channels: 1-457
# Telescope: SWIFT Instrument: XRT Channel Type: PI
# Exposure Time: 240.7 sec
# No response loaded.
#
#***Warning! One or more spectra are missing responses,
# and are not suitable for fit.
#
!XSPEC12>backgrnd 1 xrt_wt_nopu_back.pha;
#Net count rate (cts/s) for Spectrum:1 2.018e+01 +/- 2.911e-01 (99.2 % total)
#***Warning! One or more spectra are missing responses,
# and are not suitable for fit.
#
!XSPEC12>response 1 /local/data/raid2/caldb/data/swift/xrt/cpf/rmf/swxwt0to2s0_20010101v011.rmf;
#***Warning: Detected response matrix energy bin value = 0 (or neg).
# XSPEC will instead use small finite value (response file will not be altered).
#Response successfully loaded.
#
!XSPEC12>arf 1 xrt_wt_nopu_src.arf;
#Arf successfully loaded.
#
!XSPEC12>setplot energy;
#
!XSPEC12>query no;
#
!XSPEC12>ignore 1:**-0.3 10.0-**;
# 30 channels (1-30) ignored in spectrum # 1
# 25 channels (433-457) ignored in spectrum # 1
#
#
!XSPEC12>tclout rate 1;
#
!XSPEC12>scan $xspec_tclout "%f %f" cr crerr;
#
!XSPEC12>ignore bad;
#
#ignore: 292 channels ignored from source number 1
#
!XSPEC12>model wabs*powerlaw ; 0.1 ; 2.0 ; 1.0;
#
#========================================================================
#Model wabs<1>*powerlaw<2> Source No.: 1 Active/On
#Model Model Component Parameter Unit Value
# par comp
# 1 1 wabs nH 10^22 0.100000 +/- 0.0
# 2 2 powerlaw PhoIndex 2.00000 +/- 0.0
# 3 2 powerlaw norm 1.00000 +/- 0.0
#________________________________________________________________________
#
#
# Chi-Squared = 139915.0 using 165 PHA bins.
# Reduced chi-squared = 863.6727 for 162 degrees of freedom
# Null hypothesis probability = 0.000000e+00
# Current data and model not fit yet.
#
!XSPEC12>tclout dof;
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set dof [lindex $xspec_tclout 0];
#
!XSPEC12>if {$dof <= 0} {; exit;;};
#
!XSPEC12>renorm;
#
# Chi-Squared = 433.23 using 165 PHA bins.
# Reduced chi-squared = 2.6742 for 162 degrees of freedom
# Null hypothesis probability = 1.331241e-26
# Current data and model not fit yet.
#
!XSPEC12>fit;
# renorm: no renormalization necessary
#
#Chi-Squared Lvl Par # 1 2 3
#232.782 -3 0.0495859 2.03258 0.122258
#146.112 -4 0.0427531 2.15227 0.138207
#144.783 -5 0.0487871 2.16837 0.140764
#144.783 -6 0.0486512 2.16809 0.14069
#========================================
# Variances and Principal Axes
# 1 2 3
# 3.3324E-06| -0.5426 0.0023 0.8400
# 1.8555E-05| 0.8286 -0.1625 0.5357
# 1.9741E-03| -0.1377 -0.9867 -0.0863
#----------------------------------------
#
#====================================
# Covariance Matrix
# 1 2 3
#5.114e-05 2.657e-04 3.016e-05
#2.657e-04 1.922e-03 1.664e-04
#3.016e-05 1.664e-04 2.236e-05
#------------------------------------
#
#========================================================================
#Model wabs<1>*powerlaw<2> Source No.: 1 Active/On
#Model Model Component Parameter Unit Value
# par comp
# 1 1 wabs nH 10^22 4.86512E-02 +/- 7.15147E-03
# 2 2 powerlaw PhoIndex 2.16809 +/- 4.38458E-02
# 3 2 powerlaw norm 0.140690 +/- 4.72906E-03
#________________________________________________________________________
#
#
# Chi-Squared = 144.78 using 165 PHA bins.
# Reduced chi-squared = 0.89372 for 162 degrees of freedom
# Null hypothesis probability = 8.302852e-01
#
!XSPEC12>fit;
#
#Chi-Squared Lvl Par # 1 2 3
#144.783 -3 0.0486598 2.16811 0.140695
#========================================
# Variances and Principal Axes
# 1 2 3
# 3.3282E-06| -0.5428 0.0023 0.8399
# 1.8531E-05| 0.8285 -0.1624 0.5359
# 1.9732E-03| -0.1376 -0.9867 -0.0862
#----------------------------------------
#
#====================================
# Covariance Matrix
# 1 2 3
#5.107e-05 2.654e-04 3.011e-05
#2.654e-04 1.922e-03 1.662e-04
#3.011e-05 1.662e-04 2.233e-05
#------------------------------------
#
#========================================================================
#Model wabs<1>*powerlaw<2> Source No.: 1 Active/On
#Model Model Component Parameter Unit Value
# par comp
# 1 1 wabs nH 10^22 4.86598E-02 +/- 7.14601E-03
# 2 2 powerlaw PhoIndex 2.16811 +/- 4.38373E-02
# 3 2 powerlaw norm 0.140695 +/- 4.72517E-03
#________________________________________________________________________
#
#
# Chi-Squared = 144.78 using 165 PHA bins.
# Reduced chi-squared = 0.89372 for 162 degrees of freedom
# Null hypothesis probability = 8.302852e-01
#
!XSPEC12>fit;
#
#Chi-Squared Lvl Par # 1 2 3
#144.783 2 0.0486598 2.16811 0.140695
#========================================
# Variances and Principal Axes
# 1 2 3
# 3.3285E-06| -0.5428 0.0023 0.8399
# 1.8532E-05| 0.8285 -0.1624 0.5359
# 1.9733E-03| -0.1376 -0.9867 -0.0862
#----------------------------------------
#
#====================================
# Covariance Matrix
# 1 2 3
#5.107e-05 2.654e-04 3.012e-05
#2.654e-04 1.922e-03 1.662e-04
#3.012e-05 1.662e-04 2.233e-05
#------------------------------------
#
#========================================================================
#Model wabs<1>*powerlaw<2> Source No.: 1 Active/On
#Model Model Component Parameter Unit Value
# par comp
# 1 1 wabs nH 10^22 4.86598E-02 +/- 7.14635E-03
# 2 2 powerlaw PhoIndex 2.16811 +/- 4.38381E-02
# 3 2 powerlaw norm 0.140695 +/- 4.72543E-03
#________________________________________________________________________
#
#
# Chi-Squared = 144.78 using 165 PHA bins.
# Reduced chi-squared = 0.89372 for 162 degrees of freedom
# Null hypothesis probability = 8.302852e-01
#
!XSPEC12>error maximum 100 1;
# Parameter Confidence Range (2.706)
# 1 0.0372475 0.0607537 (-0.0114123,0.0120939)
#
!XSPEC12>error maximum 100 2;
# Parameter Confidence Range (2.706)
# 2 2.09539 2.2442 (-0.0727247,0.0760822)
#
!XSPEC12>error maximum 100 3;
# Parameter Confidence Range (2.706)
# 3 0.133132 0.148758 (-0.00756311,0.0080634)
#
!XSPEC12>fit;
#
#Chi-Squared Lvl Par # 1 2 3
#144.783 2 0.0486598 2.16811 0.140695
#========================================
# Variances and Principal Axes
# 1 2 3
# 3.3285E-06| -0.5428 0.0023 0.8399
# 1.8532E-05| 0.8285 -0.1624 0.5359
# 1.9733E-03| -0.1376 -0.9867 -0.0862
#----------------------------------------
#
#====================================
# Covariance Matrix
# 1 2 3
#5.107e-05 2.654e-04 3.012e-05
#2.654e-04 1.922e-03 1.662e-04
#3.012e-05 1.662e-04 2.233e-05
#------------------------------------
#
#========================================================================
#Model wabs<1>*powerlaw<2> Source No.: 1 Active/On
#Model Model Component Parameter Unit Value
# par comp
# 1 1 wabs nH 10^22 4.86598E-02 +/- 7.14635E-03
# 2 2 powerlaw PhoIndex 2.16811 +/- 4.38381E-02
# 3 2 powerlaw norm 0.140695 +/- 4.72543E-03
#________________________________________________________________________
#
#
# Chi-Squared = 144.78 using 165 PHA bins.
# Reduced chi-squared = 0.89372 for 162 degrees of freedom
# Null hypothesis probability = 8.302852e-01
#
!XSPEC12>open xrt_wt_nopu_pha.dat w
#
!XSPEC12>set fname [open xrt_wt_nopu_pha.dat w];
#
!XSPEC12>format "%-16s%-11.2e%-11.2e%s" "rate" $cr $crerr "\[counts/s (0.3-10.0 keV)\]"
#
!XSPEC12>puts $fname [format "%-16s%-11.2e%-11.2e%s" "rate" $cr $crerr "\[counts/s (0.3-10.0 keV)\]"];
#
!XSPEC12>puts $fname "";
#
!XSPEC12>tclout model;
#
!XSPEC12>set model $xspec_tclout;
#
!XSPEC12>puts $fname "$model";
#
!XSPEC12>tclout modpar;
#
!XSPEC12>set npar $xspec_tclout;
#
!XSPEC12>for {set i 1} {$i <= $npar} {incr i} {; tclout pinfo $i; set pname [lindex $xspec_tclout 0]; tclout param $i; set param [lindex $xspec_tclout 0]; tclout error $i; scan $xspec_tclout "%f %f" paramlo paramhi; if {$paramlo == 0 ;; $paramhi == 0} {; set errlo 0; set errhi 0; } else {; set errlo [expr $paramlo - $param]; set errhi [expr $paramhi - $param]; }; switch $i {; 1 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"]}; 2 {puts $fname [format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi]}; 3 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[cnts/cm^2/s/keV\]"]}; };};
#
!XSPEC12>set i 1
#
!XSPEC12>tclout pinfo $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set pname [lindex $xspec_tclout 0]
#
!XSPEC12>tclout param $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set param [lindex $xspec_tclout 0]
#
!XSPEC12>tclout error $i
#
!XSPEC12>scan $xspec_tclout "%f %f" paramlo paramhi
#
!XSPEC12>if {$paramlo == 0 ;; $paramhi == 0} {; set errlo 0; set errhi 0; } else {; set errlo [expr $paramlo - $param]; set errhi [expr $paramhi - $param]; }
#
!XSPEC12>expr $paramlo - $param
#
!XSPEC12>set errlo [expr $paramlo - $param]
#
!XSPEC12>expr $paramhi - $param
#
!XSPEC12>set errhi [expr $paramhi - $param]
#
!XSPEC12>switch $i {; 1 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"]}; 2 {puts $fname [format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi]}; 3 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[cnts/cm^2/s/keV\]"]}; }
#
!XSPEC12>format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"
#
!XSPEC12>puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"]
#
!XSPEC12>incr i
#
!XSPEC12>tclout pinfo $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set pname [lindex $xspec_tclout 0]
#
!XSPEC12>tclout param $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set param [lindex $xspec_tclout 0]
#
!XSPEC12>tclout error $i
#
!XSPEC12>scan $xspec_tclout "%f %f" paramlo paramhi
#
!XSPEC12>if {$paramlo == 0 ;; $paramhi == 0} {; set errlo 0; set errhi 0; } else {; set errlo [expr $paramlo - $param]; set errhi [expr $paramhi - $param]; }
#
!XSPEC12>expr $paramlo - $param
#
!XSPEC12>set errlo [expr $paramlo - $param]
#
!XSPEC12>expr $paramhi - $param
#
!XSPEC12>set errhi [expr $paramhi - $param]
#
!XSPEC12>switch $i {; 1 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"]}; 2 {puts $fname [format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi]}; 3 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[cnts/cm^2/s/keV\]"]}; }
#
!XSPEC12>format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi
#
!XSPEC12>puts $fname [format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi]
#
!XSPEC12>incr i
#
!XSPEC12>tclout pinfo $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set pname [lindex $xspec_tclout 0]
#
!XSPEC12>tclout param $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set param [lindex $xspec_tclout 0]
#
!XSPEC12>tclout error $i
#
!XSPEC12>scan $xspec_tclout "%f %f" paramlo paramhi
#
!XSPEC12>if {$paramlo == 0 ;; $paramhi == 0} {; set errlo 0; set errhi 0; } else {; set errlo [expr $paramlo - $param]; set errhi [expr $paramhi - $param]; }
#
!XSPEC12>expr $paramlo - $param
#
!XSPEC12>set errlo [expr $paramlo - $param]
#
!XSPEC12>expr $paramhi - $param
#
!XSPEC12>set errhi [expr $paramhi - $param]
#
!XSPEC12>switch $i {; 1 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"]}; 2 {puts $fname [format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi]}; 3 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[cnts/cm^2/s/keV\]"]}; }
#
!XSPEC12>format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[cnts/cm^2/s/keV\]"
#
!XSPEC12>puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[cnts/cm^2/s/keV\]"]
#
!XSPEC12>incr i
#
!XSPEC12>set norm $param;
#
!XSPEC12>set normerrlo $errlo;
#
!XSPEC12>set normerrhi $errhi;
#
!XSPEC12>tclout stat;
#
!XSPEC12>set chisq $xspec_tclout;
#
!XSPEC12>tclout dof;
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set dof [lindex $xspec_tclout 0];
#
!XSPEC12>tclout stat;
#
!XSPEC12>expr $chisq/$dof
#
!XSPEC12>set rchisq [expr $chisq/$dof];
#
!XSPEC12>format "%-16s%.2f" "chisq" $chisq
#
!XSPEC12>puts $fname [format "%-16s%.2f" "chisq" $chisq];
#
!XSPEC12>format "%-16s%d" "dof" $dof
#
!XSPEC12>puts $fname [format "%-16s%d" "dof" $dof];
#
!XSPEC12>format "%-16s%.2f" "rchisq" $rchisq
#
!XSPEC12>puts $fname [format "%-16s%.2f" "rchisq" $rchisq];
#
!XSPEC12>puts $fname "";
#
!XSPEC12>setplot command scr white;
#
!XSPEC12>setplot command col 2 on 2;
#
!XSPEC12>setplot command label top GRB060526;
#
!XSPEC12>setplot command label file;
#
!XSPEC12>setplot command r y2 -3.99 3.99;
#
!XSPEC12>setplot command time off;
#
!XSPEC12>cpd xrt_wt_nopu_pha.gif/gif;
#
!XSPEC12>plot ldata delchi;
#
!XSPEC12>model wabs*pegpwrlw ; 0.1 ; 2.0 ; 0.3 ; 10.0 ; 1.0;
#
#========================================================================
#Model wabs<1>*pegpwrlw<2> Source No.: 1 Active/On
#Model Model Component Parameter Unit Value
# par comp
# 1 1 wabs nH 10^22 0.100000 +/- 0.0
# 2 2 pegpwrlw PhoIndex 2.00000 +/- 0.0
# 3 2 pegpwrlw eMin keV 0.300000 frozen
# 4 2 pegpwrlw eMax keV 10.0000 frozen
# 5 2 pegpwrlw norm 1.00000 +/- 0.0
#________________________________________________________________________
#
#
# Chi-Squared = 4636.48 using 165 PHA bins.
# Reduced chi-squared = 28.6203 for 162 degrees of freedom
# Null hypothesis probability = 0.000000e+00
# Current data and model not fit yet.
#
!XSPEC12>renorm;
#
# Chi-Squared = 433.23 using 165 PHA bins.
# Reduced chi-squared = 2.6742 for 162 degrees of freedom
# Null hypothesis probability = 1.331365e-26
# Current data and model not fit yet.
#
!XSPEC12>fit;
#
#Chi-Squared Lvl Par # 1 2 5
#237.073 -3 0.0495231 2.03193 671.783
#146.102 -4 0.042724 2.15287 721.306
#144.783 -5 0.0488041 2.16846 731.399
#144.783 -6 0.0486536 2.16811 731.102
#========================================
# Variances and Principal Axes
# 1 2 5
# 7.8857E-06| -0.9941 0.1081 0.0002
# 1.5327E-03| 0.1081 0.9941 -0.0012
# 2.7282E+02| 0.0003 0.0012 1.0000
#----------------------------------------
#
#====================================
# Covariance Matrix
# 1 2 3
#5.116e-05 2.658e-04 8.335e-02
#2.658e-04 1.923e-03 3.337e-01
#8.335e-02 3.337e-01 2.728e+02
#------------------------------------
#
#========================================================================
#Model wabs<1>*pegpwrlw<2> Source No.: 1 Active/On
#Model Model Component Parameter Unit Value
# par comp
# 1 1 wabs nH 10^22 4.86536E-02 +/- 7.15252E-03
# 2 2 pegpwrlw PhoIndex 2.16811 +/- 4.38519E-02
# 3 2 pegpwrlw eMin keV 0.300000 frozen
# 4 2 pegpwrlw eMax keV 10.0000 frozen
# 5 2 pegpwrlw norm 731.102 +/- 16.5174
#________________________________________________________________________
#
#
# Chi-Squared = 144.78 using 165 PHA bins.
# Reduced chi-squared = 0.89372 for 162 degrees of freedom
# Null hypothesis probability = 8.302858e-01
#
!XSPEC12>fit;
#
#Chi-Squared Lvl Par # 1 2 5
#144.782 -3 0.0486612 2.16813 731.119
#========================================
# Variances and Principal Axes
# 1 2 5
# 7.8708E-06| -0.9941 0.1081 0.0002
# 1.5335E-03| 0.1081 0.9941 -0.0012
# 2.7229E+02| 0.0003 0.0012 1.0000
#----------------------------------------
#
#====================================
# Covariance Matrix
# 1 2 3
#5.106e-05 2.654e-04 8.312e-02
#2.654e-04 1.922e-03 3.325e-01
#8.312e-02 3.325e-01 2.723e+02
#------------------------------------
#
#========================================================================
#Model wabs<1>*pegpwrlw<2> Source No.: 1 Active/On
#Model Model Component Parameter Unit Value
# par comp
# 1 1 wabs nH 10^22 4.86612E-02 +/- 7.14594E-03
# 2 2 pegpwrlw PhoIndex 2.16813 +/- 4.38380E-02
# 3 2 pegpwrlw eMin keV 0.300000 frozen
# 4 2 pegpwrlw eMax keV 10.0000 frozen
# 5 2 pegpwrlw norm 731.119 +/- 16.5013
#________________________________________________________________________
#
#
# Chi-Squared = 144.78 using 165 PHA bins.
# Reduced chi-squared = 0.89372 for 162 degrees of freedom
# Null hypothesis probability = 8.302896e-01
#
!XSPEC12>fit;
#
#Chi-Squared Lvl Par # 1 2 5
#144.782 -3 0.0486584 2.16811 731.115
#========================================
# Variances and Principal Axes
# 1 2 5
# 7.8715E-06| -0.9941 0.1081 0.0002
# 1.5337E-03| 0.1081 0.9941 -0.0012
# 2.7232E+02| 0.0003 0.0012 1.0000
#----------------------------------------
#
#====================================
# Covariance Matrix
# 1 2 3
#5.107e-05 2.655e-04 8.313e-02
#2.655e-04 1.922e-03 3.326e-01
#8.313e-02 3.326e-01 2.723e+02
#------------------------------------
#
#========================================================================
#Model wabs<1>*pegpwrlw<2> Source No.: 1 Active/On
#Model Model Component Parameter Unit Value
# par comp
# 1 1 wabs nH 10^22 4.86584E-02 +/- 7.14651E-03
# 2 2 pegpwrlw PhoIndex 2.16811 +/- 4.38417E-02
# 3 2 pegpwrlw eMin keV 0.300000 frozen
# 4 2 pegpwrlw eMax keV 10.0000 frozen
# 5 2 pegpwrlw norm 731.115 +/- 16.5021
#________________________________________________________________________
#
#
# Chi-Squared = 144.78 using 165 PHA bins.
# Reduced chi-squared = 0.89372 for 162 degrees of freedom
# Null hypothesis probability = 8.302898e-01
#
!XSPEC12>error maximum 100 1;
# Parameter Confidence Range (2.706)
# 1 0.0372482 0.0607524 (-0.0114102,0.012094)
#
!XSPEC12>error maximum 100 2;
# Parameter Confidence Range (2.706)
# 2 2.09539 2.24419 (-0.0727128,0.0760821)
#
!XSPEC12>error maximum 100 5;
# Parameter Confidence Range (2.706)
# 5 704.997 758.83 (-26.1185,27.7152)
#
!XSPEC12>fit;
#
#Chi-Squared Lvl Par # 1 2 5
#144.782 1 0.0486584 2.16811 731.115
#========================================
# Variances and Principal Axes
# 1 2 5
# 7.8714E-06| -0.9941 0.1081 0.0002
# 1.5337E-03| 0.1081 0.9941 -0.0012
# 2.7227E+02| 0.0003 0.0012 1.0000
#----------------------------------------
#
#====================================
# Covariance Matrix
# 1 2 3
#5.107e-05 2.654e-04 8.311e-02
#2.654e-04 1.922e-03 3.325e-01
#8.311e-02 3.325e-01 2.723e+02
#------------------------------------
#
#========================================================================
#Model wabs<1>*pegpwrlw<2> Source No.: 1 Active/On
#Model Model Component Parameter Unit Value
# par comp
# 1 1 wabs nH 10^22 4.86584E-02 +/- 7.14598E-03
# 2 2 pegpwrlw PhoIndex 2.16811 +/- 4.38393E-02
# 3 2 pegpwrlw eMin keV 0.300000 frozen
# 4 2 pegpwrlw eMax keV 10.0000 frozen
# 5 2 pegpwrlw norm 731.115 +/- 16.5007
#________________________________________________________________________
#
#
# Chi-Squared = 144.78 using 165 PHA bins.
# Reduced chi-squared = 0.89372 for 162 degrees of freedom
# Null hypothesis probability = 8.302898e-01
#
!XSPEC12>tclout model;
#
!XSPEC12>set model $xspec_tclout;
#
!XSPEC12>puts $fname "$model";
#
!XSPEC12>for {set i 1} {$i <= 5} {incr i} {; tclout pinfo $i; set pname [lindex $xspec_tclout 0]; tclout param $i; set param [lindex $xspec_tclout 0]; tclout error $i; scan $xspec_tclout "%f %f" paramlo paramhi; if {$paramlo == 0 ;; $paramhi == 0} {; set errlo 0; set errhi 0; } else {; set errlo [expr $paramlo - $param]; set errhi [expr $paramhi - $param]; }; switch $i {; 1 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"]}; 2 {puts $fname [format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi]}; 3 - 4 {puts $fname [format "%-16s%-11.1f%s" $pname $param "\[keV\]"]}; 5 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^-12 ergs/cm^2/s\]"]}; };};
#
!XSPEC12>set i 1
#
!XSPEC12>tclout pinfo $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set pname [lindex $xspec_tclout 0]
#
!XSPEC12>tclout param $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set param [lindex $xspec_tclout 0]
#
!XSPEC12>tclout error $i
#
!XSPEC12>scan $xspec_tclout "%f %f" paramlo paramhi
#
!XSPEC12>if {$paramlo == 0 ;; $paramhi == 0} {; set errlo 0; set errhi 0; } else {; set errlo [expr $paramlo - $param]; set errhi [expr $paramhi - $param]; }
#
!XSPEC12>expr $paramlo - $param
#
!XSPEC12>set errlo [expr $paramlo - $param]
#
!XSPEC12>expr $paramhi - $param
#
!XSPEC12>set errhi [expr $paramhi - $param]
#
!XSPEC12>switch $i {; 1 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"]}; 2 {puts $fname [format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi]}; 3 - 4 {puts $fname [format "%-16s%-11.1f%s" $pname $param "\[keV\]"]}; 5 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^-12 ergs/cm^2/s\]"]}; }
#
!XSPEC12>format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"
#
!XSPEC12>puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"]
#
!XSPEC12>incr i
#
!XSPEC12>tclout pinfo $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set pname [lindex $xspec_tclout 0]
#
!XSPEC12>tclout param $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set param [lindex $xspec_tclout 0]
#
!XSPEC12>tclout error $i
#
!XSPEC12>scan $xspec_tclout "%f %f" paramlo paramhi
#
!XSPEC12>if {$paramlo == 0 ;; $paramhi == 0} {; set errlo 0; set errhi 0; } else {; set errlo [expr $paramlo - $param]; set errhi [expr $paramhi - $param]; }
#
!XSPEC12>expr $paramlo - $param
#
!XSPEC12>set errlo [expr $paramlo - $param]
#
!XSPEC12>expr $paramhi - $param
#
!XSPEC12>set errhi [expr $paramhi - $param]
#
!XSPEC12>switch $i {; 1 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"]}; 2 {puts $fname [format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi]}; 3 - 4 {puts $fname [format "%-16s%-11.1f%s" $pname $param "\[keV\]"]}; 5 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^-12 ergs/cm^2/s\]"]}; }
#
!XSPEC12>format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi
#
!XSPEC12>puts $fname [format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi]
#
!XSPEC12>incr i
#
!XSPEC12>tclout pinfo $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set pname [lindex $xspec_tclout 0]
#
!XSPEC12>tclout param $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set param [lindex $xspec_tclout 0]
#
!XSPEC12>tclout error $i
#
!XSPEC12>scan $xspec_tclout "%f %f" paramlo paramhi
#
!XSPEC12>if {$paramlo == 0 ;; $paramhi == 0} {; set errlo 0; set errhi 0; } else {; set errlo [expr $paramlo - $param]; set errhi [expr $paramhi - $param]; }
#
!XSPEC12>set errlo 0
#
!XSPEC12>set errhi 0
#
!XSPEC12>switch $i {; 1 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"]}; 2 {puts $fname [format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi]}; 3 - 4 {puts $fname [format "%-16s%-11.1f%s" $pname $param "\[keV\]"]}; 5 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^-12 ergs/cm^2/s\]"]}; }
#
!XSPEC12>format "%-16s%-11.1f%s" $pname $param "\[keV\]"
#
!XSPEC12>puts $fname [format "%-16s%-11.1f%s" $pname $param "\[keV\]"]
#
!XSPEC12>incr i
#
!XSPEC12>tclout pinfo $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set pname [lindex $xspec_tclout 0]
#
!XSPEC12>tclout param $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set param [lindex $xspec_tclout 0]
#
!XSPEC12>tclout error $i
#
!XSPEC12>scan $xspec_tclout "%f %f" paramlo paramhi
#
!XSPEC12>if {$paramlo == 0 ;; $paramhi == 0} {; set errlo 0; set errhi 0; } else {; set errlo [expr $paramlo - $param]; set errhi [expr $paramhi - $param]; }
#
!XSPEC12>set errlo 0
#
!XSPEC12>set errhi 0
#
!XSPEC12>switch $i {; 1 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"]}; 2 {puts $fname [format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi]}; 3 - 4 {puts $fname [format "%-16s%-11.1f%s" $pname $param "\[keV\]"]}; 5 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^-12 ergs/cm^2/s\]"]}; }
#
!XSPEC12>format "%-16s%-11.1f%s" $pname $param "\[keV\]"
#
!XSPEC12>puts $fname [format "%-16s%-11.1f%s" $pname $param "\[keV\]"]
#
!XSPEC12>incr i
#
!XSPEC12>tclout pinfo $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set pname [lindex $xspec_tclout 0]
#
!XSPEC12>tclout param $i
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set param [lindex $xspec_tclout 0]
#
!XSPEC12>tclout error $i
#
!XSPEC12>scan $xspec_tclout "%f %f" paramlo paramhi
#
!XSPEC12>if {$paramlo == 0 ;; $paramhi == 0} {; set errlo 0; set errhi 0; } else {; set errlo [expr $paramlo - $param]; set errhi [expr $paramhi - $param]; }
#
!XSPEC12>expr $paramlo - $param
#
!XSPEC12>set errlo [expr $paramlo - $param]
#
!XSPEC12>expr $paramhi - $param
#
!XSPEC12>set errhi [expr $paramhi - $param]
#
!XSPEC12>switch $i {; 1 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^22\]"]}; 2 {puts $fname [format "%-16s%-11.2f%+-11.2f%+.2f" $pname $param $errlo $errhi]}; 3 - 4 {puts $fname [format "%-16s%-11.1f%s" $pname $param "\[keV\]"]}; 5 {puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^-12 ergs/cm^2/s\]"]}; }
#
!XSPEC12>format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^-12 ergs/cm^2/s\]"
#
!XSPEC12>puts $fname [format "%-16s%-11.2e%+-11.2e%+-11.2e%s" $pname $param $errlo $errhi "\[10^-12 ergs/cm^2/s\]"]
#
!XSPEC12>incr i
#
!XSPEC12>expr $param * 1.0e-12
#
!XSPEC12>set fl [expr $param * 1.0e-12];
#
!XSPEC12>expr $errlo * 1.0e-12
#
!XSPEC12>set flerrlo [expr $errlo * 1.0e-12];
#
!XSPEC12>expr $errhi * 1.0e-12
#
!XSPEC12>set flerrhi [expr $errhi * 1.0e-12];
#
!XSPEC12>tclout stat;
#
!XSPEC12>set chisq $xspec_tclout;
#
!XSPEC12>tclout dof;
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>set dof [lindex $xspec_tclout 0];
#
!XSPEC12>tclout stat;
#
!XSPEC12>expr $chisq/$dof
#
!XSPEC12>set rchisq [expr $chisq/$dof];
#
!XSPEC12>format "%-16s%.2f" "chisq" $chisq
#
!XSPEC12>puts $fname [format "%-16s%.2f" "chisq" $chisq];
#
!XSPEC12>format "%-16s%d" "dof" $dof
#
!XSPEC12>puts $fname [format "%-16s%d" "dof" $dof];
#
!XSPEC12>format "%-16s%.2f" "rchisq" $rchisq
#
!XSPEC12>puts $fname [format "%-16s%.2f" "rchisq" $rchisq];
#
!XSPEC12>puts $fname "";
#
!XSPEC12>expr $norm/$cr
#
!XSPEC12>set cv [expr $norm/$cr];
#
!XSPEC12>expr $normerrlo/$cr
#
!XSPEC12>set cverrlo [expr $normerrlo/$cr];
#
!XSPEC12>expr $normerrhi/$cr
#
!XSPEC12>set cverrhi [expr $normerrhi/$cr];
#
!XSPEC12>format "%-38s%-11.2e%+-11.2e%+-11.2e%s" "conversion_factor (norm_at_1keV)" $cv $cverrlo $cverrhi "\[cm^2/keV\]"
#
!XSPEC12>puts $fname [format "%-38s%-11.2e%+-11.2e%+-11.2e%s" "conversion_factor (norm_at_1keV)" $cv $cverrlo $cverrhi "\[cm^2/keV\]"];
#
!XSPEC12>expr $fl/$cr
#
!XSPEC12>set cv [expr $fl/$cr];
#
!XSPEC12>expr $flerrlo/$cr
#
!XSPEC12>set cverrlo [expr $flerrlo/$cr];
#
!XSPEC12>expr $flerrhi/$cr
#
!XSPEC12>set cverrhi [expr $flerrhi/$cr];
#
!XSPEC12>format "%-38s%-11.2e%+-11.2e%+-11.2e%s" "conversion_factor (unabsorbed_flux)" $cv $cverrlo $cverrhi "\[ergs/cm^2/count\]"
#
!XSPEC12>puts $fname [format "%-38s%-11.2e%+-11.2e%+-11.2e%s" "conversion_factor (unabsorbed_flux)" $cv $cverrlo $cverrhi "\[ergs/cm^2/count\]"];
#
!XSPEC12>dummyrsp 0.01 200.0;
#***Warning! One or more spectra are missing responses,
# and are not suitable for fit.
#
!XSPEC12>flux 14.0 194.0;
# Model Flux 0.0053308 photons (3.1051e-10 ergs/cm^2/s) range (14.000 - 194.00 keV)
#
!XSPEC12>tclout flux;
#
!XSPEC12>lindex $xspec_tclout 0
#
!XSPEC12>format "%.2e" [lindex $xspec_tclout 0]
#
!XSPEC12>set fl [format "%.2e" [lindex $xspec_tclout 0]];
#
!XSPEC12>format "%-38s%-11.2e%s" "extrapolated_flux (14-194 keV)" $fl "\[ergs/cm^2/s\]"
#
!XSPEC12>puts $fname [format "%-38s%-11.2e%s" "extrapolated_flux (14-194 keV)" $fl "\[ergs/cm^2/s\]"];
#
!XSPEC12>close $fname;
#
!XSPEC12>log none;
#