54 #ifndef G4UrbanAdjointMscModel_h
55 #define G4UrbanAdjointMscModel_h 1
101 G4double& currentMinimalStep)
override;
227 }
while ((temptlimit <
tlimit-delta) ||
228 (temptlimit >
tlimit+delta));
242 G4double facz = 0.990395+w*(-0.168386+w*0.093286) ;
248 coeffc1 = 2.3785 - Z13*(4.1981e-1 - Z13*6.3100e-2);
249 coeffc2 = 4.7526e-1 + Z13*(1.7694 - Z13*3.3885e-1);
250 coeffc3 = 2.3683e-1 - Z13*(1.8111 - Z13*3.2774e-1);
251 coeffc4 = 1.7888e-2 + Z13*(1.9659e-2 - Z13*2.6664e-3);
266 G4double a = (2.*xmeanth+9.*x2meanth-3.)/(2.*xmeanth-3.*x2meanth+1.);
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 G4Exp(G4double initial_x)
Exponential Function double precision.
G4double SampleCosineTheta(G4double trueStepLength, G4double KineticEnergy)
std::vector< ExP01TrackerHit * > a
G4LossTableManager * theManager
G4ParticleChangeForMSC * fParticleChange
const G4ParticleDefinition * particle
const G4ParticleDefinition * positron
virtual G4double ComputeGeomPathLength(G4double truePathLength) override
void SetNewDisplacementFlag(G4bool)
G4double currentKinEnergy
const G4String & GetParticleName() const
G4UrbanAdjointMscModel(const G4String &nam="UrbanMsc")
G4double GetPDGCharge() const
virtual ~G4UrbanAdjointMscModel()
static const G4double emax
virtual G4double ComputeTruePathLengthLimit(const G4Track &track, G4double ¤tMinimalStep) override
G4double GetPDGMass() const
G4double G4Log(G4double x)
virtual void StartTracking(G4Track *) override
CLHEP::HepRandomEngine * rndmEngineMod
G4double currentRadLength
void SampleDisplacementNew(G4double sinTheta, G4double phi)
static G4Electron * Electron()
ThreeVector shoot(const G4int Ap, const G4int Af)
G4int currentMaterialIndex
G4double Randomizetlimit()
void SetParticle(const G4ParticleDefinition *)
G4UrbanAdjointMscModel & operator=(const G4UrbanAdjointMscModel &right)=delete
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
G4double ComputeTheta0(G4double truePathLength, G4double KineticEnergy)
const G4MaterialCutsCouple * couple
virtual G4ThreeVector & SampleScattering(const G4ThreeVector &, G4double safety) override
G4bool latDisplasmentbackup
virtual G4double ComputeTrueStepLength(G4double geomStepLength) override
void SampleDisplacement(G4double sinTheta, G4double phi)
static constexpr double eplus