| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
>> 26 // This code plementation is the intellectual property of
>> 27 // the GEANT4 collaboration.
>> 28 //
>> 29 // By copying, distributing or modifying the Program (or any work
>> 30 // based on the Program) you indicate your acceptance of this statement,
>> 31 // and all its terms.
>> 32 //
>> 33 // $Id: G4ExcitedMesonConstructor.hh 67971 2013-03-13 10:13:24Z gcosmo $
>> 34 //
>> 35 //
26 // ------------------------------------------- 36 // --------------------------------------------------------------
27 // GEANT 4 class implementation file << 37 // GEANT 4 class implementation file
28 // 38 //
29 // History: first implementation, based o 39 // History: first implementation, based on object model of
30 // 10 oct 1998 H.Kurashige 40 // 10 oct 1998 H.Kurashige
31 // ------------------------------------------- 41 // ---------------------------------------------------------------
32 #ifndef G4ExcitedMesonConstructor_h 42 #ifndef G4ExcitedMesonConstructor_h
33 #define G4ExcitedMesonConstructor_h 1 43 #define G4ExcitedMesonConstructor_h 1
34 44
35 #include "globals.hh" 45 #include "globals.hh"
36 class G4DecayTable; << 46 #include "G4ios.hh"
>> 47 class G4DecayTable;
37 48
38 class G4ExcitedMesonConstructor 49 class G4ExcitedMesonConstructor
39 { << 50 {
40 // This class is a utility class for const << 51 //This class is a utility class for construction
41 // short lived particles << 52 //short lived particles
42 public: <<
43 enum <<
44 { <<
45 NMultiplets = 10 <<
46 }; <<
47 <<
48 enum <<
49 { <<
50 NMesonTypes = 5 <<
51 }; <<
52 <<
53 enum <<
54 { <<
55 NumberOfDecayModes = 19 <<
56 }; <<
57 53
58 public: 54 public:
59 G4ExcitedMesonConstructor(G4int nStates = << 55 G4ExcitedMesonConstructor(G4int nStates = 0, G4int isoSpin=0);
60 virtual ~G4ExcitedMesonConstructor() = def << 56 virtual ~G4ExcitedMesonConstructor();
61 << 57
62 virtual void Construct(G4int indexOfState << 58 public:
63 << 59 virtual void Construct(G4int indexOfState = -1);
>> 60
64 protected: 61 protected:
65 void ConstructMesons(G4int indexOfState, G 62 void ConstructMesons(G4int indexOfState, G4int indexOfType);
66 << 63
67 G4String GetName(G4int iIso3, G4int iState 64 G4String GetName(G4int iIso3, G4int iState, G4int idxType);
68 G4double GetCharge(G4int iIsoSpin3); 65 G4double GetCharge(G4int iIsoSpin3);
69 G4int GetEncoding(G4int iIsoSpin3, G4int i << 66 G4int GetEncoding(G4int iIsoSpin3, G4int idxState, G4int idxType);
70 G4int GetQuarkContents(G4int iQ, G4int iIs << 67 G4int GetQuarkContents(G4int iQ, G4int iIso3, G4int iType);
71 68
72 protected: << 69 public:
73 G4DecayTable* CreateDecayTable(const G4Str << 70 enum { NMultiplets = 10 };
74 << 71 protected:
75 G4DecayTable* AddKPiMode(G4DecayTable* tab << 72 enum {
76 G4int iType); << 73 N11P1 = 0, N13P0 = 1, N13P1 = 2, N13P2 = 3,
77 G4DecayTable* AddKStarPiMode(G4DecayTable* << 74 N11D2 = 4, N13D1 = 5, N13D3 = 6,
78 G4int iIso3, << 75 N21S0 = 7, N23S1 = 8, N23P2 = 9
79 G4DecayTable* AddKStar2PiMode(G4DecayTable <<
80 G4int iIso3, <<
81 G4DecayTable* AddKRhoMode(G4DecayTable* ta <<
82 G4int iType); <<
83 G4DecayTable* AddKTwoPiMode(G4DecayTable* <<
84 G4int iType); <<
85 G4DecayTable* AddKOmegaMode(G4DecayTable* <<
86 G4int iType); <<
87 G4DecayTable* AddKEtaMode(G4DecayTable* ta <<
88 G4int iType); <<
89 G4DecayTable* AddPiGammaMode(G4DecayTable* <<
90 G4int iIso3, <<
91 G4DecayTable* AddRhoGammaMode(G4DecayTable <<
92 G4int iIso3, <<
93 G4DecayTable* Add2PiMode(G4DecayTable* tab <<
94 G4int iIso); <<
95 G4DecayTable* AddPiRhoMode(G4DecayTable* t <<
96 G4int iIso); <<
97 G4DecayTable* AddPiEtaMode(G4DecayTable* t <<
98 G4int iIso); <<
99 G4DecayTable* AddPiF2Mode(G4DecayTable* ta <<
100 G4int iIso); <<
101 G4DecayTable* AddPiF0Mode(G4DecayTable* ta <<
102 G4int iIso); <<
103 G4DecayTable* AddPiA2Mode(G4DecayTable* ta <<
104 G4int iIso); <<
105 G4DecayTable* Add3PiMode(G4DecayTable* tab <<
106 G4int iIso); <<
107 G4DecayTable* Add4PiMode(G4DecayTable* tab <<
108 G4int iIso); <<
109 G4DecayTable* AddKKStarMode(G4DecayTable* <<
110 G4int iIso); <<
111 G4DecayTable* Add2PiEtaMode(G4DecayTable* <<
112 G4int iIso); <<
113 G4DecayTable* AddRhoEtaMode(G4DecayTable* <<
114 G4int iIso); <<
115 G4DecayTable* Add2PiRhoMode(G4DecayTable* <<
116 G4int iIso); <<
117 G4DecayTable* Add2PiOmegaMode(G4DecayTable <<
118 G4int iIso3, <<
119 G4DecayTable* AddPiOmegaMode(G4DecayTable* <<
120 G4int iIso3, <<
121 G4DecayTable* Add2EtaMode(G4DecayTable* ta <<
122 G4int iIso); <<
123 G4DecayTable* Add2KMode(G4DecayTable* tabl <<
124 G4int iIso); <<
125 G4DecayTable* Add2KPiMode(G4DecayTable* ta <<
126 G4int iIso); <<
127 <<
128 G4bool Exist(G4int idxState, G4int idxType <<
129 G4double GetCharge(G4int iIsoSpin3, G4int <<
130 <<
131 protected: <<
132 enum <<
133 { <<
134 N11P1 = 0, <<
135 N13P0 = 1, <<
136 N13P1 = 2, <<
137 N13P2 = 3, <<
138 N11D2 = 4, <<
139 N13D1 = 5, <<
140 N13D3 = 6, <<
141 N21S0 = 7, <<
142 N23S1 = 8, <<
143 N23P2 = 9 <<
144 }; <<
145 <<
146 enum <<
147 { <<
148 TPi = 0, <<
149 TEta = 1, <<
150 TEtaPrime = 2, <<
151 TK = 3, <<
152 TAntiK = 4 <<
153 }; 76 };
>> 77
>> 78 public:
>> 79 enum { NMesonTypes = 5 };
>> 80 protected:
>> 81 enum { TPi=0, TEta=1, TEtaPrime=2, TK=3, TAntiK=4 };
154 82
>> 83 protected:
155 const G4String type; 84 const G4String type;
156 const G4int leptonNumber{0}; << 85 const G4int leptonNumber;
157 const G4int baryonNumber{0}; << 86 const G4int baryonNumber;
158 87
159 static const char* name[NMultiplets][NMeso << 88 G4bool Exist(G4int idxState, G4int idxType);
160 static const G4double mass[NMultiplets][NM << 89 G4double GetCharge(G4int iIsoSpin3, G4int idxType);
161 static const G4double massKdiff[NMultiplet << 90 static const char* name[ NMultiplets ][ NMesonTypes ];
162 static const G4double width[NMultiplets][N << 91 static const G4double mass[ NMultiplets ][ NMesonTypes ];
163 static const G4double widthKdiff[NMultiple << 92 static const G4double massKdiff[ NMultiplets ];
164 static const G4int iIsoSpin[NMesonTypes]; << 93 static const G4double width[ NMultiplets ][ NMesonTypes ];
165 static const G4int iSpin[NMultiplets]; << 94 static const G4double widthKdiff[ NMultiplets ];
166 static const G4int iParity[NMultiplets]; << 95 static const G4int iIsoSpin[ NMesonTypes ];
167 static const G4int iGParity[NMultiplets][N << 96 static const G4int iSpin[ NMultiplets ];
168 static const G4int iChargeConjugation[NMul << 97 static const G4int iParity[ NMultiplets ];
169 static const G4int encodingOffset[NMultipl << 98 static const G4int iGParity[ NMultiplets ][ NMesonTypes ];
170 << 99 static const G4int iChargeConjugation[ NMultiplets ];
171 enum << 100 static const G4int encodingOffset[ NMultiplets ];
172 { << 101
173 MPiGamma = 0, << 102 public:
174 MRhoGamma = 1, << 103 enum { NumberOfDecayModes = 19 };
175 M2Pi = 2, << 104 protected:
176 MPiRho = 3, << 105 enum { MPiGamma = 0, MRhoGamma=1, M2Pi=2, MPiRho=3,
177 M3Pi = 4, << 106 M3Pi= 4, MPiEta=5, M4Pi=6, MKKStar=7,
178 MPiEta = 5, << 107 M2PiEta=8, MRhoEta=9, M2PiRho=10, M2PiOmega=11,
179 M4Pi = 6, << 108 M2Eta=12, M2K=13, M2KPi=14, MPiOmega=15,
180 MKKStar = 7, << 109 MPiF2=16, MPiF0=17, MPiA2=18 };
181 M2PiEta = 8, << 110 enum { MKPi = 0, MKStarPi=1, MKRho=2, MKOmega=3,
182 MRhoEta = 9, << 111 MKStar2Pi=4, MKTwoPi=5, MKEta=6};
183 M2PiRho = 10, << 112
184 M2PiOmega = 11, << 113
185 M2Eta = 12, << 114 static const G4double bRatio[ NMultiplets ][ NMesonTypes ][ NumberOfDecayModes];
186 M2K = 13, << 115
187 M2KPi = 14, << 116 G4DecayTable* CreateDecayTable(const G4String&,
188 MPiOmega = 15, << 117 G4int , G4int, G4int);
189 MPiF2 = 16, << 118
190 MPiF0 = 17, << 119 G4DecayTable* AddKPiMode( G4DecayTable* table, const G4String& name,
191 MPiA2 = 18 << 120 G4double br, G4int iIso3, G4int iType);
192 }; << 121 G4DecayTable* AddKStarPiMode( G4DecayTable* table, const G4String& name,
193 enum << 122 G4double br, G4int iIso3, G4int iType);
194 { << 123 G4DecayTable* AddKStar2PiMode( G4DecayTable* table, const G4String& name,
195 MKPi = 0, << 124 G4double br, G4int iIso3, G4int iType);
196 MKStarPi = 1, << 125 G4DecayTable* AddKRhoMode( G4DecayTable* table, const G4String& name,
197 MKRho = 2, << 126 G4double br, G4int iIso3, G4int iType);
198 MKOmega = 3, << 127 G4DecayTable* AddKTwoPiMode( G4DecayTable* table, const G4String& name,
199 MKStar2Pi = 4, << 128 G4double br, G4int iIso3, G4int iType);
200 MKTwoPi = 5, << 129 G4DecayTable* AddKOmegaMode( G4DecayTable* table, const G4String& name,
201 MKEta = 6 << 130 G4double br, G4int iIso3, G4int iType);
202 }; << 131 G4DecayTable* AddKEtaMode( G4DecayTable* table, const G4String& name,
>> 132 G4double br, G4int iIso3, G4int iType);
>> 133 G4DecayTable* AddPiGammaMode( G4DecayTable* table, const G4String& name,
>> 134 G4double br, G4int iIso3,G4int iIso);
>> 135 G4DecayTable* AddRhoGammaMode( G4DecayTable* table, const G4String& name,
>> 136 G4double br, G4int iIso3, G4int iIso);
>> 137 G4DecayTable* Add2PiMode( G4DecayTable* table, const G4String& name,
>> 138 G4double br, G4int iIso3, G4int iIso);
>> 139 G4DecayTable* AddPiRhoMode( G4DecayTable* table, const G4String& name,
>> 140 G4double br, G4int iIso3, G4int iIso);
>> 141 G4DecayTable* AddPiEtaMode( G4DecayTable* table, const G4String& name,
>> 142 G4double br, G4int iIso3, G4int iIso);
>> 143 G4DecayTable* AddPiF2Mode( G4DecayTable* table, const G4String& name,
>> 144 G4double br, G4int iIso3, G4int iIso);
>> 145 G4DecayTable* AddPiF0Mode( G4DecayTable* table, const G4String& name,
>> 146 G4double br, G4int iIso3, G4int iIso);
>> 147 G4DecayTable* AddPiA2Mode( G4DecayTable* table, const G4String& name,
>> 148 G4double br, G4int iIso3, G4int iIso);
>> 149 G4DecayTable* Add3PiMode( G4DecayTable* table, const G4String& name,
>> 150 G4double br, G4int iIso3, G4int iIso);
>> 151 G4DecayTable* Add4PiMode( G4DecayTable* table, const G4String& name,
>> 152 G4double br, G4int iIso3, G4int iIso);
>> 153 G4DecayTable* AddKKStarMode( G4DecayTable* table, const G4String& name,
>> 154 G4double br, G4int iIso3, G4int iIso);
>> 155 G4DecayTable* Add2PiEtaMode( G4DecayTable* table, const G4String& name,
>> 156 G4double br, G4int iIso3, G4int iIso);
>> 157 G4DecayTable* AddRhoEtaMode( G4DecayTable* table, const G4String& name,
>> 158 G4double br, G4int iIso3, G4int iIso);
>> 159 G4DecayTable* Add2PiRhoMode( G4DecayTable* table, const G4String& name,
>> 160 G4double br, G4int iIso3, G4int iIso);
>> 161 G4DecayTable* Add2PiOmegaMode( G4DecayTable* table, const G4String& name,
>> 162 G4double br, G4int iIso3, G4int iIso);
>> 163 G4DecayTable* AddPiOmegaMode( G4DecayTable* table, const G4String& name,
>> 164 G4double br, G4int iIso3, G4int iIso);
>> 165 G4DecayTable* Add2EtaMode( G4DecayTable* table, const G4String& name,
>> 166 G4double br, G4int iIso3, G4int iIso);
>> 167 G4DecayTable* Add2KMode( G4DecayTable* table, const G4String& name,
>> 168 G4double br, G4int iIso3, G4int iIso);
>> 169 G4DecayTable* Add2KPiMode( G4DecayTable* table, const G4String& name,
>> 170 G4double br, G4int iIso3, G4int iIso);
>> 171
>> 172
203 173
204 static const G4double bRatio[NMultiplets][ <<
205 }; 174 };
206 175
207 inline G4String G4ExcitedMesonConstructor::Get << 176
>> 177 inline
>> 178 G4String G4ExcitedMesonConstructor::GetName(G4int iIso3,
>> 179 G4int iState,
>> 180 G4int iType)
208 { 181 {
209 G4String particle = name[iState][iType]; 182 G4String particle = name[iState][iType];
210 if (iType == TPi) { 183 if (iType == TPi) {
211 if (iIso3 == +2) { << 184 if ( iIso3 == +2 ){
212 particle += "+"; 185 particle += "+";
213 } << 186 } else if ( iIso3 == -2 ){
214 else if (iIso3 == -2) { <<
215 particle += "-"; 187 particle += "-";
216 } << 188 } else {
217 else { <<
218 particle += "0"; 189 particle += "0";
219 } 190 }
220 } << 191 } else if (iType == TK) {
221 else if (iType == TK) { << 192 if ( iIso3 == +1 ){
222 if (iIso3 == +1) { <<
223 particle += "+"; 193 particle += "+";
224 } << 194 } else if ( iIso3 == -1 ){
225 else if (iIso3 == -1) { <<
226 particle += "0"; 195 particle += "0";
227 } 196 }
228 } << 197 } else if (iType == TAntiK) {
229 else if (iType == TAntiK) { << 198 if ( iIso3 == +1 ){
230 if (iIso3 == +1) { <<
231 particle += "0"; 199 particle += "0";
232 particle = "anti_" + particle; 200 particle = "anti_" + particle;
233 } << 201 } else if ( iIso3 == -1 ){
234 else if (iIso3 == -1) { <<
235 particle += "-"; 202 particle += "-";
236 } 203 }
237 } << 204 }
238 return particle; 205 return particle;
239 } 206 }
240 207
241 #endif 208 #endif
242 209