Geant4  10.03.p01
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Member Functions | List of all members
G4GammaConversionToMuons Class Reference

#include <G4GammaConversionToMuons.hh>

Inheritance diagram for G4GammaConversionToMuons:
Inheritance graph
[legend]
Collaboration diagram for G4GammaConversionToMuons:
Collaboration graph
[legend]

Public Member Functions

 G4GammaConversionToMuons (const G4String &processName="GammaToMuPair", G4ProcessType type=fElectromagnetic)
 
 ~G4GammaConversionToMuons ()
 
G4bool IsApplicable (const G4ParticleDefinition &) override
 
void BuildPhysicsTable (const G4ParticleDefinition &) override
 
void PrintInfoDefinition ()
 
void SetCrossSecFactor (G4double fac)
 
G4double GetCrossSecFactor ()
 
G4double GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition) override
 
G4double GetCrossSectionPerAtom (const G4DynamicParticle *aDynamicGamma, G4Element *anElement)
 
G4VParticleChange * PostStepDoIt (const G4Track &aTrack, const G4Step &aStep) override
 
virtual G4double ComputeCrossSectionPerAtom (G4double GammaEnergy, G4double AtomicZ, G4double AtomicA)
 
G4double ComputeMeanFreePath (G4double GammaEnergy, G4Material *aMaterial)
 
- Public Member Functions inherited from G4VDiscreteProcess
 G4VDiscreteProcess (const G4String &, G4ProcessType aType=fNotDefined)
 
 G4VDiscreteProcess (G4VDiscreteProcess &)
 
virtual ~G4VDiscreteProcess ()
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
 
virtual G4double AlongStepGetPhysicalInteractionLength (const G4Track &, G4double, G4double, G4double &, G4GPILSelection *)
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
 
virtual G4VParticleChange * AtRestDoIt (const G4Track &, const G4Step &)
 
virtual G4VParticleChange * AlongStepDoIt (const G4Track &, const G4Step &)
 
- Public Member Functions inherited from G4VProcess
 G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)
 
 G4VProcess (const G4VProcess &right)
 
virtual ~G4VProcess ()
 
G4int operator== (const G4VProcess &right) const
 
G4int operator!= (const G4VProcess &right) const
 
G4double GetCurrentInteractionLength () const
 
void SetPILfactor (G4double value)
 
G4double GetPILfactor () const
 
G4double AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)
 
G4double AtRestGPIL (const G4Track &track, G4ForceCondition *condition)
 
G4double PostStepGPIL (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
 
virtual void PreparePhysicsTable (const G4ParticleDefinition &)
 
virtual G4bool StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
 
virtual G4bool RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
 
const G4String & GetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)
 
const G4String & GetProcessName () const
 
G4ProcessType GetProcessType () const
 
void SetProcessType (G4ProcessType)
 
G4int GetProcessSubType () const
 
void SetProcessSubType (G4int)
 
virtual void StartTracking (G4Track *)
 
virtual void EndTracking ()
 
virtual void SetProcessManager (const G4ProcessManager *)
 
virtual const G4ProcessManager * GetProcessManager ()
 
virtual void ResetNumberOfInteractionLengthLeft ()
 
G4double GetNumberOfInteractionLengthLeft () const
 
G4double GetTotalNumberOfInteractionLengthTraversed () const
 
G4bool isAtRestDoItIsEnabled () const
 
G4bool isAlongStepDoItIsEnabled () const
 
G4bool isPostStepDoItIsEnabled () const
 
virtual void DumpInfo () const
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
virtual void SetMasterProcess (G4VProcess *masterP)
 
const G4VProcess * GetMasterProcess () const
 
virtual void BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
 
virtual void PrepareWorkerPhysicsTable (const G4ParticleDefinition &)
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VProcess
static const G4String & GetProcessTypeName (G4ProcessType)
 
- Protected Member Functions inherited from G4VProcess
void SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)
 
void ClearNumberOfInteractionLengthLeft ()
 
- Protected Attributes inherited from G4VProcess
const G4ProcessManager * aProcessManager
 
