Geant4  10.02.p03
Public Member Functions | Private Member Functions | Private Attributes | List of all members
G4EmCaptureCascade Class Reference

#include <G4EmCaptureCascade.hh>

Inheritance diagram for G4EmCaptureCascade:
Collaboration diagram for G4EmCaptureCascade:

Public Member Functions

 G4EmCaptureCascade ()
 
 ~G4EmCaptureCascade ()
 
G4HadFinalState * ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
 
void ModelDescription (std::ostream &outFile) const
 
- Public Member Functions inherited from G4HadronicInteraction
 G4HadronicInteraction (const G4String &modelName="HadronicModel")
 
virtual ~G4HadronicInteraction ()
 
virtual G4double SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
 
virtual G4bool IsApplicable (const G4HadProjectile &, G4Nucleus &)
 
G4double GetMinEnergy () const
 
G4double GetMinEnergy (const G4Material *aMaterial, const G4Element *anElement) const
 
void SetMinEnergy (G4double anEnergy)
 
void SetMinEnergy (G4double anEnergy, const G4Element *anElement)
 
void SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)
 
G4double GetMaxEnergy () const
 
G4double GetMaxEnergy (const G4Material *aMaterial, const G4Element *anElement) const
 
void SetMaxEnergy (const G4double anEnergy)
 
void SetMaxEnergy (G4double anEnergy, const G4Element *anElement)
 
void SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)
 
const G4HadronicInteraction * GetMyPointer () const
 
virtual G4int GetVerboseLevel () const
 
virtual void SetVerboseLevel (G4int value)
 
const G4String & GetModelName () const
 
void DeActivateFor (const G4Material *aMaterial)
 
void ActivateFor (const G4Material *aMaterial)
 
void DeActivateFor (const G4Element *anElement)
 
void ActivateFor (const G4Element *anElement)
 
G4bool IsBlocked (const G4Material *aMaterial) const
 
G4bool IsBlocked (const G4Element *anElement) const
 
void SetRecoilEnergyThreshold (G4double val)
 
G4double GetRecoilEnergyThreshold () const
 
G4bool operator== (const G4HadronicInteraction &right) const
 
G4bool operator!= (const G4HadronicInteraction &right) const
 
virtual const std::pair< G4double, G4double > GetFatalEnergyCheckLevels () const
 
virtual std::pair< G4double, G4double > GetEnergyMomentumCheckLevels () const
 
void SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 

Private Member Functions

void AddNewParticle (G4ParticleDefinition *aParticle, G4double kinEnergy)
 
G4EmCaptureCascade & operator= (const G4EmCaptureCascade &right)
 
 G4EmCaptureCascade (const G4EmCaptureCascade &)
 

Private Attributes

G4HadFinalState result
 
G4ParticleDefinition * theElectron
 
G4ParticleDefinition * theGamma
 
G4double fMuMass
 
G4double fTime
 
G4double fLevelEnergy [14]
 
G4double fKLevelEnergy [93]
 

Additional Inherited Members

- Protected Member Functions inherited from G4HadronicInteraction
void SetModelName (const G4String &nam)
 
G4bool IsBlocked () const
 
void Block ()
 
- Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState theParticleChange
 
G4int verboseLevel
 
G4double theMinEnergy
 
G4double theMaxEnergy
 
G4bool isBlocked
 

Detailed Description

Definition at line 67 of file G4EmCaptureCascade.hh.

Constructor & Destructor Documentation

◆ G4EmCaptureCascade() [1/2]

G4EmCaptureCascade::G4EmCaptureCascade ( )

Definition at line 56 of file G4EmCaptureCascade.cc.

