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 G4double	CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const =0
virtual const G4CrossSectionVector *	GetComponents () const =0
virtual G4String	Name () const =0
virtual void	Print () const
virtual void	PrintAll (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
virtual G4bool	IsValid (G4double e) const
virtual G4double	HighLimit () const
virtual G4double	LowLimit () 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