G4VParticleChange * pParticleChange
 
G4ParticleChange aParticleChange
 
G4double theNumberOfInteractionLengthLeft
 
G4double currentInteractionLength
 
G4double theInitialNumberOfInteractionLength
 
G4String theProcessName
 
G4String thePhysicsTableFileName
 
G4ProcessType theProcessType
 
G4int theProcessSubType
 
G4double thePILfactor
 
G4bool enableAtRestDoIt
 
G4bool enableAlongStepDoIt
 
G4bool enablePostStepDoIt
 
G4int verboseLevel
 

Detailed Description

Definition at line 61 of file G4GammaConversionToMuons.hh.

Constructor & Destructor Documentation

G4GammaConversionToMuons::G4GammaConversionToMuons ( const G4String &  processName = "GammaToMuPair",
G4ProcessType  type = fElectromagnetic 
)
explicit

Definition at line 55 of file G4GammaConversionToMuons.cc.

57  : G4VDiscreteProcess (processName, type),
58  Mmuon(G4MuonPlus::MuonPlus()->GetPDGMass()),
59  Rc(elm_coupling/Mmuon),
60  LowestEnergyLimit (4.*Mmuon), // 4*Mmuon
61  HighestEnergyLimit(1e21*eV), // ok to 1e21eV=1e12GeV, then LPM suppression
62  CrossSecFactor(1.)
63 {
64  SetProcessSubType(fGammaConversionToMuMu);
65  MeanFreePath = DBL_MAX;
66 }
G4MuonPlus::MuonPlus
static G4MuonPlus * MuonPlus()
Definition: G4MuonPlus.cc:99
fGammaConversionToMuMu
Definition: G4EmProcessSubType.hh:61
eV
static constexpr double eV
Definition: G4SIunits.hh:215
G4VProcess::SetProcessSubType
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:432
CLHEP::elm_coupling
static constexpr double elm_coupling
Definition: PhysicalConstants.h:95
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83

Here is the call graph for this function:

G4GammaConversionToMuons::~G4GammaConversionToMuons ( )

Definition at line 72 of file G4GammaConversionToMuons.cc.

73 {}

Member Function Documentation

void G4GammaConversionToMuons::BuildPhysicsTable ( const G4ParticleDefinition &  )
overridevirtual

Reimplemented from G4VProcess.

Definition at line 85 of file G4GammaConversionToMuons.cc.

87 { //here no tables, just calling PrintInfoDefinition
88  PrintInfoDefinition();
89 }

Here is the call graph for this function:

G4double G4GammaConversionToMuons::ComputeCrossSectionPerAtom ( G4double  GammaEnergy,
G4double  AtomicZ,
G4double  AtomicA 
)
virtual

Definition at line 152 of file G4GammaConversionToMuons.cc.

