54 #ifndef G4UrbanMscModel_h
55 #define G4UrbanMscModel_h 1
99 G4double& currentMinimalStep)
override;
222 G4double facz = 0.990395+w*(-0.168386+w*0.093286) ;
228 coeffc1 = 2.3785 - Z13*(4.1981e-1 - Z13*6.3100e-2);
229 coeffc2 = 4.7526e-1 + Z13*(1.7694 - Z13*3.3885e-1);
230 coeffc3 = 2.3683e-1 - Z13*(1.8111 - Z13*3.2774e-1);
231 coeffc4 = 1.7888e-2 + Z13*(1.9659e-2 - Z13*2.6664e-3);
246 G4double a = (2.*xmeanth+9.*x2meanth-3.)/(2.*xmeanth-3.*x2meanth+1.);
252 ? -1.+2.*
G4Exp(
G4Log(rdm)/(a+1.)) : -1.+2.*rdm;
G4bool latDisplasmentbackup
G4UrbanMscModel(const G4String &nam="UrbanMsc")
G4double Randomizetlimit()
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
T max(const T t1, const T t2)
brief Return the largest of the two arguments
std::vector< ExP01TrackerHit * > a
const G4ParticleDefinition * positron
virtual G4ThreeVector & SampleScattering(const G4ThreeVector &, G4double safety) override
virtual void StartTracking(G4Track *) override
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
G4LossTableManager * theManager
void SampleDisplacement(G4double sinTheta, G4double phi)
G4double ComputeTheta0(G4double truePathLength, G4double KineticEnergy)
G4double GetPDGCharge() const
static const G4double emax
G4double GetPDGMass() const
G4double G4Log(G4double x)
const G4MaterialCutsCouple * couple
G4double SampleCosineTheta(G4double trueStepLength, G4double KineticEnergy)
G4ParticleChangeForMSC * fParticleChange
virtual ~G4UrbanMscModel()
void SetParticle(const G4ParticleDefinition *)
virtual G4double ComputeGeomPathLength(G4double truePathLength) override
CLHEP::HepRandomEngine * rndmEngineMod
ThreeVector shoot(const G4int Ap, const G4int Af)
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *particle, G4double KineticEnergy, G4double AtomicNumber, G4double AtomicWeight=0., G4double cut=0., G4double emax=DBL_MAX) override
G4double SimpleScattering(G4double xmeanth, G4double x2meanth)
G4double currentRadLength
virtual G4double ComputeTruePathLengthLimit(const G4Track &track, G4double ¤tMinimalStep) override
G4double currentKinEnergy
void SampleDisplacementNew(G4double sinTheta, G4double phi)
G4int currentMaterialIndex
G4UrbanMscModel & operator=(const G4UrbanMscModel &right)=delete
virtual G4double ComputeTrueStepLength(G4double geomStepLength) override
const G4ParticleDefinition * particle
static constexpr double eplus