70 0.7628, 0.8983, 0.9801 };
73 0.1569, 0.1112, 0.0506 };
83 limitKinEnergy(100.*
keV),
84 logLimitKinEnergy(
G4Log(limitKinEnergy)),
85 twoln10(2.0*
G4Log(10.0)),
140 if(cutEnergy < maxEnergy) {
143 G4double energy2 = totEnergy*totEnergy;
144 G4double beta2 = kineticEnergy*(kineticEnergy + 2.0*
mass)/energy2;
146 cross = 1.0/cutEnergy - 1.0/maxEnergy - beta2*
G4Log(maxEnergy/cutEnergy)/tmax
147 + 0.5*(maxEnergy - cutEnergy)/energy2;
154 G4double logstep = logtmax - logtmin;
157 for (
G4int ll=0; ll<8; ll++)
162 dcross +=
wgi[ll]*(1.0/ep - beta2/tmax + 0.5*ep/energy2)*a1*(a3 - a1);
188 (p,kineticEnergy,cutEnergy,maxEnergy);
203 (p,kineticEnergy,cutEnergy,maxEnergy);
227 -(1.0 + cutEnergy/tmax)*beta2;
230 G4double del = 0.5*cutEnergy/totEnergy;
242 if (dedx < 0.0) dedx = 0.0 ;
250 G4double ftot2= 0.5/(totEnergy*totEnergy);
252 for (
G4int ll=0; ll<8; ll++)
257 dloss +=
wgi[ll]*(1.0 - beta2*ep/tmax + ep*ep*ftot2)*a1*(a3 - a1);
280 if(minKinEnergy >= maxKinEnergy) {
return; }
284 G4double etot2 = totEnergy*totEnergy;
285 G4double beta2 = kineticEnergy*(kineticEnergy + 2.0*
mass)/etot2;
298 deltaKinEnergy = minKinEnergy*maxKinEnergy
299 /(minKinEnergy*(1.0 - q) + maxKinEnergy*q);
302 f = 1.0 - beta2*deltaKinEnergy/tmax
303 + 0.5*deltaKinEnergy*deltaKinEnergy/etot2;
312 G4cout <<
"G4MuBetheBlochModel::SampleSecondary Warning! "
313 <<
"Majorant " << grej <<
" < "
314 << f <<
" for edelta= " << deltaKinEnergy
315 <<
" tmin= " << minKinEnergy <<
" max= " << maxKinEnergy
323 G4double totalMomentum = totEnergy*sqrt(beta2);
325 (deltaMomentum * totalMomentum);
327 G4double sint = sqrt(1.0 - cost*cost);
331 G4ThreeVector deltaDirection(sint*cos(phi),sint*sin(phi), cost) ;
336 kineticEnergy -= deltaKinEnergy;
338 direction = dir.
unit();
344 deltaDirection,deltaKinEnergy);
345 vdp->push_back(delta);
G4ParticleDefinition * theElectron
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4double HighOrderCorrections(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy, G4double cutEnergy)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
static constexpr double keV
G4double MaxSecondaryKinEnergy(const G4DynamicParticle *dynParticle)
static constexpr double twopi_mc2_rcl2
G4double ShellCorrection(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
const G4ThreeVector & GetMomentumDirection() const
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kineticEnergy, G4double Z, G4double A, G4double cutEnergy, G4double maxEnergy) override
G4double DensityCorrection(G4double x)
Hep3Vector & rotateUz(const Hep3Vector &)
G4double logLimitKinEnergy
void SetProposedMomentumDirection(const G4ThreeVector &dir)
G4double G4Log(G4double x)
G4ParticleChangeForLoss * GetParticleChangeForLoss()
G4EmCorrections * EmCorrections()
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy) override
static constexpr double electron_mass_c2
static constexpr double twopi
virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition *, G4double kinEnergy) override
G4double GetMeanExcitationEnergy() const
virtual G4double MinEnergyCut(const G4ParticleDefinition *, const G4MaterialCutsCouple *) override
static G4Electron * Electron()
void SetParticle(const G4ParticleDefinition *p)
virtual G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy) override
void SetProposedKineticEnergy(G4double proposedKinEnergy)
G4IonisParamMat * GetIonisation() const
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
static G4LossTableManager * Instance()
const G4Material * GetMaterial() const
static constexpr double fine_structure_const
G4double GetElectronDensity() const
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy) override
G4ParticleChangeForLoss * fParticleChange
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4MuBetheBlochModel(const G4ParticleDefinition *p=nullptr, const G4String &nam="MuBetheBloch")
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override