158 {
159  if(Egam <= LowestEnergyLimit) return 0.0 ; // below threshold return 0
160 
161  G4int Z = G4lrint(ZZ);
162  G4double CrossSection = 0.0;
163  G4NistManager* nist = G4NistManager::Instance();
164 
165  G4double PowThres,Ecor,B,Dn,Zthird,Winfty,WMedAppr,
166  Wsatur,sigfac;
167 
168  if(Z==1) // special case of Hydrogen
169  { B=202.4;
170  Dn=1.49;
171  }
172  else
173  { B=183.;
174  Dn=1.54*nist->GetA27(Z);
175  }
176  Zthird=1./nist->GetZ13(Z); // Z**(-1/3)
177  Winfty=B*Zthird*Mmuon/(Dn*electron_mass_c2);
178  WMedAppr=1./(4.*Dn*sqrte*Mmuon);
179  Wsatur=Winfty/WMedAppr;
180  sigfac=4.*fine_structure_const*Z*Z*Rc*Rc;
181  PowThres=1.479+0.00799*Dn;
182  Ecor=-18.+4347./(B*Zthird);
183 
184  G4double CorFuc=1.+.04*G4Log(1.+Ecor/Egam);
185  //G4double Eg=pow(1.-4.*Mmuon/Egam,PowThres)*pow( pow(Wsatur,PowSat)+
186  // pow(Egam,PowSat),1./PowSat); // threshold and saturation
187  G4double Eg=G4Exp(G4Log(1.-4.*Mmuon/Egam)*PowThres)*
188  G4Exp(G4Log( G4Exp(G4Log(Wsatur)*PowSat)+G4Exp(G4Log(Egam)*PowSat))/PowSat);
189  CrossSection=7./9.*sigfac*G4Log(1.+WMedAppr*CorFuc*Eg);
190  CrossSection*=CrossSecFactor; // increase the CrossSection by (by default 1)
191  return CrossSection;
192 }
G4NistManager::GetA27
G4double GetA27(G4int Z) const
Definition: G4NistManager.hh:586
B
double B(double temperature)
Definition: G4DNAElectronHoleRecombination.cc:69
PowSat
static const G4double PowSat
Definition: G4GammaConversionToMuons.cc:53
G4int
int G4int
Definition: G4Types.hh:78
G4NistManager::Instance
static G4NistManager * Instance()
Definition: G4NistManager.cc:68
CLHEP::electron_mass_c2
static constexpr double electron_mass_c2
Definition: PhysicalConstants.h:79
G4NistManager::GetZ13
G4double GetZ13(G4double Z) const
Definition: G4NistManager.hh:572
sqrte
static const G4double sqrte
Definition: G4GammaConversionToMuons.cc:52
G4Log
G4double G4Log(G4double x)
Definition: G4Log.hh:230
G4Exp
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
Definition: G4Exp.hh:183
G4lrint
int G4lrint(double ad)
Definition: templates.hh:163
G4double
double G4double
Definition: G4Types.hh:76
CLHEP::fine_structure_const
static constexpr double fine_structure_const
Definition: PhysicalConstants.h:96
Z
G4int Z
Definition: G4ParticleHPJENDLHEData.cc:262

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4GammaConversionToMuons::ComputeMeanFreePath ( G4double  GammaEnergy,
G4Material *  aMaterial 
)

Definition at line 114 of file G4GammaConversionToMuons.cc.

118 {
119  const G4ElementVector* theElementVector = aMaterial->GetElementVector();
120  const G4double* NbOfAtomsPerVolume = aMaterial->GetVecNbOfAtomsPerVolume();
121 
122  G4double SIGMA = 0 ;
123 
124  for ( size_t i=0 ; i < aMaterial->GetNumberOfElements(); ++i)
125  {
126  G4double AtomicZ = (*theElementVector)[i]->GetZ();
127  G4double AtomicA = (*theElementVector)[i]->GetA()/(g/mole);
128  SIGMA += NbOfAtomsPerVolume[i] *
129  ComputeCrossSectionPerAtom(GammaEnergy,AtomicZ,AtomicA);
130  }
131  return SIGMA > DBL_MIN ? 1./SIGMA : DBL_MAX;
132 }
G4ElementVector
std::vector< G4Element * > G4ElementVector
Definition: G4ElementVector.hh:44
g
static constexpr double g
Definition: G4SIunits.hh:183
G4Material::GetElementVector
const G4ElementVector * GetElementVector() const
Definition: G4Material.hh:190
G4Material::GetVecNbOfAtomsPerVolume
const G4double * GetVecNbOfAtomsPerVolume() const
Definition: G4Material.hh:206
G4GammaConversionToMuons::ComputeCrossSectionPerAtom
virtual G4double ComputeCrossSectionPerAtom(G4double GammaEnergy, G4double AtomicZ, G4double AtomicA)
Definition: G4GammaConversionToMuons.cc:152
DBL_MIN
#define DBL_MIN
Definition: templates.hh:75
G4Material::GetNumberOfElements
size_t GetNumberOfElements() const
Definition: G4Material.hh:186
G4double
double G4double
Definition: G4Types.hh:76
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
mole
static constexpr double mole
Definition: G4SIunits.hh:286

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4GammaConversionToMuons::GetCrossSecFactor ( )
inline

