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;
108 ppData.push_back(
ppPDGFit[i] * GeV);
109 pnData.push_back(
npPDGFit[i] * GeV);
111 KpData.push_back(
KpPDGFit[i] * GeV);
112 KnData.push_back(
KnPDGFit[i] * GeV);
116 for (i=2; i<
nFit; i++)
130 xMap[piPlusp] = pipData;
131 xMap[piMinusp] = pipData;
132 xMap[KPlusp] = KpData;
133 xMap[KPlusn] = KnData;
134 xMap[KMinusp] = KpData;
135 xMap[KMinusn] = KnData;
171 if ( (enc1 < 0 && enc2 >0) || (enc2 < 0 && enc1 >0) ) coeff = 1.;
174 std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> trkPair(def1,def2);
177 trkPair = std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *>(def2,def1);
179 std::vector<G4double>
data;
181 if (
xMap.find(trkPair) !=
xMap.end())
184 PairDoubleMap::const_iterator iter;
185 for (iter =
xMap.begin(); iter !=
xMap.end(); ++iter)
187 std::pair<const G4ParticleDefinition *,const G4ParticleDefinition *> thePair = (*iter).first;
188 if (thePair == trkPair)
190 data = (*iter).second;
205 if (sqrtS < eMinFit || sqrtS > eMaxFit)
207 G4cout <<
"WARNING! G4XPDGTotal::PDGTotal extrapolating cross section at "
209 <<
" GeV outside the PDG fit range "
210 << eMinFit /
GeV <<
" - " << eMaxFit /
GeV <<
" GeV " <<
G4endl;
223 G4cout <<
"WARNING! G4XPDGTotal::PDGTotal "
224 << name1 <<
"-" << name2
225 <<
" total cross section: Ecm "
226 << sqrtS / GeV <<
" GeV, negative cross section "
G4bool operator==(const G4XPDGTotal &right) const
static const G4double ppPDGFit[5]
G4bool operator!=(const G4XPDGTotal &right) const
const G4LorentzVector & Get4Momentum() const
static const G4double KnPDGFit[5]
static const G4double npPDGFit[5]
virtual G4bool IsValid(G4double e) const
const G4String & GetParticleName() const
const G4ParticleDefinition * GetDefinition() const
static G4PionPlus * PionPlusDefinition()
G4double GetPDGMass() const
virtual G4double CrossSection(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
const XML_Char const XML_Char * data
virtual G4String Name() const
static constexpr double millibarn
static G4Pow * GetInstance()
static const G4double gammagammaPDGFit[5]
static G4Proton * ProtonDefinition()
static const G4double gammapPDGFit[5]
G4double powA(G4double A, G4double y) const
static G4PionMinus * PionMinusDefinition()
static const G4double _lowLimit
double epsilon(double density, double temperature)
static G4KaonMinus * KaonMinusDefinition()
static G4Gamma * GammaDefinition()
static const G4double KpPDGFit[5]
std::map< G4pDefPair, std::vector< G4double >, std::less< G4pDefPair > > xMap
static const G4double _highLimit
G4double PDGTotal(G4double coeff, G4double ecm) const
static G4Neutron * NeutronDefinition()
G4int GetPDGEncoding() const
G4GLOB_DLL std::ostream G4cout
static G4KaonPlus * KaonPlusDefinition()
static const G4double pipPDGFit[5]
static constexpr double GeV
G4bool InLimits(G4double e, G4double eLow, G4double eHigh) const