XSPEC version: 12.11.0 Build Date/Time: Mon Mar 29 15:43:38 2021 XSPEC12># Load spectra XSPEC12>@load_spectra !XSPEC12>data epoch01_grouped_spec_A.pha 1 spectrum in use Spectral Data File: epoch01_grouped_spec_A.pha Spectrum 1 Net count rate (cts/s) for Spectrum:1 1.847e-03 +/- 4.107e-04 (31.0 % total) Assigned to Data Group 1 and Plot Group 1 Noticed Channels: 1-3933 Telescope: NuSTAR Instrument: FPMA Channel Type: PI Exposure Time: 4.027e+04 sec Using fit statistic: chi Using Background File nu90701312002A01_bk.pha Background Exposure Time: 4.027e+04 sec Using Response (RMF) File nu90701312002A01_sr.rmf for Source 1 Using Auxiliary Response (ARF) File nu90701312002A01_sr.arf !XSPEC12>data 2:2 epoch01_grouped_spec_B.pha 2 spectra in use Spectral Data File: epoch01_grouped_spec_B.pha Spectrum 2 Net count rate (cts/s) for Spectrum:2 2.227e-03 +/- 4.296e-04 (34.4 % total) Assigned to Data Group 2 and Plot Group 2 Noticed Channels: 1-3781 Telescope: NuSTAR Instrument: FPMB Channel Type: PI Exposure Time: 3.97e+04 sec Using fit statistic: chi Using Background File nu90701312002B01_bk.pha Background Exposure Time: 3.97e+04 sec Using Response (RMF) File nu90701312002B01_sr.rmf for Source 1 Using Auxiliary Response (ARF) File nu90701312002B01_sr.arf !XSPEC12>ignore bad ignore: 3929 channels ignored from source number 1 ignore: 3774 channels ignored from source number 2 !XSPEC12>ignore **-3.5 78.-** 48 channels (1-48) ignored in spectrum # 1 2187 channels (1747-3933) ignored in spectrum # 1 !XSPEC12>ignore 2: **-3.5 78.-** 48 channels (1-48) ignored in spectrum # 2 2187 channels (1595-3781) ignored in spectrum # 2 XSPEC12># Load model XSPEC12>@const_phabs_pow !XSPEC12>method leven 1000 0.01 !XSPEC12>abund aspl Solar Abundance Vector set to aspl: Asplund M, Grevesse N., Sauval A.J. & Scott P., 2009, ARAA, 47, 481 !XSPEC12>xsect vern Cross Section Table set to vern: Verner, Ferland, Korista, and Yakovlev 1996 !XSPEC12>cosmo 70 0 0.73 !XSPEC12>xset delta 0.01 Parameter delta values for fits will now be 0.01 * parValue !XSPEC12>xset APEC_TRACE_ABUND Fe !XSPEC12>systematic 0 Model systematic error set to 0 !XSPEC12>model constant*phabs*powerlaw ; 1 -0.01 0 0 1e+10 1e+10 ; 0.493154 0.001 0 0 100000 1e+06 ; 3.36982 0.01 -3 -2 9 10 ; 0.000544694 0.01 0 0 1e+20 1e+24 ; 1 -0.01 0 0 1e+10 1e+10 ; = p2 ; = p3 ; = p4 ======================================================================== Model constant<1>*phabs<2>*powerlaw<3> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp Data group: 1 1 1 constant factor 1.00000 frozen 2 2 phabs nH 10^22 0.493154 +/- 0.0 3 3 powerlaw PhoIndex 3.36982 +/- 0.0 4 3 powerlaw norm 5.44694E-04 +/- 0.0 Data group: 2 5 1 constant factor 1.00000 frozen 6 2 phabs nH 10^22 0.493154 = p2 7 3 powerlaw PhoIndex 3.36982 = p3 8 3 powerlaw norm 5.44694E-04 = p4 ________________________________________________________________________ Fit statistic : Chi-Squared 2.48 using 3 bins. Chi-Squared 4.75 using 6 bins. Total fit statistic 7.23 with 6 d.o.f. Test statistic : Chi-Squared 7.23 using 9 bins. Null hypothesis probability of 3.00e-01 with 6 degrees of freedom Current data and model not fit yet. !XSPEC12>bayes off Fit statistic : Chi-Squared 2.48 using 3 bins. Chi-Squared 4.75 using 6 bins. Total fit statistic 7.23 with 6 d.o.f. Test statistic : Chi-Squared 7.23 using 9 bins. Null hypothesis probability of 3.00e-01 with 6 degrees of freedom Current data and model not fit yet. XSPEC12># Abundances and APEC_TRACE_ABUND are in the model, XSPEC12># but we can override here XSPEC12>#abund aspl XSPEC12>#xset APEC_TRACE_ABUND Fe XSPEC12># Refit the model XSPEC12>fit 1000 Parameters Chi-Squared |beta|/N Lvl 2:nH 3:PhoIndex 4:norm 7.22298 0.0420218 -1 0.0435215 3.37726 0.000549217 ======================================== Variances and Principal Axes 2 3 4 8.6125E-09| -0.0000 -0.0010 1.0000 5.2002E-01| -0.0347 0.9994 0.0010 2.2218E+03| -0.9994 -0.0347 -0.0000 ---------------------------------------- ==================================== Covariance Matrix 1 2 3 2.219e+03 7.692e+01 1.028e-01 7.692e+01 3.187e+00 4.079e-03 1.028e-01 4.079e-03 5.282e-06 ------------------------------------ ======================================================================== Model constant<1>*phabs<2>*powerlaw<3> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp Data group: 1 1 1 constant factor 1.00000 frozen 2 2 phabs nH 10^22 4.35215E-02 +/- 47.1073 3 3 powerlaw PhoIndex 3.37726 +/- 1.78520 4 3 powerlaw norm 5.49217E-04 +/- 2.29826E-03 Data group: 2 5 1 constant factor 1.00000 frozen 6 2 phabs nH 10^22 4.35215E-02 = p2 7 3 powerlaw PhoIndex 3.37726 = p3 8 3 powerlaw norm 5.49217E-04 = p4 ________________________________________________________________________ Fit statistic : Chi-Squared 2.51 using 3 bins. Chi-Squared 4.71 using 6 bins. Total fit statistic 7.22 with 6 d.o.f. Test statistic : Chi-Squared 7.22 using 9 bins. Null hypothesis probability of 3.01e-01 with 6 degrees of freedom XSPEC12>save model tmpmodel XSPEC12># Set plot parameters XSPEC12>setplot energy XSPEC12>setplot command lwidth 2 1 XSPEC12>setplot command lwidth 2 on 0..20 2 XSPEC12>setplot device const_phabs_pow.ps/cps XSPEC12>plot ldata delchi XSPEC12># Print the number of counts XSPEC12>show rates Spectral Data File: epoch01_grouped_spec_A.pha Assigned to Data Group 1 and Plot Group 1 Net count rate (cts/s) for Spectrum:1 5.929e-04 +/- 1.806e-04 (49.7 % total) Spectral data counts: 48 Model predicted rate: 8.51541E-04 Spectral Data File: epoch01_grouped_spec_B.pha Assigned to Data Group 2 and Plot Group 2 Net count rate (cts/s) for Spectrum:2 1.360e-03 +/- 2.484e-04 (60.0 % total) Spectral data counts: 90 Model predicted rate: 8.28644E-04 XSPEC12># Flux determination the new way XSPEC12> XSPEC12># freeze all the parameters XSPEC12>for {set i 1} {$i < 100} {incr i} { freeze $i } parameter 1 is already frozen. Fit statistic : Chi-Squared 2.51 using 3 bins. Chi-Squared 4.71 using 6 bins. Total fit statistic 7.22 with 7 d.o.f. Test statistic : Chi-Squared 7.22 using 9 bins. Null hypothesis probability of 4.06e-01 with 7 degrees of freedom Current data and model not fit yet. Fit statistic : Chi-Squared 2.51 using 3 bins. Chi-Squared 4.71 using 6 bins. Total fit statistic 7.22 with 8 d.o.f. Test statistic : Chi-Squared 7.22 using 9 bins. Null hypothesis probability of 5.13e-01 with 8 degrees of freedom Current data and model not fit yet. Fit statistic : Chi-Squared 2.51 using 3 bins. Chi-Squared 4.71 using 6 bins. Total fit statistic 7.22 with 9 d.o.f. Test statistic : Chi-Squared 7.22 using 9 bins. Null hypothesis probability of 6.14e-01 with 9 degrees of freedom Current data and model not fit yet. parameter 5 is already frozen. parameter 6 is linked and cannot be frozen. parameter 7 is linked and cannot be frozen. parameter 8 is linked and cannot be frozen. ***XSPEC Error: Invalid Range: Range value 9 is outside the allowed range limits (1-8) XSPEC12># edit model adding cflux() XSPEC12>editmod constant*phabs*cflux(powerlaw) Input parameter value, delta, min, bot, top, and max values for ... 0.5 -0.1( 0.005) 0 0 1e+06 1e+06 3:data group 1::cflux:Emin>3.5 10 -0.1( 0.1) 0 0 1e+06 1e+06 4:data group 1::cflux:Emax>78 -12 0.01( 0.12) -100 -100 100 100 5:data group 1::cflux:lg10Flux> Input parameter value, delta, min, bot, top, and max values for ... 0.5 -0.1( 0.005) 0 0 1e+06 1e+06 10:data group 2::cflux:Emin> 10 -0.1( 0.1) 0 0 1e+06 1e+06 11:data group 2::cflux:Emax> -12 0.01( 0.12) -100 -100 100 100 12:data group 2::cflux:lg10Flux> Fit statistic : Chi-Squared 1500.68 using 3 bins. Chi-Squared 677.48 using 6 bins. Total fit statistic 2178.15 with 8 d.o.f. Test statistic : Chi-Squared 2178.15 using 9 bins. Null hypothesis probability of 0.00e+00 with 8 degrees of freedom Current data and model not fit yet. ======================================================================== Model constant<1>*phabs<2>*cflux<3>*powerlaw<4> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp Data group: 1 1 1 constant factor 1.00000 frozen 2 2 phabs nH 10^22 4.35215E-02 frozen 3 3 cflux Emin keV 3.50000 frozen 4 3 cflux Emax keV 78.0000 frozen 5 3 cflux lg10Flux cgs -12.0000 +/- 0.0 6 4 powerlaw PhoIndex 3.37726 frozen 7 4 powerlaw norm 5.49217E-04 frozen Data group: 2 8 1 constant factor 1.00000 frozen 9 2 phabs nH 10^22 4.35215E-02 = p2 10 3 cflux Emin keV 3.50000 = p3 11 3 cflux Emax keV 78.0000 = p4 12 3 cflux lg10Flux cgs -12.0000 = p5 13 4 powerlaw PhoIndex 3.37726 = p6 14 4 powerlaw norm 5.49217E-04 = p7 ________________________________________________________________________ XSPEC12># fit the model to compute the unabsorbed flux XSPEC12>fit 1000 Warning: renorm - no variable model to allow renormalization Parameters Chi-Squared |beta|/N Lvl 5:lg10Flux 261.929 5709.22 -3 -12.3799 29.588 816.129 -4 -12.6919 8.12489 117.974 -5 -12.8815 7.22559 15.341 -6 -12.9405 7.22002 1.05157 -7 -12.9458 ==================== Variances and Principal Axes 5 5.0436E-03| 1.0000 -------------------- ============ Covariance Matrix 1 5.044e-03 ------------ ======================================================================== Model constant<1>*phabs<2>*cflux<3>*powerlaw<4> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp Data group: 1 1 1 constant factor 1.00000 frozen 2 2 phabs nH 10^22 4.35215E-02 frozen 3 3 cflux Emin keV 3.50000 frozen 4 3 cflux Emax keV 78.0000 frozen 5 3 cflux lg10Flux cgs -12.9458 +/- 7.10180E-02 6 4 powerlaw PhoIndex 3.37726 frozen 7 4 powerlaw norm 5.49217E-04 frozen Data group: 2 8 1 constant factor 1.00000 frozen 9 2 phabs nH 10^22 4.35215E-02 = p2 10 3 cflux Emin keV 3.50000 = p3 11 3 cflux Emax keV 78.0000 = p4 12 3 cflux lg10Flux cgs -12.9458 = p5 13 4 powerlaw PhoIndex 3.37726 = p6 14 4 powerlaw norm 5.49217E-04 = p7 ________________________________________________________________________ Fit statistic : Chi-Squared 2.64 using 3 bins. Chi-Squared 4.58 using 6 bins. Total fit statistic 7.22 with 8 d.o.f. Test statistic : Chi-Squared 7.22 using 9 bins. Null hypothesis probability of 5.13e-01 with 8 degrees of freedom XSPEC12># restore the original model XSPEC12>editmod constant*phabs*powerlaw Fit statistic : Chi-Squared 2.51 using 3 bins. Chi-Squared 4.71 using 6 bins. Total fit statistic 7.22 with 9 d.o.f. Test statistic : Chi-Squared 7.22 using 9 bins. Null hypothesis probability of 6.14e-01 with 9 degrees of freedom Current data and model not fit yet. ======================================================================== Model constant<1>*phabs<2>*powerlaw<3> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp Data group: 1 1 1 constant factor 1.00000 frozen 2 2 phabs nH 10^22 4.35215E-02 frozen 3 3 powerlaw PhoIndex 3.37726 frozen 4 3 powerlaw norm 5.49217E-04 frozen Data group: 2 5 1 constant factor 1.00000 frozen 6 2 phabs nH 10^22 4.35215E-02 = p2 7 3 powerlaw PhoIndex 3.37726 = p3 8 3 powerlaw norm 5.49217E-04 = p4 ________________________________________________________________________ XSPEC12># compute the absorbed flux XSPEC12>editmod cflux(constant*phabs*powerlaw) Input parameter value, delta, min, bot, top, and max values for ... 0.5 -0.1( 0.005) 0 0 1e+06 1e+06 1:data group 1::cflux:Emin>3.5 10 -0.1( 0.1) 0 0 1e+06 1e+06 2:data group 1::cflux:Emax>78 -12 0.01( 0.12) -100 -100 100 100 3:data group 1::cflux:lg10Flux> Input parameter value, delta, min, bot, top, and max values for ... 0.5 -0.1( 0.005) 0 0 1e+06 1e+06 8:data group 2::cflux:Emin> 10 -0.1( 0.1) 0 0 1e+06 1e+06 9:data group 2::cflux:Emax> -12 0.01( 0.12) -100 -100 100 100 10:data group 2::cflux:lg10Flux> Fit statistic : Chi-Squared 1504.34 using 3 bins. Chi-Squared 679.32 using 6 bins. Total fit statistic 2183.67 with 8 d.o.f. Test statistic : Chi-Squared 2183.67 using 9 bins. Null hypothesis probability of 0.00e+00 with 8 degrees of freedom Current data and model not fit yet. ======================================================================== Model cflux<1>(constant<2>*phabs<3>*powerlaw<4>) Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp Data group: 1 1 1 cflux Emin keV 3.50000 frozen 2 1 cflux Emax keV 78.0000 frozen 3 1 cflux lg10Flux cgs -12.0000 +/- 0.0 4 2 constant factor 1.00000 frozen 5 3 phabs nH 10^22 4.35215E-02 frozen 6 4 powerlaw PhoIndex 3.37726 frozen 7 4 powerlaw norm 5.49217E-04 frozen Data group: 2 8 1 cflux Emin keV 3.50000 = p1 9 1 cflux Emax keV 78.0000 = p2 10 1 cflux lg10Flux cgs -12.0000 = p3 11 2 constant factor 1.00000 frozen 12 3 phabs nH 10^22 4.35215E-02 = p5 13 4 powerlaw PhoIndex 3.37726 = p6 14 4 powerlaw norm 5.49217E-04 = p7 ________________________________________________________________________ XSPEC12> XSPEC12>fit 1000 Warning: renorm - no variable model to allow renormalization Parameters Chi-Squared |beta|/N Lvl 3:lg10Flux 262.629 5722.9 -3 -12.3799 29.662 818.066 -4 -12.6921 8.12923 118.255 -5 -12.8819 7.22564 15.3828 -6 -12.9410 7.22002 1.05585 -7 -12.9463 ==================== Variances and Principal Axes 3 5.0430E-03| 1.0000 -------------------- ============ Covariance Matrix 1 5.043e-03 ------------ ======================================================================== Model cflux<1>(constant<2>*phabs<3>*powerlaw<4>) Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp Data group: 1 1 1 cflux Emin keV 3.50000 frozen 2 1 cflux Emax keV 78.0000 frozen 3 1 cflux lg10Flux cgs -12.9463 +/- 7.10144E-02 4 2 constant factor 1.00000 frozen 5 3 phabs nH 10^22 4.35215E-02 frozen 6 4 powerlaw PhoIndex 3.37726 frozen 7 4 powerlaw norm 5.49217E-04 frozen Data group: 2 8 1 cflux Emin keV 3.50000 = p1 9 1 cflux Emax keV 78.0000 = p2 10 1 cflux lg10Flux cgs -12.9463 = p3 11 2 constant factor 1.00000 frozen 12 3 phabs nH 10^22 4.35215E-02 = p5 13 4 powerlaw PhoIndex 3.37726 = p6 14 4 powerlaw norm 5.49217E-04 = p7 ________________________________________________________________________ Fit statistic : Chi-Squared 2.64 using 3 bins. Chi-Squared 4.58 using 6 bins. Total fit statistic 7.22 with 8 d.o.f. Test statistic : Chi-Squared 7.22 using 9 bins. Null hypothesis probability of 5.13e-01 with 8 degrees of freedom XSPEC12># restore the original model XSPEC12>editmod constant*phabs*powerlaw Fit statistic : Chi-Squared 2.51 using 3 bins. Chi-Squared 4.71 using 6 bins. Total fit statistic 7.22 with 9 d.o.f. Test statistic : Chi-Squared 7.22 using 9 bins. Null hypothesis probability of 6.14e-01 with 9 degrees of freedom Current data and model not fit yet. ======================================================================== Model constant<1>*phabs<2>*powerlaw<3> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp Data group: 1 1 1 constant factor 1.00000 frozen 2 2 phabs nH 10^22 4.35215E-02 frozen 3 3 powerlaw PhoIndex 3.37726 frozen 4 3 powerlaw norm 5.49217E-04 frozen Data group: 2 5 1 constant factor 1.00000 frozen 6 2 phabs nH 10^22 4.35215E-02 = p2 7 3 powerlaw PhoIndex 3.37726 = p3 8 3 powerlaw norm 5.49217E-04 = p4 ________________________________________________________________________ XSPEC12> XSPEC12># Flux determination the old way XSPEC12> XSPEC12># Absorved flux over the NuSTAR band XSPEC12>flux 3.5 78.0 Spectrum Number: 1 Data Group Number: 1 Model Flux 1.1733e-05 photons (1.1209e-13 ergs/cm^2/s) range (3.5000 - 78.000 keV) Spectrum Number: 2 Data Group Number: 2 Model Flux 1.1733e-05 photons (1.1209e-13 ergs/cm^2/s) range (3.5000 - 78.000 keV) XSPEC12># Absorbed flux over a wide band XSPEC12>energies extend low 0.1 Models will use response energies extended to: Low: 0.1 in 200 log bins Fit statistic : Chi-Squared 2.51 using 3 bins. Chi-Squared 4.71 using 6 bins. Total fit statistic 7.22 with 9 d.o.f. Test statistic : Chi-Squared 7.22 using 9 bins. Null hypothesis probability of 6.14e-01 with 9 degrees of freedom Current data and model not fit yet. XSPEC12>energies extend high 500 Models will use response energies extended to: Low: 0.1 in 200 log bins High: 500.0 in 200 log bins Fit statistic : Chi-Squared 2.51 using 3 bins. Chi-Squared 4.71 using 6 bins. Total fit statistic 7.22 with 9 d.o.f. Test statistic : Chi-Squared 7.22 using 9 bins. Null hypothesis probability of 6.14e-01 with 9 degrees of freedom Current data and model not fit yet. XSPEC12>flux 0.3 78.0 Spectrum Number: 1 Data Group Number: 1 Model Flux 0.0025615 photons (2.4218e-12 ergs/cm^2/s) range (0.30000 - 78.000 keV) Spectrum Number: 2 Data Group Number: 2 Model Flux 0.0025615 photons (2.4218e-12 ergs/cm^2/s) range (0.30000 - 78.000 keV) XSPEC12># Unabsorbed flux over the NuSTAR band XSPEC12>newpar 2 0.0 Fit statistic : Chi-Squared 2.52 using 3 bins. Chi-Squared 4.70 using 6 bins. Total fit statistic 7.22 with 9 d.o.f. Test statistic : Chi-Squared 7.22 using 9 bins. Null hypothesis probability of 6.14e-01 with 9 degrees of freedom Current data and model not fit yet. XSPEC12> XSPEC12>flux 3.5 78.0 Spectrum Number: 1 Data Group Number: 1 Model Flux 1.1751e-05 photons (1.1222e-13 ergs/cm^2/s) range (3.5000 - 78.000 keV) Spectrum Number: 2 Data Group Number: 2 Model Flux 1.1751e-05 photons (1.1222e-13 ergs/cm^2/s) range (3.5000 - 78.000 keV) XSPEC12># Unabsorbed flux over a wide band XSPEC12>flux 0.3 78.0 Spectrum Number: 1 Data Group Number: 1 Model Flux 0.0040434 photons (3.353e-12 ergs/cm^2/s) range (0.30000 - 78.000 keV) Spectrum Number: 2 Data Group Number: 2 Model Flux 0.0040434 photons (3.353e-12 ergs/cm^2/s) range (0.30000 - 78.000 keV) XSPEC12>exit XSPEC: quit