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,
83 for(
G4int i=1; i<nlevels; ++i) {
84 G4int z1 = listK[idx];
119 G4double e = 13.6 *
eV * Z * Z * mass/ electron_mass_c2;
123 for(
G4int i = 2; i < 15; ++i) {
143 if((nAuger < nElec) && ((pGamma + 10000.0) *
G4UniformRand() < 10000.0) ) {
155 G4int iLevel = nLevel - 1 ;
156 if(var > 10.0) iLevel -=
G4int(var-10.0) + 1;
157 if( iLevel < 0 ) iLevel = 0;
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";
G4double fKLevelEnergy[93]
static G4double GetNuclearMass(const G4double A, const G4double Z)
void ModelDescription(std::ostream &outFile) const
G4ParticleDefinition * theGamma
void SetStatusChange(G4HadFinalStateStatus aS)
double A(double temperature)
G4double GetGlobalTime() const
G4double fLevelEnergy[14]
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
void AddNewParticle(G4ParticleDefinition *aParticle, G4double kinEnergy)
G4double GetPDGMass() const
static G4MuonMinus * MuonMinus()
void SetLocalEnergyDeposit(G4double aE)
static G4Electron * Electron()
G4ParticleDefinition * theElectron