pro qe_1000nm

   !P.CHARSIZE=8
    !P.CHARTHICK=10.
    !X.THICK=20.
    !Y.THICK=20.
    !P.THICK=80.
    symsize=10.
    filenamecharsize=8.
    !P.TICKLEN=0.03

    !X.MINOR=2
    ; define forground and background colors
    white=255
    black=0
	yellow=120

	txtfile='qe_1000nm.txt'

 	readcol,txtfile,temp,QE
 	set_plot, 'z'
    device, set_resolution=[4800,3200]
    device, set_font='Courier'
    temps=temp(sort(temp))
    QEs=QE(sort(temp))
    plot,temps,QEs,position=[0.15,0.15,0.8,0.95],xrange=[160,240],$
    	psym=2,xtitle='Temperature (K)',ystyle=8,$
         ytitle='QE (normalized fraction)',title='QE vs. Temperature',$
         yrange=[0,max(QE)+0.1],background=white,color=yellow

    jpgimg=tvrd()
    write_jpeg,'qe_1000nm.jpg', congrid(jpgimg, 1280, 1024, /interp, /center), quality=100

end