G4XPDGElastic Class Reference

#include <G4XPDGElastic.hh>

Inheritance 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

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 54 of file G4XPDGElastic.hh.

Constructor & Destructor Documentation

G4XPDGElastic()

G4XPDGElastic::G4XPDGElastic

(

)

Definition at line 76 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;
}

~G4XPDGElastic()

G4XPDGElastic::~G4XPDGElastic ( )

virtual

Definition at line 139 of file G4XPDGElastic.cc.

{ }

Member Function Documentation

CrossSection()

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

virtual

Implements G4VCrossSectionSource.

Definition at line 155 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;
}

Referenced by G4XMesonBaryonElastic::CrossSection().

GetComponents()

virtual const G4CrossSectionVector * G4XPDGElastic::GetComponents ( ) const

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 68 of file G4XPDGElastic.hh.

{ return 0; }

IsValid()

G4bool G4XPDGElastic::IsValid ( G4double e ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 260 of file G4XPDGElastic.cc.

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

LowLimit()

virtual G4double G4XPDGElastic::LowLimit ( ) const

inlinevirtual

Reimplemented from G4VCrossSectionSource.

Definition at line 74 of file G4XPDGElastic.hh.

{ return _lowLimit; }

Name()

G4String G4XPDGElastic::Name ( ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 253 of file G4XPDGElastic.cc.

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

operator!=()

G4bool G4XPDGElastic::operator!=

(

const G4XPDGElastic &

right

)

const

Definition at line 149 of file G4XPDGElastic.cc.

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

operator==()

G4bool G4XPDGElastic::operator==

(

const G4XPDGElastic &

right

)

const

Definition at line 143 of file G4XPDGElastic.cc.

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

---

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

• G4XPDGElastic.hh
• G4XPDGElastic.cc