56 LowestKineticEnergy(1*
GeV), HighestKineticEnergy(1*
PeV),
96 for (
G4int j = 0; j < nEl; j++) {
97 Z =
G4lrint((*theElementTable)[j]->GetZ());
119 {0.0199,0.1017,0.2372,0.4083,0.5917,0.7628,0.8983,0.9801};
121 {0.0506,0.1112,0.1569,0.1813,0.1813,0.1569,0.1112,0.0506};
128 if (KineticEnergy <=
CutFixed)
return CrossSection;
131 G4double epmax = KineticEnergy + Mass - 0.5*proton_mass_c2;
132 if (epmax <= epmin)
return CrossSection;
142 for (
G4int l = 0; l < kkk; ++l) {
144 for (
G4int ll = 0; ll < 8; ++ll) {
152 CrossSection *= hhh ;
153 if (CrossSection < 0.) { CrossSection = 0.; }
167 static const G4double coeffn = fine_structure_const/
pi;
170 G4double TotalEnergy = KineticEnergy + ParticleMass;
174 if ((epsilon >= TotalEnergy - 0.5*proton_mass_c2) ||
175 (epsilon <=
CutFixed) ) {
return DCrossSection; }
184 G4double mass2 = ParticleMass*ParticleMass;
186 G4double up = TotalEnergy*TotalEnergy*v1/mass2*(1.+mass2*v2/(alam2*v1));
187 G4double down = 1.+epsilon/alam*(1.+alam/(2.*proton_mass_c2)+epsilon/alam);
189 DCrossSection = coeffn*aeff*sigph/epsilon*
190 (-v1+(v1+0.5*v2*(1.+2.*mass2/alam2))*
G4Log(up/down));
192 if (DCrossSection < 0.) { DCrossSection = 0.; }
193 return DCrossSection;
G4double HighestKineticEnergy
G4double GetKineticEnergy() const
static const double microbarn
G4double ComputeDDMicroscopicCrossSection(G4double incidentKE, G4double Z, G4double A, G4double epsilon)
G4double ComputeMicroscopicCrossSection(G4double incidentKE, G4double A)
static G4NistManager * Instance()
static size_t GetNumberOfElements()
G4double GetElementCrossSection(const G4DynamicParticle *particle, G4int Z, const G4Material *)
void PutValue(size_t index, G4double theValue)
G4double Energy(size_t index) const
G4bool IsElementApplicable(const G4DynamicParticle *particle, G4int Z, const G4Material *)
G4double Value(G4double theEnergy, size_t &lastidx) const
static const G4double A[nN]
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
virtual ~G4KokoulinMuonNuclearXS()
G4double GetPDGMass() const
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double energy(const ThreeVector &p, const G4double m)
G4KokoulinMuonNuclearXS()
G4double GetAtomicMassAmu(const G4String &symb) const
void BuildCrossSectionTable()
static G4MuonMinus * MuonMinus()
void BuildPhysicsTable(const G4ParticleDefinition &)
static G4PhysicsVector * theCrossSection[MAXZMUN]
std::vector< G4Element * > G4ElementTable
static G4ElementTable * GetElementTable()
G4double LowestKineticEnergy
G4bool isInitialized()
Check if the generator is initialized.