#include <G4XNNTotalLowE.hh>

Inheritance diagram for G4XNNTotalLowE:

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

	G4XNNTotalLowE ()

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

virtual void	Print () const

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

virtual	~G4XNNTotalLowE ()

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 Types
typedef std::map< const G4ParticleDefinition , G4LowEXsection , std::less< const G4ParticleDefinition * > >	LowEMap

Private Member Functions
	G4XNNTotalLowE (const G4XNNTotalLowE &right)

const G4XNNTotalLowE &	operator= (const G4XNNTotalLowE &right)

Private Attributes
std::map< const G4ParticleDefinition , G4LowEXsection , std::less< const G4ParticleDefinition * > >	theCrossSections

Static Private Attributes
static const G4double	npTot [29]

static const G4double	ppTot [29]

static const G4double	ss [29]

static const G4int	tableSize = 29

Detailed Description

Definition at line 40 of file G4XNNTotalLowE.hh.

Member Typedef Documentation

◆ LowEMap

typedef std::map<const G4ParticleDefinition *, G4LowEXsection*, std::less<const G4ParticleDefinition *> > G4XNNTotalLowE::LowEMap

private

Definition at line 72 of file G4XNNTotalLowE.hh.

Constructor & Destructor Documentation

◆ G4XNNTotalLowE() [1/2]

G4XNNTotalLowE::G4XNNTotalLowE ( )

Definition at line 53 of file G4XNNTotalLowE.cc.

{ 
    
  G4LowEXsection* pp = new G4LowEXsection();
  G4LowEXsection* np = new G4LowEXsection();
 
  G4int i;
  for (i=0; i<tableSize; i++)
    {
      std::pair<double,double> it;
      it.first=ss[i];
      it.second=ppTot[i]; pp->push_back(it);
      it.second=npTot[i]; np->push_back(it);
    }
  theCrossSections[G4Proton::ProtonDefinition()] = pp;
  theCrossSections[G4Neutron::NeutronDefinition()] = np;
}

References G4Neutron::NeutronDefinition(), npTot, G4InuclParticleNames::pp, ppTot, G4Proton::ProtonDefinition(), ss, tableSize, and theCrossSections.

◆ ~G4XNNTotalLowE()

G4XNNTotalLowE::~G4XNNTotalLowE ( )

virtual

Definition at line 72 of file G4XNNTotalLowE.cc.

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

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

◆ G4XNNTotalLowE() [2/2]

G4XNNTotalLowE::G4XNNTotalLowE ( const G4XNNTotalLowE & right )

private

Member Function Documentation

◆ CrossSection()

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

virtual

Implements G4VCrossSectionSource.

Definition at line 79 of file G4XNNTotalLowE.cc.

{
  G4double result = 0.;
  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
  const G4ParticleDefinition * key = FindKeyParticle(trk1,trk2);
 
  if (theCrossSections.find(key)!= theCrossSections.end())
  {
    LowEMap::const_iterator iter;
    for (iter = theCrossSections.begin(); iter != theCrossSections.end(); ++iter)
    {
      if ((*iter).first == key)
      {
    result = (*iter).second->CrossSection(sqrtS);
      }
    }
  }
  else
  {
    throw G4HadronicException(__FILE__, __LINE__, "G4XNNTotalLowE: particle key out of range");
  }
 
  return result;
}

References G4VCrossSectionSource::FindKeyParticle(), G4KineticTrack::Get4Momentum(), and theCrossSections.

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

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 50 of file G4XNNTotalLowE.hh.

50{ return 0; }

◆ HighLimit()

virtual G4double G4XNNTotalLowE::HighLimit ( ) const

inlinevirtual

Reimplemented from G4VCrossSectionSource.

Definition at line 54 of file G4XNNTotalLowE.hh.

54{ return 3.*CLHEP::GeV; }

CLHEP::GeV

static constexpr double GeV

Definition: SystemOfUnits.h:184

References CLHEP::GeV.

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

◆ IsValid()

