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