66 static const G4int nlevels = 28;
67 static const G4int listK[nlevels] = {
68 1, 2, 4, 6, 8, 11, 14, 17, 18, 21, 24,
69 26, 29, 32, 38, 40, 41, 44, 49, 53, 55,
70 60, 65, 70, 75, 81, 85, 92};
71 static const G4double listKEnergy[nlevels] = {
72 0.00275, 0.011, 0.043, 0.098, 0.173, 0.326,
73 0.524, 0.765, 0.853, 1.146, 1.472,
74 1.708, 2.081, 2.475, 3.323, 3.627,
75 3.779, 4.237, 5.016, 5.647, 5.966,
76 6.793, 7.602, 8.421, 9.249, 10.222,
82 for(
G4int i=1; i<nlevels; ++i) {
122 for(
G4int i=1; i<14; ++i) {
123 fLevelEnergy[i] = e/(
G4double)((i+1)*(i+1));
142 if((nAuger < nElec) && ((pGamma + 10000.0) *
G4UniformRand() < 10000.0) ) {
144 deltaE = fLevelEnergy[nLevel-1] - fLevelEnergy[nLevel];
154 G4int iLevel = nLevel - 1 ;
155 if(var > 10.0) iLevel -=
G4int(var-10.0) + 1;
156 if( iLevel < 0 ) iLevel = 0;
157 deltaE = fLevelEnergy[iLevel] - fLevelEnergy[nLevel];
164 }
while( nLevel > 0 );
174 outFile <<
"Simulation of electromagnetic cascade from capture level"
175 <<
" to K-shell of the mesonic atom\n."
176 <<
"Probabilities of gamma and Auger transitions from\n"
177 <<
" N.C.Mukhopadhyay Phys. Rep. 30 (1977) 1.\n";