; hsi_brad_qe.pro
;
; Compare the old 'N19961106 HSI QE with
; Brad's newest...
;

; Current CIP HSI QE:
cip_hsi = rdb_read(!DDHXDSDIR+'/hsi_qe_N19961106.rdb')

; Brad's new QE
brad_hsi = rdb_read(!DDHXDSDIR+'/hsi_qe_N19990530.rdb')

; and UVIS transmission
uvis = rdb_read(!DDHXDSDIR+'/hsi_uvis_N19980529.rdb')
; convert Energy to keV
uvis.energy = uvis.energy/1000.0

; We'll plot the raw the MCP QEs so need the UVIS interpolated to
; the same energy grids...
brad_uvis = INTERPOL(uvis.trans, uvis.energy, brad_hsi.energy)
cip_uvis = INTERPOL(uvis.trans, uvis.energy, cip_hsi.energy)


; Make two sets of three-panel plots:
; - one linear-Energy from 0.0 to 0.6 keV
; - one log-Energy from 0.06 to 10.0 keV

pre_print_portrait
!p.multi=[0,1,3]

pltxrange = [0.05,0.6]

plot, brad_hsi.energy, brad_hsi.QE, PSYM=-6, LINESTYLE=2, $
	XRANGE=pltxrange, XSTYLE=1, $
	XTITLE="Energy (keV)", YTITLE="HSI QE", $
	YRANGE=[0.0, 0.30], CHARSIZE=1.8, $
	TITLE="HSI QE: 'N19961106 - dotted; 'N19990530 - dashed"
oplot, cip_hsi.energy, cip_hsi.QE, PSYM=-4, LINESTYLE=1

plot, uvis.energy, uvis.trans, $
	XRANGE=pltxrange, XSTYLE=1, CHARSIZE=1.8, $
	XTITLE="Energy (keV)", YTITLE="UVIS Transmission", $
	TITLE="HSI UVIS Transmission, 'N19980529"
oplot, cip_hsi.energy, cip_uvis, PSYM=4
oplot, brad_hsi.energy, brad_uvis, PSYM=6


plot, brad_hsi.energy, brad_hsi.QE/brad_uvis, PSYM=-6, LINESTYLE=2, $
	XRANGE=pltxrange, XSTYLE=1, $
	XTITLE="Energy (keV)", YTITLE="Raw-MCP QE", $
	CHARSIZE=1.8, YRANGE=[0.0,0.35], YSTYLE=1, $
	TITLE="HSI raw-MCP QE: 'N19961106 - dotted; 'N19990530 - dashed"
oplot, cip_hsi.energy, cip_hsi.QE/cip_uvis, PSYM=-4, LINESTYLE=1

device,/close
$cp idl.ps hsi_new_qe_lowE.ps

pre_print_portrait
!p.multi=[0,1,3]

pltxrange = [0.06,10.0]

plot_oi, brad_hsi.energy, brad_hsi.QE, PSYM=-6, LINESTYLE=2, $
	XRANGE=pltxrange, XSTYLE=1, $
	XTITLE="Energy (keV)", YTITLE="HSI QE", $
	YRANGE=[0.0, 0.30], CHARSIZE=1.8, $
	TITLE="HSI QE: 'N19961106 - dotted; 'N19990530 - dashed"
oplot, cip_hsi.energy, cip_hsi.QE, PSYM=-4, LINESTYLE=1

plot_oi, uvis.energy, uvis.trans, $
	XRANGE=pltxrange, XSTYLE=1, CHARSIZE=1.8, $
	XTITLE="Energy (keV)", YTITLE="UVIS Transmission", $
	TITLE="HSI UVIS Transmission, 'N19980529"
oplot, cip_hsi.energy, cip_uvis, PSYM=4
oplot, brad_hsi.energy, brad_uvis, PSYM=6


plot_oi, brad_hsi.energy, brad_hsi.QE/brad_uvis, PSYM=-6, LINESTYLE=2, $
	XRANGE=pltxrange, XSTYLE=1, $
	XTITLE="Energy (keV)", YTITLE="Raw-MCP QE", $
	CHARSIZE=1.8, YRANGE=[0.0,0.35], YSTYLE=1, $
	TITLE="HSI raw-MCP QE: 'N19961106 - dotted; 'N19990530 - dashed"
oplot, cip_hsi.energy, cip_hsi.QE/cip_uvis, PSYM=-4, LINESTYLE=1

device,/close
$cp idl.ps hsi_new_qe.ps

set_plot, 'X'
!p.multi=0