#include <G4XPDGElastic.hh>

Inheritance diagram for G4XPDGElastic:

Collaboration diagram for G4XPDGElastic:

Public Member Functions
	G4XPDGElastic ()

virtual	~G4XPDGElastic ()

G4bool	operator== (const G4XPDGElastic &right) const

G4bool	operator!= (const G4XPDGElastic &right) const

virtual G4double	CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

virtual const G4CrossSectionVector *	GetComponents () const

virtual G4String	Name () const

virtual G4bool	IsValid (G4double e) const

virtual G4double	LowLimit () const

Public Member Functions inherited from G4VCrossSectionSource
	G4VCrossSectionSource ()

virtual	~G4VCrossSectionSource ()

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 G4double	HighLimit () const

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

Private Member Functions
	G4XPDGElastic (const G4XPDGElastic &right)

const G4XPDGElastic &	operator= (const G4XPDGElastic &right)

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

Static Private Attributes
static const G4double	_lowLimit = 5. * GeV

static const G4double	_highLimit = DBL_MAX

static const G4int	nPar = 7

static const G4double	pPiPlusPDGFit [7] = { 2., 200., 0., 11.4, -0.4, 0.079, 0. }

static const G4double	pPiMinusPDGFit [7] = { 2., 360., 1.76, 11.2, -0.64, 0.043, 0. }

static const G4double	pKPlusPDGFit [7] = { 2., 175., 5.0, 8.1, -1.8, 0.16, -1.3 }

static const G4double	pKMinusPDGFit [7] = { 2., 175., 7.3, 0., 0., 0.29, -2.40 }

static const G4double	ppPDGFit [7] = { 2., 2100., 11.9, 26.9, -1.21, 0.169, -1.85 }

static const G4double	ppbarPDGFit [7] = { 5., 1730000., 10.2, 52.7, -1.16, 0.125, -1.28 }

static const G4double	npbarPDGFit [7] = { 1.1, 5.55, 36.5, 0., 0., 0., -11.9 }

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionSource
G4bool	InLimits (G4double e, G4double eLow, G4double eHigh) const