Definition at line 87 of file G4GammaConversionToMuons.hh.

87 { return CrossSecFactor;}
G4double G4GammaConversionToMuons::GetCrossSectionPerAtom ( const G4DynamicParticle *  aDynamicGamma,
G4Element *  anElement 
)

Definition at line 136 of file G4GammaConversionToMuons.cc.

141 {
142  G4double GammaEnergy = aDynamicGamma->GetKineticEnergy();
143  G4double AtomicZ = anElement->GetZ();
144  G4double AtomicA = anElement->GetN();
145  G4double crossSection =
146  ComputeCrossSectionPerAtom(GammaEnergy,AtomicZ,AtomicA);
147  return crossSection;
148 }
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4Element::GetN
G4double GetN() const
Definition: G4Element.hh:135
G4Element::GetZ
G4double GetZ() const
Definition: G4Element.hh:131
G4GammaConversionToMuons::ComputeCrossSectionPerAtom
virtual G4double ComputeCrossSectionPerAtom(G4double GammaEnergy, G4double AtomicZ, G4double AtomicA)
Definition: G4GammaConversionToMuons.cc:152
G4double
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

G4double G4GammaConversionToMuons::GetMeanFreePath ( const G4Track &  aTrack,
G4double  previousStepSize,
G4ForceCondition *  condition 
)
overridevirtual

Implements G4VDiscreteProcess.

Definition at line 93 of file G4GammaConversionToMuons.cc.

99 {
100  const G4DynamicParticle* aDynamicGamma = aTrack.GetDynamicParticle();
101  G4double GammaEnergy = aDynamicGamma->GetKineticEnergy();
102  G4Material* aMaterial = aTrack.GetMaterial();
103 
104  if (GammaEnergy <= LowestEnergyLimit)
105  MeanFreePath = DBL_MAX;
106  else
107  MeanFreePath = ComputeMeanFreePath(GammaEnergy,aMaterial);
108 
109  return MeanFreePath;
110 }
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4Track::GetDynamicParticle
const G4DynamicParticle * GetDynamicParticle() const
G4GammaConversionToMuons::ComputeMeanFreePath
G4double ComputeMeanFreePath(G4double GammaEnergy, G4Material *aMaterial)
Definition: G4GammaConversionToMuons.cc:114
G4Track::GetMaterial
G4Material * GetMaterial() const
G4double
double G4double
Definition: G4Types.hh:76
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83

Here is the call graph for this function:

G4bool G4GammaConversionToMuons::IsApplicable ( const G4ParticleDefinition &  particle)
overridevirtual

Reimplemented from G4VProcess.

Definition at line 77 of file G4GammaConversionToMuons.cc.

79 {
80  return ( &particle == G4Gamma::Gamma() );
81 }
G4Gamma::Gamma
static G4Gamma * Gamma()
Definition: G4Gamma.cc:86

Here is the call graph for this function:

G4VParticleChange * G4GammaConversionToMuons::PostStepDoIt ( const G4Track &  aTrack,
const G4Step &  aStep 
)
overridevirtual

Reimplemented from G4VDiscreteProcess.

Definition at line 206 of file G4GammaConversionToMuons.cc.

