G4NeutronHPJENDLHEData Class Reference

#include <G4NeutronHPJENDLHEData.hh>

Inheritance diagram for G4NeutronHPJENDLHEData:

Public Member Functions
	G4NeutronHPJENDLHEData ()
	G4NeutronHPJENDLHEData (G4String, G4ParticleDefinition *)
	~G4NeutronHPJENDLHEData ()
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	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual void	CrossSectionDescription (std::ostream &) const
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 46 of file G4NeutronHPJENDLHEData.hh.

Constructor & Destructor Documentation

G4NeutronHPJENDLHEData() [1/2]

G4NeutronHPJENDLHEData::G4NeutronHPJENDLHEData

(

)

Definition at line 58 of file G4NeutronHPJENDLHEData.cc.

{
   ;
}

G4NeutronHPJENDLHEData() [2/2]

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

Definition at line 65 of file G4NeutronHPJENDLHEData.cc.

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

~G4NeutronHPJENDLHEData()

G4NeutronHPJENDLHEData::~G4NeutronHPJENDLHEData

(

)

Definition at line 74 of file G4NeutronHPJENDLHEData.cc.

{
   ; 
   //delete theCrossSections;
}

Member Function Documentation

BuildPhysicsTable()

void G4NeutronHPJENDLHEData::BuildPhysicsTable ( const G4ParticleDefinition &  aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 82 of file G4NeutronHPJENDLHEData.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/";
   G4NeutronHPNames 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 const G4ElementTable *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;
 
       }
 
   }
 
}

Referenced by G4NeutronHPJENDLHEData().

DumpPhysicsTable()

void G4NeutronHPJENDLHEData::DumpPhysicsTable ( const G4ParticleDefinition &  aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 219 of file G4NeutronHPJENDLHEData.cc.

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

GetCrossSection()

G4double G4NeutronHPJENDLHEData::GetCrossSection	(	const G4DynamicParticle *	aP,
		const G4Element *	anE,
		G4double	aT
	)

Definition at line 228 of file G4NeutronHPJENDLHEData.cc.

{
 
   // Primary energy >20MeV
   // Thus
   // Not take account of Doppler broadening 
   // also
   // Not take account of Target thermal motions
 
   G4double result = 0;
 
   G4double ek = aP->GetKineticEnergy();
 
   G4int nIso = anE->GetNumberOfIsotopes();
   G4int Z = static_cast<G4int> ( anE->GetZ() );
   if ( nIso!=0 )
   {
      for ( G4int i1 = 0; i1 < nIso; i1++ )
      {
 
         G4int A = anE->GetIsotope(i1)->GetN();
         G4double frac = anE->GetRelativeAbundanceVector()[ i1 ];   // This case do NOT request "*perCent".    
 
         result += frac * getXSfromThisIsotope( Z , A , ek );
         //G4cout << reactionType << " XS in barn " << Z << " " << A << " " << frac << " " << getXSfromThisIsotope( Z , A , ek )/barn << G4endl; 
 
      }
   }
   else
   {
 
      G4StableIsotopes theStableOnes;
      G4int first = theStableOnes.GetFirstIsotope( Z );
      for ( G4int i1 = 0; i1 < theStableOnes.GetNumberOfIsotopes( static_cast<G4int>(anE->GetZ() ) ); i1++)
      {
         
         G4int A = theStableOnes.GetIsotopeNucleonCount( first+i1 );
         G4double frac = theStableOnes.GetAbundance( first+i1 )*perCent;  // This case request "*perCent". 
 
         result += frac * getXSfromThisIsotope( Z , A , ek );
         //G4cout << reactionType << " XS in barn " << Z << " " << A << " " << frac << " " << getXSfromThisIsotope( Z , A , ek )/barn << G4endl; 
          
      }
   }
   return result;
 
}

IsApplicable()

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

Definition at line 39 of file G4NeutronHPJENDLHEData.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;
 
}

IsZAApplicable()

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

inline

Definition at line 58 of file G4NeutronHPJENDLHEData.hh.

   { return false; }

---

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

• G4NeutronHPJENDLHEData.hh
• G4NeutronHPJENDLHEData.cc