G4ParticleHPJENDLHEData Class Reference

#include <G4ParticleHPJENDLHEData.hh>

Inheritance diagram for G4ParticleHPJENDLHEData:

Collaboration diagram for G4ParticleHPJENDLHEData:

Public Member Functions
	G4ParticleHPJENDLHEData ()
	G4ParticleHPJENDLHEData (G4String, G4ParticleDefinition *)
	~G4ParticleHPJENDLHEData ()
G4bool	IsApplicable (const G4DynamicParticle *, const G4Element *)
G4bool	IsZAApplicable (const G4DynamicParticle *, G4double, G4double)
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, G4double aT)
void	BuildPhysicsTable (const G4ParticleDefinition &)
void	DumpPhysicsTable (const G4ParticleDefinition &)
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
virtual G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=0, const G4Material *mat=0)
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=0)
G4double	ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=0)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
virtual G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy)
virtual void	CrossSectionDescription (std::ostream &) const
virtual G4int	GetVerboseLevel () const
virtual void	SetVerboseLevel (G4int value)
G4double	GetMinKinEnergy () const
void	SetMinKinEnergy (G4double value)
G4double	GetMaxKinEnergy () const
void	SetMaxKinEnergy (G4double value)
const G4String &	GetName () const

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Detailed Description

Definition at line 48 of file G4ParticleHPJENDLHEData.hh.

Constructor & Destructor Documentation

G4ParticleHPJENDLHEData::G4ParticleHPJENDLHEData

(

)

Definition at line 60 of file G4ParticleHPJENDLHEData.cc.

{
  for ( std::map< G4int , std::map< G4int , G4PhysicsVector* >* >::iterator itZ = mIsotope.begin();
        itZ != mIsotope.end(); ++itZ ) {
    std::map< G4int , G4PhysicsVector* >* pointer_map = itZ->second;
    if ( pointer_map ) {
      for ( std::map< G4int , G4PhysicsVector* >::iterator itA = pointer_map->begin();
            itA != pointer_map->end() ; ++itA ) {
        G4PhysicsVector* pointerPhysicsVector = itA->second;
        if ( pointerPhysicsVector ) {
          delete pointerPhysicsVector;
          itA->second = NULL;
        }
      }
      delete pointer_map;
      itZ->second = NULL;
    }
  }
  mIsotope.clear();
}

G4ParticleHPJENDLHEData::G4ParticleHPJENDLHEData	(	G4String	reaction,
		G4ParticleDefinition *	pd
	)

Definition at line 83 of file G4ParticleHPJENDLHEData.cc.

:G4VCrossSectionDataSet( "JENDLHE"+reaction+"CrossSection" )
{
   reactionName = reaction;
   BuildPhysicsTable( *pd );
}

Here is the call graph for this function:

G4ParticleHPJENDLHEData::~G4ParticleHPJENDLHEData

(

)

Definition at line 92 of file G4ParticleHPJENDLHEData.cc.

{
   ; 
   //delete theCrossSections;
}

Member Function Documentation

