77 :
fName(name), fSymbol(symbol)
82 ed <<
"Fail to create G4Element " << name
83 <<
" Z= " << zeff <<
" < 1 !";
88 ed <<
"G4Element Warning: " << name <<
" Z= " << zeff
89 <<
" A= " << aeff/(
g/
mole);
103 ed <<
"Fail to create G4Element " << name
104 <<
" with Z= " << zeff <<
" N= " <<
fNeff
105 <<
" N < Z is not allowed" <<
G4endl;
130 :
fName(name),fSymbol(symbol)
134 size_t n = size_t(nIsotopes);
138 ed <<
"Fail to create G4Element " << name
139 <<
" <" << symbol <<
"> with " << nIsotopes
156 ed <<
"Fail to add Isotope to G4Element " <<
fName
169 ed <<
"Fail to add Isotope Z= " << iz <<
" to G4Element " <<
fName
181 ed <<
"Fail to add Isotope Z= " << iz <<
" to G4Element " <<
fName
182 <<
" - more isotopes than declaired ";
195 if(wtSum > 0.0) {
fAeff /= wtSum; }
246 : fZeff(0), fNeff(0), fAeff(0)
292 static const G4double k1 = 0.0083 , k2 = 0.20206 ,k3 = 0.0020 , k4 = 0.0369 ;
297 fCoulomb = (k1*az4 + k2 + 1./(1.+az2))*az2 - (k3*az4 + k4)*az4;
308 static const G4double Lrad_light[] = {5.31 , 4.79 , 4.74 , 4.71} ;
309 static const G4double Lprad_light[] = {6.144 , 5.621 , 5.805 , 5.924} ;
317 if (iz <= 3) { Lrad = Lrad_light[iz] ; Lprad = Lprad_light[iz] ; }
318 else { Lrad = log184 - logZ3 ; Lprad = log1194 - 2*logZ3;}
333 const std::vector<G4String> elmnames =
340 for(
G4int i=0; i<
n; ++i) {
347 for(
G4int i=0; i<
n; ++i) {
351 std::ostringstream strm;
353 (*theIsotopeVector)[idx] =
new G4Isotope(strm.str(),
Z,
N, 0.0, 0);
359 if(xsum != 0.0 && xsum != 1.0) {
371 ed <<
"Invalid argument " << i <<
" in for G4Element " <<
fName
372 <<
" with Z= " <<
fZeff
386 ed <<
"Invalid argument " << i <<
" for G4Element " <<
fName
387 <<
" with Z= " <<
fZeff
389 G4Exception(
"G4Element::GetNbOfShellElectrons()",
"mat016",
423 G4cout <<
"\n---> warning from G4Element::GetElement(). The element: "
424 << elementName <<
" does not exist in the table. Return NULL pointer."
434 std::ios::fmtflags mode = flux.flags();
435 flux.setf(std::ios::fixed,std::ios::floatfield);
439 <<
" Element: " << element->
fName <<
" (" << element->
fSymbol <<
")"
440 <<
" Z = " << std::setw(4) << std::setprecision(1) << element->
fZeff
441 <<
" N = " << std::setw(5) << std::setprecision(1)
443 <<
" A = " << std::setw(6) << std::setprecision(3)
449 <<
" abundance: " << std::setw(6) << std::setprecision(3)
452 flux.precision(prec);
453 flux.setf(mode,std::ios::floatfield);
470 flux <<
"\n***** Table : Nb of elements = " << ElementTable.size()
473 for (
size_t i=0; i<ElementTable.size(); i++) flux << ElementTable[i]
G4double GetAtomicShell(G4int index) const
void ComputeCoulombFactor()
void InitializePointers()
std::ostringstream G4ExceptionDescription
G4IonisParamElm * fIonisation
std::ostream & operator<<(std::ostream &, const BasicVector3D< float > &)
static constexpr double alpha_rcl2
const G4String & GetName() const
static G4Element * GetElement(G4String name, G4bool warning=true)
static constexpr double perCent
G4int GetNumberOfNistIsotopes(G4int Z) const
static constexpr double perMillion
G4double G4Log(G4double x)
G4double * fRelativeAbundanceVector
G4Element(const G4String &name, const G4String &symbol, G4double Zeff, G4double Aeff)
G4IsotopeVector * theIsotopeVector
G4int * fNbOfShellElectrons
void AddIsotope(G4Isotope *isotope, G4double RelativeAbundance)
static constexpr double g
const std::vector< G4String > & GetNistElementNames() const
std::vector< G4Isotope * > G4IsotopeVector
G4double GetIsotopeAbundance(G4int Z, G4int N) const
std::vector< G4Element * > G4ElementTable
void ComputeLradTsaiFactor()
static G4int GetNumberOfShells(G4int Z)
static G4int GetNumberOfElectrons(G4int Z, G4int SubshellNb)
void AddNaturalIsotopes()
static G4ElementTable theElementTable
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
G4int GetNbOfShellElectrons(G4int index) const
static G4double GetBindingEnergy(G4int Z, G4int SubshellNb)
G4GLOB_DLL std::ostream G4cout
void ComputeDerivedQuantities()
static size_t GetNumberOfElements()
static constexpr double mole
static constexpr double fine_structure_const
static G4NistManager * Instance()
static G4ElementTable * GetElementTable()
G4int GetNistFirstIsotopeN(G4int Z) const