const G4ParticleDefinition *	FindLightParticle (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

G4double	FcrossX (G4double e, G4double e0, G4double sigma, G4double eParam, G4double power) const

const G4ParticleDefinition *	FindKeyParticle (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

Detailed Description

Definition at line 55 of file G4XPDGElastic.hh.

Member Typedef Documentation

◆ PairDoubleMap

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

private

Definition at line 100 of file G4XPDGElastic.hh.

Constructor & Destructor Documentation

◆ G4XPDGElastic() [1/2]

G4XPDGElastic::G4XPDGElastic ( )

Definition at line 77 of file G4XPDGElastic.cc.

 {
   const G4ParticleDefinition * proton = G4Proton::ProtonDefinition();
   const G4ParticleDefinition * neutron = G4Neutron::NeutronDefinition();
   const G4ParticleDefinition * piPlus = G4PionPlus::PionPlusDefinition();
   const G4ParticleDefinition * piMinus = G4PionMinus::PionMinusDefinition();
   const G4ParticleDefinition * KPlus = G4KaonPlus::KaonPlusDefinition();
   const G4ParticleDefinition * KMinus = G4KaonMinus::KaonMinusDefinition();
   const G4ParticleDefinition * antiproton = G4AntiProton::AntiProtonDefinition();
   
   std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> pp(proton,proton);
   std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> pn(proton,neutron);
   std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> piPlusp(piPlus,proton);
   std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> piMinusp(piMinus,proton);
   std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> KPlusp(KPlus,proton);
   std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> KMinusp(KMinus,proton);
   std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> nn(neutron,neutron);
   std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> ppbar(proton,antiproton);
   std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> npbar(antiproton,neutron);
 
   std::vector<G4double> ppData;
   std::vector<G4double> pPiPlusData;
   std::vector<G4double> pPiMinusData;
   std::vector<G4double> pKPlusData;
   std::vector<G4double> pKMinusData;
   std::vector<G4double> ppbarData;
   std::vector<G4double> npbarData;
 
   G4int i;
   for (i=0; i<2; i++) 
     {      
       ppData.push_back(ppPDGFit[i] * GeV); 
       pPiPlusData.push_back(pPiPlusPDGFit[i] * GeV); 
       pPiMinusData.push_back(pPiMinusPDGFit[i] * GeV); 
       pKPlusData.push_back(pKPlusPDGFit[i] * GeV); 
       pKMinusData.push_back(pKMinusPDGFit[i] * GeV); 
       ppbarData.push_back(ppbarPDGFit[i] * GeV);
       npbarData.push_back(npbarPDGFit[i] * GeV);
     }
 
   for (i=2; i<nPar; i++) 
     {      
       ppData.push_back(ppPDGFit[i]); 
       pPiPlusData.push_back(pPiPlusPDGFit[i]); 
       pPiMinusData.push_back(pPiMinusPDGFit[i]); 
       pKPlusData.push_back(pKPlusPDGFit[i]); 
       pKMinusData.push_back(pKMinusPDGFit[i]); 
       ppbarData.push_back(ppbarPDGFit[i]);
       npbarData.push_back(npbarPDGFit[i]);
     }
 
   xMap[nn] = ppData;
   xMap[pp] = ppData;
   xMap[pn] = ppData;
   xMap[piPlusp] = pPiPlusData;
   xMap[piMinusp] = pPiMinusData;
   xMap[KPlusp] = pKPlusData;
   xMap[KMinusp] = pKMinusData;
   xMap[ppbar] = ppbarData;
   xMap[npbar] = npbarData;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ~G4XPDGElastic()

G4XPDGElastic::~G4XPDGElastic ( )

virtual

Definition at line 140 of file G4XPDGElastic.cc.

141 { }

◆ G4XPDGElastic() [2/2]

G4XPDGElastic::G4XPDGElastic ( const G4XPDGElastic & right )

private

Member Function Documentation

◆ CrossSection()

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

virtual

Implements G4VCrossSectionSource.

Definition at line 156 of file G4XPDGElastic.cc.

 {
   // Elastic Cross-section fit, 1994 Review of Particle Properties, (1994), 1
 
   G4double sigma = 0.;
 
   const G4ParticleDefinition* def1 = trk1.GetDefinition();
   const G4ParticleDefinition* def2 = trk2.GetDefinition();
   
   G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
   G4double m_1 = def1->GetPDGMass();
   G4double m_2 = def2->GetPDGMass();
   G4double m_max = std::max(m_1,m_2);
   //  if (m1 > m) m = m1;
   
   G4double pLab = 0.;
 
   if (m_max > 0. && sqrtS > (m_1 + m_2)) 
     {
       pLab = std::sqrt( (sqrtS*sqrtS - (m_1+m_2)*(m_1+m_2) ) * (sqrtS*sqrtS - (m_1-m_2)*(m_1-m_2)) ) / (2*m_max);
       
       // The PDG fit formula requires p in GeV/c
       
       // 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; 
       G4double pMinFit = 0.;
       G4double pMaxFit = 0.;
       G4double aFit = 0.;
       G4double bFit = 0.;
       G4double cFit = 0.;
       G4double dFit = 0.;
       G4double nFit = 0.;
 
       // Debug
 //      G4cout << "Map has " << xMap.size() << " elements" << G4endl;
       // Debug end
  
       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;
           pMinFit = data[0];
           pMaxFit = data[1];
           aFit = data[2];
           bFit = data[3];
           cFit = data[5];
           dFit = data[6];
           nFit = data[4];
           
           if (pLab < pMinFit) return 0.0;
           if (pLab > pMaxFit )
             G4cout << "WARNING! G4XPDGElastic::PDGElastic " 
                << trk1.GetDefinition()->GetParticleName() << "-" 
                << trk2.GetDefinition()->GetParticleName() 
                << " elastic cross section: momentum " 
                << pLab / GeV << " GeV outside valid fit range " 
                << pMinFit /GeV << "-" << pMaxFit / GeV
                << G4endl;
           
           pLab /= GeV;
           if (pLab > 0.) 
             {
               G4double logP = G4Log(pLab);
               sigma = aFit + bFit * G4Pow::GetInstance()->powA(pLab, nFit) + cFit * logP*logP + dFit * logP;
               sigma = sigma * millibarn;
             }
         }
         }
     }
       else
     {
       G4cout << "G4XPDGElastic::CrossSection - Data not found in Map" << G4endl;
     }
     }
   
   if (sigma < 0.)
     {
       G4cout << "WARNING! G4XPDGElastic::PDGElastic "      
          << def1->GetParticleName() << "-" << def2->GetParticleName()
          << " elastic cross section: momentum " 
          << pLab << " GeV, negative cross section " 
          << sigma / millibarn << " mb set to 0" << G4endl;
       sigma = 0.;
     }
   
   return sigma;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetComponents()

virtual const G4CrossSectionVector* G4XPDGElastic::GetComponents ( ) const

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 69 of file G4XPDGElastic.hh.

69 { return 0; }

Here is the call graph for this function:

◆ IsValid()

G4bool G4XPDGElastic::IsValid ( G4double e ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 261 of file G4XPDGElastic.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

◆ LowLimit()

virtual G4double G4XPDGElastic::LowLimit ( ) const

inlinevirtual

Reimplemented from G4VCrossSectionSource.

Definition at line 75 of file G4XPDGElastic.hh.

75 { return _lowLimit; }

G4XPDGElastic::_lowLimit

static const G4double _lowLimit

Definition: G4XPDGElastic.hh:86

Here is the call graph for this function:

◆ Name()

G4String G4XPDGElastic::Name ( ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 254 of file G4XPDGElastic.cc.

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

Here is the caller graph for this function:

◆ operator!=()

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

Definition at line 150 of file G4XPDGElastic.cc.

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

◆ operator=()

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

private

Here is the caller graph for this function:

◆ operator==()

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

Definition at line 144 of file G4XPDGElastic.cc.

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

Member Data Documentation

◆ _highLimit

const G4double G4XPDGElastic::_highLimit = DBL_MAX

staticprivate

Definition at line 87 of file G4XPDGElastic.hh.

◆ _lowLimit

const G4double G4XPDGElastic::_lowLimit = 5. * GeV

staticprivate

Definition at line 86 of file G4XPDGElastic.hh.

◆ nPar

const G4int G4XPDGElastic::nPar = 7

staticprivate

Definition at line 89 of file G4XPDGElastic.hh.

◆ npbarPDGFit

const G4double G4XPDGElastic::npbarPDGFit = { 1.1, 5.55, 36.5, 0., 0., 0., -11.9 }

staticprivate

Definition at line 96 of file G4XPDGElastic.hh.

◆ pKMinusPDGFit

const G4double G4XPDGElastic::pKMinusPDGFit = { 2., 175., 7.3, 0., 0., 0.29, -2.40 }

staticprivate

Definition at line 93 of file G4XPDGElastic.hh.

◆ pKPlusPDGFit

const G4double G4XPDGElastic::pKPlusPDGFit = { 2., 175., 5.0, 8.1, -1.8, 0.16, -1.3 }

staticprivate

Definition at line 92 of file G4XPDGElastic.hh.

◆ ppbarPDGFit

const G4double G4XPDGElastic::ppbarPDGFit = { 5., 1730000., 10.2, 52.7, -1.16, 0.125, -1.28 }

staticprivate

Definition at line 95 of file G4XPDGElastic.hh.

◆ pPiMinusPDGFit

const G4double G4XPDGElastic::pPiMinusPDGFit = { 2., 360., 1.76, 11.2, -0.64, 0.043, 0. }

staticprivate

Definition at line 91 of file G4XPDGElastic.hh.

◆ pPiPlusPDGFit

const G4double G4XPDGElastic::pPiPlusPDGFit = { 2., 200., 0., 11.4, -0.4, 0.079, 0. }

staticprivate

Definition at line 90 of file G4XPDGElastic.hh.

◆ ppPDGFit

const G4double G4XPDGElastic::ppPDGFit = { 2., 2100., 11.9, 26.9, -1.21, 0.169, -1.85 }

staticprivate

Definition at line 94 of file G4XPDGElastic.hh.

◆ xMap

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

private

Definition at line 98 of file G4XPDGElastic.hh.

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

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ PairDoubleMap

Constructor & Destructor Documentation

◆ G4XPDGElastic() [1/2]

◆ ~G4XPDGElastic()

◆ G4XPDGElastic() [2/2]

Member Function Documentation

◆ CrossSection()

◆ GetComponents()

◆ IsValid()

◆ LowLimit()

◆ Name()

◆ operator!=()

◆ operator=()

◆ operator==()

Member Data Documentation

◆ _highLimit

◆ _lowLimit

◆ nPar

◆ npbarPDGFit

◆ pKMinusPDGFit

◆ pKPlusPDGFit

◆ ppbarPDGFit

◆ pPiMinusPDGFit

◆ pPiPlusPDGFit

◆ ppPDGFit

◆ xMap