79 singleScatteringMode(false)
143 G4Exception(
"G4WentzelVIRelModel::ComputeCrossSectionPerAtom",
"em0011",
176 G4double tlimit = currentMinimalStep;
241 if(rcut > rlimit) { rlimit =
std::min(rlimit, rcut*sqrt(rlimit/rcut)); }
243 if(rlimit < tlimit) { tlimit = rlimit; }
318 static const G4double singleScatLimit = 1.0e-7;
365 }
else if(
xtsec > 0.0) {
415 static const G4double thinlimit = 0.1;
471 for (; i<nelm; ++i) {
if(
xsecn[i] >= qsec) {
break; } }
485 if(x2 <= 0.0) { --nMscSteps; }
503 if(z0 > 0.01) { zzz =
G4Exp(-1.0/z0); }
509 if(cost > 1.0) { cost = 1.0; }
510 else if(cost < -1.0) { cost =-1.0; }
511 sint = sqrt((1.0 - cost)*(1.0 + cost));
517 temp.set(vx1,vy1,cost);
529 if(d >= 0.0) { dz = sqrt(d) - r; }
530 else { dx = dy = dz = 0.0; }
539 }
while (0 < nMscSteps);
541 dir.rotateUz(oldDirection);
572 const G4double* theAtomNumDensityVector =
589 for (
G4int i=0; i<nelm; ++i) {
605 if(nucsec > 0.0) { esec /= nucsec; }
static G4Pow * GetInstance()
void StartTracking(G4Track *)
ThreeVector shoot(const G4int Ap, const G4int Af)
G4double ComputeNuclearCrossSection(G4double CosThetaMin, G4double CosThetaMax)
virtual G4double ComputeTruePathLengthLimit(const G4Track &track, G4double ¤tMinimalStep)
G4ThreeVector SampleSingleScattering(G4double CosThetaMin, G4double CosThetaMax, G4double elecRatio=0.0)
static G4LossTableManager * Instance()
virtual G4ThreeVector & SampleScattering(const G4ThreeVector &, G4double safety)
std::vector< G4Element * > G4ElementVector
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
const G4DynamicParticle * GetDynamicParticle() const
G4double GetProductionCut(G4int index) const
G4StepStatus GetStepStatus() const
G4double ConvertTrueToGeom(G4double &tLength, G4double &gLength)
G4WentzelVIRelModel(const G4String &nam="WentzelVIUni")
const G4MaterialCutsCouple * GetMaterialCutsCouple() const
virtual G4double ComputeGeomPathLength(G4double truePathLength)
G4ParticleDefinition * GetDefinition() const
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double KineticEnergy, G4double AtomicNumber, G4double AtomicWeight=0., G4double cut=DBL_MAX, G4double emax=DBL_MAX)
G4WentzelVIRelXSection * wokvi
G4double ComputeTransportCrossSectionPerAtom(G4double CosThetaMax)
const G4Step * GetStep() const
const G4ElementVector * GetElementVector() const
static G4NistManager * Instance()
G4double SetupKinematic(G4double kinEnergy, const G4Material *mat)
G4double SetupTarget(G4int Z, G4double cut=DBL_MAX)
void DefineMaterial(const G4MaterialCutsCouple *)
G4double ComputeXSectionPerVolume()
G4bool singleScatteringMode
G4StepPoint * GetPreStepPoint() const
G4double GetEnergy(const G4ParticleDefinition *part, G4double range, const G4MaterialCutsCouple *couple)
const G4double * GetVecNbOfAtomsPerVolume() const
const G4MaterialCutsCouple * CurrentCouple() const
const G4ThreeVector & GetPosition() const
G4double GetRange(const G4ParticleDefinition *part, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
G4ParticleChangeForMSC * GetParticleChangeForMSC(const G4ParticleDefinition *p=0)
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
std::vector< G4double > prob
G4double ComputeGeomLimit(const G4Track &, G4double &presafety, G4double limit)
G4ThreeVector fDisplacement
G4double GetTransportMeanFreePath(const G4ParticleDefinition *part, G4double kinEnergy)
G4double GetRadlen() const
const G4DataVector * currentCuts
G4NistManager * fNistManager
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
G4double ComputeSafety(const G4ThreeVector &position, G4double limit)
G4int currentMaterialIndex
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4LossTableManager * theManager
void ProposeMomentumDirection(const G4ThreeVector &Pfinal)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
virtual ~G4WentzelVIRelModel()
virtual G4double ComputeTrueStepLength(G4double geomStepLength)
G4double PolarAngleLimit() const
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double GetSafety() const
G4MscStepLimitType steppingAlgorithm
std::vector< G4double > xsecn
G4ParticleChangeForMSC * fParticleChange
size_t GetNumberOfElements() const
const G4ParticleDefinition * particle
G4double ComputeElectronCrossSection(G4double CosThetaMin, G4double CosThetaMax)
const G4MaterialCutsCouple * currentCouple
void Initialise(const G4ParticleDefinition *, G4double CosThetaLim)
void SetupParticle(const G4ParticleDefinition *)
const G4Material * currentMaterial
G4ProductionCuts * GetProductionCuts() const