Randall Smith
sherpa
-----------------------------------------------------
Welcome to Sherpa: CXC's Modeling and Fitting Program
-----------------------------------------------------
Version: 2.0 (2 June 2000)
Type HELP for help options.
Type EXIT, QUIT, or BYE to leave the program.
Notes:
Temporary files for visualization will be written to the directory:
/tmp
To change this so that these files are not deleted when you exit Sherpa,
edit the variable ASCDS_WORK_PATH in your $HOME/.cxcds.*sh setup script.
Solar abundances set to Anders & Grevesse
sherpa> capella = readfile("acisf01318N002_pha2.fits")
Warning: could not find SYS_ERR column
sherpa> print(capella)
filename = acisf01318N002_pha2.fits
path = /data/rsmith2/capella/ObsID1318/tertiary/demo/
filter = NULL
ncols = 11
nrows = 12
bgscale = 1
areascale = 1
exptime = 26699.7
channels = Float_Type[98304]
counts = Float_Type[98304]
grouping = Integer_Type[12,8192]
binhi = Float_Type[12,8192]
binlo = Float_Type[12,8192]
errors = Float_Type[98304]
background = String_Type[12]
arf = String_Type[12]
response = String_Type[12]
numgroups = Integer_Type[12]
numchans = Integer_Type[12]
sherpa> data capella
sherpa> analysis
Analysis Space for Dataset 1: Channel
Analysis Space for Dataset 2: Channel
Analysis Space for Dataset 3: Channel
Analysis Space for Dataset 4: Channel
Analysis Space for Dataset 5: Channel
Analysis Space for Dataset 6: Channel
Analysis Space for Dataset 7: Channel
Analysis Space for Dataset 8: Channel
Analysis Space for Dataset 9: Channel
Analysis Space for Dataset 10: Channel
Analysis Space for Dataset 11: Channel
Analysis Space for Dataset 12: Channel
sherpa> analysis wave
sherpa> rsp[hegm1]
sherpa> rsp[hegp1]
sherpa> rsp[megm1]
sherpa> rsp[megp1]
sherpa> hegm1.arf = heg_m1.arf
The inferred file type is ARF. If this is not what you want, please
specify the type explicitly in the data command.
sherpa> hegp1.arf = heg_p1.arf
The inferred file type is ARF. If this is not what you want, please
specify the type explicitly in the data command.
sherpa> megm1.arf = meg_m1.arf
The inferred file type is ARF. If this is not what you want, please
specify the type explicitly in the data command.
sherpa> megp1.arf = meg_p1.arf
The inferred file type is ARF. If this is not what you want, please
specify the type explicitly in the data command.
sherpa> instrument 3 = hegm1
sherpa> instrument 4 = hegp1
sherpa> instrument 9 = megm1
sherpa> instrument 10= megp1
sherpa> import("guide")
GUIDE Initialized
WARNING: GUIDE is BETA software.
WARNING: USE AT YOUR OWN RISK.
sherpa> set plot noerrorbars
sherpa> lp 4 data 3 data 4 data 9 data 10
sherpa> ignore allsets all
sherpa> notice allsets wave 8.2:8.6
sherpa> lp 4 data 3 data 4 data 9 data 10
sherpa> source 3,4 = ngauss[hg1] + poly[hp1]
sherpa> source 9,10 = ngauss[mg1] + poly[mp1]
sherpa> show source 3
Source 3: (hg1 + hp1)
normgauss1d[hg1] (integrate: on)
Param Type Value Min Max Units
----- ---- ----- --- --- -----
1 fwhm thawed 0.0242 0.0002 2.4239
2 pos thawed 8.4217 8.1992 8.5992
3 ampl thawed 0.0003 3.2096e-06 0.0321
polynom1d[hp1] (integrate: on)
Param Type Value Min Max Units
----- ---- ----- --- --- -----
1 c0 thawed 0.0062 -0 0.0124
2 c1 frozen 0 -3.1097 3.1097
3 c2 frozen 0 -7.7738 7.7738
4 c3 frozen 0 -0 0.0124
5 c4 frozen 0 -0 0.0124
6 c5 frozen 0 -0 0.0124
7 c6 frozen 0 -0 0.0124
8 c7 frozen 0 -0 0.0124
9 c8 frozen 0 -0 0.0124
10 offset frozen 0 -8.1992 8.5992
sherpa> hg1.fwhm = 0.012
sherpa> freeze hg1.fwhm
sherpa> mg1.fwhm = 0.023
sherpa> freeze mg1.fwhm
sherpa> mg1.ampl => hg1.ampl
sherpa> fit 3,4,9,10
powll: v1.2
powll: initial function value = 3.46138E+04
powll: converged to minimum = 1.78342E+02 at iteration = 7
powll: final function value = 1.78342E+02
hg1.pos 8.42229
hg1.ampl 0.000161188
hp1.c0 0.000334997
mg1.pos 8.41981
mp1.c0 0.000382111
sherpa> identify(8.42)
Found 8 lines.
Lambda -- Ion UpperLev LowerLev Emis(ph cm^3/s) @ Peak Temp
8.4084 -- Fe XXII 158 -> 7, 1.427e-20 @ logT = 7.10
8.4144 -- Fe XXII 157 -> 6, 8.474e-20 @ logT = 7.10
8.4157 -- Mg XII 1014 -> 15, 1.083e-20 @ logT = 6.90
8.4160 -- Mg XII 1055 -> 13, 5.761e-19 @ logT = 6.90
8.4181 -- Mg XII 1022 -> 9, 4.523e-19 @ logT = 6.90
8.4193 -- Mg XII 4 -> 1, 4.770e-17 @ logT = 6.90
8.4209 -- Fe XXII 177 -> 8, 1.283e-18 @ logT = 7.10
8.4247 -- Mg XII 3 -> 1, 2.266e-17 @ logT = 6.90
sherpa> describe(12,12,4,1)
Ion Mg XII, energy level 1:
electron configuration: 1S1
: 1s 2S0.5
1
energy above ground state (eV): 0.000000
Quantum state: n=1, l=0, s=2.000000, degeneracy=2.000000
Data source: Sampson,Goett, & Clark,1983,ADNDT,29,467
-------------------------------------------
Ion Mg XII, energy level 4:
electron configuration: 2P1
: 2p 2P1.5
1
energy above ground state (eV): 1469.430054
Quantum state: n=2, l=1, s=2.000000, degeneracy=4.000000
Data source: Sampson,Goett, & Clark,1983,ADNDT,29,467
-------------------------------------------
Interactions between levels 1 and 4 in Mg XII:
Electron collisional transitions between 1 and 4 exist.
Reference is P.Desai:Merged Sampson,Kiselius
No file with Z=12 and type=5 in filemap.
read_pcoll: No proton collision data file for Mg XII in file map.
Radiative transitions between 1 and 4 exist
Observed wavelength is 8.419300 Angstroms
Theoretical wavelength is 8.438330 Angstroms
Transition rate (A value) is 6.391090e+12 s^-1
Reference is Fernley,Taylor,Seaton1987,JPhysB,20,6457
sherpa> hg1 lineid "APECline(12,12,4,1)+APECline(12,12,3,1)"
sherpa> hg1 filter "ignore allsets all; notice allsets wave 8.2:8.6"
sherpa> write mdl "MgXII_MDL.fits"
WARNING: Dataset 1 will not be included in the fit.
WARNING: Dataset 6 will not be included in the fit.
sherpa> prism MgXII_MDL.fits