00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036 #ifndef G4ExcitedNucleonConstructor_h
00037 #define G4ExcitedNucleonConstructor_h 1
00038
00039 #include "globals.hh"
00040 #include "G4ios.hh"
00041 #include "G4ExcitedBaryonConstructor.hh"
00042
00043 class G4ExcitedNucleonConstructor: public G4ExcitedBaryonConstructor
00044 {
00045
00046
00047
00048 public:
00049 G4ExcitedNucleonConstructor();
00050 virtual ~G4ExcitedNucleonConstructor();
00051
00052 protected:
00053 virtual G4int GetEncoding(G4int iIsoSpin3, G4int idxState);
00054
00055 protected:
00056 virtual G4bool Exist( G4int ){return true;}
00057
00058 virtual G4int GetQuarkContents(G4int, G4int);
00059 virtual G4String GetName(G4int iIso3, G4int iState);
00060 virtual G4String GetMultipletName(G4int iState);
00061 virtual G4double GetMass( G4int state, G4int iso);
00062 virtual G4double GetWidth( G4int state, G4int iso);
00063 virtual G4int GetiSpin(G4int iState);
00064 virtual G4int GetiParity(G4int iState);
00065 virtual G4int GetEncodingOffset(G4int iState);
00066
00067 virtual G4DecayTable* CreateDecayTable(const G4String& name,
00068 G4int iIso3, G4int iState,
00069 G4bool fAnti = false);
00070 private:
00071 G4DecayTable* AddNGammaMode( G4DecayTable* table, const G4String& name,
00072 G4double br, G4int iIso3, G4bool fAnti);
00073 G4DecayTable* AddNPiMode( G4DecayTable* table, const G4String& name,
00074 G4double br, G4int iIso3, G4bool fAnti);
00075 G4DecayTable* AddNEtaMode( G4DecayTable* table, const G4String& name,
00076 G4double br, G4int iIso3, G4bool fAnti);
00077 G4DecayTable* AddNOmegaMode( G4DecayTable* table, const G4String& name,
00078 G4double br, G4int iIso3, G4bool fAnti);
00079 G4DecayTable* AddNRhoMode( G4DecayTable* table, const G4String& name,
00080 G4double br, G4int iIso3, G4bool fAnti);
00081 G4DecayTable* AddN2PiMode( G4DecayTable* table, const G4String& name,
00082 G4double br, G4int iIso3, G4bool fAnti);
00083 G4DecayTable* AddDeltaPiMode( G4DecayTable* table, const G4String& name,
00084 G4double br, G4int iIso3, G4bool fAnti);
00085 G4DecayTable* AddNStarPiMode( G4DecayTable* table, const G4String& name,
00086 G4double br, G4int iIso3, G4bool fAnti);
00087 G4DecayTable* AddLambdaKMode( G4DecayTable* table, const G4String& name,
00088 G4double br, G4int iIso3, G4bool fAnti);
00089
00090 public:
00091 enum { NStates = 15 };
00092 private:
00093 enum { NucleonIsoSpin = 1 };
00094
00095 private:
00096 static const char* name[ NStates ];
00097 static const G4double mass[ NStates ];
00098 static const G4double width[ NStates ];
00099 static const G4int iSpin[ NStates ];
00100 static const G4int iParity[ NStates ];
00101 static const G4int encodingOffset[ NStates ];
00102
00103 public:
00104 enum { NumberOfDecayModes = 9 };
00105 private:
00106 enum { NGamma=0, NPi=1, NEta=2, NOmega=3, NRho=4,
00107 N2Pi=5, DeltaPi=6, NStarPi=7, LambdaK = 8 };
00108 private:
00109 static const G4double bRatio[ NStates ][ NumberOfDecayModes];
00110 };
00111
00112 inline
00113 G4double G4ExcitedNucleonConstructor::GetMass(G4int iState, G4int)
00114 {
00115 return mass[iState];
00116 }
00117
00118 inline
00119 G4double G4ExcitedNucleonConstructor::GetWidth(G4int iState, G4int)
00120 {
00121 return width[iState];
00122 }
00123
00124 inline
00125 G4int G4ExcitedNucleonConstructor::GetiSpin(G4int iState)
00126 {
00127 return iSpin[iState];
00128 }
00129
00130 inline
00131 G4int G4ExcitedNucleonConstructor::GetiParity(G4int iState)
00132 {
00133 return iParity[iState];
00134 }
00135
00136 inline
00137 G4int G4ExcitedNucleonConstructor::GetEncodingOffset(G4int iState)
00138 {
00139 return encodingOffset[iState];
00140 }
00141
00142 inline
00143 G4int G4ExcitedNucleonConstructor::GetQuarkContents(G4int iQ, G4int iIso3)
00144 {
00145
00146
00147
00148 G4int quark=0;
00149 if ( iQ == 0 ){
00150
00151 quark = 2;
00152 } else if ( iQ == 2 ){
00153
00154 quark = 1;
00155 } else {
00156 if ( iIso3 == -1 ){
00157
00158 quark = 1;
00159 } else {
00160
00161 quark = 2;
00162 }
00163 }
00164 return quark;
00165 }
00166
00167 inline
00168 G4String G4ExcitedNucleonConstructor::GetMultipletName(G4int iState)
00169 {
00170 return name[iState];
00171 }
00172
00173 inline
00174 G4String G4ExcitedNucleonConstructor::GetName(G4int iIso3, G4int iState)
00175 {
00176 G4String particle = name[iState];
00177 if ( iIso3 == -1 ){
00178 particle += "0";
00179 } else {
00180 particle += "+";
00181 }
00182 return particle;
00183 }
00184 #endif
00185
00186
00187
00188
00189
00190
00191
00192
00193
00194