G4ParticleHPFFFissionFS Class Reference

#include <G4ParticleHPFFFissionFS.hh>

Inheritance diagram for G4ParticleHPFFFissionFS:

Collaboration diagram for G4ParticleHPFFFissionFS:

Public Member Functions
	G4ParticleHPFFFissionFS ()
	~G4ParticleHPFFFissionFS ()
void	Init (G4double A, G4double Z, G4int M, G4String &dirName, G4String &aFSType, G4ParticleDefinition *)
G4DynamicParticleVector *	ApplyYourself (G4int nNeutrons)
G4ParticleHPFinalState *	New ()
void	GetAFissionFragment (G4double, G4int &, G4int &, G4int &)
Public Member Functions inherited from G4ParticleHPFissionBaseFS
	G4ParticleHPFissionBaseFS ()
virtual	~G4ParticleHPFissionBaseFS ()
void	Init (G4double A, G4double Z, G4int M, G4String &dirName, G4String &bit, G4ParticleDefinition *)
G4DynamicParticleVector *	ApplyYourself (G4int Prompt)
virtual G4double	GetXsec (G4double anEnergy)
virtual G4ParticleHPVector *	GetXsec ()
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 *projectile)
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 (G4double anA, G4double aZ, G4int aM, G4ParticleHPDataUsed used)
void	SetProjectile (G4ParticleDefinition *projectile)

Private Member Functions
G4HadFinalState *	ApplyYourself (const G4HadProjectile &)

Private Attributes
std::map< G4int, std::map< G4double, std::map< G4int, G4double > *> *>	FissionProductYieldData
std::map< G4int, std::map< G4double, G4int > *>	mMTInterpolation

Additional Inherited Members
Protected Member Functions inherited from G4ParticleHPFinalState
void	adjust_final_state (G4LorentzVector)
G4bool	DoNotAdjustFinalState ()
Protected Attributes inherited from G4ParticleHPFinalState
G4bool	hasXsec
G4bool	hasFSData
G4bool	hasAnyData
G4ParticleHPNames	theNames
G4Cache< G4HadFinalState *>	theResult
G4double	theBaseA
G4double	theBaseZ
G4int	theBaseM
G4int	theNDLDataZ
G4int	theNDLDataA
G4int	theNDLDataM
G4ParticleDefinition *	theProjectile

Detailed Description

Definition at line 37 of file G4ParticleHPFFFissionFS.hh.

Constructor & Destructor Documentation

G4ParticleHPFFFissionFS()

G4ParticleHPFFFissionFS::G4ParticleHPFFFissionFS ( )

inline

Definition at line 41 of file G4ParticleHPFFFissionFS.hh.

{ hasXsec = false; }

Here is the call graph for this function:

Here is the caller graph for this function:

~G4ParticleHPFFFissionFS()

G4ParticleHPFFFissionFS::~G4ParticleHPFFFissionFS

(

)

Definition at line 36 of file G4ParticleHPFFFissionFS.cc.

{
    std::map<G4int,std::map<G4double,std::map<G4int,G4double >* >* >::iterator 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 ) {
            std::map<G4double,std::map<G4int,G4double>*>::iterator 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();
    }
    
    std::map< G4int , std::map< G4double , G4int >* >::iterator ii = mMTInterpolation.begin();
    while ( ii != mMTInterpolation.end() ) { // Loop checking, 11.05.2015, T. Koi
        delete ii->second;
        mMTInterpolation.erase(ii);
        ii = mMTInterpolation.begin();
    }
}

Here is the caller graph for this function:

Member Function Documentation

ApplyYourself() [1/2]

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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

ApplyYourself() [2/2]

G4HadFinalState* G4ParticleHPFFFissionFS::ApplyYourself ( const G4HadProjectile &  )

inlineprivatevirtual

Reimplemented from G4ParticleHPFinalState.

Definition at line 58 of file G4ParticleHPFFFissionFS.hh.

{ return NULL; }

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;

   std::map< G4double , std::map< G4int , G4double >* >* mEnergyFSPData = FissionProductYieldData.find( 454 )->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->begin()->first;
   }
   else
   {
      //Find closest energy point
      G4double Dmin=DBL_MAX; 
      G4int i = 0;
      for ( std::map< G4double , std::map< G4int , G4double >* >::iterator it = mEnergyFSPData->begin() ; 
      it != mEnergyFSPData->end() ; it++ )
      {
         G4double e = (it->first);
         G4double d = std::fabs ( energy - e ); 
         if ( d < Dmin ) 
         {
            Dmin = d;
            key_energy = e;
         }
         i++;
      }
   }

   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 ( std::map<G4int,G4double>::iterator it = mFSPYieldData->begin() ; it != mFSPYieldData->end() ; 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;
}

Here is the caller graph for this function:

Init()

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

virtual

Implements G4ParticleHPFinalState.

Definition at line 64 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
      std::map< G4double , std::map< G4int , G4double >* >* mEnergyFSPData = new std::map< G4double , std::map< G4int , G4double >* >;

      std::map< G4double , G4int >* 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
         std::map<G4int,G4double>* 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();
}

Here is the call graph for this function:

Here is the caller graph for this function:

New()

G4ParticleHPFinalState* G4ParticleHPFFFissionFS::New ( )

inlinevirtual

Implements G4ParticleHPFinalState.

Definition at line 48 of file G4ParticleHPFFFissionFS.hh.

      {
         G4ParticleHPFFFissionFS * theNew = new G4ParticleHPFFFissionFS;
         return theNew;
      }

Here is the call graph for this function:

Member Data Documentation

FissionProductYieldData

std::map< G4int , std::map< G4double , std::map< G4int , G4double >* >* > G4ParticleHPFFFissionFS::FissionProductYieldData

private

Definition at line 61 of file G4ParticleHPFFFissionFS.hh.

mMTInterpolation

std::map< G4int , std::map< G4double , G4int >* > G4ParticleHPFFFissionFS::mMTInterpolation

private

Definition at line 62 of file G4ParticleHPFFFissionFS.hh.

---

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

• G4ParticleHPFFFissionFS.hh
• G4ParticleHPFFFissionFS.cc