PRO JAN20overlay, PS=ps, TIMES=times, BEAM=beam, REVERSE=reverse, ELABELS=elabels
;
; Prepares an overlay plot of TRACE, 250-500 keV contours and 2.2 MeV circles.
;
; 13-May-05		Initial version, adapted from NOV02overlay.pro
; 17-May-05	gh	Incorporated TRACE pointing corrections
; 28-Nov-05 gh	Updated to reflect November 26 mapping results and optimized for poster.
; 29-Nov-05 gh	Added /REVERSE switch to show TRACE image as negative
; 1-Dec-05 gh	Restore 250-500 keV centroid error circle.
;  2-Dec-05 gh	Add /ELABELS keyword to insert date and energies
;
; beam = FWHM beam diameter to be overlaid
;
; See Sam Krucker's emails of 3 Dec 2003 and 25 May 2004.
;
cs = 1.4		; 1.0 or 1.4, depending on how plot looks
DEFAULT, elabels, 0
DEFAULT, notitle 0
IF ELABELS EQ 0 THEN notitle = 1
;
IF KEYWORD_SET(ps) THEN BEGIN
	SET_PLOT, 'ps'
	DEVICE, /COLOR, BITS_PER_PIXEL=8
ENDIF
;
; Plot TRACE map
fits2map, 'c:\HESSI_DATA\SYNOPTIC\trac_1600___a2_20050120_065230.fts', tmap
tmap.xc=tmap.xc-9.					; Correct TRACE pointing per Sam Krucker's 5/15/05 email
tmap.yc=tmap.yc+9.5
LOADCT, 8							; GREEN/WHITE LINEAR
IF KEYWORD_SET(reverse) THEN tmap.data = -tmap.data				; reverse to make a white background
PLOT_MAP, tmap, FOV=[3.,3.], NOTITLE=notitle
;
; Overplot 250-500  electron map, with footpoint centroids marked.
fits2map, 'JAN20-Q.fits', xmap
maxx	= MAX(xmap.data)
LOADCT, 6							; PRISM
PLOT_MAP, xmap, /OVER, /CONT, LEVELS=[0.3,.5,.7,.9]*maxx, LCOLOR=64, CTHICK=3
OPLOT, [821.4, 816.7], [239.4, 272.4], PSYM=1, THICK=3, COLOR=64					; peak locations
OPLOT_C, 818.4, 249.2, 1.3, 1.3, COLOR=64, THICK=7									; centroid circle
;
; Overplot 2.2MeV centroid locations
OPLOT_C, 806.9, 267.5,  4.1,  4.1, COLOR=192, THICK=7
IF elabels NE 0 THEN BEGIN
	XYOUTS, 795., 345., ' 250-500  keV', COLOR=64,  CHARTHICK=2, CHARSIZE=1.4
	XYOUTS, 795., 335., '2215-2231 keV', COLOR=192, CHARTHICK=2, CHARSIZE=1.4
ENDIF
;
;
OPLOT_C, 750,200, 17.5,17.5, COLOR=192
LOADCT, 8
IF KEYWORD_SET(reverse) THEN beamcolor= 0 ELSE beamcolor=256
XYOUTS, 735, 201, '2.2MeV', COLOR=beamcolor, CHARSIZE=1.2
XYOUTS, 735, 191, ' BEAM',  COLOR=beamcolor, CHARSIZE=1.2
;OPLOT_C, 750, 200, 91.5, 91.5, COLOR=192, LINESTYLE=1
IF KEYWORD_SET(ps) THEN BEGIN
	DEVICE,/CLOSE
	SET_PLOT, 'win'
ENDIF
END