G4LENDCrossSection Class Reference

#include <G4LENDCrossSection.hh>

Inheritance diagram for G4LENDCrossSection:

Collaboration diagram for G4LENDCrossSection:

Public Member Functions
	G4LENDCrossSection (const G4String name="")
	~G4LENDCrossSection ()
G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int, G4int, const G4Element *, const G4Material *)
G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int, G4int, const G4Isotope *, const G4Element *, const G4Material *)
void	BuildPhysicsTable (const G4ParticleDefinition &)
void	DumpPhysicsTable (const G4ParticleDefinition &)
void	ChangeDefaultEvaluation (G4String name_tmp)
void	AllowNaturalAbundanceTarget ()
void	AllowAnyCandidateTarget ()
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, 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 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

Protected Member Functions
void	create_used_target_map ()
virtual G4double	getLENDCrossSection (G4GIDI_target *, G4double, G4double)
G4double	GetUltraLowEnergyExtrapolatedXS (G4double, G4double, G4double, G4double, G4double)
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)

Protected Attributes
G4ParticleDefinition *	proj
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Private Member Functions
void	recreate_used_target_map ()

Private Attributes
std::map< G4int, G4LENDUsedTarget *>	usedTarget_map
G4String	default_evaluation
G4bool	allow_nat
G4bool	allow_any
G4LENDManager *	lend_manager

Detailed Description

Definition at line 57 of file G4LENDCrossSection.hh.

Constructor & Destructor Documentation

G4LENDCrossSection()

G4LENDCrossSection::G4LENDCrossSection

(

const G4String

name = ""

)

Definition at line 84 of file G4LENDCrossSection.cc.

:G4VCrossSectionDataSet( nam )
{

   proj = NULL; //will be set in an inherited class
   //default_evaluation = "endl99";
   //default_evaluation = "ENDF.B-VII.0";
   default_evaluation = "ENDF/BVII.1";

   allow_nat = false;
   allow_any = false;

   SetMinKinEnergy(  0*MeV );
   SetMaxKinEnergy( 20*MeV );

   lend_manager = G4LENDManager::GetInstance(); 

}

Here is the call graph for this function:

~G4LENDCrossSection()

G4LENDCrossSection::~G4LENDCrossSection

(

)

Definition at line 103 of file G4LENDCrossSection.cc.

{

   for ( std::map< G4int , G4LENDUsedTarget* >::iterator 
         it = usedTarget_map.begin() ; it != usedTarget_map.end() ; it ++ )
   { 
      delete it->second;  
   }

}

Member Function Documentation

AllowAnyCandidateTarget()

void G4LENDCrossSection::AllowAnyCandidateTarget ( )

inline

Definition at line 96 of file G4LENDCrossSection.hh.

{ allow_any = true; };

Here is the caller graph for this function:

AllowNaturalAbundanceTarget()

void G4LENDCrossSection::AllowNaturalAbundanceTarget ( )

inline

Definition at line 95 of file G4LENDCrossSection.hh.

{ allow_nat = true; };

BuildPhysicsTable()

void G4LENDCrossSection::BuildPhysicsTable ( const G4ParticleDefinition &  )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 114 of file G4LENDCrossSection.cc.

{
   create_used_target_map();
}

Here is the call graph for this function:

ChangeDefaultEvaluation()

void G4LENDCrossSection::ChangeDefaultEvaluation ( G4String  name_tmp )

inline

Definition at line 94 of file G4LENDCrossSection.hh.

{ default_evaluation = name_tmp; };

Here is the caller graph for this function:

create_used_target_map()

void G4LENDCrossSection::create_used_target_map ( )

protected

Definition at line 276 of file G4LENDCrossSection.cc.

