From 91f81425751adda276dd7430b95be09e6df790e8 Mon Sep 17 00:00:00 2001 From: Rafael Pezzi Date: Fri, 9 Sep 2011 15:28:18 -0300 Subject: [PATCH] =?UTF-8?q?Incluida=20geometria=20do=20experimento=20na=20?= =?UTF-8?q?compara=C3=A7=C3=A3o=20da=20perda=20de=20Energia=20de=20Bessel?= =?UTF-8?q?=20e=20Gaussiana?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- gaussbesseltest.py | 28 ++++++++++++++++++++-------- 1 file changed, 20 insertions(+), 8 deletions(-) diff --git a/gaussbesseltest.py b/gaussbesseltest.py index bf03c93..5627ffe 100644 --- a/gaussbesseltest.py +++ b/gaussbesseltest.py @@ -1,28 +1,40 @@ # 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) -- 2.22.0