71 twoln10(2.0*
G4Log(10.0))
211 return material->GetElectronDensity()
247 if(eloss >= preKinEnergy) {
return; }
268 res += (
G4Log(x)*(tau + 1.)*(tau + 1.)/(tau * (tau + 2.0)) + 1.0 - x)
290 if(res > preKinEnergy) { res = preKinEnergy; }
291 else if(res < 0.0) { res = eloss; }
301 vector<G4DynamicParticle*>* vdp,
312 if(minKinEnergy >= maxKinEnergy) {
return; }
318 G4double etot2 = totEnergy*totEnergy;
319 G4double beta2 = kineticEnergy*(kineticEnergy + 2.0*
mass)/etot2;
324 if( 0.0 <
spin ) { fmax += 0.5*maxKinEnergy*maxKinEnergy/etot2; }
332 deltaKinEnergy = minKinEnergy*maxKinEnergy
333 /(minKinEnergy*(1.0 - rndm[0]) + maxKinEnergy*rndm[0]);
335 f = 1.0 - beta2*deltaKinEnergy/tmax;
337 f1 = 0.5*deltaKinEnergy*deltaKinEnergy/etot2;
342 }
while( fmax*rndm[1] > f);
354 grej *= (1.0 +
magMoment2*(x2 - f1/f)/(1.0 + x2));
357 G4cout <<
"### G4LindhardSorensenIonModel WARNING: grej= " << grej
359 <<
" Ekin(MeV)= " << kineticEnergy
360 <<
" delEkin(MeV)= " << deltaKinEnergy
363 if(rndmEngineMod->
flat() > grej) {
return; }
383 G4double sint = std::sqrt((1.0 - cost)*(1.0 + cost));
387 deltaDirection.
set(sint*std::cos(phi),sint*std::sin(phi), cost) ;
405 vdp->push_back(delta);
408 kineticEnergy -= deltaKinEnergy;
410 finalP = finalP.
unit();
static const G4double emax
G4double G4Log(G4double x)
static constexpr double GeV
G4GLOB_DLL std::ostream G4cout
void set(double x, double y, double z)
Hep3Vector & rotateUz(const Hep3Vector &)
virtual void flatArray(const int size, double *vect)=0
void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy) override
virtual G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy) override
const G4ThreeVector & GetMomentumDirection() const
G4ParticleDefinition * GetDefinition() const
G4double GetKineticEnergy() const
G4ThreeVector GetMomentum() const
G4double GetTotalMomentum() const
static G4Electron * Electron()
G4double EffectiveChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
G4double BarkasCorrection(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
G4double GetParticleCharge(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
G4double EffectiveChargeCorrection(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
G4double GetDEDX(const G4Material *, const G4int Z, const G4double e, const G4double loge) const
G4double GetMeanExcitationEnergy() const
G4double GetDeltaL(G4int Z, G4double gamma) const
G4BetheBlochModel * fBBModel
~G4LindhardSorensenIonModel() override
static std::vector< G4float > * fact[MAXZION]
G4double MaxSecondaryEnergy(const G4ParticleDefinition *, G4double kinEnergy) override
G4ParticleChangeForLoss * fParticleChange
static G4LindhardSorensenData * lsdata
static G4IonICRU73Data * fIonData
G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kineticEnergy, G4double Z, G4double A, G4double cutEnergy, G4double maxEnergy) override
G4BraggIonModel * fBraggModel
G4double GetChargeSquareRatio() const
G4double MinEnergyCut(const G4ParticleDefinition *, const G4MaterialCutsCouple *couple) override
G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
void SetParticle(const G4ParticleDefinition *p)
void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
G4LindhardSorensenIonModel(const G4ParticleDefinition *p=nullptr, const G4String &nam="LindhardSorensen")
void CorrectionsAlongStep(const G4MaterialCutsCouple *couple, const G4DynamicParticle *dp, const G4double &length, G4double &eloss) override
const G4ParticleDefinition * particle
G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy) override
G4ParticleDefinition * theElectron
G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy) override
G4double GetParticleCharge(const G4ParticleDefinition *p, const G4Material *mat, G4double kineticEnergy) override
void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy) override
static G4LossTableManager * Instance()
G4EmCorrections * EmCorrections()
const G4Material * GetMaterial() const
G4ProductionCuts * GetProductionCuts() const
G4IonisParamMat * GetIonisation() const
G4double GetElectronDensity() const
G4double GetA27(G4int Z) const
static G4NistManager * Instance()
void SetProposedKineticEnergy(G4double proposedKinEnergy)
void SetProposedMomentumDirection(const G4ThreeVector &dir)
G4double GetPDGMagneticMoment() const
G4int GetAtomicNumber() const
G4double GetPDGMass() const
G4double GetPDGCharge() const
const G4String & GetParticleName() const
G4double GetPDGSpin() const
G4double GetProductionCut(G4int index) const
virtual G4ThreeVector & SampleDirection(const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, const G4Material *)=0
virtual void SetParticleAndCharge(const G4ParticleDefinition *, G4double q2)
G4VEmFluctuationModel * GetModelOfFluctuations()
G4VEmAngularDistribution * GetAngularDistribution()
G4int SelectRandomAtomNumber(const G4Material *)
void SetDeexcitationFlag(G4bool val)
void SetAngularDistribution(G4VEmAngularDistribution *)
G4bool UseAngularGeneratorFlag() const
G4ParticleChangeForLoss * GetParticleChangeForLoss()
static constexpr double eplus
static constexpr double electron_mass_c2
static constexpr double proton_mass_c2
static constexpr double GeV
static constexpr double hbar_Planck
static constexpr double twopi
static constexpr double MeV
static constexpr double c_squared
static constexpr double twopi_mc2_rcl2
T max(const T t1, const T t2)
brief Return the largest of the two arguments
T min(const T t1, const T t2)
brief Return the smallest of the two arguments