Geant4  10.02.p03
Public Member Functions | Private Types | Private Member Functions | Private Attributes | Static Private Attributes | List of all members
G4XPDGTotal Class Reference

#include <G4XPDGTotal.hh>

Inheritance diagram for G4XPDGTotal:
Collaboration 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 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

G4double PDGTotal (G4double coeff, G4double ecm) const
 
 G4XPDGTotal (const G4XPDGTotal &right)
 
const G4XPDGTotal & operator= (const G4XPDGTotal &right)
 

Private Attributes

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

Static Private Attributes

static const G4double _lowLimit = 3. * GeV
 
static const G4double _highLimit = DBL_MAX
 
static const G4int nFit = 5
 
static const G4double ppPDGFit [5] = { 3., 40000., 18.256, 60.19, 33.43 }
 
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 KpPDGFit [5] = { 3., 40., 10.376, 15.57, 13.19 }
 
static const G4double KnPDGFit [5] = { 3., 40., 10.376, 14.29, 7.38 }
 
static const G4double gammapPDGFit [5] = { 3., 300., 0.0577, 0.1171, 0. }
 
static const G4double gammagammaPDGFit [5] = { 3., 300., 0.000156, 0.00032, 0. }
 

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.

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.

71 {
72  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> pp(G4Proton::ProtonDefinition(),
73  G4Proton::ProtonDefinition());
74  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> pn(G4Proton::ProtonDefinition(),
75  G4Neutron::NeutronDefinition());
76  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> piPlusp(G4PionPlus::PionPlusDefinition(),
77  G4Proton::ProtonDefinition());
78  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> piMinusp(G4PionMinus::PionMinusDefinition(),
79  G4Proton::ProtonDefinition());
80  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> KPlusp(G4KaonPlus::KaonPlusDefinition(),
81  G4Proton::ProtonDefinition());
82  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> KPlusn(G4KaonPlus::KaonPlusDefinition(),
83  G4Neutron::NeutronDefinition());
84  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> KMinusp(G4KaonMinus::KaonMinusDefinition(),
85  G4Proton::ProtonDefinition());
86  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> KMinusn(G4KaonMinus::KaonMinusDefinition(),
87  G4Neutron::NeutronDefinition());
88  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> gp(G4Gamma::GammaDefinition(),
89  G4Proton::ProtonDefinition());
90  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> gg(G4Gamma::GammaDefinition(),
91  G4Gamma::GammaDefinition());
92  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> nn(G4Neutron::NeutronDefinition(),
93  G4Neutron::NeutronDefinition());
94 
95  std::vector<G4double> nnData;
96  std::vector<G4double> ppData;
97  std::vector<G4double> pnData;
98  std::vector<G4double> pipData;
99  std::vector<G4double> KpData;
100  std::vector<G4double> KnData;
101  std::vector<G4double> gpData;
102  std::vector<G4double> ggData;
103 
104  G4int i;
105  for (i=0; i<2; i++)
106  {
107  nnData.push_back(ppPDGFit[i] * GeV);
108  ppData.push_back(ppPDGFit[i] * GeV);
109  pnData.push_back(npPDGFit[i] * GeV);
110  pipData.push_back(pipPDGFit[i] * GeV);
111  KpData.push_back(KpPDGFit[i] * GeV);
112  KnData.push_back(KnPDGFit[i] * GeV);
113  gpData.push_back(gammapPDGFit[i] * GeV);
114  ggData.push_back(gammagammaPDGFit[i] * GeV);
115  }
116  for (i=2; i<nFit; i++)
117  {
118  nnData.push_back(ppPDGFit[i]);
119  ppData.push_back(ppPDGFit[i]);
120  pnData.push_back(npPDGFit[i]);
121  pipData.push_back(pipPDGFit[i]);
122  KpData.push_back(KpPDGFit[i]);
123  KnData.push_back(KnPDGFit[i]);
124  gpData.push_back(gammapPDGFit[i]);
125  ggData.push_back(gammagammaPDGFit[i]);
126  }
127 
128  xMap[pp] = ppData;
129  xMap[pn] = pnData;
130  xMap[piPlusp] = pipData;
131  xMap[piMinusp] = pipData;
132  xMap[KPlusp] = KpData;
133  xMap[KPlusn] = KnData;
134  xMap[KMinusp] = KpData;
135  xMap[KMinusn] = KnData;
136  xMap[gp] = gpData;
137  xMap[gg] = ggData;
138  xMap[nn] = nnData;
139 }
G4XPDGTotal::nFit
static const G4int nFit
Definition: G4XPDGTotal.hh:90
G4KaonPlus::KaonPlusDefinition
static G4KaonPlus * KaonPlusDefinition()
Definition: G4KaonPlus.cc:108
G4XPDGTotal::npPDGFit
static const G4double npPDGFit[5]
Definition: G4XPDGTotal.hh:92
G4XPDGTotal::KpPDGFit
static const G4double KpPDGFit[5]
Definition: G4XPDGTotal.hh:94
G4KaonMinus::KaonMinusDefinition
static G4KaonMinus * KaonMinusDefinition()
Definition: G4KaonMinus.cc:108
G4Proton::ProtonDefinition
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
G4XPDGTotal::KnPDGFit
static const G4double KnPDGFit[5]
Definition: G4XPDGTotal.hh:95
G4int
int G4int
Definition: G4Types.hh:78
G4PionMinus::PionMinusDefinition
static G4PionMinus * PionMinusDefinition()
Definition: G4PionMinus.cc:93
G4PionPlus::PionPlusDefinition
static G4PionPlus * PionPlusDefinition()
Definition: G4PionPlus.cc:93
GeV
static const double GeV
Definition: G4SIunits.hh:214
G4XPDGTotal::gammagammaPDGFit
static const G4double gammagammaPDGFit[5]
Definition: G4XPDGTotal.hh:97
G4XPDGTotal::ppPDGFit
static const G4double ppPDGFit[5]
Definition: G4XPDGTotal.hh:91
G4XPDGTotal::gammapPDGFit
static const G4double gammapPDGFit[5]
Definition: G4XPDGTotal.hh:96
G4XPDGTotal::pipPDGFit
static const G4double pipPDGFit[5]
Definition: G4XPDGTotal.hh:93
G4XPDGTotal::xMap
std::map< G4pDefPair, std::vector< G4double >, std::less< G4pDefPair > > xMap
Definition: G4XPDGTotal.hh:99
G4Neutron::NeutronDefinition
static G4Neutron * NeutronDefinition()
Definition: G4Neutron.cc:99
G4Gamma::GammaDefinition
static G4Gamma * GammaDefinition()
Definition: G4Gamma.cc:81
Here is the call graph for this function:
Here is the caller graph for this function:

◆ ~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.

160 {
161  G4double sigma = 0.;
162 
163  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
164 
165  const G4ParticleDefinition* def1 = trk1.GetDefinition();
166  const G4ParticleDefinition* def2 = trk2.GetDefinition();
167 
168  G4double enc1 = def1->GetPDGEncoding();
169  G4double enc2 = def2->GetPDGEncoding();
170  G4double coeff = -1.;
171  if ( (enc1 < 0 && enc2 >0) || (enc2 < 0 && enc1 >0) ) coeff = 1.;
172 
173  // Order the pair: first is the lower mass particle, second is the higher mass one
174  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> trkPair(def1,def2);
175 
176  if (def1->GetPDGMass() > def2->GetPDGMass())
177  trkPair = std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *>(def2,def1);
178 
179  std::vector<G4double> data;
180 
181  if (xMap.find(trkPair) != xMap.end())
182  {
183 
184  PairDoubleMap::const_iterator iter;
185  for (iter = xMap.begin(); iter != xMap.end(); ++iter)
186  {
187  std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> thePair = (*iter).first;
188  if (thePair == trkPair)
189  {
190  data = (*iter).second;
191 
192  G4double eMinFit = data[0];
193  G4double eMaxFit = data[1];
194  G4double xFit = data[2];
195  G4double y1Fit = data[3];
196  G4double y2Fit = data[4];
197 
198  // Total Cross-section fit, 1998 Review of Particle Properties, European Phys. J. 3(1998), 1
199 
200  // Parameters from the PDG fit
201  static const G4double epsilon = 0.095;
202  static const G4double eta1 = -0.34;
203  static const G4double eta2 = -0.55;
204 
205  if (sqrtS < eMinFit || sqrtS > eMaxFit)
206  {
207  G4cout << "WARNING! G4XPDGTotal::PDGTotal extrapolating cross section at "
208  << sqrtS / GeV
209  << " GeV outside the PDG fit range "
210  << eMinFit / GeV << " - " << eMaxFit / GeV << " GeV " << G4endl;
211  }
212 
213  G4double S = (sqrtS * sqrtS) / (GeV*GeV);
214 
215  sigma = ( (xFit * G4Pow::GetInstance()->powA(S,epsilon)) +
216  (y1Fit * G4Pow::GetInstance()->powA(S,eta1)) +
217  (coeff * y2Fit * G4Pow::GetInstance()->powA(S,eta2)) ) * millibarn;
218 
219  if (sigma < 0.)
220  {
221  G4String name1 = def1->GetParticleName();
222  G4String name2 = def2->GetParticleName();
223  G4cout << "WARNING! G4XPDGTotal::PDGTotal "
224  << name1 << "-" << name2
225  << " total cross section: Ecm "
226  << sqrtS / GeV << " GeV, negative cross section "
227  << sigma / millibarn << " mb set to 0" << G4endl;
228  sigma = 0.;
229  }
230  }
231  }
232  }
233  return sigma;
234 }
G4KineticTrack::GetDefinition
const G4ParticleDefinition * GetDefinition() const
Definition: G4KineticTrack.hh:210
G4Pow::GetInstance
static G4Pow * GetInstance()
Definition: G4Pow.cc:55
S
double S(double temp)
Definition: G4DNAElectronHoleRecombination.cc:76
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:120
G4cout
G4GLOB_DLL std::ostream G4cout
GeV
static const double GeV
Definition: G4SIunits.hh:214
millibarn
static const double millibarn
Definition: G4SIunits.hh:105
G4endl
#define G4endl
Definition: G4ios.hh:61
G4Pow::powA
G4double powA(G4double A, G4double y) const
Definition: G4Pow.hh:259
G4double
double G4double
Definition: G4Types.hh:76
G4XPDGTotal::xMap
std::map< G4pDefPair, std::vector< G4double >, std::less< G4pDefPair > > xMap
Definition: G4XPDGTotal.hh:99
epsilon
double epsilon(double density, double temperature)
Definition: G4DNAElectronHoleRecombination.cc:95
G4KineticTrack::Get4Momentum
const G4LorentzVector & Get4Momentum() const
Definition: G4KineticTrack.hh:240
Here is the call graph for this function:
Here is the caller graph for this function:

◆ GetComponents()

virtual const G4CrossSectionVector* G4XPDGTotal::GetComponents ( ) const
inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 69 of file G4XPDGTotal.hh.

69 { return 0; }
Here is the call graph for this function:

◆ IsValid()

G4bool G4XPDGTotal::IsValid ( G4double  e) const
virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 244 of file G4XPDGTotal.cc.

245 {
246  G4bool answer = InLimits(e,_lowLimit,_highLimit);
247 
248  return answer;
249 }
G4XPDGTotal::_highLimit
static const G4double _highLimit
Definition: G4XPDGTotal.hh:88
G4bool
bool G4bool
Definition: G4Types.hh:79
G4VCrossSectionSource::InLimits
G4bool InLimits(G4double e, G4double eLow, G4double eHigh) const
Definition: G4VCrossSectionSource.cc:134
G4XPDGTotal::_lowLimit
static const G4double _lowLimit
Definition: G4XPDGTotal.hh:87
Here is the call graph for this function:
Here is the caller graph for this function:

◆ LowLimit()

virtual G4double G4XPDGTotal::LowLimit ( ) const
inlinevirtual

Reimplemented from G4VCrossSectionSource.

Definition at line 75 of file G4XPDGTotal.hh.

75 { return _lowLimit; }
G4XPDGTotal::_lowLimit
static const G4double _lowLimit
Definition: G4XPDGTotal.hh:87
Here is the call graph for this function:

◆ Name()

G4String G4XPDGTotal::Name ( ) const
virtual

Implements G4VCrossSectionSource.

Definition at line 237 of file G4XPDGTotal.cc.

238 {
239  G4String name = "PDGTotal ";
240  return name;
241 }
name
G4String name
Definition: TRTMaterials.hh:40
Here is the caller graph for this function:

◆ operator!=()

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

Definition at line 152 of file G4XPDGTotal.cc.

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

◆ operator=()

const G4XPDGTotal& G4XPDGTotal::operator= ( const G4XPDGTotal &  right)
private
Here is the caller graph for this function:

◆ operator==()

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

Definition at line 146 of file G4XPDGTotal.cc.

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

◆ PDGTotal()

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

Definition at line 252 of file G4XPDGTotal.cc.

253 {
254  return 0.;
255 }
Here is the caller graph for this function:

Member Data Documentation

◆ _highLimit

const G4double G4XPDGTotal::_highLimit = DBL_MAX
staticprivate

Definition at line 88 of file G4XPDGTotal.hh.

◆ _lowLimit

const G4double G4XPDGTotal::_lowLimit = 3. * GeV
staticprivate

Definition at line 87 of file G4XPDGTotal.hh.

◆ gammagammaPDGFit

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

Definition at line 97 of file G4XPDGTotal.hh.

◆ gammapPDGFit

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

Definition at line 96 of file G4XPDGTotal.hh.

◆ KnPDGFit

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

Definition at line 95 of file G4XPDGTotal.hh.

◆ KpPDGFit

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

Definition at line 94 of file G4XPDGTotal.hh.

◆ nFit

const G4int G4XPDGTotal::nFit = 5
staticprivate

Definition at line 90 of file G4XPDGTotal.hh.

◆ npPDGFit

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

Definition at line 92 of file G4XPDGTotal.hh.

◆ pipPDGFit

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

Definition at line 93 of file G4XPDGTotal.hh.

◆ ppPDGFit

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

Definition at line 91 of file G4XPDGTotal.hh.

◆ xMap

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

Definition at line 99 of file G4XPDGTotal.hh.

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