{

   lend_manager->RequestChangeOfVerboseLevel( verboseLevel );

   size_t numberOfElements = G4Element::GetNumberOfElements();
   static const G4ElementTable* theElementTable = G4Element::GetElementTable();

   for ( size_t i = 0 ; i < numberOfElements ; ++i )
   {

      const G4Element* anElement = (*theElementTable)[i];
      G4int numberOfIsotope = anElement->GetNumberOfIsotopes(); 

      if ( numberOfIsotope > 0 )
      {
      // User Defined Abundances   
         for ( G4int i_iso = 0 ; i_iso < numberOfIsotope ; i_iso++ )
         {
            G4int iZ = anElement->GetIsotope( i_iso )->GetZ();
            G4int iA = anElement->GetIsotope( i_iso )->GetN();
            G4int iIsomer = anElement->GetIsotope( i_iso )->Getm();

            //G4LENDUsedTarget* aTarget = new G4LENDUsedTarget ( G4Neutron::Neutron() , default_evaluation , iZ , iA );  
            G4LENDUsedTarget* aTarget = new G4LENDUsedTarget ( proj , default_evaluation , iZ , iA );  
            if ( allow_nat == true ) aTarget->AllowNat();
            if ( allow_any == true ) aTarget->AllowAny();
            usedTarget_map.insert( std::pair< G4int , G4LENDUsedTarget* > ( lend_manager->GetNucleusEncoding( iZ , iA , iIsomer ) , aTarget ) );
         }
      }
      else
      {
      // Natural Abundances   
         G4NistElementBuilder* nistElementBuild = lend_manager->GetNistElementBuilder();
         G4int iZ = int ( anElement->GetZ() );
         //G4cout << nistElementBuild->GetNumberOfNistIsotopes( int ( anElement->GetZ() ) ) << G4endl;
         G4int numberOfNistIso = nistElementBuild->GetNumberOfNistIsotopes( int ( anElement->GetZ() ) ); 

         for ( G4int ii = 0 ; ii < numberOfNistIso ; ii++ )
         {
            //G4cout << nistElementBuild->GetIsotopeAbundance( iZ , nistElementBuild->GetNistFirstIsotopeN( iZ ) + i ) << G4endl;
            if ( nistElementBuild->GetIsotopeAbundance( iZ , nistElementBuild->GetNistFirstIsotopeN( iZ ) + ii ) > 0 )
            {
               G4int iMass = nistElementBuild->GetNistFirstIsotopeN( iZ ) + ii;  
               //G4cout << iZ << " " << nistElementBuild->GetNistFirstIsotopeN( iZ ) + i << " " << nistElementBuild->GetIsotopeAbundance ( iZ , iMass ) << G4endl;  
               G4int iIsomer = 0; 

               G4LENDUsedTarget* aTarget = new G4LENDUsedTarget ( proj , default_evaluation , iZ , iMass );  
               if ( allow_nat == true ) aTarget->AllowNat();
               if ( allow_any == true ) aTarget->AllowAny();
               usedTarget_map.insert( std::pair< G4int , G4LENDUsedTarget* > ( lend_manager->GetNucleusEncoding( iZ , iMass , iIsomer ) , aTarget ) );

            }

         }
      }
   }

   G4cout << "Dump UsedTarget for " << GetName() << G4endl;
   //G4cout << "Requested Evaluation, Z , A -> Actual Evaluation, Z , A(0=Nat) , Pointer of Target" << G4endl;
   G4cout << "Requested Evaluation, Z , A -> Actual Evaluation, Z , A(0=Nat) " << G4endl;
   for ( std::map< G4int , G4LENDUsedTarget* >::iterator 
         it = usedTarget_map.begin() ; it != usedTarget_map.end() ; it ++ )
   {
      G4cout 
         << " " << it->second->GetWantedEvaluation() 
         << ", " << it->second->GetWantedZ() 
         << ", " << it->second->GetWantedA() 
         << " -> " << it->second->GetActualEvaluation() 
         << ", " << it->second->GetActualZ() 
         << ", " << it->second->GetActualA() 
         //<< ", " << it->second->GetTarget() 
         << G4endl; 
   } 

}

Here is the call graph for this function:

Here is the caller graph for this function:

DumpPhysicsTable()

