========================== Some example SPEX sessions ========================== This document currently contains four examples: 1) An AGN with a simple power law spectrum 2) An AGN with a warm absorber 3) An AGN with a warm absorber and relativistically-broadened accretion disc emission lines. 4) An AGN with a two phase warm absorber (e.g. NGC 3783) and three narrow emission lines. ----- 1) An AGN with a simple power law spectrum Assume you have a spectrum ex_1.spo and a response matrix ex_1.res. The following set of commands will read in the data, plot it, construct a model, then calculate and plot the model on the data. Model component 3 serves to take account of the RGS instrumental oxygen feature. data ex_1 ex_1 #read in data and response distance 58.2 mpc #set distance plot type data #plot data plot ux a plot uy a plot device xw plot comp pow #set up model comp absm comp slab comp relation 1 2,3 model show #show model par 1 2 nh value 8.7e-4 #set galactic hydrogen column par 1 2 nh status f par 1 1 gam value 1.655 #set power law par 1 1 gam status f par 1 1 norm value 52872244.72 par 1 1 norm status f par 1 3 o1 value 21.53 #instrumental oxygen feature par 1 3 o1 status f par show #show model parameters calculate #calculate model plot #plot data and model --- 2) An AGN with a warm absorber These commands follow the form above, but this time the model is of an AGN with a two-phase outflowing warm absorber. The hotter phase contains ions such as O VII, N VII, C VI and L-shell iron and is blueshifted at 500 km/s. The cooler phase, which is not outflowing, contains an Unresolved Transition Array (UTA) of M-shell iron. Again, the RGS instrumental oxygen feature is taken account of. data ex_2 ex_2 #read in data and response distance 58.2 mpc #set distance plot type data #plot data plot ux a plot uy a plot device xw plot comp pow #set up model comp absm comp reds comp slab comp slab comp slab comp relation 1 4,6,3,2,5 par 1 2 nh value 8.7e-4 #set galactic hydrogen column par 1 2 nh status f par 1 1 gam value 1.655 #set power law par 1 1 gam status f par 1 1 norm value 52872244.72 par 1 1 norm status f par 1 3 z value 0.0097 #set cosmological redshift par 1 3 z status f par 1 4 o8 val 20.72996 #oxygen in hot phase par 1 4 o8 sta f par 1 4 o7 val 21.58468 par 1 4 o7 sta f par 1 4 o6 val 19.95547 par 1 4 o6 sta f par 1 4 o1 val 1 par 1 4 o1 sta f par 1 4 n7 val 20.64288 #nitrogen in hot phase par 1 4 n7 sta f par 1 4 n6 val 20.54342 par 1 4 n6 sta f par 1 4 n5 val 1 par 1 4 n5 sta f par 1 4 n1 val 1 par 1 4 n1 sta f par 1 4 c6 val 20.94342 #carbon in hot phase par 1 4 c6 sta f par 1 4 c5 val 21.88750 par 1 4 c5 sta f par 1 4 c4 val 1 par 1 4 c4 sta f par 1 4 c1 val 1 par 1 4 c1 sta f par 1 4 fe17 val 21.71937 #iron in hot phase par 1 4 fe17 sta f par 1 4 fe18 val 21.72341 par 1 4 fe18 sta f par 1 4 fe19 val 21.53607 par 1 4 fe19 sta f par 1 4 fe20 val 21.61097 par 1 4 fe20 sta f par 1 4 fe21 val 21.79189 par 1 4 fe21 sta f par 1 4 fe22 val 20.28482 par 1 4 fe22 sta f par 1 4 fe23 val 20.79127 par 1 4 fe23 sta f par 1 4 fe24 val 12.90610 par 1 4 fe24 sta f par 1 4 v val 100 #set velocity broadening par 1 4 v sta f par 1 4 zv val -500 #blueshift of hot phase par 1 4 zv sta f par 1 5 o1 val 21.53 #instrumental oxygen edge par 1 5 o1 sta f par 1 6 fe07 val 20.86841 #UTA in cold phase par 1 6 fe07 sta f par 1 6 fe08 val 21.00819 par 1 6 fe08 sta f par 1 6 fe09 val 20.44001 par 1 6 fe09 sta f par 1 6 fe10 val 20.94657 par 1 6 fe10 sta f par 1 6 fe11 val 20.70722 par 1 6 fe11 sta f par 1 6 fe12 val 20.14890 par 1 6 fe12 sta f par 1 6 fe14 val 20.40249 par 1 6 fe14 sta f par 1 6 v val 100 #velocity broadening of cold phase par 1 6 v sta f par 1 6 zv val 0 #blueshift of cold phase par 1 6 zv sta f calculate #calculate model plot #plot model --- 3) An AGN with a warm absorber and relativistically-broadened accretion disc emission lines. This example uses the same structure and warm absorber model as in the previous example, but this time there are relativistically broadened oxygen and nitrogen emission lines as well. Note that the Laor model is multiplicative: in this case it is convolved with a delta function. data ex_3 ex_3 #read in data and response distance 58.2 mpc #set distance plot type data #plot data plot ux a plot uy a plot device xw plot comp pow #set up model comp absm comp reds comp slab comp slab comp slab comp delt comp delt comp laor comp laor comp relation 1 4,6,3,2,5 comp relation 7 9,4,6,3,2,5 comp relation 8 10,4,6,3,2,5 par 1 2 nh value 8.7e-4 #set galactic hydrogen column par 1 2 nh status f par 1 1 gam value 1.655 #set power law par 1 1 gam status f par 1 1 norm value 52872244.72 par 1 1 norm status f par 1 3 z value 0.0097 #set cosmological redshift par 1 3 z status f par 1 4 o8 val 20.72996 #oxygen in hot phase par 1 4 o8 sta f par 1 4 o7 val 21.58468 par 1 4 o7 sta f par 1 4 o6 val 19.95547 par 1 4 o6 sta f par 1 4 o1 val 1 par 1 4 o1 sta f par 1 4 n7 val 20.64288 #nitrogen in hot phase par 1 4 n7 sta f par 1 4 n6 val 20.54342 par 1 4 n6 sta f par 1 4 n5 val 1 par 1 4 n5 sta f par 1 4 n1 val 1 par 1 4 n1 sta f par 1 4 c6 val 20.94342 #carbon in hot phase par 1 4 c6 sta f par 1 4 c5 val 21.88750 par 1 4 c5 sta f par 1 4 c4 val 1 par 1 4 c4 sta f par 1 4 c1 val 1 par 1 4 c1 sta f par 1 4 fe17 val 21.71937 #iron in hot phase par 1 4 fe17 sta f par 1 4 fe18 val 21.72341 par 1 4 fe18 sta f par 1 4 fe19 val 21.53607 par 1 4 fe19 sta f par 1 4 fe20 val 21.61097 par 1 4 fe20 sta f par 1 4 fe21 val 21.79189 par 1 4 fe21 sta f par 1 4 fe22 val 20.28482 par 1 4 fe22 sta f par 1 4 fe23 val 20.79127 par 1 4 fe23 sta f par 1 4 fe24 val 12.90610 par 1 4 fe24 sta f par 1 4 v val 100 #velocity broadening in hot phase par 1 4 v sta f par 1 4 zv val -500 #blueshift of hot phase par 1 4 zv sta f par 1 5 o1 val 21.53 #instrumental oxygen feature par 1 5 o1 sta f par 1 6 fe07 val 20.86841 #UTA in cold phase par 1 6 fe07 sta f par 1 6 fe08 val 21.00819 par 1 6 fe08 sta f par 1 6 fe09 val 20.44001 par 1 6 fe09 sta f par 1 6 fe10 val 20.94657 par 1 6 fe10 sta f par 1 6 fe11 val 20.70722 par 1 6 fe11 sta f par 1 6 fe12 val 20.14890 par 1 6 fe12 sta f par 1 6 fe14 val 20.40249 par 1 6 fe14 sta f par 1 6 v val 100 #velocity broadening of cold phase par 1 6 v sta f par 1 6 zv val 0 #blueshift of cold phase par 1 6 zv sta f par 1 7 e val 0.654 #delta function for O VIII line par 1 7 e sta f par 1 7 norm val 1.5e7 par 1 7 norm sta f par 1 8 e val 0.500 #delta function for N VII line par 1 8 e sta f par 1 8 norm val 1.5e7 par 1 8 norm sta f par 1 9 q val 4 #relativistic effect on O VIII line par 1 9 q sta f par 1 9 i val 40 par 1 9 i sta f par 1 10 q val 4 #relativistic effect on N VII line par 1 10 q sta f par 1 10 i val 40 par 1 10 i sta f calculate #calculate model plot #plot data and model --- 4) NGC 3783 2 zone warm absorber model with three narrow emission lines data ngc3783_v52_nov ngc3783_v52_nov distance 0.00973 z plot type data plot ux a plot uy fa plot device xw p rx 6 38 p ry 0 15 comp pow comp absm comp reds comp xabs comp xabs comp gaus comp gaus comp gaus comp relation 1 4 , 5 , 3 , 2 comp relation 6 4 , 5 , 3 , 2 comp relation 7 4 , 5 , 3 , 2 comp relation 8 4 , 5 , 3 , 2 par 1 2 nh value 8.7e-4 par 1 2 nh status f par 1 1 gam value 1.54 par 1 1 gam status f par 1 1 norm value 71483945.64 par 1 1 norm status f par 1 3 z value 0.00973 par 1 3 z status f par 1 4 nh val 2.8429955E-02 par 1 4 nh sta f par 1 4 v val 300 par 1 4 zv val -1000 par 1 4 xil val 2.4 par 1 4 xil sta f par 1 4 02 val 1 par 1 4 06 val 0.8 par 1 4 07 val 0.2 par 1 4 08 val 0.7453385 par 1 4 10 val 0.9365730 par 1 4 11 val 1 par 1 4 12 val 0.8689545 par 1 4 14 val 2.5 par 1 4 16 val 1 par 1 4 18 val 0.7 par 1 4 26 val 0.9424475 par 1 5 nh val 5.3898338E-04 par 1 5 nh sta f par 1 5 v val 300 par 1 5 zv val -1000 par 1 5 xil val 0.3 par 1 5 xil sta f par 1 5 02 val 1 par 1 5 06 val 1 par 1 5 07 val 1 par 1 5 08 val 1 par 1 5 10 val 1 par 1 5 11 val 1 par 1 5 12 val 1 par 1 5 14 val 1 par 1 5 16 val 1 par 1 5 18 val 1 par 1 5 26 val 10.7 par 1 6 e val 0.5757136 par 1 6 e sta f par 1 6 fwhm val 0 par 1 6 fwhm sta f par 1 6 norm val 1.7894856E+08 par 1 6 norm sta f par 1 7 e val 0.5703016 par 1 7 e sta f par 1 7 fwhm val 0 par 1 7 fwhm sta f par 1 7 norm val 684000.0 par 1 7 norm sta f par 1 8 e val 0.5627100 par 1 8 e sta f par 1 8 fwhm val 1.8611701E-03 par 1 8 fwhm sta f par 1 8 norm val 1064890. par 1 8 norm sta f calculate plot ------------------