#include <G4XnpTotalLowE.hh>

Inheritance diagram for G4XnpTotalLowE:

Public Member Functions
virtual G4double	CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

	G4XnpTotalLowE ()

virtual const G4CrossSectionVector *	GetComponents () const

virtual G4double	HighLimit () const

virtual G4bool	IsValid (G4double e) const

virtual G4double	LowLimit () const

virtual G4String	Name () const

G4bool	operator!= (const G4VCrossSectionSource &right) const

G4bool	operator!= (const G4XnpTotalLowE &right) const

G4bool	operator== (const G4VCrossSectionSource &right) const

G4bool	operator== (const G4XnpTotalLowE &right) const

virtual void	Print () const

virtual void	PrintAll (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

virtual	~G4XnpTotalLowE ()

Protected Member Functions
G4double	FcrossX (G4double e, G4double e0, G4double sigma, G4double eParam, G4double power) const

const G4ParticleDefinition *	FindKeyParticle (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

const G4ParticleDefinition *	FindLightParticle (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

G4bool	InLimits (G4double e, G4double eLow, G4double eHigh) const

Private Member Functions
	G4XnpTotalLowE (const G4XnpTotalLowE &right)

const G4XnpTotalLowE &	operator= (const G4XnpTotalLowE &right)

Private Attributes
G4double	_eMax

G4double	_eMin

G4PhysicsVector *	_sigma

Static Private Attributes
static const G4double	_eMinTable = 1.8964808

static const G4double	_eStepLog = 0.01

static const G4double	_highLimit = 3.*GeV

static const G4double	_lowLimit = 0.

static const G4double	_sigmaTable [101]

static const G4int	_tableSize = 101

Detailed Description

Definition at line 51 of file G4XnpTotalLowE.hh.

Constructor & Destructor Documentation

◆ G4XnpTotalLowE() [1/2]

G4XnpTotalLowE::G4XnpTotalLowE ( )

Definition at line 88 of file G4XnpTotalLowE.cc.

{ 
  // Cross-sections are available in the range (_eMin,_eMax)
 
  _eMin = _eMinTable * GeV;
  _eMin = G4Exp(G4Log(_eMinTable)-_eStepLog)*GeV;
  _eMax = G4Exp(G4Log(_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 G4PhysicsLogVector(_eMin,_eMax,_tableSize);
  G4int i;
  for (i=0; i<_tableSize; i++)
    {
      G4double value = _sigmaTable[i] * millibarn;
      _sigma->PutValue(i,value);
    }
}

References _eMax, _eMin, _eMinTable, _eStepLog, _highLimit, _lowLimit, _sigma, _sigmaTable, _tableSize, G4Exp(), G4Log(), GeV, millibarn, and G4PhysicsVector::PutValue().

◆ ~G4XnpTotalLowE()

G4XnpTotalLowE::~G4XnpTotalLowE ( )

virtual

Definition at line 114 of file G4XnpTotalLowE.cc.

{
  if (_sigma) delete _sigma;
  _sigma=0;
}

References _sigma.

◆ G4XnpTotalLowE() [2/2]

G4XnpTotalLowE::G4XnpTotalLowE ( const G4XnpTotalLowE & right )

private

Member Function Documentation

◆ CrossSection()

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

virtual

Implements G4VCrossSectionSource.

Definition at line 133 of file G4XnpTotalLowE.cc.

{
  G4double sigma = 0.;
  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
  G4bool dummy = false;
 
  const G4ParticleDefinition* proton = G4Proton::ProtonDefinition();
  const G4ParticleDefinition* neutron = G4Neutron::NeutronDefinition();
 
  const G4ParticleDefinition* def1 = trk1.GetDefinition();
  const 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;
}

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

◆ FcrossX()

G4double G4VCrossSectionSource::FcrossX	(	G4double	e,
		G4double	e0,
		G4double	sigma,
		G4double	eParam,
		G4double	power
	)		const

protectedinherited

Definition at line 174 of file G4VCrossSectionSource.cc.

{
  G4double result = 0.;
 
  G4double denom = eParam*eParam + (e-e0)*(e-e0);
  if (denom > 0.) 
  {
    G4double value = (2.* eParam * sigma * (e-e0) / denom) * G4Pow::GetInstance()->powA(((e0 + eParam) / e), power);
    result = std::max(0., value);
  }
  return result;
}     

References G4Pow::GetInstance(), G4INCL::Math::max(), and G4Pow::powA().

◆ FindKeyParticle()

const G4ParticleDefinition * G4VCrossSectionSource::FindKeyParticle	(	const G4KineticTrack &	trk1,
		const G4KineticTrack &	trk2
	)		const

protectedinherited

Definition at line 49 of file G4VCrossSectionSource.cc.

{
  const G4ParticleDefinition * result;
  
  const G4ParticleDefinition * p1 = trk1.GetDefinition();
  const G4ParticleDefinition * p2 = trk2.GetDefinition();
  
  if( (p1==G4Proton::Proton() && p2==G4Proton::Proton() ) ||
      (p1==G4Neutron::Neutron() && p2==G4Neutron::Neutron()) )
  {
    result = G4Proton::Proton();
  }
  else if( (p1==G4Neutron::Neutron() && p2==G4Proton::Proton()) ||
           (p2==G4Neutron::Neutron() && p1==G4Proton::Proton()) )
  {
    result = G4Neutron::Neutron();
  }
  else
  {
    throw G4HadronicException(__FILE__, __LINE__, "G4VCrossSectionSource: unklnown particles in FindKeyParticle");
  }
  return result;
}

References G4KineticTrack::GetDefinition(), G4Neutron::Neutron(), and G4Proton::Proton().

Referenced by G4XNNElasticLowE::CrossSection(), and G4XNNTotalLowE::CrossSection().

◆ FindLightParticle()

const G4ParticleDefinition * G4VCrossSectionSource::FindLightParticle	(	const G4KineticTrack &	trk1,
		const G4KineticTrack &	trk2
	)		const

protectedinherited

Definition at line 158 of file G4VCrossSectionSource.cc.

{
  G4double mass1 = trk1.GetDefinition()->GetPDGMass();
  G4double mass2 = trk2.GetDefinition()->GetPDGMass();
  if (mass1 < mass2)
    {
      return trk1.GetDefinition();
    }
  else
    {
      return trk2.GetDefinition();
    }
}

References G4KineticTrack::GetDefinition(), and G4ParticleDefinition::GetPDGMass().

Referenced by G4XMesonBaryonElastic::CrossSection().

◆ GetComponents()

virtual const G4CrossSectionVector * G4XnpTotalLowE::GetComponents ( ) const

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 65 of file G4XnpTotalLowE.hh.

65{return 0; }

◆ HighLimit()

virtual G4double G4XnpTotalLowE::HighLimit ( ) const

inlinevirtual

Reimplemented from G4VCrossSectionSource.

Definition at line 73 of file G4XnpTotalLowE.hh.

73{ return _highLimit; }

References _highLimit.

◆ InLimits()

G4bool G4VCrossSectionSource::InLimits	(	G4double	e,
		G4double	eLow,
		G4double	eHigh
	)		const

protectedinherited

Definition at line 133 of file G4VCrossSectionSource.cc.

{
  G4bool answer = false;
  if (e >= eLow && e <= eHigh) answer = true;
  return answer;
}

Referenced by G4XAnnihilationChannel::IsValid(), G4XAqmElastic::IsValid(), G4XAqmTotal::IsValid(), G4XMesonBaryonElastic::IsValid(), G4XNNElasticLowE::IsValid(), G4XnpElasticLowE::IsValid(), IsValid(), G4XPDGElastic::IsValid(), and G4XPDGTotal::IsValid().

◆ IsValid()

G4bool G4XnpTotalLowE::IsValid ( G4double e ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 185 of file G4XnpTotalLowE.cc.

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

References _highLimit, _lowLimit, and G4VCrossSectionSource::InLimits().

◆ LowLimit()

G4double G4VCrossSectionSource::LowLimit ( ) const

virtualinherited

Reimplemented in G4XPDGElastic, and G4XPDGTotal.

Definition at line 140 of file G4VCrossSectionSource.cc.

{
  return 0.; 
}

Referenced by G4CrossSectionPatch::CrossSection(), G4VCrossSectionSource::IsValid(), and G4CrossSectionPatch::Transition().

◆ Name()

G4String G4XnpTotalLowE::Name ( ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 178 of file G4XnpTotalLowE.cc.

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

References G4InuclParticleNames::name().

Referenced by Print().

◆ operator!=() [1/2]

G4bool G4VCrossSectionSource::operator!= ( const G4VCrossSectionSource & right ) const

inherited

Definition at line 80 of file G4VCrossSectionSource.cc.

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

◆ operator!=() [2/2]

G4bool G4XnpTotalLowE::operator!= ( const G4XnpTotalLowE & right ) const

Definition at line 127 of file G4XnpTotalLowE.cc.

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

◆ operator=()

const G4XnpTotalLowE & G4XnpTotalLowE::operator= ( const G4XnpTotalLowE & right )

private

◆ operator==() [1/2]

G4bool G4VCrossSectionSource::operator== ( const G4VCrossSectionSource & right ) const

inherited

Definition at line 74 of file G4VCrossSectionSource.cc.

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

◆ operator==() [2/2]

G4bool G4XnpTotalLowE::operator== ( const G4XnpTotalLowE & right ) const

Definition at line 121 of file G4XnpTotalLowE.cc.

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

◆ Print()

void G4XnpTotalLowE::Print ( ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 159 of file G4XnpTotalLowE.cc.

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

References _sigma, _tableSize, G4cout, G4endl, G4PhysicsVector::GetLowEdgeEnergy(), G4PhysicsVector::GetValue(), GeV, millibarn, Name(), and G4VCrossSectionSource::Print().

◆ PrintAll()

void G4VCrossSectionSource::PrintAll	(	const G4KineticTrack &	trk1,
		const G4KineticTrack &	trk2
	)		const

virtualinherited

Definition at line 107 of file G4VCrossSectionSource.cc.

{
  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
  G4double sigma = CrossSection(trk1,trk2) / millibarn;
  G4cout << "---- " << Name() << ": "
     << "Ecm = " << sqrtS / GeV << " GeV -  " 
     << " Cross section = " << sigma << " mb "
     << G4endl;
 
  G4int nComponents = 0;
  const G4CrossSectionVector* components = GetComponents();
  if (components != 0)
    {
      nComponents = components->size();
    }
  G4int i;
  for (i=0; i<nComponents; i++)
    {
      G4cout << "* Component " << i << ": ";
      G4CrossSectionSourcePtr componentPtr = (*components)[i];
      G4VCrossSectionSource* component = componentPtr();
      component->PrintAll(trk1,trk2);
    }
}

References G4VCrossSectionSource::CrossSection(), G4cout, G4endl, G4KineticTrack::Get4Momentum(), G4VCrossSectionSource::GetComponents(), GeV, millibarn, G4VCrossSectionSource::Name(), and G4VCrossSectionSource::PrintAll().

Referenced by G4VCrossSectionSource::PrintAll().

Field Documentation

◆ _eMax

G4double G4XnpTotalLowE::_eMax

private

Definition at line 93 of file G4XnpTotalLowE.hh.

Referenced by CrossSection(), and G4XnpTotalLowE().

◆ _eMin

G4double G4XnpTotalLowE::_eMin

private

Definition at line 92 of file G4XnpTotalLowE.hh.

Referenced by CrossSection(), and G4XnpTotalLowE().

◆ _eMinTable

const G4double G4XnpTotalLowE::_eMinTable = 1.8964808

staticprivate

Definition at line 88 of file G4XnpTotalLowE.hh.

Referenced by G4XnpTotalLowE().

◆ _eStepLog

const G4double G4XnpTotalLowE::_eStepLog = 0.01

staticprivate

Definition at line 89 of file G4XnpTotalLowE.hh.

Referenced by G4XnpTotalLowE().

◆ _highLimit

const G4double G4XnpTotalLowE::_highLimit = 3.*GeV

staticprivate

Definition at line 85 of file G4XnpTotalLowE.hh.

Referenced by G4XnpTotalLowE(), HighLimit(), and IsValid().

◆ _lowLimit

const G4double G4XnpTotalLowE::_lowLimit = 0.

staticprivate

Definition at line 84 of file G4XnpTotalLowE.hh.

Referenced by G4XnpTotalLowE(), and IsValid().

◆ _sigma

G4PhysicsVector* G4XnpTotalLowE::_sigma

private

Definition at line 91 of file G4XnpTotalLowE.hh.

Referenced by CrossSection(), G4XnpTotalLowE(), Print(), and ~G4XnpTotalLowE().

◆ _sigmaTable

const G4double G4XnpTotalLowE::_sigmaTable

staticprivate

Initial value:

= 
{ 
0,  
20, 93.38, 55.26, 44.50, 41.33, 38.48, 37.20, 35.98,
02, 34.47, 34.37, 34.67, 35.23, 35.97, 36.75, 37.37,
77, 38.03, 38.40, 38.83, 39.26, 39.67, 40.06, 40.45,
79, 41.06, 41.31, 41.52, 41.70, 41.81, 41.87, 41.98,
12, 42.29, 42.55, 42.82, 43.01, 43.12, 43.16, 43.14,
06, 42.95, 42.81, 42.67, 42.54, 42.45, 42.38, 42.33,
30, 42.29, 42.28, 42.26, 42.24, 42.21, 42.17, 42.14,
10, 42.07, 42.06, 42.05, 42.04, 42.03, 42.02, 42.00,
97, 41.94, 41.89, 41.84, 41.79, 41.73, 41.67, 41.61,
55, 41.49, 41.44, 41.38, 41.34, 41.31, 41.29, 41.28,
27, 41.28, 41.30, 41.33, 41.36, 41.40, 41.44, 41.49,
50, 41.51, 41.51, 41.51, 41.52, 41.51, 41.51, 41.50,
50, 41.49, 41.47, 41.46
}

Definition at line 86 of file G4XnpTotalLowE.hh.

Referenced by G4XnpTotalLowE().

◆ _tableSize

const G4int G4XnpTotalLowE::_tableSize = 101

staticprivate

Definition at line 87 of file G4XnpTotalLowE.hh.

Referenced by G4XnpTotalLowE(), and Print().

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

source/processes/hadronic/models/im_r_matrix/include/G4XnpTotalLowE.hh
source/processes/hadronic/models/im_r_matrix/src/G4XnpTotalLowE.cc

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4XnpTotalLowE() [1/2]

◆ ~G4XnpTotalLowE()

◆ G4XnpTotalLowE() [2/2]

Member Function Documentation

◆ CrossSection()

◆ FcrossX()

◆ FindKeyParticle()

◆ FindLightParticle()

◆ GetComponents()

◆ HighLimit()

◆ InLimits()

◆ IsValid()

◆ LowLimit()

◆ Name()

◆ operator!=() [1/2]

◆ operator!=() [2/2]

◆ operator=()

◆ operator==() [1/2]

◆ operator==() [2/2]

◆ Print()

◆ PrintAll()

Field Documentation

◆ _eMax

◆ _eMin

◆ _eMinTable

◆ _eStepLog

◆ _highLimit

◆ _lowLimit

◆ _sigma

◆ _sigmaTable

◆ _tableSize