120 if(!
dedx0) {
dedx0 =
new std::vector<G4double>; }
125 if(n < numOfCouples) {
dedx0->resize(numOfCouples); }
128 for(
G4int i=0; i<numOfCouples; ++i) {
133 G4double vF = electron_Compton_length*pow(3*
pi*
pi*eDensity,0.3333333333);
135 (
G4Log(2*vF/fine_structure_const) - 0.5)/vF;
174 dedx = (kapa1*dedx1 + kapa2*dedx2)/(kapa1 + kapa2);
194 G4double dedx = log(2.0 * electron_mass_c2 * bg2 / eexc) - 0.5;
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);
273 siga = (invbeta2 - 0.5) * twopi_mc2_rcl2 * tmax * length
ThreeVector shoot(const G4int Ap, const G4int Af)
G4IonisParamMat * GetIonisation() const
G4double LowEnergyLimit() const
virtual G4double Dispersion(const G4Material *, const G4DynamicParticle *, G4double tmax, G4double length)
G4double GetAdensity() const
G4ParticleChangeForLoss * GetParticleChangeForLoss()
G4double GetKineticEnergy() const
G4double HighEnergyLimit() const
G4double GetX1density() const
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
double B(double temperature)
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
static const G4double bg2lim
G4mplIonisationModel(G4double mCharge, const G4String &nam="mplIonisation")
void SetHighEnergyLimit(G4double)
G4double pi_hbarc2_over_mc2
const G4MaterialCutsCouple * CurrentCouple() const
size_t GetTableSize() const
G4double GetElectronDensity() const
G4double ComputeDEDXAhlen(const G4Material *material, G4double bg2)
static const G4double beta2lim
void SetParticle(const G4ParticleDefinition *p)
virtual ~G4mplIonisationModel()
static const G4double emax
G4double G4Log(G4double x)
static G4ProductionCutsTable * GetProductionCutsTable()
G4double GetX0density() const
G4double GetPDGMass() const
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
G4double GetCdensity() const
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4ParticleChangeForLoss * fParticleChange
const G4double x[NPOINTSGL]
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double GetMeanExcitationEnergy() const
const G4ParticleDefinition * monopole
static const G4double twoln10
static std::vector< G4double > * dedx0
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy)
void SetLowEnergyLimit(G4double)
G4double GetMdensity() const
virtual G4double SampleFluctuations(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmax, G4double length, G4double meanLoss)
const G4Material * GetMaterial() const