66 d1= 2.7965e-1*
barn, d2=-1.8300e-1*
barn, d3= 6.7527 *
barn, d4=-1.9798e+1*
barn,
67 e1= 1.9756e-5*
barn, e2=-1.0205e-2*
barn, e3=-7.3913e-2*
barn, e4= 2.7079e-2*
barn,
68 f1=-3.9178e-7*
barn, f2= 6.8241e-5*
barn, f3= 6.0480e-5*
barn, f4= 3.0274e-4*
barn;
111 if ( Z < 0.9999 )
return xSection;
112 if ( GammaEnergy < 0.1*
keV )
return xSection;
115 static const G4double a = 20.0 ,
b = 230.0 ,
c = 440.0;
117 G4double p1Z = Z*(d1 + e1*Z + f1*Z*Z), p2Z = Z*(d2 + e2*Z + f2*Z*Z),
118 p3Z = Z*(d3 + e3*Z + f3*Z*Z), p4Z = Z*(d4 + e4*Z + f4*Z*Z);
121 if (Z < 1.5) T0 = 40.0*
keV;
124 xSection = p1Z*
G4Log(1.+2.*X)/X
125 + (p2Z + p3Z*X + p4Z*X*X)/(1. + a*X +
b*X*X +
c*X*X*X);
128 if (GammaEnergy < T0) {
132 + (p2Z + p3Z*X + p4Z*X*X)/(1. + a*X +
b*X*X +
c*X*X*X);
133 G4double c1 = -T0*(sigma-xSection)/(xSection*dT0);
135 if (Z > 1.5) c2 = 0.375-0.0556*
G4Log(Z);
137 xSection *=
G4Exp(-y*(c1+c2*y));
174 G4double epsilon, epsilonsq, onecost, sint2, greject ;
179 G4double alpha2 = 0.5*(1.- epsilon0sq);
184 epsilonsq = epsilon*epsilon;
188 epsilon = sqrt(epsilonsq);
191 onecost = (1.- epsilon)/(epsilon*E0_m);
192 sint2 = onecost*(2.-onecost);
193 greject = 1. - epsilon*sint2/(1.+ epsilonsq);
201 if(sint2 < 0.0) { sint2 = 0.0; }
210 G4ThreeVector gamDirection1(sinTeta*cos(Phi), sinTeta*sin(Phi), cosTeta);
211 gamDirection1.
rotateUz(gamDirection0);
212 G4double gamEnergy1 = epsilon*gamEnergy0;
226 G4double eKinEnergy = gamEnergy0 - gamEnergy1;
229 G4ThreeVector eDirection = gamEnergy0*gamDirection0 - gamEnergy1*gamDirection1;
230 eDirection = eDirection.
unit();
234 fvect->push_back(dp);
G4ParticleChangeForGamma * fParticleChange
G4KleinNishinaCompton(const G4ParticleDefinition *p=0, const G4String &nam="Klein-Nishina")
G4double GetKineticEnergy() const
void InitialiseElementSelectors(const G4ParticleDefinition *, const G4DataVector &)
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
G4double lowestGammaEnergy
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A, G4double cut, G4double emax)
const G4ThreeVector & GetMomentumDirection() const
Hep3Vector & rotateUz(const Hep3Vector &)
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
static G4Electron * Electron()
void SetProposedKineticEnergy(G4double proposedKinEnergy)
G4ParticleDefinition * theElectron
void ProposeTrackStatus(G4TrackStatus status)
virtual ~G4KleinNishinaCompton()
virtual void InitialiseLocal(const G4ParticleDefinition *, G4VEmModel *masterModel)
G4ParticleDefinition * theGamma
G4ParticleChangeForGamma * GetParticleChangeForGamma()
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)