#include <G4XPDGTotal.hh>

Inheritance diagram for G4XPDGTotal:

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

	G4XPDGTotal ()

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

G4bool	operator== (const G4VCrossSectionSource &right) const

G4bool	operator== (const G4XPDGTotal &right) const

virtual void	Print () const

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

virtual	~G4XPDGTotal ()

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< G4pDefPair, std::vector< G4double >, std::less< G4pDefPair > >	PairDoubleMap

Private Member Functions
	G4XPDGTotal (const G4XPDGTotal &right)

const G4XPDGTotal &	operator= (const G4XPDGTotal &right)

G4double	PDGTotal (G4double coeff, G4double ecm) const

Private Attributes
std::map< G4pDefPair, std::vector< G4double >, std::less< G4pDefPair > >	xMap

Static Private Attributes
static const G4double	_highLimit = DBL_MAX

static const G4double	_lowLimit = 3. * GeV

static const G4double	gammagammaPDGFit [5] = { 3., 300., 0.000156, 0.00032, 0. }

static const G4double	gammapPDGFit [5] = { 3., 300., 0.0577, 0.1171, 0. }

static const G4double	KnPDGFit [5] = { 3., 40., 10.376, 14.29, 7.38 }

static const G4double	KpPDGFit [5] = { 3., 40., 10.376, 15.57, 13.19 }

static const G4int	nFit = 5

static const G4double	npPDGFit [5] = { 3., 40., 18.256, 61.14, 29.80 }

static const G4double	pipPDGFit [5] = { 3., 40., 11.568, 27.55, 5.62 }

static const G4double	ppPDGFit [5] = { 3., 40000., 18.256, 60.19, 33.43 }

Detailed Description

Definition at line 55 of file G4XPDGTotal.hh.

Member Typedef Documentation

◆ PairDoubleMap

typedef std::map<G4pDefPair, std::vector<G4double>, std::less<G4pDefPair> > G4XPDGTotal::PairDoubleMap

private

Definition at line 100 of file G4XPDGTotal.hh.

Constructor & Destructor Documentation

◆ G4XPDGTotal() [1/2]

G4XPDGTotal::G4XPDGTotal ( )

Definition at line 70 of file G4XPDGTotal.cc.

{
  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> pp(G4Proton::ProtonDefinition(),
                    G4Proton::ProtonDefinition());
  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> pn(G4Proton::ProtonDefinition(),
                    G4Neutron::NeutronDefinition());
  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> piPlusp(G4PionPlus::PionPlusDefinition(),
                     G4Proton::ProtonDefinition());
  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> piMinusp(G4PionMinus::PionMinusDefinition(),
                      G4Proton::ProtonDefinition());
  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> KPlusp(G4KaonPlus::KaonPlusDefinition(),
                    G4Proton::ProtonDefinition());
  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> KPlusn(G4KaonPlus::KaonPlusDefinition(),
                    G4Neutron::NeutronDefinition());
  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> KMinusp(G4KaonMinus::KaonMinusDefinition(),
                     G4Proton::ProtonDefinition());
  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> KMinusn(G4KaonMinus::KaonMinusDefinition(),
                     G4Neutron::NeutronDefinition());
  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> gp(G4Gamma::GammaDefinition(),
                    G4Proton::ProtonDefinition());
  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> gg(G4Gamma::GammaDefinition(),
                    G4Gamma::GammaDefinition());
  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> nn(G4Neutron::NeutronDefinition(),
                    G4Neutron::NeutronDefinition());
 
  std::vector<G4double> nnData;
  std::vector<G4double> ppData;
  std::vector<G4double> pnData;
  std::vector<G4double> pipData;
  std::vector<G4double> KpData;
  std::vector<G4double> KnData;
  std::vector<G4double> gpData;
  std::vector<G4double> ggData;
 
  G4int i;
  for (i=0; i<2; i++) 
    {   
      nnData.push_back(ppPDGFit[i] * GeV); 
      ppData.push_back(ppPDGFit[i] * GeV); 
      pnData.push_back(npPDGFit[i] * GeV); 
      pipData.push_back(pipPDGFit[i] * GeV); 
      KpData.push_back(KpPDGFit[i] * GeV); 
      KnData.push_back(KnPDGFit[i] * GeV); 
      gpData.push_back(gammapPDGFit[i] * GeV); 
      ggData.push_back(gammagammaPDGFit[i] * GeV); 
    }
  for (i=2; i<nFit; i++) 
    {   
      nnData.push_back(ppPDGFit[i]); 
      ppData.push_back(ppPDGFit[i]); 
      pnData.push_back(npPDGFit[i]); 
      pipData.push_back(pipPDGFit[i]); 
      KpData.push_back(KpPDGFit[i]); 
      KnData.push_back(KnPDGFit[i]); 
      gpData.push_back(gammapPDGFit[i]); 
      ggData.push_back(gammagammaPDGFit[i]); 
    }
 
  xMap[pp] = ppData;
  xMap[pn] = pnData;
  xMap[piPlusp] = pipData;
  xMap[piMinusp] = pipData;
  xMap[KPlusp] = KpData;
  xMap[KPlusn] = KnData;
  xMap[KMinusp] = KpData;
  xMap[KMinusn] = KnData;
  xMap[gp] = gpData;
  xMap[gg] = ggData;
  xMap[nn] = nnData;
}

References G4Gamma::GammaDefinition(), gammagammaPDGFit, gammapPDGFit, GeV, G4KaonMinus::KaonMinusDefinition(), G4KaonPlus::KaonPlusDefinition(), KnPDGFit, KpPDGFit, G4Neutron::NeutronDefinition(), nFit, G4InuclParticleNames::nn, npPDGFit, G4PionMinus::PionMinusDefinition(), G4PionPlus::PionPlusDefinition(), pipPDGFit, G4InuclParticleNames::pn, G4InuclParticleNames::pp, ppPDGFit, G4Proton::ProtonDefinition(), and xMap.

◆ ~G4XPDGTotal()

G4XPDGTotal::~G4XPDGTotal ( )

virtual

Definition at line 142 of file G4XPDGTotal.cc.

143{ }

◆ G4XPDGTotal() [2/2]

G4XPDGTotal::G4XPDGTotal ( const G4XPDGTotal & right )

private

Member Function Documentation

◆ CrossSection()

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

virtual

Implements G4VCrossSectionSource.

Definition at line 158 of file G4XPDGTotal.cc.

{
  G4double sigma = 0.;
 
  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
  
  const G4ParticleDefinition* def1 = trk1.GetDefinition();
  const G4ParticleDefinition* def2 = trk2.GetDefinition();
 
  G4double enc1 = def1->GetPDGEncoding();
  G4double enc2 = def2->GetPDGEncoding();
  G4double coeff = -1.;
  if ( (enc1 < 0 && enc2 >0) || (enc2 < 0 && enc1 >0) ) coeff = 1.;
 
  // Order the pair: first is the lower mass particle, second is the higher mass one
  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> trkPair(def1,def2);
 
  if (def1->GetPDGMass() > def2->GetPDGMass())
    trkPair = std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *>(def2,def1);
 
  std::vector<G4double> data;   
      
  if (xMap.find(trkPair) != xMap.end())
    {
 
      PairDoubleMap::const_iterator iter;
      for (iter = xMap.begin(); iter != xMap.end(); ++iter)
    {
      std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> thePair = (*iter).first;
      if (thePair == trkPair)
        {
          data = (*iter).second;
      
          G4double eMinFit = data[0];
          G4double eMaxFit = data[1];
          G4double xFit = data[2];
          G4double y1Fit = data[3];
          G4double y2Fit = data[4];
          
          // Total Cross-section fit, 1998 Review of Particle Properties, European Phys. J. 3(1998), 1
          
          // Parameters from the PDG fit
          static const G4double epsilon = 0.095;
          static const G4double eta1 = -0.34;
          static const G4double eta2 = -0.55;
          
          if (sqrtS < eMinFit || sqrtS > eMaxFit)
          {
              G4cout << "WARNING! G4XPDGTotal::PDGTotal extrapolating cross section at " 
               << sqrtS / GeV 
               << " GeV outside the PDG fit range "
               << eMinFit / GeV << " - " << eMaxFit / GeV << " GeV " << G4endl;
          }
          
          G4double S = (sqrtS * sqrtS) / (GeV*GeV);
          
          sigma = ( (xFit * G4Pow::GetInstance()->powA(S,epsilon)) + 
            (y1Fit * G4Pow::GetInstance()->powA(S,eta1)) + 
            (coeff * y2Fit * G4Pow::GetInstance()->powA(S,eta2)) ) * millibarn;
          
          if (sigma < 0.)
        {
          G4String name1 = def1->GetParticleName();
          G4String name2 = def2->GetParticleName();
          G4cout << "WARNING! G4XPDGTotal::PDGTotal "      
             << name1 << "-" << name2 
             << " total cross section: Ecm " 
             << sqrtS / GeV << " GeV, negative cross section " 
             << sigma / millibarn << " mb set to 0" << G4endl;
          sigma = 0.;
        }
        }
    }
    }
  return sigma;
}