void G4LENDCrossSection::DumpPhysicsTable ( const G4ParticleDefinition &  aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 119 of file G4LENDCrossSection.cc.

{

  if ( &aP != proj ) 
     throw G4HadronicException(__FILE__, __LINE__, "Attempt to use LEND data for particles other than neutrons!!!");  

   G4cout << G4endl;
   G4cout << "Dump Cross Sections of " << GetName() << G4endl;
   G4cout << "(Pointwise cross-section at 300 Kelvin.)" << G4endl;
   G4cout << G4endl;

   G4cout << "Target informaiton " << G4endl;

   for ( std::map< G4int , G4LENDUsedTarget* >::iterator 
         it = usedTarget_map.begin() ; it != usedTarget_map.end() ; it ++ )
   {
      G4cout 
         << "Wanted " << it->second->GetWantedEvaluation() 
         << ", Z= " << it->second->GetWantedZ() 
         << ", A= " << it->second->GetWantedA() 
         << "; Actual " << it->second->GetActualEvaluation() 
         << ", Z= " << it->second->GetActualZ() 
         << ", A= " << it->second->GetActualA() 
         << ", " << it->second->GetTarget() 
         << G4endl; 

      G4int ie = 0;

      G4GIDI_target* aTarget = it->second->GetTarget();
      G4double aT = 300;
      for ( ie = 0 ; ie < 130 ; ie++ )
      {
         G4double ke = 1.0e-5 * G4Pow::GetInstance()->powA ( 10.0 , ie/10.0 ) *eV;

         if ( ke < 20*MeV )
         {
            G4cout << "  " << GetName() << ", cross section at " << ke/eV << " [eV] = " << getLENDCrossSection ( aTarget , ke , aT )/barn << " [barn] " << G4endl;
         }
      }
      G4cout << G4endl;

   }

}

Here is the call graph for this function:

GetIsoCrossSection()

G4double G4LENDCrossSection::GetIsoCrossSection ( const G4DynamicParticle *  dp, G4int  iZ, G4int  iA, const G4Isotope *  isotope, const G4Element *  , const G4Material *  material  )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 57 of file G4LENDCrossSection.cc.

{

   G4double xs = 0.0;
   G4double ke = dp->GetKineticEnergy();
   G4double temp = material->GetTemperature();
   G4int iM = isotope->Getm();

   G4GIDI_target* aTarget = usedTarget_map.find( lend_manager->GetNucleusEncoding( iZ , iA , iM ) )->second->GetTarget();

   xs = getLENDCrossSection ( aTarget , ke , temp );

   return xs;
}

Here is the call graph for this function:

getLENDCrossSection()

virtual G4double G4LENDCrossSection::getLENDCrossSection ( G4GIDI_target *  , G4double  , G4double    )

inlineprotectedvirtual

Reimplemented in G4LENDInelasticCrossSection, G4LENDCaptureCrossSection, G4LENDElasticCrossSection, and G4LENDFissionCrossSection.

Definition at line 129 of file G4LENDCrossSection.hh.

{ return 0.0; };

Here is the call graph for this function:

Here is the caller graph for this function:

GetUltraLowEnergyExtrapolatedXS()

G4double G4LENDCrossSection::GetUltraLowEnergyExtrapolatedXS ( G4double  x1, G4double  x2, G4double  y1, G4double  y2, G4double  ke  )

protected

Definition at line 354 of file G4LENDCrossSection.cc.

{
   //XS propotinal to 1/v at low energy -> 1/root(E) 
   //XS = a * 1/root(E) + b  
   G4double a = ( y2 - y1 ) / ( 1/std::sqrt(x2) - 1/std::sqrt(x1) );
   G4double b = y1 - a * 1/std::sqrt(x1);
   G4double result = a * 1/std::sqrt(ke) + b;
   return result;
}

Here is the caller graph for this function:

IsIsoApplicable()

G4bool G4LENDCrossSection::IsIsoApplicable ( const G4DynamicParticle *  dp, G4int  , G4int  , const G4Element *  , const G4Material *    )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 47 of file G4LENDCrossSection.cc.

{
   G4double eKin = dp->GetKineticEnergy();
   if ( dp->GetDefinition() != proj ) return false;
   if ( eKin > GetMaxKinEnergy() || eKin < GetMinKinEnergy() ) return false;

   return true;
}

Here is the call graph for this function:

recreate_used_target_map()

void G4LENDCrossSection::recreate_used_target_map ( )

private

Definition at line 262 of file G4LENDCrossSection.cc.

{
   for ( std::map< G4int , G4LENDUsedTarget* >::iterator 
         it = usedTarget_map.begin() ; it != usedTarget_map.end() ; it ++ )
   { 
      delete it->second;  
   }
   usedTarget_map.clear();

   create_used_target_map();
}

Here is the call graph for this function:

Member Data Documentation

allow_any

G4bool G4LENDCrossSection::allow_any

private

Definition at line 120 of file G4LENDCrossSection.hh.

allow_nat

G4bool G4LENDCrossSection::allow_nat

private

Definition at line 119 of file G4LENDCrossSection.hh.

default_evaluation

G4String G4LENDCrossSection::default_evaluation

private

Definition at line 118 of file G4LENDCrossSection.hh.

lend_manager

G4LENDManager* G4LENDCrossSection::lend_manager

private

Definition at line 122 of file G4LENDCrossSection.hh.

proj

G4ParticleDefinition* G4LENDCrossSection::proj

protected

Definition at line 127 of file G4LENDCrossSection.hh.

usedTarget_map

std::map< G4int , G4LENDUsedTarget* > G4LENDCrossSection::usedTarget_map

private

Definition at line 96 of file G4LENDCrossSection.hh.

---

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

• G4LENDCrossSection.hh
• G4LENDCrossSection.cc