50 #ifndef G4PairProductionRelModel_h
51 #define G4PairProductionRelModel_h 1
67 const G4String& nam =
"BetheHeitlerLPM");
221 screenVal = 21.12 - 4.184*std::log(delta+0.952);
223 screenVal = 20.868 - delta*(3.242 - 0.625*delta);
235 screenVal = 21.12 - 4.184*std::log(delta+0.952);
237 screenVal = 20.209 - delta*(1.930 + 0.086*delta);
249 if (ScreenVariable > 1.)
250 screenVal = 42.24 - 8.368*std::log(ScreenVariable+0.952);
252 screenVal = 42.392 - ScreenVariable*(7.796 - 1.961*ScreenVariable);
266 if (ScreenVariable > 1.)
267 screenVal = 42.24 - 8.368*std::log(ScreenVariable+0.952);
269 screenVal = 41.405 - ScreenVariable*(5.828 - 0.8945*ScreenVariable);
281 return std::exp( (42.24-FZ)/8.368 ) + 0.952;
286 return 4.*136./
z13*(CLHEP::electron_mass_c2/k);
G4double Phi1(G4double delta) const
G4bool use_completescreening
G4double ComputeXSectionPerAtom(G4double totalEnergy, G4double Z)
void SetLPMconstant(G4double val)
G4double Phi2(G4double delta) const
G4PairProductionRelModel(const G4ParticleDefinition *p=0, const G4String &nam="BetheHeitlerLPM")
G4double GetZ13(G4double Z)
G4double GetfCoulomb() const
static const G4double preS1
static const G4double Finel_light[5]
G4ParticleDefinition * thePositron
static const G4double facFel
static const G4double xgi[8]
G4PairProductionRelModel & operator=(const G4PairProductionRelModel &right)
G4double ScreenFunction1(G4double ScreenVariable)
virtual void SetupForMaterial(const G4ParticleDefinition *, const G4Material *, G4double)
virtual void InitialiseLocal(const G4ParticleDefinition *, G4VEmModel *masterModel)
static const G4double wgi[8]
static const G4double logTwo
virtual ~G4PairProductionRelModel()
G4double DeltaMin(G4double) const
void CalcLPMFunctions(G4double k, G4double eplus)
static const G4double facFinel
G4double ScreenFunction2(G4double ScreenVariable)
G4double LPMconstant() const
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
G4ParticleChangeForGamma * fParticleChange
G4ParticleDefinition * theGamma
void SetCurrentElement(G4double)
G4double DeltaMax() const
G4double GetLOGZ(G4int Z)
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
G4double ComputeRelDXSectionPerAtom(G4double eplusEnergy, G4double totalEnergy, G4double Z)
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0., G4double cut=0., G4double emax=DBL_MAX)
G4ParticleDefinition * theElectron
static const G4double Fel_light[5]
G4double ComputeDXSectionPerAtom(G4double eplusEnergy, G4double totalEnergy, G4double Z)
const G4Element * GetCurrentElement() const