79 else if (secKinetic <= 200.*
eV)
88 cosTheta = std::sqrt(1.-sin2O);
100 if (secKinetic>100*
eV) cosTheta = std::sqrt(secKinetic / maxSecKinetic);
104 G4double sint = sqrt((1.0 - cosTheta)*(1.0 + cosTheta));
void set(double x, double y, double z)
const G4ThreeVector & GetMomentumDirection() const
virtual G4ThreeVector & SampleDirection(const G4DynamicParticle *dp, G4double kinEnergyFinal, G4int Z, const G4Material *mat=0)
Hep3Vector & rotateUz(const Hep3Vector &)
G4double GetPDGMass() const
G4ParticleDefinition * GetDefinition() const
static constexpr double electron_mass_c2
virtual G4ThreeVector & SampleDirectionForShell(const G4DynamicParticle *dp, G4double kinEnergyFinal, G4int Z, G4int shellIdx, const G4Material *mat=0)
static constexpr double twopi
static constexpr double eV
static G4Electron * Electron()
void PrintGeneratorInformation() const
G4ThreeVector fLocalDirection
G4double GetKineticEnergy() const
virtual ~G4DNABornAngle()
G4DNABornAngle(const G4String &name="")
const G4ParticleDefinition * fElectron