;pro map_sep27

goto,maps_8

obs_time_interval=['2002/09/27 19:34:00','2002/09/27 19:42:00']
map_time_interval=['2002/09/27 19:35:00','2002/09/27 19:36:00']
map_time_interval=['2002/09/27 19:36:00','2002/09/27 19:37:00']
                   
xy_offset=[-585,123] ;[0,0]

; Keep the same for each flare:
    energy_band=[12.,25.]
    energy_band=[6.,25.]

    det_index_mask=[0,0,1,1,1,1,1,1,0] ;[0,0,0,0,0,0,1,1,1] ; latter for full disk
    pixel_size=[2,2] ;[32,32] use the latter for full disk
    image_dim=[64,64]

algorithm = ['back_proj','mem_vis','clean','mem_sato','forward_fit','pixons']


image_algorithm=algorithm[0]
if (!d.name NE 'X') then image_algorithm=algorithm[0]
xyoffset=xy_offset

obj_im = hsi_image() 
  obj_im -> set,use_auto_time_bin=0
  obj_im -> set,obs_time_interval=obs_time_interval
  obj_im -> set,energy_band=energy_band ;
  obj_im -> set,image_algorithm=image_algorithm
  obj_im -> set,pixel_size=pixel_size 
  obj_im -> set,det_index_mask=det_index_mask ;[0,0,0,0,0,0,1,1,1] ;
  obj_im -> set,time_bin_def=[1,1,1,1,1,1,1,1,1] 
  obj_im -> set,image_dim=image_dim
if image_algorithm eq 'clean' then $
  obj_im -> set,clean_niter=333

istart=anytim(map_time_interval[0])-anytim(obs_time_interval[0])
iend=anytim(map_time_interval[1])-anytim(obs_time_interval[0])

obj_im -> set,time_range=[istart,iend] ; if this had been set to anytim units, then
                                ; obs_time_interval will be set equal
                                ; to time_range (see
                                ; http://hessi.ssl.berkeley.edu/
                                ;~cmj/hessi/doc.html#Introduction

p=obj_im->get()  ; get the parameters
print,'xyoffset=',p.xyoffset
print,'obs_time_interval=',anytim(p.obs_time_interval,/yo)
print,'time_range=',p.time_range

print,'det_index_mask=',p.det_index_mask
print,'energy_band=',p.energy_band
print,'pixel_size=',p.pixel_size
print,'image_dim=',p.image_dim
print,'time_bin_def=',p.time_bin_def
print,'time_range=',p.time_range
print,'use_auto_time_bin=',p.use_auto_time_bin

print,'Getting the map data...'
im=obj_im -> getdata(obs_time_interval=obs_time_interval, $
                     XYOFFSET=xyoffset,                   $
                     energy_band=energy_band,             $
                     use_auto_time_bin=0,                 $
                     det_index_mask=[0,0,1,1,1,1,1,1,0] ) 
;                     time_bin_def=[1,1,1,1,1,1,1,1,1] ) 

if (!d.name EQ 'X') then obj_im->plotman

; Getting calibrated eventlist...'
;cbe=obj_im->getdata(class_name='hsi_calib_eventlist',$
;          energy_band=energy_band,det_index_mask=[1,1,1,1,1,1,1,1,1],$
;          OBS_TIME_INTERVAL=obs_time_interval, $
;          XYOFFSET=xyoffset,                   $
;          use_auto_time_bin=0,                 $
;          time_bin_def=[1,1,1,1,1,1,1,1,1])

maps_8:

;istart= 0.0000000   iend=       60.000000 ; 1934-1935
istart=420   &  iend=  480 ; 1941-1942

energy_band=[12,25]
for j=0,7 do begin

  istart=istart-60 & iend=iend-60
  obj_im -> set,time_range=[istart,iend] 
  p=obj_im->get()  ; get the parameters
  print,'time_range=',p.time_range
  print,'det_index_mask=',p.det_index_mask
  print,'energy_band=',energy_band[0],energy_band[1]

  im=obj_im -> getdata(obs_time_interval=obs_time_interval, $
                     XYOFFSET=xyoffset,                   $
                     energy_band=energy_band,             $
                     det_index_mask=[0,0,1,1,1,1,1,1,0] ) 

  obj_im->plotman
endfor




end