Incluida geometria do experimento na comparação da perda de Energia de Bessel e Gaussiana

gaussbesseltest.py

 # Teste para a perda de energia baseada na funcao modificada de Bessel do primeiro tipo de ordem 1.
-# Apresenta os graficos para a perda de energia para algumas espessuras para feixe com incidencia normal,
+# Apresenta os graficos para a perda de energia para algumas espessuras na geometria definida,
 # comparadas com as correspondentes perdas de energias Gaussianas nas mesmas espessuras.
 # Referencia: R. P. Pezzi, et al. Applied Physics Letters, 92 164102 (2008)
 
 from pylab import *
 from scipy.special import i1
-
+PI=math.pi
 
 dedx=192. # eV/nm
 dw2dx=20740. # eV^2/nm
-e = arange(0,6000)
-PI=math.pi
+Theta_in=0.
+Theta_out=70.
+Theta_s = PI - (Theta_in + Theta_out)*PI/180
+mt = 178.49
+mi = 1
+k = ( (sqrt (mt**2+mi**2*(sin(Theta_s)**2)) + mi*cos(Theta_s) ) / (mi+mt) )**2	
+
+
+dedxef=dedx*(k/cos(Theta_in*PI/180.) + 1/cos(Theta_out*PI/180.))
+dw2dxef=dw2dx*(k**2./cos(Theta_in*PI/180.)+1./cos(Theta_out*PI/180.))
+e = arange(0,23000,20)
+print dedx
+print dedxef/dedx
+print (k/cos(Theta_in*PI/180.) + 1/cos(Theta_out*PI/180.))
 
 for x in [1,5,10,15,20,25]:
     # Bessel:
-    alpha=dedx*(2./dw2dx)
-    m=alpha*dedx
+    alpha=dedxef*(2./dw2dxef)
+    m=alpha*dedxef
     lbd=m*x*alpha
     bessel = lbd*exp(-m*x-alpha*e)*i1(2.*sqrt(lbd*e))/(sqrt(lbd*e))
     plt.plot(e,bessel)
 
     # Gaussiana:
-    Em=x*dedx
-    sigma=x*dw2dx
+    Em=x*dedxef
+    sigma=x*dw2dxef
     gaussian = exp((-(e-Em)**2)/(2.*sigma))/sqrt(2.*PI*sigma) 
     plt.plot(e,gaussian)