#include <G4NeutronHPCaptureData.hh>

Inheritance diagram for G4NeutronHPCaptureData:

Public Member Functions
	G4NeutronHPCaptureData ()

	~G4NeutronHPCaptureData ()

G4bool	IsIsoApplicable (const G4DynamicParticle , G4int, G4int, const G4Element , const G4Material *)

G4double	GetIsoCrossSection (const G4DynamicParticle , G4int, G4int, const G4Isotope , const G4Element , const G4Material )

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , G4double aT)

void	BuildPhysicsTable (const G4ParticleDefinition &)

void	DumpPhysicsTable (const G4ParticleDefinition &)

void	IgnoreOnFlightDopplerBroadening ()

void	EnableOnFlightDopplerBroadening ()

G4int	GetVerboseLevel () const

void	SetVerboseLevel (G4int)

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

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 49 of file G4NeutronHPCaptureData.hh.

Constructor & Destructor Documentation

G4NeutronHPCaptureData::G4NeutronHPCaptureData ( )

Definition at line 46 of file G4NeutronHPCaptureData.cc.

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

 :G4VCrossSectionDataSet("NeutronHPCaptureXS")
 {
    SetMinKinEnergy( 0*MeV );                                   
    SetMaxKinEnergy( 20*MeV );                                   
 
    ke_cache = 0.0;
    xs_cache = 0.0;
    element_cache = NULL;
    material_cache = NULL;
 
    theCrossSections = 0;
    onFlightDB = true;
 
    //BuildPhysicsTable(*G4Neutron::Neutron());
 }

G4NeutronHPCaptureData::~G4NeutronHPCaptureData ( )

Definition at line 63 of file G4NeutronHPCaptureData.cc.

References G4PhysicsTable::clearAndDestroy().

 {
    if ( theCrossSections != 0 ) theCrossSections->clearAndDestroy();
 
    delete theCrossSections;
 }

Member Function Documentation

