source: trunk/source/particles/shortlived/include/G4ExcitedMesonConstructor.hh@ 1199

Last change on this file since 1199 was 1196, checked in by garnier, 16 years ago

update CVS release candidate geant4.9.3.01
File size: 8.8 KB
Line 
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// 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,v 1.8 2006/06/29 19:26:36 gunter Exp $
34// GEANT4 tag $Name: geant4-09-03-cand-01 $
35//
36//
37// --------------------------------------------------------------
38// GEANT 4 class implementation file
39//
40// History: first implementation, based on object model of
41// 10 oct 1998 H.Kurashige
42// ---------------------------------------------------------------
43#ifndef G4ExcitedMesonConstructor_h
44#define G4ExcitedMesonConstructor_h 1
45
46#include "globals.hh"
47#include "G4ios.hh"
48class G4DecayTable;
49
50class G4ExcitedMesonConstructor
51{
52 //This class is a utility class for construction
53 //short lived particles
54
55 public:
56 G4ExcitedMesonConstructor(G4int nStates = 0, G4int isoSpin=0);
57 virtual ~G4ExcitedMesonConstructor();
58
59 public:
60 virtual void Construct(G4int indexOfState = -1);
61
62 protected:
63 void ConstructMesons(G4int indexOfState, G4int indexOfType);
64
65 G4String GetName(G4int iIso3, G4int iState, G4int idxType);
66 G4double GetCharge(G4int iIsoSpin3);
67 G4int GetEncoding(G4int iIsoSpin3, G4int idxState, G4int idxType);
68 G4int GetQuarkContents(G4int iQ, G4int iIso3, G4int iType);
69
70 public:
71 enum { NMultiplets = 10 };
72 protected:
73 enum {
74 N11P1 = 0, N13P0 = 1, N13P1 = 2, N13P2 = 3,
75 N11D2 = 4, N13D1 = 5, N13D3 = 6,
76 N21S0 = 7, N23S1 = 8, N23P2 = 9
77 };
78
79 public:
80 enum { NMesonTypes = 5 };
81 protected:
82 enum { TPi=0, TEta=1, TEtaPrime=2, TK=3, TAntiK=4 };
83
84 protected:
85 const G4String type;
86 const G4int leptonNumber;
87 const G4int baryonNumber;
88
89 G4bool Exist(G4int idxState, G4int idxType);
90 G4double GetCharge(G4int iIsoSpin3, G4int idxType);
91 static const char* name[ NMultiplets ][ NMesonTypes ];
92 static const G4double mass[ NMultiplets ][ NMesonTypes ];
93 static const G4double massKdiff[ NMultiplets ];
94 static const G4double width[ NMultiplets ][ NMesonTypes ];
95 static const G4double widthKdiff[ NMultiplets ];
96 static const G4int iIsoSpin[ NMesonTypes ];
97 static const G4int iSpin[ NMultiplets ];
98 static const G4int iParity[ NMultiplets ];
99 static const G4int iGParity[ NMultiplets ][ NMesonTypes ];
100 static const G4int iChargeConjugation[ NMultiplets ];
101 static const G4int encodingOffset[ NMultiplets ];
102
103 public:
104 enum { NumberOfDecayModes = 19 };
105 protected:
106 enum { MPiGamma = 0, MRhoGamma=1, M2Pi=2, MPiRho=3,
107 M3Pi= 4, MPiEta=5, M4Pi=6, MKKStar=7,
108 M2PiEta=8, MRhoEta=9, M2PiRho=10, M2PiOmega=11,
109 M2Eta=12, M2K=13, M2KPi=14, MPiOmega=15,
110 MPiF2=16, MPiF0=17, MPiA2=18 };
111 enum { MKPi = 0, MKStarPi=1, MKRho=2, MKOmega=3,
112 MKStar2Pi=4, MKTwoPi=5, MKEta=6};
113
114
115 static const G4double bRatio[ NMultiplets ][ NMesonTypes ][ NumberOfDecayModes];
116
117 G4DecayTable* CreateDecayTable(const G4String&,
118 G4int , G4int, G4int);
119
120 G4DecayTable* AddKPiMode( G4DecayTable* table, const G4String& name,
121 G4double br, G4int iIso3, G4int iType);
122 G4DecayTable* AddKStarPiMode( G4DecayTable* table, const G4String& name,
123 G4double br, G4int iIso3, G4int iType);
124 G4DecayTable* AddKStar2PiMode( G4DecayTable* table, const G4String& name,
125 G4double br, G4int iIso3, G4int iType);
126 G4DecayTable* AddKRhoMode( G4DecayTable* table, const G4String& name,
127 G4double br, G4int iIso3, G4int iType);
128 G4DecayTable* AddKTwoPiMode( G4DecayTable* table, const G4String& name,
129 G4double br, G4int iIso3, G4int iType);
130 G4DecayTable* AddKOmegaMode( G4DecayTable* table, const G4String& name,
131 G4double br, G4int iIso3, G4int iType);
132 G4DecayTable* AddKEtaMode( G4DecayTable* table, const G4String& name,
133 G4double br, G4int iIso3, G4int iType);
134 G4DecayTable* AddPiGammaMode( G4DecayTable* table, const G4String& name,
135 G4double br, G4int iIso3,G4int iIso);
136 G4DecayTable* AddRhoGammaMode( G4DecayTable* table, const G4String& name,
137 G4double br, G4int iIso3, G4int iIso);
138 G4DecayTable* Add2PiMode( G4DecayTable* table, const G4String& name,
139 G4double br, G4int iIso3, G4int iIso);
140 G4DecayTable* AddPiRhoMode( G4DecayTable* table, const G4String& name,
141 G4double br, G4int iIso3, G4int iIso);
142 G4DecayTable* AddPiEtaMode( G4DecayTable* table, const G4String& name,
143 G4double br, G4int iIso3, G4int iIso);
144 G4DecayTable* AddPiF2Mode( G4DecayTable* table, const G4String& name,
145 G4double br, G4int iIso3, G4int iIso);
146 G4DecayTable* AddPiF0Mode( G4DecayTable* table, const G4String& name,
147 G4double br, G4int iIso3, G4int iIso);
148 G4DecayTable* AddPiA2Mode( G4DecayTable* table, const G4String& name,
149 G4double br, G4int iIso3, G4int iIso);
150 G4DecayTable* Add3PiMode( G4DecayTable* table, const G4String& name,
151 G4double br, G4int iIso3, G4int iIso);
152 G4DecayTable* Add4PiMode( G4DecayTable* table, const G4String& name,
153 G4double br, G4int iIso3, G4int iIso);
154 G4DecayTable* AddKKStarMode( G4DecayTable* table, const G4String& name,
155 G4double br, G4int iIso3, G4int iIso);
156 G4DecayTable* Add2PiEtaMode( G4DecayTable* table, const G4String& name,
157 G4double br, G4int iIso3, G4int iIso);
158 G4DecayTable* AddRhoEtaMode( G4DecayTable* table, const G4String& name,
159 G4double br, G4int iIso3, G4int iIso);
160 G4DecayTable* Add2PiRhoMode( G4DecayTable* table, const G4String& name,
161 G4double br, G4int iIso3, G4int iIso);
162 G4DecayTable* Add2PiOmegaMode( G4DecayTable* table, const G4String& name,
163 G4double br, G4int iIso3, G4int iIso);
164 G4DecayTable* AddPiOmegaMode( G4DecayTable* table, const G4String& name,
165 G4double br, G4int iIso3, G4int iIso);
166 G4DecayTable* Add2EtaMode( G4DecayTable* table, const G4String& name,
167 G4double br, G4int iIso3, G4int iIso);
168 G4DecayTable* Add2KMode( G4DecayTable* table, const G4String& name,
169 G4double br, G4int iIso3, G4int iIso);
170 G4DecayTable* Add2KPiMode( G4DecayTable* table, const G4String& name,
171 G4double br, G4int iIso3, G4int iIso);
172
173
174
175};
176
177
178inline
179 G4String G4ExcitedMesonConstructor::GetName(G4int iIso3,
180 G4int iState,
181 G4int iType)
182{
183 G4String particle = name[iState][iType];
184 if (iType == TPi) {
185 if ( iIso3 == +2 ){
186 particle += "+";
187 } else if ( iIso3 == -2 ){
188 particle += "-";
189 } else {
190 particle += "0";
191 }
192 } else if (iType == TK) {
193 if ( iIso3 == +1 ){
194 particle += "+";
195 } else if ( iIso3 == -1 ){
196 particle += "0";
197 }
198 } else if (iType == TAntiK) {
199 if ( iIso3 == +1 ){
200 particle += "0";
201 particle = "anti_" + particle;
202 } else if ( iIso3 == -1 ){
203 particle += "-";
204 }
205 }
206 return particle;
207}
208
209#endif
Note: See TracBrowser for help on using the repository browser.