Two-Second Maps

IDL SCRIPT USED FOR SIMULATION

  ; 2 Gaussian components for sim_timerange=[0,1]:
  gg1 = {gaussian_source_str}
    gg1.amplitude = 0.5               ; 1.0 for sim_timerange=[1,2]
    gg1.xysigma = [1., 1.]
    gg1.xypos = [-25.,10.]
    gg1.tilt_angle_deg = 0.
  gg2 = {gaussian_source_str}
    gg2.amplitude = 1.0               ; 0.5 for sim_timerange=[1,2]
    gg2.xysigma = [1., 1.]
    gg2.xypos = [25.,-10.]
    gg2.tilt_angle_deg = 0.

  model = [gg1,gg2]
  obj_im = hsi_image() 
  obj_im->set,pixel_size=2
  obj_im->set,det_index=[0,1,1,1,1,1,1,1,0]
  obj_im->set,time_range=[0,1]    
; obj_im->set,sim_time_range=[1,2]  

  im = obj_im -> getdata(/plot,sim_model=model,sim_photons=6000,$
  sim_energy_band=[6,100],energy_band=[6,100],$
  /force) 

; Get the calibrated eventlists if you want them:
; cbe=obj_im->getdata(class_name='hsi_calib_eventlist')

: plot if you desire:
; !p.multi=[0,1,4]
; !x.style=1 & !P.CHARSIZE=2
; tj=['1)','2)','3)','4)']
; tj=['5)','6)','7)','8)']
; ti="OVER RESOLVED COUNT PROFILE FOR 15'' SOURCE (DET="
; for j=0,3 do plot,(*cbe(j+4,0)).count,psym=10,$
;  TIT=TI+TJ[j],xtIT='TIME BIN',YTIT='COUNTS/S'

; Use plotman to plot if graphic control is needed:
  obj_im->plotman

 end