Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/materials/src/G4AtomicShells.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /materials/src/G4AtomicShells.cc (Version 11.3.0) and /materials/src/G4AtomicShells.cc (Version 9.2.p3)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
                                                   >>  25 //
                                                   >>  26 //
                                                   >>  27 // $Id: G4AtomicShells.cc,v 1.7 2006/10/17 15:15:46 vnivanch Exp $
                                                   >>  28 // GEANT4 tag $Name: geant4-09-02-patch-03 $
                                                   >>  29 //
                                                   >>  30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
 25                                                    31 
 26 // 24-04-98, implementation of the first versi     32 // 24-04-98, implementation of the first version, V. Grichine
 27 // 16-11-98, GetBindingEnergy(Z,ShellNb), M.Ma     33 // 16-11-98, GetBindingEnergy(Z,ShellNb), M.Maire
 28 // 19-07-04, add a protection in GetNumberOfSh     34 // 19-07-04, add a protection in GetNumberOfShells(), mma
 29 // 11-02-05, GetNumberOfElectrons(Z,ShellNb),      35 // 11-02-05, GetNumberOfElectrons(Z,ShellNb), V.Ivanchenko
 30 // 30-04-10, added fIndexOfShells, V.Ivanchenk << 
 31 // 12-10-22, data added for 101 >= Z <=104 by  << 
 32 //           Electron Binding energy from Tabl << 
 33 //           Electron Configuration from CRC H << 
 34 //           and Physics 8th Ed                << 
 35                                                    36 
 36 #include "G4AtomicShells.hh"                   <<  37 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
 37                                                    38 
 38 #include "G4SystemOfUnits.hh"                  <<  39 
                                                   >>  40 #include "G4AtomicShells.hh"
                                                   >>  41 #include <assert.h>
 39                                                    42 
 40 //....oooOO0OOooo........oooOO0OOooo........oo     43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
 41                                                    44 
 42 // clang-format off                            <<  45 const G4int
 43 const G4int G4AtomicShells::fNumberOfShells[10 <<  46 G4AtomicShells::fNumberOfShells[101] =
 44 {                                                  47 {
 45   1,   // nonexistent zero element             <<  48  0 ,  // nonexisting zero element
 46   1,  1,  2,  2,  3,    3,  4,  4,  3,  4, //  << 
 47   5,  5,  6,  6,  6,    6,  6,  7,  8,  8, //  << 
 48   9,  9,  9,  9,  9,    9,  9, 10, 10, 10, //  << 
 49  11, 11, 11, 11, 11,   12, 13, 13, 14, 14, //  << 
 50  14, 14, 14, 14, 14,   15, 15, 15, 16, 16, //  << 
 51  16, 16, 16, 17, 18,   18, 19, 19, 19, 19, //  << 
 52  19, 19, 19, 20, 19,   19, 19, 19, 19, 20, //  << 
 53  21, 21, 21, 21, 21,   21, 21, 21, 22, 22, //  << 
 54  23, 23, 23, 23, 24,   24, 25, 25, 26, 26, //  << 
 55  27, 27, 27, 26, 26,   27, 27, 26, 26, 26, //  << 
 56  27, 27, 28, 28                            //  << 
 57 };                                             << 
 58                                                    49 
 59 // The total shell number is:                  <<  50  1 ,  1 ,  2 ,  2 ,  3 ,    3 ,  4 ,  4 ,  3 ,  4 ,  //  1 - 10
 60 // 1 + G4AtomicShells::TotalNumberOfShells(104 << 
 61                                                    51 
 62 //....oooOO0OOooo........oooOO0OOooo........oo <<  52  5 ,  5 ,  6 ,  6 ,  6 ,    6 ,  6 ,  7 ,  8 ,  8 ,  // 11 - 20
                                                   >>  53 
                                                   >>  54  9 ,  9 ,  9 ,  9 ,  9 ,    9 ,  9 , 10 , 10 , 10 ,  // 21 - 30
                                                   >>  55 
                                                   >>  56 11 , 11 , 11 , 11 , 11 ,   12 , 13 , 13 , 14 , 14 ,  // 31 - 40
                                                   >>  57 
                                                   >>  58 14 , 14 , 14 , 14 , 14 ,   15 , 15 , 15 , 16 , 16 ,  // 41 - 50
                                                   >>  59 
                                                   >>  60 // ----------------------------------------------------------
                                                   >>  61 
                                                   >>  62 16 , 16 , 16 , 17 , 18 ,   18 , 19 , 19 , 19 , 19 ,  // 51 - 60
                                                   >>  63 
                                                   >>  64 19 , 19 , 19 , 20 , 19 ,   19 , 19 , 19 , 19 , 20 ,  // 61 - 70
                                                   >>  65 
                                                   >>  66 21 , 21 , 21 , 21 , 21 ,   21 , 21 , 21 , 22 , 22 ,  // 71 - 80
                                                   >>  67 
                                                   >>  68 23 , 23 , 23 , 23 , 24 ,   24 , 25 , 25 , 26 , 26 ,  // 81 - 90
                                                   >>  69 
                                                   >>  70 27 , 27 , 27 , 26 , 26 ,   27 , 27 , 26 , 26 , 26    // 91 - 100
 63                                                    71 
 64 const G4int G4AtomicShells::fIndexOfShells[105 << 
 65 {                                              << 
 66     0,   // nonexistent zero element           << 
 67     1,    2,    3,    5,    7,     10,   13,   << 
 68    28,   33,   38,   44,   50,     56,   62,   << 
 69    91,  100,  109,  118,  127,    136,  145,   << 
 70   184,  195,  206,  217,  228,    239,  251,   << 
 71   305,  319,  333,  347,  361,    375,  390,   << 
 72   452,  468,  484,  500,  517,    535,  553,   << 
 73   629,  648,  667,  686,  706,    725,  744,   << 
 74   821,  842,  863,  884,  905,    926,  947,   << 
 75  1033, 1056, 1079, 1102, 1125,   1149, 1173, 1 << 
 76  1275, 1302, 1329, 1356, 1382,   1408, 1435, 1 << 
 77  1540, 1567, 1594, 1622                        << 
 78 };                                                 72 };
 79                                                    73 
                                                   >>  74 // The total shell number is:
                                                   >>  75 // 1 + G4AtomicShells::TotalNumberOfShells(100) = 1 + 1539 = 1540 
                                                   >>  76 
 80 //....oooOO0OOooo........oooOO0OOooo........oo     77 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
 81                                                    78 
 82 const G4double G4AtomicShells::fBindingEnergie <<  79 const G4double
                                                   >>  80 G4AtomicShells::fBindingEnergies[1540] =
 83 {                                                  81 {
 84   0.0,   // Nonexisting zero element           <<  82   0.0 ,   // Nonexisting zero element
 85                                                    83 
 86   //  H  ------------------------------------- <<  84   //  H  ---------------------------------------------------------
 87   13.60,                                       <<  85   13.60 ,
 88                                                    86 
 89   //  He ------------------------------------- <<  87   //  He ---------------------------------------------------------
 90   24.59,                                       <<  88   24.59 ,
 91                                                    89 
 92   //  Li  ------------------------------------ <<  90   //  Li  --------------------------------------------------------
 93   58.0,  5.39,                                 <<  91   58.0 ,  5.39 ,
 94                                                    92 
 95   //  Be  ------------------------------------ <<  93   //  Be  --------------------------------------------------------
 96   115.0, 9.32,                                 <<  94   115.0 , 9.32 ,
 97                                                    95 
 98   //  B   ------------------------------------ <<  96   //  B   --------------------------------------------------------
 99   192.0, 12.93, 8.3,                           <<  97   192.0 , 12.93 , 8.3 ,
100                                                    98 
101   //  C  ------------------------------------- <<  99   //  C  ---------------------------------------------------------
102   288.0, 16.59, 11.26,                         << 100   288.0 , 16.59 , 11.26 ,
103                                                   101 
104   //  N  -------------------------------------    102   //  N  ---------------------------------------------------------
105   403.0, 37.3, 20.33, 14.53,                   << 103   403.0 , 37.3 , 20.33 , 14.53 ,
106                                                   104 
107   //  O  -------------------------------------    105   //  O  ---------------------------------------------------------
108   543.1, 41.6, 28.48, 13.62,                   << 106   543.1 , 41.6 , 28.48 , 13.62 ,
109                                                   107 
110   //  F  -------------------------------------    108   //  F  ---------------------------------------------------------
111   696.7, 37.85, 17.42,                         << 109   696.7 , 37.85 , 17.42 ,
112                                                   110 
113   //  Ne -------------------------------------    111   //  Ne ---------------------------------------------------------
114   870.1, 48.47, 21.66, 21.56,                  << 112   870.1 , 48.47 , 21.66 , 21.56 ,
115                                                   113 
116   //  Na -------------------------------------    114   //  Na ---------------------------------------------------------
117   1075.0, 66.0, 34.0, 34.0, 5.14,              << 115   1075.0 , 66.0 , 34.0 , 34.0 , 5.14 ,
118                                                   116 
119   //  Mg -------------------------------------    117   //  Mg ---------------------------------------------------------
120   1308.0, 92.0, 54.0, 54.0, 7.65,              << 118   1308.0 , 92.0 , 54.0 , 54.0 , 7.65 ,
121                                                   119 
122   //  Al -------------------------------------    120   //  Al ---------------------------------------------------------
123   1564.0, 121., 77.0, 77.0, 10.62, 5.99,       << 121   1564.0 , 121. , 77.0 , 77.0 , 10.62 , 5.99 ,
124                                                   122 
125   //  Si -------------------------------------    123   //  Si ---------------------------------------------------------
126   1844.0, 154.0, 104.0, 104.0, 13.46, 8.15,    << 124   1844.0 , 154.0 , 104.0 , 104.0 , 13.46 , 8.15 ,
127                                                   125 
128   //  P  -------------------------------------    126   //  P  ---------------------------------------------------------
129   2148.0, 191.0, 135.0, 134.0, 16.15, 10.49,   << 127   2148.0 , 191.0 , 135.0 , 134.0 , 16.15 , 10.49 ,
130                                                   128 
131   //  S  -------------------------------------    129   //  S  ---------------------------------------------------------
132   2476.0, 232.0, 170.0, 168.0, 20.20, 10.36,   << 130   2476.0 , 232.0 , 170.0 , 168.0 , 20.20 , 10.36 ,
133                                                   131 
134   //  Cl -------------------------------------    132   //  Cl ---------------------------------------------------------
135   2829.0, 277.0, 208.0, 206.0, 24.54, 12.97,   << 133   2829.0 , 277.0 , 208.0 , 206.0 , 24.54 , 12.97 ,
136                                                   134 
137   //  Ar -------------------------------------    135   //  Ar ---------------------------------------------------------
138   3206.3, 326.5, 250.6, 248.5, 29.24, 15.94, 1 << 136   3206.3 , 326.5 , 250.6 , 248.5 , 29.24 , 15.94 , 15.76 ,
139                                                   137 
140   //  K  -------------------------------------    138   //  K  ---------------------------------------------------------
141   3610.0, 381.0, 299.0, 296.0, 37.0, 19.0, 18. << 139   3610.0 , 381.0 , 299.0 , 296.0 , 37.0 , 19.0 , 18.7 , 4.34 ,
142                                                   140 
143   //  Ca -------------------------------------    141   //  Ca ---------------------------------------------------------
144   4041.0, 441.0, 353.0, 349.0, 46.0, 28.0, 28. << 142   4041.0 , 441.0 , 353.0 , 349.0 , 46.0 , 28.0 , 28.0 , 6.11 ,
145                                                   143 
146   //  Sc -------------------------------------    144   //  Sc ---------------------------------------------------------
147   4494.0, 503.0, 408.0, 403.0, 55.0, 33.0, 33. << 145   4494.0 , 503.0 , 408.0 , 403.0 , 55.0 , 33.0 , 33.0 , 8.0 , 6.54 ,
148                                                   146 
149   //  Ti -------------------------------------    147   //  Ti ---------------------------------------------------------
150   4966.0, 567.0, 465.0, 459.0, 64.0, 39.0, 38. << 148   4966.0 , 567.0 , 465.0 , 459.0 , 64.0 , 39.0 , 38.0 , 8.0 , 6.82 ,
151                                                   149 
152   //  V  -------------------------------------    150   //  V  ---------------------------------------------------------
153   5465.0, 633.0, 525.0, 518.0, 72.0, 44.0, 43. << 151   5465.0 , 633.0 , 525.0 , 518.0 , 72.0 , 44.0 , 43.0 , 8.0 , 6.74 ,
154                                                   152 
155   //  Cr -------------------------------------    153   //  Cr ---------------------------------------------------------
156   5989.0, 702.0, 589.0, 580.0, 80.0, 49.0, 48. << 154   5989.0 , 702.0 , 589.0 , 580.0 , 80.0 , 49.0 , 48.0 , 8.25 , 6.77 ,
157                                                   155 
158   //  Mn -------------------------------------    156   //  Mn ---------------------------------------------------------
159   6539.0, 755.0, 656.0, 645.0, 89.0, 55.0, 53. << 157   6539.0 , 755.0 , 656.0 , 645.0 , 89.0 , 55.0 , 53.0 , 9.0 , 7.43 ,
160                                                   158 
161   //  Fe -------------------------------------    159   //  Fe ---------------------------------------------------------
162   7112.0, 851.0, 726.0, 713.0, 98.0, 61.0, 59. << 160   7112.0 , 851.0 , 726.0 , 713.0 , 98.0 , 61.0 , 59.0 , 9.0 , 7.87 ,
163                                                   161 
164   //  Co -------------------------------------    162   //  Co ---------------------------------------------------------
165   7709.0, 931.0, 800.0, 785.0, 107.0, 68.0, 66 << 163   7709.0 , 931.0 , 800.0 , 785.0 , 107.0 , 68.0 , 66.0 , 9.0 , 7.86 ,
166                                                   164 
167   //  Ni -------------------------------------    165   //  Ni ---------------------------------------------------------
168   8333.0, 1015.0, 877.0, 860.0, 117.0, 75.0, 7 << 166   8333.0 , 1015.0 , 877.0 , 860.0 , 117.0 , 75.0 , 73.0 , 10.0 , 10.0 ,
169      7.64,                                     << 167      7.64 ,
170                                                   168 
171   //  Cu -------------------------------------    169   //  Cu ---------------------------------------------------------
172   8979.0, 1103.0, 958.0, 938.0, 127.0, 82.0, 8 << 170   8979.0 , 1103.0 , 958.0 , 938.0 , 127.0 , 82.0 , 80.0 , 11.0 , 10.4 ,
173      7.73,                                     << 171      7.73 ,
174                                                   172 
175   //  Zn -------------------------------------    173   //  Zn ---------------------------------------------------------
176   9659.0, 1198.0, 1047.0, 1024.0, 141.0, 94.0, << 174   9659.0 , 1198.0 , 1047.0 , 1024.0 , 141.0 , 94.0 , 91.0 , 12.0 , 11.2 ,
177      9.39,                                     << 175      9.39 ,
178                                                   176 
179   //  Ga -------------------------------------    177   //  Ga ---------------------------------------------------------
180   10367.0, 1302.0, 1146.0, 1119.0, 162.0, 111. << 178   10367.0 , 1302.0 , 1146.0 , 1119.0 , 162.0 , 111.0 , 107.0 , 21.0 ,
181      20.0,   11.0,    6.0,                     << 179      20.0 ,   11.0 ,    6.0 ,
182                                                   180 
183   //  Ge -------------------------------------    181   //  Ge ---------------------------------------------------------
184   11103.0, 1413.0, 1251.0, 1220.0, 184.0, 130. << 182   11103.0 , 1413.0 , 1251.0 , 1220.0 , 184.0 , 130.0 , 125.0 , 33.0 ,
185      32.0,   14.3,    7.9,                     << 183      32.0 ,   14.3 ,    7.9 ,
186                                                   184 
187   //  As -------------------------------------    185   //  As ---------------------------------------------------------
188   11867.0, 1531.0, 1362.0, 1327.0, 208.0, 151. << 186   11867.0 , 1531.0 , 1362.0 , 1327.0 , 208.0 , 151.0 , 145.0 , 46.0 ,
189      45.0,   17.0,    9.81,                    << 187      45.0 ,   17.0 ,    9.81 ,
190                                                   188 
191   //  Se -------------------------------------    189   //  Se ---------------------------------------------------------
192   12658.0, 1656.0, 1479.0, 1439.0, 234.0, 173. << 190   12658.0 , 1656.0 , 1479.0 , 1439.0 , 234.0 , 173.0 , 166.0 , 61.0 ,
193      60.0,   20.15,   9.75,                    << 191      60.0 ,   20.15 ,   9.75 ,
194                                                   192 
195   //  Br -------------------------------------    193   //  Br ---------------------------------------------------------
196   13474.0, 1787.0, 1602.0, 1556.0, 262.0, 197. << 194   13474.0 , 1787.0 , 1602.0 , 1556.0 , 262.0 , 197.0 , 189.0 , 77.0 ,
197      76.0,   23.8,   11.85,                    << 195      76.0 ,   23.8 ,   11.85 ,
198                                                   196 
199   //  Kr -------------------------------------    197   //  Kr ---------------------------------------------------------
200   14326.0, 1924.6, 1730.9, 1678.4, 292.8, 222. << 198   14326.0 , 1924.6 , 1730.9 , 1678.4 , 292.8 , 222.2 , 214.4 , 95.0 ,
201      93.8,   27.51,  14.65,  14.0,             << 199      93.8 ,   27.51 ,  14.65 ,  14.0 ,
202                                                   200 
203   //  Rb -------------------------------------    201   //  Rb ---------------------------------------------------------
204   15200.0, 2068.0, 1867.0, 1807.0, 325.0, 251. << 202   15200.0 , 2068.0 , 1867.0 , 1807.0 , 325.0 , 251.0 , 242.0 , 116.0 ,
205     114.0,   32.0,   16.0,   15.3,   4.18,     << 203     114.0 ,   32.0 ,   16.0 ,   15.3 ,   4.18 ,
206                                                   204 
207   //  Sr -------------------------------------    205   //  Sr ---------------------------------------------------------
208   16105.0, 2219.0, 2010.0, 1943.0, 361.0, 283. << 206   16105.0 , 2219.0 , 2010.0 , 1943.0 , 361.0 , 283.0 , 273.0 , 139.0 ,
209     137.0,   40.0,   23.0,   22.0,   5.69,     << 207     137.0 ,   40.0 ,   23.0 ,   22.0 ,   5.69 ,
210                                                   208 
211   //  Y  -------------------------------------    209   //  Y  ---------------------------------------------------------
212   17038.0, 2375.0, 2158.0, 2083.0, 397.0, 315. << 210   17038.0 , 2375.0 , 2158.0 , 2083.0 , 397.0 , 315.0 , 304.0 , 163.0 ,
213     161.0,   48.0,   30.0,   29.0,   6.48,   6 << 211     161.0 ,   48.0 ,   30.0 ,   29.0 ,   6.38 ,  6.48 ,
214                                                   212 
215   //  Zr -------------------------------------    213   //  Zr ---------------------------------------------------------
216   17998.0, 2536.0, 2311.0, 2227.0, 434.0, 348. << 214   17998.0 , 2536.0 , 2311.0 , 2227.0 , 434.0 , 348.0 , 335.0 , 187.0 ,
217     185.0,   56.0,   35.0,   33.0,   8.61,  6. << 215     185.0 ,   56.0 ,   35.0 ,   33.0 ,   8.61 ,  6.84 ,
218                                                   216 
219   //  Nb -------------------------------------    217   //  Nb ----------------------------------------------------------
220   18986.0, 2702.0, 2469.0, 2375.0, 472.0, 382. << 218   18986.0 , 2702.0 , 2469.0 , 2375.0 , 472.0 , 382.0 , 367.0 , 212.0 ,
221     209.0,   62.0,   40.0,   38.0,   7.17,  6. << 219     209.0 ,   62.0 ,   40.0 ,   38.0 ,   7.17 ,  6.88 ,
222                                                   220 
223   //  Mo -------------------------------------    221   //  Mo -----------------------------------------------------------
224   20000.0, 2872.0, 2632.0, 2527.0, 511.0, 416. << 222   20000.0 , 2872.0 , 2632.0 , 2527.0 , 511.0 , 416.0 , 399.0 , 237.0 ,
225     234.0,   68.0,   45.0,   42.0,   8.56,  7. << 223     234.0 ,   68.0 ,   45.0 ,   42.0 ,   8.56 ,  7.10 ,
226                                                   224 
227   //  Tc -------------------------------------    225   //  Tc -----------------------------------------------------------
228   21044.0, 3048.0, 2800.0, 2683.0, 551.0, 451. << 226   21044.0 , 3048.0 , 2800.0 , 2683.0 , 551.0 , 451.0 , 432.0 , 263.0 ,
229     259.0,   74.0,   49.0,   45.0,   8.6,   7. << 227     259.0 ,   74.0 ,   49.0 ,   45.0 ,   8.6 ,   7.28 ,
230                                                   228 
231   //  Ru -------------------------------------    229   //  Ru -----------------------------------------------------------
232   22117.0, 3230.0, 2973.0, 2844.0, 592.0, 488. << 230   22117.0 , 3230.0 , 2973.0 , 2844.0 , 592.0 , 488.0 , 466.0 , 290.0 ,
233     286.0,   81.0,   53.0,   49.0,   8.5,   7. << 231     286.0 ,   81.0 ,   53.0 ,   49.0 ,   8.5 ,   7.37 ,
234                                                   232 
235   //  Rh -------------------------------------    233   //  Rh -----------------------------------------------------------
236   23220.0, 3418.0, 3152.0, 3010.0, 634.0, 526. << 234   23220.0 , 3418.0 , 3152.0 , 3010.0 , 634.0 , 526.0 , 501.0 , 318.0 ,
237     313.0,   87.0,   58.0,   53.0,   9.56,  7. << 235     313.0 ,   87.0 ,   58.0 ,   53.0 ,   9.56 ,  7.46 ,
238                                                   236 
239   //  Pd -------------------------------------    237   //  Pd -----------------------------------------------------------
240   24350.0, 3611.0, 3337.0, 3180.0, 677.0, 565. << 238   24350.0 , 3611.0 , 3337.0 , 3180.0 , 677.0 , 565.0 , 537.0 , 347.0 ,
241     342.0,   93.0,   63.0,   57.0,   8.78,  8. << 239     342.0 ,   93.0 ,   63.0 ,   57.0 ,   8.78 ,  8.34 ,  7.52 ,
242                                                   240 
243   //  Ag -------------------------------------    241   //  Ag -----------------------------------------------------------
244   25514.0, 3812.0, 3530.0, 3357.0, 724.0, 608. << 242   25514.0 , 3812.0 , 3530.0 , 3357.0 , 724.0 , 608.0 , 577.0 , 379.0 ,
245     373.0,  101.0,   69.0,   63.0,  11.0,  10. << 243     373.0 ,  101.0 ,   69.0 ,   63.0 ,  11.0 ,  10.0 ,   7.58 ,
246                                                   244 
247   //  Cd -------------------------------------    245   //  Cd ------------------------------------------------------------
248   26711.0, 4022.0, 3732.0, 3542.0, 775.0, 655. << 246   26711.0 , 4022.0 , 3732.0 , 3542.0 , 775.0 , 655.0 , 621.0 , 415.0 ,
249     408.0,  112.0,   78.0,   71.0,  14.0,  13. << 247     408.0 ,  112.0 ,   78.0 ,   71.0 ,  14.0 ,  13.0 ,   8.99 ,
250                                                   248 
251   //  In -------------------------------------    249   //  In ------------------------------------------------------------
252   27940.0, 4242.0, 3943.0, 3735.0, 830.0, 707. << 250   27940.0 , 4242.0 , 3943.0 , 3735.0 , 830.0 , 707.0 , 669.0 , 455.0 ,
253     447.0,  126.0,   90.0,   82.0,  21.0,  20. << 251     447.0 ,  126.0 ,   90.0 ,   82.0 ,  21.0 ,  20.0 ,  10.0 ,   5.79 ,
254                                                   252 
255   //  Sn -------------------------------------    253   //  Sn ------------------------------------------------------------
256   29200.0, 4469.0, 4160.0, 3933.0, 888.0, 761. << 254   29200.0 , 4469.0 , 4160.0 , 3933.0 , 888.0 , 761.0 , 719.0 , 497.0 ,
257     489.0,  141.0,  102.0,   93.0,  29.0,  28. << 255     489.0 ,  141.0 ,  102.0 ,   93.0 ,  29.0 ,  28.0 ,  12.0 ,   7.34 ,
258                                                   256 
259   //  Sb -------------------------------------    257   //  Sb ------------------------------------------------------------
260   30419.0, 4698.0, 4385.0, 4137.0, 949.0, 817. << 258   30419.0 , 4698.0 , 4385.0 , 4137.0 , 949.0 , 817.0 , 771.0 , 542.0 ,
261     533.0,  157.0,  114.0,  104.0,  38.0,  37. << 259     533.0 ,  157.0 ,  114.0 ,  104.0 ,  38.0 ,  37.0 ,  15.0 ,   8.64 ,
262                                                   260 
263   //  Te -------------------------------------    261   //  Te ------------------------------------------------------------
264   31814.0, 4939.0, 4612.0, 4347.0, 1012.0, 876 << 262   31814.0 , 4939.0 , 4612.0 , 4347.0 , 1012.0 , 876.0 , 825.0 , 589.0 ,
265     578.0,  174.0,  127.0,  117.0,   48.0,  46 << 263     578.0 ,  174.0 ,  127.0 ,  117.0 ,   48.0 ,  46.0 ,  17.84 ,  9.01 ,
266                                                   264 
267   //  I  -------------------------------------    265   //  I  ------------------------------------------------------------
268   33169.0, 5188.0, 4852.0, 4557.0, 1078.0, 937 << 266   33169.0 , 5188.0 , 4852.0 , 4557.0 , 1078.0 , 937.0 , 881.0 , 638.0 ,
269     626.0,  193.0,  141.0,  131.0,   58.0,  56 << 267     626.0 ,  193.0 ,  141.0 ,  131.0 ,   58.0 ,  56.0 ,  20.61 , 10.45 ,
270                                                   268 
271   //  Xe -------------------------------------    269   //  Xe ------------------------------------------------------------
272   34570.0, 5460.0, 5110.0, 4790.0, 1148.7, 100 << 270   34570.0 , 5460.0 , 5110.0 , 4790.0 , 1148.7 , 1002.1 , 940.6 , 689.0 ,
273     676.4,  213.2,  157.0,  145.5,   69.5,   6 << 271     676.4 ,  213.2 ,  157.0 ,  145.5 ,   69.5 ,   67.5 ,  23.39 , 13.43 ,
274      12.13,                                    << 272      12.13 ,
275                                                   273 
276   //  Cs -------------------------------------    274   //  Cs ------------------------------------------------------------
277   35985.0, 5714.0, 5359.0, 5012.0, 1220.0, 106 << 275   35985.0 , 5714.0 , 5359.0 , 5012.0 , 1220.0 , 1068.0 , 1000.0 , 742.0 ,
278     728.0,  233.0,  174.0,  164.0,   81.0,   7 << 276     728.0 ,  233.0 ,  174.0 ,  164.0 ,   81.0 ,   79.0 ,   25.0 ,  14.0 ,
279      12.3,    3.89,                            << 277      12.3 ,    3.89 ,
280                                                   278 
281   //  Ba -------------------------------------    279   //  Ba -------------------------------------------------------------
282   37441.0, 5989.0, 5624.0, 5247.0, 1293.0, 113 << 280   37441.0 , 5989.0 , 5624.0 , 5247.0 , 1293.0 , 1138.0 , 1063.0 , 797.0 ,
283     782.0,  254.0,  193.0,  181.0,   94.0,   9 << 281     782.0 ,  254.0 ,  193.0 ,  181.0 ,   94.0 ,   92.0 ,   31.0 ,  18.0 ,
284      16.0,    5.21,                            << 282      16.0 ,    5.21 ,
285                                                   283 
286   //  La -------------------------------------    284   //  La -------------------------------------------------------------
287   38925.0, 6266.0, 5891.0, 5483.0, 1365.0, 120 << 285   38925.0 , 6266.0 , 5891.0 , 5483.0 , 1365.0 , 1207.0 , 1124.0 , 851.0 ,
288     834.0,  273.0,  210.0,  196.0,  105.0,  10 << 286     834.0 ,  273.0 ,  210.0 ,  196.0 ,  105.0 ,  103.0 ,   36.0 ,  22.0 ,
289      19.0,    5.75,   5.58,                    << 287      19.0 ,    5.75 ,   5.58 ,
290                                                   288 
291   //  Ce -------------------------------------    289   //  Ce -------------------------------------------------------------
292   40443.0, 6548.0, 6164.0, 5723.0, 1437.0, 127 << 290   40443.0 , 6548.0 , 6164.0 , 5723.0 , 1437.0 , 1275.0 , 1184.0 , 903.0 ,
293     885.0,  291.0,  225.0,  209.0,  114.0,  11 << 291     885.0 ,  291.0 ,  225.0 ,  209.0 ,  114.0 ,  111.0 ,    6.0 ,  39.0 ,
294      22.0,    6.0,    5.65,                    << 292      25.0 ,   22.0 ,    5.65 ,
295                                                   293 
296   //  Pr -------------------------------------    294   //  Pr -------------------------------------------------------------
297   41991.0, 6835.0, 6440.0, 5964.0, 1509.0, 134 << 295   41991.0 , 6835.0 , 6440.0 , 5964.0 , 1509.0 , 1342.0 , 1244.0 , 954.0 ,
298     934.0,  307.0,  238.0,  220.0,  121.0,  11 << 296     934.0 ,  307.0 ,  238.0 ,  220.0 ,  121.0 ,  117.0 ,    6.0 ,  41.0 ,
299      27.0,   24.0,    6.0,  5.42 ,             << 297      27.0 ,   24.0 ,    5.42 ,
300                                                   298 
301   //  Nd -------------------------------------    299   //  Nd -------------------------------------------------------------
302   43569.0, 7126.0, 6722.0, 6208.0, 1580.0, 140 << 300   43569.0 , 7126.0 , 6722.0 , 6208.0 , 1580.0 , 1408.0 , 1303.0 , 1005.0 ,
303     983.0,  321.0,  250.0,  230.0,  126.0,  12 << 301     983.0 ,  321.0 ,  250.0 ,  230.0 ,  126.0 ,  122.0 ,    6.0 ,   42.0 ,
304      28.0,   25.0,    6.0,  5.49,              << 302      28.0 ,   25.0 ,    5.49 ,
305                                                   303 
306   //  Pm -------------------------------------    304   //  Pm -------------------------------------------------------------
307   45184.0, 7428.0, 7013.0, 6459.0, 1653.0, 147 << 305   45184.0 , 7428.0 , 7013.0 , 6459.0 , 1653.0 , 1476.0 , 1362.0 , 1057.0 ,
308    1032.0,  325.0,  261.0,  240.0,  131.0,  12 << 306    1032.0 ,  325.0 ,  261.0 ,  240.0 ,  131.0 ,  127.0 ,    6.0 ,   43.0 ,
309      28.0,   25.0,    6.0,   5.55,             << 307      28.0 ,   25.0 ,    5.55 ,
310                                                   308 
311   //  Sm -------------------------------------    309   //  Sm -------------------------------------------------------------
312   46834.0, 7737.0, 7312.0, 6716.0, 1728.0, 154 << 310   46834.0 , 7737.0 , 7312.0 , 6716.0 , 1728.0 , 1546.0 , 1422.0 , 1110.0 ,
313    1083.0,  349.0,  273.0,  251.0,  137.0,  13 << 311    1083.0 ,  349.0 ,  273.0 ,  251.0 ,  137.0 ,  132.0 ,    6.0 ,   44.0 ,
314      29.0,   25.0,    6.0,    5.63,            << 312      29.0 ,   25.0 ,    5.63 ,
315                                                   313 
316   //  Eu -------------------------------------    314   //  Eu -------------------------------------------------------------
317   48519.0, 8052.0, 7617.0, 6977.0, 1805.0, 161 << 315   48519.0 , 8052.0 , 7617.0 , 6977.0 , 1805.0 , 1618.0 , 1484.0 , 1164.0 ,
318    1135.0,  364.0,  286.0,  262.0,  143.0,  13 << 316    1135.0 ,  364.0 ,  286.0 ,  262.0 ,  143.0 ,  137.0 ,    6.0 ,   45.0 ,
319      30.0,   26.0,    6.0,    5.68,            << 317      30.0 ,   26.0 ,    5.68 ,
320                                                   318 
321   //  Gd -------------------------------------    319   //  Gd -------------------------------------------------------------
322   50239.0, 8376.0, 7930.0, 7243.0, 1884.0, 169 << 320   50239.0 , 8376.0 , 7930.0 , 7243.0 , 1884.0 , 1692.0 , 1547.0 , 1220.0 ,
323    1189.0,  380.0,  300.0,  273.0,  150.0,  14 << 321    1189.0 ,  380.0 ,  300.0 ,  273.0 ,  150.0 ,  143.0 ,    6.0 ,   46.0 ,
324      31.0,   27.0,    6.16,   6.0,    6.0,     << 322      31.0 ,   27.0 ,    6.0 ,    6.16 ,
325                                                   323 
326   //  Tb -------------------------------------    324   //  Tb -------------------------------------------------------------
327   51996.0, 8708.0, 8252.0, 7514.0, 1965.0, 176 << 325   51996.0 , 8708.0 , 8252.0 , 7514.0 , 1965.0 , 1768.0 , 1612.0 , 1277.0 ,
328    1243.0,  398.0,  315.0,  285.0,  157.0,  15 << 326    1243.0 ,  398.0 ,  315.0 ,  285.0 ,  157.0 ,  150.0 ,    6.0 ,   48.0 ,
329      32.0,   28.0,    6.0,    5.85,            << 327      32.0 ,   28.0 ,    5.85 ,
330                                                   328 
331   //  Dy -------------------------------------    329   //  Dy -------------------------------------------------------------
332   53789.0, 9046.0, 8581.0, 7790.0, 2048.0, 184 << 330   53789.0 , 9046.0 , 8581.0 , 7790.0 , 2048.0 , 1846.0 , 1678.0 , 1335.0 ,
333    1298.0,  416.0,  331.0,  297.0,  164.0,  15 << 331    1298.0 ,  416.0 ,  331.0 ,  297.0 ,  164.0 ,  157.0 ,    6.0 ,   50.0 ,
334      33.0,   28.0,    6.0,    5.93,            << 332      33.0 ,   28.0 ,    5.93 ,
335                                                   333 
336   //  Ho -------------------------------------    334   //  Ho -------------------------------------------------------------
337   55618.0, 9394.0, 8918.0, 8071.0, 2133.0, 192 << 335   55618.0 , 9394.0 , 8918.0 , 8071.0 , 2133.0 , 1926.0 , 1746.0 , 1395.0 ,
338    1354.0,  434.0,  348.0,  310.0,  172.0,  16 << 336    1354.0 ,  434.0 ,  348.0 ,  310.0 ,  172.0 ,  164.0 ,    6.0 ,   52.0 ,
339      34.0,   29.0,    6.02,   6.0,             << 337      34.0 ,   29.0 ,    6.02 ,
340                                                   338 
341   //  Er -------------------------------------    339   //  Er -------------------------------------------------------------
342   57486.0, 9751.0, 9264.0, 8358.0, 2220.0, 200 << 340   57486.0 , 9751.0 , 9264.0 , 8358.0 , 2220.0 , 2008.0 , 1815.0 , 1456.0 ,
343    1412.0,  452.0,  365.0,  323.0,  181.0,  17 << 341    1412.0 ,  452.0 ,  365.0 ,  323.0 ,  181.0 ,  172.0 ,    6.0 ,   54.0 ,
344      35.0,   30.0,    6.10,   6.0,             << 342      35.0 ,   30.0 ,    6.10 ,
345                                                   343 
346   //  Tu -------------------------------------    344   //  Tu -------------------------------------------------------------
347   59390.0, 10116.0, 9617.0, 8648.0, 2309.0, 20 << 345   59390.0 , 10116.0 , 9617.0 , 8648.0 , 2309.0 , 2092.0 , 1885.0 , 1518.0 ,
348    1471.0,   471.0,  382.0,  336.0,  190.0,  1 << 346    1471.0 ,   471.0 ,  382.0 ,  336.0 ,  190.0 ,  181.0 ,    7.0 ,   56.0 ,
349      36.0,    30.0,    7.0,    6.18,           << 347      36.0 ,    30.0 ,    6.18 ,
350                                                   348 
351   //  Yb -------------------------------------    349   //  Yb -------------------------------------------------------------
352   61332.0, 10486.0, 9978.0, 8944.0, 2401.0, 21 << 350   61332.0 , 10486.0 , 9978.0 , 8944.0 , 2401.0 , 2178.0 , 1956.0 , 1580.0 ,
353    1531.0,   490.0,  399.0,  349.0,  200.0,  1 << 351    1531.0 ,   490.0 ,  399.0 ,  349.0 ,  200.0 ,  190.0 ,    8.0 ,    7.0 ,
354      58.0,    37.0,   31.0,    8.0,    7.0,    << 352      58.0 ,    37.0 ,   31.0 ,    6.25 ,
355                                                   353 
356   //  Lu -------------------------------------    354   //  Lu -------------------------------------------------------------
357   63314.0, 10870.0, 10349.0, 9244.0, 2499.0, 2 << 355   63314.0 , 10870.0 , 10349.0 , 9244.0 , 2499.0 , 2270.0 , 2032.0 , 1647.0 ,
358    1596.0,   514.0,   420.0,  366.0,  213.0,   << 356    1596.0 ,   514.0 ,   420.0 ,  366.0 ,  213.0 ,  202.0 ,   13.0 ,   12.0 ,
359      62.0,    39.0,    32.0,   13.0,   12.0,   << 357      62.0 ,    39.0 ,    32.0 ,    6.6 ,    7.0 ,
360                                                   358 
361   //  Hf -------------------------------------    359   //  Hf -------------------------------------------------------------
362   65351.0, 11271.0, 10739.0, 9561.0, 2604.0, 2 << 360   65351.0 , 11271.0 , 10739.0 , 9561.0 , 2604.0 , 2369.0 , 2113.0 , 1720.0 ,
363    1665.0,   542.0,   444.0,  386.0,  229.0,   << 361    1665.0 ,   542.0 ,   444.0 ,  386.0 ,  229.0 ,  217.0 ,   21.0 ,   20.0 ,
364      68.0,    43.0,    35.0,   21.0,   20.0,   << 362      68.0 ,    43.0 ,    35.0 ,    7.0 ,    7.5 ,
365                                                   363 
366   //  Ta -------------------------------------    364   //  Ta -------------------------------------------------------------
367   67416.0, 11682.0, 11136.0, 9881.0, 2712.0, 2 << 365   67416.0 , 11682.0 , 11136.0 , 9881.0 , 2712.0 , 2472.0 , 2197.0 , 1796.0 ,
368    1737.0,   570.0,   469.0,  407.0,  245.0,   << 366    1737.0 ,   570.0 ,   469.0 ,  407.0 ,  245.0 ,  232.0 ,   30.0 ,   28.0 ,
369      74.0,    47.0,    38.0,   30.0,   28.0,   << 367      74.0 ,    47.0 ,    38.0 ,    8.3 ,    7.9 ,
370                                                   368 
371   //  W  -------------------------------------    369   //  W  -------------------------------------------------------------
372   69525.0, 12100.0, 11544.0, 10207.0, 2823.0,  << 370   69525.0 , 12100.0 , 11544.0 , 10207.0 , 2823.0 , 2577.0 , 2283.0 , 1874.0 ,
373    1811.0,   599.0,   495.0,   428.0,  261.0,  << 371    1811.0 ,   599.0 ,   495.0 ,   428.0 ,  261.0 ,  248.0 ,   38.0 ,   36.0 ,
374      80.0,    51.0,    41.0,    38.0,   36.0,  << 372      80.0 ,    51.0 ,    41.0 ,     9.0 ,    8.0 ,
375                                                   373 
376   //  Re -------------------------------------    374   //  Re -------------------------------------------------------------
377   71676.0, 12527.0, 11959.0, 10535.0, 2937.0,  << 375   71676.0 , 12527.0 , 11959.0 , 10535.0 , 2937.0 , 2686.0 , 2371.0 , 1953.0 ,
378    1887.0,   629.0,   522.0,   450.0,  278.0,  << 376    1887.0 ,   629.0 ,   522.0 ,   450.0 ,  278.0 ,  264.0 ,   47.0 ,   45.0 ,
379      86.0,    56.0,    47.0,   45.0,    45.0,  << 377      86.0 ,    56.0 ,    45.0 ,     9.6 ,    7.9 ,
380                                                   378 
381   //  Os -------------------------------------    379   //  Os -------------------------------------------------------------
382   73871.0, 12968.0, 12385.0, 10871.0, 3054.0,  << 380   73871.0 , 12968.0 , 12385.0 , 10871.0 , 3054.0 , 2797.0 , 2461.0 , 2035.0 ,
383    1964.0,   660.0,   551.0,   473.0,  295.0,  << 381    1964.0 ,   660.0 ,   551.0 ,   473.0 ,  295.0 ,  280.0 ,   56.0 ,   54.0 ,
384      92.0,    61.0,    56.0,    54.0,   49.0,  << 382      92.0 ,    61.0 ,    49.0 ,     9.6 ,    8.5 ,
385                                                   383 
386   //  Ir -------------------------------------    384   //  Ir -------------------------------------------------------------
387   76111.0, 13419.0, 12824.0, 11215.0, 3175.0,  << 385   76111.0 , 13419.0 , 12824.0 , 11215.0 , 3175.0 , 2912.0 , 2554.0 , 2119.0 ,
388    2044.0,   693.0,   581.0,   497.0,  314.0,  << 386    2044.0 ,   693.0 ,   581.0 ,   497.0 ,  314.0 ,  298.0 ,   67.0 ,   64.0 ,
389      99.0,    67.0,    66.0,    64.0,   53.0,  << 387      99.0 ,    66.0 ,    53.0 ,     9.6 ,    9.1 ,
390                                                   388 
391   //  Pt -------------------------------------    389   //  Pt -------------------------------------------------------------
392   78395.0, 13880.0, 13273.0, 11564.0, 3300.0,  << 390   78395.0 , 13880.0 , 13273.0 , 11564.0 , 3300.0 , 3030.0 , 2649.0 , 2206.0 ,
393    2126.0,   727.0,   612.0,   522.0,  335.0,  << 391    2126.0 ,   727.0 ,   612.0 ,   522.0 ,  335.0 ,  318.0 ,   78.0 ,   75.0 ,
394     106.0,    78.0,    75.0,    71.0,   57.0,  << 392     106.0 ,    71.0 ,    57.0 ,     9.6 ,    9.0 ,
395                                                   393 
396   //  Au -------------------------------------    394   //  Au -------------------------------------------------------------
397   80725.0, 14353.0, 13734.0, 11919.0, 3430.0,  << 395   80725.0 , 14353.0 , 13734.0 , 11919.0 , 3430.0 , 3153.0 , 2748.0 , 2295.0 ,
398    2210.0,   764.0,   645.0,   548.0,  357.0,  << 396    2210.0 ,   764.0 ,   645.0 ,   548.0 ,  357.0 ,  339.0 ,   91.0 ,   87.0 ,
399     114.0,    91.0,    87.0,    76.0,   61.0,  << 397     114.0 ,    76.0 ,    61.0 ,    12.5 ,   11.1 ,    9.23 ,
400                                                   398 
401   //  Hg -------------------------------------    399   //  Hg -------------------------------------------------------------
402   83102.0, 14839.0, 14209.0, 12284.0, 3567.0,  << 400   83102.0 , 14839.0 , 14209.0 , 12284.0 , 3567.0 , 3283.0 , 2852.0 , 2390.0 ,
403    2300.0,   806.0,   683.0,   579.0,  382.0,  << 401    2300.0 ,   806.0 ,   683.0 ,   579.0 ,  382.0 ,  363.0 ,  107.0 ,  103.0 ,
404     125.0,   107.0,   103.0,    85.0,   68.0,  << 402     125.0 ,    85.0 ,    68.0 ,    14.0 ,   12.0 ,   10.4 ,
405                                                   403 
406   //  Tl -------------------------------------    404   //  Tl -------------------------------------------------------------
407   85530.0, 15347.0, 14698.0, 12658.0, 3710.0,  << 405   85530.0 , 15347.0 , 14698.0 , 12658.0 , 3710.0 , 3420.0 , 2961.0 , 2490.0 ,
408    2394.0,   852.0,   726.0,   615.0,  411.0,  << 406    2394.0 ,   852.0 ,   726.0 ,   615.0 ,  411.0 ,  391.0 ,  127.0 ,  123.0 ,
409     139.0,   127.0,   123.0,    98.0,   79.0,  << 407     139.0 ,    98.0 ,    79.0 ,    21.0 ,   19.0 ,    8.0 ,    6.11 ,
410       6.11,                                    << 
411                                                   408 
412   //  Pb -------------------------------------    409   //  Pb -------------------------------------------------------------
413   88005.0, 15861.0, 15200.0, 13055.0, 3857.0,  << 410   88005.0 , 15861.0 , 15200.0 , 13055.0 , 3857.0 , 3560.0 , 3072.0 , 2592.0 ,
414    2490.0,   899.0,   769.0,   651.0,  441.0,  << 411    2490.0 ,   899.0 ,   769.0 ,   651.0 ,  441.0 ,  419.0 ,  148.0 ,  144.0 ,
415     153.0,   148.0,   144.0,   111.0,   90.0,  << 412     153.0 ,   111.0 ,    90.0 ,    27.0 ,   25.0 ,   10.0 ,    7.42 ,
416       7.42,                                    << 
417                                                   413 
418   //  Bi -------------------------------------    414   //  Bi -------------------------------------------------------------
419   90526.0, 16388.0, 15711.0, 13419.0, 4007.0,  << 415   90526.0 , 16388.0 , 15711.0 , 13419.0 , 4007.0 , 3704.0 , 3185.0 , 2696.0 ,
420    2588.0,   946.0,   813.0,   687.0,  472.0,  << 416    2588.0 ,   946.0 ,   813.0 ,   687.0 ,  472.0 ,  448.0 ,  170.0 ,  165.0 ,
421     167.0,   165.0,   125.0,   101.0,   34.0,  << 417     167.0 ,   125.0 ,   101.0 ,    34.0 ,   32.0 ,   12.0 ,    7.29 ,
422                                                   418 
423   //  Po -------------------------------------    419   //  Po -------------------------------------------------------------
424   93105.0, 16939.0, 16244.0, 13814.0, 4161.0,  << 420   93105.0 , 16939.0 , 16244.0 , 13814.0 , 4161.0 , 3852.0 , 3301.0 , 2802.0 ,
425    2687.0,   994.0,   858.0,   724.0,  503.0,  << 421    2687.0 ,   994.0 ,   858.0 ,   724.0 ,  503.0 ,  478.0 ,  193.0 ,  187.0 ,
426     181.0,   139.0,   112.0,    41.0,   38.0,  << 422     181.0 ,   139.0 ,   112.0 ,    41.0 ,   38.0 ,   15.0 ,    8.43 ,
427                                                   423 
428   //  At -------------------------------------    424   //  At -------------------------------------------------------------
429   95730.0, 17493.0, 16785.0, 14214.0, 4320.0,  << 425   95730.0 , 17493.0 , 16785.0 , 14214.0 , 4320.0 , 4005.0 , 3420.0 , 2910.0 ,
430    2788.0,  1044.0,   904.0,   761.0,  535.0,  << 426    2788.0 ,  1044.0 ,   904.0 ,   761.0 ,  535.0 ,  508.0 ,  217.0 ,  211.0 ,
431     196.0,   153.0,   123.0,    48.0,   44.0,  << 427     196.0 ,   153.0 ,   123.0 ,    48.0 ,   44.0 ,   19.0 ,   11.0 ,    9.3 ,
432                                                   428 
433   //  Rn -------------------------------------    429   //  Rn -------------------------------------------------------------
434   98404.0, 18049.0, 17337.0, 14619.0, 4483.0,  << 430   98404.0 , 18049.0 , 17337.0 , 14619.0 , 4483.0 , 4162.0 , 3452.0 , 3109.0 ,
435    2890.0,  1096.0,   951.0,   798.0,  567.0,  << 431    2890.0 ,  1096.0 ,   951.0 ,   798.0 ,  567.0 ,  538.0 ,  242.0 ,  235.0 ,
436     212.0,   167.0,   134.0,    55.0,   51.0,  << 432     212.0 ,   167.0 ,   134.0 ,    55.0 ,   51.0 ,   24.0 ,   14.0 ,   10.7 ,
437                                                   433 
438   //  Fr -------------------------------------    434   //  Fr -------------------------------------------------------------
439  101137.0, 18639.0, 17907.0, 15031.0, 4652.0,  << 435  101137.0 , 18639.0 , 17907.0 , 15031.0 , 4652.0 , 4324.0 , 3666.0 , 3134.0 ,
440    2998.0,  1153.0,  1003.0,   839.0,  603.0,  << 436    2998.0 ,  1153.0 ,  1003.0 ,   839.0 ,  603.0 ,  572.0 ,  268.0 ,  260.0 ,
441     231.0,   183.0,   147.0,    65.0,   61.0,  << 437     231.0 ,   183.0 ,   147.0 ,    65.0 ,   61.0 ,   33.0 ,   19.0 ,   14.0 ,
442       4.0,                                     << 438       4.0 ,
443                                                   439 
444   //  Ra -------------------------------------    440   //  Ra -------------------------------------------------------------
445  103922.0, 19237.0, 18484.0, 15444.0, 4822.0,  << 441  103922.0 , 19237.0 , 18484.0 , 15444.0 , 4822.0 , 4491.0 , 3793.0 , 3254.0 ,
446    3111.0,  1214.0,  1060.0,   884.0,  642.0,  << 442    3111.0 ,  1214.0 ,  1060.0 ,   884.0 ,  642.0 ,  609.0 ,  296.0 ,  287.0 ,
447     253.0,   201.0,   161.0,    77.0,   73.0,  << 443     253.0 ,   201.0 ,   161.0 ,    77.0 ,   73.0 ,   40.0 ,   25.0 ,   19.0 ,
448       5.28,                                    << 444       5.28 ,
449                                                   445 
450   //  Ac -------------------------------------    446   //  Ac -------------------------------------------------------------
451  106755.0, 19840.0, 19083.0, 15871.0, 5002.0,  << 447  106755.0 , 19840.0 , 19083.0 , 15871.0 , 5002.0 , 4656.0 , 3921.0 , 3374.0 ,
452    3223.0,  1274.0,  1116.0,   928.0,  680.0,  << 448    3223.0 ,  1274.0 ,  1116.0 ,   928.0 ,  680.0 ,  645.0 ,  322.0 ,  313.0 ,
453     274.0,   218.0,   174.0,    88.0,   83.0,  << 449     274.0 ,   218.0 ,   174.0 ,    88.0 ,   83.0 ,   45.0 ,   29.0 ,   22.0 ,
454       6.3,     5.7,                            << 450       5.7 ,     6.3 ,
455                                                   451 
456   //  Th -------------------------------------    452   //  Th -------------------------------------------------------------
457  109651.0, 20472.0, 19693.0, 16300.0, 5182.0,  << 453  109651.0 , 20472.0 , 19693.0 , 16300.0 , 5182.0 , 4830.0 , 4049.0 , 3494.0 ,
458    3335.0,  1333.0,  1171.0,   970.0,  717.0,  << 454    3335.0 ,  1333.0 ,  1171.0 ,   970.0 ,  717.0 ,  679.0 ,  347.0 ,  338.0 ,
459     293.0,   233.0,   185.0,    97.0,   91.0,  << 455     293.0 ,   233.0 ,   185.0 ,    97.0 ,   91.0 ,   50.0 ,   33.0 ,   25.0 ,
460       6.0,     6.0,                            << 456       6.0 ,     6.0 ,
461                                                   457 
462   //  Pa -------------------------------------    458   //  Pa -------------------------------------------------------------
463  112601.0, 21105.0, 20314.0, 16733.0, 5367.0,  << 459  112601.0 , 21105.0 , 20314.0 , 16733.0 , 5367.0 , 5001.0 , 4178.0 , 3613.0 ,
464    3446.0,  1390.0,  1225.0,  1011.0,  752.0,  << 460    3446.0 ,  1390.0 ,  1225.0 ,  1011.0 ,  752.0 ,  712.0 ,  372.0 ,  362.0 ,
465     312.0,   248.0,   195.0,   104.0,   97.0,  << 461     312.0 ,   248.0 ,   195.0 ,   104.0 ,   97.0 ,    6.0 ,   50.0 ,   32.0 ,
466      24.0,     6.0,     6.0,     6.0,          << 462      24.0 ,     6.0 ,     6.0 ,
467                                                   463 
468   //  U  -------------------------------------    464   //  U  -------------------------------------------------------------
469  115606.0, 21757.0, 20948.0, 17166.0, 5548.0,  << 465  115606.0 , 21757.0 , 20948.0 , 17166.0 , 5548.0 , 5182.0 , 4308.0 , 3733.0 ,
470    3557.0,  1446.0,  1278.0,  1050.0,  785.0,  << 466    3557.0 ,  1446.0 ,  1278.0 ,  1050.0 ,  785.0 ,  743.0 ,  396.0 ,  386.0 ,
471     329.0,   261.0,   203.0,   110.0,  101.0,  << 467     329.0 ,   261.0 ,   203.0 ,   110.0 ,  101.0 ,    6.0 ,   52.0 ,   34.0 ,
472      24.0,     6.1,     6.0,     6.0,          << 468      24.0 ,     6.1 ,     6.0 ,
473                                                   469 
474   //  Np -------------------------------------    470   //  Np -------------------------------------------------------------
475  118678.0, 22426.0, 21600.0, 17610.0, 5723.0,  << 471  118678.0 , 22426.0 , 21600.0 , 17610.0 , 5723.0 , 5366.0 , 4440.0 , 3854.0 ,
476    3669.0,  1504.0,  1331.0,  1089.0,  819.0,  << 472    3669.0 ,  1504.0 ,  1331.0 ,  1089.0 ,  819.0 ,  774.0 ,  421.0 ,  410.0 ,
477     346.0,   274.0,   211.0,   116.0,  106.0,  << 473     346.0 ,   274.0 ,   211.0 ,   116.0 ,  106.0 ,    6.0 ,   54.0 ,   35.0 ,
478      25.0,     6.0,     6.0,     6.0,          << 474      25.0 ,     6.0 ,     6.0 ,
479                                                   475 
480   //  Pu -------------------------------------    476   //  Pu -------------------------------------------------------------
481  121818.0, 23097.0, 22266.0, 18056.0, 5933.0,  << 477  121818.0 , 23097.0 , 22266.0 , 18056.0 , 5933.0 , 5541.0 , 4557.0 , 3977.0 ,
482    3783.0,  1563.0,  1384.0,  1128.0,  853.0,  << 478    3783.0 ,  1563.0 ,  1384.0 ,  1128.0 ,  853.0 ,  805.0 ,  446.0 ,  434.0 ,
483     356.0,   287.0,   219.0,   122.0,  111.0,  << 479     356.0 ,   287.0 ,   219.0 ,   122.0 ,  111.0 ,    6.0 ,   53.0 ,   34.0 ,
484      23.0,     6.0,     6.0,                   << 480      23.0 ,     6.0 ,
485                                                   481 
486   //  Am -------------------------------------    482   //  Am -------------------------------------------------------------
487  125027.0, 23773.0, 22944.0, 18504.0, 6121.0,  << 483  125027.0 , 23773.0 , 22944.0 , 18504.0 , 6121.0 , 5710.0 , 4667.0 , 4102.0 ,
488    3898.0,  1623.0,  1439.0,  1167.0,  887.0,  << 484    3898.0 ,  1623.0 ,  1439.0 ,  1167.0 ,  887.0 ,  836.0 ,  467.0 ,  452.0 ,
489     355.0,   301.0,   220.0,   123.0,  112.0,  << 485     355.0 ,   301.0 ,   220.0 ,   123.0 ,  112.0 ,    6.0 ,   54.0 ,   44.0 ,
490      36.0,     6.0,     6.0,                   << 486      36.0 ,     6.0 ,
491                                                   487 
492   //  Cm -------------------------------------    488   //  Cm -------------------------------------------------------------
493  128220.0, 24460.0, 23779.0, 18930.0, 6288.0,  << 489  128220.0 , 24460.0 , 23779.0 , 18930.0 , 6288.0 , 5895.0 , 4797.0 , 4236.0 ,
494    4014.0,  1664.0,  1493.0,  1194.0,  919.0,  << 490    4014.0 ,  1664.0 ,  1493.0 ,  1194.0 ,  919.0 ,  864.0 ,  494.0 ,  479.0 ,
495     384.0,   314.0,   239.0,   126.0,  119.0,  << 491     384.0 ,   314.0 ,   239.0 ,   126.0 ,  119.0 ,   11.0 ,   60.0 ,   39.0 ,
496      27.0,    11.0,     5.0,     6.0,          << 492      27.0 ,     5.0 ,     6.0 ,
497                                                   493 
498   //  Bk -------------------------------------    494   //  Bk -------------------------------------------------------------
499  131590.0, 25275.0, 24385.0, 19452.0, 6556.0,  << 495  131590.0 , 25275.0 , 24385.0 , 19452.0 , 6556.0 , 6147.0 , 4977.0 , 4366.0 ,
500    4133.0,  1729.0,  1554.0,  1236.0,  955.0,  << 496    4133.0 ,  1729.0 ,  1554.0 ,  1236.0 ,  955.0 ,  898.0 ,  520.0 ,  504.0 ,
501     401.0,   329.0,   248.0,   142.0,  124.0,  << 497     401.0 ,   329.0 ,   248.0 ,   142.0 ,  124.0 ,   12.0 ,   63.0 ,   41.0 ,
502      27.0,    12.0,     6.0,     4.0,          << 498      27.0 ,     4.0 ,     6.0 ,
503                                                   499 
504   //  Cf -------------------------------------    500   //  Cf -------------------------------------------------------------
505  135960.0, 26110.0, 25250.0, 19930.0, 6754.0,  << 501  135960.0 , 26110.0 , 25250.0 , 19930.0 , 6754.0 , 6359.0 , 5109.0 , 4492.0 ,
506    4247.0,  1789.0,  1610.0,  1273.0,  987.0,  << 502    4247.0 ,  1789.0 ,  1610.0 ,  1273.0 ,  987.0 ,  925.0 ,  546.0 ,  529.0 ,
507     412.0,   338.0,   251.0,   142.0,  129.0,  << 503     412.0 ,   338.0 ,   251.0 ,   142.0 ,  129.0 ,    9.0 ,   61.0 ,   39.0 ,
508      25.0,     9.0,     6.0,                   << 504      25.0 ,     6.0 ,
509                                                << 505 
510   //  Es ------------------------------------- << 506   //  Es -------------------------------------------------------------
511  139490.0, 26900.0, 26020.0, 20410.0, 6977.0,  << 507  139490.0 , 26900.0 , 26020.0 , 20410.0 , 6977.0 , 6754.0 , 5252.0 , 4630.0 ,
512    4369.0,  1857.0,  1674.0,  1316.0, 1024.0,  << 508    4369.0 ,  1857.0 ,  1674.0 ,  1316.0 , 1024.0 ,  959.0 ,  573.0 ,  554.0 ,
513     429.0,   353.0,   260.0,   148.0,  135.0,  << 509     429.0 ,   353.0 ,   260.0 ,   148.0 ,  135.0 ,    9.0 ,   63.0 ,   40.0 ,
514      25.0,     9.0,     6.0,                   << 510      25.0 ,     6.0 ,
515                                                << 511 
516   //  Fm ------------------------------------- << 512   //  Fm -------------------------------------------------------------
517  143090.0, 27700.0, 26810.0, 20900.0, 7205.0,  << 513  143090.0 , 27700.0 , 26810.0 , 20900.0 , 7205.0 , 6793.0 , 5397.0 , 4766.0 ,
518    4498.0,  1933.0,  1746.0,  1366.0, 1068.0,  << 514    4498.0 ,  1933.0 ,  1746.0 ,  1366.0 , 1068.0 , 1000.0 ,  606.0 ,  587.0 ,
519     453.0,   375.0,   275.0,   160.0,  145.0,  << 515     453.0 ,   375.0 ,   275.0 ,   160.0 ,  145.0 ,   15.0 ,   69.0 ,   45.0 ,
520      29.0,    15.0,     7.0,                   << 516      29.0 ,     7.0 ,
521                                                << 517 
522   // Md -------------------------------------- << 
523  146526.0, 28387.0, 27438.0, 21356.0, 7440.0,  << 
524    4615.0,  2024.0,  1816.0,  1424.0, 1105.0,  << 
525     471.0,   389.0,   272.0,   154.0,  137.0,  << 
526      37.0,    17.0,     5.9,                   << 
527                                                << 
528   // No -------------------------------------- << 
529  149208.0, 29221.0, 28255.0, 21851.0, 7678.0,  << 
530    4741.0,  2097.0,  1885.0,  1469.0, 1145.0,  << 
531     490.0,   406.0,   280.0,   161.0,  142.0,  << 
532      38.0,    18.0,     6.0,                   << 
533                                                << 
534   // Lr -------------------------------------- << 
535  152970.0, 30083.0, 29103.0, 22359.0, 7930.0,  << 
536    4876.0,  2180.0,  1963.0,  1523.0, 1192.0,  << 
537     516.0,   429.0,   296.0,   174.0,  154.0,  << 
538      44.0,    21.0,     3.9,     6.9,          << 
539                                                << 
540   // Rf -------------------------------------- << 
541  156288.0, 30881.0, 29986.0, 22907.0, 8161.0,  << 
542    5014.0,  2237.0,  2035.0,  1554.0, 1233.0,  << 
543     535.0,   448.0,   319.0,   190.0,  171.0,  << 
544      55.0,    33.0,     5.0,     7.5           << 
545 };                                                518 };
546                                                   519 
547 const G4int G4AtomicShells::fNumberOfElectrons << 520 const G4int
                                                   >> 521 G4AtomicShells::fNumberOfElectrons[1540] =
548 {                                                 522 {
549   0,   // nonexistent zero element             << 523   0,   // Nonexisting zero element
550   //  H  -------------------------------------    524   //  H  ---------------------------------------------------------
551   1,                                              525   1,
552   //  He -------------------------------------    526   //  He ---------------------------------------------------------
553   2,                                              527   2,
554   //  Li  ------------------------------------    528   //  Li  --------------------------------------------------------
555   2, 1,                                           529   2, 1,
556   //  Be  ------------------------------------    530   //  Be  --------------------------------------------------------
557   2, 2,                                           531   2, 2,
558   //  B   ------------------------------------    532   //  B   --------------------------------------------------------
559   2, 2, 1,                                        533   2, 2, 1,
560   //  C  -------------------------------------    534   //  C  ---------------------------------------------------------
561   2, 2, 2,                                        535   2, 2, 2,
562   //  N  -------------------------------------    536   //  N  ---------------------------------------------------------
563   2, 2, 2, 1,                                     537   2, 2, 2, 1,
564   //  O  -------------------------------------    538   //  O  ---------------------------------------------------------
565   2, 2, 2, 2,                                     539   2, 2, 2, 2,
566   //  F  -------------------------------------    540   //  F  ---------------------------------------------------------
567   2, 2, 5,                                        541   2, 2, 5,
568   //  Ne -------------------------------------    542   //  Ne ---------------------------------------------------------
569   2, 2, 2, 4,                                     543   2, 2, 2, 4,
570   //  Na -------------------------------------    544   //  Na ---------------------------------------------------------
571   2, 2, 2, 4, 1,                                  545   2, 2, 2, 4, 1,
572   //  Mg -------------------------------------    546   //  Mg ---------------------------------------------------------
573   2, 2, 2, 4, 2,                                  547   2, 2, 2, 4, 2,
574   //  Al -------------------------------------    548   //  Al ---------------------------------------------------------
575   2, 2, 2, 4, 2, 1,                               549   2, 2, 2, 4, 2, 1,
576   //  Si -------------------------------------    550   //  Si ---------------------------------------------------------
577   2, 2, 2, 4, 2, 2,                               551   2, 2, 2, 4, 2, 2,
578   //  P  -------------------------------------    552   //  P  ---------------------------------------------------------
579   2, 2, 2, 4, 2, 3,                               553   2, 2, 2, 4, 2, 3,
580   //  S  -------------------------------------    554   //  S  ---------------------------------------------------------
581   2, 2, 2, 4, 2, 4,                               555   2, 2, 2, 4, 2, 4,
582   //  Cl -------------------------------------    556   //  Cl ---------------------------------------------------------
583   2, 2, 2, 4, 2, 5,                               557   2, 2, 2, 4, 2, 5,
584   //  Ar -------------------------------------    558   //  Ar ---------------------------------------------------------
585   2, 2, 2, 4, 2, 2, 4,                            559   2, 2, 2, 4, 2, 2, 4,
586   //  K  -------------------------------------    560   //  K  ---------------------------------------------------------
587   2, 2, 2, 4, 2, 2, 4, 1,                         561   2, 2, 2, 4, 2, 2, 4, 1,
588   //  Ca -------------------------------------    562   //  Ca ---------------------------------------------------------
589   2, 2, 2, 4, 2, 2, 4, 2,                         563   2, 2, 2, 4, 2, 2, 4, 2,
590   //  Sc -------------------------------------    564   //  Sc ---------------------------------------------------------
591   2, 2, 2, 4, 2, 2, 4, 1, 2,                      565   2, 2, 2, 4, 2, 2, 4, 1, 2,
592   //  Ti -------------------------------------    566   //  Ti ---------------------------------------------------------
593   2, 2, 2, 4, 2, 2, 4, 2, 2,                      567   2, 2, 2, 4, 2, 2, 4, 2, 2,
594   //  V  -------------------------------------    568   //  V  ---------------------------------------------------------
595   2, 2, 2, 4, 2, 2, 4, 3, 2,                      569   2, 2, 2, 4, 2, 2, 4, 3, 2,
596   //  Cr -------------------------------------    570   //  Cr ---------------------------------------------------------
597   2, 2, 2, 4, 2, 2, 4, 4, 2,                      571   2, 2, 2, 4, 2, 2, 4, 4, 2,
598   //  Mn -------------------------------------    572   //  Mn ---------------------------------------------------------
599   2, 2, 2, 4, 2, 2, 4, 5, 2,                      573   2, 2, 2, 4, 2, 2, 4, 5, 2,
600   //  Fe -------------------------------------    574   //  Fe ---------------------------------------------------------
601   2, 2, 2, 4, 2, 2, 4, 6, 2,                      575   2, 2, 2, 4, 2, 2, 4, 6, 2,
602   //  Co -------------------------------------    576   //  Co ---------------------------------------------------------
603   2, 2, 2, 4, 2, 2, 4, 7, 2,                      577   2, 2, 2, 4, 2, 2, 4, 7, 2,
604   //  Ni -------------------------------------    578   //  Ni ---------------------------------------------------------
605   2, 2, 2, 4, 2, 2, 4, 4, 4, 2,                   579   2, 2, 2, 4, 2, 2, 4, 4, 4, 2,
606   //  Cu -------------------------------------    580   //  Cu ---------------------------------------------------------
607   2, 2, 2, 4, 2, 2, 4, 4, 5, 2,                   581   2, 2, 2, 4, 2, 2, 4, 4, 5, 2,
608   //  Zn -------------------------------------    582   //  Zn ---------------------------------------------------------
609   2, 2, 2, 4, 2, 2, 4, 4, 6, 2,                   583   2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
610   //  Ga -------------------------------------    584   //  Ga ---------------------------------------------------------
611   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 1,                585   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 1,
612   //  Ge -------------------------------------    586   //  Ge ---------------------------------------------------------
613   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2,                587   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2,
614   //  As -------------------------------------    588   //  As ---------------------------------------------------------
615   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 3,                589   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 3,
616   //  Se -------------------------------------    590   //  Se ---------------------------------------------------------
617   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 4,                591   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 4,
618   //  Br -------------------------------------    592   //  Br ---------------------------------------------------------
619   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 5,                593   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 5,
620   //  Kr -------------------------------------    594   //  Kr ---------------------------------------------------------
621   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4,             595   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4,
622   //  Rb -------------------------------------    596   //  Rb ---------------------------------------------------------
623   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 1,          597   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 1,
624   //  Sr -------------------------------------    598   //  Sr ---------------------------------------------------------
625   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 2,          599   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 2,
626   //  Y  -------------------------------------    600   //  Y  ---------------------------------------------------------
627   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 2, 1,    << 601   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 1, 2,
628   //  Zr -------------------------------------    602   //  Zr ---------------------------------------------------------
629   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 2, 2,       603   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 2, 2,
630   //  Nb -------------------------------------    604   //  Nb ----------------------------------------------------------
631   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 3, 2,       605   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 3, 2,
632   //  Mo -------------------------------------    606   //  Mo -----------------------------------------------------------
633   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 2,       607   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 2,
634   //  Tc -------------------------------------    608   //  Tc -----------------------------------------------------------
635   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 5, 2,       609   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 5, 2,
636   //  Ru -------------------------------------    610   //  Ru -----------------------------------------------------------
637   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 6, 2,       611   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 6, 2,
638   //  Rh -------------------------------------    612   //  Rh -----------------------------------------------------------
639   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 7, 2,       613   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 7, 2,
640   //  Pd -------------------------------------    614   //  Pd -----------------------------------------------------------
641   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 4, 2,    615   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 4, 2,
642   //  Ag -------------------------------------    616   //  Ag -----------------------------------------------------------
643   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 5, 2,    617   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 5, 2,
644   //  Cd -------------------------------------    618   //  Cd ------------------------------------------------------------
645   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,    619   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,
646   //  In -------------------------------------    620   //  In ------------------------------------------------------------
647   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,    621   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 1,
648   //  Sn -------------------------------------    622   //  Sn ------------------------------------------------------------
649   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,    623   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2,
650   //  Sb -------------------------------------    624   //  Sb ------------------------------------------------------------
651   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,    625   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 3,
652   //  Te -------------------------------------    626   //  Te ------------------------------------------------------------
653   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,    627   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 4,
654   //  I  -------------------------------------    628   //  I  ------------------------------------------------------------
655   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,    629   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 5,
656   //  Xe -------------------------------------    630   //  Xe ------------------------------------------------------------
657   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,    631   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4,
658   //  Cs -------------------------------------    632   //  Cs ------------------------------------------------------------
659   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,    633   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 1,
660   //  Ba -------------------------------------    634   //  Ba -------------------------------------------------------------
661   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,    635   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 2,
662   //  La -------------------------------------    636   //  La -------------------------------------------------------------
663   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,    637   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 1, 2,
664   //  Ce -------------------------------------    638   //  Ce -------------------------------------------------------------
665   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 639   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 2, 4, 2,
666   //  Pr -------------------------------------    640   //  Pr -------------------------------------------------------------
667   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 641   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 3, 2, 2, 4, 2,
668   //  Nd -------------------------------------    642   //  Nd -------------------------------------------------------------
669   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 643   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 4, 2, 2, 4, 2,
670   //  Pm -------------------------------------    644   //  Pm -------------------------------------------------------------
671   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 645   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 5, 2, 2, 4, 2,
672   //  Sm -------------------------------------    646   //  Sm -------------------------------------------------------------
673   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 647   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 2, 2, 4, 2,
674   //  Eu -------------------------------------    648   //  Eu -------------------------------------------------------------
675   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 649   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 7, 2, 2, 4, 2,
676   //  Gd -------------------------------------    650   //  Gd -------------------------------------------------------------
677   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 651   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 7, 2, 2, 4, 1, 2,
678   //  Tb -------------------------------------    652   //  Tb -------------------------------------------------------------
679   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 653   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 9, 2, 2, 4, 2,
680   //  Dy -------------------------------------    654   //  Dy -------------------------------------------------------------
681   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 655   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6,10, 2, 2, 4, 2,
682   //  Ho -------------------------------------    656   //  Ho -------------------------------------------------------------
683   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 657   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6,11, 2, 2, 4, 2,
684   //  Er -------------------------------------    658   //  Er -------------------------------------------------------------
685   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 659   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6,12, 2, 2, 4, 2,
686   //  Tu -------------------------------------    660   //  Tu -------------------------------------------------------------
687   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 661   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6,13, 2, 2, 4, 2,
688   //  Yb -------------------------------------    662   //  Yb -------------------------------------------------------------
689   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 663   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 2,
690   //  Lu -------------------------------------    664   //  Lu -------------------------------------------------------------
691   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 665   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 1, 2,
692   //  Hf -------------------------------------    666   //  Hf -------------------------------------------------------------
693   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 667   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 2, 2,
694   //  Ta -------------------------------------    668   //  Ta -------------------------------------------------------------
695   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 669   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 3, 2,
696   //  W  -------------------------------------    670   //  W  -------------------------------------------------------------
697   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 671   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 2,
698   //  Re -------------------------------------    672   //  Re -------------------------------------------------------------
699   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 673   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 5, 2,
700   //  Os -------------------------------------    674   //  Os -------------------------------------------------------------
701   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 675   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 6, 2,
702   //  Ir -------------------------------------    676   //  Ir -------------------------------------------------------------
703   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 677   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 7, 2,
704   //  Pt -------------------------------------    678   //  Pt -------------------------------------------------------------
705   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 679   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 9, 1,
706   //  Au -------------------------------------    680   //  Au -------------------------------------------------------------
707   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 681   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 1,
708   //  Hg -------------------------------------    682   //  Hg -------------------------------------------------------------
709   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 683   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2,
710   //  Tl -------------------------------------    684   //  Tl -------------------------------------------------------------
711   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 685   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 1,
712   //  Pb -------------------------------------    686   //  Pb -------------------------------------------------------------
713   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, << 687   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
714   //  Bi -------------------------------------    688   //  Bi -------------------------------------------------------------
715   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 689   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 3,
716   //  Po -------------------------------------    690   //  Po -------------------------------------------------------------
717   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6,    691   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 4,
718   //  At -------------------------------------    692   //  At -------------------------------------------------------------
719   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6,    693   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 3,
720   //  Rn -------------------------------------    694   //  Rn -------------------------------------------------------------
721   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6,    695   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4,
722   //  Fr -------------------------------------    696   //  Fr -------------------------------------------------------------
723   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6,    697   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 1,
724   //  Ra -------------------------------------    698   //  Ra -------------------------------------------------------------
725   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6,    699   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 2,
726   //  Ac -------------------------------------    700   //  Ac -------------------------------------------------------------
727   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 701   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 1, 2,
728   //  Th -------------------------------------    702   //  Th -------------------------------------------------------------
729   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6,    703   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 2, 2,
730   //  Pa -------------------------------------    704   //  Pa -------------------------------------------------------------
731   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 705   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 2, 4, 1, 2,
732   //  U  -------------------------------------    706   //  U  -------------------------------------------------------------
733   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 707   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 3, 2, 2, 4, 1, 2,
734   //  Np -------------------------------------    708   //  Np -------------------------------------------------------------
735   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 709   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 4, 2, 2, 4, 1, 2,
736   //  Pu -------------------------------------    710   //  Pu -------------------------------------------------------------
737   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 711   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 6, 2, 2, 4, 2,
738   //  Am -------------------------------------    712   //  Am -------------------------------------------------------------
739   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 713   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 7, 2, 2, 4, 2,
740   //  Cm -------------------------------------    714   //  Cm -------------------------------------------------------------
741   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 715   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 7, 2, 2, 4, 1, 2,
742   //  Bk -------------------------------------    716   //  Bk -------------------------------------------------------------
743   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 717   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 8, 2, 2, 4, 1, 2,
744   //  Cf -------------------------------------    718   //  Cf -------------------------------------------------------------
745   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 719   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6,10, 2, 2, 4, 2,
746   //  Es -------------------------------------    720   //  Es -------------------------------------------------------------
747   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 721   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6,11, 2, 2, 4, 2,
748   //  Fm -------------------------------------    722   //  Fm -------------------------------------------------------------
749   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 723   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6,12, 2, 2, 4, 2
750   //  Md ------------------------------------- << 
751   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 
752   //  No ------------------------------------- << 
753   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 
754   //  Lr ------------------------------------- << 
755   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 
756   //  Rf ------------------------------------- << 
757   2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, << 
758 };                                                724 };
759 // clang-format on                             << 
760                                                   725 
761 //....oooOO0OOooo........oooOO0OOooo........oo    726 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
762                                                   727 
763 G4int G4AtomicShells::GetNumberOfShells(G4int  << 728 G4int
                                                   >> 729 G4AtomicShells::GetNumberOfShells(G4int Z)  
764 {                                                 730 {
765 #ifdef G4VERBOSE                               << 731   assert (Z>=1 && Z<=101);
766   if (Z < 0 || Z > 104) {                      << 
767     Z = PrintErrorZ(Z, "GetNumberOfShells");   << 
768   }                                            << 
769 #endif                                         << 
770   return fNumberOfShells[Z];                      732   return fNumberOfShells[Z];
771 }                                                 733 }
772                                                   734 
773 //....oooOO0OOooo........oooOO0OOooo........oo    735 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
774                                                   736 
775 G4double G4AtomicShells::GetBindingEnergy(G4in << 737 G4double 
776 {                                              << 738 G4AtomicShells::GetBindingEnergy(G4int Z, G4int ShellNb)
777 #ifdef G4VERBOSE                               << 
778   if (Z < 0 || Z > 104) {                      << 
779     Z = PrintErrorZ(Z, "GetBindingEnergy");    << 
780   }                                            << 
781   if (ShellNb < 0 || ShellNb >= fNumberOfShell << 
782     ShellNb = PrintErrorShell(Z, ShellNb, "Get << 
783   }                                            << 
784 #endif                                         << 
785   return fBindingEnergies[fIndexOfShells[Z] +  << 
786 }                                              << 
787                                                << 
788 //....oooOO0OOooo........oooOO0OOooo........oo << 
789                                                << 
790 G4int G4AtomicShells::GetNumberOfElectrons(G4i << 
791 {                                                 739 {
792 #ifdef G4VERBOSE                               << 740   assert (Z>=1 && Z<=101 && ShellNb<fNumberOfShells[Z]);
793   if (Z < 0 || Z > 104) {                      << 741   
794     Z = PrintErrorZ(Z, "GetNumberOfElectrons") << 742   G4int indice = 1;
795   }                                            << 743   for (G4int z = 1 ; z < Z ; z++) indice += fNumberOfShells[z];
796   if (ShellNb < 0 || ShellNb >= fNumberOfShell << 744   indice += ShellNb;  
797     ShellNb = PrintErrorShell(Z, ShellNb, "Get << 745   return fBindingEnergies[indice]*eV;
798   }                                            << 
799 #endif                                         << 
800   return fNumberOfElectrons[fIndexOfShells[Z]  << 
801 }                                                 746 }
802                                                   747 
803 //....oooOO0OOooo........oooOO0OOooo........oo    748 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
804                                                   749 
805 G4int G4AtomicShells::GetNumberOfFreeElectrons << 750 G4int
                                                   >> 751 G4AtomicShells::GetNumberOfElectrons(G4int Z, G4int ShellNb)
806 {                                                 752 {
807 #ifdef G4VERBOSE                               << 753   assert (Z>=1 && Z<=101 && ShellNb<fNumberOfShells[Z]);
808   if (Z < 0 || Z > 104) {                      << 754   
809     Z = PrintErrorZ(Z, "GetNumberOfFreeElectro << 755   G4int indice = 1;
810   }                                            << 756   for (G4int z = 1 ; z < Z ; z++) indice += fNumberOfShells[z];
811 #endif                                         << 757   indice += ShellNb;  
812   G4int idx = fIndexOfShells[Z];               << 758   return fNumberOfElectrons[indice];
813   G4int idxmax = idx + fNumberOfShells[Z];     << 
814   G4int n = 0;                                 << 
815   for (G4int i = idx; i < idxmax; ++i) {       << 
816     if (fBindingEnergies[i] * CLHEP::eV <= th) << 
817       n += fNumberOfElectrons[i];              << 
818     }                                          << 
819   }                                            << 
820   return n;                                    << 
821 }                                                 759 }
822                                                   760 
823 //....oooOO0OOooo........oooOO0OOooo........oo    761 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
824                                                   762 
825 G4double G4AtomicShells::GetTotalBindingEnergy << 763 G4double G4AtomicShells::GetTotalBindingEnergy (G4int Z)
826 {                                                 764 {
827 #ifdef G4VERBOSE                               << 765   assert (Z>=1 && Z<=101);
828   if (Z < 0 || Z > 104) {                      << 766   
829     Z = PrintErrorZ(Z, "GetTotalBindingEnergy" << 767   G4int idx = 1;
830   }                                            << 768   for (G4int z = 1 ; z < Z ; z++) idx += fNumberOfShells[z];
831 #endif                                         << 
832   G4int idx = fIndexOfShells[Z];               << 
833   G4int idxmax = idx + fNumberOfShells[Z];     << 
834   G4double energy = 0.0;                          769   G4double energy = 0.0;
835   for (G4int i = idx; i < idxmax; ++i) {       << 770   G4int idxmax = idx +  fNumberOfShells[Z];
836     energy += fBindingEnergies[i];             << 771   for (G4int i=idx; i<idxmax; i++) {energy += fBindingEnergies[i];}
837   }                                            << 772   return energy*eV;
838   return energy * CLHEP::eV;                   << 
839 }                                                 773 }
840                                                   774 
841 //....oooOO0OOooo........oooOO0OOooo........oo    775 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
842                                                << 
843 #ifdef G4VERBOSE                               << 
844                                                << 
845 G4int G4AtomicShells::PrintErrorZ(G4int Z, con << 
846 {                                              << 
847   G4String sss = "G4AtomicShells::" + ss + "() << 
848   G4ExceptionDescription ed;                   << 
849   ed << "Atomic number out of range Z= " << Z; << 
850   G4Exception(sss, "mat060", FatalException, e << 
851   return 1;                                    << 
852 }                                              << 
853                                                << 
854 //....oooOO0OOooo........oooOO0OOooo........oo << 
855                                                << 
856 G4int G4AtomicShells::PrintErrorShell(G4int Z, << 
857 {                                              << 
858   G4String sss = "G4AtomicShells::" + ss + "() << 
859   G4ExceptionDescription ed;                   << 
860   ed << "Shell number out of range Nshell= " < << 
861   G4Exception(sss, "mat061", FatalException, e << 
862   return 0;                                    << 
863 }                                              << 
864 #endif                                         << 
865                                                   776