53 Verbosity_(G4FFGDefaultValues::
Verbosity),
54 YieldType_(WhichYield)
67 Verbosity_(Verbosity),
68 YieldType_(WhichYield)
80 Verbosity_(Verbosity),
81 YieldType_(WhichYield)
109 }
catch (std::exception &
e)
146 return NumberOfElements;
155 if(WhichYield >= 0 && WhichYield < YieldContainerTable_->G4GetNumberOfElements())
180 if(!dataStream.good())
185 "Fission product data not available");
189 throw std::exception();
410 G4int currentEnergy = 0;
423 std::vector< G4double > projectileEnergies;
424 std::map< const G4int, std::pair< std::vector< G4double >, std::vector< G4double > > > intermediateData;
425 std::map< const G4int, std::pair< std::vector< G4double >, std::vector< G4double > > >::iterator dataIterator;
427 while(dataStream.good())
429 dataStream >> MT >> MF >> dummy >> blockCount;
433 for(
G4int b = 0; b < blockCount; ++b)
435 dataStream >> incidentEnergy >> itemCount >> interpolation;
436 maxSize = maxSize >= itemCount ? maxSize : itemCount;
441 projectileEnergies.push_back(incidentEnergy);
442 currentEnergy = projectileEnergies.size() - 1;
449 for(
G4int i = 0; i < itemCount; ++i)
451 dataStream >> isotope >> metastate >> yield;
455 identifier = isotope * 10 + metastate;
457 dataIterator = intermediateData.insert(std::make_pair(
460 std::vector< G4double >(projectileEnergies.size(), 0.0),
461 std::vector< G4double >(projectileEnergies.size(), 0.0)))).
first;
463 if(dataIterator->second.first.size() < projectileEnergies.size())
465 dataIterator->second.first.resize(projectileEnergies.size());
466 dataIterator->second.second.resize(projectileEnergies.size());
469 dataIterator->second.first[currentEnergy] = yield;
470 dataIterator->second.second[currentEnergy] =
error;
495 for(dataIterator = intermediateData.begin(); dataIterator != intermediateData.end(); ++dataIterator)
497 identifier = dataIterator->first;
498 metastate = identifier % 10;
513 "Unsupported metastable state supplied in fission yield data. Defaulting to the ground state");
519 NewProduct = (identifier - metastate) / 10;
523 if(energyGroup < (
signed)dataIterator->second.first.size())
525 yield = dataIterator->second.first[energyGroup];
526 error = dataIterator->second.second[energyGroup];
533 NewYield[energyGroup] = yield;
534 NewError[energyGroup] =
error;
static G4ParticleHPManager * GetInstance()
void ReadInData(std::istringstream &dataStream)
G4ENDFTapeRead(G4String FileLocation, G4String FileName, G4FFGEnumerations::YieldType WhichYield, G4FFGEnumerations::FissionCause WhichCause)
G4ENDFYieldDataContainer * G4GetYield(G4int WhichYield)
#define G4FFG_FUNCTIONLEAVE__
#define G4FFG_DATA_FUNCTIONENTER__
T * G4GetNewContainer(void)
G4long G4GetNumberOfElements(void)
void SetYieldError(G4double *YieldError)
T * G4GetContainer(unsigned int WhichContainer)
void SetMetaState(G4FFGEnumerations::MetaState MetaState)
void GetDataStream(G4String, std::istringstream &iss)
void SetProduct(G4int Product)
void SetYieldProbability(G4double *YieldProbability)
void G4SetVerbosity(G4int WhatVerbosity)
G4int G4GetNumberOfEnergyGroups(void)
G4int G4GetNumberOfFissionProducts(void)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
G4TableTemplate< G4ENDFYieldDataContainer > * YieldContainerTable_
void Initialize(G4String dataFile)
G4double * G4GetEnergyGroupValues(void)
#define G4FFG_DATA_FUNCTIONLEAVE__
G4double * EnergyGroupValues_
static PROLOG_HANDLER error
#define G4FFG_FUNCTIONENTER__
const G4FFGEnumerations::YieldType YieldType_