70 0.7628, 0.8983, 0.9801 };
73 0.1569, 0.1112, 0.0506 };
83 limitKinEnergy(100.*
keV),
84 logLimitKinEnergy(
G4Log(limitKinEnergy)),
88 alphaprime(fine_structure_const/
twopi)
120 G4double tmax = 2.0*electron_mass_c2*tau*(tau + 2.) /
145 if(cutEnergy < maxEnergy) {
148 G4double energy2 = totEnergy*totEnergy;
149 G4double beta2 = kineticEnergy*(kineticEnergy + 2.0*
mass)/energy2;
151 cross = 1.0/cutEnergy - 1.0/maxEnergy - beta2*
G4Log(maxEnergy/cutEnergy)/tmax
152 + 0.5*(maxEnergy - cutEnergy)/energy2;
159 G4double logstep = logtmax - logtmin;
162 for (
G4int ll=0; ll<8; ll++)
167 dcross +=
wgi[ll]*(1.0/ep - beta2/tmax + 0.5*ep/energy2)*a1*(a3 - a1);
173 cross *= twopi_mc2_rcl2/beta2;
193 (p,kineticEnergy,cutEnergy,maxEnergy);
208 (p,kineticEnergy,cutEnergy,maxEnergy);
231 G4double dedx =
G4Log(2.0*electron_mass_c2*bg2*cutEnergy/eexc2)
232 -(1.0 + cutEnergy/tmax)*beta2;
235 G4double del = 0.5*cutEnergy/totEnergy;
247 if (dedx < 0.0) dedx = 0.0 ;
255 G4double ftot2= 0.5/(totEnergy*totEnergy);
257 for (
G4int ll=0; ll<8; ll++)
262 dloss +=
wgi[ll]*(1.0 - beta2*ep/tmax + ep*ep*ftot2)*a1*(a3 - a1);
267 dedx *= twopi_mc2_rcl2*eDensity/beta2;
285 if(minKinEnergy >= maxKinEnergy) {
return; }
289 G4double etot2 = totEnergy*totEnergy;
290 G4double beta2 = kineticEnergy*(kineticEnergy + 2.0*
mass)/etot2;
303 deltaKinEnergy = minKinEnergy*maxKinEnergy
304 /(minKinEnergy*(1.0 - q) + maxKinEnergy*q);
307 f = 1.0 - beta2*deltaKinEnergy/tmax
308 + 0.5*deltaKinEnergy*deltaKinEnergy/etot2;
317 G4cout <<
"G4MuBetheBlochModel::SampleSecondary Warning! "
318 <<
"Majorant " << grej <<
" < "
319 << f <<
" for edelta= " << deltaKinEnergy
320 <<
" tmin= " << minKinEnergy <<
" max= " << maxKinEnergy
327 sqrt(deltaKinEnergy * (deltaKinEnergy + 2.0*electron_mass_c2));
328 G4double totalMomentum = totEnergy*sqrt(beta2);
329 G4double cost = deltaKinEnergy * (totEnergy + electron_mass_c2) /
330 (deltaMomentum * totalMomentum);
332 G4double sint = sqrt(1.0 - cost*cost);
336 G4ThreeVector deltaDirection(sint*cos(phi),sint*sin(phi), cost) ;
338 deltaDirection.rotateUz(direction);
341 kineticEnergy -= deltaKinEnergy;
342 G4ThreeVector dir = totalMomentum*direction - deltaMomentum*deltaDirection;
343 direction = dir.unit();
349 deltaDirection,deltaKinEnergy);
350 vdp->push_back(delta);
G4ParticleChangeForLoss * fParticleChange
G4double MaxSecondaryKinEnergy(const G4DynamicParticle *dynParticle)
G4IonisParamMat * GetIonisation() const
G4double HighOrderCorrections(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy, G4double cutEnergy)
static G4LossTableManager * Instance()
G4ParticleChangeForLoss * GetParticleChangeForLoss()
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
virtual G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
virtual G4double MinEnergyCut(const G4ParticleDefinition *, const G4MaterialCutsCouple *)
static const G4double bg2lim
G4ParticleDefinition * theElectron
G4GLOB_DLL std::ostream G4cout
G4double GetElectronDensity() const
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kineticEnergy, G4double Z, G4double A, G4double cutEnergy, G4double maxEnergy)
G4EmCorrections * EmCorrections()
const G4ThreeVector & GetMomentumDirection() const
static const double twopi
void SetProposedKineticEnergy(G4double proposedKinEnergy)
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy)
void SetProposedMomentumDirection(const G4ThreeVector &dir)
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4double ShellCorrection(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
G4double DensityCorrection(G4double x)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
const G4double x[NPOINTSGL]
virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition *, G4double kinEnergy)
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double GetMeanExcitationEnergy() const
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
static G4Electron * Electron()
static const G4double twoln10
G4MuBetheBlochModel(const G4ParticleDefinition *p=0, const G4String &nam="MuBetheBloch")
void SetParticle(const G4ParticleDefinition *p)
G4double logLimitKinEnergy
virtual ~G4MuBetheBlochModel()
const G4Material * GetMaterial() const