52 "1-",
"1+",
"2-",
"2+",
"3-",
"3+",
"4-",
"4+",
"5-",
"5+"};
55 : fMinProbability(1.e-8),fVerbose(0)
58 char* directory = getenv(
"G4LEVELGAMMADATA");
63 "Environment variable G4LEVELGAMMADATA is not defined");
68 for(
G4int i=0; i<10; ++i) {
fICC[i] = 0.0; }
97 std::ostringstream ss;
98 ss <<
"/z" << Z <<
".a" <<
A;
119 std::ifstream infile(filename, std::ios::in);
122 if (!infile.is_open()) {
124 G4cout <<
" G4LevelReader: fail open file for Z= "
126 <<
" <" << filename <<
">" <<
G4endl;
132 G4cout <<
"G4LevelReader: open file for Z= "
134 <<
" <" << filename <<
">" <<
G4endl;
157 for(
G4int i=0; i<10; ++i) {
166 G4cout <<
"Line #" << nline <<
" in file <" << filename
167 <<
" is corrupted Z= " << Z <<
" A= " << A <<
G4endl;
180 <<
" Ntransitions= " << nn <<
" fNorm1= " <<
fNorm1
189 for(
size_t i=0; i<
nn; ++i) {
220 if(!end) { next =
true; }
250 G4float ediffMin = std::abs(efinal);
259 if(ediff < ediffMin) {
273 for(; tnum<10; ++tnum) {
if(
fTrans[tnum] ==
fPol) {
break; } }
284 G4cout <<
"G4LevelReader: Line #" << nline <<
" Z= " << Z <<
" A= "
285 <<
" this file is too long - stop loading" <<
G4endl;
297 G4cout <<
"Reader: new manager for Z= " << Z <<
" A= " << A
298 <<
" Nlevels= " <<
vEnergy.size() <<
" E[0]= "
313 if(dataFile.fail()) { okay =
false; }
314 else { x = strtod(buffer,
nullptr); }
323 if(dataFile.fail()) { okay =
false; }
330 std::vector<G4float>* vec =
nullptr;
357 if(LL < 3) {
for(
G4int i=LL+1; i<=4; ++i) {
fICC[i] = 0.0; } }
358 if(M < 5) {
for(
G4int i=M+4; i<=8; ++i) {
fICC[i] = 0.0; } }
359 if(N < 1) {
fICC[9] = 0.0; }
360 for(
G4int i=0; i<10; ++i) {
366 vec =
new std::vector<G4float>;
368 for(
G4int i=0; i<9; ++i) {
371 vec->push_back(1.0f);
static G4Pow * GetInstance()
std::vector< G4int > vSpin
std::vector< const G4NucLevel * > vLevel
const std::vector< G4float > * NormalizedICCProbability(G4int Z)
std::vector< G4float > vTransEnergy
const G4LevelManager * MakeLevelManager(G4int Z, G4int A, const G4String &filename)
G4double logZ(G4int Z) const
std::vector< G4float > vTime
G4GLOB_DLL std::ostream G4cout
double A(double temperature)
G4bool ReadDataItem(std::istream &dataFile, G4double &x)
std::vector< G4float > vTimeg
static const double second
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
static G4String fTrans[10]
std::vector< const std::vector< G4float > * > vShellProbability
T max(const T t1, const T t2)
brief Return the largest of the two arguments
std::vector< G4int > vTrans
const G4double x[NPOINTSGL]
static const G4int LL[nN]
std::vector< G4float > vGammaCumProbability
const G4LevelManager * CreateLevelManager(G4int Z, G4int A)
std::vector< G4float > vGammaECumProbability
G4float MaxLevelEnergy() const
std::vector< G4float > vEnergy
std::vector< G4float > vGammaProbability