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};
101 lowestKinEnergy(1.0*
GeV),
102 minThreshold(0.9*
keV)
112 for(
G4int i=1; i<93; ++i) {
187 const G4double* theAtomicNumDensityVector =
196 dedx += loss*theAtomicNumDensityVector[i];
199 if(dedx < 0.) dedx = 0.;
210 static const G4int k2=5;
218 if(vcut>vmax) { bbb = vmax; }
220 if(kkk < 1) { kkk = 1; }
225 for(
G4int l=0; l<kkk; l++)
227 for(
G4int i=0; i<6; i++)
235 loss *=hhh*totalEnergy ;
249 static const G4int k2 = 4;
252 if(cut >= tkin)
return cross;
260 if(kkk < 1) { kkk = 1; }
266 for(
G4int l=0; l<kkk; l++)
268 for(
G4int i=0; i<6; i++)
293 if(gammaEnergy > tkin) {
return dxsection; }
297 G4double delta = 0.5*mass*mass*v/(E-gammaEnergy) ;
301 if(iz < 1) { iz = 1; }
302 else if(iz > 92) { iz = 92; }
319 G4double fn=
G4Log(rab1/(dnstar*(electron_mass_c2+rab0*rab1))*
320 (mass+delta*(dnstar*sqrte-2.))) ;
321 if(fn <0.) { fn = 0.; }
325 if(gammaEnergy<epmax1)
328 fe=
G4Log(rab2*mass/((1.+delta*
rmass/(electron_mass_c2*sqrte))*
329 (electron_mass_c2+rab0*rab2))) ;
333 dxsection =
coeff*(1.-v*(1. - 0.75*v))*Z*(fn*Z +
fe)/gammaEnergy;
352 if (cut >= tmax)
return cross;
355 if(tmax < kineticEnergy) {
364 std::vector<G4DynamicParticle*>* vdp,
375 if(tmin >= tmax)
return;
386 G4double totalMomentum = sqrt(kineticEnergy*(kineticEnergy + 2.0*mass));
417 G4double dirx = sint*cos(phi), diry = sint*sin(phi), dirz = cos(theta) ;
420 gDirection.rotateUz(partDirection);
422 partDirection *= totalMomentum;
423 partDirection -= gEnergy*gDirection;
424 partDirection = partDirection.unit();
427 kineticEnergy -= gEnergy;
434 vdp->push_back(aGamma);
G4MuBremsstrahlungModel(const G4ParticleDefinition *p=0, const G4String &nam="MuBrem")
G4ParticleChangeForLoss * GetParticleChangeForLoss()
std::vector< G4Element * > G4ElementVector
virtual void InitialiseLocal(const G4ParticleDefinition *, G4VEmModel *masterModel)
G4ParticleDefinition * theGamma
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
void InitialiseElementSelectors(const G4ParticleDefinition *, const G4DataVector &)
G4double HighEnergyLimit() const
G4double ComputeMicroscopicCrossSection(G4double tkin, G4double Z, G4double cut)
G4double GetZ13(G4double Z)
void SetParticle(const G4ParticleDefinition *)
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
const G4ElementVector * GetElementVector() const
static G4NistManager * Instance()
virtual G4double ComputeDMicroscopicCrossSection(G4double tkin, G4double Z, G4double gammaEnergy)
static const G4double xgi[6]
virtual G4double MinEnergyCut(const G4ParticleDefinition *, const G4MaterialCutsCouple *)
static const G4double wgi[6]
static const G4double sqrte
const G4ParticleDefinition * particle
const G4ThreeVector & GetMomentumDirection() const
static const double twopi
void SetProposedKineticEnergy(G4double proposedKinEnergy)
virtual ~G4MuBremsstrahlungModel()
std::vector< G4EmElementSelector * > * GetElementSelectors()
void SetProposedMomentumDirection(const G4ThreeVector &dir)
const G4double * GetAtomicNumDensityVector() const
G4double G4Log(G4double x)
virtual G4double MinPrimaryEnergy(const G4Material *, const G4ParticleDefinition *, G4double)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy)
G4double ComputMuBremLoss(G4double Z, G4double tkin, G4double cut)
void SetElementSelectors(std::vector< G4EmElementSelector * > *)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4ParticleChangeForLoss * fParticleChange
const G4double x[NPOINTSGL]
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kineticEnergy, G4double Z, G4double A, G4double cutEnergy, G4double maxEnergy)
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
size_t GetNumberOfElements() const
G4ThreeVector G4ParticleMomentum
const G4Element * SelectRandomAtom(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)