212 {
213  aParticleChange.Initialize(aTrack);
214  G4Material* aMaterial = aTrack.GetMaterial();
215 
216  // current Gamma energy and direction, return if energy too low
217  const G4DynamicParticle *aDynamicGamma = aTrack.GetDynamicParticle();
218  G4double Egam = aDynamicGamma->GetKineticEnergy();
219  if (Egam <= LowestEnergyLimit) {
220  return G4VDiscreteProcess::PostStepDoIt(aTrack,aStep);
221  }
222  G4ParticleMomentum GammaDirection = aDynamicGamma->GetMomentumDirection();
223 
224  // select randomly one element constituting the material
225  const G4Element* anElement = SelectRandomAtom(aDynamicGamma, aMaterial);
226  G4int Z = G4lrint(anElement->GetZ());
227  G4NistManager* nist = G4NistManager::Instance();
228 
229  G4double B,Dn;
230  G4double A027 = nist->GetA27(Z);
231 
232  if(Z==1) // special case of Hydrogen
233  { B=202.4;
234  Dn=1.49;
235  }
236  else
237  { B=183.;
238  Dn=1.54*A027;
239  }
240  G4double Zthird=1./nist->GetZ13(Z); // Z**(-1/3)
241  G4double Winfty=B*Zthird*Mmuon/(Dn*electron_mass_c2);
242  G4double C1Num=0.35*A027;
243  G4double C1Num2=C1Num*C1Num;
244  G4double C2Term2=electron_mass_c2/(183.*Zthird*Mmuon);
245 
246  G4double GammaMuonInv=Mmuon/Egam;
247  G4double sqrtx=sqrt(.25-GammaMuonInv);
248  G4double xmax=.5+sqrtx;
249  G4double xmin=.5-sqrtx;
250 
251  // generate xPlus according to the differential cross section by rejection
252  G4double Ds2=(Dn*sqrte-2.);
253  G4double sBZ=sqrte*B*Zthird/electron_mass_c2;
254  G4double LogWmaxInv=1./G4Log(Winfty*(1.+2.*Ds2*GammaMuonInv)
255  /(1.+2.*sBZ*Mmuon*GammaMuonInv));
256  G4double xPlus,xMinus,xPM,result,W;
257  G4int nn = 0;
258  const G4int nmax = 1000;
259  do
260  { xPlus=xmin+G4UniformRand()*(xmax-xmin);
261  xMinus=1.-xPlus;
262  xPM=xPlus*xMinus;
263  G4double del=Mmuon*Mmuon/(2.*Egam*xPM);
264  W=Winfty*(1.+Ds2*del/Mmuon)/(1.+sBZ*del);
265  if(W<=1. || nn > nmax) { break; } // to avoid negative cross section at xmin
266  G4double xxp=1.-4./3.*xPM; // the main xPlus dependence
267  result=xxp*G4Log(W)*LogWmaxInv;
268  if(result>1.) {
269  G4cout << "G4GammaConversionToMuons::PostStepDoIt WARNING:"
270  << " in dSigxPlusGen, result=" << result << " > 1" << G4endl;
271  }
272  ++nn;
273  if(nn >= nmax) { break; }
274  }
275  // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
276  while (G4UniformRand() > result);
277 
278  // now generate the angular variables via the auxilary variables t,psi,rho
279  G4double t;
280  G4double psi;
281  G4double rho;
282 
283  G4double thetaPlus,thetaMinus,phiHalf; // final angular variables
284  nn = 0;
285  do // t, psi, rho generation start (while angle < pi)
286  {
287  //generate t by the rejection method
288  G4double C1=C1Num2* GammaMuonInv/xPM;
289  G4double f1_max=(1.-xPM) / (1.+C1);
290  G4double f1; // the probability density
291  do
292  {
293  ++nn;
294  t=G4UniformRand();
295  f1=(1.-2.*xPM+4.*xPM*t*(1.-t)) / (1.+C1/(t*t));
296  if(f1<0 || f1> f1_max) // should never happend
297  {
298  G4cout << "G4GammaConversionToMuons::PostStepDoIt WARNING:"
299  << "outside allowed range f1=" << f1 << " is set to zero"
300  << G4endl;
301  f1 = 0.0;
302  }
303  if(nn > nmax) { break; }
304  }
305  // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
306  while ( G4UniformRand()*f1_max > f1);
307  // generate psi by the rejection method
308  G4double f2_max=1.-2.*xPM*(1.-4.*t*(1.-t));
309 
310  // long version
311  G4double f2;
312  do
313  {
314  ++nn;
315  psi=2.*pi*G4UniformRand();
316  f2=1.-2.*xPM+4.*xPM*t*(1.-t)*(1.+cos(2.*psi));
317  if(f2<0 || f2> f2_max) // should never happend
318  {
319  G4cout << "G4GammaConversionToMuons::PostStepDoIt WARNING:"
320  << "outside allowed range f2=" << f2 << " is set to zero"
321  << G4endl;
322  f2 = 0.0;
323  }
324  if(nn >= nmax) { break; }
325  }
326  // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
327  while ( G4UniformRand()*f2_max > f2);
328 
329  // generate rho by direct transformation
330  G4double C2Term1=GammaMuonInv/(2.*xPM*t);
331  G4double C2=4./sqrt(xPM)*pow(C2Term1*C2Term1+C2Term2*C2Term2,2.);
332  G4double rhomax=1.9/A027*(1./t-1.);
333  G4double beta=G4Log( (C2+rhomax*rhomax*rhomax*rhomax)/C2 );
334  rho=G4Exp(G4Log(C2 *( G4Exp(beta*G4UniformRand())-1. ))*0.25);
335 
336  //now get from t and psi the kinematical variables
337  G4double u=sqrt(1./t-1.);
338  G4double xiHalf=0.5*rho*cos(psi);
339  phiHalf=0.5*rho/u*sin(psi);
340 
341  thetaPlus =GammaMuonInv*(u+xiHalf)/xPlus;
342  thetaMinus=GammaMuonInv*(u-xiHalf)/xMinus;
343 
344  // protection against infinite loop
345  if(nn > nmax) {
346  if(std::abs(thetaPlus)>pi) { thetaPlus = 0.0; }
347  if(std::abs(thetaMinus)>pi) { thetaMinus = 0.0; }
348  }
349 
350  // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
351  } while ( std::abs(thetaPlus)>pi || std::abs(thetaMinus) >pi);
352 
353  // now construct the vectors
354  // azimuthal symmetry, take phi0 at random between 0 and 2 pi
355  G4double phi0=2.*pi*G4UniformRand();
356  G4double EPlus=xPlus*Egam;
357  G4double EMinus=xMinus*Egam;
358 
359  // mu+ mu- directions for gamma in z-direction
360  G4ThreeVector MuPlusDirection ( sin(thetaPlus) *cos(phi0+phiHalf),
361  sin(thetaPlus) *sin(phi0+phiHalf), cos(thetaPlus) );
362  G4ThreeVector MuMinusDirection (-sin(thetaMinus)*cos(phi0-phiHalf),
363  -sin(thetaMinus) *sin(phi0-phiHalf), cos(thetaMinus) );
364  // rotate to actual gamma direction
365  MuPlusDirection.rotateUz(GammaDirection);
366  MuMinusDirection.rotateUz(GammaDirection);
367  aParticleChange.SetNumberOfSecondaries(2);
368  // create G4DynamicParticle object for the particle1
369  G4DynamicParticle* aParticle1= new G4DynamicParticle(
370  G4MuonPlus::MuonPlus(),MuPlusDirection,EPlus-Mmuon);
371  aParticleChange.AddSecondary(aParticle1);
372  // create G4DynamicParticle object for the particle2
373  G4DynamicParticle* aParticle2= new G4DynamicParticle(
374  G4MuonMinus::MuonMinus(),MuMinusDirection,EMinus-Mmuon);
375  aParticleChange.AddSecondary(aParticle2);
376  //
377  // Kill the incident photon
378  //
379  aParticleChange.ProposeMomentumDirection( 0., 0., 0. ) ;
380  aParticleChange.ProposeEnergy( 0. ) ;
381  aParticleChange.ProposeTrackStatus( fStopAndKill ) ;
382  // Reset NbOfInteractionLengthLeft and return aParticleChange
383  return G4VDiscreteProcess::PostStepDoIt( aTrack, aStep );
384 }
result
G4double G4ParticleHPJENDLHEData::G4double result
Definition: G4ParticleHPJENDLHEData.cc:257
C2
const double C2
Definition: MCGIDI_KalbachMann.cc:15
G4MuonPlus::MuonPlus
static G4MuonPlus * MuonPlus()
Definition: G4MuonPlus.cc:99
C1
const double C1
Definition: MCGIDI_KalbachMann.cc:15
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4Track::GetDynamicParticle
const G4DynamicParticle * GetDynamicParticle() const
G4Element::GetZ
G4double GetZ() const
Definition: G4Element.hh:131
B
double B(double temperature)
Definition: G4DNAElectronHoleRecombination.cc:69
G4int
int G4int
Definition: G4Types.hh:78
G4NistManager::Instance
static G4NistManager * Instance()
Definition: G4NistManager.cc:68
CLHEP::electron_mass_c2
static constexpr double electron_mass_c2
Definition: PhysicalConstants.h:79
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:97
G4cout
G4GLOB_DLL std::ostream G4cout
sqrte
static const G4double sqrte
Definition: G4GammaConversionToMuons.cc:52
G4DynamicParticle::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4Track::GetMaterial
G4Material * GetMaterial() const
G4Log
G4double G4Log(G4double x)
Definition: G4Log.hh:230
G4Exp
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
Definition: G4Exp.hh:183
G4ParticleChange::Initialize
virtual void Initialize(const G4Track &)
Definition: G4ParticleChange.cc:228
G4lrint
int G4lrint(double ad)
Definition: templates.hh:163
G4VParticleChange::SetNumberOfSecondaries
void SetNumberOfSecondaries(G4int totSecondaries)
G4ParticleChange::ProposeEnergy
void ProposeEnergy(G4double finalEnergy)
G4VProcess::aParticleChange
G4ParticleChange aParticleChange
Definition: G4VProcess.hh:289
fStopAndKill
Definition: G4TrackStatus.hh:56
G4ParticleChange::AddSecondary
void AddSecondary(G4Track *aSecondary)
Definition: G4ParticleChange.cc:218
G4MuonMinus::MuonMinus
static G4MuonMinus * MuonMinus()
Definition: G4MuonMinus.cc:100
G4endl
#define G4endl
Definition: G4ios.hh:61
pi
static constexpr double pi
Definition: G4SIunits.hh:75
G4ParticleChange::ProposeMomentumDirection
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
G4double
double G4double
Definition: G4Types.hh:76
G4VParticleChange::ProposeTrackStatus
void ProposeTrackStatus(G4TrackStatus status)
G4VDiscreteProcess::PostStepDoIt
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
Definition: G4VDiscreteProcess.cc:112
Z
G4int Z
Definition: G4ParticleHPJENDLHEData.cc:262