57  : G4HadronicInteraction("emCaptureCascade")
58 {
59  theElectron = G4Electron::Electron();
60  theGamma = G4Gamma::Gamma();
61  fMuMass = G4MuonMinus::MuonMinus()->GetPDGMass();
62  fTime = 0.0;
63 
64  // Calculate the Energy of K Mesoatom Level for this Element using
65  // the Energy of Hydrogen Atom taken into account finite size of the
66  // nucleus
67  const G4int nlevels = 28;
68  const G4int listK[nlevels] = {
69  1, 2, 4, 6, 8, 11, 14, 17, 18, 21, 24,
70  26, 29, 32, 38, 40, 41, 44, 49, 53, 55,
71  60, 65, 70, 75, 81, 85, 92};
72  const G4double listKEnergy[nlevels] = {
73  0.00275, 0.011, 0.043, 0.098, 0.173, 0.326,
74  0.524, 0.765, 0.853, 1.146, 1.472,
75  1.708, 2.081, 2.475, 3.323, 3.627,
76  3.779, 4.237, 5.016, 5.647, 5.966,
77  6.793, 7.602, 8.421, 9.249, 10.222,
78  10.923,11.984};
79 
80  fKLevelEnergy[0] = 0.0;
81  fKLevelEnergy[1] = listKEnergy[0];
82  G4int idx = 1;
83  for(G4int i=1; i<nlevels; ++i) {
84  G4int z1 = listK[idx];
85  G4int z2 = listK[i];
86  if(z1+1 < z2) {
87  G4double dz = G4double(z2 - z1);
88  G4double y1 = listKEnergy[idx]/G4double(z1*z1);
89  G4double y2 = listKEnergy[i]/G4double(z2*z2);
90  for(G4int z=z1+1; z<z2; ++z) {
91  fKLevelEnergy[z] = (y1 + (y2 - y1)*(z - z1)/dz)*z*z;
92  }
93  }
94  fKLevelEnergy[z2] = listKEnergy[i];
95  idx = i;
96  }
97  for( G4int i = 0; i<14; ++i) { fLevelEnergy[i] = 0.0; }
98 }
G4EmCaptureCascade::fKLevelEnergy
G4double fKLevelEnergy[93]
Definition: G4EmCaptureCascade.hh:95
y2
Double_t y2[nxs]
Definition: Hadr00/scripts/Style.C:21
y1
Double_t y1[nxs]
Definition: Hadr00/scripts/Style.C:20
G4EmCaptureCascade::theGamma
G4ParticleDefinition * theGamma
Definition: G4EmCaptureCascade.hh:91
G4int
int G4int
Definition: G4Types.hh:78
G4HadronicInteraction::G4HadronicInteraction
G4HadronicInteraction(const G4String &modelName="HadronicModel")
Definition: G4HadronicInteraction.cc:42
G4Gamma::Gamma
static G4Gamma * Gamma()
Definition: G4Gamma.cc:86
G4MuonMinus::MuonMinus
static G4MuonMinus * MuonMinus()
Definition: G4MuonMinus.cc:100
G4Electron::Electron
static G4Electron * Electron()
Definition: G4Electron.cc:94
G4double
double G4double
Definition: G4Types.hh:76
G4EmCaptureCascade::theElectron
G4ParticleDefinition * theElectron
Definition: G4EmCaptureCascade.hh:90
Here is the call graph for this function:

◆ ~G4EmCaptureCascade()

G4EmCaptureCascade::~G4EmCaptureCascade ( )

Definition at line 102 of file G4EmCaptureCascade.cc.

103 {}

◆ G4EmCaptureCascade() [2/2]

G4EmCaptureCascade::G4EmCaptureCascade ( const G4EmCaptureCascade &  )
private

Member Function Documentation

◆ AddNewParticle()

void G4EmCaptureCascade::AddNewParticle ( G4ParticleDefinition *  aParticle,
G4double  kinEnergy 
)
inlineprivate

Definition at line 101 of file G4EmCaptureCascade.hh.

103 {
104  G4DynamicParticle* dp = new G4DynamicParticle(aParticle,
105  G4RandomDirection(),
106  kinEnergy);
107  G4HadSecondary hs(dp);
108  hs.SetTime(fTime);
109  result.AddSecondary(hs);
110 }
G4EmCaptureCascade::result
G4HadFinalState result
Definition: G4EmCaptureCascade.hh:89
G4RandomDirection
G4ThreeVector G4RandomDirection()
Definition: G4RandomDirection.hh:56
G4HadFinalState::AddSecondary
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
Definition: G4HadFinalState.hh:61
Here is the call graph for this function:
Here is the caller graph for this function:

◆ ApplyYourself()

G4HadFinalState * G4EmCaptureCascade::ApplyYourself ( const G4HadProjectile &  aTrack,
G4Nucleus &  targetNucleus 
)
virtual

Implements G4HadronicInteraction.

Definition at line 108 of file G4EmCaptureCascade.cc.

