#include <G4XnpElasticLowE.hh>

Inheritance diagram for G4XnpElasticLowE:

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

	G4XnpElasticLowE ()

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 G4XnpElasticLowE &right) const

G4bool	operator== (const G4VCrossSectionSource &right) const

G4bool	operator== (const G4XnpElasticLowE &right) const

virtual void	Print () const

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

virtual	~G4XnpElasticLowE ()

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
	G4XnpElasticLowE (const G4XnpElasticLowE &right)

const G4XnpElasticLowE &	operator= (const G4XnpElasticLowE &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 G4XnpElasticLowE.hh.

Constructor & Destructor Documentation

◆ G4XnpElasticLowE() [1/2]

G4XnpElasticLowE::G4XnpElasticLowE ( )

Definition at line 88 of file G4XnpElasticLowE.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
 
  if (_eMin < _lowLimit)
    throw G4HadronicException(__FILE__, __LINE__, "G4XnpElasticLowE::G4XnpElasticLowE - Low energy limit not valid");
    
  if (_highLimit > _eMax)
    throw G4HadronicException(__FILE__, __LINE__, "G4XnpElasticLowE::G4XnpElasticLowE - 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().

◆ ~G4XnpElasticLowE()

G4XnpElasticLowE::~G4XnpElasticLowE ( )

virtual

Definition at line 114 of file G4XnpElasticLowE.cc.

{
   delete _sigma;
   _sigma = 0;
}

References _sigma.

◆ G4XnpElasticLowE() [2/2]

G4XnpElasticLowE::G4XnpElasticLowE ( const G4XnpElasticLowE & right )

private

Member Function Documentation

◆ CrossSection()

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

virtual

Implements G4VCrossSectionSource.

Definition at line 133 of file G4XnpElasticLowE.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 * G4XnpElasticLowE::GetComponents ( ) const

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 65 of file G4XnpElasticLowE.hh.

65{ return 0; }

◆ HighLimit()

virtual G4double G4XnpElasticLowE::HighLimit ( ) const

inlinevirtual

Reimplemented from G4VCrossSectionSource.

Definition at line 73 of file G4XnpElasticLowE.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(), IsValid(), G4XnpTotalLowE::IsValid(), G4XPDGElastic::IsValid(), and G4XPDGTotal::IsValid().

◆ IsValid()

G4bool G4XnpElasticLowE::IsValid ( G4double e ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 184 of file G4XnpElasticLowE.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 G4XnpElasticLowE::Name ( ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 177 of file G4XnpElasticLowE.cc.

{
  G4String name("npElasticLowE");
  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 G4XnpElasticLowE::operator!= ( const G4XnpElasticLowE & right ) const

Definition at line 127 of file G4XnpElasticLowE.cc.

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

◆ operator=()

const G4XnpElasticLowE & G4XnpElasticLowE::operator= ( const G4XnpElasticLowE & 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 G4XnpElasticLowE::operator== ( const G4XnpElasticLowE & right ) const

Definition at line 121 of file G4XnpElasticLowE.cc.

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

◆ Print()

void G4XnpElasticLowE::Print ( ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 159 of file G4XnpElasticLowE.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 G4XnpElasticLowE::_eMax

private

Definition at line 93 of file G4XnpElasticLowE.hh.

Referenced by CrossSection(), and G4XnpElasticLowE().

◆ _eMin

G4double G4XnpElasticLowE::_eMin

private

Definition at line 92 of file G4XnpElasticLowE.hh.

Referenced by CrossSection(), and G4XnpElasticLowE().

◆ _eMinTable

const G4double G4XnpElasticLowE::_eMinTable = 1.8964808

staticprivate

Definition at line 88 of file G4XnpElasticLowE.hh.

Referenced by G4XnpElasticLowE().

◆ _eStepLog

const G4double G4XnpElasticLowE::_eStepLog = 0.01

staticprivate

Definition at line 89 of file G4XnpElasticLowE.hh.

Referenced by G4XnpElasticLowE().

◆ _highLimit

const G4double G4XnpElasticLowE::_highLimit = 3.*GeV

staticprivate

Definition at line 85 of file G4XnpElasticLowE.hh.

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

◆ _lowLimit

const G4double G4XnpElasticLowE::_lowLimit = 0.

staticprivate

Definition at line 84 of file G4XnpElasticLowE.hh.

Referenced by G4XnpElasticLowE(), and IsValid().

◆ _sigma

G4PhysicsVector* G4XnpElasticLowE::_sigma

private

Definition at line 91 of file G4XnpElasticLowE.hh.

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

◆ _sigmaTable

const G4double G4XnpElasticLowE::_sigmaTable

staticprivate

Initial value:

= 
{ 
0,  
20, 93.38, 55.26, 44.50, 41.33, 38.48, 37.20, 35.98,
02, 34.47, 32.48, 30.76, 29.46, 28.53, 27.84, 27.20,
53, 25.95, 25.59, 25.46, 25.00, 24.49, 24.08, 23.86,
17, 22.70, 21.88, 21.48, 20.22, 19.75, 18.97, 18.39,
98, 17.63, 17.21, 16.72, 16.68, 16.58, 16.42, 16.22,
98, 15.71, 15.42, 15.14, 14.87, 14.65, 14.44, 14.26,
10, 13.95, 13.80, 13.64, 13.47, 13.29, 13.09, 12.89,
68, 12.47, 12.27, 12.06, 11.84, 11.76, 11.69, 11.60,
50, 11.41, 11.29, 11.17, 11.06, 10.93, 10.81, 10.68,
56, 10.44, 10.33, 10.21, 10.12, 10.03,  9.96,  9.89,
83,  9.80,  9.77,  9.75,  9.74,  9.74,  9.74,  9.76,
73,  9.70,  9.68,  9.65,  9.63,  9.60,  9.57,  9.55,
52,  9.49,  9.46,  9.43
}

Definition at line 86 of file G4XnpElasticLowE.hh.

Referenced by G4XnpElasticLowE().

◆ _tableSize

const G4int G4XnpElasticLowE::_tableSize = 101

staticprivate

Definition at line 87 of file G4XnpElasticLowE.hh.

Referenced by G4XnpElasticLowE(), and Print().

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

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

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4XnpElasticLowE() [1/2]

◆ ~G4XnpElasticLowE()

◆ G4XnpElasticLowE() [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