66 HighestEnergyLimit = 1000*
TeV;
101 CrossSecFactor = fac;
102 G4cout <<
"The cross section for AnnihiToMuPair is artificially "
103 <<
"increased by the CrossSecFactor=" << CrossSecFactor <<
G4endl;
114 static const G4double Sig0 =
pi*Rmuon*Rmuon/3.;
117 if (Epos < LowestEnergyLimit)
return CrossSection;
119 G4double xi = LowestEnergyLimit/Epos;
120 G4double SigmaEl = Sig0*xi*(1.+xi/2.)*sqrt(1.-xi);
121 CrossSection = SigmaEl*
Z;
137 G4double AtomicZ = (*theElementVector)[i]->GetZ();
138 SIGMA += NbOfAtomsPerVolume[i] *
160 if(CurrentSigma >
DBL_MIN) mfp = 1.0/(CurrentSigma*CrossSecFactor);
189 if (Epos < LowestEnergyLimit) {
195 G4double xi = LowestEnergyLimit/Epos;
210 G4double Ecm = sqrt(0.5*Mele*(Epos+Mele));
211 G4double Pcm = sqrt(Ecm*Ecm-Mmuon*Mmuon);
212 G4double beta = sqrt((Epos-Mele)/(Epos+Mele));
218 G4double EmuPlus = gamma*( Ecm+cost*beta*Pcm);
219 G4double EmuMinus = gamma*( Ecm-cost*beta*Pcm);
220 G4double PmuPlusZ = gamma*(beta*Ecm+cost* Pcm);
221 G4double PmuMinusZ = gamma*(beta*Ecm-cost* Pcm);
227 G4double PmuPlus = sqrt(Pt*Pt+PmuPlusZ *PmuPlusZ );
228 G4double PmuMinus = sqrt(Pt*Pt+PmuMinusZ*PmuMinusZ);
233 MuPlusDirection ( PmuPlusX/PmuPlus, PmuPlusY/PmuPlus, PmuPlusZ/PmuPlus );
235 MuMinusDirection(PmuMinusX/PmuMinus,PmuMinusY/PmuMinus,PmuMinusZ/PmuMinus);
239 MuPlusDirection.
rotateUz(PositronDirection);
240 MuMinusDirection.
rotateUz(PositronDirection);
264 G4String comments =
"e+e->mu+mu- annihilation, atomic e- at rest, SubType=.";
267 G4cout <<
" threshold at " << LowestEnergyLimit/
GeV <<
" GeV"
268 <<
" good description up to "
269 << HighestEnergyLimit/
TeV <<
" TeV for all Z." <<
G4endl;
static G4MuonPlus * MuonPlus()
std::vector< G4Element * > G4ElementVector
void PrintInfoDefinition()
G4double GetKineticEnergy() const
const G4DynamicParticle * GetDynamicParticle() const
G4double ComputeCrossSectionPerAtom(G4double PositronEnergy, G4double AtomicZ)
G4bool IsApplicable(const G4ParticleDefinition &)
G4AnnihiToMuPair(const G4String &processName="AnnihiToMuPair", G4ProcessType type=fElectromagnetic)
const G4ElementVector * GetElementVector() const
const G4double * GetVecNbOfAtomsPerVolume() const
G4GLOB_DLL std::ostream G4cout
G4VParticleChange * PostStepDoIt(const G4Track &aTrack, const G4Step &aStep)
const G4ThreeVector & GetMomentumDirection() const
Hep3Vector & rotateUz(const Hep3Vector &)
void SetProcessSubType(G4int)
const G4String & GetProcessName() const
G4Material * GetMaterial() const
virtual void Initialize(const G4Track &)
static G4Positron * Positron()
void BuildPhysicsTable(const G4ParticleDefinition &)
G4double GetMeanFreePath(const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *)
G4double GetPDGMass() const
void SetNumberOfSecondaries(G4int totSecondaries)
void ProposeEnergy(G4double finalEnergy)
G4ParticleChange aParticleChange
void AddSecondary(G4Track *aSecondary)
static G4MuonMinus * MuonMinus()
size_t GetNumberOfElements() const
G4double CrossSectionPerVolume(G4double PositronEnergy, const G4Material *)
void SetCrossSecFactor(G4double fac)
void ProposeTrackStatus(G4TrackStatus status)
G4int GetProcessSubType() const
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)