Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros
G4ExcitedSigmaConstructor.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//
27// $Id$
28//
29//
30// --------------------------------------------------------------
31// GEANT 4 class implementation file
32//
33// History: first implementation, based on object model of
34// 10 oct 1998 H.Kurashige
35// ---------------------------------------------------------------
36
37
38#include "G4ExcitedSigmaConstructor.hh"
39
40#include "G4SystemOfUnits.hh"
41#include "G4ParticleDefinition.hh"
42#include "G4ParticleTable.hh"
43#include "G4ShortLivedTable.hh"
44#include "G4PhaseSpaceDecayChannel.hh"
45#include "G4VDecayChannel.hh"
46#include "G4DecayTable.hh"
47
48G4ExcitedSigmaConstructor::G4ExcitedSigmaConstructor():
49 G4ExcitedBaryonConstructor(NStates, SigmaIsoSpin)
50{
51
52}
53
54G4ExcitedSigmaConstructor::~G4ExcitedSigmaConstructor()
55{
56}
57
58G4DecayTable* G4ExcitedSigmaConstructor::CreateDecayTable(
59 const G4String& parentName,
60 G4int iIso3,
61 G4int iState,
62 G4bool fAnti)
63{
64
65 // create decay table
66 G4DecayTable* decayTable = new G4DecayTable();
67
68 G4double br;
69 if ( (br=bRatio[iState][NK]) >0.0) {
70 AddNKMode( decayTable, parentName, br, iIso3, fAnti);
71 }
72
73 if ( (br=bRatio[iState][NKStar]) >0.0) {
74 AddNKStarMode( decayTable, parentName, br, iIso3, fAnti);
75 }
76
77 if ( (br=bRatio[iState][SigmaPi]) >0.0) {
78 AddSigmaPiMode( decayTable, parentName, br, iIso3, fAnti);
79 }
80
81 if ( (br=bRatio[iState][SigmaStarPi]) >0.0) {
82 AddSigmaStarPiMode( decayTable, parentName, br, iIso3, fAnti);
83 }
84
85 if ( (br=bRatio[iState][LambdaPi]) >0.0) {
86 AddLambdaPiMode( decayTable, parentName, br, iIso3, fAnti);
87 }
88
89 if ( (br=bRatio[iState][SigmaEta]) >0.0) {
90 AddSigmaEtaMode( decayTable, parentName, br, iIso3, fAnti);
91 }
92
93 if ( (br=bRatio[iState][LambdaStarPi]) >0.0) {
94 AddLambdaStarPiMode( decayTable, parentName, br, iIso3, fAnti);
95 }
96
97 if ( (br=bRatio[iState][DeltaK]) >0.0) {
98 AddDeltaKMode( decayTable, parentName, br, iIso3, fAnti);
99 }
100
101 return decayTable;
102}
103
104G4DecayTable* G4ExcitedSigmaConstructor::AddSigmaEtaMode(
105 G4DecayTable* decayTable, const G4String& nameParent,
106 G4double br, G4int iIso3, G4bool fAnti)
107{
108 G4VDecayChannel* mode;
109 //
110 G4String daughterH;
111 if (iIso3== +2) {
112 daughterH = "sigma+";
113 } else if (iIso3== 0) {
114 daughterH = "sigma0";
115 } else if (iIso3== -2) {
116 daughterH = "sigma-";
117 }
118 if (fAnti) daughterH = "anti_" + daughterH;
119
120 // create decay channel [parent BR #daughters]
121 mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
122 daughterH,"eta");
123 // add decay table
124 decayTable->Insert(mode);
125
126 return decayTable;
127}
128
129G4DecayTable* G4ExcitedSigmaConstructor::AddNKMode(
130 G4DecayTable* decayTable, const G4String& nameParent,
131 G4double br, G4int iIso3, G4bool fAnti)
132
133
134{
135 G4VDecayChannel* mode;
136
137 G4String daughterN;
138 G4String daughterK;
139 G4double r = 0.;
140
141 // ------------ N K- ------------
142 // determine daughters
143 if (iIso3== +2) {
144 r=0.;
145 } else if (iIso3== 0) {
146 daughterN = "proton";
147 r = br/2.;
148 } else if (iIso3== -2) {
149 daughterN = "neutron";
150 r = br;
151 }
152 if (!fAnti) {
153 daughterK = "kaon-";
154 } else {
155 daughterK = "kaon+";
156 }
157 if (fAnti) daughterN = "anti_" + daughterN;
158 if (r>0.) {
159 // create decay channel [parent BR #daughters]
160 mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
161 daughterN,daughterK);
162 // add decay table
163 decayTable->Insert(mode);
164 }
165
166 // ------------ N K0 ------------
167 // determine daughters
168 if (iIso3== +2) {
169 daughterN = "proton";
170 r=br;
171 } else if (iIso3== 0) {
172 daughterN = "neutron";
173 r = br/2.;
174 } else if (iIso3== -2) {
175 r = 0.;
176 }
177 if (!fAnti) {
178 daughterK = "anti_kaon0";
179 } else {
180
181 daughterK = "kaon0";
182 }
183 if (fAnti) daughterN = "anti_" + daughterN;
184 if (r>0.) {
185 // create decay channel [parent BR #daughters]
186 mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
187 daughterN,daughterK);
188 // add decay table
189 decayTable->Insert(mode);
190 }
191
192 return decayTable;
193}
194
195G4DecayTable* G4ExcitedSigmaConstructor::AddDeltaKMode(
196 G4DecayTable* decayTable, const G4String& nameParent,
197 G4double br, G4int iIso3, G4bool fAnti)
198{
199 G4VDecayChannel* mode;
200
201 G4String daughterN;
202 G4String daughterK;
203 G4double r = 0.;
204
205 // ------------ N K- ------------
206 // determine daughters
207 if (iIso3== +2) {
208 daughterN = "delta++";
209 r=0.75*br;
210 } else if (iIso3== 0) {
211 daughterN = "delta+";
212 r = br/2.;
213 } else if (iIso3== -2) {
214 daughterN = "delta0";
215 r = 0.25*br;
216 }
217 if (!fAnti) {
218 daughterK = "kaon-";
219 } else {
220 daughterK = "kaon+";
221 }
222 if (fAnti) daughterN = "anti_" + daughterN;
223 if (r>0.) {
224 // create decay channel [parent BR #daughters]
225 mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
226 daughterN,daughterK);
227 // add decay table
228 decayTable->Insert(mode);
229 }
230
231 // ------------ N K0 ------------
232 // determine daughters
233 if (iIso3== +2) {
234 daughterN = "delta+";
235 r=0.25*br;
236 } else if (iIso3== 0) {
237 daughterN = "delta0";
238 r = br/2.;
239 } else if (iIso3== -2) {
240 daughterN = "delta-";
241 r=0.75*br;
242 }
243 if (!fAnti) {
244 daughterK = "anti_kaon0";
245 } else {
246 daughterK = "kaon0";
247 }
248 if (fAnti) daughterN = "anti_" + daughterN;
249 if (r>0.) {
250
251
252 // create decay channel [parent BR #daughters]
253 mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
254 daughterN,daughterK);
255 // add decay table
256 decayTable->Insert(mode);
257 }
258
259 return decayTable;
260}
261
262
263G4DecayTable* G4ExcitedSigmaConstructor::AddNKStarMode(
264 G4DecayTable* decayTable, const G4String& nameParent,
265 G4double br, G4int iIso3, G4bool fAnti)
266{
267 G4VDecayChannel* mode;
268
269 G4String daughterN;
270 G4String daughterK;
271 G4double r = 0.;
272
273 // ------------ N K- ------------
274 // determine daughters
275 if (iIso3== +2) {
276
277 r=0.;
278 } else if (iIso3== 0) {
279 daughterN = "proton";
280 r = br/2.;
281 } else if (iIso3== -2) {
282 daughterN = "neutron";
283 r = br;
284 }
285 if (!fAnti) {
286 daughterK = "k_star-";
287 } else {
288 daughterK = "k_star+";
289 }
290 if (fAnti) daughterN = "anti_" + daughterN;
291 if (r>0.) {
292 // create decay channel [parent BR #daughters]
293 mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
294 daughterN,daughterK);
295 // add decay table
296 decayTable->Insert(mode);
297 }
298
299 // ------------ N K0 ------------
300
301 // determine daughters
302 if (iIso3== +2) {
303 daughterN = "proton";
304 r=br;
305 } else if (iIso3== 0) {
306 daughterN = "neutron";
307 r = br/2.;
308 } else if (iIso3== -2) {
309 r = 0.;
310 }
311 if (!fAnti) {
312 daughterK = "anti_k_star0";
313 } else {
314 daughterK = "k_star0";
315 }
316 if (fAnti) daughterN = "anti_" + daughterN;
317 // create decay channel [parent BR #daughters]
318 if (r>0.) {
319 mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
320 daughterN,daughterK);
321 // add decay table
322 decayTable->Insert(mode);
323 }
324
325 return decayTable;
326}
327
328G4DecayTable* G4ExcitedSigmaConstructor::AddSigmaPiMode(
329 G4DecayTable* decayTable, const G4String& nameParent,
330 G4double br, G4int iIso3, G4bool fAnti)
331{
332 G4VDecayChannel* mode;
333
334 G4String daughterSigma;
335 G4String daughterPi;
336 G4double r = 0.;
337
338 // ------------ Sigma+ pi - ------------
339 // determine daughters
340 if (iIso3== +2) {
341 r = 0.;
342 } else if (iIso3== 0) {
343 daughterSigma = "sigma+";
344 r = br/2.;
345 } else if (iIso3== -2) {
346 daughterSigma = "sigma0";
347 r = br/2.;
348 }
349 if (!fAnti) {
350 daughterPi = "pi-";
351 } else {
352 daughterPi = "pi+";
353 }
354 if (fAnti) daughterSigma = "anti_" + daughterSigma;
355 if (r>0.) {
356 // create decay channel [parent BR #daughters]
357 mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
358 daughterSigma,daughterPi);
359 // add decay table
360 decayTable->Insert(mode);
361 }
362 // ------------ Sigma0 Pi0 ------------
363 // determine daughters
364 if (iIso3== +2) {
365 daughterSigma = "sigma+";
366 r = br/2.;
367 } else if (iIso3== 0) {
368 r = 0.;
369 } else if (iIso3== -2) {
370 daughterSigma = "sigma-";
371
372
373 r = br/2.;
374 }
375 daughterPi = "pi0";
376 if (fAnti) daughterSigma = "anti_" + daughterSigma;
377 if (r>0.) {
378 // create decay channel [parent BR #daughters]
379 mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
380 daughterSigma,daughterPi);
381 // add decay table
382 decayTable->Insert(mode);
383 }
384
385 // ------------ Sigma- pi + ------------
386 // determine daughters
387 if (iIso3== +2) {
388 daughterSigma = "sigma0";
389 r = br/2.;
390 } else if (iIso3== 0) {
391 daughterSigma = "sigma-";
392 r = br/2.;
393 } else if (iIso3== -2) {
394 r = 0.;
395 }
396 if (!fAnti) {
397 daughterPi = "pi+";
398 } else {
399 daughterPi = "pi-";
400 }
401 if (fAnti) daughterSigma = "anti_" + daughterSigma;
402 if (r>0.) {
403 // create decay channel [parent BR #daughters]
404 mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
405 daughterSigma,daughterPi);
406 // add decay table
407 decayTable->Insert(mode);
408 }
409
410 return decayTable;
411}
412
413
414G4DecayTable* G4ExcitedSigmaConstructor::AddSigmaStarPiMode(
415 G4DecayTable* decayTable, const G4String& nameParent,
416 G4double br, G4int iIso3, G4bool fAnti)
417{
418 G4VDecayChannel* mode;
419
420
421 G4String daughterSigma;
422 G4String daughterPi;
423 G4double r = 0.;
424
425 // ------------ Sigma+ pi - ------------
426 // determine daughters
427 if (iIso3== +2) {
428 r = 0.;
429 } else if (iIso3== 0) {
430 daughterSigma = "sigma(1385)+";
431 r = br/2.;
432 } else if (iIso3== -2) {
433 daughterSigma = "sigma(1385)0";
434 r = br/2.;
435 }
436 if (!fAnti) {
437 daughterPi = "pi-";
438 } else {
439 daughterPi = "pi+";
440 }
441 if (fAnti) daughterSigma = "anti_" + daughterSigma;
442 if (r>0.) {
443 // create decay channel [parent BR #daughters]
444 mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
445 daughterSigma,daughterPi);
446 // add decay table
447 decayTable->Insert(mode);
448 }
449 // ------------ Sigma0 Pi0 ------------
450 // determine daughters
451 if (iIso3== +2) {
452 daughterSigma = "sigma(1385)+";
453 r = br/2.;
454 } else if (iIso3== 0) {
455 r = 0.;
456 } else if (iIso3== -2) {
457 daughterSigma = "sigma(1385)-";
458 r = br/2.;
459 }
460 daughterPi = "pi0";
461 if (fAnti) daughterSigma = "anti_" + daughterSigma;
462 if (r>0.) {
463 // create decay channel [parent BR #daughters]
464 mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
465 daughterSigma,daughterPi);
466
467 // add decay table
468 decayTable->Insert(mode);
469 }
470
471 // ------------ Sigma- pi + ------------
472 // determine daughters
473 if (iIso3== +2) {
474 daughterSigma = "sigma(1385)0";
475 r = br/2.;
476 } else if (iIso3== 0) {
477 daughterSigma = "sigma(1385)-";
478 r = br/2.;
479 } else if (iIso3== -2) {
480 r = 0.;
481 }
482 if (!fAnti) {
483 daughterPi = "pi+";
484 } else {
485 daughterPi = "pi-";
486 }
487 if (fAnti) daughterSigma = "anti_" + daughterSigma;
488 if (r>0.) {
489 // create decay channel [parent BR #daughters]
490
491 mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
492 daughterSigma,daughterPi);
493 // add decay table
494 decayTable->Insert(mode);
495 }
496
497 return decayTable;
498}
499
500G4DecayTable* G4ExcitedSigmaConstructor::AddLambdaPiMode(
501 G4DecayTable* decayTable, const G4String& nameParent,
502 G4double br, G4int iIso3, G4bool fAnti)
503{
504 G4VDecayChannel* mode;
505
506 G4String daughterLambda = "lambda";
507 G4String daughterPi;
508
509 // determine daughters
510 if (iIso3== +2) {
511 if (!fAnti) {
512 daughterPi = "pi+";
513 } else {
514 daughterPi = "pi-";
515 }
516 } else if (iIso3== 0) {
517 daughterPi = "pi0";
518 } else if (iIso3== -2) {
519 if (!fAnti) {
520 daughterPi = "pi-";
521 } else {
522 daughterPi = "pi+";
523 }
524 }
525 if (fAnti) daughterLambda = "anti_" + daughterLambda;
526 // create decay channel [parent BR #daughters]
527 mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
528 daughterLambda, daughterPi);
529 // add decay table
530 decayTable->Insert(mode);
531
532 return decayTable;
533}
534
535G4DecayTable* G4ExcitedSigmaConstructor::AddLambdaStarPiMode(
536 G4DecayTable* decayTable, const G4String& nameParent,
537 G4double br, G4int iIso3, G4bool fAnti)
538{
539 G4VDecayChannel* mode;
540
541 G4String daughterLambda = "lambda(1405)";
542 G4String daughterPi;
543
544 // determine daughters
545 if (iIso3== +2) {
546 if (!fAnti) {
547 daughterPi = "pi+";
548 } else {
549 daughterPi = "pi-";
550 }
551 } else if (iIso3== 0) {
552 daughterPi = "pi0";
553 } else if (iIso3== -2) {
554 if (!fAnti) {
555 daughterPi = "pi-";
556 } else {
557 daughterPi = "pi+";
558 }
559 }
560
561 if (fAnti) daughterLambda = "anti_" + daughterLambda;
562 // create decay channel [parent BR #daughters]
563 mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
564 daughterLambda,daughterPi);
565 // add decay table
566 decayTable->Insert(mode);
567
568 return decayTable;
569}
570
571
572G4double G4ExcitedSigmaConstructor::GetMass(G4int iState, G4int iso3)
573{
574 G4double fmass = mass[iState];
575 if (iState == 0 ) {
576 if (iso3== +2) fmass -= 0.9*MeV; // sigma+
577 else if (iso3== -2) fmass += 3.5*MeV; // sigma-
578 }
579 return fmass;
580}
581
582G4double G4ExcitedSigmaConstructor::GetWidth(G4int iState, G4int iso3)
583{
584 G4double fw=width[iState];
585 if (iState == 0 ) {
586 if (iso3== +2) fw = 35.8*MeV; // sigma+
587 else if (iso3== -2) fw = 39.4*MeV; // sigma-
588 }
589 return fw;
590}
591
592const char* G4ExcitedSigmaConstructor::name[] = {
593 "sigma(1385)","sigma(1660)","sigma(1670)","sigma(1750)","sigma(1775)",
594 "sigma(1915)","sigma(1940)","sigma(2030)"
595};
596
597const G4double G4ExcitedSigmaConstructor::mass[] = {
598 1.3837*GeV, 1.660*GeV, 1.670*GeV, 1.750*GeV, 1.775*GeV,
599 1.915*GeV, 1.940*GeV, 2.030*GeV
600};
601
602const G4double G4ExcitedSigmaConstructor::width[] = {
603 36.0*MeV, 100.0*MeV, 60.0*MeV, 90.0*MeV, 120.0*MeV,
604 120.0*MeV, 220.0*MeV, 180.0*MeV
605};
606
607const G4int G4ExcitedSigmaConstructor::iSpin[] = {
608 3, 1, 3, 1, 5,
609 5, 3, 7
610};
611
612const G4int G4ExcitedSigmaConstructor::iParity[] = {
613 +1, +1, -1, -1, -1,
614 +1, -1, +1
615};
616
617
618const G4int G4ExcitedSigmaConstructor::encodingOffset[] = {
619 0, 10000, 10000, 20000, 0,
620 10000, 20000, 0
621};
622
623const G4double G4ExcitedSigmaConstructor::bRatio[ G4ExcitedSigmaConstructor::NStates ][ G4ExcitedSigmaConstructor::NumberOfDecayModes] =
624{
625 { 0.0, 0.0, 0.12, 0.0, 0.88, 0.0, 0.0, 0.0},
626 { 0.30, 0.0, 0.35, 0.0, 0.35, 0.0, 0.0, 0.0},
627 { 0.15, 0.0, 0.70, 0.0, 0.15, 0.0, 0.0, 0.0},
628 { 0.40, 0.0, 0.05, 0.0, 0.0, 0.55, 0.0, 0.0},
629 { 0.40, 0.0, 0.04, 0.10, 0.23, 0.0, 0.23, 0.0},
630 { 0.15, 0.0, 0.40, 0.05, 0.40, 0.0, 0.0, 0.0},
631 { 0.10, 0.15, 0.15, 0.15, 0.15, 0.0, 0.15, 0.15},
632 { 0.20, 0.04, 0.10, 0.10, 0.20, 0.0, 0.18, 0.18}
633
634};
635
636
637
638
639
640
641
642
643
644
645
G4double
double G4double
Definition: G4Types.hh:64
G4int
int G4int
Definition: G4Types.hh:66
G4bool
bool G4bool
Definition: G4Types.hh:67
G4DecayTable::Insert
void Insert(G4VDecayChannel *aChannel)
Definition: G4DecayTable.cc:60
G4ExcitedSigmaConstructor::CreateDecayTable
virtual G4DecayTable * CreateDecayTable(const G4String &name, G4int iIso3, G4int iState, G4bool fAnti=false)
Definition: G4ExcitedSigmaConstructor.cc:58
G4ExcitedSigmaConstructor::GetWidth
virtual G4double GetWidth(G4int state, G4int iso)
Definition: G4ExcitedSigmaConstructor.cc:582
G4ExcitedSigmaConstructor::GetMass
virtual G4double GetMass(G4int state, G4int iso)
Definition: G4ExcitedSigmaConstructor.cc:572