Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4QParticle.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$
28 //
29 // ---------------- G4QParticle ----------------
30 // by Mikhail Kossov, Sept 1999.
31 // class for Particles in the CHIPS Model
32 // -------------------------------------------------------------------
33 // 1 2 3 4 5 6 7 8 9
34 //34567890123456789012345678901234567890123456789012345678901234567890123456789012345678901
35 // -----------------------------------------------------------------------
36 // Short description: The G4QParticle is a part of the CHIPS World. It is
37 // characterized by the quark content, spin, mass, width and a vector of
38 // the decay channels (G4QDecayCannelVector).
39 // -----------------------------------------------------------------------
40 //#define debug
41 //#define pdebug
42 
43 #include "G4QParticleVector.hh"
44 
45 G4QParticle::G4QParticle()
46 {
47 #ifdef debug
48  G4cout<<"G4QParticle::Constr: Default constructor is called"<<G4endl;
49 #endif
50 }
51 
52 G4QParticle::G4QParticle(G4int thePDG)
53 {
54  aQPDG = G4QPDGCode(thePDG);
55  aQuarkCont = aQPDG.GetQuarkContent();
56  aDecay = InitDecayVector(aQPDG.GetQCode());
57 }
58 
59 G4QParticle::G4QParticle(G4bool f, G4int theQCode)
60 {
61  aQPDG = G4QPDGCode(f,theQCode);
62 #ifdef debug
63  G4cout<<"G4QParticle::Constr: PDG="<<aQPDG.GetPDGCode()<<G4endl;
64 #endif
65  aQuarkCont = aQPDG.GetQuarkContent();
66  aDecay = InitDecayVector(theQCode);
67 }
68 
69 G4QParticle::G4QParticle(const G4QParticle& right)
70 {
71  aQPDG = right.aQPDG;
72  //aDecay (Vector)
73  G4int nD = right.aDecay.size();
74  if(nD) for(G4int id=0; id<nD; id++)
75  {
76  G4QDecayChan* curD = new G4QDecayChan(right.aDecay[id]);
77  aDecay.push_back(curD);
78  }
79 
80  aQuarkCont = right.aQuarkCont;
81 }
82 
83 G4QParticle::G4QParticle(G4QParticle* right)
84 {
85  aQPDG = right->aQPDG;
86  //aDecay (Vector)
87  G4int nD = right->aDecay.size();
88  if(nD) for(G4int id=0; id<nD; id++)
89  {
90  G4QDecayChan* curD = new G4QDecayChan(right->aDecay[id]);
91  aDecay.push_back(curD);
92  }
93 
94  aQuarkCont = right->aQuarkCont;
95 }
96 
97 G4QParticle::~G4QParticle()
98 {
99  G4int nDC=aDecay.size();
100  //G4cout<<"G4QParticle::Destructor: Before nDC="<<nDC<<G4endl; // TMP
101  if(nDC) std::for_each(aDecay.begin(), aDecay.end(), DeleteQDecayChan());
102  //G4cout<<"G4QParticle::Destructor: After"<<G4endl; // TMP
103  aDecay.clear();
104 }
105 
106 // Assignment operator
107 const G4QParticle& G4QParticle::operator=(const G4QParticle &right)
108 {
109  if(this != &right) // Beware of self assignment
110  {
111  aQPDG = right.aQPDG;
112  //aDecay (Vector)
113  G4int iD = aDecay.size();
114  if(iD) for(G4int jd=0; jd<iD; jd++) delete aDecay[jd];
115  aDecay.clear();
116  G4int nD = right.aDecay.size();
117  if(nD) for(G4int id=0; id<nD; id++)
118  {
119  G4QDecayChan* curD = new G4QDecayChan(right.aDecay[id]);
120  aDecay.push_back(curD);
121  }
122 
123  aQuarkCont = right.aQuarkCont;
124  }
125  return *this;
126 }
127 
128 // Standard output for QParticle
129 std::ostream& operator<<(std::ostream& lhs, G4QParticle& rhs)
130 {
131  G4QPDGCode rhsQPDG = rhs.GetQPDG();
132  lhs << G4endl << "Particle with PDG=" << rhsQPDG << ", Spin=" << rhs.GetSpin()
133  << ", mass=" << rhs.GetMass() << ", width=" << rhs.GetWidth() << G4endl;
134  lhs<<" Quark Content of the Particle="<<rhs.GetQContent()<<", Decay Channels:"<<G4endl;
135  G4QDecayChanVector DCV = rhs.GetDecayVector();
136  G4int n = DCV.size();
137  for (int i=0; i<n; i++)
138  {
139  lhs << DCV[i]->GetDecayChanLimit() << "PDG codes";
140  G4QPDGCodeVector PCV=DCV[i]->GetVecOfSecHadrons();
141  G4int m_value = PCV.size();
142  for (int j=0; j<m_value; j++)
143  {
144  if(!j) lhs << ":";
145  else lhs << ",";
146  lhs << PCV[j]->GetPDGCode() ;
147  }
148  }
149  return lhs;
150 }
151 
152 // Initialize the PDG-Particle by QCode @@ Can be improved, using PDG.DATA file
153 G4QDecayChanVector G4QParticle::InitDecayVector(G4int nQ)
154 {
155  //static G4int nP = 486; // Up to A=80
156  //static const G4int nP = 494; // Up to A=80(?) "Isonuclear revision"
157  static const G4int nP = 512; // A<57 "Leptons/Hypernuclei" G4QCHIPSWorld::GetParticles(!)
158  //static G4QDecayChanVector* DecayDB = new G4QDecayChanVector[nP];
159  static G4QDecayChanVector DecayDB[nP];
160  static int limit= 0;
161  if(nQ>=limit && nQ<nP)
162  {
163  //*** Secondary PDG-particles should be ordered in a Channel by increasing width***!***
164  //** Channels should be ordered by increasing minimum mass of the secondary particles**
165  //if(limit<= 0 && nQ>= 0)DecayDB[ 0] = 0; // e- (11)
166  //if(limit<= 1 && nQ>= 1)DecayDB[ 1] = 0; // nu_e (12)
167  //if(limit<= 2 && nQ>= 2)DecayDB[ 2] = 0; // mu- (13)
168  //if(limit<= 3 && nQ>= 3)DecayDB[ 3] = 0; // mu_e (14)
169  //if(limit<= 4 && nQ>= 4)DecayDB[ 4] = 0; // tau- (15)
170  //if(limit<= 5 && nQ>= 5)DecayDB[ 5] = 0; // nu_tau (16)
171  //if(limit<= 6 && nQ>= 6)DecayDB[ 6] = 0; // gamma (22)
172  if(limit<= 7 && nQ>= 7) // Z0 (23)
173  {
174  DecayDB[ 7].push_back(new G4QDecayChan(.036, 11, -11));
175  DecayDB[ 7].push_back(new G4QDecayChan(.073, 13, -13));
176  DecayDB[ 7].push_back(new G4QDecayChan(.107, 15, -15));
177  DecayDB[ 7].push_back(new G4QDecayChan(.174, 12, -12)); // @@ Fake invisible decays
178  DecayDB[ 7].push_back(new G4QDecayChan(.240, 14, -14));
179  DecayDB[ 7].push_back(new G4QDecayChan(.307, 16, -16));
180  DecayDB[ 7].push_back(new G4QDecayChan(.400,2112,-2112)); // @@ Fake Hadronic decays
181  DecayDB[ 7].push_back(new G4QDecayChan(.500,2212,-2212)); // @@ Need heavy quarks
182  DecayDB[ 7].push_back(new G4QDecayChan(.600,2212,-2212, 111));
183  DecayDB[ 7].push_back(new G4QDecayChan(.700,2112,-2112, 111));
184  DecayDB[ 7].push_back(new G4QDecayChan(.800,2212,-2112,-211));
185  DecayDB[ 7].push_back(new G4QDecayChan(.990,2112,-2212, 211));
186  DecayDB[ 7].push_back(new G4QDecayChan(1.00,2112,-3122, 111));
187  }
188  if(limit<= 8 && nQ>= 8) // W- (24) @@ Update HadronicDecays
189  {
190  DecayDB[ 8].push_back(new G4QDecayChan(.107, 11, -12));
191  DecayDB[ 8].push_back(new G4QDecayChan(.214, 13, -14));
192  DecayDB[ 8].push_back(new G4QDecayChan(.321, 15, -16));
193  DecayDB[ 8].push_back(new G4QDecayChan(.421,2112,-2212)); // @@ Fake Hadronic decays
194  DecayDB[ 8].push_back(new G4QDecayChan(.521,2112,-2112,-211));
195  DecayDB[ 8].push_back(new G4QDecayChan(.621,2212,-2212,-211));
196  DecayDB[ 8].push_back(new G4QDecayChan(.721,3122,-3122,-211));
197  DecayDB[ 8].push_back(new G4QDecayChan(.821,2112,-2212, 111));
198  DecayDB[ 8].push_back(new G4QDecayChan(.921,3122,-2212, 111));
199  DecayDB[ 8].push_back(new G4QDecayChan(1.00,2112,-3122,-211));
200  }
201  //if(limit<= 9 && nQ>= 9)DecayDB[ 9] = 0; // H0 (25)
202  //if(limit<= 10 && nQ>= 10)DecayDB[ 10] = 0; // H- (37)
203  if(limit<= 11 && nQ>= 11) // Low sigma=pi,pi S-wave : f_0 (800)
204  {
205  DecayDB[ 11].push_back(new G4QDecayChan(.333,211,-211));
206  DecayDB[ 11].push_back(new G4QDecayChan(1.00,111, 111));
207  }
208  if(limit<= 12 && nQ>= 12) // Midle Regeon-Pomeron : f_0 (980)
209  {
210  DecayDB[ 12].push_back(new G4QDecayChan(.333,211,-211));
211  DecayDB[ 12].push_back(new G4QDecayChan(1.00,111, 111));
212  }
213  if(limit<= 13 && nQ>= 13) // High Regeon-Pomeron : f_0 (1500)
214  {
215  DecayDB[ 13].push_back(new G4QDecayChan(.019,221, 331));
216  DecayDB[ 13].push_back(new G4QDecayChan(.070,221, 221));
217  DecayDB[ 13].push_back(new G4QDecayChan(.113,311,-311));
218  DecayDB[ 13].push_back(new G4QDecayChan(.156,321,-321));
219  DecayDB[ 13].push_back(new G4QDecayChan(.578,211,-211)); //@@ include 4pi decays
220  DecayDB[ 13].push_back(new G4QDecayChan(1.00,111, 111));
221  }
222  //if(limit<= 14 && nQ>= 14)DecayDB[ 14].push_back(new G4QDecayChan(1.00,22,22));//Pi0
223  //if(limit<= 15 && nQ>= 15)DecayDB[ 15] = 0; // Pi +
224  if(limit<= 16 && nQ>= 16) // eta
225  {
226  DecayDB[ 16].push_back(new G4QDecayChan(.226,211,-211,111));
227  DecayDB[ 16].push_back(new G4QDecayChan(.551,111, 111,111));
228  DecayDB[ 16].push_back(new G4QDecayChan(.598,211,-211, 22));
229  DecayDB[ 16].push_back(new G4QDecayChan(.606, 11, -11, 22)); //@@ .002 (pi+)(pi-)2gam
230  DecayDB[ 16].push_back(new G4QDecayChan(1.00, 22, 22));
231  }
232  //if(limit<= 17 && nQ>= 17) // K 0 (K_short - probab 1/2) @@@@@@@@@@@@
233  //{
234  // DecayDB[ 17].push_back(new G4QDecayChan(.6861,211,-211));
235  // DecayDB[ 17].push_back(new G4QDecayChan(1.00, 111, 111));
236  //}
237  //if(limit<= 18 && nQ>= 18)DecayDB[ 8] = 0; // K +
238  if(limit<= 19 && nQ>= 19) // eta'
239  {
240  DecayDB[ 19].push_back(new G4QDecayChan(.443,211,-211,221));
241  DecayDB[ 19].push_back(new G4QDecayChan(.652,111, 111,221));
242  DecayDB[ 19].push_back(new G4QDecayChan(.947, 22, 223));
243  DecayDB[ 19].push_back(new G4QDecayChan(.949,111, 111,111));
244  DecayDB[ 19].push_back(new G4QDecayChan(.979, 22, 113));
245  DecayDB[ 19].push_back(new G4QDecayChan(1.00, 22, 22));
246  }
247  //if(limit<= 20 && nQ>= 20)DecayDB[ 20] = 0; // n
248  //if(limit<= 21 && nQ>= 21)DecayDB[ 21] = 0; // p
249  //if(limit<= 22 && nQ>= 22) // Lambda =--=>> all week decays are closed at this time
250  //{
251  // DecayDB[ 22].push_back(new G4QDecayChan(.640,2212,-211));
252  // DecayDB[ 22].push_back(new G4QDecayChan(1.00,2112, 111));
253  //}
254  //if(limit<= 23 &&nQ>=23)DecayDB[23].push_back(new G4QDecayChan(1.,2112,-211));//Sigma-
255  if(limit<= 24 &&nQ>=24)DecayDB[24].push_back(new G4QDecayChan(1.,3122,22));//Sigma0(EM)
256  //if(limit<= 25 && nQ>= 25) // Sigma +
257  //{
258  // DecayDB[ 25].push_back(new G4QDecayChan(.484,2112, 211));
259  // DecayDB[ 25].push_back(new G4QDecayChan(1.00,2212, 111));
260  //}
261  //if(limit<= 26 && nQ>=26)DecayDB[26].push_back(new G4QDecayChan(1.,3122,-211));// Ksi-
262  //if(limit<= 27 && nQ>=27)DecayDB[27].push_back(new G4QDecayChan(1.,3122, 111));// Ksi0
263  if(limit<= 28 && nQ>= 28)DecayDB[ 28].push_back(new G4QDecayChan(1., 211,-211));// rho0
264  if(limit<= 29 && nQ>= 29)DecayDB[ 29].push_back(new G4QDecayChan(1., 211, 111));// rho+
265  if(limit<= 30 && nQ>= 30) // omega
266  {
267  DecayDB[ 30].push_back(new G4QDecayChan(.891, 211,-211,111));
268  DecayDB[ 30].push_back(new G4QDecayChan(.908, 211,-211));
269  DecayDB[ 30].push_back(new G4QDecayChan(.997, 22, 111));
270  DecayDB[ 30].push_back(new G4QDecayChan(.998, 11, -11, 111)); //@@NeedsMoreAccurate
271  DecayDB[ 30].push_back(new G4QDecayChan(.998, 13, -13, 111));
272  DecayDB[ 30].push_back(new G4QDecayChan(1.00, 22, 221));
273  }
274  if(limit<= 31 && nQ>= 31) // K* 0
275  {
276  DecayDB[ 31].push_back(new G4QDecayChan(.667,-211, 321));
277  DecayDB[ 31].push_back(new G4QDecayChan(1.00, 111, 311));
278  }
279  if(limit<= 32 && nQ>= 32) // K* +
280  {
281  DecayDB[ 32].push_back(new G4QDecayChan(.667, 211, 311));
282  DecayDB[ 32].push_back(new G4QDecayChan(1.00, 111, 321));
283  }
284  if(limit<= 33 && nQ>= 33) // phi
285  {
286  DecayDB[ 33].push_back(new G4QDecayChan(.491, 311,-311));
287  DecayDB[ 33].push_back(new G4QDecayChan(.831, 321,-321));
288  DecayDB[ 33].push_back(new G4QDecayChan(.844, 22, 221));
289  DecayDB[ 33].push_back(new G4QDecayChan(.846, 22, 111));
290  DecayDB[ 33].push_back(new G4QDecayChan(.897, 211,-213));
291  DecayDB[ 33].push_back(new G4QDecayChan(.948,-211, 213));
292  DecayDB[ 33].push_back(new G4QDecayChan(1.00, 111, 113));
293  }
294  if(limit<= 34 && nQ>= 34)DecayDB[34].push_back(new G4QDecayChan(1.,2112,-211));//Delta-
295  if(limit<= 35 && nQ>= 35) // Delta 0
296  {
297  DecayDB[ 35].push_back(new G4QDecayChan(.333,2212,-211));
298  DecayDB[ 35].push_back(new G4QDecayChan(1.00,2112, 111));
299  }
300  if(limit<= 36 && nQ>= 36) // Delta +
301  {
302  DecayDB[ 36].push_back(new G4QDecayChan(.333,2112, 211));
303  DecayDB[ 36].push_back(new G4QDecayChan(1.00,2212, 111));
304  }
305  if(limit<= 37 && nQ>= 37)DecayDB[37].push_back(new G4QDecayChan(1.,2212,211));//Delta++
306  if(limit<= 38 && nQ>= 38) // Lambda* (1520)
307  {
308  DecayDB[ 38].push_back(new G4QDecayChan(.225,3112,-311));
309  DecayDB[ 38].push_back(new G4QDecayChan(.450,3222,-321));
310  DecayDB[ 38].push_back(new G4QDecayChan(.453,3112,211,111));
311  DecayDB[ 38].push_back(new G4QDecayChan(.456,3212,211,-211));
312  DecayDB[ 38].push_back(new G4QDecayChan(.459,3212,111,111));
313  DecayDB[ 38].push_back(new G4QDecayChan(.462,3222,-211,111));
314  DecayDB[ 38].push_back(new G4QDecayChan(.512,3122,211,-211));
315  DecayDB[ 38].push_back(new G4QDecayChan(.562,3122,111,111));
316  DecayDB[ 38].push_back(new G4QDecayChan(.702,3222,-211));
317  DecayDB[ 38].push_back(new G4QDecayChan(.842,3212, 111));
318  DecayDB[ 38].push_back(new G4QDecayChan(.982,3112, 211));
319  DecayDB[ 38].push_back(new G4QDecayChan(1.00,3122, 22));
320  }
321  if(limit<= 39 && nQ>= 39) // Sigma* -
322  {
323  DecayDB[ 39].push_back(new G4QDecayChan(.060,3112, 111));
324  DecayDB[ 39].push_back(new G4QDecayChan(.120,3212,-211));
325  DecayDB[ 39].push_back(new G4QDecayChan(1.00,3122,-211));
326  }
327  if(limit<= 40 && nQ>= 40) // Sigma* 0
328  {
329  DecayDB[ 40].push_back(new G4QDecayChan(.040,3112, 211));
330  DecayDB[ 40].push_back(new G4QDecayChan(.080,3222,-211));
331  DecayDB[ 40].push_back(new G4QDecayChan(.120,3212, 111));
332  DecayDB[ 40].push_back(new G4QDecayChan(1.00,3122, 111));
333  }
334  if(limit<= 41 && nQ>= 41) // Sigma* +
335  {
336  DecayDB[ 41].push_back(new G4QDecayChan(.060,3212, 211));
337  DecayDB[ 41].push_back(new G4QDecayChan(.120,3222, 111));
338  DecayDB[ 41].push_back(new G4QDecayChan(1.00,3122, 211));
339  }
340  if(limit<= 42 && nQ>= 42) // Ksi* -
341  {
342  DecayDB[ 42].push_back(new G4QDecayChan(.667,3322,-211));
343  DecayDB[ 42].push_back(new G4QDecayChan(1.00,3312, 111));
344  }
345  if(limit<= 43 && nQ>= 43) // Ksi* 0
346  {
347  DecayDB[ 43].push_back(new G4QDecayChan(.667,3312, 211));
348  DecayDB[ 43].push_back(new G4QDecayChan(1.00,3322, 111));
349  }
350  //if(limit<= 44 && nQ>= 44) // OMEGA - (Weak)
351  //{
352  // DecayDB[ 44].push_back(new G4QDecayChan(.678,3122, 321));
353  // DecayDB[ 44].push_back(new G4QDecayChan(.914,3322,-211));
354  // DecayDB[ 44].push_back(new G4QDecayChan(1.00,3312, 111));
355  //}
356  if(limit<= 45 && nQ>= 45) // a_2 0
357  {
358  DecayDB[ 45].push_back(new G4QDecayChan(.070, 211,-211,223));
359  DecayDB[ 45].push_back(new G4QDecayChan(.106, 111, 111,223));
360  DecayDB[ 45].push_back(new G4QDecayChan(.131, 321,-321));
361  DecayDB[ 45].push_back(new G4QDecayChan(.156, 311,-311));
362  DecayDB[ 45].push_back(new G4QDecayChan(.301, 111, 221));
363  DecayDB[ 45].push_back(new G4QDecayChan(.534,-211, 213));
364  DecayDB[ 45].push_back(new G4QDecayChan(.767, 211,-213));
365  DecayDB[ 45].push_back(new G4QDecayChan(1.00, 111, 113));
366  }
367  if(limit<= 46 && nQ>= 46) // a_2 +
368  {
369  DecayDB[ 46].push_back(new G4QDecayChan(.106,111,211,223));
370  DecayDB[ 46].push_back(new G4QDecayChan(.156, 321,-311));
371  DecayDB[ 46].push_back(new G4QDecayChan(.301, 211, 221));
372  DecayDB[ 46].push_back(new G4QDecayChan(.651, 211, 113));
373  DecayDB[ 46].push_back(new G4QDecayChan(1.00, 111, 213));
374  }
375  if(limit<= 47 && nQ>= 47) // f_2 0
376  {
377  DecayDB[ 47].push_back(new G4QDecayChan(.005, 221, 221));
378  DecayDB[ 47].push_back(new G4QDecayChan(.028, 311,-311));
379  DecayDB[ 47].push_back(new G4QDecayChan(.051, 321,-321));
380  DecayDB[ 47].push_back(new G4QDecayChan(.123, 111, 113));
381  DecayDB[ 47].push_back(new G4QDecayChan(.126, 111, 221));
382  DecayDB[ 47].push_back(new G4QDecayChan(.152, 211,-211,113));
383  DecayDB[ 47].push_back(new G4QDecayChan(.717, 211,-211));
384  DecayDB[ 47].push_back(new G4QDecayChan(1.00, 111, 111));
385  }
386  if(limit<= 48 && nQ>= 48) // K_2 0
387  {
388  DecayDB[ 48].push_back(new G4QDecayChan(.028, 311, 223));
389  DecayDB[ 48].push_back(new G4QDecayChan(.074, 211,-211,313));
390  DecayDB[ 48].push_back(new G4QDecayChan(.143,111,-211,323));
391  DecayDB[ 48].push_back(new G4QDecayChan(.166,111, 111,313));
392  DecayDB[ 48].push_back(new G4QDecayChan(.190,-211, 323));
393  DecayDB[ 48].push_back(new G4QDecayChan(.314, 111, 313));
394  DecayDB[ 48].push_back(new G4QDecayChan(.357, 311, 113));
395  DecayDB[ 48].push_back(new G4QDecayChan(.500, 321,-213));
396  DecayDB[ 48].push_back(new G4QDecayChan(.750,-211, 321));
397  DecayDB[ 48].push_back(new G4QDecayChan(1.00, 111, 311));
398  }
399  if(limit<= 49 && nQ>= 49) // K_2 +
400  {
401  DecayDB[ 49].push_back(new G4QDecayChan(.028, 321, 223));
402  DecayDB[ 49].push_back(new G4QDecayChan(.074,211,-211,323));
403  DecayDB[ 49].push_back(new G4QDecayChan(.143,111, 211,313));
404  DecayDB[ 49].push_back(new G4QDecayChan(.166,111, 111,323));
405  DecayDB[ 49].push_back(new G4QDecayChan(.190, 211, 313));
406  DecayDB[ 49].push_back(new G4QDecayChan(.314, 111, 323));
407  DecayDB[ 49].push_back(new G4QDecayChan(.357, 311, 213));
408  DecayDB[ 49].push_back(new G4QDecayChan(.500, 321, 113));
409  DecayDB[ 49].push_back(new G4QDecayChan(.750, 211, 311));
410  DecayDB[ 49].push_back(new G4QDecayChan(1.00, 111, 321));
411  }
412  if(limit<= 50 && nQ>= 50) // f_2' 0
413  {
414  DecayDB[ 50].push_back(new G4QDecayChan(.103, 221, 221));
415  DecayDB[ 50].push_back(new G4QDecayChan(.547, 311,-311));
416  DecayDB[ 50].push_back(new G4QDecayChan(.991, 321,-321));
417  DecayDB[ 50].push_back(new G4QDecayChan(.997, 211,-211));
418  DecayDB[ 50].push_back(new G4QDecayChan(1.00, 111, 111));
419  }
420  if(limit<= 51 && nQ>= 51) // N_5/2 0
421  {
422  DecayDB[ 51].push_back(new G4QDecayChan(.040, 211, 1114));
423  DecayDB[ 51].push_back(new G4QDecayChan(.080, 111, 2114));
424  DecayDB[ 51].push_back(new G4QDecayChan(.120,-211, 2214));
425  DecayDB[ 51].push_back(new G4QDecayChan(.180, 2112, 113));
426  DecayDB[ 51].push_back(new G4QDecayChan(.210, 2212,-213));
427  DecayDB[ 51].push_back(new G4QDecayChan(.340, 2112, 110));
428  DecayDB[ 51].push_back(new G4QDecayChan(.780, 2212,-211));
429  DecayDB[ 51].push_back(new G4QDecayChan(1.00, 2112, 111));
430  }
431  if(limit<= 52 && nQ>= 52) // N_5/2 +
432  {
433  DecayDB[ 52].push_back(new G4QDecayChan(.040,-211, 2224));
434  DecayDB[ 52].push_back(new G4QDecayChan(.080, 211, 2114));
435  DecayDB[ 52].push_back(new G4QDecayChan(.120, 111, 2214));
436  DecayDB[ 52].push_back(new G4QDecayChan(.180, 2112, 213));
437  DecayDB[ 52].push_back(new G4QDecayChan(.210, 2212, 113));
438  DecayDB[ 52].push_back(new G4QDecayChan(.340, 2212, 229));
439  DecayDB[ 52].push_back(new G4QDecayChan(.780, 2112, 211));
440  DecayDB[ 52].push_back(new G4QDecayChan(1.00, 2212, 111));
441  }
442  if(limit<= 53 && nQ>= 53) // LAMBDA_5/2
443  {
444  DecayDB[ 53].push_back(new G4QDecayChan(.350, 2112,-311));
445  DecayDB[ 53].push_back(new G4QDecayChan(.700, 2212,-321));
446  DecayDB[ 53].push_back(new G4QDecayChan(.740, 211, 3114));
447  DecayDB[ 53].push_back(new G4QDecayChan(.780,-211, 3224));
448  DecayDB[ 53].push_back(new G4QDecayChan(.820, 111, 3214));
449  DecayDB[ 53].push_back(new G4QDecayChan(.880, 3112, 211));
450  DecayDB[ 53].push_back(new G4QDecayChan(.940, 3222,-211));
451  DecayDB[ 53].push_back(new G4QDecayChan(1.00, 3212, 111));
452  }
453  if(limit<= 54 && nQ>= 54) // SIGMA_5/2 -
454  {
455  DecayDB[ 54].push_back(new G4QDecayChan(.600, 2112,-321));
456  DecayDB[ 54].push_back(new G4QDecayChan(.660,-211, 3214));
457  DecayDB[ 54].push_back(new G4QDecayChan(.720, 111, 3114));
458  DecayDB[ 54].push_back(new G4QDecayChan(.810, 3212,-211));
459  DecayDB[ 54].push_back(new G4QDecayChan(.900, 3112, 111));
460  DecayDB[ 54].push_back(new G4QDecayChan(1.00, 3122,-211));
461  }
462  if(limit<= 55 && nQ>= 55) // SIGMA_5/2 0
463  {
464  DecayDB[ 55].push_back(new G4QDecayChan(.300, 2112,-311));
465  DecayDB[ 55].push_back(new G4QDecayChan(.600, 2212,-321));
466  DecayDB[ 55].push_back(new G4QDecayChan(.640, 211, 3114));
467  DecayDB[ 55].push_back(new G4QDecayChan(.680,-211, 3224));
468  DecayDB[ 55].push_back(new G4QDecayChan(.720, 111, 3214));
469  DecayDB[ 55].push_back(new G4QDecayChan(.780, 3112, 211));
470  DecayDB[ 55].push_back(new G4QDecayChan(.840, 3222,-211));
471  DecayDB[ 55].push_back(new G4QDecayChan(.900, 3212, 111));
472  DecayDB[ 55].push_back(new G4QDecayChan(1.00, 3122, 111));
473  }
474  if(limit<= 56 && nQ>= 56) // SIGMA_5/2 +
475  {
476  DecayDB[ 56].push_back(new G4QDecayChan(.600, 2212,-311));
477  DecayDB[ 56].push_back(new G4QDecayChan(.660, 211, 3214));
478  DecayDB[ 56].push_back(new G4QDecayChan(.720, 111, 3224));
479  DecayDB[ 56].push_back(new G4QDecayChan(.810, 3212, 211));
480  DecayDB[ 56].push_back(new G4QDecayChan(.900, 3222, 111));
481  DecayDB[ 56].push_back(new G4QDecayChan(1.00, 3122, 211));
482  }
483  if(limit<= 57 && nQ>= 57) // KSI_5/2 -
484  {
485  DecayDB[ 57].push_back(new G4QDecayChan(.400, 3112,-311));
486  DecayDB[ 57].push_back(new G4QDecayChan(.800, 3212,-321));
487  DecayDB[ 57].push_back(new G4QDecayChan(1.00, 3122,-321));
488  }
489  if(limit<= 58 && nQ>= 58) // KSI_5/2 0
490  {
491  DecayDB[ 58].push_back(new G4QDecayChan(.400, 3212,-311));
492  DecayDB[ 58].push_back(new G4QDecayChan(.800, 3222,-321));
493  DecayDB[ 58].push_back(new G4QDecayChan(1.00, 3122,-311));
494  }
495  if(limit<= 59 && nQ>= 59) // rho_3 0
496  {
497  DecayDB[ 59].push_back(new G4QDecayChan(.019,311,-313));
498  DecayDB[ 59].push_back(new G4QDecayChan(.038,321,-323));
499  DecayDB[ 59].push_back(new G4QDecayChan(.046,311,-311));
500  DecayDB[ 59].push_back(new G4QDecayChan(.054,321,-321));
501  DecayDB[ 59].push_back(new G4QDecayChan(.224,111, 223));
502  DecayDB[ 59].push_back(new G4QDecayChan(.404,111,-211,213));
503  DecayDB[ 59].push_back(new G4QDecayChan(.584,111, 211,-213));
504  DecayDB[ 59].push_back(new G4QDecayChan(.764,111, 111,113));
505  DecayDB[ 59].push_back(new G4QDecayChan(1.00,211,-211));
506  }
507  if(limit<= 60 && nQ>= 60) // rho_3 +
508  {
509  DecayDB[ 60].push_back(new G4QDecayChan(.019, 321,-313));
510  DecayDB[ 60].push_back(new G4QDecayChan(.038,-311, 323));
511  DecayDB[ 60].push_back(new G4QDecayChan(.054, 321,-311));
512  DecayDB[ 60].push_back(new G4QDecayChan(.224, 211, 223));
513  DecayDB[ 60].push_back(new G4QDecayChan(.404,211,-211,213));
514  DecayDB[ 60].push_back(new G4QDecayChan(.584,211,211,-213));
515  DecayDB[ 60].push_back(new G4QDecayChan(.764,211,111,113));
516  DecayDB[ 60].push_back(new G4QDecayChan(1.00, 211, 111));
517  }
518  if(limit<= 61 && nQ>= 61) // omega_3
519  {
520  DecayDB[ 61].push_back(new G4QDecayChan(.020,211,-211,223));
521  DecayDB[ 61].push_back(new G4QDecayChan(.040,111, 111,223));
522  DecayDB[ 61].push_back(new G4QDecayChan(.060, 211,-213));
523  DecayDB[ 61].push_back(new G4QDecayChan(.080,-211, 213));
524  DecayDB[ 61].push_back(new G4QDecayChan(1.00, 111, 113));
525  }
526  if(limit<= 62 && nQ>= 62) // K_3 0
527  {
528  DecayDB[ 62].push_back(new G4QDecayChan(.030, 111, 315));
529  DecayDB[ 62].push_back(new G4QDecayChan(.060,-211, 325));
530  DecayDB[ 62].push_back(new G4QDecayChan(.340, 311, 331));
531  DecayDB[ 62].push_back(new G4QDecayChan(.400, 111, 313));
532  DecayDB[ 62].push_back(new G4QDecayChan(.520,-211, 323));
533  DecayDB[ 62].push_back(new G4QDecayChan(.620, 311, 113));
534  DecayDB[ 62].push_back(new G4QDecayChan(.820, 321,-213));
535  DecayDB[ 62].push_back(new G4QDecayChan(.940,-211, 321));
536  DecayDB[ 62].push_back(new G4QDecayChan(1.00, 111, 311));
537  }
538  if(limit<= 63 && nQ>= 63) // K_3 +
539  {
540  DecayDB[ 63].push_back(new G4QDecayChan(.030, 211, 315));
541  DecayDB[ 63].push_back(new G4QDecayChan(.060, 111, 325));
542  DecayDB[ 63].push_back(new G4QDecayChan(.340, 321, 331));
543  DecayDB[ 63].push_back(new G4QDecayChan(.400, 211, 313));
544  DecayDB[ 63].push_back(new G4QDecayChan(.520, 111, 323));
545  DecayDB[ 63].push_back(new G4QDecayChan(.620, 311, 213));
546  DecayDB[ 63].push_back(new G4QDecayChan(.820, 321, 113));
547  DecayDB[ 63].push_back(new G4QDecayChan(.940, 211, 311));
548  DecayDB[ 63].push_back(new G4QDecayChan(1.00, 111, 321));
549  }
550  if(limit<= 64 && nQ>= 64) // phi_3
551  {
552  DecayDB[ 64].push_back(new G4QDecayChan(.250, 321,-321));
553  DecayDB[ 64].push_back(new G4QDecayChan(.500, 311,-311));
554  DecayDB[ 64].push_back(new G4QDecayChan(.625, 321,-323));
555  DecayDB[ 64].push_back(new G4QDecayChan(.750,-321, 323));
556  DecayDB[ 64].push_back(new G4QDecayChan(.875, 311,-313));
557  DecayDB[ 64].push_back(new G4QDecayChan(1.00,-311, 313));
558  }
559  if(limit<= 65 && nQ>= 65) // DELTA_7/2 -
560  {
561  DecayDB[ 65].push_back(new G4QDecayChan(.200, 2112,-213 ));
562  DecayDB[ 65].push_back(new G4QDecayChan(.320,-211 , 2114));
563  DecayDB[ 65].push_back(new G4QDecayChan(.500, 111 , 1114));
564  DecayDB[ 65].push_back(new G4QDecayChan(1.00, 2112,-211 ));
565  }
566  if(limit<= 66 && nQ>= 66) // DELTA_7/2 0
567  {
568  DecayDB[ 66].push_back(new G4QDecayChan(.133, 2112, 113 ));
569  DecayDB[ 66].push_back(new G4QDecayChan(.200, 2212,-213 ));
570  DecayDB[ 66].push_back(new G4QDecayChan(.360,-211 , 2214));
571  DecayDB[ 66].push_back(new G4QDecayChan(.480, 211 , 1114));
572  DecayDB[ 66].push_back(new G4QDecayChan(.500, 111 , 2114));
573  DecayDB[ 66].push_back(new G4QDecayChan(.666, 2212,-211 ));
574  DecayDB[ 66].push_back(new G4QDecayChan(1.00, 2112, 111 ));
575  }
576  if(limit<= 67 && nQ>= 67) // DELTA_7/2 +
577  {
578  DecayDB[ 67].push_back(new G4QDecayChan(.133, 2112, 213 ));
579  DecayDB[ 67].push_back(new G4QDecayChan(.200, 2212, 113 ));
580  DecayDB[ 67].push_back(new G4QDecayChan(.360,-211 , 2224));
581  DecayDB[ 67].push_back(new G4QDecayChan(.480, 211 , 2114));
582  DecayDB[ 67].push_back(new G4QDecayChan(.500, 111 , 2214));
583  DecayDB[ 67].push_back(new G4QDecayChan(.666, 2112, 211 ));
584  DecayDB[ 67].push_back(new G4QDecayChan(1.00, 2212, 111 ));
585  }
586  if(limit<= 68 && nQ>= 68) // DELTA_7/2 ++
587  {
588  DecayDB[ 68].push_back(new G4QDecayChan(.200, 2212, 213 ));
589  DecayDB[ 68].push_back(new G4QDecayChan(.320, 211 , 2214));
590  DecayDB[ 68].push_back(new G4QDecayChan(.500, 111 , 2224));
591  DecayDB[ 68].push_back(new G4QDecayChan(1.00, 2212, 211 ));
592  }
593  if(limit<= 69 && nQ>= 69) // LAMBDA_7/2
594  {
595  DecayDB[ 69].push_back(new G4QDecayChan(.160, 3122, 223 ));
596  DecayDB[ 69].push_back(new G4QDecayChan(.260, 2112,-313 ));
597  DecayDB[ 69].push_back(new G4QDecayChan(.360, 2212,-323 ));
598  DecayDB[ 69].push_back(new G4QDecayChan(.400, 3312, 321 ));
599  DecayDB[ 69].push_back(new G4QDecayChan(.440, 3322, 311 ));
600  DecayDB[ 69].push_back(new G4QDecayChan(.480, 3122, 221 ));
601  DecayDB[ 69].push_back(new G4QDecayChan(.520, 2112,-311 ));
602  DecayDB[ 69].push_back(new G4QDecayChan(.560, 2212,-321 ));
603  DecayDB[ 69].push_back(new G4QDecayChan(.600, 3112, 211 ));
604  DecayDB[ 69].push_back(new G4QDecayChan(.800, 3222,-211 ));
605  DecayDB[ 69].push_back(new G4QDecayChan(1.00, 3212, 111 ));
606  }
607  if(limit<= 70 && nQ>= 70) // SIGMA_7/2 -
608  {
609  DecayDB[ 70].push_back(new G4QDecayChan(.030, 2112,-323 ));
610  DecayDB[ 70].push_back(new G4QDecayChan(.165,-311 , 1114));
611  DecayDB[ 70].push_back(new G4QDecayChan(.210,-321 , 2114));
612  DecayDB[ 70].push_back(new G4QDecayChan(.390,-211 , 3124));
613  DecayDB[ 70].push_back(new G4QDecayChan(.450,-211 , 3214));
614  DecayDB[ 70].push_back(new G4QDecayChan(.510, 111 , 3114));
615  DecayDB[ 70].push_back(new G4QDecayChan(.540, 311 , 3314));
616  DecayDB[ 70].push_back(new G4QDecayChan(.570,-211 , 3212));
617  DecayDB[ 70].push_back(new G4QDecayChan(.600, 111 , 3112));
618  DecayDB[ 70].push_back(new G4QDecayChan(.780,-321 , 2112));
619  DecayDB[ 70].push_back(new G4QDecayChan(1.00,-211 , 3122));
620  }
621  if(limit<= 71 && nQ>= 71) // SIGMA_7/2 0
622  {
623  DecayDB[ 71].push_back(new G4QDecayChan(.015, 2112,-313 ));
624  DecayDB[ 71].push_back(new G4QDecayChan(.030, 2212,-321 ));
625  DecayDB[ 71].push_back(new G4QDecayChan(.120,-311 , 2114));
626  DecayDB[ 71].push_back(new G4QDecayChan(.210,-321 , 2214));
627  DecayDB[ 71].push_back(new G4QDecayChan(.390, 111 , 3124));
628  DecayDB[ 71].push_back(new G4QDecayChan(.450,-211 , 3224));
629  DecayDB[ 71].push_back(new G4QDecayChan(.510, 211 , 3114));
630  DecayDB[ 71].push_back(new G4QDecayChan(.525, 311 , 3324));
631  DecayDB[ 71].push_back(new G4QDecayChan(.540, 321 , 3314));
632  DecayDB[ 71].push_back(new G4QDecayChan(.570,-211 , 3222));
633  DecayDB[ 71].push_back(new G4QDecayChan(.600, 211 , 3112));
634  DecayDB[ 71].push_back(new G4QDecayChan(.690,-311 , 2112));
635  DecayDB[ 71].push_back(new G4QDecayChan(.780,-321 , 2212));
636  DecayDB[ 71].push_back(new G4QDecayChan(1.00, 111 , 3122));
637  }
638  if(limit<= 72 && nQ>= 72) // SIGMA_7/2 +
639  {
640  DecayDB[ 72].push_back(new G4QDecayChan(.030, 2212,-313 ));
641  DecayDB[ 72].push_back(new G4QDecayChan(.165,-321 , 2224));
642  DecayDB[ 72].push_back(new G4QDecayChan(.210,-311 , 2214));
643  DecayDB[ 72].push_back(new G4QDecayChan(.390, 211 , 3124));
644  DecayDB[ 72].push_back(new G4QDecayChan(.450, 211 , 3214));
645  DecayDB[ 72].push_back(new G4QDecayChan(.510, 111 , 3224));
646  DecayDB[ 72].push_back(new G4QDecayChan(.540, 321 , 3324));
647  DecayDB[ 72].push_back(new G4QDecayChan(.570, 211 , 3212));
648  DecayDB[ 72].push_back(new G4QDecayChan(.600, 111 , 3222));
649  DecayDB[ 72].push_back(new G4QDecayChan(.780,-311 , 2212));
650  DecayDB[ 72].push_back(new G4QDecayChan(1.00, 211 , 3122));
651  }
652  if(limit<= 73 && nQ>= 73) // KSI_7/2 -
653  {
654  DecayDB[ 73].push_back(new G4QDecayChan(.400, 3112,-311));
655  DecayDB[ 73].push_back(new G4QDecayChan(.800, 3212,-321));
656  DecayDB[ 73].push_back(new G4QDecayChan(1.00, 3122,-321));
657  }
658  if(limit<= 74 && nQ>= 74) // KSI_7/2 0
659  {
660  DecayDB[ 74].push_back(new G4QDecayChan(.400, 3212,-311));
661  DecayDB[ 74].push_back(new G4QDecayChan(.800, 3222,-321));
662  DecayDB[ 74].push_back(new G4QDecayChan(1.00, 3122,-311));
663  }
664  if(limit<= 75 && nQ>= 75) // OMEGA_7/2 -
665  {
666  DecayDB[ 75].push_back(new G4QDecayChan(.250,-311 , 3314));
667  DecayDB[ 75].push_back(new G4QDecayChan(.500,-321 , 3324));
668  DecayDB[ 75].push_back(new G4QDecayChan(.750, 3312,-313 ));
669  DecayDB[ 75].push_back(new G4QDecayChan(1.00, 3322,-323 ));
670  }
671  if(limit<= 76 && nQ>= 76) // a_4 0
672  {
673  DecayDB[ 76].push_back(new G4QDecayChan(.200, 311,-311));
674  DecayDB[ 76].push_back(new G4QDecayChan(.400, 321,-321));
675  DecayDB[ 76].push_back(new G4QDecayChan(.600, 111, 221));
676  DecayDB[ 76].push_back(new G4QDecayChan(1.00, 111, 211,-211));
677  }
678  if(limit<= 77 && nQ>= 77) // a_4 +
679  {
680  DecayDB[ 77].push_back(new G4QDecayChan(.400, 321,-311));
681  DecayDB[ 77].push_back(new G4QDecayChan(.600, 211, 221));
682  DecayDB[ 77].push_back(new G4QDecayChan(.800, 211, 211,-211));
683  DecayDB[ 77].push_back(new G4QDecayChan(1.00, 211, 111, 111));
684  }
685  if(limit<= 78 && nQ>= 78) // f_4 0
686  {
687  DecayDB[ 78].push_back(new G4QDecayChan(.020, 333, 333));
688  DecayDB[ 78].push_back(new G4QDecayChan(.340, 223, 223));
689  DecayDB[ 78].push_back(new G4QDecayChan(.350, 221, 221));
690  DecayDB[ 78].push_back(new G4QDecayChan(.360, 311,-311));
691  DecayDB[ 78].push_back(new G4QDecayChan(.370, 321,-321));
692  DecayDB[ 78].push_back(new G4QDecayChan(.670, 213,-213));
693  DecayDB[ 78].push_back(new G4QDecayChan(.820, 113, 113));
694  DecayDB[ 78].push_back(new G4QDecayChan(.940, 211,-211));
695  DecayDB[ 78].push_back(new G4QDecayChan(1.00, 111, 111));
696  }
697  if(limit<= 79 && nQ>= 79) // K_4 0
698  {
699  DecayDB[ 79].push_back(new G4QDecayChan(.060, 333, 313));
700  DecayDB[ 79].push_back(new G4QDecayChan(.260, 223, 313));
701  DecayDB[ 79].push_back(new G4QDecayChan(.380, 313, 113));
702  DecayDB[ 79].push_back(new G4QDecayChan(.400, 323,-213));
703  DecayDB[ 79].push_back(new G4QDecayChan(.500, 223, 311));
704  DecayDB[ 79].push_back(new G4QDecayChan(.550, 311, 113));
705  DecayDB[ 79].push_back(new G4QDecayChan(.600, 321,-213));
706  DecayDB[ 79].push_back(new G4QDecayChan(.700, 311, 113));
707  DecayDB[ 79].push_back(new G4QDecayChan(.800, 321,-213));
708  DecayDB[ 79].push_back(new G4QDecayChan(.900, 311, 111));
709  DecayDB[ 79].push_back(new G4QDecayChan(1.00, 321,-211));
710  }
711  if(limit<= 80 && nQ>= 80) // K_4 +
712  {
713  DecayDB[ 80].push_back(new G4QDecayChan(.060, 333, 323));
714  DecayDB[ 80].push_back(new G4QDecayChan(.260, 223, 323));
715  DecayDB[ 80].push_back(new G4QDecayChan(.380, 313, 213));
716  DecayDB[ 80].push_back(new G4QDecayChan(.400, 323, 113));
717  DecayDB[ 80].push_back(new G4QDecayChan(.500, 223, 321));
718  DecayDB[ 80].push_back(new G4QDecayChan(.550, 321, 113));
719  DecayDB[ 80].push_back(new G4QDecayChan(.600, 311, 213));
720  DecayDB[ 80].push_back(new G4QDecayChan(.700, 311, 211));
721  DecayDB[ 80].push_back(new G4QDecayChan(.800, 321, 111));
722  DecayDB[ 80].push_back(new G4QDecayChan(.900, 311, 211));
723  DecayDB[ 80].push_back(new G4QDecayChan(1.00, 321, 111));
724  }
725  if(limit<=81&&nQ>=81)DecayDB[81].push_back(new G4QDecayChan(1., 333,333));//phi_4(2300)
726  if(limit<=82&&nQ>=82)DecayDB[82].push_back(new G4QDecayChan(1.,2212, 2224));//pDelta++
727  if(limit<=83&&nQ>=83)DecayDB[83].push_back(new G4QDecayChan(1.,2112, 1114));//nDelta-
728  if(limit<=84&&nQ>=84)DecayDB[84].push_back(new G4QDecayChan(1.,2224, 2224));//D++D++
729  if(limit<=85&&nQ>=85)DecayDB[85].push_back(new G4QDecayChan(1.,1114, 1114));//Del-Del-
730  if(limit<=86&&nQ>=86)DecayDB[86].push_back(new G4QDecayChan(1.,2212,2212,2224));//ppD++
731  if(limit<=87&&nQ>=87)DecayDB[87].push_back(new G4QDecayChan(1.,2112,2112,1114));//nnD-
732  if(limit<=88&&nQ>=88)DecayDB[88].push_back(new G4QDecayChan(1.,2212,2224,2224));//p2D++
733  if(limit<=89&&nQ>=89)DecayDB[89].push_back(new G4QDecayChan(1.,2112,1114,1114));//nD-D-
734  //if(limit<=90&&nQ>=90)DecayDB[90] = 0; // neutron (n) as a quark exchange fragment
735  //if(limit<=91&&nQ>=91)DecayDB[91] = 0; // proton (p) as a quark exchange fragment
736  //if(limit<=92&&nQ>=92)DecayDB[92] = 0; // neutron (L/Sigma0) as aQuarkExchangeFragment
737  //if(limit<=93&&nQ>=93)DecayDB[93] = 0; // neutron (Sigma-) as a quarkExchangeFragment
738  //if(limit<=94&&nQ>=94)DecayDB[94] = 0; // neutron (Sigma+) as a quarkExchangeFragment
739  //if(limit<=95&&nQ>=95)DecayDB[95] = 0; // neutron (Xi-) as a quark exchange fragment
740  //if(limit<=96&&nQ>=96)DecayDB[96] = 0; // neutron (Xi0) as a quark exchange fragment
741  //if(limit<=97&&nQ>=97)DecayDB[97] = 0; // neutron (Omega-) as a quarkExchangeFragment
742  if(limit<=98&&nQ>=98)DecayDB[98].push_back(new G4QDecayChan(1.,2112, 2112)); //nn (src)
743  if(limit<=99&&nQ>=99)DecayDB[99].push_back(new G4QDecayChan(1.,2212, 2112)); //d/pn(?)
744  if(limit<=100&&nQ>=100)DecayDB[100].push_back(new G4QDecayChan(1.,2212,2212));//pp(src)
745  if(limit<=101&&nQ>=101)DecayDB[101].push_back(new G4QDecayChan(1.,3122,2112));//Ln
746  if(limit<=102&&nQ>=102)DecayDB[102].push_back(new G4QDecayChan(1.,3122,2212));//Lp
747  if(limit<=104&&nQ>=104)DecayDB[104].push_back(new G4QDecayChan(1.,3112,2112));//nSig-
748  if(limit<=103&&nQ>=103)DecayDB[103].push_back(new G4QDecayChan(1.,3122,3122));//LL
749  if(limit<=105&&nQ>=105)DecayDB[105].push_back(new G4QDecayChan(1.,3222,2212));//pSig+
750  //if(limit<=106&&nQ>=106)DecayDB[106] = 0; // t
751  //if(limit<=107&&nQ>=107)DecayDB[107] = 0; // He3
752  //Lnn=>Lambda+n+n decay to avoid the final state Hypernucleus
753  if(limit<=108&&nQ>=108)DecayDB[108].push_back(new G4QDecayChan(1.,3122,2112,2112));
754  if(limit<=109&&nQ>=109)DecayDB[109].push_back(new G4QDecayChan(1.,3122,90001001));// Ld
755  //Lpp=>Lambda+p+p decay to avoid the final state Hypernucleus
756  if(limit<=110&&nQ>=110)DecayDB[110].push_back(new G4QDecayChan(1.,3122,2212,2212));
757  //LLn=>Lambda+Lambda+n decay to avoid the final state Hypernucleus
758  if(limit<=111&&nQ>=111)DecayDB[111].push_back(new G4QDecayChan(1.,3122,3122,2112));
759  //LLp=>Lambda+Lambda+p decay to avoid the final state Hypernucleus
760  if(limit<=112&&nQ>=112)DecayDB[112].push_back(new G4QDecayChan(1.,3122,3122,2212));
761  // nnSigma-=>n+n+Sigma- decay to avoid the final state Hypernucleus
762  if(limit<=113&&nQ>=113)DecayDB[113].push_back(new G4QDecayChan(1.,2112,2112,3112));
763  // ------- Nuclear fragments
764  //if(limit<=114 && nQ>=114)
765  //{
766  // if(limit<114) limit=101;
767  // for (int i=limit; i<nQ; i++) DecayDB[i] = 0;
768  //}
769  //Update the limit
770  limit=nQ+1;
771 #ifdef debug
772  G4cout<<"G4QParticle::InitDecayVector: limit is set to "<<limit<<G4endl;
773 #endif
774  }
775  //if(!nQ)G4cout<<"G4QParticle::InitDecayVector:Q=0,nD="<<DecayDB[abs(nQ)].size()<<G4endl;
776  return DecayDB[std::abs(nQ)];
777 }
778 
779 // Initialize the Particle by a Q Code
780 void G4QParticle::InitQParticle(G4int theQCode)
781 {
782  aQPDG.InitByQCode(theQCode);
783  aQuarkCont = aQPDG.GetQuarkContent();
784  aDecay = InitDecayVector(theQCode);
785  //if(!theQCode)G4cout<<"G4QPar::InitQP:PDG="<<GetPDGCode()<<",n="<<aDecay.size()<<G4endl;
786 }
787 
788 // Initialize the Particle by a PDG Code
789 void G4QParticle::InitPDGParticle(G4int thePDGCode)
790 {
791  aQPDG = G4QPDGCode(thePDGCode);
792  aQuarkCont = aQPDG.GetQuarkContent();
793  aDecay = InitDecayVector(aQPDG.GetQCode());
794 }
Generated on Sat May 25 2013 14:33:47 for Geant4 by   doxygen 1.8.4