Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4ResonanceNames.cc
Go to the documentation of this file.
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26#include "globals.hh"
27#include "G4ResonanceNames.hh"
28#include "G4VCrossSectionSource.hh"
29#include "G4KineticTrack.hh"
30#include "G4CrossSectionVector.hh"
31#include "G4ParticleDefinition.hh"
32#include "G4ShortLivedConstructor.hh"
33#include "G4ShortLivedTable.hh"
34#include "G4ParticleTable.hh"
35#include "G4PionPlus.hh"
36
37G4ResonanceNames::G4ResonanceNames()
38{
39 // Excited N resonances
40
41 nameNstar.push_back("N(1440)+");
42 nameNstar.push_back("N(1440)0");
43 nameNstar.push_back("N(1520)+");
44 nameNstar.push_back("N(1520)0");
45 nameNstar.push_back("N(1535)+");
46 nameNstar.push_back("N(1535)0");
47 nameNstar.push_back("N(1650)+");
48 nameNstar.push_back("N(1650)0");
49 nameNstar.push_back("N(1675)+");
50 nameNstar.push_back("N(1675)0");
51 nameNstar.push_back("N(1680)+");
52 nameNstar.push_back("N(1680)0");
53 nameNstar.push_back("N(1700)+");
54 nameNstar.push_back("N(1700)0");
55 nameNstar.push_back("N(1710)+");
56 nameNstar.push_back("N(1710)0");
57 nameNstar.push_back("N(1720)+");
58 nameNstar.push_back("N(1720)0");
59 nameNstar.push_back("N(1900)+");
60 nameNstar.push_back("N(1900)0");
61 nameNstar.push_back("N(1990)+");
62 nameNstar.push_back("N(1990)0");
63 nameNstar.push_back("N(2090)+");
64 nameNstar.push_back("N(2090)0");
65 nameNstar.push_back("N(2190)+");
66 nameNstar.push_back("N(2190)0");
67 nameNstar.push_back("N(2220)+");
68 nameNstar.push_back("N(2220)0");
69 nameNstar.push_back("N(2250)+");
70 nameNstar.push_back("N(2250)0");
71
72
73 // Delta
74
75 G4String d1232Minus("delta-");
76 G4String d1232Zero("delta0");
77 G4String d1232Plus("delta+");
78 G4String d1232PlusPlus("delta++");
79 nameDelta.push_back(d1232Minus);
80 nameDelta.push_back(d1232Zero);
81 nameDelta.push_back(d1232Plus);
82 nameDelta.push_back(d1232PlusPlus);
83
84
85 // Excited Delta resonances
86
87 nameDeltastar.push_back("delta(1600)+");
88 nameDeltastar.push_back("delta(1600)++");
89 nameDeltastar.push_back("delta(1600)-");
90 nameDeltastar.push_back("delta(1600)0");
91 nameDeltastar.push_back("delta(1620)+");
92 nameDeltastar.push_back("delta(1620)++");
93 nameDeltastar.push_back("delta(1620)-");
94 nameDeltastar.push_back("delta(1620)0");
95 nameDeltastar.push_back("delta(1700)+");
96 nameDeltastar.push_back("delta(1700)++");
97 nameDeltastar.push_back("delta(1700)-");
98 nameDeltastar.push_back("delta(1700)0");
99 nameDeltastar.push_back("delta(1900)+");
100 nameDeltastar.push_back("delta(1900)++");
101 nameDeltastar.push_back("delta(1900)-");
102 nameDeltastar.push_back("delta(1900)0");
103 nameDeltastar.push_back("delta(1905)+");
104 nameDeltastar.push_back("delta(1905)++");
105 nameDeltastar.push_back("delta(1905)-");
106 nameDeltastar.push_back("delta(1905)0");
107 nameDeltastar.push_back("delta(1910)+");
108 nameDeltastar.push_back("delta(1910)++");
109 nameDeltastar.push_back("delta(1910)-");
110 nameDeltastar.push_back("delta(1910)0");
111 nameDeltastar.push_back("delta(1920)+");
112 nameDeltastar.push_back("delta(1920)++");
113 nameDeltastar.push_back("delta(1920)-");
114 nameDeltastar.push_back("delta(1920)0");
115 nameDeltastar.push_back("delta(1930)+");
116 nameDeltastar.push_back("delta(1930)++");
117 nameDeltastar.push_back("delta(1930)-");
118 nameDeltastar.push_back("delta(1930)0");
119 nameDeltastar.push_back("delta(1950)+");
120 nameDeltastar.push_back("delta(1950)++");
121 nameDeltastar.push_back("delta(1950)-");
122 nameDeltastar.push_back("delta(1950)0");
123
124
125 // Lambda
126
127 nameLambda.push_back("lambda");
128 nameLambda.push_back("lambda(1405)");
129 nameLambda.push_back("lambda(1520)");
130 nameLambda.push_back("lambda(1600)");
131 nameLambda.push_back("lambda(1670)");
132 nameLambda.push_back("lambda(1690)");
133 nameLambda.push_back("lambda(1800)");
134 nameLambda.push_back("lambda(1810)");
135 nameLambda.push_back("lambda(1820)");
136 nameLambda.push_back("lambda(1830)");
137 nameLambda.push_back("lambda(1890)");
138 nameLambda.push_back("lambda(2100)");
139 nameLambda.push_back("lambda(2110)");
140
141
142 // Sigma
143
144 nameSigma.push_back("sigma(1385)+");
145 nameSigma.push_back("sigma(1385)-");
146 nameSigma.push_back("sigma(1385)0");
147 nameSigma.push_back("sigma(1660)+");
148 nameSigma.push_back("sigma(1660)-");
149 nameSigma.push_back("sigma(1660)0");
150 nameSigma.push_back("sigma(1670)+");
151 nameSigma.push_back("sigma(1670)-");
152 nameSigma.push_back("sigma(1670)0");
153 nameSigma.push_back("sigma(1750)+");
154 nameSigma.push_back("sigma(1750)-");
155 nameSigma.push_back("sigma(1750)0");
156 nameSigma.push_back("sigma(1775)+");
157 nameSigma.push_back("sigma(1775)-");
158 nameSigma.push_back("sigma(1775)0");
159 nameSigma.push_back("sigma(1915)+");
160 nameSigma.push_back("sigma(1915)-");
161 nameSigma.push_back("sigma(1915)0");
162 nameSigma.push_back("sigma(1940)+");
163 nameSigma.push_back("sigma(1940)-");
164 nameSigma.push_back("sigma(1940)0");
165 nameSigma.push_back("sigma(2030)+");
166 nameSigma.push_back("sigma(2030)-");
167 nameSigma.push_back("sigma(2030)0");
168
169
170 // Xi
171
172 nameXi.push_back("xi(1530)-");
173 nameXi.push_back("xi(1530)0");
174 nameXi.push_back("xi(1690)-");
175 nameXi.push_back("xi(1690)0");
176 nameXi.push_back("xi(1820)-");
177 nameXi.push_back("xi(1820)0");
178 nameXi.push_back("xi(1950)-");
179 nameXi.push_back("xi(1950)0");
180 nameXi.push_back("xi(2030)-");
181 nameXi.push_back("xi(2030)0");
182
183
184 G4ShortLivedConstructor ShortLived;
185 ShortLived.ConstructParticle();
186 G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
187
188 size_t i;
189
190 // Fill a map with the lowest resonance for each category
191 for (i=0; i<nameNstar.size(); i++)
192 {
193 lowResMap[nameNstar[i]] = particleTable->FindParticle("N(1440)0");
194 }
195
196 for (i=0; i<nameDeltastar.size(); i++)
197 {
198 lowResMap[nameDeltastar[i]] = particleTable->FindParticle("delta0");
199 }
200
201 for (i=0; i<nameDelta.size(); i++)
202 {
203 lowResMap[nameDelta[i]] = particleTable->FindParticle("delta0");
204 }
205
206 for (i=0; i<nameLambda.size(); i++)
207 {
208 lowResMap[nameLambda[i]] = particleTable->FindParticle("lambda");
209 }
210
211 for (i=0; i<nameSigma.size(); i++)
212 {
213 lowResMap[nameSigma[i]] = particleTable->FindParticle("sigma0");
214 }
215
216 shortMap["N(1440)0"] = "N(1440)";
217 shortMap["N(1440)+"] = "N(1440)";
218
219 shortMap["N(1520)0"] = "N(1520)";;
220 shortMap["N(1520)+"] = "N(1520)";
221
222 shortMap["N(1535)0"] = "N(1535)";
223 shortMap["N(1535)+"] = "N(1535)";
224
225 shortMap["N(1650)0"] = "N(1650)";
226 shortMap["N(1650)+"] = "N(1650)";
227
228 shortMap["N(1675)0"] = "N(1675)";
229 shortMap["N(1675)+"] = "N(1675)";
230
231 shortMap["N(1680)0"] = "N(1680)";
232 shortMap["N(1680)+"] = "N(1680)";
233
234 shortMap["N(1700)0"] = "N(1700)";
235 shortMap["N(1700)+"] = "N(1700)";
236
237 shortMap["N(1710)0"] = "N(1710)";
238 shortMap["N(1710)+"] = "N(1710)";
239
240 shortMap["N(1720)0"] = "N(1720)";
241 shortMap["N(1720)+"] = "N(1720)";
242
243 shortMap["N(1900)0"] = "N(1900)";
244 shortMap["N(1900)+"] = "N(1900)";
245
246 shortMap["N(1990)0"] = "N(1990)";
247 shortMap["N(1990)+"] = "N(1990)";
248
249 shortMap["N(2090)0"] = "N(2090)";
250 shortMap["N(2090)+"] = "N(2090)";
251
252 shortMap["N(2190)0"] = "N(2190)";
253 shortMap["N(2190)+"] = "N(2190)";
254
255 shortMap["N(2220)0"] = "N(2220)";
256 shortMap["N(2220)+"] = "N(2220)";
257
258 shortMap["N(2250)0"] = "N(2250)";
259 shortMap["N(2250)+"] = "N(2250)";
260
261
262 // Excited Delta
263
264 shortMap["delta(1600)-"] = "delta(1600)";
265 shortMap["delta(1600)0"] = "delta(1600)";
266 shortMap["delta(1600)+"] = "delta(1600)";
267 shortMap["delta(1600)++"] = "delta(1600)";
268
269 shortMap["delta(1620)-"] = "delta(1620)";
270 shortMap["delta(1620)0"] = "delta(1620)";
271 shortMap["delta(1620)+"] = "delta(1620)";
272 shortMap["delta(1620)++"] = "delta(1620)";
273
274 shortMap["delta(1700)-"] = "delta(1700)";
275 shortMap["delta(1700)0"] = "delta(1700)";
276 shortMap["delta(1700)+"] = "delta(1700)";
277 shortMap["delta(1700)++"] = "delta(1700)";
278
279 shortMap["delta(1900)-"] = "delta(1900)";
280 shortMap["delta(1900)0"] = "delta(1900)";
281 shortMap["delta(1900)+"] = "delta(1900)";
282 shortMap["delta(1900)++"] = "delta(1900)";
283
284 shortMap["delta(1905)-"] = "delta(1905)";
285 shortMap["delta(1905)0"] = "delta(1905)";
286 shortMap["delta(1905)+"] = "delta(1905)";
287 shortMap["delta(1905)++"] = "delta(1905)";
288
289 shortMap["delta(1910)-"] = "delta(1910)";
290 shortMap["delta(1910)0"] = "delta(1910)";
291 shortMap["delta(1910)+"] = "delta(1910)";
292 shortMap["delta(1910)++"] = "delta(1910)";
293
294 shortMap["delta(1920)-"] = "delta(1920)";
295 shortMap["delta(1920)0"] = "delta(1920)";
296 shortMap["delta(1920)+"] = "delta(1920)";
297 shortMap["delta(1920)++"] = "delta(1920)";
298
299 shortMap["delta(1930)-"] = "delta(1930)";
300 shortMap["delta(1930)0"] = "delta(1930)";
301 shortMap["delta(1930)+"] = "delta(1930)";
302 shortMap["delta(1930)++"] = "delta(1930)";
303
304 shortMap["delta(1950)-"] = "delta(1950)";
305 shortMap["delta(1950)0"] = "delta(1950)";
306 shortMap["delta(1950)+"] = "delta(1950)";
307 shortMap["delta(1950)++"] = "delta(1950)";
308
309 // Delta
310
311 shortMap["delta-"] = "delta";
312 shortMap["delta0"] = "delta";
313 shortMap["delta+"] = "delta";
314 shortMap["delta++"] = "delta";
315
316}
317
318
319G4ResonanceNames::~G4ResonanceNames()
320{ }
321
322
323G4bool G4ResonanceNames::operator==(const G4ResonanceNames &right) const
324{
325 return(this == (G4ResonanceNames*) &right);
326}
327
328
329G4bool G4ResonanceNames::operator!=(const G4ResonanceNames &right) const
330{
331 return (this != (G4ResonanceNames*) &right);
332}
333
334
335G4double G4ResonanceNames::MinMass(const G4String& name)
336{
337 // Cut, from UrQMD (reference still to be detailed)
338 const G4double coeff = 0.001;
339
340 G4double lowMass = 0.;
341
342 G4ParticleDefinition* def = 0;
343
344 if (lowResMap.find(name) != lowResMap.end())
345 {
346 def = lowResMap[name];
347 }
348 else
349 {
350 G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
351 def = particleTable->FindParticle(name);
352 }
353 if (def != 0)
354 {
355 G4ParticleDefinition* pion = G4PionPlus::PionPlusDefinition();
356 lowMass = (1. + coeff) * def->GetPDGMass() + pion->GetPDGMass();
357 }
358 else
359 {
360 G4cout << "G4ResonanceNames::MinMass - " << name << " not found" << G4endl;
361 throw G4HadronicException(__FILE__, __LINE__, "G4ResonanceNames::MinMass - resonance name not found");
362 }
363 return lowMass;
364}
365
366
367const G4String G4ResonanceNames::ShortName(const G4String& name)
368{
369 G4String shortName = "";
370 if (shortMap.find(name) != shortMap.end())
371 {
372 shortName = shortMap[name];
373 }
374 return shortName;
375}
G4double
double G4double
Definition: G4Types.hh:64
G4bool
bool G4bool
Definition: G4Types.hh:67
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cout
G4DLLIMPORT std::ostream G4cout
G4ParticleTable::FindParticle
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
Definition: G4ParticleTable.cc:472
G4ParticleTable::GetParticleTable
static G4ParticleTable * GetParticleTable()
Definition: G4ParticleTable.cc:63
G4PionPlus::PionPlusDefinition
static G4PionPlus * PionPlusDefinition()
Definition: G4PionPlus.cc:93
G4ResonanceNames::operator==
G4bool operator==(const G4ResonanceNames &right) const
Definition: G4ResonanceNames.cc:323
G4ResonanceNames::MinMass
G4double MinMass(const G4String &name)
Definition: G4ResonanceNames.cc:335
G4ResonanceNames::operator!=
G4bool operator!=(const G4ResonanceNames &right) const
Definition: G4ResonanceNames.cc:329
G4ResonanceNames::ShortName
const G4String ShortName(const G4String &name)
Definition: G4ResonanceNames.cc:367