52#ifndef G4UrbanMscModel_h
53#define G4UrbanMscModel_h 1
96 G4double& currentMinimalStep)
override;
116 void InitialiseModelCache();
147 G4double tlimitminfix,tlimitminfix2;
181 G4double coeffc1, coeffc2, coeffc3, coeffc4;
185 static std::vector<mscData*> msc;
193 G4bool latDisplasmentbackup;
208 ChargeSquare = charge*charge;
214inline G4double G4UrbanMscModel::Randomizetlimit()
217 if(tlimit > tlimitmin)
219 res = G4RandGauss::shoot(rndmEngineMod,tlimit,0.1*(tlimit-tlimitmin));
220 res = std::max(res, tlimitmin);
232 G4double a = (2.*xmeanth+9.*x2meanth-3.)/(2.*xmeanth-3.*x2meanth+1.);
238 ? -1.+2.*
G4Exp(
G4Log(rdm)/(a+1.)) : -1.+2.*rdm;
242inline G4double G4UrbanMscModel::ComputeStepmin()
246 G4double rat = currentKinEnergy*invmev;
247 return lambda0*1.e-3/(2.e-3+rat*(msc[idx]->stepmina+msc[idx]->stepminb*rat));
250inline G4double G4UrbanMscModel::ComputeTlimitmin()
252 G4double x = (particle == positron) ?
253 0.7*msc[idx]->sqrtZ*stepmin : 0.87*msc[idx]->Z23*stepmin;
254 if(currentKinEnergy < tlow) { x *= 0.5*(1.+currentKinEnergy/tlow); }
255 return std::max(x, tlimitminfix);
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4double G4Log(G4double x)
G4double GetPDGMass() const
G4double GetPDGCharge() const
G4ThreeVector & SampleScattering(const G4ThreeVector &, G4double safety) override
G4double ComputeTrueStepLength(G4double geomStepLength) override
G4double ComputeTheta0(G4double truePathLength, G4double KineticEnergy)
void StartTracking(G4Track *) override
~G4UrbanMscModel() override
void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
G4double ComputeGeomPathLength(G4double truePathLength) override
G4UrbanMscModel & operator=(const G4UrbanMscModel &right)=delete
G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *particle, G4double KineticEnergy, G4double AtomicNumber, G4double AtomicWeight=0., G4double cut=0., G4double emax=DBL_MAX) override
G4double ComputeTruePathLengthLimit(const G4Track &track, G4double ¤tMinimalStep) override
G4UrbanMscModel(const G4UrbanMscModel &)=delete