119 static const G4double a = 20.0 , b = 230.0 , c = 440.0;
122 p3Z = Z*(
d3 +
e3*Z +
f3*Z*Z), p4Z = Z*(
d4 +
e4*Z +
f4*Z*Z);
125 if (Z < 1.5) { T0 = 40.0*
keV; }
127 G4double X =
max(GammaEnergy, T0) / electron_mass_c2;
128 xSection = p1Z*
G4Log(1.+2.*X)/X
129 + (p2Z + p3Z*X + p4Z*X*X)/(1. + a*X + b*X*X + c*X*X*X);
132 if (GammaEnergy < T0) {
134 X = (T0+
dT0) / electron_mass_c2 ;
136 + (p2Z + p3Z*X + p4Z*X*X)/(1. + a*X + b*X*X + c*X*X*X);
137 G4double c1 = -T0*(sigma-xSection)/(xSection*dT0);
139 if (Z > 1.5) { c2 = 0.375-0.0556*
G4Log(Z); }
141 xSection *=
G4Exp(-y*(c1+c2*y));
150 std::vector<G4DynamicParticle*>* fvect,
166 G4double E0_m = gamEnergy0 / electron_mass_c2 ;
179 G4double alpha2 = alpha1 + 0.5*(1.- epsilon0sq);
181 CLHEP::HepRandomEngine* rndmEngineMod = G4Random::getTheEngine();
191 rndmEngineMod->flatArray(3, rndm);
193 if ( alpha1 > alpha2*rndm[0] ) {
194 epsilon =
G4Exp(-alpha1*rndm[1]);
198 epsilonsq = epsilon0sq + (1.- epsilon0sq)*rndm[1];
199 epsilon = sqrt(epsilonsq);
202 onecost = (1.-
epsilon)/(epsilon*E0_m);
203 sint2 = onecost*(2.-onecost);
204 greject = 1. - epsilon*sint2/(1.+ epsilonsq);
207 }
while (greject < rndm[2]);
213 if(sint2 < 0.0) { sint2 = 0.0; }
222 G4ThreeVector gamDirection1(sinTeta*cos(Phi), sinTeta*sin(Phi), cosTeta);
223 gamDirection1.rotateUz(gamDirection0);
224 G4double gamEnergy1 = epsilon*gamEnergy0;
239 G4double eKinEnergy = gamEnergy0 - gamEnergy1;
242 G4ThreeVector eDirection = gamEnergy0*gamDirection0 - gamEnergy1*gamDirection1;
243 eDirection = eDirection.unit();
247 fvect->push_back(dp);
static c2_factory< G4double > c2
G4double LowEnergyLimit() const
G4ParticleChangeForGamma * fParticleChange
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
void InitialiseElementSelectors(const G4ParticleDefinition *, const G4DataVector &)
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
static const G4int nlooplim
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A, G4double cut, G4double emax)
const G4ThreeVector & GetMomentumDirection() const
static const double twopi
G4double lowestSecondaryEnergy
std::vector< G4EmElementSelector * > * GetElementSelectors()
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
void SetElementSelectors(std::vector< G4EmElementSelector * > *)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4KleinNishinaCompton(const G4ParticleDefinition *p=nullptr, const G4String &nam="Klein-Nishina")
static G4Electron * Electron()
void SetProposedKineticEnergy(G4double proposedKinEnergy)
G4ParticleDefinition * theElectron
void ProposeTrackStatus(G4TrackStatus status)
virtual ~G4KleinNishinaCompton()
virtual void InitialiseLocal(const G4ParticleDefinition *, G4VEmModel *masterModel)
double epsilon(double density, double temperature)
G4ParticleDefinition * theGamma
static const G4double dT0
G4ParticleChangeForGamma * GetParticleChangeForGamma()
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)