77 beta2lim(betalim*betalim),
78 bg2lim(beta2lim*(1.0 + beta2lim))
120 if(!
dedx0) {
dedx0 =
new std::vector<G4double>; }
125 if(n < numOfCouples) {
dedx0->resize(numOfCouples); }
128 for(
G4int i=0; i<numOfCouples; ++i) {
174 dedx = (kapa1*dedx1 + kapa2*dedx2)/(kapa1 + kapa2);
198 if(
nmpl > 1) k = 0.346;
201 const G4double B[7] = { 0.0, 0.248, 0.672, 1.022, 1.243, 1.464, 1.685};
203 dedx += 0.5 * k - B[
nmpl];
210 if ( x < x1den ) deltam += aden * pow((x1den-x), mden);
211 dedx -= 0.5 * deltam;
217 if (dedx < 0.0) dedx = 0.;
242 G4double twomeanLoss = meanLoss + meanLoss;
244 if(twomeanLoss < siga) {
248 x = (loss - meanLoss)/siga;
255 }
while (0.0 > loss || loss > twomeanLoss);
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double pi_hbarc2_over_mc2
void SetHighEnergyLimit(G4double)
static constexpr double twopi_mc2_rcl2
static constexpr double hbarc
virtual G4double SampleFluctuations(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmax, G4double length, G4double meanLoss) override
G4double GetX1density() const
G4mplIonisationModel(G4double mCharge, const G4String &nam="mplIonisation")
static const G4double emax
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy) override
G4double GetPDGMass() const
G4double G4Log(G4double x)
G4double LowEnergyLimit() const
static constexpr double electron_Compton_length
G4ParticleChangeForLoss * GetParticleChangeForLoss()
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy) override
static constexpr double g
static constexpr double cm2
void SetParticle(const G4ParticleDefinition *p)
static constexpr double electron_mass_c2
G4ParticleChangeForLoss * fParticleChange
G4double GetAdensity() const
G4double GetMeanExcitationEnergy() const
size_t GetTableSize() const
ThreeVector shoot(const G4int Ap, const G4int Af)
G4double GetX0density() const
const G4MaterialCutsCouple * CurrentCouple() const
const G4ParticleDefinition * monopole
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
static G4ProductionCutsTable * GetProductionCutsTable()
virtual G4double Dispersion(const G4Material *, const G4DynamicParticle *, G4double tmax, G4double length) override
G4double GetCdensity() const
void SetLowEnergyLimit(G4double)
G4IonisParamMat * GetIonisation() const
G4double GetMdensity() const
static std::vector< G4double > * dedx0
G4double GetKineticEnergy() const
virtual ~G4mplIonisationModel()
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
const G4Material * GetMaterial() const
static constexpr double pi
static constexpr double fine_structure_const
G4double HighEnergyLimit() const
double B(double temperature)
static constexpr double GeV
G4double GetElectronDensity() const
G4double ComputeDEDXAhlen(const G4Material *material, G4double bg2)
T min(const T t1, const T t2)
brief Return the smallest of the two arguments