#include <G4LENDCrossSection.hh>

Inheritance 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

Detailed Description

Definition at line 57 of file G4LENDCrossSection.hh.

Constructor & Destructor Documentation

G4LENDCrossSection::G4LENDCrossSection ( const G4String name = "" )

Definition at line 83 of file G4LENDCrossSection.cc.

References G4LENDManager::GetInstance(), python.hepunit::MeV, proj, G4VCrossSectionDataSet::SetMaxKinEnergy(), and G4VCrossSectionDataSet::SetMinKinEnergy().

 :G4VCrossSectionDataSet( nam )
 {
 
    proj = NULL; //will be set in an inherited class
    //default_evaluation = "endl99";
    default_evaluation = "ENDF.B-VII.0";
 
    allow_nat = false;
    allow_any = false;
 
    SetMinKinEnergy(  0*MeV );
    SetMaxKinEnergy( 20*MeV );
 
    lend_manager = G4LENDManager::GetInstance(); 
 
 }

G4LENDCrossSection::~G4LENDCrossSection ( )

Definition at line 101 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

void G4LENDCrossSection::AllowAnyCandidateTarget ( )

inline

Definition at line 96 of file G4LENDCrossSection.hh.

Referenced by G4NeutronLENDBuilder::Build(), and G4HadronElasticPhysicsLEND::ConstructProcess().

96 { allow_any = true; };

void G4LENDCrossSection::AllowNaturalAbundanceTarget ( )

inline

Definition at line 95 of file G4LENDCrossSection.hh.

95 { allow_nat = true; };

void G4LENDCrossSection::BuildPhysicsTable ( const G4ParticleDefinition & )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 112 of file G4LENDCrossSection.cc.

References create_used_target_map().

 {
    create_used_target_map();
 }

void G4LENDCrossSection::ChangeDefaultEvaluation ( G4String name_tmp )

inline

Definition at line 94 of file G4LENDCrossSection.hh.

Referenced by G4NeutronLENDBuilder::Build(), and G4HadronElasticPhysicsLEND::ConstructProcess().

94 { default_evaluation = name_tmp; };

void G4LENDCrossSection::create_used_target_map ( )

protected

Definition at line 274 of file G4LENDCrossSection.cc.

References G4LENDUsedTarget::AllowAny(), G4LENDUsedTarget::AllowNat(), G4cout, G4endl, G4Element::GetElementTable(), G4Element::GetIsotope(), G4NistElementBuilder::GetIsotopeAbundance(), G4Isotope::Getm(), G4Isotope::GetN(), G4VCrossSectionDataSet::GetName(), G4LENDManager::GetNistElementBuilder(), G4NistElementBuilder::GetNistFirstIsotopeN(), G4LENDManager::GetNucleusEncoding(), G4Element::GetNumberOfElements(), G4Element::GetNumberOfIsotopes(), G4NistElementBuilder::GetNumberOfNistIsotopes(), G4Isotope::GetZ(), G4Element::GetZ(), int(), proj, G4LENDManager::RequestChangeOfVerboseLevel(), and G4VCrossSectionDataSet::verboseLevel.

Referenced by BuildPhysicsTable().

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

void G4LENDCrossSection::DumpPhysicsTable ( const G4ParticleDefinition & aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 117 of file G4LENDCrossSection.cc.

References python.hepunit::barn, python.hepunit::eV, G4cout, G4endl, getLENDCrossSection(), G4VCrossSectionDataSet::GetName(), python.hepunit::MeV, and proj.

 {
 
   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 * std::pow ( 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;
 
    }
 
 }

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 56 of file G4LENDCrossSection.cc.

References G4DynamicParticle::GetKineticEnergy(), getLENDCrossSection(), G4Isotope::Getm(), G4LENDManager::GetNucleusEncoding(), G4Material::GetTemperature(), and python.hepunit::second.

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

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

inlineprotectedvirtual

Definition at line 129 of file G4LENDCrossSection.hh.

Referenced by DumpPhysicsTable(), and GetIsoCrossSection().

129 { return 0.0; };

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

protected

Definition at line 351 of file G4LENDCrossSection.cc.

References test::a, and test::b.

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

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 46 of file G4LENDCrossSection.cc.

References G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetKineticEnergy(), G4VCrossSectionDataSet::GetMaxKinEnergy(), G4VCrossSectionDataSet::GetMinKinEnergy(), and proj.

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

Field Documentation

G4ParticleDefinition* G4LENDCrossSection::proj

protected

Definition at line 127 of file G4LENDCrossSection.hh.

Referenced by create_used_target_map(), DumpPhysicsTable(), G4LENDCaptureCrossSection::G4LENDCaptureCrossSection(), G4LENDCrossSection(), G4LENDElasticCrossSection::G4LENDElasticCrossSection(), G4LENDFissionCrossSection::G4LENDFissionCrossSection(), G4LENDInelasticCrossSection::G4LENDInelasticCrossSection(), and IsIsoApplicable().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Field Documentation