XSPEC12>model (apec+apec)*wabs Input parameter value, delta, min, bot, top, and max values for ... 1 0.01 0.008 0.008 64 64 1:apec:kT>0.16 1 -0.001 0 0 5 5 2:apec:Abundanc> 0 -0.01 0 0 10 10 3:apec:redshift> 1 0.01 0 0 1e+24 1e+24 4:apec:norm>.0007 1 0.01 0.008 0.008 64 64 5:apec:kT>.39 1 -0.001 0 0 5 5 6:apec:Abundanc> 0 -0.01 0 0 10 10 7:apec:redshift> 1 0.01 0 0 1e+24 1e+24 8:apec:norm>0.0016 1 0.001 0 0 100000 1e+06 9:wabs:nH>.018 ======================================================================== Model (apec<1> + apec<2>)wabs<3> Source No.: 1 Active/Off Model Model Component Parameter Unit Value par comp 1 1 apec kT keV 0.160000 +/- 0.0 2 1 apec Abundanc 1.00000 frozen 3 1 apec redshift 0.0 frozen 4 1 apec norm 7.00000E-04 +/- 0.0 5 2 apec kT keV 0.390000 +/- 0.0 6 2 apec Abundanc 1.00000 frozen 7 2 apec redshift 0.0 frozen 8 2 apec norm 1.60000E-03 +/- 0.0 9 3 wabs nH 10^22 1.80000E-02 +/- 0.0 ________________________________________________________________________ XSPEC12>fakeit none For fake spectrum #1 response file is needed: meg_p1.rmf ...and ancillary file: acisf02575MEG_1_garf.fits Use counting statistics in creating fake data? (y): Input optional fake file prefix: tcar_v1_mp Fake data file name (tcar_v1_mpmeg_p1.fak): tcar_v1_mp.fak Exposure time, correction norm (1.00000, 1.00000): 200000. No background will be applied to fake spectrum #1 Warning: RMF CHANTYPE keyword is not consistent with spectrum OK, plot looks decent. Let's start again, and make a data file with all four spectrum MEG and HEG +/- 1. XSPEC12>model (apec+apec)*wabs Input parameter value, delta, min, bot, top, and max values for ... 1 0.01 0.008 0.008 64 64 1:apec:kT>.16 1 -0.001 0 0 5 5 2:apec:Abundanc> 0 -0.01 0 0 10 10 3:apec:redshift> 1 0.01 0 0 1e+24 1e+24 4:apec:norm>.0007 1 0.01 0.008 0.008 64 64 5:apec:kT>.39 1 -0.001 0 0 5 5 6:apec:Abundanc> 0 -0.01 0 0 10 10 7:apec:redshift> 1 0.01 0 0 1e+24 1e+24 8:apec:norm>.0016 1 0.001 0 0 100000 1e+06 9:wabs:nH>.018 ======================================================================== Model (apec<1> + apec<2>)wabs<3> Source No.: 1 Active/Off Model Model Component Parameter Unit Value par comp 1 1 apec kT keV 0.160000 +/- 0.0 2 1 apec Abundanc 1.00000 frozen 3 1 apec redshift 0.0 frozen 4 1 apec norm 7.00000E-04 +/- 0.0 5 2 apec kT keV 0.390000 +/- 0.0 6 2 apec Abundanc 1.00000 frozen 7 2 apec redshift 0.0 frozen 8 2 apec norm 1.60000E-03 +/- 0.0 9 3 wabs nH 10^22 1.80000E-02 +/- 0.0 ________________________________________________________________________ XSPEC12>fakeit 4 For fake spectrum #1 response file is needed: meg_n1.rmf ...and ancillary file: acisf02575MEG_-1_garf.fits For fake spectrum #2 response file is needed: meg_p1.rmf ...and ancillary file: acisf02575MEG_1_garf.fits For fake spectrum #3 response file is needed: heg_m1.rmf ...and ancillary file: acis640HEG_-1_garf.fits For fake spectrum #4 response file is needed: heg_p1.rmf ...and ancillary file: acis640HEG_1_garf.fits Use counting statistics in creating fake data? (y): Input optional fake file prefix: Fake data file name (meg_n1.fak): all4_v1.fak Exposure time, correction norm (1.00000, 1.00000): 200000. Fake data file name (meg_p1.fak): Exposure time, correction norm (1.00000, 1.00000): 200000. Fake data file name (heg_m1.fak): Exposure time, correction norm (1.00000, 1.00000): 200000. Fake data file name (heg_p1.fak): Exposure time, correction norm (1.00000, 1.00000): 200000. No background will be applied to fake spectrum #1 Warning: RMF CHANTYPE keyword is not consistent with spectrum No background will be applied to fake spectrum #2 Warning: RMF CHANTYPE keyword is not consistent with spectrum No background will be applied to fake spectrum #3 Warning: RMF CHANTYPE keyword is not consistent with spectrum No background will be applied to fake spectrum #4 Warning: RMF CHANTYPE keyword is not consistent with spectrum XSPEC12>@setup !XSPEC12>setplot wave; !XSPEC12>query yes; !XSPEC12>statistic cstat; C-statistic = 6701.37 using 32768 PHA bins and 32763 degrees of freedom. Valid fit does not exist. !XSPEC12>cpd /xw; !XSPEC12>setplot command log off; loading the faked datasets back in: !XSPEC12>data 1 meg_n1.fak; Warning: RMF CHANTYPE keyword is not consistent with spectrum 1 spectrum in use Source File: meg_n1.fak Net count rate (cts/s) for Spectrum:1 1.491e-02 +/- 2.731e-04 Assigned to Data Group 1 and Plot Group 1 Noticed Channels: 1-8192 Telescope: CHANDRA Instrument: ACIS Channel Type: PHA Exposure Time: 2e+05 sec Using Response (RMF) File meg_n1.rmf for Source 1 Using Auxiliary Response (ARF) File acisf02575MEG_-1_garf.fits !XSPEC12>data 2 meg_p1.fak; Warning: RMF CHANTYPE keyword is not consistent with spectrum 2 spectra in use Source File: meg_p1.fak Net count rate (cts/s) for Spectrum:2 7.545e-03 +/- 1.942e-04 Assigned to Data Group 1 and Plot Group 2 Noticed Channels: 1-8192 Telescope: CHANDRA Instrument: ACIS Channel Type: PHA Exposure Time: 2e+05 sec Using Response (RMF) File meg_p1.rmf for Source 1 Using Auxiliary Response (ARF) File acisf02575MEG_1_garf.fits !XSPEC12>data 3 heg_m1.fak; Warning: RMF CHANTYPE keyword is not consistent with spectrum 3 spectra in use Source File: heg_m1.fak Net count rate (cts/s) for Spectrum:3 1.990e-03 +/- 9.975e-05 Assigned to Data Group 1 and Plot Group 3 Noticed Channels: 1-8192 Telescope: CHANDRA Instrument: ACIS Channel Type: PHA Exposure Time: 2e+05 sec Using Response (RMF) File heg_m1.rmf for Source 1 Using Auxiliary Response (ARF) File acis640HEG_-1_garf.fits !XSPEC12>data 4 heg_p1.fak; Warning: RMF CHANTYPE keyword is not consistent with spectrum 4 spectra in use Source File: heg_p1.fak Net count rate (cts/s) for Spectrum:4 1.445e-03 +/- 8.500e-05 Assigned to Data Group 1 and Plot Group 4 Noticed Channels: 1-8192 Telescope: CHANDRA Instrument: ACIS Channel Type: PHA Exposure Time: 2e+05 sec Using Response (RMF) File heg_p1.rmf for Source 1 Using Auxiliary Response (ARF) File acis640HEG_1_garf.fits and fitting a 2T + wabs apec model to the four faked datasets: Model (apec<1> + apec<2>)wabs<3> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 apec kT keV 0.182424 +/- 6388.45 2 1 apec Abundanc 1.00000 frozen 3 1 apec redshift 0.0 frozen 4 1 apec norm 1.02179E-03 +/- 348.250 5 2 apec kT keV 0.244428 +/- 4065.34 6 2 apec Abundanc 1.00000 frozen 7 2 apec redshift 0.0 frozen 8 2 apec norm 1.31947E-03 +/- 339.148 9 3 wabs nH 10^22 2.67781E-02 +/- 1.04741E-02 ________________________________________________________________________ C-statistic = 6698.82 using 32768 PHA bins and 32763 degrees of freedom. ----- 9Mar09 - starting some new simulations: *first, loading faked data and writing out individual, four column, data only files (for input to bigplot) There are 687 HEG counts and 4492 MEG counts Looking at some individual lines: O VIII Ly-alpha at 18.97 Fitting the continuum first: 18.7-18.88; 19.08-19.20 ======================================================================== Model powerlaw<1> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 powerlaw PhoIndex 2.00000 frozen 2 1 powerlaw norm 1.97019E-04 +/- 1.45064E-04 ________________________________________________________________________ C-statistic = 39.83 using 352 PHA bins and 351 degrees of freedom. XSPEC12>plot XSPEC12>goodness 1000 0.20% of realizations are < best fit statistic 39.83 XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 0.000131 0.000280 (-0.000065,0.000085) fitting the line on 18.93-19.01 770 MEG counts; 41 HEG counts ======================================================================== Model wgauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 wgauss sigma_l "mA" 3.15665 +/- 0.592017 2 1 wgauss delta_l "mA" 0.320322 +/- 0.320572 3 1 wgauss waveleng "A" 18.9690 frozen 4 1 wgauss norm 5.34708E-04 +/- 3.56942E-05 5 2 powerlaw PhoIndex 2.00000 frozen 6 2 powerlaw norm 1.97000E-04 frozen ________________________________________________________________________ C-statistic = 49.54 using 92 PHA bins and 89 degrees of freedom. XSPEC12>plot XSPEC12>goodness 100 51.00% of realizations are < best fit statistic 49.54 XSPEC12>error 1. 1 Parameter Confidence Range (1.0000000) 1 2.397473 3.787992 (-0.759180,0.631338) XSPEC12>error 2.7 1 Parameter Confidence Range (2.700000) 1 1.595471 4.162059 (-1.561182,1.005406) XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) 2 0.000691 0.640371 (-0.319632,0.320048) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 -0.204962 0.846859 (-0.525285,0.526536) XSPEC12>error 1. 4 Parameter Confidence Range (1.000000) 4 0.000516 0.000553 (-0.000019,0.000019) XSPEC12>error 2.7 4 Parameter Confidence Range (2.700000) 4 0.000504 0.000566 (-0.000030,0.000032) XSPEC12>iplot PLT> wdata /home/cohen/xspec/2009/props/thetaCar/OVIII_1897_bestfit_wgauss.dat *Mg XI fir complex fitting the continuum on 9.0-9.1; 9.4-9.5 Model powerlaw<1> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 powerlaw PhoIndex 2.00000 frozen 2 1 powerlaw norm 6.98086E-05 +/- 2.62022E-05 ________________________________________________________________________ C-statistic = 102.97 using 232 PHA bins and 231 degrees of freedom. XSPEC12>plot XSPEC12>goodness 100 nosim 51.00% of realizations are < best fit statistic 102.97 XSPEC12>error 1. 2 Parameter Confidence Range (1.0000000) 2 0.000057 0.000084 (-0.000012,0.000014) ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 1.47232 +/- 0.502874 2 1 hegauss G 0.728848 +/- 0.232053 3 1 hegauss sigma_v "km/s" 36.0290 +/- 15.6795 4 1 hegauss delta_v "km/s" 37.1995 +/- 15.6875 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 4.26194E-06 +/- 1.39853E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 6.98000E-05 frozen ________________________________________________________________________ C-statistic = 243.69 using 356 PHA bins and 351 degrees of freedom. XSPEC12>plot XSPEC12>error 1. 1 Parameter Confidence Range (1.000000) 1 0.922152 1.893120 (-0.550004,0.420965) XSPEC12>goodness 1000 nosim 94.80% of realizations are < best fit statistic 243.69 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2009/props/thetaCar/MgXI_hegauss_fit_v1.dat XSPEC12>goodness 1000 nosim 94.80% of realizations are < best fit statistic 243.69 XSPEC12>error 1. 1 Parameter Confidence Range (1.000000) 1 0.922152 1.893120 (-0.550004,0.420965) XSPEC12>error 2.7 1 Parameter Confidence Range (2.7000000) 1 0.747941 2.364554 (-0.724048,0.892564) XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) 2 0.580914 0.924531 (-0.147916,0.195700) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 0.398280 1.045274 (-0.330545,0.316449) -- better fit found while using the "error" cmd: ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 1.30850 +/- 0.464000 2 1 hegauss G 0.715194 +/- 0.160623 3 1 hegauss sigma_v "km/s" 38.7428 +/- 60.4509 4 1 hegauss delta_v "km/s" 49.2454 +/- 19.6841 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 4.13835E-06 +/- 8.72256E-07 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 6.98000E-05 frozen ________________________________________________________________________ Using energies from responses. C-statistic = 243.04 using 356 PHA bins and 351 degrees of freedom. -- but then "fit": ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 1.30851 +/- 0.445964 2 1 hegauss G 0.715195 +/- 0.219772 3 1 hegauss sigma_v "km/s" 38.7428 +/- 19.3348 4 1 hegauss delta_v "km/s" 49.2454 +/- 19.3384 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 4.25947E-06 +/- 1.34475E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 6.98000E-05 frozen ________________________________________________________________________ C-statistic = 243.24 using 356 PHA bins and 351 degrees of freedom. XSPEC12>iplot PLT> wdata /home/cohen/xspec/2009/props/thetaCar/MgXI_hegauss_fit_v1b.dat Hmm... maybe we should *exclude the HEG spectra* XSPEC12>ignore 3-4:0.-** XSPEC12>fit ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 1.05945 +/- 1.57147 2 1 hegauss G 0.277884 +/- 0.530302 3 1 hegauss sigma_v "km/s" 33.5621 +/- 556.711 4 1 hegauss delta_v "km/s" 34.0893 +/- 180.486 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.17349E-06 +/- 2.57075E-05 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 6.98000E-05 frozen ________________________________________________________________________ C-statistic = 107.15 using 118 PHA bins and 113 degrees of freedom. XSPEC12>goodness 100 nosim 32.00% of realizations are < best fit statistic 107.15 XSPEC12>error 1. 1 Parameter Confidence Range (1.0000000) 1 0.737728 1.701682 (-0.322278,0.641676) XSPEC12>error 2.7 1 Parameter Confidence Range (2.700000) 1 0.563241 2.244754 (-0.485967,1.195547) XSPEC12>steppar 2 0.2 0.8 30 C-Statistic Delta G C-Statistic 2 107.16551 0.07242 0 0.20000 106.53548 -0.55761 1 0.22000 106.64677 -0.44632 2 0.24000 106.95029 -0.14279 3 0.26000 107.43725 0.34416 4 0.28000 108.06146 0.96838 5 0.30000 108.82803 1.73494 6 0.32000 109.70913 2.61604 7 0.34000 110.68287 3.58978 8 0.36000 111.73364 4.64055 9 0.38000 112.84835 5.75526 10 0.40000 114.01597 6.92288 11 0.42000 113.28554 6.19245 12 0.44000 113.12049 6.02740 13 0.46000 113.07763 5.98454 14 0.48000 113.06729 5.97420 15 0.50000 113.08622 5.99314 16 0.52000 113.13521 6.04212 17 0.54000 112.82394 5.73086 18 0.56000 112.84376 5.75067 19 0.58000 112.89655 5.80347 20 0.60000 112.97883 5.88574 21 0.62000 113.08743 5.99435 22 0.64000 113.21962 6.12654 23 0.66000 113.37301 6.27992 24 0.68000 113.54548 6.45239 25 0.70000 113.73520 6.64212 26 0.72000 113.93117 6.83808 27 0.74000 114.15995 7.06686 28 0.76000 114.40386 7.31077 29 0.78000 114.66158 7.56850 30 0.80000 Argh! New best-fit at an unrealistically small value of G Note: there are 111 counts in the MEG spectrum in this complex. OK - at this point, let's remake the fake spectrum, but using hegauss + pow; that way, at least we'll know the input R, G, etc. (instead of relying on the apec output); but we'll use the normalization here. New xspec session: XSPEC12>model hegauss + pow ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/Off Model Model Component Parameter Unit Value par comp 1 1 hegauss R 2.70000 frozen 2 1 hegauss G 0.710000 frozen 3 1 hegauss sigma_v "km/s" 0.0 frozen 4 1 hegauss delta_v "km/s" 0.0 frozen 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.20000E-06 +/- 0.0 7 2 powerlaw PhoIndex 2.00000 +/- 0.0 8 2 powerlaw norm 7.00000E-05 +/- 0.0 ________________________________________________________________________ XSPEC12>fakeit 2 For fake spectrum #1 response file is needed: /home/cohen/xspec/2009/props/thetaCar/meg_n1.rmf ...and ancillary file: /home/cohen/xspec/2009/props/thetaCar/acisf02575MEG_-1_garf.fits For fake spectrum #2 response file is needed: /home/cohen/xspec/2009/props/thetaCar/meg_p1.rmf ...and ancillary file: /home/cohen/xspec/2009/props/thetaCar/acisf02575MEG_1_garf.fits Use counting statistics in creating fake data? (y): Input optional fake file prefix: /home/cohen/xspec/2009/props/thetaCar/mgxi Fake data file name (/home/cohen/xspec/2009/props/thetaCar/mgxi/home/cohen/xspec/2009/props/thetaCar/meg_n1.fak): mgxi_meg_n1.fak Exposure time, correction norm (1.00000, 1.00000): 200000. Fake data file name (/home/cohen/xspec/2009/props/thetaCar/mgxi/home/cohen/xspec/2009/props/thetaCar/meg_p1.fak): mgxi_meg_p1.fak Exposure time, correction norm (1.00000, 1.00000): 200000. No background will be applied to fake spectrum #1 Warning: RMF CHANTYPE keyword is not consistent with spectrum No background will be applied to fake spectrum #2 Warning: RMF CHANTYPE keyword is not consistent with spectrum (Note: accidentally wrote fake files into local models directory; had to move them to working directory - don't forget to give the full path for output files too) OK; while still in the same xspec session, I fit the faked data: ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 3.65427 +/- 1.32730 2 1 hegauss G 0.747751 +/- 0.138609 3 1 hegauss sigma_v "km/s" 0.0 frozen 4 1 hegauss delta_v "km/s" 0.0 frozen 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 6.85808E-06 +/- 1.20784E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 7.00000E-05 frozen ________________________________________________________________________ C-statistic = 125.58 using 118 PHA bins and 115 degrees of freedom. XSPEC12>plot XSPEC12>error 1. 1 Parameter Confidence Range (1) 1 2.58694 5.34323 (-1.06655,1.68974) XSPEC12>error 2.7 1 Parameter Confidence Range (2.7) 1 2.10083 7.11515 (-1.55341,3.4609) XSPEC12>error 4. 1 Parameter Confidence Range (4) 1 1.87619 8.34396 (-1.77759,4.69018) XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 0.620926 0.899788 (-0.126875,0.151987) XSPEC12>error 2.7 2 Parameter Confidence Range (2.7) 2 0.548409 1.01203 (-0.199343,0.264282) and wrote out the fake data and refitted model: XSPEC12>iplot PLT> wdata /home/cohen/xspec/2009/props/thetaCar/mgxi_hegauss_refit_v1.dat For future fitting of this faked dataset, see: load_mgxi_fak.xcm Note: 163 counts in MEG dataset ----- 12mar09 - making two more fake datasets, with R=0.1 and with R=1.2 XSPEC12>model hegauss + pow ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/Off Model Model Component Parameter Unit Value par comp 1 1 hegauss R 0.100000 frozen 2 1 hegauss G 0.710000 frozen 3 1 hegauss sigma_v "km/s" 0.0 frozen 4 1 hegauss delta_v "km/s" 0.0 frozen 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.20000E-06 +/- 0.0 7 2 powerlaw PhoIndex 2.00000 +/- 0.0 8 2 powerlaw norm 7.00000E-05 +/- 0.0 ________________________________________________________________________ XSPEC12>fakeit 4 For fake spectrum #1 response file is needed: /home/cohen/xspec/2009/props/thetaCar/meg_n1.rmf ...and ancillary file: /home/cohen/xspec/2009/props/thetaCar/acisf02575MEG_-1_garf.fits For fake spectrum #2 response file is needed: /home/cohen/xspec/2009/props/thetaCar/meg_p1.rmf ...and ancillary file: /home/cohen/xspec/2009/props/thetaCar/acisf02575MEG_1_garf.fits For fake spectrum #3 response file is needed: /home/cohen/xspec/2009/props/thetaCar/heg_m1.rmf ...and ancillary file: /home/cohen/xspec/2009/props/thetaCar/acis640HEG_-1_garf.fits For fake spectrum #4 response file is needed: /home/cohen/xspec/2009/props/thetaCar/heg_p1.rmf ...and ancillary file: /home/cohen/xspec/2009/props/thetaCar/acis640HEG_1_garf.fits Use counting statistics in creating fake data? (y): y Input optional fake file prefix: Enter new name (/home/cohen/xspec/2009/props/thetaCar/meg_n1.fak): R01_meg_n1.fak Exposure time, correction norm (1.00000, 1.00000): 200000. Fake data file name (/home/cohen/xspec/2009/props/thetaCar/meg_p1.fak): R01_meg_p1.fak Exposure time, correction norm (1.00000, 1.00000): 200000. Fake data file name (/home/cohen/xspec/2009/props/thetaCar/heg_m1.fak): R01_heg_n1.fak Exposure time, correction norm (1.00000, 1.00000): 200000. Fake data file name (/home/cohen/xspec/2009/props/thetaCar/heg_p1.fak): R01_heg_p1.fak Exposure time, correction norm (1.00000, 1.00000): 200000. refitting the R = 0.1 faked data: ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 7.96537E-02 +/- 5.10237E-02 2 1 hegauss G 0.591709 +/- 0.146161 3 1 hegauss sigma_v "km/s" 31.3589 +/- 6.18107E-02 4 1 hegauss delta_v "km/s" 16.3161 +/- 14.1106 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.81212E-06 +/- 2.29417E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 7.00000E-05 frozen ________________________________________________________________________ C-statistic = 191.39 using 356 PHA bins and 351 degrees of freedom. XSPEC12>goodness 1000 nosim 11.20% of realizations are < best fit statistic 191.39 XSPEC12>error 1. 1 Parameter Confidence Range (1.0000000) 1 0.040506 0.131472 (-0.039141,0.051825) XSPEC12>error 2.7 1 Parameter Confidence Range (2.700000) 1 0.020524 0.172331 (-0.059122,0.092685) XSPEC12>error 4. 1 Parameter Confidence Range (4.000000) 1 0.010989 0.198047 (-0.068656,0.118402) XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) 2 0.505104 0.687569 (-0.086685,0.095780) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 0.455850 0.757880 (-0.135924,0.166106) XSPEC12>error 4. 2 Parameter Confidence Range (4.000000) 2 0.430414 0.799279 (-0.161349,0.207516) XSPEC12>error 1. 3 Parameter Confidence Range (1.000000) Apparent non-monotonicity in statistic space detected. Current bracket values 31.361339, 31.299529 and delta stat 0.000000, 34.536965 but latest trial 31.311677 gives 34.538015 Suggest that you check this result using the steppar command. 3 31.330434 49.852226 (-0.031073,18.490719) XSPEC12>error 1. 4 Parameter Confidence Range (1.000000) ***Warning: Number of trials exceeded before convergence. Current trial values 26.908635, 26.908796 and delta statistic 0.874812, 1.338940 ***Warning: Identical values of the parameter give different values of the statistic. Please check your result for the high end of the confidence range. 4 1.534259 26.908767 (-14.852693,10.521816) XSPEC12>error 1. 6 Parameter Confidence Range (1.000000) 6 0.000007 0.000008 (-0.000001,0.000001) XSPEC12>error 2.7 6 Parameter Confidence Range (2.700000) 6 0.000007 0.000009 (-0.000001,0.000001) erXSPEC12>error 4.0 6 Parameter Confidence Range (4.000000) 6 0.000007 0.000009 (-0.000001,0.000001) XSPEC12>steppar 6 6.e-6 9.6e-6 30 C-Statistic Delta norm C-Statistic 6 203.26892 11.88332 0 0.00001 201.62563 10.24004 1 0.00001 200.12203 8.73643 2 0.00001 198.75492 7.36933 3 0.00001 197.51960 6.13400 4 0.00001 196.40928 5.02368 5 0.00001 195.42149 4.03589 6 0.00001 194.55114 3.16555 7 0.00001 193.79391 2.40831 8 0.00001 193.14716 1.76157 9 0.00001 192.60687 1.22128 10 0.00001 192.16971 0.78412 11 0.00001 191.83251 0.44692 12 0.00001 191.59226 0.20666 13 0.00001 191.44589 0.06030 14 0.00001 191.39005 0.00445 15 0.00001 191.42285 0.03726 16 0.00001 191.54139 0.15580 17 0.00001 191.74250 0.35690 18 0.00001 192.02318 0.63759 19 0.00001 192.16899 0.78339 20 0.00001 192.43888 1.05329 21 0.00001 192.77854 1.39295 22 0.00001 193.18957 1.80397 23 0.00001 193.66713 2.28154 24 0.00001 194.21675 2.83116 25 0.00001 194.83374 3.44814 26 0.00001 195.51544 4.12985 27 0.00001 196.26158 4.87598 28 0.00001 197.06797 5.68237 29 0.00001 197.93418 6.54858 30 0.00001 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2009/props/thetaCar/R01_hegauss.dat OK... that's all we need: the 4 faked spectra, a fit to them, confidence limits on that fit, and an ascii file containing four spectra plus the best-fit model. Moving on to the R=1.2 model (admittedly, a moderately arbitrary, and somewhat extreme choice; should probably do R=0.6 or something, too). ...in same xspec session, just defining a new model: ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 1.20000 frozen 2 1 hegauss G 0.710000 frozen 3 1 hegauss sigma_v "km/s" 0.0 frozen 4 1 hegauss delta_v "km/s" 0.0 frozen 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.20000E-06 +/- 0.0 7 2 powerlaw PhoIndex 2.00000 +/- 0.0 8 2 powerlaw norm 7.00000E-05 +/- 0.0 ________________________________________________________________________ C-statistic = 259.99 using 356 PHA bins and 353 degrees of freedom. Valid fit does not exist. XSPEC12>fakeit 4 Use counting statistics in creating fake data? (y): Input optional fake file prefix: Fake data file name (R01_meg_n1.fak.fak): R12_meg_n1.fak Exposure time, correction norm (200000., 1.00000): 200000. Fake data file name (R01_meg_p1.fak.fak): F12_meg_p1.fak Exposure time, correction norm (200000., 1.00000): 200000. Fake data file name (R01_heg_n1.fak.fak): R12_heg_n1.fak Exposure time, correction norm (200000., 1.00000): 200000. Fake data file name (R01_heg_p1.fak.fak): R12_heg_p1.fak Exposure time, correction norm (200000., 1.00000): 200000. No background will be applied to fake spectrum #1 No background will be applied to fake spectrum #2 No background will be applied to fake spectrum #3 No background will be applied to fake spectrum #4 ...written into local models folder; change one wrong name, and mv files ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 1.25511 +/- 0.335066 2 1 hegauss G 0.558550 +/- 0.133774 3 1 hegauss sigma_v "km/s" 30.5806 +/- 8.18227 4 1 hegauss delta_v "km/s" 12.9500 +/- 8.18774 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.19385E-06 +/- 1.96127E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 7.00000E-05 frozen ________________________________________________________________________ C-statistic = 180.83 using 356 PHA bins and 351 degrees of freedom. XSPEC12>goodness 1000 nosim 2.30% of realizations are < best fit statistic 180.83 XSPEC12>error 1. 1 Parameter Confidence Range (1.0000000) 1 0.973537 1.666874 (-0.283031,0.410306) XSPEC12>error 2.7 1 Parameter Confidence Range (2.700000) 1 0.820745 1.989120 (-0.437155,0.731220) XSPEC12>error 4.0 1 Parameter Confidence Range (4.000000) 1 0.746059 2.198747 (-0.513053,0.939636) XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) 2 0.469743 0.649486 (-0.087886,0.091856) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 0.421634 0.718663 (-0.135812,0.161218) XSPEC12>error 4.0 2 Parameter Confidence Range (4.000000) 2 0.396797 0.759914 (-0.160479,0.202638) XSPEC12>steppar 6 6.6e-6 8.0e-6 70 C-Statistic Delta norm C-Statistic 6 182.34653 1.53682 0 0.00001 182.24960 1.43989 1 0.00001 182.15860 1.34889 2 0.00001 182.07074 1.26103 3 0.00001 181.98586 1.17614 4 0.00001 181.90406 1.09435 5 0.00001 181.82535 1.01564 6 0.00001 181.74970 0.93999 7 0.00001 181.67709 0.86738 8 0.00001 181.60751 0.79780 9 0.00001 181.54094 0.73123 10 0.00001 181.47736 0.66765 11 0.00001 181.41676 0.60705 12 0.00001 181.35912 0.54941 13 0.00001 181.30442 0.49471 14 0.00001 181.25265 0.44294 15 0.00001 181.20379 0.39408 16 0.00001 181.15783 0.34811 17 0.00001 181.11474 0.30503 18 0.00001 181.07452 0.26480 19 0.00001 181.03714 0.22743 20 0.00001 181.00260 0.19289 21 0.00001 180.97087 0.16116 22 0.00001 180.94195 0.13224 23 0.00001 180.91374 0.10403 24 0.00001 180.89038 0.08067 25 0.00001 180.86978 0.06007 26 0.00001 180.85192 0.04221 27 0.00001 180.83679 0.02708 28 0.00001 180.82438 0.01466 29 0.00001 180.81466 0.00495 30 0.00001 180.80762 -0.00209 31 0.00001 180.80326 -0.00645 32 0.00001 180.80156 -0.00815 33 0.00001 180.80250 -0.00721 34 0.00001 180.80607 -0.00364 35 0.00001 180.81225 0.00254 36 0.00001 180.82103 0.01132 37 0.00001 180.83240 0.02269 38 0.00001 180.84634 0.03663 39 0.00001 180.86284 0.05313 40 0.00001 180.88189 0.07218 41 0.00001 180.90344 0.09373 42 0.00001 180.92758 0.11786 43 0.00001 180.95422 0.14451 44 0.00001 180.98335 0.17364 45 0.00001 181.01497 0.20526 46 0.00001 181.04905 0.23934 47 0.00001 181.08560 0.27588 48 0.00001 181.12459 0.31487 49 0.00001 181.16601 0.35630 50 0.00001 181.20986 0.40015 51 0.00001 181.25612 0.44641 52 0.00001 181.30478 0.49507 53 0.00001 181.35583 0.54612 54 0.00001 181.40926 0.59955 55 0.00001 181.46506 0.65535 56 0.00001 181.52321 0.71350 57 0.00001 181.58371 0.77400 58 0.00001 181.64654 0.83683 59 0.00001 181.71170 0.90199 60 0.00001 181.77917 0.96946 61 0.00001 181.84894 1.03923 62 0.00001 181.92100 1.11129 63 0.00001 181.99534 1.18563 64 0.00001 182.07192 1.26221 65 0.00001 182.15076 1.34105 66 0.00001 182.23183 1.42212 67 0.00001 182.31514 1.50543 68 0.00001 182.40066 1.59095 69 0.00001 182.48840 1.67869 70 0.00001 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2009/props/thetaCar/R12_hegauss.dat OK, now we'll do a simulation with R=0.6: ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 0.600000 frozen 2 1 hegauss G 0.710000 frozen 3 1 hegauss sigma_v "km/s" 0.0 frozen 4 1 hegauss delta_v "km/s" 0.0 frozen 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.20000E-06 +/- 0.0 7 2 powerlaw PhoIndex 2.00000 +/- 0.0 8 2 powerlaw norm 7.00000E-05 +/- 0.0 ________________________________________________________________________ XSPEC12>fakeit 4 Use counting statistics in creating fake data? (y): Input optional fake file prefix: R06_meg_n1.fak Fake data file name (R06_meg_n1.fakR12_meg_n1.fak.fak): R06_meg_n1.fak Exposure time, correction norm (200000., 1.00000): 200000. Fake data file name (R06_meg_n1.fakF12_meg_p1.fak.fak): R06_meg_p1.fak Exposure time, correction norm (200000., 1.00000): 200000. Fake data file name (R06_meg_n1.fakR12_heg_n1.fak.fak): R06_heg_n1.fak Exposure time, correction norm (200000., 1.00000): 200000. Fake data file name (R06_meg_n1.fakR12_heg_p1.fak.fak): R06_heg_p1.fak Exposure time, correction norm (200000., 1.00000): 200000. No background will be applied to fake spectrum #1 No background will be applied to fake spectrum #2 No background will be applied to fake spectrum #3 No background will be applied to fake spectrum #4 ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 0.442948 +/- 0.105362 2 1 hegauss G 0.755422 +/- 0.161550 3 1 hegauss sigma_v "km/s" 27.7652 +/- 10.6296 4 1 hegauss delta_v "km/s" 0.458424 +/- 10.6455 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 8.75684E-06 +/- 2.03421E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 7.00000E-05 frozen ________________________________________________________________________ C-statistic = 242.56 using 356 PHA bins and 351 degrees of freedom. XSPEC12>goodness 1000 nosim 88.20% of realizations are < best fit statistic 242.56 XSPEC12>error 1. 1 Parameter Confidence Range (1.0000000) 1 0.340871 0.550673 (-0.101385,0.108417) XSPEC12>error 2.7 1 Parameter Confidence Range (2.700000) 1 0.289554 0.637375 (-0.152074,0.195747) XSPEC12>error 4.0 1 Parameter Confidence Range (4.000000) 1 0.263523 0.690724 (-0.177536,0.249665) XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) 2 0.661072 0.877012 (-0.095096,0.120844) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 0.602701 0.959536 (-0.153613,0.203222) XSPEC12>error 4. 2 Parameter Confidence Range (4.000000) 2 0.572342 1.008387 (-0.184106,0.251939) XSPEC12>steppar 6 6.6e-6 9.6e-6 30 C-Statistic Delta norm C-Statistic 6 257.39795 14.84048 0 0.00001 255.91944 13.36198 1 0.00001 254.53133 11.97387 2 0.00001 253.23116 10.67369 3 0.00001 252.01655 9.45909 4 0.00001 250.88524 8.32777 5 0.00001 249.83503 7.27757 6 0.00001 248.86384 6.30638 7 0.00001 247.96965 5.41219 8 0.00001 247.15052 4.59306 9 0.00001 246.40457 3.84711 10 0.00001 245.73002 3.17256 11 0.00001 245.12515 2.56769 12 0.00001 244.58819 2.03073 13 0.00001 244.11754 1.56007 14 0.00001 243.71165 1.15418 15 0.00001 243.36903 0.81156 16 0.00001 243.08825 0.53079 17 0.00001 242.86503 0.30757 18 0.00001 242.70332 0.14585 19 0.00001 242.59999 0.04253 20 0.00001 242.55325 -0.00421 21 0.00001 242.56192 0.00445 22 0.00001 242.62474 0.06728 23 0.00001 242.74060 0.18314 24 0.00001 242.90838 0.35091 25 0.00001 243.12700 0.56953 26 0.00001 243.39542 0.83796 27 0.00001 243.71265 1.15518 28 0.00001 244.07769 1.52023 29 0.00001 244.48960 1.93214 30 0.00001 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2009/props/thetaCar/R06_hegauss.dat OK, fine. Maurice suggests it would be very good if we could get away with smaller exposure times. So... using the slightly optimistic normalization of 7.2e-6 for the Mg XI complex, let's simulate these same three cases using only 100 ks exposure: ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 0.100000 frozen 2 1 hegauss G 0.710000 frozen 3 1 hegauss sigma_v "km/s" 0.0 frozen 4 1 hegauss delta_v "km/s" 0.0 frozen 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.20000E-06 +/- 0.0 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 7.00000E-05 frozen ________________________________________________________________________ Using energies from responses. XSPEC12>fakeit 4 Use counting statistics in creating fake data? (y): Input optional fake file prefix: Fake data file name (R06_meg_n1.fak.fak): R01_meg_n1_100ks.fak Exposure time, correction norm (200000., 1.00000): 100000. Fake data file name (R06_meg_p1.fak.fak): R01_meg_p1_100ks.fak Exposure time, correction norm (200000., 1.00000): 100000. Fake data file name (R06_heg_n1.fak.fak): R01_heg_n1_100ks.fak Exposure time, correction norm (200000., 1.00000): 100000. Fake data file name (R06_heg_p1.fak.fak): R01_heg_p1_100ks.fak Exposure time, correction norm (200000., 1.00000): 100000. No background will be applied to fake spectrum #1 No background will be applied to fake spectrum #2 No background will be applied to fake spectrum #3 No background will be applied to fake spectrum #4 ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 0.156425 +/- 9.17325E-02 2 1 hegauss G 0.497701 +/- 0.110154 3 1 hegauss sigma_v "km/s" 74.8486 +/- 40.9976 4 1 hegauss delta_v "km/s" 46.4888 +/- 22.8159 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 8.21119E-06 +/- 1.70588E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 7.00000E-05 frozen ________________________________________________________________________ C-statistic = 164.40 using 356 PHA bins and 351 degrees of freedom. XSPEC12>error 1. 1 Parameter Confidence Range (1.000000) 1 0.079275 0.266836 (-0.077162,0.110398) XSPEC12>error 2.7 1 Parameter Confidence Range (2.700000) 1 0.042242 0.358945 (-0.114196,0.202507) XSPEC12>error 4. 1 Parameter Confidence Range (4.000000) 1 0.025297 0.420203 (-0.131141,0.263766) XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) 2 0.398083 0.618085 (-0.099629,0.120374) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 0.343095 0.708843 (-0.154617,0.211131) XSPEC12>error 4. 2 Parameter Confidence Range (4.000000) 2 0.315073 0.763053 (-0.182639,0.265342) XSPEC12>steppar 6 6.e-6 10.e-6 40 C-Statistic Delta norm C-Statistic 6 172.76120 8.36219 0 0.00001 171.94515 7.54614 1 0.00001 171.17686 6.77785 2 0.00001 170.45475 6.05574 3 0.00001 169.77770 5.37869 4 0.00001 169.14421 4.74520 5 0.00001 168.55311 4.15410 6 0.00001 167.99667 3.59766 7 0.00001 167.48548 3.08647 8 0.00001 167.01498 2.61597 9 0.00001 166.57990 2.18089 10 0.00001 166.18382 1.78482 11 0.00001 165.82545 1.42645 12 0.00001 165.50345 1.10444 13 0.00001 165.21692 0.81792 14 0.00001 164.96646 0.56745 15 0.00001 164.75043 0.35143 16 0.00001 164.56800 0.16899 17 0.00001 164.41834 0.01934 18 0.00001 164.30069 -0.09832 19 0.00001 164.21428 -0.18473 20 0.00001 164.15838 -0.24063 21 0.00001 164.13229 -0.26672 22 0.00001 164.13533 -0.26367 23 0.00001 164.16685 -0.23215 24 0.00001 164.22622 -0.17279 25 0.00001 164.31281 -0.08619 26 0.00001 164.42605 0.02704 27 0.00001 164.56535 0.16634 28 0.00001 164.73012 0.33111 29 0.00001 164.91812 0.51911 30 0.00001 165.11243 0.71342 31 0.00001 165.34292 0.94391 32 0.00001 165.59372 1.19472 33 0.00001 165.86988 1.47088 34 0.00001 166.16683 1.76782 35 0.00001 166.48698 2.08797 36 0.00001 166.82797 2.42897 37 0.00001 167.19133 2.79233 38 0.00001 167.57343 3.17442 39 0.00001 167.97826 3.57925 40 0.00001 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2009/props/thetaCar/R01_hegauss_100ks.dat Next, R=0.6 and 100 ks: ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 0.600000 frozen 2 1 hegauss G 0.710000 frozen 3 1 hegauss sigma_v "km/s" 0.0 frozen 4 1 hegauss delta_v "km/s" 0.0 frozen 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.20000E-06 +/- 0.0 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 7.00000E-05 frozen ________________________________________________________________________ XSPEC12>fakeit 4 Use counting statistics in creating fake data? (y): Input optional fake file prefix: Fake data file name (R01_meg_n1_100ks.fak.fak): R06_meg_n1_100ks.fak Exposure time, correction norm (100000., 1.00000): 100000. Fake data file name (R01_meg_p1_100ks.fak.fak): R06_meg_p1_100ks.fak Exposure time, correction norm (100000., 1.00000): 100000. Fake data file name (R01_heg_n1_100ks.fak.fak): R06_heg_n1_100ks.fak Exposure time, correction norm (100000., 1.00000): 100000. Fake data file name (R01_heg_p1_100ks.fak.fak): R06_heg_p1_100ks.fak Exposure time, correction norm (100000., 1.00000): 100000. No background will be applied to fake spectrum #1 No background will be applied to fake spectrum #2 No background will be applied to fake spectrum #3 No background will be applied to fake spectrum #4 ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 1.23683 +/- 0.470736 2 1 hegauss G 0.706617 +/- 0.160505 3 1 hegauss sigma_v "km/s" 30.2621 +/- 39.0698 4 1 hegauss delta_v "km/s" 11.3116 +/- 34.2984 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.70745E-06 +/- 1.74504E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 7.00000E-05 frozen ________________________________________________________________________ C-statistic = 156.73 using 356 PHA bins and 351 degrees of freedom. ...hmmm, that's a pretty big R... XSPEC12>goodness 1000 nosim 16.00% of realizations are < best fit statistic 156.73 XSPEC12>error 1. 1 Parameter Confidence Range (1.0000000) 1 0.864824 1.766316 (-0.381529,0.519963) XSPEC12>error 2.7 1 Parameter Confidence Range (2.700000) 1 0.671164 2.214840 (-0.575198,0.968478) XSPEC12>error 4. 1 Parameter Confidence Range (4.000000) 1 0.586831 2.525132 (-0.659910,1.278391) *Maybe the very low count HEG data is fucking things up; let's look just at the MEG: ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 1.73540 +/- 9.14602 2 1 hegauss G 0.660956 +/- 1.30122 3 1 hegauss sigma_v "km/s" 28.1970 +/- 386.574 4 1 hegauss delta_v "km/s" 1.76902 +/- 543.953 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.81665E-06 +/- 1.23529E-05 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 7.00000E-05 frozen ________________________________________________________________________ C-statistic = 83.25 using 118 PHA bins and 113 degrees of freedom. XSPEC12>plot XSPEC12>error 1. 1 Parameter Confidence Range (1.000000) Apparent non-monotonicity in statistic space detected. Current bracket values 0.951223, 0.948793 and delta stat 0.956316, 1.018293 but latest trial 0.949814 gives 0.894657 Suggest that you check this result using the steppar command. 1 0.950008 2.769020 (-0.788521,1.030491) So, no, that's not the problem. Maybe we've just got a statistical outlier simulation. Let's rerun fakeit. ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 0.600000 frozen 2 1 hegauss G 0.710000 frozen 3 1 hegauss sigma_v "km/s" 0.0 frozen 4 1 hegauss delta_v "km/s" 0.0 frozen 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.20000E-06 +/- 0.0 7 2 powerlaw PhoIndex 2.00000 +/- 0.0 8 2 powerlaw norm 7.00000E-05 +/- 0.0 ________________________________________________________________________ XSPEC12>fakeit 4 Use counting statistics in creating fake data? (y): Input optional fake file prefix: Fake data file name (R06_meg_n1_100ks.fak.fak): R06_meg_n1_100ks.fak Exposure time, correction norm (100000., 1.00000): 100000. Fake data file name (R06_meg_p1_100ks.fak.fak): R06_meg_p1_100ks.fak Exposure time, correction norm (100000., 1.00000): 100000. Fake data file name (R06_heg_n1_100ks.fak.fak): R06_heg_n1_100ks.fak Exposure time, correction norm (100000., 1.00000): 100000. Fake data file name (R06_heg_p1_100ks.fak.fak): R06_heg_p1_100ks.fak Exposure time, correction norm (100000., 1.00000): 100000. No background will be applied to fake spectrum #1 No background will be applied to fake spectrum #2 No background will be applied to fake spectrum #3 No background will be applied to fake spectrum #4 (Overwrote the older four *fak) ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 0.809900 +/- 0.334343 2 1 hegauss G 0.597144 +/- 0.207342 3 1 hegauss sigma_v "km/s" 26.4560 +/- 6.70747E-02 4 1 hegauss delta_v "km/s" -5.52108 +/- 31.3044 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.90983E-06 +/- 2.97003E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 7.00000E-05 frozen ________________________________________________________________________ C-statistic = 151.67 using 356 PHA bins and 351 degrees of freedom. XSPEC12>goodness 1000 nosim 10.80% of realizations are < best fit statistic 151.67 XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) 2 0.482898 0.753502 (-0.115945,0.154660) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 0.416489 0.875463 (-0.182695,0.276279) XSPEC12>error 4. 2 Parameter Confidence Range (4.000000) ***Warning: Identical values of the parameter give different values of the statistic. Please check your result for the high end of the confidence range. 2 0.383134 0.947682 (-0.216362,0.348186) XSPEC12>steppar 6 6.e-6 10.e-6 40 C-Statistic Delta norm C-Statistic 6 157.99172 6.32459 0 0.00001 157.28387 5.61674 1 0.00001 156.62482 4.95769 2 0.00001 156.01505 4.34792 3 0.00001 155.44902 3.78189 4 0.00001 154.92921 3.26208 5 0.00001 154.45224 2.78511 6 0.00001 154.01694 2.34980 7 0.00001 153.62215 1.95502 8 0.00001 153.26679 1.59966 9 0.00001 152.94980 1.28267 10 0.00001 152.67017 1.00304 11 0.00001 152.42692 0.75979 12 0.00001 152.21912 0.55199 13 0.00001 152.04586 0.37873 14 0.00001 151.90629 0.23916 15 0.00001 151.79958 0.13244 16 0.00001 151.72466 0.05753 17 0.00001 151.68061 0.01347 18 0.00001 151.66756 0.00043 19 0.00001 151.68445 0.01731 20 0.00001 151.73051 0.06337 21 0.00001 151.80505 0.13792 22 0.00001 151.90743 0.24030 23 0.00001 152.03704 0.36991 24 0.00001 152.19327 0.52614 25 0.00001 152.37555 0.70842 26 0.00001 152.58331 0.91618 27 0.00001 152.81601 1.14888 28 0.00001 153.07312 1.40599 29 0.00001 153.35414 1.68700 30 0.00001 153.65856 1.99143 31 0.00001 153.98591 2.31878 32 0.00001 154.33574 2.66861 33 0.00001 154.70758 3.04045 34 0.00001 155.10100 3.43387 35 0.00001 155.51558 3.84845 36 0.00001 155.95091 4.28378 37 0.00001 156.40659 4.73945 38 0.00001 156.88222 5.21509 39 0.00001 157.37745 5.71032 40 0.00001 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2009/props/thetaCar/R06_hegauss_100ks.dat Finally, making the R=1.2 simulation with 100 ks exposure: ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 1.20000 frozen 2 1 hegauss G 0.710000 frozen 3 1 hegauss sigma_v "km/s" 0.0 frozen 4 1 hegauss delta_v "km/s" 0.0 frozen 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 7.20000E-06 +/- 0.0 7 2 powerlaw PhoIndex 2.00000 +/- 0.0 8 2 powerlaw norm 7.00000E-05 +/- 0.0 ________________________________________________________________________ XSPEC12>fakeit 4 Use counting statistics in creating fake data? (y): Input optional fake file prefix: Fake data file name (R06_meg_n1_100ks.fak.fak): R12_meg_n1_100ks.fak Exposure time, correction norm (100000., 1.00000): 100000. Fake data file name (R06_meg_p1_100ks.fak.fak): R12_meg_p1_100ks.fak Exposure time, correction norm (100000., 1.00000): 100000. Fake data file name (R06_heg_n1_100ks.fak.fak): R12_heg_n1_100ks.fak Exposure time, correction norm (100000., 1.00000): 100000. Fake data file name (R06_heg_p1_100ks.fak.fak): R12_heg_p1_100ks.fak Exposure time, correction norm (100000., 1.00000): 100000. No background will be applied to fake spectrum #1 No background will be applied to fake spectrum #2 No background will be applied to fake spectrum #3 No background will be applied to fake spectrum #4 ======================================================================== Model hegauss<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hegauss R 1.39442 +/- 0.706745 2 1 hegauss G 0.612716 +/- 0.296251 3 1 hegauss sigma_v "km/s" 21.5021 +/- 9.69979E-02 4 1 hegauss delta_v "km/s" -27.7718 +/- 20.9051 5 1 hegauss Z 12.0000 frozen 6 1 hegauss norm 6.44729E-06 +/- 3.72349E-06 7 2 powerlaw PhoIndex 2.00000 frozen 8 2 powerlaw norm 7.00000E-05 frozen ________________________________________________________________________ C-statistic = 151.62 using 356 PHA bins and 351 degrees of freedom. XSPEC12>goodness 1000 nosim 18.40% of realizations are < best fit statistic 151.62 XSPEC12>error 1. 1 Parameter Confidence Range (1.0000000) 1 0.808537 2.098413 (-0.583131,0.706745) XSPEC12>error 2.7 1 Parameter Confidence Range (2.700000) 1 0.581733 2.821118 (-0.807212,1.432174) XSPEC12>error 4. 1 Parameter Confidence Range (4.000000) 1 0.489866 3.403532 (-0.896400,2.017267) XSPEC12>error 1. 2 Parameter Confidence Range (1.000000) 2 0.477878 0.823690 (-0.135834,0.209977) XSPEC12>error 2.7 2 Parameter Confidence Range (2.700000) 2 0.405540 0.989624 (-0.208243,0.375842) XSPEC12>error 4. 2 Parameter Confidence Range (4.000000) 2 0.368350 1.095895 (-0.245995,0.481550) XSPEC12>steppar 6 5.e-6 9.e-6 40 C-Statistic Delta norm C-Statistic 6 156.37034 4.76521 0 0.00001 155.71193 4.10680 1 0.00001 155.11048 3.50535 2 0.00001 154.56240 2.95727 3 0.00001 154.06797 2.46284 4 0.00001 153.62280 2.01767 5 0.00001 153.22725 1.62212 6 0.00001 152.87859 1.27346 7 0.00001 152.57530 0.97017 8 0.00001 152.31594 0.71081 9 0.00001 152.09912 0.49399 10 0.00001 151.92353 0.31840 11 0.00001 151.78782 0.18269 12 0.00001 151.69050 0.08537 13 0.00001 151.63039 0.02526 14 0.00001 151.60641 0.00128 15 0.00001 151.59833 -0.00680 16 0.00001 151.63553 0.03040 17 0.00001 151.70459 0.09946 18 0.00001 151.77532 0.17019 19 0.00001 151.89992 0.29479 20 0.00001 152.02087 0.41574 21 0.00001 152.19643 0.59130 22 0.00001 152.36154 0.75641 23 0.00001 152.58394 0.97881 24 0.00001 152.78910 1.18397 25 0.00001 153.05690 1.45177 26 0.00001 153.35098 1.74585 27 0.00001 153.67055 2.06542 28 0.00001 154.00006 2.39493 29 0.00001 154.36160 2.75647 30 0.00001 154.74677 3.14164 31 0.00001 155.13707 3.53194 32 0.00001 155.56296 3.95783 33 0.00001 156.01123 4.40610 34 0.00001 156.48137 4.87624 35 0.00001 156.97289 5.36776 36 0.00001 157.48530 5.88017 37 0.00001 158.01816 6.41303 38 0.00001 158.57102 6.96589 39 0.00001 159.03657 7.43144 40 0.00001 XSPEC12>iplot PLT> wdata /home/cohen/xspec/2009/props/thetaCar/R12_hegauss_100ks.dat