73 alpha2(fine_structure_const*fine_structure_const)
81 G4double p0 = electron_mass_c2*classic_electr_radius;
94 for(
G4int j=1; j<100; ++j) {
200 return cosTetMaxNuc2;
209 if(cosTMax >= 1.0) {
return xSection; }
225 y = x2*(1.0 - 1.3333333*x + 3*x2);
226 if(0.0 < factB) { y -= fb*x2*x*(0.6666667 -
x); }
229 xlog =
G4Log(1.0 + x);
231 if(0.0 < factB) { y -= fb*(x + x1 - 2*xlog); }
237 G4cout <<
"G4WentzelOKandVIxSection::ComputeTransportCrossSectionPerAtom"
238 <<
" scattering on e- <0"
241 <<
" e(MeV)= " <<
tkin <<
" p(MeV/c)= " << sqrt(
mom2)
244 G4cout <<
" 1-costm= " << 1.0-costm <<
" screenZ= " <<
screenZ
264 y = x2*(1.0 - 1.3333333*x + 3*x2);
265 if(0.0 < factB) { y -= fb*x2*x*(0.6666667 -
x); }
268 xlog =
G4Log(1.0 + x);
270 if(0.0 < factB) { y -= fb*(x + x1 - 2*xlog); }
276 G4cout <<
"G4WentzelOKandVIxSection::ComputeTransportCrossSectionPerAtom"
277 <<
" scattering on nucleus <0"
283 <<
" x= " <<
" x1= " << x1 <<
G4endl;
309 temp.set(0.0,0.0,1.0);
310 CLHEP::HepRandomEngine* rndmEngineMod = G4Random::getTheEngine();
315 if(elecRatio > 0.0) {
316 if(rndmEngineMod->flat() <= elecRatio) {
329 /((1.0 + z1*
factD)*fm*fm);
331 if(rndmEngineMod->flat() <= grej) {
334 if(cost > 1.0) { cost = 1.0; }
335 else if(cost < -1.0) { cost =-1.0; }
336 G4double sint = sqrt((1.0 - cost)*(1.0 + cost));
339 temp.set(sint*cos(phi),sint*sin(phi),cost);
353 G4double tmax = 2.0*electron_mass_c2*tau*(tau + 2.)/
354 (1.0 + 2.0*ratio*(tau + 1.0) + ratio*ratio);
363 G4double mom21 = t*(t + 2.0*electron_mass_c2);
const G4ParticleDefinition * theElectron
const G4ParticleDefinition * theProton
static G4Pow * GetInstance()
void SetupParticle(const G4ParticleDefinition *)
static G4double ScreenRSquareElec[100]
void SetTargetMass(G4double value)
G4WentzelOKandVIxSection(G4bool combined=true)
CLHEP::Hep3Vector G4ThreeVector
static G4double FormFactor[100]
G4double ComputeSecondTransportMoment(G4double CosThetaMax)
const G4ParticleDefinition * particle
void ComputeMaxElectronScattering(G4double cut)
static G4NistManager * Instance()
G4double SetupTarget(G4int Z, G4double cut=DBL_MAX)
const G4String & GetParticleName() const
virtual ~G4WentzelOKandVIxSection()
G4double ComputeTransportCrossSectionPerAtom(G4double CosThetaMax)
G4GLOB_DLL std::ostream G4cout
G4double Z13(G4int Z) const
static const double twopi
static G4Proton * Proton()
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
static G4Positron * Positron()
G4double GetPDGMass() const
const G4ParticleDefinition * thePositron
T max(const T t1, const T t2)
brief Return the largest of the two arguments
const G4Material * currentMaterial
const G4double x[NPOINTSGL]
static G4EmParameters * Instance()
G4double GetAtomicMassAmu(const G4String &symb) const
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4NistManager * fNistManager
G4double Z23(G4int Z) const
G4double GetPDGSpin() const
static G4Electron * Electron()
static G4double ScreenRSquare[100]
static const double eplus
G4double GetPDGCharge() const
static const double fermi
void Initialise(const G4ParticleDefinition *, G4double CosThetaLim)
G4double FactorForAngleLimit() const
G4ThreeVector & SampleSingleScattering(G4double CosThetaMin, G4double CosThetaMax, G4double elecRatio=0.0)