source: trunk/source/particles/shortlived/include/G4ExcitedDeltaConstructor.hh@ 850

Last change on this file since 850 was 850, checked in by garnier, 17 years ago

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