100 if ( !getenv(
"G4DPMJET2_5DATA") )
102 G4cout <<
"ENVIRONMENT VARIABLE G4DPMJET2_5DATA NOT SET " <<
G4endl;
104 "Please setenv G4DPMJET2_5DATA to point to the dpmjet2.5 data files.");
151 std::ostringstream os;
159 std::ostringstream os;
194 std::ifstream glauberFile(filename);
205 G4int asci = glauberFile.peek();
206 if (asci >= 48 && asci <= 57) i = asci - 48;
207 else if (asci >= 65 && asci <= 90) i = asci - 55;
209 G4cerr <<
"ERROR IN G4GlaubAADataSetHandler::LoadGlauberDataReturnPtr"
212 G4cerr <<
"IDENTIFIED AS GLAUBER FILE TYPE " <<i <<
G4endl;
218 glauberFile >>(*glauberData);
223 G4int AP1 = glauberData->GetAP();
224 G4int AT1 = glauberData->GetAT();
225 G4bool found = AP1==AP && AT1==AT;
237 theIndex.insert (G4GlaubAADataSetIndex::value_type(n,glauberData));
241 G4cout <<
"****************************************"
242 <<
"****************************************"
244 G4cout <<
"In G4GlaubAADataSetHandler::LoadGlauberDataReturnPtr"
246 G4cout <<
"LOADED GLAUBER DATA SET FOR PROJECTILE A = " <<AP
247 <<
" & TARGET A = " <<AT <<
G4endl;
249 G4cout <<
"****************************************"
250 <<
"****************************************"
257 G4cout <<
"****************************************"
258 <<
"****************************************"
260 G4cout <<
"In G4GlaubAADataSetHandler::LoadGlauberDataReturnPtr"
262 G4cout <<
"LOADED GLAUBER DATA SET FOR PROJECTILE A = " <<AP
263 <<
" & TARGET A = " <<AT <<
" TEMPORARILY" <<
G4endl;
264 G4cout <<
"****************************************"
265 <<
"****************************************"
271 G4cerr <<
"WARNING IN G4GlaubAADataSetHandler::LoadGlauberData"
273 G4cerr <<
"GLAUBER FILE " <<filename <<
" LOOKED LIKE IT SHOULD CONTAIN"
275 G4cerr <<
"DATA FOR AP = " <<AP <<
" AND AT = " <<AT <<
G4endl;
276 G4cerr <<
"BUT CONTAINED AP = " <<AP1 <<
" AND AT = " <<AT1 <<
G4endl;
300 for (G4GlaubAADataSetIndex::iterator it=
theIndex.begin(); it!=
theIndex.end();
315 const G4int AT)
const
318 G4GlaubAADataSetIndex::const_iterator it =
theIndex.find(glauberID);
323 std::ifstream glauberFile;
324 glauberFile.open(filename);
348 G4cout <<
"WARNING: G4GlaubAADataSetHandler::CheckIfSpace:"
350 G4cout <<
"MAXIMUM NUMBER OF GLAUBER DATASETS REACHED : "
367 G4cerr <<
"ERROR IN G4GlaubAADataSetHandler::SetMaxGlauberDataSets"
369 G4cerr <<
"MAXIMUM NUMBER OF GLAUBER DATA SETS WOULD BE EXCEEDED IF YOU"
371 G4cerr <<
"SET THIS VALUE. ATTEMPTED TO SET TO " <<n
396 G4GlaubAADataSetIndex::iterator it =
theIndex.find(glauberID);
454 if (ppn1*
GeV > 1.0*
eV) ppn = ppn1;
467 if (ppn1*
GeV > 1.0*
eV) ppn = ppn1;
G4bool IsGlauberDataSetAvailable(const G4int AP, const G4int AT) const
G4GlaubAADataSetHandler()
void SetMaxGlauberDataSets(const G4int n)
G4int GetIndexID(const G4int AP, const G4int AT) const
Definition of the G4ParamType1GlaubAADataSet class.
G4GlaubAADataSet * theCurrentGlauberDataSet
void ResetCurrentGlauberDataSet()
G4String GetStringID(const G4int AP, const G4int AT) const
G4bool usingLocalGlauberDataSet
G4String glauberDataSetDir
G4GlaubAADataSet * GetCurrentGlauberDataSet() const
Definition of the G4FullGlaubAADataSet class.
virtual ~G4GlaubAADataSetHandler()
struct ccdpm25dtumat dtumat_
G4String GetProjectileStringID(const G4int AP) const
G4GLOB_DLL std::ostream G4cout
G4VGlauberDataSet * LoadGlauberDataReturnPtr(const G4int AP, const G4int AT)
static G4GlaubAADataSetHandler * instance
G4double GetValueN(const G4double v, const G4double ppn1=0.0)
G4double GetValueM(const G4double v, const G4double ppn1=0.0)
G4int UnloadAllGlauberData()
G4String GetTargetStringID(const G4int AT) const
void SetVerboseLevel(const G4int i)
G4bool SetCurrentGlauberDataSet(const G4int AP, const G4int AT, const G4double ppn=0.0)
Definition of the G4DPMJET2_5Interface class.
G4bool CheckIfSpace() const
Definition of the G4GlaubAADataSetHandler class.
static G4GlaubAADataSetHandler * getInstance()
G4GlaubAADataSetIndex theIndex
G4GLOB_DLL std::ostream G4cerr