Fit a warmabs model to a spectrum in XPEC
In this example we apply warmabs model to the Chandra HETG spectrum of NGC 37843. We use the canned spectrum of observation ID 2090 from TGCat. We use default population files in this example and run XSPEC from the command line:
XSPEC12>data meg_1.pha.gz
XSPEC12>notice 0.9-1.5
XSPEC12>lmod xstarmod \$WARMABS_DATA
XSPEC12>model warmabs*powerlaw
Input parameter value, delta, min, bot, top, and max values for ...
0 0.01( 0.01) -3 -3 2 2
1:warmabs:column>0.35
0 0.231( 0.231) -4 -4 5 5
2:warmabs:rlogxi>2.36
1 -0.01( 0.01) 0 0 1000 1000
3:warmabs:Cabund>0
1 -0.01( 0.01) 0 0 1000 1000
4:warmabs:Nabund>0
1 -0.01( 0.01) 0 0 1000 1000
5:warmabs:Oabund>0
1 -0.01( 0.01) 0 0 1000 1000
6:warmabs:Fabund>0
1 -0.01( 0.01) 0 0 1000 1000
7:warmabs:Neabund>1
1 -0.01( 0.01) 0 0 1000 1000
8:warmabs:Naabund>1
1 -0.01( 0.01) 0 0 1000 1000
9:warmabs:Mgabund>1
1 -0.01( 0.01) 0 0 1000 1000
10:warmabs:Alabund>1
1 -0.01( 0.01) 0 0 1000 1000
11:warmabs:Siabund>1
1 -0.01( 0.01) 0 0 1000 1000
12:warmabs:Pabund>0
1 -0.01( 0.01) 0 0 1000 1000
13:warmabs:Sabund>1
1 -0.01( 0.01) 0 0 1000 1000
14:warmabs:Clabund>0
1 -0.01( 0.01) 0 0 1000 1000
15:warmabs:Arabund>1
1 -0.01( 0.01) 0 0 1000 1000
16:warmabs:Kabund>0
1 -0.01( 0.01) 0 0 1000 1000
17:warmabs:Caabund>1
1 -0.01( 0.01) 0 0 1000 1000
18:warmabs:Scabund>0
1 -0.01( 0.01) 0 0 1000 1000
19:warmabs:Tiabund>0
1 -0.01( 0.01) 0 0 1000 1000
20:warmabs:Vabund>0
1 -0.01( 0.01) 0 0 1000 1000
21:warmabs:Crabund>0
1 -0.01( 0.01) 0 0 1000 1000
22:warmabs:Mnabund>0
1 -0.01( 0.01) 0 0 1000 1000
23:warmabs:Feabund>1
1 -0.01( 0.01) 0 0 1000 1000
24:warmabs:Coabund>0
1 -0.01( 0.01) 0 0 1000 1000
25:warmabs:Niabund>0
1 -0.01( 0.01) 0 0 1000 1000
26:warmabs:Cuabund>0
1 -0.01( 0.01) 0 0 1000 1000
27:warmabs:Znabund>0
0 -0.01( 0.01) 0 0 1 1
28:warmabs:write_outfile>0
0 -0.01( 0.01) 0 0 10000 10000
29:warmabs:outfile_idx>0
0 -0.01( 0.01) 0 0 10000 10000
30:warmabs:vturb>100
0 -0.01( 0.01) 0 0 10 10
31:warmabs:Redshift>9.73000E-03
1 0.01( 0.01) -3 -2 9 10
32:powerlaw:PhoIndex>0.107390
1 0.01( 0.01) 0 0 1e+20 1e+24
33:powerlaw:norm>4.31010E-03
========================================================================
Model warmabs<1>*powerlaw<2> Source No.: 1 Active/On
Model Model Component Parameter Unit Value
par comp
1 1 warmabs column cm^-2 0.350000 +/- 0.0
2 1 warmabs rlogxi 2.36000 +/- 0.0
3 1 warmabs Cabund 0.0 frozen
4 1 warmabs Nabund 0.0 frozen
5 1 warmabs Oabund 0.0 frozen
6 1 warmabs Fabund 0.0 frozen
7 1 warmabs Neabund 1.00000 frozen
8 1 warmabs Naabund 1.00000 frozen
9 1 warmabs Mgabund 1.00000 frozen
10 1 warmabs Alabund 1.00000 frozen
11 1 warmabs Siabund 1.00000 frozen
12 1 warmabs Pabund 0.0 frozen
13 1 warmabs Sabund 1.00000 frozen
14 1 warmabs Clabund 0.0 frozen
15 1 warmabs Arabund 1.00000 frozen
16 1 warmabs Kabund 0.0 frozen
17 1 warmabs Caabund 1.00000 frozen
18 1 warmabs Scabund 0.0 frozen
19 1 warmabs Tiabund 0.0 frozen
20 1 warmabs Vabund 0.0 frozen
21 1 warmabs Crabund 0.0 frozen
22 1 warmabs Mnabund 0.0 frozen
23 1 warmabs Feabund 1.00000 frozen
24 1 warmabs Coabund 0.0 frozen
25 1 warmabs Niabund 0.0 frozen
26 1 warmabs Cuabund 0.0 frozen
27 1 warmabs Znabund 0.0 frozen
28 1 warmabs write_outfile 0.0 frozen
29 1 warmabs outfile_idx 0.0 frozen
30 1 warmabs vturb km/s 100.000 frozen
31 1 warmabs Redshift 9.73000E-03 frozen
32 2 powerlaw PhoIndex 0.107390 +/- 0.0
33 2 powerlaw norm 4.31010E-03 +/- 0.0
________________________________________________________________________
doing setup: 0 1 0 0 0
photemis/warmabs v2.47
xstar version 2.58e
Loading Atomic Database...
Loading Atomic Database...
Atomic Data Version: 2021-08-03T18:36:21
number of pointers= 1216828
number of reals= 199199716
number of integers= 6205342
number of chars= 753844
initializng database...
initializng database...
number of lines= 736256
number of rrcs= 301301
max levels per ion: 744
habund 0.1000E+01 0.1000E+01
heabund 0.1000E+01 0.1000E+00
liabund 0.0000E+00 0.0000E+00
beabund 0.0000E+00 0.0000E+00
babund 0.0000E+00 0.0000E+00
cabund 0.1000E+01 0.3700E-03
nabund 0.1000E+01 0.1100E-03
oabund 0.1000E+01 0.6800E-03
fabund 0.1000E+01 0.3980E-07
neabund 0.1000E+01 0.2800E-04
naabund 0.1000E+01 0.1780E-05
mgabund 0.1000E+01 0.3500E-04
alabund 0.1000E+01 0.2450E-05
siabund 0.1000E+01 0.3500E-04
pabund 0.1000E+01 0.3310E-06
sabund 0.1000E+01 0.1600E-04
clabund 0.1000E+01 0.3980E-06
arabund 0.1000E+01 0.4500E-05
kabund 0.1000E+01 0.8910E-07
caabund 0.1000E+01 0.2100E-05
scabund 0.1000E+01 0.1660E-08
tiabund 0.1000E+01 0.1350E-06
vabund 0.1000E+01 0.2510E-07
crabund 0.1000E+01 0.7080E-06
mnabund 0.1000E+01 0.2510E-06
feabund 0.1000E+01 0.2500E-04
coabund 0.1000E+01 0.1260E-06
niabund 0.1000E+01 0.2000E-05
cuabund 0.1000E+01 0.3160E-07
znabund 0.1000E+01 0.1580E-07
Reading pop files
using files xout_detail.fits etc. 219
Fit statistic : Chi-Squared 29732.67 using 8192 bins.
Test statistic : Chi-Squared 29732.67 using 8192 bins.
Null hypothesis probability of 0.00e+00 with 8188 degrees of freedom
Current data and model not fit yet.
XSPEC12>plot
