| Date 08/09/2010 |
|
||||||
|
| Script S6_2_10.m | |||
|
|||
|
%==============================================
%compton 3 %============================================== % %speed light in vacuum (in m/sec), Planck constant (in Joulexsec) %rest mass of the electron (in kg) c=2.998e008; h=6.626e-034; m0=9.109e-031; %constant hc hc=h*c; %Compton constant const=h/(c*m0) %wavelength of the incident beam in meters lambda=12.0*1.0e-012; %scattering angle of the electrons in radians fi=59.6*pi/180; cosfi=cos(fi); cosfi2=cosfi^2; %energy of the incident photons in Joule Wf=hc/lambda; Wf2=Wf^2; %and in keV WfkeV=(Wf/1.6)*1.0e+016 %preliminary calculi m0c2=m0*c^2; num=2*Wf2*m0c2*cosfi2; den1=(m0c2+Wf)^2; den2=Wf2*cosfi2; den=den1+den2; %kinetic energy of electrons in Joule T=num/den; %and in keV Tkev=(T/1.6)*1.0e+016 %energy of the scattered photons %W'f in keV WpfkeV=WfkeV-Tkev %and in Joule Wpf=(WfkeV-Tkev)*(1.6*1.0e-016); %wavelength lambdap of the scattered photons in meters lambdap=hc/Wpf; %and in picometers lambdap_pico=lambdap*(1.0e+012) % const1=(lambdap-lambda)/const coseno=1-const1; %angle of the scattered photons in degrees tetag=(180/pi)*acos(coseno) %============================================== % |
|||
| Top | |||