void G4ParticleHPJENDLHEData::BuildPhysicsTable ( const G4ParticleDefinition &  aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 100 of file G4ParticleHPJENDLHEData.cc.

{

//   if ( &aP != G4Neutron::Neutron() ) 
//      throw G4HadronicException(__FILE__, __LINE__, "Attempt to use NeutronHP data for particles other than neutrons!!!");  
   particleName = aP.GetParticleName();

   G4String baseName = getenv( "G4NEUTRONHPDATA" );
   G4String dirName = baseName+"/JENDL_HE/"+particleName+"/"+reactionName ;
   G4String aFSType = "/CrossSection/";
   G4ParticleHPNames theNames; 

   G4String filename;

// Create JENDL_HE data 
// Create map element or isotope  

   size_t numberOfElements = G4Element::GetNumberOfElements();
   //theCrossSections = new G4PhysicsTable( numberOfElements );

   // make a PhysicsVector for each element

   static G4ThreadLocal G4ElementTable *theElementTable  = 0 ; if (!theElementTable) theElementTable= G4Element::GetElementTable();
   vElement.clear();
   vElement.resize( numberOfElements );
   for ( size_t i = 0; i < numberOfElements; ++i )
   {

      G4Element* theElement = (*theElementTable)[i];
      vElement[i] = false;

      // isotope
      G4int nIso = (*theElementTable)[i]->GetNumberOfIsotopes();
      G4int Z = static_cast<G4int> ((*theElementTable)[i]->GetZ());
      if ( nIso!=0 )
      {
         G4bool found_at_least_one = false; 
         for ( G4int i1 = 0; i1 < nIso; i1++ )
         {
             G4int A = theElement->GetIsotope(i1)->GetN();

             if ( isThisNewIsotope( Z , A ) ) 
             { 

                std::stringstream ss; 
                ss << dirName << aFSType << Z << "_" << A << "_" << theNames.GetName( Z-1 );
                filename = ss.str();
                std::fstream file;
                file.open ( filename , std::fstream::in );
                G4int dummy;
                file >> dummy;
                if ( file.good() ) 
                {

                   //G4cout << "Found file for Z=" << Z << ", A=" << A << ", as " << filename << G4endl;
                   found_at_least_one = true;

                   // read the file
                   G4PhysicsVector* aPhysVec = readAFile ( &file );

                   //Regist 

                   registAPhysicsVector( Z , A , aPhysVec );

                }
                else 
                { 
                   //G4cout << "No file for "<< reactionType << " Z=" << Z << ", A=" << A << G4endl;
                }

                file.close();

             }
             else
             {
                found_at_least_one = TRUE;
             }
          }

          if ( found_at_least_one ) vElement[i] = true;

       }
       else
       {
          G4StableIsotopes theStableOnes;
          G4int first = theStableOnes.GetFirstIsotope( Z );
          G4bool found_at_least_one = FALSE; 
          for ( G4int i1 = 0; i1 < theStableOnes.GetNumberOfIsotopes( static_cast<G4int>(theElement->GetZ() ) ); i1++)
          {
             G4int A = theStableOnes.GetIsotopeNucleonCount( first+i1 );

             if ( isThisNewIsotope( Z , A ) ) 
             {

                std::stringstream ss; 
                ss << dirName << aFSType << Z << "_" << A << "_" << theNames.GetName( Z-1 );
                filename = ss.str();

                std::fstream file;
                file.open ( filename , std::fstream::in );
                G4int dummy;
                file >> dummy;
                if ( file.good() ) 
                {
                   //G4cout << "Found file for Z=" << Z << ", A=" << A << ", as " << filename << G4endl;
                   found_at_least_one = TRUE;
                   //Read the file

                   G4PhysicsVector* aPhysVec = readAFile ( &file );

                   //Regist the PhysicsVector
                   registAPhysicsVector( Z , A , aPhysVec );

                }
                else 
                { 
                   //G4cout << "No file for "<< reactionType << " Z=" << Z << ", A=" << A << G4endl;
                }

                file.close();
             }
             else
             {
                found_at_least_one = TRUE;
             }
          }

          if ( found_at_least_one ) vElement[i] = true;

       }

   }

}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4ParticleHPJENDLHEData::DumpPhysicsTable ( const G4ParticleDefinition &  aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 237 of file G4ParticleHPJENDLHEData.cc.

{
  if(&aP!=G4Neutron::Neutron()) 
     throw G4HadronicException(__FILE__, __LINE__, "Attempt to use NeutronHP data for particles other than neutrons!!!");  
//  G4cout << "G4ParticleHPJENDLHEData::DumpPhysicsTable still to be implemented"<<G4endl;
}

Here is the call graph for this function:

G4double G4ParticleHPJENDLHEData::GetCrossSection	(	const G4DynamicParticle *	,
		const G4Element *	,
		G4double	aT
	)

G4bool G4ParticleHPJENDLHEData::IsApplicable	(	const G4DynamicParticle *	aP,
		const G4Element *	anE
	)

Definition at line 41 of file G4ParticleHPJENDLHEData.cc.

{

   G4bool result = true;
   G4double eKin = aP->GetKineticEnergy();
   //if(eKin>20*MeV||aP->GetDefinition()!=G4Neutron::Neutron()) result = false;
   if ( eKin < 20*MeV || 3*GeV < eKin || aP->GetDefinition()!=G4Neutron::Neutron() ) 
   {
      result = false;
   } 
// Element Check 
   else if ( !(vElement[ anE->GetIndex() ]) ) result = false;

   return result;

}

Here is the call graph for this function:

G4bool G4ParticleHPJENDLHEData::IsZAApplicable ( const G4DynamicParticle *  , G4double  , G4double    )

inline

Definition at line 60 of file G4ParticleHPJENDLHEData.hh.

   { return false; }

---

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

• geant4.10.03.p01/source/processes/hadronic/models/particle_hp/include/G4ParticleHPJENDLHEData.hh
• geant4.10.03.p01/source/processes/hadronic/models/particle_hp/src/G4ParticleHPJENDLHEData.cc