Six Components

SCRIPT SIMULATING 6-COMPONENT SOURCE

; sim_multiple_src
gg1 = {gaussian_source_str}
  gg1.amplitude = 1.
  gg1.xysigma = [1., 1.]
  gg1.xypos = [-25,-5.]
  gg1.tilt_angle_deg = 0.
 gg2 = {gaussian_source_str}
  gg2.amplitude = 1.
  gg2.xysigma = [1., 1.]
  gg2.xypos = [25.,5.]
  gg2.tilt_angle_deg = 0.
 gg3 = {gaussian_source_str}
  gg3.amplitude = 0.5
  gg3.xysigma = [1., 1.]
  gg3.xypos = [-20.,-5.]
  gg3.tilt_angle_deg = 0.
 gg4 = {gaussian_source_str}
  gg4.amplitude = 0.5
  gg4.xysigma = [1., 1.]
  gg4.xypos = [-15.,5.]
  gg4.tilt_angle_deg = 0.

 gg5 = {gaussian_source_str}
  gg5.amplitude = 0.5
  gg5.xysigma = [1., 1.]
  gg5.xypos = [15.,5.]
  gg5.tilt_angle_deg = 0.

 gg6 = {gaussian_source_str}
  gg6.amplitude = 0.5
  gg6.xysigma = [1., 1.]
  gg6.xypos = [20.,10.]
  gg6.tilt_angle_deg = 0.

  
  
  
model = [gg1,gg2,gg3,gg4,gg5,gg6]
obj_im = hsi_image() 
obj_im->set,pixel_size=4
obj_im->set,det_index=[0,1,1,1,1,1,1,1,0]

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

cbe=obj_im->getdata(class_name='hsi_calib_eventlist')

!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'

obj_im->plotman

end