| 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 G4ExcitedDeltaConstructor_h
33 #define G4ExcitedDeltaConstructor_h 1
34
35 #include "G4ExcitedBaryonConstructor.hh"
36 #include "globals.hh"
37
38 class G4ExcitedDeltaConstructor : public G4ExcitedBaryonConstructor
39 {
40 // This class is a utility class for construction
41 // short lived particles
42 public:
43 enum
44 {
45 NStates = 9
46 };
47
48 enum
49 {
50 NumberOfDecayModes = 5
51 };
52
53 public:
54 G4ExcitedDeltaConstructor();
55 ~G4ExcitedDeltaConstructor() override = default;
56
57 protected:
58 G4int GetEncoding(G4int iIsoSpin3, G4int idxState) override;
59
60 G4bool Exist(G4int) override { return true; }
61
62 G4int GetQuarkContents(G4int, G4int) override;
63 G4String GetName(G4int iIso3, G4int iState) override;
64 G4String GetMultipletName(G4int iState) override;
65
66 G4double GetMass(G4int state, G4int iso) override;
67 G4double GetWidth(G4int state, G4int iso) override;
68 G4int GetiSpin(G4int iState) override;
69 G4int GetiParity(G4int iState) override;
70 G4int GetEncodingOffset(G4int iState) override;
71
72 G4DecayTable* CreateDecayTable(const G4String& name, G4int iIso3, G4int iState,
73 G4bool fAnti = false) override;
74
75 private:
76 G4DecayTable* AddNGammaMode(G4DecayTable* table, const G4String& name, G4double br, G4int iIso3,
77 G4bool fAnti);
78 G4DecayTable* AddNPiMode(G4DecayTable* table, const G4String& name, G4double br, G4int iIso3,
79 G4bool fAnti);
80 G4DecayTable* AddNRhoMode(G4DecayTable* table, const G4String& name, G4double br, G4int iIso3,
81 G4bool fAnti);
82 G4DecayTable* AddDeltaPiMode(G4DecayTable* table, const G4String& name, G4double br,
83 G4int iIso3, G4bool fAnti);
84 G4DecayTable* AddNStarPiMode(G4DecayTable* table, const G4String& name, G4double br,
85 G4int iIso3, G4bool fAnti);
86
87 private:
88 enum
89 {
90 DeltaIsoSpin = 3
91 };
92
93 static const char* name[NStates];
94 static const G4double mass[NStates];
95 static const G4double width[NStates];
96 static const G4int iSpin[NStates];
97 static const G4int iParity[NStates];
98 static const G4int encodingOffset[NStates];
99
100 enum
101 {
102 NGamma = 0,
103 NPi = 1,
104 NRho = 2,
105 DeltaPi = 3,
106 NStarPi = 4
107 };
108
109 static const G4double bRatio[NStates][NumberOfDecayModes];
110 };
111
112 inline G4double G4ExcitedDeltaConstructor::GetMass(G4int iState, G4int)
113 {
114 return mass[iState];
115 }
116
117 inline G4double G4ExcitedDeltaConstructor::GetWidth(G4int iState, G4int)
118 {
119 return width[iState];
120 }
121
122 inline G4int G4ExcitedDeltaConstructor::GetiSpin(G4int iState)
123 {
124 return iSpin[iState];
125 }
126
127 inline G4int G4ExcitedDeltaConstructor::GetiParity(G4int iState)
128 {
129 return iParity[iState];
130 }
131
132 inline G4int G4ExcitedDeltaConstructor::GetEncodingOffset(G4int iState)
133 {
134 return encodingOffset[iState];
135 }
136
137 inline G4int G4ExcitedDeltaConstructor::GetQuarkContents(G4int iQ, G4int iIso3)
138 {
139 // Quark contents
140 // iIso3 = +3 : uuu
141 // iIso3 = +1 : uud
142 // iIso3 = -1 : udd
143 // iIso3 = -3 : ddd
144 G4int quark = 0;
145 if (iQ == 0) {
146 if (iIso3 == -3) {
147 // d-quark
148 quark = 1;
149 }
150 else {
151 // u-quark
152 quark = 2;
153 }
154 }
155 else if (iQ == 2) {
156 if (iIso3 == +3) {
157 // u-quark
158 quark = 2;
159 }
160 else {
161 // d-quark
162 quark = 1;
163 }
164 }
165 else {
166 if ((iIso3 == -1) || (iIso3 == -3)) {
167 // d-quark
168 quark = 1;
169 }
170 else {
171 // u-quark
172 quark = 2;
173 }
174 }
175 return quark;
176 }
177
178 inline G4String G4ExcitedDeltaConstructor::GetMultipletName(G4int iState)
179 {
180 return name[iState];
181 }
182
183 inline G4String G4ExcitedDeltaConstructor::GetName(G4int iIso3, G4int iState)
184 {
185 G4String particle = name[iState];
186 if (iIso3 == -3) {
187 particle += "-";
188 }
189 else if (iIso3 == -1) {
190 particle += "0";
191 }
192 else if (iIso3 == +1) {
193 particle += "+";
194 }
195 else {
196 particle += "++";
197 }
198 return particle;
199 }
200 #endif
201