G4ParticleHPFFFissionFS Class Reference

#include <G4ParticleHPFFFissionFS.hh>

Inheritance diagram for G4ParticleHPFFFissionFS:

Public Member Functions
	G4ParticleHPFFFissionFS ()
	~G4ParticleHPFFFissionFS () override
void	Init (G4double A, G4double Z, G4int M, G4String &dirName, G4String &aFSType, G4ParticleDefinition *) override
G4DynamicParticleVector *	ApplyYourself (G4int nNeutrons)
G4ParticleHPFinalState *	New () override
void	GetAFissionFragment (G4double, G4int &, G4int &, G4int &)
Public Member Functions inherited from G4ParticleHPFissionBaseFS
	G4ParticleHPFissionBaseFS ()
	~G4ParticleHPFissionBaseFS () override
void	Init (G4double A, G4double Z, G4int M, G4String &dirName, G4String &bit, G4ParticleDefinition *) override
G4DynamicParticleVector *	ApplyYourself (G4int Prompt)
G4double	GetXsec (G4double anEnergy) override
G4ParticleHPVector *	GetXsec () override
void	SetNeutronRP (const G4ReactionProduct &aNeutron)
void	SetTarget (const G4ReactionProduct &aTarget)
Public Member Functions inherited from G4ParticleHPFinalState
	G4ParticleHPFinalState ()
virtual	~G4ParticleHPFinalState ()
void	Init (G4double A, G4double Z, G4String &dirName, G4String &aFSType, G4ParticleDefinition *p)
G4bool	HasXsec ()
G4bool	HasFSData ()
G4bool	HasAnyData ()
void	SetA_Z (G4double anA, G4double aZ, G4int aM=0)
G4double	GetZ ()
G4double	GetN ()
G4double	GetA ()
G4int	GetM ()
void	SetAZMs (G4ParticleHPDataUsed used)
void	SetAZMs (G4double anA, G4double aZ, G4int aM, G4ParticleHPDataUsed used)
void	SetProjectile (G4ParticleDefinition *projectile)
G4ParticleHPFinalState &	operator= (const G4ParticleHPFinalState &right)=delete
	G4ParticleHPFinalState (const G4ParticleHPFinalState &)=delete

Additional Inherited Members
Protected Member Functions inherited from G4ParticleHPFinalState
void	adjust_final_state (G4LorentzVector)
Protected Attributes inherited from G4ParticleHPFinalState
G4ParticleDefinition *	theProjectile {nullptr}
G4ParticleHPManager *	fManager
G4IonTable *	ionTable
G4int	theBaseA {0}
G4int	theBaseZ {0}
G4int	theBaseM {0}
G4int	theNDLDataZ {0}
G4int	theNDLDataA {0}
G4int	theNDLDataM {0}
G4int	secID {-1}
G4bool	hasXsec {true}
G4bool	hasFSData {true}
G4bool	hasAnyData {true}
G4ParticleHPNames	theNames
G4Cache< G4HadFinalState * >	theResult

Detailed Description

Definition at line 37 of file G4ParticleHPFFFissionFS.hh.

Constructor & Destructor Documentation

G4ParticleHPFFFissionFS()

G4ParticleHPFFFissionFS::G4ParticleHPFFFissionFS ( )

inline

Definition at line 40 of file G4ParticleHPFFFissionFS.hh.

{ hasXsec = false; }

Referenced by New().

~G4ParticleHPFFFissionFS()

G4ParticleHPFFFissionFS::~G4ParticleHPFFFissionFS ( )

override

Definition at line 37 of file G4ParticleHPFFFissionFS.cc.

