54 #ifndef G4UrbanAdjointMscModel_h
55 #define G4UrbanAdjointMscModel_h 1
101 G4double& currentMinimalStep)
override;
121 inline void UpdateCache();
150 G4double tlimitminfix,tlimitminfix2;
174 G4int currentMaterialIndex;
179 G4double coeffc1,coeffc2,coeffc3,coeffc4;
184 G4bool latDisplasmentbackup ;
199 displacementFlag = val;
208 if (p1 != particle) {
212 ChargeSquare = charge*charge;
218 inline G4double G4UrbanAdjointMscModel::Randomizetlimit()
221 if(tlimit > tlimitmin)
227 }
while ((temptlimit < tlimit-delta) ||
228 (temptlimit > tlimit+delta));
230 else { temptlimit = tlimitmin; }
237 inline void G4UrbanAdjointMscModel::UpdateCache()
242 G4double facz = 0.990395+w*(-0.168386+w*0.093286) ;
243 coeffth1 = facz*(1. - 8.7780e-2/Zeff);
244 coeffth2 = facz*(4.0780e-2 + 1.7315e-4*Zeff);
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.);
271 if(rndmEngineMod->
flat() < prob) {
274 cth = -1.+2.*rndmEngineMod->
flat();
ThreeVector shoot(const G4int Ap, const G4int Af)
virtual G4double ComputeTruePathLengthLimit(const G4Track &track, G4double ¤tMinimalStep) override
std::vector< ExP01TrackerHit * > a
G4UrbanAdjointMscModel(const G4String &nam="UrbanMsc")
virtual ~G4UrbanAdjointMscModel()
void SetNewDisplacementFlag(G4bool)
const G4String & GetParticleName() const
virtual G4double ComputeTrueStepLength(G4double geomStepLength) override
virtual void StartTracking(G4Track *) override
static constexpr double eplus
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *particle, G4double KineticEnergy, G4double AtomicNumber, G4double AtomicWeight=0., G4double cut=0., G4double emax=DBL_MAX) override
static const G4double emax
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4double GetPDGMass() const
G4double ComputeTheta0(G4double truePathLength, G4double KineticEnergy)
virtual G4ThreeVector & SampleScattering(const G4ThreeVector &, G4double safety) override
static G4Electron * Electron()
virtual G4double ComputeGeomPathLength(G4double truePathLength) override
G4double GetPDGCharge() const
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override