55const G4double G4XPDGTotal::ppPDGFit[5] = { 3., 40000., 18.256, 60.19, 33.43 };
57const G4double G4XPDGTotal::npPDGFit[5] = { 3., 40., 18.256, 61.14, 29.80 };
59const G4double G4XPDGTotal::pipPDGFit[5] = { 3., 40., 11.568, 27.55, 5.62 };
61const G4double G4XPDGTotal::KpPDGFit[5] = { 3., 40., 10.376, 15.57, 13.19 };
63const G4double G4XPDGTotal::KnPDGFit[5] = { 3., 40., 10.376, 14.29, 7.38 };
65const G4double G4XPDGTotal::gammapPDGFit[5] = { 3., 300., 0.0577, 0.1171, 0. };
67const G4double G4XPDGTotal::gammagammaPDGFit[5] = { 3., 300., 0.000156, 0.00032, 0. };
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;
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);
116 for (i=2; i<nFit; i++)
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]);
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;
213 G4double S = (sqrtS * sqrtS) / (GeV*GeV);
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;
G4double GetPDGMass() const
G4int GetPDGEncoding() const
const G4String & GetParticleName() const