{
  auto it = FissionProductYieldData.begin();
  while (it != FissionProductYieldData.end()) {  // Loop checking, 11.05.2015, T. Koi
    std::map<G4double, std::map<G4int, G4double>*>* firstLevel = it->second;
    if (firstLevel != nullptr) {
      auto it2 = firstLevel->begin();
      while (it2 != firstLevel->end()) {  // Loop checking, 11.05.2015, T. Koi
        delete it2->second;
        it2->second = 0;
        firstLevel->erase(it2);
        it2 = firstLevel->begin();
      }
    }
    delete firstLevel;
    it->second = 0;
    FissionProductYieldData.erase(it);
    it = FissionProductYieldData.begin();
  }
 
  auto ii = mMTInterpolation.begin();
  while (ii != mMTInterpolation.end()) {  // Loop checking, 11.05.2015, T. Koi
    delete ii->second;
    mMTInterpolation.erase(ii);
    ii = mMTInterpolation.begin();
  }
}

Member Function Documentation

ApplyYourself()

G4DynamicParticleVector * G4ParticleHPFFFissionFS::ApplyYourself

(

G4int

nNeutrons

)

Definition at line 151 of file G4ParticleHPFFFissionFS.cc.

{
  G4DynamicParticleVector* aResult;
  //    G4cout <<"G4ParticleHPFFFissionFS::ApplyYourself +"<<G4endl;
  aResult = G4ParticleHPFissionBaseFS::ApplyYourself(nNeutrons);
  return aResult;
}

GetAFissionFragment()

void G4ParticleHPFFFissionFS::GetAFissionFragment	(	G4double	energy,
		G4int &	fragZ,
		G4int &	fragA,
		G4int &	fragM )

Definition at line 159 of file G4ParticleHPFFFissionFS.cc.

{
  // G4cout << "G4ParticleHPFFFissionFS::GetAFissionFragment " << G4endl;
 
  G4double rand = G4UniformRand();
  // G4cout << rand << G4endl;
 
  auto ptr = FissionProductYieldData.find(454);
  if (ptr == FissionProductYieldData.end())
    return;
 
  auto mEnergyFSPData = ptr->second;
 
  // It is not clear that the treatment of the scheme 2 on two-dimensional interpolation.
  // So, here just use the closest energy point array of yield data.
  // TK120531
  G4double key_energy = DBL_MAX;
  if (mEnergyFSPData->size() == 1) {
    key_energy = mEnergyFSPData->cbegin()->first;
  }
  else {
    // Find closest energy point
    G4double Dmin = DBL_MAX;
    for (auto it = mEnergyFSPData->cbegin(); it != mEnergyFSPData->cend(); ++it) {
      G4double e = (it->first);
      G4double d = std::fabs(energy - e);
      if (d < Dmin) {
        Dmin = d;
        key_energy = e;
      }
    }
  }
 
  std::map<G4int, G4double>* mFSPYieldData = (*mEnergyFSPData)[key_energy];
 
  G4int ifrag = 0;
  G4double ceilling =
    mFSPYieldData->rbegin()->second;  // Because of numerical accuracy, this is not always 2
  for (auto it = mFSPYieldData->cbegin(); it != mFSPYieldData->cend(); ++it) {
    // if ( ( rand - it->second/ceilling ) < 1.0e-6 ) std::cout << rand - it->second/ceilling <<
    // std::endl;
    if (rand <= it->second / ceilling) {
      // G4cout << it->first << " " << it->second/ceilling << G4endl;
      ifrag = it->first;
      break;
    }
  }
 
  fragZ = ifrag / 100000;
  fragA = (ifrag % 100000) / 100;
  fragM = (ifrag % 100);
 
  // G4cout << fragZ << " " << fragA << " " << fragM << G4endl;
}

Referenced by G4ParticleHPFissionFS::ApplyYourself().

Init()

void G4ParticleHPFFFissionFS::Init ( G4double A, G4double Z, G4int M, G4String & dirName, G4String & aFSType, G4ParticleDefinition *  )

overridevirtual

Implements G4ParticleHPFinalState.

Definition at line 65 of file G4ParticleHPFFFissionFS.cc.

