65 #ifdef G4MULTITHREADED
147 #ifdef G4MULTITHREADED
148 G4MUTEXLOCK(&G4WentzelOKandVIxSection::WentzelOKandVIxSectionMutex);
161 for(
G4int j=2; j<100; ++j) {
169 #ifdef G4MULTITHREADED
170 G4MUTEXUNLOCK(&G4WentzelOKandVIxSection::WentzelOKandVIxSectionMutex);
255 return cosTetMaxNuc2;
264 if(cosTMax >= 1.0) {
return xSection; }
274 xSection = x2*((1.0 - 1.3333333*x + 3*
x2) - fb*x*(0.6666667 - x));
278 xSection = xlog - x1 - fb*(x + x1 - 2*xlog);
284 G4cout <<
"G4WentzelOKandVIxSection::ComputeTransportCrossSectionPerAtom"
285 <<
" scattering on e- <0"
287 G4cout <<
"cross= " << xSection
288 <<
" e(MeV)= " <<
tkin <<
" p(MeV/c)= " << sqrt(
mom2)
291 G4cout <<
" 1-costm= " << 1.0-costm <<
" screenZ= " <<
screenZ
311 y = x2*((1.0 - 1.3333333*x + 3*
x2) - fb*x*(0.6666667 - x));
315 y = xlog - x1 - fb*(x + x1 - 2*xlog);
321 G4cout <<
"G4WentzelOKandVIxSection::ComputeTransportCrossSectionPerAtom"
322 <<
" scattering on nucleus <0"
360 if(elecRatio > 0.0) {
361 if(rndmEngineMod->
flat() <= elecRatio) {
378 fm =
G4Exp(-2*formf*z1);
392 *fm*fm/(1.0 + z1*
factD);
399 if(cost > 1.0) { cost = 1.0; }
400 else if(cost < -1.0) { cost =-1.0; }
401 G4double sint = sqrt((1.0 - cost)*(1.0 + cost));
404 temp.
set(sint*cos(phi),sint*sin(phi),cost);
419 (1.0 + 2.0*ratio*(tau + 1.0) + ratio*ratio);
void set(double x, double y, double z)
G4double ComputeSecondTransportMoment(G4double CosThetaMax)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
const G4ParticleDefinition * thePositron
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4ThreeVector & SampleSingleScattering(G4double CosThetaMin, G4double CosThetaMax, G4double elecRatio)
std::vector< ExP01TrackerHit * > a
G4double Z13(G4int Z) const
G4double A13(G4double A) const
static constexpr double MeV
const G4Material * currentMaterial
const G4ParticleDefinition * theProton
Float_t x1[n_points_granero]
static constexpr double hbarc
const G4String & GetParticleName() const
G4double GetPDGCharge() const
static G4Proton * Proton()
G4double GetAtomicMassAmu(const G4String &symb) const
G4double GetPDGMass() const
static constexpr double classic_electr_radius
G4double G4Log(G4double x)
static constexpr double proton_mass_c2
virtual ~G4WentzelOKandVIxSection()
static constexpr double fermi
G4double GetPDGSpin() const
static constexpr double amu_c2
#define G4MUTEX_INITIALIZER
static G4Pow * GetInstance()
const G4ParticleDefinition * particle
void ComputeMaxElectronScattering(G4double cut)
void Initialise(const G4ParticleDefinition *, G4double cosThetaLim)
G4double FlatFormfactor(G4double x)
void SetTargetMass(G4double value)
G4NuclearFormfactorType fNucFormfactor
static constexpr double electron_mass_c2
static constexpr double twopi
static G4Positron * Positron()
static constexpr double eV
void SetupParticle(const G4ParticleDefinition *)
G4WentzelOKandVIxSection(G4bool comb=true)
static G4Electron * Electron()
#define G4MUTEXUNLOCK(mutex)
G4double ScreeningFactor() const
G4double SetupTarget(G4int Z, G4double cut)
static constexpr double eplus
void Initialise(const G4ParticleDefinition *, G4double CosThetaLim)
G4double GetInvA23() const
G4NuclearFormfactorType NuclearFormfactorType() const
const G4ParticleDefinition * theElectron
G4double GetA27(G4int Z) const
static G4double ScreenRSquareElec[100]
G4IonisParamMat * GetIonisation() const
static G4double FormFactor[100]
#define G4MUTEXLOCK(mutex)
G4ScreeningMottCrossSection * fMottXSection
G4GLOB_DLL std::ostream G4cout
G4double RatioMottRutherfordCosT(G4double)
G4NistManager * fNistManager
G4double FactorForAngleLimit() const
void SetupKinematic(G4double kinEnergy, G4double Z)
G4double ComputeTransportCrossSectionPerAtom(G4double CosThetaMax)
static constexpr double fine_structure_const
Float_t x2[n_points_geant4]
static G4EmParameters * Instance()
static constexpr double pi
G4double Z23(G4int Z) const
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
static G4NistManager * Instance()
G4double SetupKinematic(G4double kinEnergy, const G4Material *mat)
static G4double ScreenRSquare[100]