void G4NeutronHPCaptureData::BuildPhysicsTable ( const G4ParticleDefinition & aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 109 of file G4NeutronHPCaptureData.cc.

References G4PhysicsTable::clearAndDestroy(), G4cout, G4endl, G4ThreadLocal, G4Element::GetElementTable(), G4Element::GetNumberOfElements(), G4NeutronHPData::Instance(), G4NeutronHPData::MakePhysicsVector(), G4Neutron::Neutron(), and G4PhysicsTable::push_back().

 {
   if(&aP!=G4Neutron::Neutron()) 
      throw G4HadronicException(__FILE__, __LINE__, "Attempt to use NeutronHP data for particles other than neutrons!!!");  
 
 //080428
    if ( getenv( "G4NEUTRONHP_NEGLECT_DOPPLER" ) ) 
    {
       G4cout << "Find environment variable of \"G4NEUTRONHP_NEGLECT_DOPPLER\"." << G4endl;
       G4cout << "On the fly Doppler broadening will be neglect in the cross section calculation of capture reaction of neutrons (<20MeV)." << G4endl;
       onFlightDB = false;
    }
   
   size_t numberOfElements = G4Element::GetNumberOfElements();
   // G4cout << "CALLED G4NeutronHPCaptureData::BuildPhysicsTable "<<numberOfElements<<G4endl;
    // TKDB
    //if ( theCrossSections == 0 ) theCrossSections = new G4PhysicsTable( numberOfElements );
    if ( theCrossSections == NULL ) 
       theCrossSections = new G4PhysicsTable( numberOfElements );
    else
       theCrossSections->clearAndDestroy();
 
   // make a PhysicsVector for each element
 
   static G4ThreadLocal G4ElementTable *theElementTable  = 0 ; if (!theElementTable) theElementTable= G4Element::GetElementTable();
   for( size_t i=0; i<numberOfElements; ++i )
   {
      if(getenv("CaptureDataIndexDebug"))
      {
        G4int index_debug = ((*theElementTable)[i])->GetIndex();
        G4cout << "IndexDebug "<< i <<" "<<index_debug<<G4endl;
      }
      G4PhysicsVector* physVec = G4NeutronHPData::
       Instance()->MakePhysicsVector((*theElementTable)[i], this);
      theCrossSections->push_back(physVec);
   }
 }

void G4NeutronHPCaptureData::DumpPhysicsTable ( const G4ParticleDefinition & aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 147 of file G4NeutronHPCaptureData.cc.

References python.hepunit::barn, python.hepunit::eV, G4cout, G4endl, G4ThreadLocal, G4Element::GetElementTable(), G4Element::GetNumberOfElements(), python.hepunit::MeV, and G4Neutron::Neutron().

 {
   if(&aP!=G4Neutron::Neutron()) 
      throw G4HadronicException(__FILE__, __LINE__, "Attempt to use NeutronHP data for particles other than neutrons!!!");  
 
 //
 // Dump element based cross section
 // range 10e-5 eV to 20 MeV
 // 10 point per decade
 // in barn
 //
 
    G4cout << G4endl;
    G4cout << G4endl;
    G4cout << "Capture Cross Section of Neutron HP"<< G4endl;
    G4cout << "(Pointwise cross-section at 0 Kelvin.)" << G4endl;
    G4cout << G4endl;
    G4cout << "Name of Element" << G4endl;
    G4cout << "Energy[eV]  XS[barn]" << G4endl;
    G4cout << G4endl;
 
    size_t numberOfElements = G4Element::GetNumberOfElements();
    static G4ThreadLocal G4ElementTable *theElementTable  = 0 ; if (!theElementTable) theElementTable= G4Element::GetElementTable();
 
    for ( size_t i = 0 ; i < numberOfElements ; ++i )
    {
 
       G4cout << (*theElementTable)[i]->GetName() << G4endl;
 
       G4int ie = 0;
 
       for ( ie = 0 ; ie < 130 ; ie++ )
       {
          G4double eKinetic = 1.0e-5 * std::pow ( 10.0 , ie/10.0 ) *eV;
          G4bool outOfRange = false;
 
          if ( eKinetic < 20*MeV )
          {
             G4cout << eKinetic/eV << " " << (*((*theCrossSections)(i))).GetValue(eKinetic, outOfRange)/barn << G4endl;
          }
 
       }
 
       G4cout << G4endl;
    }
 
 
 //  G4cout << "G4NeutronHPCaptureData::DumpPhysicsTable still to be implemented"<<G4endl;
 }

void G4NeutronHPCaptureData::EnableOnFlightDopplerBroadening ( )

inline

Definition at line 81 of file G4NeutronHPCaptureData.hh.

81 { onFlightDB = true; };

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

Definition at line 200 of file G4NeutronHPCaptureData.cc.

References buffer, G4DynamicParticle::GetDefinition(), G4Element::GetIndex(), G4DynamicParticle::GetKineticEnergy(), G4ReactionProduct::GetKineticEnergy(), G4ReactionProduct::GetMass(), G4DynamicParticle::GetMomentum(), G4ReactionProduct::GetMomentum(), G4Element::GetN(), G4NucleiProperties::GetNuclearMass(), G4ParticleDefinition::GetPDGMass(), G4Nucleus::GetThermalNucleus(), G4Element::GetZ(), python.hepunit::k_Boltzmann, python.hepunit::kelvin, G4ReactionProduct::Lorentz(), CLHEP::Hep3Vector::mag(), G4INCL::Math::max(), G4Neutron::Neutron(), and G4ReactionProduct::SetMomentum().

Referenced by GetIsoCrossSection().

 {
   G4double result = 0;
   G4bool outOfRange;
   G4int index = anE->GetIndex();
 
   // prepare neutron
   G4double eKinetic = aP->GetKineticEnergy();
 
 //if ( getenv( "G4NEUTRONHP_NEGLECT_DOPPLER" ) )
 //080428
   if ( !onFlightDB )
   {
      G4double factor = 1.0;
      if ( eKinetic < aT * k_Boltzmann ) 
      {
         // below 0.1 eV neutrons 
         // Have to do some, but now just igonre.   
         // Will take care after performance check.  
         // factor = factor * targetV;
      }
      return ( (*((*theCrossSections)(index))).GetValue(eKinetic, outOfRange) )* factor; 
   }
 
   G4ReactionProduct theNeutron( aP->GetDefinition() );
   theNeutron.SetMomentum( aP->GetMomentum() );
   theNeutron.SetKineticEnergy( eKinetic );
 
   // prepare thermal nucleus
   G4Nucleus aNuc;
   G4double eps = 0.0001;
   G4double theA = anE->GetN();
   G4double theZ = anE->GetZ();
   G4double eleMass; 
   eleMass = G4NucleiProperties::GetNuclearMass( static_cast<G4int>(theA+eps) , static_cast<G4int>(theZ+eps) ) / G4Neutron::Neutron()->GetPDGMass();
   
   G4ReactionProduct boosted;
   G4double aXsection;
   
   // MC integration loop
   G4int counter = 0;
   G4double buffer = 0;
   G4int size = G4int(std::max(10., aT/60*kelvin));
   G4ThreeVector neutronVelocity = 1./G4Neutron::Neutron()->GetPDGMass()*theNeutron.GetMomentum();
   G4double neutronVMag = neutronVelocity.mag();
 
   while(counter == 0 || std::abs(buffer-result/std::max(1,counter)) > 0.03*buffer)
   {
     if(counter) buffer = result/counter;
     while (counter<size)
     {
       counter ++;
       G4ReactionProduct aThermalNuc = aNuc.GetThermalNucleus(eleMass, aT);
       boosted.Lorentz(theNeutron, aThermalNuc);
       G4double theEkin = boosted.GetKineticEnergy();
       aXsection = (*((*theCrossSections)(index))).GetValue(theEkin, outOfRange);
       // velocity correction, or luminosity factor...
       G4ThreeVector targetVelocity = 1./aThermalNuc.GetMass()*aThermalNuc.GetMomentum();
       aXsection *= (targetVelocity-neutronVelocity).mag()/neutronVMag;
       result += aXsection;
     }
     size += size;
   }
   result /= counter;
 /*
   // Checking impact of  G4NEUTRONHP_NEGLECT_DOPPLER
   G4cout << " result " << result << " " 
          << (*((*theCrossSections)(index))).GetValue(eKinetic, outOfRange) << " " 
          << (*((*theCrossSections)(index))).GetValue(eKinetic, outOfRange) /result << G4endl;
 */
   return result;
 }

G4double G4NeutronHPCaptureData::GetIsoCrossSection	(	const G4DynamicParticle *	dp,
		G4int	,
		G4int	,
		const G4Isotope *	,
		const G4Element *	element,
		const G4Material *	material
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 83 of file G4NeutronHPCaptureData.cc.

References GetCrossSection(), G4DynamicParticle::GetKineticEnergy(), and eplot::material.

 {
    if ( dp->GetKineticEnergy() == ke_cache && element == element_cache &&  material == material_cache ) return xs_cache;
 
    ke_cache = dp->GetKineticEnergy();
    element_cache = element;
    material_cache = material;
    G4double xs = GetCrossSection( dp , element , material->GetTemperature() );
    xs_cache = xs;
    return xs;
 }

G4int G4NeutronHPCaptureData::GetVerboseLevel ( ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 273 of file G4NeutronHPCaptureData.cc.

References G4NeutronHPManager::GetInstance(), and G4NeutronHPManager::GetVerboseLevel().

 {
    return G4NeutronHPManager::GetInstance()->GetVerboseLevel();
 }

void G4NeutronHPCaptureData::IgnoreOnFlightDopplerBroadening ( )

inline

Definition at line 80 of file G4NeutronHPCaptureData.hh.

80 { onFlightDB = false; };

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 70 of file G4NeutronHPCaptureData.cc.

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

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

void G4NeutronHPCaptureData::SetVerboseLevel ( G4int newValue )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 277 of file G4NeutronHPCaptureData.cc.

References G4NeutronHPManager::GetInstance(), and G4NeutronHPManager::SetVerboseLevel().

 {
    G4NeutronHPManager::GetInstance()->SetVerboseLevel(newValue);
 }

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation