G4XpipNTotal Class Reference

#include <G4XpipNTotal.hh>

Inheritance diagram for G4XpipNTotal:

Public Member Functions
virtual G4double	CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
	G4XpipNTotal ()
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	~G4XpipNTotal ()

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 Attributes
std::vector< std::pair< double, double > >	theLowEData
G4XPDGTotal	thePDGData

Detailed Description

Definition at line 35 of file G4XpipNTotal.hh.

Constructor & Destructor Documentation

G4XpipNTotal()

G4XpipNTotal::G4XpipNTotal

(

)

Definition at line 32 of file G4XpipNTotal.cc.

 {
   std::pair<double,double> it;
   it.first=1105.46; it.second= 6.4;  theLowEData.push_back(it);
   it.first=1133.28; it.second= 18.0; theLowEData.push_back(it);
   it.first=1165.12; it.second= 66.0; theLowEData.push_back(it);
   it.first=1198.8;  it.second= 160.; theLowEData.push_back(it);
   it.first=1212.51; it.second= 199.; theLowEData.push_back(it);
   it.first=1233.18; it.second= 195.; theLowEData.push_back(it);
   it.first=1301.86; it.second= 73.8; theLowEData.push_back(it);
   it.first=1368.93; it.second= 35.0; theLowEData.push_back(it);
   it.first=1433.81; it.second= 22.0; theLowEData.push_back(it);
   it.first=1496.42; it.second= 15.0; theLowEData.push_back(it);
   it.first=1556.84; it.second= 15.2; theLowEData.push_back(it);
   it.first=1615.21; it.second= 19.4; theLowEData.push_back(it);
   it.first=1671.7;  it.second= 25.;  theLowEData.push_back(it);
   it.first=1726.44; it.second= 26.;  theLowEData.push_back(it);
   it.first=1779.57; it.second= 30.;  theLowEData.push_back(it);
   it.first=1881.49; it.second= 40.;  theLowEData.push_back(it);
   it.first=1930.49; it.second= 40.;  theLowEData.push_back(it);
   it.first=2070.69; it.second= 30.;  theLowEData.push_back(it);
   it.first=2202.11; it.second= 29.;  theLowEData.push_back(it);
   it.first=2326.19; it.second= 30.7; theLowEData.push_back(it);
   it.first=2405.38; it.second= 30.7; theLowEData.push_back(it);
   it.first=2733.67; it.second= 28.3; theLowEData.push_back(it);
   it.first=3207.21; it.second= 26.5; theLowEData.push_back(it);
   it.first=4434.76; it.second= 25.0; theLowEData.push_back(it);
 } 

References theLowEData.

~G4XpipNTotal()

virtual G4XpipNTotal::~G4XpipNTotal ( )

inlinevirtual

Definition at line 39 of file G4XpipNTotal.hh.

{}

Member Function Documentation

CrossSection()

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

virtual

Implements G4VCrossSectionSource.

Definition at line 62 of file G4XpipNTotal.cc.

 {
   G4double sqrts = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
   if(sqrts > theLowEData.back().first/MeV) return thePDGData.CrossSection(trk1, trk2);
   G4double result = 0;
   size_t i(0), it(0);
   if(sqrts<theLowEData[0].first) return 0;
   for(i=0; i<theLowEData.size(); i++)
   {
     if(theLowEData[i].first/MeV>sqrts) break;
     it = i;
   }
   G4double x1 = G4Log(theLowEData[it].first);
   G4double x2 = G4Log(theLowEData[it+1].first);
   G4double y1 = G4Log(theLowEData[it].second);
   G4double y2 = G4Log(theLowEData[it+1].second);
   G4double x = G4Log(sqrts);
   G4double y = y1+(x-x1)*(y2-y1)/(x2-x1);
   result = G4Exp(y);
   return result*millibarn;
 }

References G4XPDGTotal::CrossSection(), G4Exp(), G4Log(), G4KineticTrack::Get4Momentum(), MeV, millibarn, second, theLowEData, and thePDGData.

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

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 41 of file G4XpipNTotal.hh.

{ return 0; }

HighLimit()

G4double G4VCrossSectionSource::HighLimit ( ) const

virtualinherited

Reimplemented in G4XNNElasticLowE, G4XNNTotalLowE, G4XnpElasticLowE, and G4XnpTotalLowE.

Definition at line 146 of file G4VCrossSectionSource.cc.

{
  return DBL_MAX; 
}

References DBL_MAX.

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

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 G4VCrossSectionSource::IsValid ( G4double  e ) const

virtualinherited

Reimplemented in G4CrossSectionComposite, G4CrossSectionPatch, G4XAnnihilationChannel, G4XAqmElastic, G4XAqmTotal, G4XMesonBaryonElastic, G4XNNElasticLowE, G4XNNTotalLowE, G4XnpElasticLowE, G4XnpTotalLowE, G4XPDGElastic, and G4XPDGTotal.

Definition at line 151 of file G4VCrossSectionSource.cc.

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

References G4VCrossSectionSource::HighLimit(), and G4VCrossSectionSource::LowLimit().

Referenced by G4CrossSectionComposite::CrossSection(), G4CrossSectionPatch::CrossSection(), G4CrossSectionComposite::IsValid(), and G4CrossSectionPatch::IsValid().

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()

virtual G4String G4XpipNTotal::Name ( ) const

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 42 of file G4XpipNTotal.hh.

{return "G4XpipNTotal";}

operator!=()

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

inherited

Definition at line 80 of file G4VCrossSectionSource.cc.

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

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

theLowEData

std::vector<std::pair<double,double> > G4XpipNTotal::theLowEData

private

Definition at line 45 of file G4XpipNTotal.hh.

Referenced by CrossSection(), and G4XpipNTotal().

thePDGData

G4XPDGTotal G4XpipNTotal::thePDGData

private

Definition at line 44 of file G4XpipNTotal.hh.

Referenced by CrossSection().

---

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

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