thm_sst_dist3d_16x6.pro (Documentation for /home/davin/idl/socware/)




function thm_sst_dist3d_16x6,ion=ion,elec=elec,time,probe=prb,index=index

dat = {thm_sst_dist3d_16x6}

spin_period = 3.

dim = size(/dimension,dat.data)
nenergy = dim[0]
nbins   = dim[1]

;one = replicate(1.,16)
if keyword_set(ion) then begin
    dat.theta = replicate(1,16) # [67.5,-67.5,0,0,0,0.]
    dat.dtheta = replicate(1,16) # [45.,45,90,90,90,90.]

    dat.phi   = replicate(1,16) # [0,0,45,135.,225,315]   ; rotate by 45?
    dat.dphi  = replicate(1,16) # [360.,360.,90.,90.,90.,90]
endif

if keyword_set(elec) then begin
    dat.theta = replicate(1,16) # [67.5,-67.5,0,0,0,0.]
    dat.dtheta = replicate(1,16) # [45.,45,90,90,90,90.]

    dat.phi   = replicate(1,16) # [0,0,45,135.,225,315]
    dat.dphi  = replicate(1,16) # [360.,360.,90.,90.,90.,90]

endif

;Reduced data has larger bins, this results in a larger phi and a larger integration time
dat.integ_t = dat.dphi / 360 * spin_period

if 0 then begin
  ; Warning: not final cals!
  idap_start = [12,19,26,34,44,69,103,150,215,306,506,906,2000,3000,4000,5000]  * 1.5 * 1000
  idap_width = [ 7, 7, 8,10,25,34, 47, 65, 91,200,400,3060,1000,1000,1000,1000] * 1.5 * 1000

  edap_start = [12,19,26,34,44,69,103,150,215,306,506,906,2000,3000,4000,5000]  * 1.5 * 1000
  edap_width = [ 7, 7, 8,10,25,34, 47, 65, 91,200,400,2000,1000,1000,1000,1000] * 1.5 * 1000

  energy = (2*idap_start + idap_width)/2   + 5000.       ; midpoint energy
  dat.energy = energy # replicate(1,nbins)       ; total energy width
  denergy = (idap_width)
  dat.denergy = denergy # replicate(1,nbins)
endif else begin ;not really the final cals either...but closer...
   thm_sst_energy_cal,energy=energy,denerg=denergy,inst = keyword_set(ion),probe=prb
   ;detector = [[replicate(0,16)],[replicate(1,16)],[replicate(2,16)],[replicate(3,16)]]
   ;anode reduction means that anode by anode corrections don't exactly apply.  This code uses a kluge, the average energy & denergy for angle bins that include multiple anodes 
   dat.energy  =  [energy[*,0],energy[*,1],(energy[*,2]+energy[*,3])/2.,(energy[*,2]+energy[*,3])/2.,(energy[*,2]+energy[*,3])/2.,(energy[*,2]+energy[*,3])/2.]
   dat.denergy = [denergy[*,0],denergy[*,1],(denergy[*,2]+denergy[*,3])/2.,(denergy[*,2]+denergy[*,3])/2.,(denergy[*,2]+denergy[*,3])/2.,(denergy[*,2]+denergy[*,3])/2.]
endelse


;weights = calc_omega_flt2(dat.theta,dat.phi,dat.dtheta,dat.dphi, 1.)
;dat.domega = weights[0,*,*]

dat.nenergy = nenergy
dat.nbins   = nbins
dat.bins = 1

dat.gf = replicate(1,nenergy) # [1.,1.,2.,2.,2.,2.] ;correction to notional GF for reduced data.  The correction in the spin-plane(phi) is accounted for by increasing the integration
                             ;time of the data, so this correction only includes the theta component. 
                             ;(Each of the spin plane bins contain the contribution from 2 different anodes, thus doubling the effective GF along the theta axis.)

dat.units_procedure = 'thm_sst_convert_units'

dat.geom_factor = .1


dat.eclipse_dphi = !values.d_nan

return,dat
end