80 beta2lim(betalim*betalim),
81 bg2lim(beta2lim*(1.0 + beta2lim))
91 G4cout <<
"### Monopole ionisation model with d-electron production, Gmag= "
127 if(!
dedx0) {
dedx0 =
new std::vector<G4double>; }
132 if(n < numOfCouples) {
dedx0->resize(numOfCouples); }
135 for(
G4int i=0; i<numOfCouples; ++i) {
185 dedx = (kapa1*dedx1 + kapa2*dedx2)/(kapa1 + kapa2);
207 if(
nmpl > 1) { k = 0.346; }
210 const G4double B[7] = { 0.0, 0.248, 0.672, 1.022, 1.243, 1.464, 1.685};
212 dedx += 0.5 * k - B[
nmpl];
221 if (dedx < 0.0) { dedx = 0.; }
239 if(cutEnergy < maxEnergy) {
273 if(minKinEnergy >= maxKinEnergy) {
return; }
280 G4double etot2 = totEnergy*totEnergy;
281 G4double beta2 = kineticEnergy*(kineticEnergy + 2.0*
mass)/etot2;
285 G4double deltaKinEnergy = minKinEnergy*maxKinEnergy
286 /(minKinEnergy*(1.0 - q) + maxKinEnergy*q);
289 G4double totMomentum = totEnergy*sqrt(beta2);
293 (deltaMomentum * totMomentum);
294 if(cost > 1.0) { cost = 1.0; }
296 G4double sint = sqrt((1.0 - cost)*(1.0 + cost));
300 G4ThreeVector deltaDirection(sint*cos(phi),sint*sin(phi), cost);
308 vdp->push_back(delta);
311 kineticEnergy -= deltaKinEnergy;
312 G4ThreeVector finalP = direction*totMomentum - deltaDirection*deltaMomentum;
313 finalP = finalP.
unit();
331 G4double twomeanLoss = meanLoss + meanLoss;
333 if(twomeanLoss < siga) {
337 x = (loss - meanLoss)/siga;
344 }
while (0.0 > loss || loss > twomeanLoss);
T max(const T t1, const T t2)
brief Return the largest of the two arguments
virtual G4double SampleFluctuations(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmax, G4double length, G4double meanLoss) override
G4double pi_hbarc2_over_mc2
void SetHighEnergyLimit(G4double)
virtual G4double Dispersion(const G4Material *, const G4DynamicParticle *, G4double tmax, G4double length) override
virtual G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
static constexpr double twopi_mc2_rcl2
const G4ThreeVector & GetMomentumDirection() const
static constexpr double hbarc
G4double DensityCorrection(G4double x)
Hep3Vector & rotateUz(const Hep3Vector &)
static const G4double emax
void SetProposedMomentumDirection(const G4ThreeVector &dir)
G4double GetPDGMass() const
G4double G4Log(G4double x)
G4double LowEnergyLimit() const
G4ParticleChangeForLoss * fParticleChange
static constexpr double electron_Compton_length
G4ParticleChangeForLoss * GetParticleChangeForLoss()
static constexpr double g
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy) override
void SetParticle(const G4ParticleDefinition *p)
G4ParticleDefinition * theElectron
G4ParticleDefinition * GetDefinition() const
static constexpr double cm2
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
static constexpr double electron_mass_c2
virtual ~G4mplIonisationWithDeltaModel()
static std::vector< G4double > * dedx0
static constexpr double twopi
G4double GetMeanExcitationEnergy() const
static G4Electron * Electron()
size_t GetTableSize() const
ThreeVector shoot(const G4int Ap, const G4int Af)
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kineticEnergy, G4double Z, G4double A, G4double cutEnergy, G4double maxEnergy) override
const G4MaterialCutsCouple * CurrentCouple() const
static G4ProductionCutsTable * GetProductionCutsTable()
static constexpr double eplus
const G4ParticleDefinition * monopole
void SetLowEnergyLimit(G4double)
void SetProposedKineticEnergy(G4double proposedKinEnergy)
G4IonisParamMat * GetIonisation() const
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
const G4Material * GetMaterial() const
static constexpr double pi
virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition *, G4double kinEnergy) override
static constexpr double fine_structure_const
G4double HighEnergyLimit() const
double B(double temperature)
static constexpr double GeV
G4double GetElectronDensity() const
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double ComputeDEDXAhlen(const G4Material *material, G4double bg2, G4double cut)
G4mplIonisationWithDeltaModel(G4double mCharge, const G4String &nam="mplIonisationWithDelta")
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy) override