{
  // G4cout << "G4ParticleHPFFFissionFS::Init" << G4endl;
  G4String aString = "FF";
 
  G4String tString = dirName;
  G4bool dbool;
  G4ParticleHPDataUsed aFile =
    theNames.GetName(static_cast<G4int>(A), static_cast<G4int>(Z), M, tString, aString, dbool);
  G4String filename = aFile.GetName();
  theBaseA = aFile.GetA();
  theBaseZ = aFile.GetZ();
 
  // 3456
  if (!dbool || (Z < 2.5 && (std::abs(theBaseZ - Z) > 0.0001 || std::abs(theBaseA - A) > 0.0001))) {
    hasAnyData = false;
    hasFSData = false;
    hasXsec = false;
    return;  // no data for exactly this isotope.
  }
  // std::ifstream theData(filename, std::ios::in);
  std::istringstream theData(std::ios::in);
  G4ParticleHPManager::GetInstance()->GetDataStream(filename, theData);
  G4double dummy;
  if (!theData) {
    // theData.close();
    hasFSData = false;
    hasXsec = false;
    hasAnyData = false;
    return;  // no data for this FS for this isotope
  }
 
  hasFSData = true;
  //          MT              Energy            FPS    Yield
  // std::map< int , std::map< double , std::map< int , double >* >* > FisionProductYieldData;
  while (theData.good())  // Loop checking, 11.05.2015, T. Koi
  {
    G4int iMT, iMF;
    G4int imax;
    // Reading the data
    //          MT       MF       AWR
    theData >> iMT >> iMF >> dummy;
    //         nBlock
    theData >> imax;
    // if ( !theData.good() ) continue;
    //         Ei                   FPS     Yield
    auto mEnergyFSPData = new std::map<G4double, std::map<G4int, G4double>*>;
 
    auto mInterporation = new std::map<G4double, G4int>;
    for (G4int i = 0; i <= imax; i++) {
      G4double YY = 0.0;
      G4double Ei;
      G4int jmax;
      G4int ip;
      //        energy of incidence neutron
      theData >> Ei;
      //        Number of data set followings
      theData >> jmax;
      //         interpolation scheme
      theData >> ip;
      mInterporation->insert(std::pair<G4double, G4int>(Ei * eV, ip));
      //         nNumber  nIP
      auto mFSPYieldData = new std::map<G4int, G4double>;
      for (G4int j = 0; j < jmax; j++) {
        G4int FSP;
        G4int mFSP;
        G4double Y;
        theData >> FSP >> mFSP >> Y;
        G4int k = FSP * 100 + mFSP;
        YY = YY + Y;
        // if ( iMT == 454 )G4cout << iMT << " " << i << " " << j << " " <<  k << " " << Y << " " <<
        // YY << G4endl;
        mFSPYieldData->insert(std::pair<G4int, G4double>(k, YY));
      }
      mEnergyFSPData->insert(
        std::pair<G4double, std::map<G4int, G4double>*>(Ei * eV, mFSPYieldData));
    }
 
    FissionProductYieldData.insert(
      std::pair<G4int, std::map<G4double, std::map<G4int, G4double>*>*>(iMT, mEnergyFSPData));
    mMTInterpolation.insert(std::pair<G4int, std::map<G4double, G4int>*>(iMT, mInterporation));
  }
  // theData.close();
}

Referenced by G4ParticleHPFissionFS::Init().

New()

G4ParticleHPFinalState * G4ParticleHPFFFissionFS::New ( )

inlineoverridevirtual

Implements G4ParticleHPFinalState.

Definition at line 48 of file G4ParticleHPFFFissionFS.hh.

    {
      auto theNew = new G4ParticleHPFFFissionFS;
      return theNew;
    }

---

The documentation for this class was generated from the following files:

• G4ParticleHPFFFissionFS.hh
• G4ParticleHPFFFissionFS.cc