52 const G4int G4XPDGTotal::nFit = 5;
54 const G4double G4XPDGTotal::ppPDGFit[5] = { 3., 40000., 18.256, 60.19, 33.43 };
56 const G4double G4XPDGTotal::npPDGFit[5] = { 3., 40., 18.256, 61.14, 29.80 };
58 const G4double G4XPDGTotal::pipPDGFit[5] = { 3., 40., 11.568, 27.55, 5.62 };
60 const G4double G4XPDGTotal::KpPDGFit[5] = { 3., 40., 10.376, 15.57, 13.19 };
62 const G4double G4XPDGTotal::KnPDGFit[5] = { 3., 40., 10.376, 14.29, 7.38 };
64 const G4double G4XPDGTotal::gammapPDGFit[5] = { 3., 300., 0.0577, 0.1171, 0. };
66 const G4double G4XPDGTotal::gammagammaPDGFit[5] = { 3., 300., 0.000156, 0.00032, 0. };
94 std::vector<G4double> nnData;
95 std::vector<G4double> ppData;
96 std::vector<G4double> pnData;
97 std::vector<G4double> pipData;
98 std::vector<G4double> KpData;
99 std::vector<G4double> KnData;
100 std::vector<G4double> gpData;
101 std::vector<G4double> ggData;
106 nnData.push_back(ppPDGFit[i] *
GeV);
107 ppData.push_back(ppPDGFit[i] * GeV);
108 pnData.push_back(npPDGFit[i] * GeV);
109 pipData.push_back(pipPDGFit[i] * GeV);
110 KpData.push_back(KpPDGFit[i] * GeV);
111 KnData.push_back(KnPDGFit[i] * GeV);
112 gpData.push_back(gammapPDGFit[i] * GeV);
113 ggData.push_back(gammagammaPDGFit[i] * GeV);
115 for (i=2; i<nFit; i++)
117 nnData.push_back(ppPDGFit[i]);
118 ppData.push_back(ppPDGFit[i]);
119 pnData.push_back(npPDGFit[i]);
120 pipData.push_back(pipPDGFit[i]);
121 KpData.push_back(KpPDGFit[i]);
122 KnData.push_back(KnPDGFit[i]);
123 gpData.push_back(gammapPDGFit[i]);
124 ggData.push_back(gammagammaPDGFit[i]);
129 xMap[piPlusp] = pipData;
130 xMap[piMinusp] = pipData;
131 xMap[KPlusp] = KpData;
132 xMap[KPlusn] = KnData;
133 xMap[KMinusp] = KpData;
134 xMap[KMinusn] = KnData;
170 if ( (enc1 < 0 && enc2 >0) || (enc2 < 0 && enc1 >0) ) coeff = 1.;
173 std::pair<G4ParticleDefinition *,G4ParticleDefinition *> trkPair(def1,def2);
176 trkPair = std::pair<G4ParticleDefinition *,G4ParticleDefinition *>(def2,def1);
178 std::vector<G4double>
data;
180 if (xMap.find(trkPair) != xMap.end())
183 PairDoubleMap::const_iterator iter;
184 for (iter = xMap.begin(); iter != xMap.end(); ++iter)
186 std::pair<G4ParticleDefinition *,G4ParticleDefinition *> thePair = (*iter).first;
187 if (thePair == trkPair)
189 data = (*iter).second;
204 if (sqrtS < eMinFit || sqrtS > eMaxFit)
206 G4cout <<
"WARNING! G4XPDGTotal::PDGTotal extrapolating cross section at "
208 <<
" GeV outside the PDG fit range "
209 << eMinFit /
GeV <<
" - " << eMaxFit /
GeV <<
" GeV " <<
G4endl;
214 sigma = ( (xFit * std::pow(S,epsilon)) +
215 (y1Fit * std::pow(S,eta1)) +
216 (coeff * y2Fit * std::pow(S,eta2)) ) *
millibarn;
222 G4cout <<
"WARNING! G4XPDGTotal::PDGTotal "
223 << name1 <<
"-" << name2
224 <<
" total cross section: Ecm "
225 << sqrtS / GeV <<
" GeV, negative cross section "
static G4KaonPlus * KaonPlusDefinition()
static G4KaonMinus * KaonMinusDefinition()
static G4Proton * ProtonDefinition()
G4int GetPDGEncoding() const
const G4String & GetParticleName() const
virtual G4bool IsValid(G4double e) const
G4ParticleDefinition * GetDefinition() const
G4GLOB_DLL std::ostream G4cout
static G4PionMinus * PionMinusDefinition()
static G4PionPlus * PionPlusDefinition()
G4bool operator==(const G4XPDGTotal &right) const
virtual G4String Name() const
G4double GetPDGMass() const
const G4LorentzVector & Get4Momentum() const
virtual G4double CrossSection(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
static G4Neutron * NeutronDefinition()
const XML_Char const XML_Char * data
G4bool InLimits(G4double e, G4double eLow, G4double eHigh) const
static G4Gamma * GammaDefinition()
G4bool operator!=(const G4XPDGTotal &right) const