PRO axaf1_versions

@effic_com
@axaf1_com

plot_em = 0
plot_superpos = 1

ax_sname = "FLAT 1 count/cm^2 sec bin"
ax_nmps = 4
ax_coating = ["Ni","Ni","Ir","Ir"]
ax_mps = [1,3,4,6]
ax_arcsec = [0,0,0,0]
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.1
ax_dlength = 136.9
ax_plots = [0,0,1,0,1]

spect_all = fltarr(ax_ne,12)
isa = 0

;MEG 		HEG
if(plot_em EQ 1) then pre_print
ax_gmat = ["Ag","Ag","Au","Au"]
ax_gthk = [0.5,0.5,1.0,1.0]
ax_gp = 1000.*[4.5,4.5,2.0,2.0]
axaf1
if(plot_em EQ 1) then post_print
spect_all(*,isa) = effic_total
isa = isa+1

;Thick MEG 	HEG
if(plot_em EQ 1) then pre_print
ax_gmat = ["Ag","Ag","Au","Au"]
ax_gthk = [0.5,1.0,1.0,1.0]
ax_gp = 1000.*[4.5,4.5,2.0,2.0]
axaf1
if(plot_em EQ 1) then post_print
spect_all(*,isa) = effic_total
isa = isa+1

;Thin HEG 	HEG
if(plot_em EQ 1) then pre_print
ax_gmat = ["Ag","Au","Au","Au"]
ax_gthk = [0.5,0.6,1.0,1.0]
ax_gp = 1000.*[4.5,2.0,2.0,2.0]
axaf1
if(plot_em EQ 1) then post_print
spect_all(*,isa) = effic_total
isa = isa+1

;HEG 		HEG
if(plot_em EQ 1) then pre_print
ax_gmat = ["Ag","Au","Au","Au"]
ax_gthk = [0.5,1.0,1.0,1.0]
ax_gp = 1000.*[4.5,2.0,2.0,2.0]
axaf1
if(plot_em EQ 1) then post_print
spect_all(*,isa) = effic_total
isa = isa+1

;MEG 		Thin HEG
if(plot_em EQ 1) then pre_print
ax_gmat = ["Ag","Ag","Au","Au"]
ax_gthk = [0.5,0.5,0.6,1.0]
ax_gp = 1000.*[4.5,4.5,2.0,2.0]
axaf1
if(plot_em EQ 1) then post_print
spect_all(*,isa) = effic_total
isa = isa+1

;Thick MEG 	Thin HEG
if(plot_em EQ 1) then pre_print
ax_gmat = ["Ag","Ag","Au","Au"]
ax_gthk = [0.5,1.0,0.6,1.0]
ax_gp = 1000.*[4.5,4.5,2.0,2.0]
axaf1
if(plot_em EQ 1) then post_print
spect_all(*,isa) = effic_total
isa = isa+1

;Thin HEG 	Thin HEG
if(plot_em EQ 1) then pre_print
ax_gmat = ["Ag","Au","Au","Au"]
ax_gthk = [0.5,0.6,0.6,1.0]
ax_gp = 1000.*[4.5,2.0,2.0,2.0]
axaf1
if(plot_em EQ 1) then post_print
spect_all(*,isa) = effic_total
isa = isa+1

;HEG 		Thin HEG
if(plot_em EQ 1) then pre_print
ax_gmat = ["Ag","Au","Au","Au"]
ax_gthk = [0.5,1.0,0.6,1.0]
ax_gp = 1000.*[4.5,2.0,2.0,2.0]
axaf1
if(plot_em EQ 1) then post_print
spect_all(*,isa) = effic_total
isa = isa+1

;MEG 		Thick MEG
if(plot_em EQ 1) then pre_print
ax_gmat = ["Ag","Ag","Ag","Au"]
ax_gthk = [0.5,0.5,1.0,1.0]
ax_gp = 1000.*[4.5,4.5,4.5,2.0]
axaf1
if(plot_em EQ 1) then post_print
spect_all(*,isa) = effic_total
isa = isa+1

;Thick MEG 	Thick MEG
if(plot_em EQ 1) then pre_print
ax_gmat = ["Ag","Ag","Ag","Au"]
ax_gthk = [0.5,1.0,1.0,1.0]
ax_gp = 1000.*[4.5,4.5,4.5,2.0]
axaf1
if(plot_em EQ 1) then post_print
spect_all(*,isa) = effic_total
isa = isa+1

;Thin HEG 	Thick MEG
if(plot_em EQ 1) then pre_print
ax_gmat = ["Ag","Au","Ag","Au"]
ax_gthk = [0.5,0.6,1.0,1.0]
ax_gp = 1000.*[4.5,2.0,4.5,2.0]
axaf1
if(plot_em EQ 1) then post_print
spect_all(*,isa) = effic_total
isa = isa+1

;HEG 		Thick MEG
if(plot_em EQ 1) then pre_print
ax_gmat = ["Ag","Au","Ag","Au"]
ax_gthk = [0.5,1.0,1.0,1.0]
ax_gp = 1000.*[4.5,2.0,4.5,2.0]
axaf1
if(plot_em EQ 1) then post_print
spect_all(*,isa) = effic_total
isa = isa+1

; Plot all of the curves on one plot!
IF(ax_plots(4) GE 1) THEN BEGIN
	!MTITLE = "HEG(dot) / MEG(dash) Comparison. Spectrum is "+ax_sname
	plot_oo, ax_es, total_spec, $
	XRANGE = [0.1,10.], YRANGE = [1.0, 200.0], LINESTYLE=0, $
	XTITLE = "Energy (keV)", YTITLE = "counts / sec bin"
	oplot, ax_es, heg_spec, LINESTYLE = 1
	oplot, ax_es, meg_spec, LINESTYLE = 2
; label the plot w/axaf1 info
axaf1_label,3,0.,0.57
END



RETURN
END