G4XNNElasticLowE Class Reference

#include <G4XNNElasticLowE.hh>

Inheritance diagram for G4XNNElasticLowE:

Public Member Functions
virtual G4double	CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
	G4XNNElasticLowE ()
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 G4XNNElasticLowE &right) const
G4bool	operator== (const G4VCrossSectionSource &right) const
G4bool	operator== (const G4XNNElasticLowE &right) const
virtual void	Print () const
virtual void	PrintAll (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
virtual	~G4XNNElasticLowE ()

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

Private Attributes
G4double	_eMax
G4double	_eMin
std::map< const G4ParticleDefinition *, G4PhysicsVector *, std::less< const G4ParticleDefinition * > >	xMap

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	npTable [101]
static const G4double	ppTable [101]
static const G4int	tableSize = 101

Detailed Description

Definition at line 39 of file G4XNNElasticLowE.hh.

Constructor & Destructor Documentation

G4XNNElasticLowE() [1/2]

G4XNNElasticLowE::G4XNNElasticLowE

(

)

Definition at line 89 of file G4XNNElasticLowE.cc.

{ 
  // Cross-sections are available in the range (_eMin,_eMax)
 
  _eMin = _eMinTable * GeV;
  _eMax = G4Exp(G4Log(_eMinTable) + tableSize * _eStepLog) * GeV;
  if (_eMin < _lowLimit)
    throw G4HadronicException(__FILE__, __LINE__, "G4XNNElasticLowE::G4XNNElasticLowE - Low energy limit not valid");    
  if (_highLimit > _eMax)
    throw G4HadronicException(__FILE__, __LINE__, "G4XNNElasticLowE::G4XNNElasticLowE - High energy limit not valid");    
  G4PhysicsVector* pp = new G4PhysicsLogVector(_eMin,_eMax,tableSize);
 
  _eMin = G4Exp(G4Log(_eMinTable)-_eStepLog)*GeV;
  if (_eMin < _lowLimit)
    throw G4HadronicException(__FILE__, __LINE__, "G4XNNElasticLowE::G4XNNElasticLowE - Low energy limit not valid");
  G4PhysicsVector* np = new G4PhysicsLogVector(_eMin,_eMax,tableSize);
 
  G4int i;
  for (i=0; i<tableSize; i++)
    {
      G4double value = ppTable[i] * millibarn;
      pp->PutValue(i,value);
      value = npTable[i] * millibarn;
      np->PutValue(i,value);
    }
  xMap[G4Proton::ProtonDefinition()] = pp;
  xMap[G4Neutron::NeutronDefinition()] = np;
}

References _eMax, _eMin, _eMinTable, _eStepLog, _highLimit, _lowLimit, G4Exp(), G4Log(), GeV, millibarn, G4Neutron::NeutronDefinition(), npTable, G4InuclParticleNames::pp, ppTable, G4Proton::ProtonDefinition(), G4PhysicsVector::PutValue(), tableSize, and xMap.

~G4XNNElasticLowE()

G4XNNElasticLowE::~G4XNNElasticLowE ( )

virtual

Definition at line 119 of file G4XNNElasticLowE.cc.

{
  delete xMap[G4Proton::ProtonDefinition()];
  delete xMap[G4Neutron::NeutronDefinition()];
}

References G4Neutron::NeutronDefinition(), G4Proton::ProtonDefinition(), and xMap.

G4XNNElasticLowE() [2/2]

G4XNNElasticLowE::G4XNNElasticLowE ( const G4XNNElasticLowE &  right )

private

Member Function Documentation

CrossSection()

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

virtual

Implements G4VCrossSectionSource.

Definition at line 139 of file G4XNNElasticLowE.cc.

{
  G4double sigma = 0.;
  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
  G4bool dummy = false;
 
  const G4ParticleDefinition * key = FindKeyParticle(trk1,trk2);
 
  typedef std::map <const G4ParticleDefinition *, G4PhysicsVector*, std::less<const G4ParticleDefinition *> > StringPhysMap;
 
  if (xMap.find(key)!= xMap.end())
    {
 
      StringPhysMap::const_iterator iter;
      for (iter = xMap.begin(); iter != xMap.end(); ++iter)
    {
      const G4ParticleDefinition * str = (*iter).first;
          if (str == key)
        {
          G4PhysicsVector* physVector = (*iter).second; 
          //     G4PhysicsVector* physVector = xMap[key];
          if (sqrtS >= _eMin && sqrtS <= _eMax)
        {
          sigma = physVector->GetValue(sqrtS,dummy);
        } else if ( sqrtS < _eMin )
                {
                  sigma = physVector->GetValue(_eMin,dummy);
                }
        //G4cout << " sqrtS / sigma " << sqrtS/GeV << " / " <<
        //          sigma/millibarn << G4endl;
        }
    }
    }
  return sigma;
}

References _eMax, _eMin, G4VCrossSectionSource::FindKeyParticle(), G4KineticTrack::Get4Momentum(), G4PhysicsVector::GetValue(), and xMap.

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 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 * G4XNNElasticLowE::GetComponents ( ) const

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 53 of file G4XNNElasticLowE.hh.

{ return 0; }

HighLimit()

virtual G4double G4XNNElasticLowE::HighLimit ( ) const

inlinevirtual

Reimplemented from G4VCrossSectionSource.

Definition at line 61 of file G4XNNElasticLowE.hh.

{ 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(), IsValid(), G4XnpElasticLowE::IsValid(), G4XnpTotalLowE::IsValid(), G4XPDGElastic::IsValid(), and G4XPDGTotal::IsValid().

IsValid()

G4bool G4XNNElasticLowE::IsValid ( G4double  e ) const

virtual

Reimplemented from G4VCrossSectionSource.

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

virtual

Implements G4VCrossSectionSource.

Definition at line 239 of file G4XNNElasticLowE.cc.

{
  G4String name("NNElasticLowE");
  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 G4XNNElasticLowE::operator!=

(

const G4XNNElasticLowE &

right

)

const

Definition at line 132 of file G4XNNElasticLowE.cc.

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

operator=()

const G4XNNElasticLowE & G4XNNElasticLowE::operator= ( const G4XNNElasticLowE &  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 G4XNNElasticLowE::operator==

(

const G4XNNElasticLowE &

right

)

const

Definition at line 126 of file G4XNNElasticLowE.cc.

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

Print()

void G4XNNElasticLowE::Print ( ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 176 of file G4XNNElasticLowE.cc.

{
  // Dump the pp cross-section table
 
  G4cout << Name() << ", pp cross-section: " << G4endl;
 
  G4bool dummy = false;
  G4int i;
  const G4ParticleDefinition * key = G4Proton::ProtonDefinition();
  G4PhysicsVector* pp = 0;
 
  typedef std::map <const G4ParticleDefinition *, G4PhysicsVector*, std::less<const G4ParticleDefinition *> > StringPhysMap;
  StringPhysMap::const_iterator iter;
 
  for (iter = xMap.begin(); iter != xMap.end(); ++iter)
    {
      const G4ParticleDefinition * str = (*iter).first;
      if (str == key)
    {
      pp = (*iter).second; 
    }
    }
  
  if (pp != 0)
    {   
      for (i=0; i<tableSize; i++)
    {
      G4double e = pp->GetLowEdgeEnergy(i);
      G4double sigma = pp->GetValue(e,dummy) / millibarn;
      G4cout << i << ") e = " << e / GeV << " GeV ---- Cross section = " << sigma << " mb " << G4endl;
    }
    }
  
  // Dump the np cross-section table
 
  G4cout << Name() << ", np cross-section: " << G4endl;
 
  key = G4Neutron::NeutronDefinition();
  G4PhysicsVector* np = 0;
  for (iter = xMap.begin(); iter != xMap.end(); ++iter)
    {
      const G4ParticleDefinition * str = (*iter).first;
      if (str == key)
    {
      np = (*iter).second; 
    }
    }
  
  //  G4PhysicsVector* np = xMap[G4Neutron::NeutronDefinition()->GetParticleName()];
  
  if (np != 0)
    {   
      for (i=0; i<tableSize; i++)
    {
      G4double e = np->GetLowEdgeEnergy(i);
      G4double sigma = np->GetValue(e,dummy) / millibarn;
      G4cout << i << ") e = " << e / GeV << " GeV ---- Cross section = " << sigma << " mb " << G4endl;
    }
    }
  G4VCrossSectionSource::Print();
}

References G4cout, G4endl, G4PhysicsVector::GetLowEdgeEnergy(), G4PhysicsVector::GetValue(), GeV, millibarn, Name(), G4Neutron::NeutronDefinition(), G4InuclParticleNames::pp, G4VCrossSectionSource::Print(), G4Proton::ProtonDefinition(), tableSize, and xMap.

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 G4XNNElasticLowE::_eMax

private

Definition at line 83 of file G4XNNElasticLowE.hh.

Referenced by CrossSection(), and G4XNNElasticLowE().

_eMin

G4double G4XNNElasticLowE::_eMin

private

Definition at line 82 of file G4XNNElasticLowE.hh.

Referenced by CrossSection(), and G4XNNElasticLowE().

_eMinTable

const G4double G4XNNElasticLowE::_eMinTable = 1.8964808

staticprivate

Definition at line 77 of file G4XNNElasticLowE.hh.

Referenced by G4XNNElasticLowE().

_eStepLog

const G4double G4XNNElasticLowE::_eStepLog = 0.01

staticprivate

Definition at line 78 of file G4XNNElasticLowE.hh.

Referenced by G4XNNElasticLowE().

_highLimit

const G4double G4XNNElasticLowE::_highLimit = 3.*GeV

staticprivate

Definition at line 73 of file G4XNNElasticLowE.hh.

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

_lowLimit

const G4double G4XNNElasticLowE::_lowLimit = 0.

staticprivate

Definition at line 72 of file G4XNNElasticLowE.hh.

Referenced by G4XNNElasticLowE(), and IsValid().

npTable

const G4double G4XNNElasticLowE::npTable

staticprivate

Initial value:

= 
{ 
  1500.00, 
  248.20, 93.38, 55.26, 44.50, 41.33, 38.48, 37.20, 35.98,
  35.02, 34.47, 32.48, 30.76, 29.46, 28.53, 27.84, 27.20,
  26.53, 25.95, 25.59, 25.46, 25.00, 24.49, 24.08, 23.86,
  23.17, 22.70, 21.88, 21.48, 20.22, 19.75, 18.97, 18.39,
  17.98, 17.63, 17.21, 16.72, 16.68, 16.58, 16.42, 16.22,
  15.98, 15.71, 15.42, 15.14, 14.87, 14.65, 14.44, 14.26,
  14.10, 13.95, 13.80, 13.64, 13.47, 13.29, 13.09, 12.89,
  12.68, 12.47, 12.27, 12.06, 11.84, 11.76, 11.69, 11.60,
  11.50, 11.41, 11.29, 11.17, 11.06, 10.93, 10.81, 10.68,
  10.56, 10.44, 10.33, 10.21, 10.12, 10.03,  9.96,  9.89,
  9.83,  9.80,  9.77,  9.75,  9.74,  9.74,  9.74,  9.76,
  9.73,  9.70,  9.68,  9.65,  9.63,  9.60,  9.57,  9.55,
  9.52,  9.49,  9.46,  9.43
}

Definition at line 75 of file G4XNNElasticLowE.hh.

Referenced by G4XNNElasticLowE().

ppTable

const G4double G4XNNElasticLowE::ppTable

staticprivate

Initial value:

= 
{ 
  60.00, 
  33.48, 26.76, 25.26, 24.55, 23.94, 23.77, 23.72, 23.98,
  25.48, 27.52, 27.72, 27.21, 25.80, 26.00, 24.32, 23.81,
  24.37, 24.36, 23.13, 22.43, 21.71, 21.01, 20.83, 20.74,
  20.25, 20.10, 20.59, 20.04, 20.83, 20.84, 21.07, 20.83,
  20.79, 21.88, 21.15, 20.92, 19.00, 18.60, 17.30, 17.00,
  16.70, 16.50, 16.20, 15.80, 15.57, 15.20, 15.00, 14.60,
  14.20, 14.00, 13.80, 13.60, 13.40, 13.20, 13.00, 12.85,
  12.70, 12.60, 12.50, 12.40, 12.30, 12.20, 12.10, 12.00,
  11.90, 11.80, 11.75, 11.70, 11.64, 11.53, 11.41, 11.31,
  11.22, 11.13, 11.05, 10.97, 10.89, 10.82, 10.75, 10.68,
  10.61, 10.54, 10.48, 10.41, 10.35, 10.28, 10.22, 10.16,
  10.13, 10.10, 10.08, 10.05, 10.02,  9.99,  9.96,  9.93,
  9.90,  9.87,  9.84,  9.80       
}

Definition at line 74 of file G4XNNElasticLowE.hh.

Referenced by G4XNNElasticLowE().

tableSize

const G4int G4XNNElasticLowE::tableSize = 101

staticprivate

Definition at line 76 of file G4XNNElasticLowE.hh.

Referenced by G4XNNElasticLowE(), and Print().

xMap

std::map<const G4ParticleDefinition *, G4PhysicsVector*, std::less<const G4ParticleDefinition *> > G4XNNElasticLowE::xMap

private

Definition at line 80 of file G4XNNElasticLowE.hh.

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

---

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

• source/processes/hadronic/models/im_r_matrix/include/G4XNNElasticLowE.hh
• source/processes/hadronic/models/im_r_matrix/src/G4XNNElasticLowE.cc