PRO JAN20_LC_COMPOSITE, BFACTOR=bfactor, DECX=decx
;
; Quickie program to generate composite Jan 20 light curve for AGU05
; BFACTOR = nc background normalization factor
;
; 30-Nov-05		Adapted from OCT20_LC_COMPOSITE
;
DEFAULT, bfactor, 1.064/6.			; corresponds to E^(-3.0) continuum count spectrum
DEFAULT, decx, 12
cs = 1.2				; Charactersize
;
GAMMA_LC_DISPLAY, 'JAN20-LC3.SAV', ENERGY=[2215,2231], RATE=rate1, /NOPLOT		; 06:40-07:10, 10s resolution, 180 PTS
GAMMA_LC_DISPLAY, 'JAN20-LC3.SAV', ENERGY=[2240,2300], RATE=rate2, /NOPLOT
GAMMA_LC_DISPLAY, 'JAN20-LC3.SAV', ENERGY=[250,500],   RATE=rate3, /NOPLOT
GAMMA_LC_DISPLAY, 'JAN20-LC3.SAV', ENERGY=[3500,7000], RATE=rate4, /NOPLOT

!P.MULTI	= [0,1,1]
ncrate 		= rate1.rates-bfactor*rate2.rates
ebrates 	= rate3.rates
nucrates	= rate4.rates
ebkgd		= MEAN([ebrates[0:11],  ebrates[120:*]])		; bkgd interval==06:40-06:42 and 07:00-07:02
nucbkgd		= MEAN([nucrates[0:11],nucrates[120:*]])
ebrates 	= ebrates - ebkgd
nucrates 	= nucrates - nucbkgd
ncrates 	= TOTAL(REFORM(ncrate, 3,60),1)/3			; 1-MINUTE TIME RESOLUTION
trel 		= (FINDGEN(180)+0.5)/6. + 40				; MINUTES AFTER 06:00
trelnc 		= (FINDGEN(60)+0.5)/2 + 40
LOADCT, 0			; BLACK AND WHITE
PLOT, trel, ebrates>0, THICK=3, XTITLE='Minutes after 06:00', TITLE='20-Jan-2005', YTITLE='Relative counts', $
		CHARTHICK=4, CHARSIZE=1.5
LOADCT, 40				; RAINBOW+WHITE
OPLOT, trel, nucrates*13, LINESTYLE=0, THICK=3, COLOR=250		; 3500-7000 in red
OPLOT, trelnc, ncrates*30, PSYM=10, THICK=3, COLOR=64			; 2215-2231 in blue
OPLOT, [46.2, 46.2], [0,40], COLOR=250, LINESTYLE=1, THICK=2
OPLOT, [47.5,47.5], [0,40], COLOR=64, LINESTYLE=1, THICK=2
LOADCT, 0
LOADCT, 0			; BLACK AND WHITE
;OPLOT, trel, ebrates>0, THICK=2
;OPLOT, [44,56], [1.5,1.5], THICK=2
END