88 std::vector<G4VEMDataSet *>::const_iterator i(
components.begin());
89 std::vector<G4VEMDataSet *>::const_iterator end(
components.end());
93 value += (*i)->FindValue(energy);
105 G4cout <<
"The data set has " << n <<
" components" <<
G4endl;
164 std::ifstream
in(fullFileName);
169 message += fullFileName;
170 message +=
"\" not found";
181 G4int shellIndex = 0;
199 if ((k%nColumns == 0) && (orig_shell_energies != 0) )
202 orig_shell_energies = 0;
204 log_shell_energies = 0;
210 if (orig_shell_energies == 0)
220 log_shell_energies->push_back(std::log10(a) + std::log10(unitEnergies));
222 else if (k%nColumns == 1)
224 orig_shell_data->push_back(a*
unitData);
225 log_shell_data->push_back(std::log10(a) + std::log10(unitData));
234 delete orig_shell_energies;
235 delete orig_shell_data;
236 delete log_shell_energies;
237 delete log_shell_data;
248 std::ifstream
in(fullFileName);
253 message += fullFileName;
254 message +=
"\" not found";
263 G4int shellIndex = 0;
279 if ((k%nColumns == 0) && (orig_shell_energies != 0) )
282 orig_shell_energies = 0;
288 if (orig_shell_energies == 0)
297 else if (k%nColumns == 1)
299 orig_shell_data->push_back(a*
unitData);
308 delete orig_shell_energies;
309 delete orig_shell_data;
319 std::ofstream out(fullFileName);
324 message += fullFileName;
341 G4DataVector::const_iterator i = energies.begin();
342 G4DataVector::const_iterator endI = energies.end();
343 G4DataVector::const_iterator j = data.begin();
355 out << ((*j)/
unitData) << std::endl;
369 out << -1.f << std::endl;
382 out << -2.f << std::endl;
400 char* path = std::getenv(
"G4LEDATA");
408 std::ostringstream fullFileName;
410 fullFileName << path <<
'/' << fileName <<
z <<
".dat";
412 return G4String(fullFileName.str().c_str());