110 {
111  result.Clear();
112  result.SetStatusChange(isAlive);
113  fTime = projectile.GetGlobalTime();
114 
115  G4int Z = targetNucleus.GetZ_asInt();
116  G4int A = targetNucleus.GetA_asInt();
117  G4double massA = G4NucleiProperties::GetNuclearMass(A, Z);
118  G4double mass = fMuMass * massA / (fMuMass + massA) ;
119  G4double e = 13.6 * eV * Z * Z * mass/ electron_mass_c2;
120 
121  // precise corrections of energy only for K-shell
122  fLevelEnergy[0] = fKLevelEnergy[Z];
123  for( G4int i = 2; i < 15; ++i) {
124  fLevelEnergy[i-1] = e/G4double(i*i);
125  }
126 
127  G4int nElec = G4int(Z);
128  G4int nAuger = 1;
129  G4int nLevel = 13;
130  G4double pGamma = Z*Z*Z*Z;
131 
132  // Capture on 14-th level
133  G4double edep = fLevelEnergy[13];
134  AddNewParticle(theElectron,edep);
135  G4double deltaE;
136 
137  // Emit new photon or electron
138  // Simplified model for probabilities
139  // N.C.Mukhopadhyay Phy. Rep. 30 (1977) 1.
140  do {
141 
142  // case of Auger electrons
143  if((nAuger < nElec) && ((pGamma + 10000.0) * G4UniformRand() < 10000.0) ) {
144  ++nAuger;
145  deltaE = fLevelEnergy[nLevel-1] - fLevelEnergy[nLevel];
146  --nLevel;
147  AddNewParticle(theElectron, deltaE);
148 
149  } else {
150 
151  // Case of photon cascade, probabilities from
152  // C.S.Wu and L.Wilets, Ann. Rev. Nuclear Sci. 19 (1969) 527.
153 
154  G4double var = (10.0 + G4double(nLevel - 1) ) * G4UniformRand();
155  G4int iLevel = nLevel - 1 ;
156  if(var > 10.0) iLevel -= G4int(var-10.0) + 1;
157  if( iLevel < 0 ) iLevel = 0;
158  deltaE = fLevelEnergy[iLevel] - fLevelEnergy[nLevel];
159  nLevel = iLevel;
160  AddNewParticle(theGamma, deltaE);
161  }
162  edep += deltaE;
163 
164  // Loop checking, 06-Aug-2015, Vladimir Ivanchenko
165  } while( nLevel > 0 );
166 
167  result.SetLocalEnergyDeposit(edep);
168  return &result;
169 }
G4EmCaptureCascade::result
G4HadFinalState result
Definition: G4EmCaptureCascade.hh:89
G4EmCaptureCascade::fKLevelEnergy
G4double fKLevelEnergy[93]
Definition: G4EmCaptureCascade.hh:95
G4NucleiProperties::GetNuclearMass
static G4double GetNuclearMass(const G4double A, const G4double Z)
Definition: G4NucleiProperties.cc:53
G4EmCaptureCascade::theGamma
G4ParticleDefinition * theGamma
Definition: G4EmCaptureCascade.hh:91
G4Nucleus::GetA_asInt
G4int GetA_asInt() const
Definition: G4Nucleus.hh:109
G4int
int G4int
Definition: G4Types.hh:78
edep
Double_t edep
Definition: brachytherapy/macro.C:11
G4HadFinalState::SetStatusChange
void SetStatusChange(G4HadFinalStateStatus aS)
Definition: G4HadFinalState.hh:67
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:97
A
double A(double temperature)
Definition: G4DNAElectronHoleRecombination.cc:59
Z
Float_t Z
Definition: advanced/microbeam/plot.C:39
python.hepunit.electron_mass_c2
float electron_mass_c2
Definition: hepunit.py:274
G4Nucleus::GetZ_asInt
G4int GetZ_asInt() const
Definition: G4Nucleus.hh:115
G4EmCaptureCascade::AddNewParticle
void AddNewParticle(G4ParticleDefinition *aParticle, G4double kinEnergy)
Definition: G4EmCaptureCascade.hh:101
eV
static const double eV
Definition: G4SIunits.hh:212
isAlive
Definition: G4HadFinalState.hh:42
G4HadFinalState::SetLocalEnergyDeposit
void SetLocalEnergyDeposit(G4double aE)
Definition: G4HadFinalState.hh:82
G4double
double G4double
Definition: G4Types.hh:76
G4EmCaptureCascade::theElectron
G4ParticleDefinition * theElectron
Definition: G4EmCaptureCascade.hh:90
Here is the call graph for this function:

◆ ModelDescription()

void G4EmCaptureCascade::ModelDescription ( std::ostream &  outFile) const
virtual

Reimplemented from G4HadronicInteraction.

Definition at line 173 of file G4EmCaptureCascade.cc.

174 {
175  outFile << "Simulation of electromagnetic cascade from capture level"
176  << " to K-shell of the mesonic atom\n."
177  << "Probabilities of gamma and Auger transitions from\n"
178  << " N.C.Mukhopadhyay Phys. Rep. 30 (1977) 1.\n";
179 }

◆ operator=()

G4EmCaptureCascade& G4EmCaptureCascade::operator= ( const G4EmCaptureCascade &  right)
private

Member Data Documentation

◆ fKLevelEnergy

G4double G4EmCaptureCascade::fKLevelEnergy[93]
private

Definition at line 95 of file G4EmCaptureCascade.hh.

◆ fLevelEnergy

G4double G4EmCaptureCascade::fLevelEnergy[14]
private

Definition at line 94 of file G4EmCaptureCascade.hh.

◆ fMuMass

G4double G4EmCaptureCascade::fMuMass
private

Definition at line 92 of file G4EmCaptureCascade.hh.

◆ fTime

G4double G4EmCaptureCascade::fTime
private

Definition at line 93 of file G4EmCaptureCascade.hh.

◆ result

G4HadFinalState G4EmCaptureCascade::result
private

Definition at line 89 of file G4EmCaptureCascade.hh.

◆ theElectron

G4ParticleDefinition* G4EmCaptureCascade::theElectron
private

Definition at line 90 of file G4EmCaptureCascade.hh.

◆ theGamma

G4ParticleDefinition* G4EmCaptureCascade::theGamma
private

Definition at line 91 of file G4EmCaptureCascade.hh.

The documentation for this class was generated from the following files: