86 {0.03377,0.16940,0.38069,0.61931,0.83060,0.96623};
88 {0.08566,0.18038,0.23396,0.23396,0.18038,0.08566};
100 fParticleChange(nullptr),
101 lowestKinEnergy(1.0*
GeV),
102 minThreshold(0.9*
keV)
112 for(
G4int i=1; i<93; ++i) {
182 const G4double* theAtomicNumDensityVector =
191 dedx += loss*theAtomicNumDensityVector[i];
194 if(dedx < 0.) dedx = 0.;
205 static const G4int k2=5;
213 if(vcut>vmax) { bbb = vmax; }
215 if(kkk < 1) { kkk = 1; }
220 for(
G4int l=0; l<kkk; l++)
222 for(
G4int i=0; i<6; i++)
230 loss *=hhh*totalEnergy ;
244 static const G4int k2 = 4;
247 if(cut >= tkin)
return cross;
255 if(kkk < 1) { kkk = 1; }
261 for(
G4int l=0; l<kkk; l++)
263 for(
G4int i=0; i<6; i++)
288 if(gammaEnergy > tkin) {
return dxsection; }
292 G4double delta = 0.5*mass*mass*v/(E-gammaEnergy) ;
296 if(iz < 1) { iz = 1; }
297 else if(iz > 92) { iz = 92; }
315 (mass+delta*(dnstar*sqrte-2.))) ;
316 if(fn <0.) { fn = 0.; }
320 if(gammaEnergy<epmax1)
328 dxsection =
coeff*(1.-v*(1. - 0.75*v))*Z*(fn*Z +
fe)/gammaEnergy;
347 if (cut >= tmax)
return cross;
350 if(tmax < kineticEnergy) {
359 std::vector<G4DynamicParticle*>* vdp,
370 if(tmin >= tmax)
return;
381 G4double totalMomentum = sqrt(kineticEnergy*(kineticEnergy + 2.0*mass));
412 G4double dirx = sint*cos(phi), diry = sint*sin(phi), dirz = cos(theta) ;
417 partDirection *= totalMomentum;
418 partDirection -= gEnergy*gDirection;
419 partDirection = partDirection.
unit();
422 kineticEnergy -= gEnergy;
429 vdp->push_back(aGamma);
void SetElementSelectors(std::vector< G4EmElementSelector * > *)
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
static constexpr double MeV
virtual void InitialiseLocal(const G4ParticleDefinition *, G4VEmModel *masterModel) override
static constexpr double keV
static const G4double xgi[6]
const G4ThreeVector & GetMomentumDirection() const
const G4ElementVector * GetElementVector() const
Hep3Vector & rotateUz(const Hep3Vector &)
const G4Element * SelectRandomAtom(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
void SetProposedMomentumDirection(const G4ThreeVector &dir)
std::vector< G4EmElementSelector * > * GetElementSelectors()
static const G4double ak1
G4double G4Log(G4double x)
G4ParticleChangeForLoss * GetParticleChangeForLoss()
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kineticEnergy, G4double Z, G4double A, G4double cutEnergy, G4double maxEnergy) override
const G4ParticleDefinition * particle
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy) override
void InitialiseElementSelectors(const G4ParticleDefinition *, const G4DataVector &)
static constexpr double electron_mass_c2
static constexpr double twopi
G4double ComputMuBremLoss(G4double Z, G4double tkin, G4double cut)
static const G4double sqrte
virtual G4double MinEnergyCut(const G4ParticleDefinition *, const G4MaterialCutsCouple *) override
static const G4double wgi[6]
virtual G4double MinPrimaryEnergy(const G4Material *, const G4ParticleDefinition *, G4double) override
virtual G4double ComputeDMicroscopicCrossSection(G4double tkin, G4double Z, G4double gammaEnergy)
std::vector< G4Element * > G4ElementVector
G4double GetA27(G4int Z) const
void SetProposedKineticEnergy(G4double proposedKinEnergy)
G4double GetKineticEnergy() const
G4ParticleDefinition * theGamma
G4double GetZ13(G4double Z) const
const G4double * GetAtomicNumDensityVector() const
void SetParticle(const G4ParticleDefinition *)
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy) override
G4double HighEnergyLimit() const
G4MuBremsstrahlungModel(const G4ParticleDefinition *p=nullptr, const G4String &nam="MuBrem")
static constexpr double GeV
G4ParticleChangeForLoss * fParticleChange
size_t GetNumberOfElements() const
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
static G4NistManager * Instance()
G4double ComputeMicroscopicCrossSection(G4double tkin, G4double Z, G4double cut)