67 static const G4int nlevels = 28;
68 static 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 static 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,
83 for(
G4int i=1; i<nlevels; ++i) {
84 G4int z1 = listK[idx];
123 for(
G4int i=1; i<14; ++i) {
124 fLevelEnergy[i] = e/(
G4double)((i+1)*(i+1));
143 if((nAuger < nElec) && ((pGamma + 10000.0) *
G4UniformRand() < 10000.0) ) {
145 deltaE = fLevelEnergy[nLevel-1] - fLevelEnergy[nLevel];
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];
165 }
while( nLevel > 0 );
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";
virtual void ModelDescription(std::ostream &outFile) const
Float_t y1[n_points_granero]
static G4MuonMinus * MuonMinus()
G4double GetPDGMass() const
virtual ~G4EmCaptureCascade()
Float_t y2[n_points_geant4]
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
G4ParticleDefinition * theGamma
G4double GetGlobalTime() const
static constexpr double electron_mass_c2
double A(double temperature)
static G4double GetNuclearMass(const G4double A, const G4double Z)
static constexpr double eV
static G4Electron * Electron()
G4double fKLevelEnergy[93]
G4double fLevelEnergy[14]
void SetLocalEnergyDeposit(G4double aE)
G4ParticleDefinition * theElectron
void AddNewParticle(G4ParticleDefinition *aParticle, G4double kinEnergy)
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
void SetStatusChange(G4HadFinalStateStatus aS)