G4XPDGTotal Class Reference

#include <G4XPDGTotal.hh>

Inheritance diagram for G4XPDGTotal:

Public Member Functions
	G4XPDGTotal ()
virtual	~G4XPDGTotal ()
G4bool	operator== (const G4XPDGTotal &right) const
G4bool	operator!= (const G4XPDGTotal &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 55 of file G4XPDGTotal.hh.

Constructor & Destructor Documentation

G4XPDGTotal()

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;
}

~G4XPDGTotal()

G4XPDGTotal::~G4XPDGTotal ( )

virtual

Definition at line 142 of file G4XPDGTotal.cc.

{ }

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;
}

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

GetComponents()

virtual const G4CrossSectionVector * G4XPDGTotal::GetComponents ( ) const

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 69 of file G4XPDGTotal.hh.

{ return 0; }

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;
}

LowLimit()

virtual G4double G4XPDGTotal::LowLimit ( ) const

inlinevirtual

Reimplemented from G4VCrossSectionSource.

Definition at line 75 of file G4XPDGTotal.hh.

{ return _lowLimit; }

Name()

G4String G4XPDGTotal::Name ( ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 237 of file G4XPDGTotal.cc.

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

operator!=()

G4bool G4XPDGTotal::operator!=

(

const G4XPDGTotal &

right

)

const

Definition at line 152 of file G4XPDGTotal.cc.

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

operator==()

G4bool G4XPDGTotal::operator==

(

const G4XPDGTotal &

right

)

const

Definition at line 146 of file G4XPDGTotal.cc.

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

---

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

• G4XPDGTotal.hh
• G4XPDGTotal.cc