Here is the call graph for this function:

void G4GammaConversionToMuons::PrintInfoDefinition ( )

Definition at line 416 of file G4GammaConversionToMuons.cc.

417 {
418  G4String comments ="gamma->mu+mu- Bethe Heitler process, SubType= ";
419  G4cout << G4endl << GetProcessName() << ": " << comments
420  << GetProcessSubType() << G4endl;
421  G4cout << " good cross section parametrization from "
422  << G4BestUnit(LowestEnergyLimit,"Energy")
423  << " to " << HighestEnergyLimit/GeV << " GeV for all Z." << G4endl;
424 }
G4BestUnit
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
Definition: G4SteppingVerbose.cc:53
G4cout
G4GLOB_DLL std::ostream G4cout
G4VProcess::GetProcessName
const G4String & GetProcessName() const
Definition: G4VProcess.hh:408
GeV
static constexpr double GeV
Definition: G4SIunits.hh:217
G4endl
#define G4endl
Definition: G4ios.hh:61
G4VProcess::GetProcessSubType
G4int GetProcessSubType() const
Definition: G4VProcess.hh:426

Here is the call graph for this function:

Here is the caller graph for this function:

void G4GammaConversionToMuons::SetCrossSecFactor ( G4double  fac)

Definition at line 196 of file G4GammaConversionToMuons.cc.

198 {
199  CrossSecFactor=fac;
200  G4cout << "The cross section for GammaConversionToMuons is artificially "
201  << "increased by the CrossSecFactor=" << CrossSecFactor << G4endl;
202 }
G4cout
G4GLOB_DLL std::ostream G4cout
fac
static const G4double fac
Definition: G4NISTStoppingData.hh:87
G4endl
#define G4endl
Definition: G4ios.hh:61
The documentation for this class was generated from the following files: