G4XnpTotalLowE Class Reference

#include <G4XnpTotalLowE.hh>

Inheritance diagram for G4XnpTotalLowE:

Public Member Functions
	G4XnpTotalLowE ()
virtual	~G4XnpTotalLowE ()
G4bool	operator== (const G4XnpTotalLowE &right) const
G4bool	operator!= (const G4XnpTotalLowE &right) const
virtual G4double	CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
virtual const G4CrossSectionVector *	GetComponents () const
virtual G4bool	IsValid (G4double e) const
virtual void	Print () const
virtual G4String	Name () const
virtual G4double	HighLimit () const

---

Detailed Description

Definition at line 51 of file G4XnpTotalLowE.hh.

---

Constructor & Destructor Documentation

G4XnpTotalLowE::G4XnpTotalLowE

(

)

Definition at line 87 of file G4XnpTotalLowE.cc.

References G4PhysicsVector::PutValue().

{ 
  // Cross-sections are available in the range (_eMin,_eMax)

  _eMin = _eMinTable * GeV;
  _eMin = std::exp(std::log(_eMinTable)-_eStepLog)*GeV;
  _eMax = std::exp(std::log(_eMinTable) + _tableSize * _eStepLog) * GeV;

  // Protections: validity limits must be compatible with available data
//  @@GF  this ought to be _lowLimit < _eMin
  if (_eMin < _lowLimit)
    throw G4HadronicException(__FILE__, __LINE__, "G4XnpTotalLowE::G4XnpTotalLowE - Low energy limit not valid");
    
  if (_highLimit > _eMax)
    throw G4HadronicException(__FILE__, __LINE__, "G4XnpTotalLowE::G4XnpTotalLowE - High energy limit not valid");
    
  _sigma = new G4PhysicsLnVector(_eMin,_eMax,_tableSize);
  G4int i;
  for (i=0; i<_tableSize; i++)
    {
      G4double value = _sigmaTable[i] * millibarn;
      _sigma->PutValue(i,value);
    }
}

G4XnpTotalLowE::~G4XnpTotalLowE

(

)

[virtual]

Definition at line 113 of file G4XnpTotalLowE.cc.

{
  delete _sigma;
}

---

Member Function Documentation

G4double G4XnpTotalLowE::CrossSection	(	const G4KineticTrack &	trk1,
		const G4KineticTrack &	trk2
	)			const [virtual]

Implements G4VCrossSectionSource.

Definition at line 131 of file G4XnpTotalLowE.cc.

References G4KineticTrack::Get4Momentum(), G4KineticTrack::GetDefinition(), G4PhysicsVector::GetValue(), neutron, G4Neutron::NeutronDefinition(), G4InuclParticleNames::proton, and G4Proton::ProtonDefinition().

{
  G4double sigma = 0.;
  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
  G4bool dummy = false;

  G4ParticleDefinition* proton = G4Proton::ProtonDefinition();
  G4ParticleDefinition* neutron = G4Neutron::NeutronDefinition();   

  G4ParticleDefinition* def1 = trk1.GetDefinition();
  G4ParticleDefinition* def2 = trk2.GetDefinition();
  if ( (def1 == proton && def2 == neutron) ||
       (def1 == neutron && def2 == proton) )
    {
      if (sqrtS >= _eMin && sqrtS <= _eMax)
        {
          sigma = _sigma->GetValue(sqrtS,dummy);
        } else if ( sqrtS < _eMin )
        {
          sigma = _sigma->GetValue(_eMin,dummy);
        }
    }

  return sigma;
}

virtual const G4CrossSectionVector* G4XnpTotalLowE::GetComponents

(

)

const [inline, virtual]

Implements G4VCrossSectionSource.

Definition at line 65 of file G4XnpTotalLowE.hh.

{return 0; }

virtual G4double G4XnpTotalLowE::HighLimit

(

)

const [inline, virtual]

Reimplemented from G4VCrossSectionSource.

Definition at line 73 of file G4XnpTotalLowE.hh.

{ return _highLimit; }

G4bool G4XnpTotalLowE::IsValid

(

G4double

e

)

const [virtual]

Reimplemented from G4VCrossSectionSource.

Definition at line 183 of file G4XnpTotalLowE.cc.

References G4VCrossSectionSource::InLimits().

{
  G4bool answer = InLimits(e,_lowLimit,_highLimit);

  return answer;
}

G4String G4XnpTotalLowE::Name

(

)

const [virtual]

Implements G4VCrossSectionSource.

Definition at line 176 of file G4XnpTotalLowE.cc.

Referenced by Print().

{
  G4String name("NNTotalLowE");
  return name;
}

G4bool G4XnpTotalLowE::operator!=

(

const G4XnpTotalLowE &

right

)

const

Definition at line 125 of file G4XnpTotalLowE.cc.

{
  return (this != (G4XnpTotalLowE *) &right);
}

G4bool G4XnpTotalLowE::operator==

(

const G4XnpTotalLowE &

right

)

const

Definition at line 119 of file G4XnpTotalLowE.cc.

{
  return (this == (G4XnpTotalLowE *) &right);
}

void G4XnpTotalLowE::Print

(

)

const [virtual]

Reimplemented from G4VCrossSectionSource.

Definition at line 157 of file G4XnpTotalLowE.cc.

References G4cout, G4endl, G4PhysicsVector::GetLowEdgeEnergy(), G4PhysicsVector::GetValue(), Name(), and G4VCrossSectionSource::Print().

{
  // Dump the cross-section table

  G4cout << Name() << "Cross-section table: " << G4endl;
  G4bool dummy = false;
  G4int i;

  for (i=0; i<_tableSize; i++)
    {
      G4double e = _sigma->GetLowEdgeEnergy(i) / GeV;
      G4double sigma = _sigma->GetValue(e,dummy) / millibarn;
      G4cout << i << ") e = " << e << " GeV ---- Cross section = " << sigma << " mb " << G4endl;
    }

  G4VCrossSectionSource::Print();
}

---

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

• G4XnpTotalLowE.hh
• G4XnpTotalLowE.cc

---