; Show HEG vs MEG resolving power
; by comparing MEG on 3 with HEG on 4 looking 
; at Cas_A_0.9_1.2_HEG (MEG) spectra
;
; Type @middle_pairs to run this
;
; All four baseline:
ax_sname = "FLAT 1 count/cm^2 sec bin"
;
ax_nmps = 4
ax_mps = [1,3,4,6]
ax_coating = ["Ir","Ir","Ir","Ir"]
ax_arcsec = [1,1,1,1]
;
ax_gmat = ["Ag","Ag","Au","Au"]
ax_gthk = [0.5,0.5,1.0,1.0]
ax_gp = [4.6,4.6,2.0,2.0]*1000.
ax_gpoly = [1.0,1.0,1.0,1.0]
ax_go = [0.88,0.88,0.85,0.85]
;
ax_detector = "ACIS - Frontside"
ax_doff = 3.3
ax_dlength = 6.*0.018*1280.
;
ax_plots = [0,0,0,0,1]



; 4 with HEG
ax_gmat = ["Ag","Ag","Au","Au"]
ax_gthk = [0.5,0.5,1.0,1.0]
ax_gp = [4.6,4.6,2.0,2.0]*1000.
ax_gpoly = [1.0,1.0,1.0,1.0]
ax_go = [0.00,0.00,0.85,0.00]

ax_sname = 'file: s_Cas_A_0.9_1.2_HEG'
axaf1
axaf1_observ, 's_Cas_A_0.9_1.2_HEG.ax1'

; 3 with MEG
ax_gmat = ["Ag","Ag","Au","Au"]
ax_gthk = [0.5,0.5,1.0,1.0]
ax_gp = [4.6,4.6,2.0,2.0]*1000.
ax_gpoly = [1.0,1.0,1.0,1.0]
ax_go = [0.00,0.88,0.00,0.00]

ax_sname = 'file: s_Cas_A_0.9_1.2_MEG'
axaf1
axaf1_observ, 's_Cas_A_0.9_1.2_MEG.ax1'