101 if ( Z < 0.9999 )
return xSection;
102 if ( GammaEnergy < 0.01*
keV )
return xSection;
105 static const G4double a = 20.0 ,
b = 230.0 ,
c = 440.0;
108 d1= 2.7965e-1*
barn, d2=-1.8300e-1*
barn, d3= 6.7527 *
barn, d4=-1.9798e+1*
barn,
109 e1= 1.9756e-5*
barn, e2=-1.0205e-2*
barn, e3=-7.3913e-2*
barn, e4= 2.7079e-2*
barn,
112 G4double p1Z = Z*(d1 + e1*Z + f1*Z*
Z), p2Z = Z*(d2 + e2*Z + f2*Z*Z),
113 p3Z = Z*(d3 + e3*Z + f3*Z*
Z), p4Z = Z*(d4 + e4*Z + f4*Z*Z);
116 if (Z < 1.5) T0 = 40.0*
keV;
119 xSection = p1Z*std::log(1.+2.*X)/X
120 + (p2Z + p3Z*X + p4Z*X*
X)/(1. + a*X +
b*X*X +
c*X*X*X);
123 if (GammaEnergy < T0) {
127 + (p2Z + p3Z*X + p4Z*X*
X)/(1. + a*X +
b*X*X +
c*X*X*X);
128 G4double c1 = -T0*(sigma-xSection)/(xSection*dT0);
130 if (Z > 1.5) c2 = 0.375-0.0556*log(Z);
132 xSection *= exp(-y*(c1+c2*y));
160 eMomDir *= std::sqrt(eMomentumC2);
182 G4double epsilon, epsilonsq, onecost, sint2, greject ;
187 G4double alpha2 = 0.5*(1.- epsilon0sq);
194 epsilonsq = epsilon*epsilon;
200 epsilon = sqrt(epsilonsq);
203 onecost = (1.- epsilon)/(epsilon*E0_m);
204 sint2 = onecost*(2.-onecost);
205 greject = 1. - epsilon*sint2/(1.+ epsilonsq);
216 G4double dirx = sinTeta*cos(Phi), diry = sinTeta*sin(Phi), dirz = cosTeta;
223 gamDirection1.
rotateUz(gamDirection0);
224 G4double gamEnergy1 = epsilon*gamEnergy0;
225 gamDirection1 *= gamEnergy1;
233 G4double eKinEnergy = gamEnergy0 - gamEnergy1;
234 G4ThreeVector eDirection = gamEnergy0*gamDirection0 - gamEnergy1*gamDirection1;
235 eDirection = eDirection.
unit();
236 G4double eFinalMom = std::sqrt(eKinEnergy*(eKinEnergy+2*electron_mass_c2));
237 eDirection *= eFinalMom;
240 gamma4vfinal.
boost(bst);
243 gamDirection1 = gamma4vfinal.
vect();
244 gamEnergy1 = gamDirection1.
mag();
245 gamDirection1 /= gamEnergy1;
254 gamDirection1 /= gamEnergy1;
269 eDirection = e4vfinal.
vect();
271 eDirection /= eFinMomMag;
273 fvect->push_back(dp);
ThreeVector shoot(const G4int Ap, const G4int Af)
Hep3Vector boostVector() const
G4ParticleChangeForGamma * fParticleChange
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A, G4double cut, G4double emax)
G4ParticleDefinition * theGamma
G4ThreeVector G4RandomDirection()
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
HepLorentzVector & boost(double, double, double)
Hep3Vector & rotateUz(const Hep3Vector &)
virtual ~G4HeatedKleinNishinaCompton()
G4HeatedKleinNishinaCompton(const G4ParticleDefinition *p=0, const G4String &nam="Heated-Klein-Nishina")
G4LorentzVector Get4Momentum() const
G4ParticleDefinition * theElectron
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double lowestGammaEnergy
G4double GetLocalEnergyDeposit() const
static G4Electron * Electron()
void SetProposedKineticEnergy(G4double proposedKinEnergy)
void ProposeTrackStatus(G4TrackStatus status)
G4ParticleChangeForGamma * GetParticleChangeForGamma()
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
CLHEP::HepLorentzVector G4LorentzVector