| Geant4 Cross Reference |
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 // GEANT 4 class implementation file
28 //
29 // History: first implementation, based on object model of
30 // 10 oct 1998 H.Kurashige
31 // ---------------------------------------------------------------
32 #ifndef G4ExcitedSigmaConstructor_h
33 #define G4ExcitedSigmaConstructor_h 1
34
35 #include "G4ExcitedBaryonConstructor.hh"
36 #include "globals.hh"
37
38 class G4ExcitedSigmaConstructor : public G4ExcitedBaryonConstructor
39 {
40 // This class is a utility class for construction
41 // short lived particles
42 public:
43 enum
44 {
45 NStates = 8
46 };
47
48 enum
49 {
50 NumberOfDecayModes = 8
51 };
52
53 public:
54 G4ExcitedSigmaConstructor();
55 ~G4ExcitedSigmaConstructor() override = default;
56
57 protected:
58 G4bool Exist(G4int) override { return true; }
59
60 G4int GetQuarkContents(G4int, G4int) override;
61 G4String GetName(G4int iIso3, G4int iState) override;
62 G4String GetMultipletName(G4int iState) override;
63 G4double GetMass(G4int state, G4int iso) override;
64 G4double GetWidth(G4int state, G4int iso) override;
65 G4int GetiSpin(G4int iState) override;
66 G4int GetiParity(G4int iState) override;
67 G4int GetEncodingOffset(G4int iState) override;
68
69 G4DecayTable* CreateDecayTable(const G4String& name, G4int iIso3, G4int iState,
70 G4bool fAnti = false) override;
71
72 private:
73 G4DecayTable* AddNKMode(G4DecayTable* table, const G4String& name, G4double br, G4int iIso3,
74 G4bool fAnti);
75 G4DecayTable* AddNKStarMode(G4DecayTable* table, const G4String& name, G4double br, G4int iIso3,
76 G4bool fAnti);
77 G4DecayTable* AddSigmaPiMode(G4DecayTable* table, const G4String& name, G4double br,
78 G4int iIso3, G4bool fAnti);
79 G4DecayTable* AddSigmaStarPiMode(G4DecayTable* table, const G4String& name, G4double br,
80 G4int iIso3, G4bool fAnti);
81 G4DecayTable* AddLambdaPiMode(G4DecayTable* table, const G4String& name, G4double br,
82 G4int iIso3, G4bool fAnti);
83 G4DecayTable* AddSigmaEtaMode(G4DecayTable* table, const G4String& name, G4double br,
84 G4int iIso3, G4bool fAnti);
85 G4DecayTable* AddLambdaStarPiMode(G4DecayTable* table, const G4String& name, G4double br,
86 G4int iIso3, G4bool fAnti);
87 G4DecayTable* AddDeltaKMode(G4DecayTable* table, const G4String& name, G4double br, G4int iIso3,
88 G4bool fAnti);
89
90 private:
91 enum
92 {
93 SigmaIsoSpin = 2
94 };
95
96 static const char* name[NStates];
97 static const G4double mass[NStates];
98 static const G4double width[NStates];
99 static const G4int iSpin[NStates];
100 static const G4int iParity[NStates];
101 static const G4int encodingOffset[NStates];
102
103 enum
104 {
105 NK = 0,
106 NKStar = 1,
107 SigmaPi = 2,
108 SigmaStarPi = 3,
109 LambdaPi = 4,
110 SigmaEta = 5,
111 LambdaStarPi = 6,
112 DeltaK = 7
113 };
114
115 static const G4double bRatio[NStates][NumberOfDecayModes];
116 };
117
118 inline G4int G4ExcitedSigmaConstructor::GetiSpin(G4int iState)
119 {
120 return iSpin[iState];
121 }
122
123 inline G4int G4ExcitedSigmaConstructor::GetiParity(G4int iState)
124 {
125 return iParity[iState];
126 }
127
128 inline G4int G4ExcitedSigmaConstructor::GetEncodingOffset(G4int iState)
129 {
130 return encodingOffset[iState];
131 }
132
133 inline G4int G4ExcitedSigmaConstructor::GetQuarkContents(G4int iQ, G4int iIso3)
134 {
135 G4int quark = 0;
136 if (iQ == 0) {
137 // s-quark
138 quark = 3;
139 }
140 else if (iQ == 1) {
141 if (iIso3 == -2) {
142 // d-quark
143 quark = 1;
144 }
145 else {
146 // u-quark
147 quark = 2;
148 }
149 }
150 else if (iQ == 2) {
151 if (iIso3 == +2) {
152 // u-quark
153 quark = 2;
154 }
155 else {
156 // d-quark
157 quark = 1;
158 }
159 }
160 return quark;
161 }
162
163 inline G4String G4ExcitedSigmaConstructor::GetMultipletName(G4int iState)
164 {
165 return name[iState];
166 }
167
168 inline G4String G4ExcitedSigmaConstructor::GetName(G4int iIso3, G4int iState)
169 {
170 G4String particle = name[iState];
171 if (iIso3 == +2) {
172 particle += "+";
173 }
174 else if (iIso3 == 0) {
175 particle += "0";
176 }
177 else if (iIso3 == -2) {
178 particle += "-";
179 }
180 return particle;
181 }
182 #endif
183