References epsilon(), G4cout, G4endl, G4KineticTrack::Get4Momentum(), G4KineticTrack::GetDefinition(), G4Pow::GetInstance(), G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetPDGEncoding(), G4ParticleDefinition::GetPDGMass(), GeV, millibarn, G4Pow::powA(), S(), and xMap.

Referenced by G4XpimNTotal::CrossSection(), and G4XpipNTotal::CrossSection().

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

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 69 of file G4XPDGTotal.hh.

69{ 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 IsValid().

◆ IsValid()

G4bool G4XPDGTotal::IsValid ( G4double e ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 244 of file G4XPDGTotal.cc.

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

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

◆ LowLimit()

virtual G4double G4XPDGTotal::LowLimit ( ) const

inlinevirtual

Reimplemented from G4VCrossSectionSource.

Definition at line 75 of file G4XPDGTotal.hh.

75{ return _lowLimit; }

References _lowLimit.

◆ Name()

G4String G4XPDGTotal::Name ( ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 237 of file G4XPDGTotal.cc.

{
  G4String name = "PDGTotal ";
  return name;
}

References G4InuclParticleNames::name().

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

Definition at line 152 of file G4XPDGTotal.cc.

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

◆ operator=()

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

Definition at line 146 of file G4XPDGTotal.cc.

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

◆ PDGTotal()

G4double G4XPDGTotal::PDGTotal	(	G4double	coeff,
		G4double	ecm
	)		const

private

Definition at line 252 of file G4XPDGTotal.cc.

{
  return 0.;
}

◆ 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

◆ _highLimit

const G4double G4XPDGTotal::_highLimit = DBL_MAX

staticprivate

Definition at line 88 of file G4XPDGTotal.hh.

Referenced by IsValid().

◆ _lowLimit

const G4double G4XPDGTotal::_lowLimit = 3. * GeV

staticprivate

Definition at line 87 of file G4XPDGTotal.hh.

Referenced by IsValid(), and LowLimit().

◆ gammagammaPDGFit

const G4double G4XPDGTotal::gammagammaPDGFit = { 3., 300., 0.000156, 0.00032, 0. }

staticprivate

Definition at line 97 of file G4XPDGTotal.hh.

Referenced by G4XPDGTotal().

◆ gammapPDGFit

const G4double G4XPDGTotal::gammapPDGFit = { 3., 300., 0.0577, 0.1171, 0. }

staticprivate

Definition at line 96 of file G4XPDGTotal.hh.

Referenced by G4XPDGTotal().

◆ KnPDGFit

const G4double G4XPDGTotal::KnPDGFit = { 3., 40., 10.376, 14.29, 7.38 }

staticprivate

Definition at line 95 of file G4XPDGTotal.hh.

Referenced by G4XPDGTotal().

◆ KpPDGFit

const G4double G4XPDGTotal::KpPDGFit = { 3., 40., 10.376, 15.57, 13.19 }

staticprivate

Definition at line 94 of file G4XPDGTotal.hh.

Referenced by G4XPDGTotal().

◆ nFit

const G4int G4XPDGTotal::nFit = 5

staticprivate

Definition at line 90 of file G4XPDGTotal.hh.

Referenced by G4XPDGTotal().

◆ npPDGFit

const G4double G4XPDGTotal::npPDGFit = { 3., 40., 18.256, 61.14, 29.80 }

staticprivate

Definition at line 92 of file G4XPDGTotal.hh.

Referenced by G4XPDGTotal().

◆ pipPDGFit

const G4double G4XPDGTotal::pipPDGFit = { 3., 40., 11.568, 27.55, 5.62 }

staticprivate

Definition at line 93 of file G4XPDGTotal.hh.

Referenced by G4XPDGTotal().

◆ ppPDGFit

const G4double G4XPDGTotal::ppPDGFit = { 3., 40000., 18.256, 60.19, 33.43 }

staticprivate

Definition at line 91 of file G4XPDGTotal.hh.

Referenced by G4XPDGTotal().

◆ xMap

std::map<G4pDefPair, std::vector<G4double>, std::less<G4pDefPair> > G4XPDGTotal::xMap

private

Definition at line 99 of file G4XPDGTotal.hh.

Referenced by CrossSection(), and G4XPDGTotal().

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

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

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ PairDoubleMap

Constructor & Destructor Documentation

◆ G4XPDGTotal() [1/2]

◆ ~G4XPDGTotal()

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

◆ PDGTotal()

◆ Print()

◆ PrintAll()

Field Documentation

◆ _highLimit

◆ _lowLimit

◆ gammagammaPDGFit

◆ gammapPDGFit

◆ KnPDGFit

◆ KpPDGFit

◆ nFit

◆ npPDGFit

◆ pipPDGFit

◆ ppPDGFit

◆ xMap