| 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 G4ExcitedXiConstructor_h
33 #define G4ExcitedXiConstructor_h 1
34
35 #include "G4ExcitedBaryonConstructor.hh"
36 #include "globals.hh"
37
38 class G4ExcitedXiConstructor : public G4ExcitedBaryonConstructor
39 {
40 // This class is a utility class for construction
41 // short lived particles
42 public:
43 enum
44 {
45 NStates = 5
46 };
47
48 enum
49 {
50 NumberOfDecayModes = 4
51 };
52
53 public:
54 G4ExcitedXiConstructor();
55 ~G4ExcitedXiConstructor() override = default;
56
57 protected:
58 G4bool Exist(G4int) override { return true; }
59
60 G4int GetQuarkContents(G4int, G4int) override;
61 G4String GetName(G4int iIso3, G4int iState) override;
62 G4String GetMultipletName(G4int iState) override;
63 G4double GetMass(G4int state, G4int iso) override;
64 G4double GetWidth(G4int state, G4int iso) override;
65 G4int GetiSpin(G4int iState) override;
66 G4int GetiParity(G4int iState) override;
67 G4int GetEncodingOffset(G4int iState) override;
68
69 G4DecayTable* CreateDecayTable(const G4String& name, G4int iIso3, G4int iState,
70 G4bool fAnti = false) override;
71
72 private:
73 G4DecayTable* AddXiPiMode(G4DecayTable* table, const G4String& name, G4double br, G4int iIso3,
74 G4bool fAnti);
75 G4DecayTable* AddXiGammaMode(G4DecayTable* table, const G4String& name, G4double br,
76 G4int iIso3, G4bool fAnti);
77 G4DecayTable* AddLambdaKMode(G4DecayTable* table, const G4String& name, G4double br,
78 G4int iIso3, G4bool fAnti);
79 G4DecayTable* AddSigmaKMode(G4DecayTable* table, const G4String& name, G4double br, G4int iIso3,
80 G4bool fAnti);
81
82 private:
83 enum
84 {
85 XiIsoSpin = 1
86 };
87
88 static const char* name[NStates];
89 static const G4double mass[NStates];
90 static const G4double width[NStates];
91 static const G4int iSpin[NStates];
92 static const G4int iParity[NStates];
93 static const G4int encodingOffset[NStates];
94
95 enum
96 {
97 XiPi = 0,
98 XiGamma = 1,
99 LambdaK = 2,
100 SigmaK = 3
101 };
102
103 static const G4double bRatio[NStates][NumberOfDecayModes];
104 };
105
106 inline G4int G4ExcitedXiConstructor::GetiSpin(G4int iState)
107 {
108 return iSpin[iState];
109 }
110
111 inline G4int G4ExcitedXiConstructor::GetiParity(G4int iState)
112 {
113 return iParity[iState];
114 }
115
116 inline G4int G4ExcitedXiConstructor::GetEncodingOffset(G4int iState)
117 {
118 return encodingOffset[iState];
119 }
120
121 inline G4int G4ExcitedXiConstructor::GetQuarkContents(G4int iQ, G4int iIso3)
122 {
123 G4int quark = 0;
124 if (iQ == 0) {
125 // s-quark
126 quark = 3;
127 }
128 else if (iQ == 1) {
129 // s-quark
130 quark = 3;
131 }
132 else if (iQ == 2) {
133 if (iIso3 == +1) {
134 // u-quark
135 quark = 2;
136 }
137 else {
138 // d-quark
139 quark = 1;
140 }
141 }
142 return quark;
143 }
144
145 inline G4String G4ExcitedXiConstructor::GetMultipletName(G4int iState)
146 {
147 return name[iState];
148 }
149
150 inline G4String G4ExcitedXiConstructor::GetName(G4int iIso3, G4int iState)
151 {
152 G4String particle = name[iState];
153 if (iIso3 == +1) {
154 particle += "0";
155 }
156 else if (iIso3 == -1) {
157 particle += "-";
158 }
159 return particle;
160 }
161 #endif
162