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