Geant4  10.00.p02
G4ExcitedBaryonConstructor.cc
Go to the documentation of this file.
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: G4ExcitedBaryonConstructor.cc 77701 2013-11-27 16:39:30Z gcosmo $
28 //
29 //
30 // --------------------------------------------------------------
31 // GEANT 4 class implementation file
32 //
33 // History: first implementation, based on object model of
34 // 10 oct 1998 H.Kurashige
35 // ---------------------------------------------------------------
36 
37 
39 
40 #include "G4SystemOfUnits.hh"
41 #include "G4ParticleDefinition.hh"
42 #include "G4ParticleTable.hh"
44 #include "G4VDecayChannel.hh"
45 #include "G4DecayTable.hh"
46 
47 
49  G4int isoSpin)
50  : NumberOfStates(nStates), iIsoSpin(isoSpin), type("baryon"),
51  iConjugation(0), iGParity(0), leptonNumber(0), baryonNumber(1)
52 {
53 }
54 
56 {
57 }
58 
60 {
61  if (idx < 0 ) {
62  for (G4int state=0; state< NumberOfStates; state +=1) {
63  ConstructParticle(state);
64  ConstructAntiParticle(state);
65  }
66  } else if (idx < NumberOfStates) {
67  ConstructParticle(idx);
69  } else {
70 #ifdef G4VERBOSE
71  if (G4ParticleTable::GetParticleTable()->GetVerboseLevel()>1) {
72  G4cerr << "G4ExcitedBaryonConstructor::Construct()";
73  G4cerr << " illegal index os state = " << idx << G4endl;
74  }
75 #endif
76  }
77 }
78 
79 
80 #include "G4ExcitedBaryons.hh"
81 
83 {
84  if (!Exist(idx) ) return;
85 
86  // Construct Resonace particles as dynamic object
87  // Arguments for constructor are as follows
88  // name mass width charge
89  // 2*spin parity C-conjugation
90  // 2*Isospin 2*Isospin3 G-parity
91  // type lepton number baryon number PDG encoding
92  // stable lifetime decay table
93 
94 
95  G4String name;
96  G4ParticleDefinition* particle;
97 
98  for (G4int iIso3 = -1*iIsoSpin; iIso3 <= iIsoSpin; iIso3 +=2) {
99  name= GetName(iIso3, idx);
100 
101  particle = new G4ExcitedBaryons(
102  name, GetMass(idx,iIso3), GetWidth(idx,iIso3), GetCharge(iIso3),
103  GetiSpin(idx), GetiParity(idx), iConjugation,
104  iIsoSpin, iIso3, iGParity,
105  type, leptonNumber, baryonNumber, GetEncoding( iIso3,idx),
106  false, 0.0, NULL
107  );
108  ((G4ExcitedBaryons*)(particle))->SetMultipletName(GetMultipletName(idx));
109  particle->SetDecayTable(CreateDecayTable( name, iIso3, idx, false));
110  }
111 }
112 
114 {
115  if (!Exist(idx) ) return;
116 
117  // Construct Resonace particles as dynamic object
118  // Arguments for constructor are as follows
119  // name mass width charge
120  // 2*spin parity C-conjugation
121  // 2*Isospin 2*Isospin3 G-parity
122  // type lepton number baryon number PDG encoding
123  // stable lifetime decay table
124 
125 
126  G4String name;
127  G4ParticleDefinition* particle;
128 
129  for (G4int iIso3 = -1*iIsoSpin; iIso3 <= iIsoSpin; iIso3 +=2) {
130  name = GetName(iIso3, idx);
131  name = "anti_" + name;
132 
133  particle = new G4ExcitedBaryons(
134  name, GetMass(idx,iIso3), GetWidth(idx,iIso3), -1.0*GetCharge(iIso3),
135  GetiSpin(idx), GetiParity(idx), iConjugation,
136  iIsoSpin, -1*iIso3, iGParity,
137  type, leptonNumber,
138  -1*baryonNumber,
139  -1*GetEncoding( iIso3,idx),
140  false, 0.0, NULL
141  );
142 
143  ((G4ExcitedBaryons*)(particle))->SetMultipletName(GetMultipletName(idx));
144  particle->SetDecayTable(CreateDecayTable( name, iIso3, idx, true));
145  }
146 
147 }
148 
150 {
151  G4double charge = 0.0;
152  static const G4double quark_charge[7] =
153  {
154  0., -1./3., +2./3., -1./3., +2./3., -1./3., +2./3.
155  };
156 
157  for (G4int idx=0; idx<3; idx+=1){
158  charge += quark_charge[GetQuarkContents(idx, iIsoSpin3)]*eplus;
159  }
160  return charge;
161 }
162 
164 {
165  G4int encoding = GetEncodingOffset(idxState);
166  encoding += 1000*GetQuarkContents(0, iIsoSpin3);
167  encoding += 100*GetQuarkContents(1, iIsoSpin3);
168  encoding += 10*GetQuarkContents(2, iIsoSpin3);
169  if (GetiSpin(idxState) <9) {
170  encoding += GetiSpin(idxState) +1;
171  } else {
172  encoding += (GetiSpin(idxState) +1)*10000000;
173  }
174  return encoding;
175 }
void SetDecayTable(G4DecayTable *aDecayTable)
G4ExcitedBaryonConstructor(G4int nStates=0, G4int isoSpin=0)
virtual G4double GetCharge(G4int iIsoSpin3)
G4String name
Definition: TRTMaterials.hh:40
virtual G4DecayTable * CreateDecayTable(const G4String &, G4int, G4int, G4bool)=0
virtual G4bool Exist(G4int)=0
virtual G4double GetMass(G4int state, G4int iso)=0
int G4int
Definition: G4Types.hh:78
virtual G4String GetName(G4int, G4int)=0
virtual void ConstructParticle(G4int indexOfState)
virtual G4String GetMultipletName(G4int)=0
virtual G4int GetiParity(G4int)=0
virtual void Construct(G4int indexOfState=-1)
#define encoding
Definition: xmlparse.cc:588
virtual G4double GetWidth(G4int state, G4int iso)=0
static G4ParticleTable * GetParticleTable()
virtual G4int GetEncodingOffset(G4int)=0
virtual G4int GetiSpin(G4int)=0
#define G4endl
Definition: G4ios.hh:61
virtual void ConstructAntiParticle(G4int indexOfState)
virtual G4int GetEncoding(G4int iIsoSpin3, G4int idxState)
double G4double
Definition: G4Types.hh:76
static const double eplus
Definition: G4SIunits.hh:178
virtual G4int GetQuarkContents(G4int, G4int)=0
G4GLOB_DLL std::ostream G4cerr