81 if(particle ==
fElectron) { maximumEnergyTransfer *= 0.5; }
83 G4double ratio = electron_mass_c2/mass;
85 maximumEnergyTransfer = 2.0*electron_mass_c2*tau*(tau + 2.) /
86 (1. + 2.0*(tau + 1.)*ratio + ratio*ratio);
90 if (secKinetic>100*
eV && secKinetic <= maximumEnergyTransfer) {
91 cosTheta = std::sqrt(secKinetic / maximumEnergyTransfer);
96 G4double sint = sqrt((1.0 - cosTheta)*(1.0 + cosTheta));
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
G4DNARuddAngle(const G4String &name="")
G4ParticleDefinition * GetDefinition() const
virtual G4ThreeVector & SampleDirection(const G4DynamicParticle *dp, G4double kinEnergyFinal, G4int Z, const G4Material *mat=0)
const G4ThreeVector & GetMomentumDirection() const
virtual ~G4DNARuddAngle()
G4double GetPDGMass() const
static G4Electron * Electron()
const G4ParticleDefinition * fElectron
G4ThreeVector fLocalDirection
virtual G4ThreeVector & SampleDirectionForShell(const G4DynamicParticle *dp, G4double kinEnergyFinal, G4int Z, G4int shellIdx, const G4Material *mat=0)
void PrintGeneratorInformation() const