G4bool G4XNNTotalLowE::IsValid ( G4double e ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 111 of file G4XNNTotalLowE.cc.

{
  G4bool result = e>0&&e<3*GeV;
  return result;
}

References GeV.

◆ 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 G4XNNTotalLowE::Name ( ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 104 of file G4XNNTotalLowE.cc.

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

References G4InuclParticleNames::name().

◆ operator!=()

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

inherited

Definition at line 80 of file G4VCrossSectionSource.cc.

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

◆ operator=()

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

private

◆ operator==()

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

inherited

Definition at line 74 of file G4VCrossSectionSource.cc.

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

◆ Print()

void G4VCrossSectionSource::Print ( ) const

virtualinherited

Reimplemented in G4XNNElasticLowE, G4XnpElasticLowE, and G4XnpTotalLowE.

Definition at line 86 of file G4VCrossSectionSource.cc.

{
  G4int nComponents = 0;
  const G4CrossSectionVector* components = GetComponents();
  if (components)
    {
      nComponents = components->size();
    }
  G4cout << "---- " << this->Name() << " ---- has " << nComponents << " components" <<G4endl;
  G4int i;
  for (i=0; i<nComponents; i++)
    {
      G4cout << "-" <<  this->Name() << " - Component " << i << ": " <<G4endl;
 
      G4CrossSectionSourcePtr componentPtr = (*components)[i];
      G4VCrossSectionSource* component = componentPtr();
      component->Print();
    }
}

References G4cout, G4endl, G4VCrossSectionSource::GetComponents(), G4VCrossSectionSource::Name(), and G4VCrossSectionSource::Print().

Referenced by G4VCollision::Print(), G4VCrossSectionSource::Print(), G4XNNElasticLowE::Print(), G4XnpElasticLowE::Print(), and G4XnpTotalLowE::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

◆ npTot

const G4double G4XNNTotalLowE::npTot

staticprivate

Initial value:

= 
{
4250, 1380, 770, 585, 465, 300, 232, 175, 140, 100, 80, 63, 50, 40, 35, 34, 34, 36.5, 37., 37.7, 38, 39, 39.8, 40.5, 40.7, 41, 41.2, 41.5, 41.5
}

Definition at line 67 of file G4XNNTotalLowE.hh.

Referenced by G4XNNTotalLowE().

◆ ppTot

const G4double G4XNNTotalLowE::ppTot

staticprivate

Initial value:

= 
{ 
  2000, 600, 250, 180, 138, 92, 75, 57, 46, 35, 29.5, 25.5, 25, 24, 23.7, 25, 29, 39, 44, 47, 48, 48, 47.5, 47, 45.6, 45, 43.3, 42.9, 42.9
}

Definition at line 65 of file G4XNNTotalLowE.hh.

Referenced by G4XNNTotalLowE().

◆ ss

const G4double G4XNNTotalLowE::ss

staticprivate

Initial value:

=
{
 1877.05, 1879.58, 1882.69, 1884.8, 1886.95, 1891.89, 1895.27, 1899.9, 1904.67, 1913.36, 1921.18, 1933.56, 1949.82, 1990.11, 2025.14, 2059.03, 2101.89, 2166.47, 2201.01, 2236.36, 2289.27, 2377.17, 2426.86, 2500.18, 2602.91, 2733.62, 2925.49, 3002.71 
}

Definition at line 66 of file G4XNNTotalLowE.hh.

Referenced by G4XNNTotalLowE().

◆ tableSize

const G4int G4XNNTotalLowE::tableSize = 29

staticprivate

Definition at line 68 of file G4XNNTotalLowE.hh.

Referenced by G4XNNTotalLowE().

◆ theCrossSections

std::map<const G4ParticleDefinition *, G4LowEXsection *, std::less<const G4ParticleDefinition *> > G4XNNTotalLowE::theCrossSections

private

Definition at line 71 of file G4XNNTotalLowE.hh.

Referenced by CrossSection(), G4XNNTotalLowE(), and ~G4XNNTotalLowE().

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

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

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ LowEMap

Constructor & Destructor Documentation

◆ G4XNNTotalLowE() [1/2]

◆ ~G4XNNTotalLowE()

◆ G4XNNTotalLowE() [2/2]

Member Function Documentation

◆ CrossSection()

◆ FcrossX()

◆ FindKeyParticle()

◆ FindLightParticle()

◆ GetComponents()

◆ HighLimit()

◆ InLimits()

◆ IsValid()

◆ LowLimit()

◆ Name()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ Print()

◆ PrintAll()

Field Documentation

◆ npTot

◆ ppTot

◆ ss

◆ tableSize

◆ theCrossSections