; JAN20GAMMA_IMG.PRO

; Script for imaging January 20, 2005 event at gamma ray energies
; Adapted from gamma_img.pro
;
;xyoffset 	= [0,0]		; default
;pixels		= [16,16]	; default
;
; JAN 20, 2005 EVENT
oim = hsi_image()
;
	;oim->set, time_range			= '2005-jan-20 '+['06:46:00', '06:47:00']		; 1 minute near peak
	;oim->set, time_range			= '2005-jan-20 '+['06:44:00', '07:02:00']		; Full 2.2 MeV time range
	;oim->set, time_range			= '2005-jan-20 '+['06:44:00', '07:26:00']		; Full flare
	;oim->set, time_range			= '2005-jan-20 '+['06:44:00', '06:54:00']		; peak
oim->set, time_range			= '2005-jan-20 '+['06:44:00', '06:56:00']		; full peak
	;oim->set, time_range			= '2005-jan-20 '+['06:56:00', '07:26:00']		; Post peak -- flat 511 period
;
	;oim-> set, energy_band			= [1000,2218]
	;oim->set, energy_band			= [1000,7000]
	;oim->set, energy_band			= [3500,7000]
	;oim->set, energy_band			= [1500,2000]
	;oim->set, energy_band			= [7000, 15000]
	;oim-> set, energy_band			= [2218,2228]		; N capture line
	;oim-> set, energy_band			= [2215,2231]		; optimized for N capture line
	;oim-> set, energy_band			= [2230,2280]		; Background for N capture line
	;oim->set, energy_band			= [2300,3500]		; Lower nuclear fraction
	;oim-> set, energy_band			= [3500, 7000]		; Broad nuclear lines
	;oim-> set, energy_band			= [1500, 2100]		; 2.2 compton edge
	;oim-> set, energy_band			= [1000, 1500]		; mid 2.2 compton edge + Brems
	;oim-> set, energy_band			= [550, 1000]		; high E bremsstrahlung
oim->set, energy_band			= [250,500]			; low E bremsstrahlung
	;oim-> set, energy_band			= [3250,6500]
	;oim-> set, energy_band			= [700,1400]
	;oim-> set, energy_band			= [1400, 2000]
	;oim-> set, energy_band			= [700, 2000]
	;oim-> set, energy_band			= [2219,2227]
	;oim-> set, energy_band			= [2228,2400]
	;oim-> set, energy_band			= [2400,4000]
	;oim-> set, energy_band			= [485,505]
	;oim-> set, energy_band			= [4000,7000]
	;oim-> set, energy_band			= [1300,1820]
	;oim-> set, energy_band			= [2000,4000,7000], eb_index=1
	;oim-> set, energy_band			= [300, 507, 509, 513, 515], eb_index=1
	;oim-> set, energy_band			= [1300, 1820, 2219, 2227], eb_index=2
	;oim-> set, energy_band			= [1356, 1373, 1618, 1639, 1761, 1783], eb_index=0
	;oim-> set, energy_band			= [508,514.]			; narrow 511 line
	;oim-> set, energy_band			= [500,520]			; full 511 line
	;oim-> set, energy_band			= [525,547]			; continuum for 511 line
;
oim-> set, det_index_mask		= [0,0,0, 1,1,1, 1,1,1]			; grids 4,5,6,7,8,9
	;oim-> set, det_index_mask			= [0,0,0, 0,0,1, 0,0,1]		; grids 6,9 only
	;oim-> set, det_index_mask			= [0,0,0, 0,1,1, 1,1,1]		; grids 5,6,7,8,9
	;oim-> set, det_index_mask			= [0,0,0, 0,0,0, 0,0,1]		; grid 9 only
	;oim-> set, det_index_mask			= [0,0,0, 0,0,1, 0,0,0]		; grids 6 only
;
oim-> set, use_auto_time_bin 	= 0
	;oim-> set, time_bin_def			= [128,128,128, 128,128,8, 128,128,32]				; grids 6,9 only
	;oim-> set, time_bin_def				= [1024,1024,1024, 1024,1024,1024, 1024,1024,32]	; grid 9 only
	;oim-> set, time_bin_def			= [1024,1024,1024, 1024,1024,10, 1024,1024,1024]	; grid 6 only
	;oim-> set, time_bin_def			= [128,128,  2,   4, 8,8,    16,32,64]
oim-> set, time_bin_def			= [128,128, 128,   4, 8,8,    16,32,64]			; grids 4,5,6,7,8,9
	;oim-> set, time_bin_def			= [128,128,128, 128,8,8,    16,32,64]				; grids 5,6,7,8,9
;
;
	;Pixels 							= [8,8]
	;pixels 							= [3,3]
pixels 							= [2,2]
	;pixels 							= [16,16]
;
xyoffset						= [806, 260]
	;image_dim						= [128,128]
image_dim						= [64,64]
;
oim-> set, xyoffset				= xyoffset
oim-> set, image_dim			= image_dim
oim-> set, pixel_size 			= pixels
oim-> set, UNIFORM_WEIGHTING 	= 1
oim-> set, NATURAL_WEIGHTING 	= 0
oim-> set, rear_segment			= 1
oim-> set, front_segment		= 0
oim-> set, use_rate				= 1
oim-> set, sum_coincidence		= 1				; ignore front/rear coincidences
oim-> set, image_algorithm		= 'clean'
oim-> set, clean_niter		= 200
	;oim-> set, image_algorithm		= 'back'
oim-> set, use_rate				= 1				; rate-based
oim-> set, use_local_average	= 1				; enable local averaging
oim-> set, auto_frequency		= 0				; do NOT use formula
oim-> set, local_average_frequency	= 1.0
oim->set, ev_filename = 'c:\hessi_data\jan20_eventlist_rearA.fits'
;
oim->plotman
;
PRINT, 'GAMMA_IMG is done.'
hsi_map_evaluator, oim, /TEXT
END