;+
;FUNCTION: e_sec_eff(Te)
;INPUT:
; Te: flt,fltarr electron temperature
;KEYWORDS
; pot flt,fltarr if fltarr, then n_elements(pot)=n_elements(te), default pot=0
; tmax flt algorithm parameter, default tmax = 2.283
; aa flt algorithm parameter, default aa = 1.35
; emax flt algorithm parameter, default emax = 325
; kk flt algorithm parameter, default kk = 2.2
; sec_max flt algorithm parameter, default sec_max=2.
; te_sec flt algorithm parameter, default te_sec=2.
;
;PURPOSE:
; Returns the fractional secondary electron flux with energy > pot
;NOTES:
;
;CREATED BY:
; J.McFadden 09-05-31
;LAST MODIFICATION:
;-
function e_sec_eff,Te,pot=pot,tmax=tmax,aa=aa,emax=emax,kk=kk,sec_max=sec_max
; Make electron secondary energy efficiency array
; formula is an offshoot of the following MCP efficiency formula
; REF: Relative electron detection efficiency of microchannel plates from 0-3 keV
; R.R. Goruganthu and W. G. Wilson, Rev. Sci. Instrum. Vol. 55, No. 12 Dec 1984
if not keyword_set(pot) then pot = 0
if not keyword_set(tmax) then tmax = 2.283
if not keyword_set(aa) then aa = 1.35
if not keyword_set(emax) then emax = 325. ; 325 eV
if not keyword_set(kk) then kk = 2.2
if not keyword_set(sec_max) then sec_max=2.
if not keyword_set(te_sec) then te_sec=2.
delta = ((Te+(pot>0.))/emax)^(1-aa)*(1-exp(-tmax*((Te+(pot>0.))/emax)^aa))/(1-exp(-tmax))
en_eff = exp(-20./Te)/exp(-20./emax)*exp(-(pot>0.)/te_sec)*exp((pot<0.)/Te)*sec_max*(1-exp(-kk*delta))/(1-exp(-kk))
return, en_eff
end