ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4ExcitedMesonConstructor.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4ExcitedMesonConstructor.cc
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 //
28 //
29 // --------------------------------------------------------------
30 // GEANT 4 class implementation file
31 // History: first implementation, based on object model of
32 // 10 oct 1998 H.Kurashige
33 // ---------------------------------------------------------------
34 // 01 Oct. 02 Fixed PDG codes for a0(1450), f0(1370), k0_star(1430)
35 
36 #include "G4ExcitedMesonConstructor.hh"
37 
38 #include "G4SystemOfUnits.hh"
39 #include "G4ParticleDefinition.hh"
40 #include "G4ParticleTable.hh"
41 #include "G4PhaseSpaceDecayChannel.hh"
42 #include "G4VDecayChannel.hh"
43 #include "G4DecayTable.hh"
44 
45 
46 G4ExcitedMesonConstructor::G4ExcitedMesonConstructor(G4int ,
47  G4int )
48  : type("meson"), leptonNumber(0), baryonNumber(0)
49 {
50 }
51 
52 G4ExcitedMesonConstructor::~G4ExcitedMesonConstructor()
53 {
54 }
55 
56 void G4ExcitedMesonConstructor::Construct(G4int idx)
57 {
58  G4int iType;
59  if (idx < 0 ) {
60  for (G4int state=0; state< NMultiplets; state +=1) {
61  for (iType = 0; iType < NMesonTypes ; iType++)
62  ConstructMesons(state, iType);
63  }
64  } else if (idx < NMultiplets ) {
65  for (iType = 0; iType < NMesonTypes ; iType++)
66  ConstructMesons(idx, iType);
67  } else {
68 #ifdef G4VERBOSE
69  if (G4ParticleTable::GetParticleTable()->GetVerboseLevel()>1) {
70  G4cerr << "G4ExcitedMesonConstructor::Construct()";
71  G4cerr << " illegal index os state = " << idx << G4endl;
72  }
73 #endif
74  }
75 }
76 
77 G4bool G4ExcitedMesonConstructor::Exist(G4int idxState, G4int idxType)
78 {
79  G4bool value = true;
80  if ( idxType == TEtaPrime ) {
81  if (idxState==N13P0) value = false;
82  if (idxState==N13D1) value = false;
83  } else if ( idxType == TPi ) {
84  if (idxState==N23P2) value = false;
85  }
86  return value;
87 }
88 
89 #include "G4ExcitedMesons.hh"
90 
91 void G4ExcitedMesonConstructor::ConstructMesons(G4int iState, G4int iType)
92 {
93  if (!Exist(iState, iType) ) return;
94 
95  // Construct Resonace particles as dynamic object
96  // Arguments for constructor are as follows
97  // name mass width
98  // charge 2*spin
99  // parity C-conjugation
100  // 2*Isospin 2*Isospin3
101  // G-parity
102  // type lepton number Baryon number
103  // PDG encoding
104  // stable lifetime decay table
105 
106 
107  G4String aName;
108  G4ExcitedMesons* particle;
109 
110  for ( G4int iIso3=(-1)*iIsoSpin[iType]; iIso3<=iIsoSpin[iType]; iIso3+=2) {
111  aName= GetName(iIso3, iState, iType);
112  G4double fmass = mass[iState][iType];
113  G4double fwidth = width[iState][iType];
114  if ( (iType== TK) || (iType==TAntiK ) ) {
115  if ( GetCharge(iIso3,iType) == 0.0) {
116  fmass += massKdiff[iState];
117  fwidth += widthKdiff[iState];
118  }
119  }
120  particle = new G4ExcitedMesons(
121  aName, fmass, fwidth,
122  GetCharge(iIso3,iType), iSpin[iState],
123  iParity[iState], iChargeConjugation[iState],
124  iIsoSpin[iType], iIso3,
125  iGParity[iState][iType],
126  type, leptonNumber, baryonNumber,
127  GetEncoding(iIso3, iState, iType),
128  false, 0.0, NULL
129  );
130 
131  if ( (iType==TEta) || (iType==TEtaPrime) || ((iType==TPi)&&(iIso3==0)) ) {
132  // set same encoding for AntiParticle
133  particle->SetAntiPDGEncoding(GetEncoding(iIso3, iState, iType));
134  }
135  particle->SetMultipletName(name[iState][iType]);
136  particle->SetDecayTable(CreateDecayTable( aName, iIso3, iState, iType));
137  }
138 }
139 
140 
141 G4int G4ExcitedMesonConstructor::GetQuarkContents(G4int iQ,
142  G4int iIso3,
143  G4int iType)
144 {
145  // Quark contents
146 
147  G4int quark=0;
148  if (iType == TPi) {
149  if ( iIso3 == 2 ){
150  if ( iQ == 0 ){ quark = 2; }
151  else { quark = 1; }
152  } else if ( iIso3 == 0 ){
153  quark = 1;
154  } else if ( iIso3 == -2 ){
155  if ( iQ == 0 ){ quark = 1; }
156  else { quark = 2; }
157  }
158  } else if (iType == TEta) {
159  quark = 2;
160 
161  } else if (iType == TEtaPrime) {
162  quark = 3;
163 
164  } else if (iType == TAntiK) {
165  if ( iIso3 == 1 ){
166  if ( iQ == 0 ){ quark = 3; }
167  else { quark = 1; }
168  } else if ( iIso3 == -1 ){
169  if ( iQ == 0 ){ quark = 3; }
170  else { quark = 2; }
171  }
172 
173  } else if (iType == TK) {
174  if ( iIso3 == 1 ){
175  if ( iQ == 0 ){ quark = 2; }
176  else { quark = 3; }
177  } else if ( iIso3 == -1 ){
178  if ( iQ == 0 ){ quark = 1; }
179  else { quark = 3; }
180  }
181 
182  }
183  return quark;
184 }
185 
186 G4double G4ExcitedMesonConstructor::GetCharge(G4int iIsoSpin3, G4int idxType )
187 {
188  static const G4double quark_charge[7] =
189  {
190  0., -1./3., +2./3., -1./3., +2./3., -1./3., +2./3.
191  };
192 
193  G4double charge = quark_charge[GetQuarkContents(0, iIsoSpin3, idxType)]*eplus;
194  charge -= quark_charge[GetQuarkContents(1, iIsoSpin3, idxType)]*eplus;
195  return charge;
196 }
197 
198 G4int G4ExcitedMesonConstructor::GetEncoding(G4int iIsoSpin3,
199  G4int idxState,
200  G4int idxType )
201 {
202  G4int encoding = encodingOffset[idxState];
203  encoding += iSpin[idxState] +1;
204  G4int iQ = 0;
205  G4int iQbar = 1;
206 
207  if ( idxType == TPi ) {
208  if (iIsoSpin3<0) {
209  iQ = 1;
210  iQbar = 0;
211  }
212  } else if ( idxType == TK ) {
213  iQ = 1;
214  iQbar = 0;
215  }
216 
217 
218  encoding += 100*GetQuarkContents(iQ, iIsoSpin3, idxType);
219  encoding += 10*GetQuarkContents(iQbar, iIsoSpin3, idxType);
220  if ( idxType == TPi ) {
221  if (iIsoSpin3<0) {
222  encoding *= -1;
223  }
224  } else if ( idxType == TAntiK ) {
225  encoding *= -1;
226  }
227 
228 // PDG2005
229 //
230  if (idxState == 9 ) {
231  if (idxType == TEta) {
232 // f2(1810) 9030225
233  encoding = 9030225;
234  } else if (idxType == TEtaPrime) {
235 // f2(2010) 9060225
236  encoding = 9060225;
237  }
238  }
239 
240 // PDG2013
241  if (idxState == 1 ) {
242  if (idxType == TEta) {
243 // f0(1370) 30221
244  encoding = 30221;
245  }
246  }
247  return encoding;
248 }
249 
250 G4DecayTable* G4ExcitedMesonConstructor::CreateDecayTable(
251  const G4String& parentName,
252  G4int iIso3,
253  G4int iState,
254  G4int iType)
255 {
256  // create decay table
257  G4DecayTable* decayTable = new G4DecayTable();
258  G4double br;
259 
260  if ((iType==TK)||(iType==TAntiK)) {
261 
262  if ( (br=bRatio[iState][iType][MKPi]) >0.0) {
263  AddKPiMode( decayTable, parentName, br, iIso3, iType );
264  }
265  if ( (br=bRatio[iState][iType][MKStarPi]) >0.0) {
266  AddKStarPiMode( decayTable, parentName, br, iIso3, iType );
267  }
268  if ( (br=bRatio[iState][iType][MKRho]) >0.0) {
269  AddKRhoMode( decayTable, parentName, br, iIso3, iType );
270  }
271  if ( (br=bRatio[iState][iType][MKOmega]) >0.0) {
272  AddKOmegaMode( decayTable, parentName, br, iIso3, iType );
273  }
274  if ( (br=bRatio[iState][iType][MKStar2Pi]) >0.0) {
275  AddKStar2PiMode( decayTable, parentName, br, iIso3, iType );
276  }
277  if ( (br=bRatio[iState][iType][MKTwoPi]) >0.0) {
278  AddKTwoPiMode( decayTable, parentName, br, iIso3, iType );
279  }
280  if ( (br=bRatio[iState][iType][MKEta]) >0.0) {
281  AddKEtaMode( decayTable, parentName, br, iIso3, iType );
282  }
283 
284  } else {
285  if ( (br=bRatio[iState][iType][MPiGamma]) >0.0) {
286  AddPiGammaMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
287  }
288  if ( (br=bRatio[iState][iType][MRhoGamma]) >0.0) {
289  AddRhoGammaMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
290  }
291  if ( (br=bRatio[iState][iType][M2Pi]) >0.0) {
292  Add2PiMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
293  }
294  if ( (br=bRatio[iState][iType][MPiRho]) >0.0) {
295  AddPiRhoMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
296  }
297  if ( (br=bRatio[iState][iType][MPiEta]) >0.0) {
298  AddPiEtaMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
299  }
300  if ( (br=bRatio[iState][iType][M3Pi]) >0.0) {
301  Add3PiMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
302  }
303  if ( (br=bRatio[iState][iType][M4Pi]) >0.0) {
304  Add4PiMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
305  }
306  if ( (br=bRatio[iState][iType][MKKStar]) >0.0) {
307  AddKKStarMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
308  }
309  if ( (br=bRatio[iState][iType][M2PiEta]) >0.0) {
310  Add2PiEtaMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
311  }
312  if ( (br=bRatio[iState][iType][MRhoEta]) >0.0) {
313  AddRhoEtaMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
314  }
315  if ( (br=bRatio[iState][iType][M2PiRho]) >0.0) {
316  Add2PiRhoMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
317  }
318  if ( (br=bRatio[iState][iType][M2PiOmega]) >0.0) {
319  Add2PiOmegaMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
320  }
321  if ( (br=bRatio[iState][iType][M2Eta]) >0.0) {
322  Add2EtaMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
323  }
324  if ( (br=bRatio[iState][iType][M2K]) >0.0) {
325  Add2KMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
326  }
327  if ( (br=bRatio[iState][iType][M2KPi]) >0.0) {
328  Add2KPiMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
329  }
330  if ( (br=bRatio[iState][iType][MPiOmega]) >0.0) {
331  AddPiOmegaMode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
332  }
333  if ( (br=bRatio[iState][iType][MPiF2]) >0.0) {
334  AddPiF2Mode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
335  }
336  if ( (br=bRatio[iState][iType][MPiF0]) >0.0) {
337  AddPiF0Mode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
338  }
339  if ( (br=bRatio[iState][iType][MPiA2]) >0.0) {
340  AddPiA2Mode( decayTable, parentName, br, iIso3, iIsoSpin[iType] );
341  }
342  }
343 
344  return decayTable;
345 }
346 
347 G4DecayTable* G4ExcitedMesonConstructor::AddKPiMode(
348  G4DecayTable* decayTable, const G4String& nameParent,
349  G4double br, G4int iIso3, G4int iType)
350 {
351  G4VDecayChannel* mode;
352  //
353  if (iIso3 == +1) {
354  if (iType == TK) {
355  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
356  "kaon+","pi0");
357  decayTable->Insert(mode);
358  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
359  "kaon0","pi+");
360  decayTable->Insert(mode);
361  }else if (iType==TAntiK) {
362  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
363  "anti_kaon0","pi0");
364  decayTable->Insert(mode);
365  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
366  "kaon-","pi+");
367  decayTable->Insert(mode);
368  }
369  } else if (iIso3 == -1) {
370  if (iType == TK) {
371  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
372  "kaon0","pi0");
373  decayTable->Insert(mode);
374  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
375  "kaon+","pi-");
376  decayTable->Insert(mode);
377 
378  }else if (iType==TAntiK) {
379  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
380  "kaon-","pi0");
381  decayTable->Insert(mode);
382  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
383  "anti_kaon0","pi-");
384  decayTable->Insert(mode);
385  }
386  }
387 
388  return decayTable;
389 }
390 G4DecayTable* G4ExcitedMesonConstructor::AddKTwoPiMode(
391  G4DecayTable* decayTable, const G4String& nameParent,
392  G4double br, G4int iIso3, G4int iType)
393 {
394  G4VDecayChannel* mode;
395  //
396  if (iIso3 == +1) {
397  if (iType == TK) {
398  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
399  "k2_star(1430)+","pi0");
400  decayTable->Insert(mode);
401  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
402  "k2_star(1430)0","pi+");
403  decayTable->Insert(mode);
404  }else if (iType==TAntiK) {
405  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
406  "anti_k2_star(1430)0","pi0");
407  decayTable->Insert(mode);
408  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
409  "k2_star(1430)-","pi+");
410  decayTable->Insert(mode);
411  }
412  } else if (iIso3 == -1) {
413  if (iType == TK) {
414  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
415  "k2_star(1430)0","pi0");
416  decayTable->Insert(mode);
417  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
418  "k2_star(1430)+","pi-");
419  decayTable->Insert(mode);
420 
421  }else if (iType==TAntiK) {
422  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
423  "k2_star(1430)-","pi0");
424  decayTable->Insert(mode);
425  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
426  "anti_k2_star(1430)0","pi-");
427  decayTable->Insert(mode);
428  }
429  }
430 
431  return decayTable;
432 }
433 
434 G4DecayTable* G4ExcitedMesonConstructor::AddKOmegaMode(
435  G4DecayTable* decayTable, const G4String& nameParent,
436  G4double br, G4int iIso3, G4int iType)
437 {
438  G4VDecayChannel* mode;
439  //
440  if (iIso3 == +1) {
441  if (iType == TK) {
442  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
443  "kaon+","omega");
444  decayTable->Insert(mode);
445  }else if (iType==TAntiK) {
446  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
447  "anti_kaon0","omega");
448  decayTable->Insert(mode);
449  }
450  } else if (iIso3 == -1) {
451  if (iType == TK) {
452  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
453  "kaon0","omega");
454  decayTable->Insert(mode);
455  }else if (iType==TAntiK) {
456  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
457  "kaon-","omega");
458  decayTable->Insert(mode);
459  }
460  }
461 
462  return decayTable;
463 }
464 
465 G4DecayTable* G4ExcitedMesonConstructor::AddKEtaMode(
466  G4DecayTable* decayTable, const G4String& nameParent,
467  G4double br, G4int iIso3, G4int iType)
468 {
469  G4VDecayChannel* mode;
470  //
471  if (iIso3 == +1) {
472  if (iType == TK) {
473  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
474  "kaon+","eta");
475  decayTable->Insert(mode);
476  }else if (iType==TAntiK) {
477  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
478  "anti_kaon0","eta");
479  decayTable->Insert(mode);
480  }
481  } else if (iIso3 == -1) {
482  if (iType == TK) {
483  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
484  "kaon0","eta");
485  decayTable->Insert(mode);
486  }else if (iType==TAntiK) {
487  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
488  "kaon-","eta");
489  decayTable->Insert(mode);
490  }
491  }
492 
493  return decayTable;
494 }
495 
496 G4DecayTable* G4ExcitedMesonConstructor::AddKRhoMode(
497  G4DecayTable* decayTable, const G4String& nameParent,
498  G4double br, G4int iIso3, G4int iType)
499 {
500  G4VDecayChannel* mode;
501  //
502  if (iIso3 == +1) {
503  if (iType == TK) {
504  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
505  "kaon+","rho0");
506  decayTable->Insert(mode);
507  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
508  "kaon0","rho+");
509  decayTable->Insert(mode);
510  }else if (iType==TAntiK) {
511  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
512  "anti_kaon0","rho0");
513  decayTable->Insert(mode);
514  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
515  "kaon-","rho+");
516  decayTable->Insert(mode);
517  }
518  } else if (iIso3 == -1) {
519  if (iType == TK) {
520  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
521  "kaon0","rho0");
522  decayTable->Insert(mode);
523  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
524  "kaon+","rho-");
525  decayTable->Insert(mode);
526 
527  }else if (iType==TAntiK) {
528  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
529  "kaon-","rho0");
530  decayTable->Insert(mode);
531  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
532  "anti_kaon0","rho-");
533  decayTable->Insert(mode);
534  }
535  }
536 
537  return decayTable;
538 }
539 
540 G4DecayTable* G4ExcitedMesonConstructor::AddKStarPiMode(
541  G4DecayTable* decayTable, const G4String& nameParent,
542  G4double br, G4int iIso3, G4int iType)
543 {
544  G4VDecayChannel* mode;
545  //
546  if (iIso3 == +1) {
547  if (iType == TK) {
548  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
549  "k_star+","pi0");
550  decayTable->Insert(mode);
551  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
552  "k_star0","pi+");
553  decayTable->Insert(mode);
554  }else if (iType==TAntiK) {
555  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
556  "anti_k_star0","pi0");
557  decayTable->Insert(mode);
558  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
559  "k_star-","pi+");
560  decayTable->Insert(mode);
561  }
562  } else if (iIso3 == -1) {
563  if (iType == TK) {
564  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
565  "k_star0","pi0");
566  decayTable->Insert(mode);
567  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
568  "k_star+","pi-");
569  decayTable->Insert(mode);
570 
571  }else if (iType==TAntiK) {
572  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 2,
573  "k_star-","pi0");
574  decayTable->Insert(mode);
575  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 2,
576  "anti_k_star0","pi-");
577  decayTable->Insert(mode);
578  }
579  }
580 
581  return decayTable;
582 }
583 
584 G4DecayTable* G4ExcitedMesonConstructor::AddKStar2PiMode(
585  G4DecayTable* decayTable, const G4String& nameParent,
586  G4double br, G4int iIso3, G4int iType)
587 {
588  // K* --> K pipi(I=1)
589  G4VDecayChannel* mode;
590  //
591  if (iIso3 == +1) {
592  if (iType == TK) {
593  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 3,
594  "k_star+","pi+","pi-");
595  decayTable->Insert(mode);
596  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 3,
597  "k_star0","pi+","pi0");
598  decayTable->Insert(mode);
599  }else if (iType==TAntiK) {
600  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 3,
601  "anti_k_star0","pi+","pi-");
602  decayTable->Insert(mode);
603  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 3,
604  "k_star-","pi+","pi0");
605  decayTable->Insert(mode);
606  }
607  } else if (iIso3 == -1) {
608  if (iType == TK) {
609  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 3,
610  "k_star0","pi+","pi-");
611  decayTable->Insert(mode);
612  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 3,
613  "k_star+","pi-","pi0");
614  decayTable->Insert(mode);
615 
616  }else if (iType==TAntiK) {
617  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 3,
618  "k_star-","pi+","pi-");
619  decayTable->Insert(mode);
620  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 3,
621  "anti_k_star0","pi-","pi0");
622  decayTable->Insert(mode);
623  }
624  }
625 
626  return decayTable;
627 }
628 
629 G4DecayTable* G4ExcitedMesonConstructor::AddPiGammaMode(
630  G4DecayTable* decayTable, const G4String& nameParent,
631  G4double br, G4int iIso3, G4int iIso)
632 {
633  if ((iIso!=2)&&(iIso!=0)) return decayTable;
634 
635  G4VDecayChannel* mode;
636  //
637  G4String daughter;
638  if (iIso3 == +2) {
639  daughter = "pi+";
640  } else if (iIso3 == 0) {
641  daughter = "pi0";
642  } else if (iIso3 ==-2) {
643  daughter = "pi-";
644  } else {
645  return decayTable;
646  }
647  // create decay channel [parent BR #daughters]
648  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
649  daughter,"gamma");
650  // add decay table
651  decayTable->Insert(mode);
652 
653  return decayTable;
654 }
655 
656 G4DecayTable* G4ExcitedMesonConstructor::AddPiOmegaMode(
657  G4DecayTable* decayTable, const G4String& nameParent,
658  G4double br, G4int iIso3, G4int iIso)
659 {
660  if ((iIso!=2)&&(iIso!=0)) return decayTable;
661 
662  G4VDecayChannel* mode;
663  //
664  G4String daughter;
665  if (iIso3 == +2) {
666  daughter = "pi+";
667  } else if (iIso3 == 0) {
668  daughter = "pi0";
669  } else if (iIso3 ==-2) {
670  daughter = "pi-";
671  } else {
672  return decayTable;
673  }
674  // create decay channel [parent BR #daughters]
675  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
676  daughter,"omega");
677  // add decay table
678  decayTable->Insert(mode);
679 
680  return decayTable;
681 }
682 
683 G4DecayTable* G4ExcitedMesonConstructor::AddRhoGammaMode(
684  G4DecayTable* decayTable, const G4String& nameParent,
685  G4double br, G4int iIso3, G4int iIso)
686 {
687  if ((iIso!=2)&&(iIso!=0)) return decayTable;
688 
689  G4VDecayChannel* mode;
690  //
691  G4String daughter;
692  if (iIso3 == +2) {
693  daughter = "rho+";
694  } else if (iIso3 == 0) {
695  daughter = "rho0";
696  } else if (iIso3 ==-2) {
697  daughter = "rho-";
698  } else {
699  return decayTable;
700  }
701  // create decay channel [parent BR #daughters]
702  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
703  daughter,"gamma");
704  // add decay table
705  decayTable->Insert(mode);
706 
707  return decayTable;
708 }
709 
710 G4DecayTable* G4ExcitedMesonConstructor::AddPiEtaMode(
711  G4DecayTable* decayTable, const G4String& nameParent,
712  G4double br, G4int iIso3, G4int iIso)
713 {
714  if ((iIso!=2)&&(iIso!=0)) return decayTable;
715 
716  G4VDecayChannel* mode;
717  //
718  G4String daughter;
719  if (iIso3 == +2) {
720  daughter = "pi+";
721  } else if (iIso3 == 0) {
722  daughter = "pi0";
723  } else if (iIso3 ==-2) {
724  daughter = "pi-";
725  } else {
726  return decayTable;
727  }
728  // create decay channel [parent BR #daughters]
729  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
730  daughter,"eta");
731  // add decay table
732  decayTable->Insert(mode);
733 
734  return decayTable;
735 }
736 
737 G4DecayTable* G4ExcitedMesonConstructor::AddRhoEtaMode(
738  G4DecayTable* decayTable, const G4String& nameParent,
739  G4double br, G4int iIso3, G4int iIso)
740 {
741  if ((iIso!=2)&&(iIso!=0)) return decayTable;
742 
743  G4VDecayChannel* mode;
744  //
745  G4String daughter;
746  if (iIso3 == +2) {
747  daughter = "rho+";
748  } else if (iIso3 == 0) {
749  daughter = "rho0";
750  } else if (iIso3 ==-2) {
751  daughter = "rho-";
752  } else {
753  return decayTable;
754  }
755  // create decay channel [parent BR #daughters]
756  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
757  daughter,"eta");
758  // add decay table
759  decayTable->Insert(mode);
760 
761  return decayTable;
762 }
763 
764 G4DecayTable* G4ExcitedMesonConstructor::AddPiF2Mode(
765  G4DecayTable* decayTable, const G4String& nameParent,
766  G4double br, G4int iIso3, G4int iIso)
767 {
768  if ((iIso!=2)&&(iIso!=0)) return decayTable;
769 
770  G4VDecayChannel* mode;
771  //
772  G4String daughter;
773  if (iIso3 == +2) {
774  daughter = "pi+";
775  } else if (iIso3 == 0) {
776  daughter = "pi0";
777  } else if (iIso3 ==-2) {
778  daughter = "pi-";
779  } else {
780  return decayTable;
781  }
782  // create decay channel [parent BR #daughters]
783  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
784  daughter,"f2(1270)");
785  // add decay table
786  decayTable->Insert(mode);
787 
788  return decayTable;
789 }
790 
791 G4DecayTable* G4ExcitedMesonConstructor::AddPiF0Mode(
792  G4DecayTable* decayTable, const G4String& nameParent,
793  G4double br, G4int iIso3, G4int iIso)
794 {
795  if ((iIso!=2)&&(iIso!=0)) return decayTable;
796 
797  G4VDecayChannel* mode;
798  //
799  G4String daughter;
800  if (iIso3 == +2) {
801  daughter = "pi+";
802  } else if (iIso3 == 0) {
803  daughter = "pi0";
804  } else if (iIso3 ==-2) {
805  daughter = "pi-";
806  } else {
807  return decayTable;
808  }
809  // create decay channel [parent BR #daughters]
810  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
811  daughter,"f0(1370)");
812  // add decay table
813  decayTable->Insert(mode);
814  return decayTable;
815 }
816 
817 
818 G4DecayTable* G4ExcitedMesonConstructor::Add2PiMode(
819  G4DecayTable* decayTable, const G4String& nameParent,
820  G4double br, G4int iIso3, G4int iIso)
821 {
822  G4VDecayChannel* mode;
823 
824  G4String daughterPi1;
825  G4String daughterPi2;
826  G4double r;
827 
828  // I = 0 states
829  if (iIso==0) {
830  if (iIso3==0) {
831  // pi+ + pi-
832  daughterPi1 = "pi+";
833  daughterPi2 = "pi-";
834  r = br*2./3.;
835  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
836  daughterPi1,daughterPi2);
837  decayTable->Insert(mode);
838 
839  // pi0 + pi0
840  daughterPi1 = "pi0";
841  daughterPi2 = "pi0";
842  r = br*1./3.;
843  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
844  daughterPi1,daughterPi2);
845  decayTable->Insert(mode);
846  }
847  } else if (iIso==2) {
848  if (iIso3==+2) {
849  // pi+ + pi0
850  daughterPi1 = "pi+";
851  daughterPi2 = "pi0";
852  r = br;
853  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
854  daughterPi1,daughterPi2);
855  // add decay table
856  decayTable->Insert(mode);
857  } else if (iIso3==0) {
858  // pi+ + pi-
859  daughterPi1 = "pi+";
860  daughterPi2 = "pi-";
861  r = br;
862  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
863  daughterPi1,daughterPi2);
864  decayTable->Insert(mode);
865  } else if (iIso3==-2) {
866  // pi- + pi0
867  daughterPi1 = "pi-";
868  daughterPi2 = "pi0";
869  r = br;
870  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
871  daughterPi1,daughterPi2);
872  decayTable->Insert(mode);
873  }
874  }
875  return decayTable;
876 }
877 
878 G4DecayTable* G4ExcitedMesonConstructor::AddPiRhoMode(
879  G4DecayTable* decayTable, const G4String& nameParent,
880  G4double br, G4int iIso3, G4int iIso)
881 {
882  G4VDecayChannel* mode;
883 
884  G4String daughterPi;
885  G4String daughterRho;
886  G4double r;
887 
888  // I = 0 states
889  if (iIso==0) {
890  if (iIso3==0) {
891  // pi+ + rho-
892  daughterPi = "pi+";
893  daughterRho = "rho-";
894  r = br/3.;
895  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
896  daughterPi,daughterRho);
897  decayTable->Insert(mode);
898 
899  // pi0 + rho0
900  daughterPi = "pi0";
901  daughterRho = "rho0";
902  r = br*1./3.;
903  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
904  daughterPi,daughterRho);
905  decayTable->Insert(mode);
906 
907  // pi- + rho+
908  daughterPi = "pi-";
909  daughterRho = "rho+";
910  r = br*1./3.;
911  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
912  daughterPi,daughterRho);
913  decayTable->Insert(mode);
914  }
915  } else if (iIso==2) {
916  if (iIso3==+2) {
917  // pi+ + rho0
918  daughterPi = "pi+";
919  daughterRho = "rho0";
920  r = br/2.;
921  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
922  daughterPi,daughterRho);
923  decayTable->Insert(mode);
924 
925  // pi0 + rho+
926  daughterPi = "pi0";
927  daughterRho = "rho+";
928  r = br/2.;
929  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
930  daughterPi,daughterRho);
931  decayTable->Insert(mode);
932  } else if (iIso3==0) {
933  // pi+ + rho-
934  daughterPi = "pi+";
935  daughterRho = "rho-";
936  r = br/2.;
937  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
938  daughterPi,daughterRho);
939  decayTable->Insert(mode);
940 
941  // pi- + rho+
942  daughterPi = "pi-";
943  daughterRho = "rho+";
944  r = br/2.;
945  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
946  daughterPi,daughterRho);
947  decayTable->Insert(mode);
948  } else if (iIso3==-2) {
949  // pi- + rho0
950  daughterPi = "pi-";
951  daughterRho = "rho0";
952  r = br/2.;
953  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
954  daughterPi,daughterRho);
955  decayTable->Insert(mode);
956 
957  // pi0 + rho-
958  daughterPi = "pi0";
959  daughterRho = "rho-";
960  r = br/2.;
961  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
962  daughterPi,daughterRho);
963  decayTable->Insert(mode);
964  }
965  }
966  return decayTable;
967 }
968 
969 G4DecayTable* G4ExcitedMesonConstructor::AddPiA2Mode(
970  G4DecayTable* decayTable, const G4String& nameParent,
971  G4double br, G4int iIso3, G4int iIso)
972 {
973  G4VDecayChannel* mode;
974 
975  G4String daughterPi;
976  G4String daughterA2;
977  G4double r;
978 
979  // I = 0 states
980  if (iIso==0) {
981  if (iIso3==0) {
982  // pi+ + a2(1320)-
983  daughterPi = "pi+";
984  daughterA2 = "a2(1320)-";
985  r = br/3.;
986  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
987  daughterPi,daughterA2);
988  decayTable->Insert(mode);
989 
990  // pi0 + a2(1320)0
991  daughterPi = "pi0";
992  daughterA2 = "a2(1320)0";
993  r = br*1./3.;
994  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
995  daughterPi,daughterA2);
996  decayTable->Insert(mode);
997 
998  // pi- + a2(1320)+
999  daughterPi = "pi-";
1000  daughterA2 = "a2(1320)+";
1001  r = br*1./3.;
1002  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
1003  daughterPi,daughterA2);
1004  decayTable->Insert(mode);
1005  }
1006  } else if (iIso==2) {
1007  if (iIso3==+2) {
1008  // pi+ + a2(1320)0
1009  daughterPi = "pi+";
1010  daughterA2 = "a2(1320)0";
1011  r = br/2.;
1012  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
1013  daughterPi,daughterA2);
1014  decayTable->Insert(mode);
1015 
1016  // pi0 + a2(1320)+
1017  daughterPi = "pi0";
1018  daughterA2 = "a2(1320)+";
1019  r = br/2.;
1020  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
1021  daughterPi,daughterA2);
1022  decayTable->Insert(mode);
1023  } else if (iIso3==0) {
1024  // pi+ + a2(1320)-
1025  daughterPi = "pi+";
1026  daughterA2 = "a2(1320)-";
1027  r = br/2.;
1028  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
1029  daughterPi,daughterA2);
1030  decayTable->Insert(mode);
1031 
1032  // pi- + a2(1320)+
1033  daughterPi = "pi-";
1034  daughterA2 = "a2(1320)+";
1035  r = br/2.;
1036  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
1037  daughterPi,daughterA2);
1038  decayTable->Insert(mode);
1039  } else if (iIso3==-2) {
1040  // pi- + a2(1320)0
1041  daughterPi = "pi-";
1042  daughterA2 = "a2(1320)0";
1043  r = br/2.;
1044  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
1045  daughterPi,daughterA2);
1046  decayTable->Insert(mode);
1047 
1048  // pi0 + a2(1320)-
1049  daughterPi = "pi0";
1050  daughterA2 = "a2(1320)-";
1051  r = br/2.;
1052  mode = new G4PhaseSpaceDecayChannel(nameParent, r, 2,
1053  daughterPi,daughterA2);
1054  decayTable->Insert(mode);
1055  }
1056  }
1057  return decayTable;
1058 }
1059 
1060 G4DecayTable* G4ExcitedMesonConstructor::Add3PiMode(
1061  G4DecayTable* decayTable, const G4String& nameParent,
1062  G4double br, G4int iIso3, G4int iIso)
1063 {
1064  G4VDecayChannel* mode;
1065 
1066  // I =0 state
1067  // This mode is X(I=0,J=1) --> pi+,pi-,pi0 mode
1068  if (iIso==0) {
1069  // pi+ + pi-
1070  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 3,
1071  "pi+","pi-","pi0");
1072  decayTable->Insert(mode);
1073  } else if (iIso==2) {
1074  // This mode is X(I=1) --> pi + pipi(I=0) mode
1075  if (iIso3==+2) {
1076  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 3,
1077  "pi+","pi0","pi0");
1078  decayTable->Insert(mode);
1079  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 3,
1080  "pi+","pi+","pi-");
1081  decayTable->Insert(mode);
1082  } else if (iIso3==0) {
1083  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 3,
1084  "pi0","pi0","pi0");
1085  decayTable->Insert(mode);
1086  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 3,
1087  "pi0","pi+","pi-");
1088  decayTable->Insert(mode);
1089  } else if (iIso3==-2) {
1090  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 3,
1091  "pi-","pi0","pi0");
1092  decayTable->Insert(mode);
1093  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 3,
1094  "pi-","pi+","pi-");
1095  decayTable->Insert(mode);
1096  }
1097  }
1098  return decayTable;
1099 }
1100 
1101 G4DecayTable* G4ExcitedMesonConstructor::Add4PiMode(
1102  G4DecayTable* decayTable, const G4String& nameParent,
1103  G4double br, G4int iIso3, G4int )
1104 {
1105  G4VDecayChannel* mode;
1106 
1107  if (iIso3==0) {
1108  // 2pi+ + 2pi-
1109  mode = new G4PhaseSpaceDecayChannel(nameParent, br/2., 4,
1110  "pi+","pi-","pi+","pi-");
1111  decayTable->Insert(mode);
1112  // pi+ + pi- + 2pi0
1113  mode = new G4PhaseSpaceDecayChannel(nameParent, br/2., 4,
1114  "pi+","pi-","pi0","pi0");
1115  decayTable->Insert(mode);
1116  } else if (iIso3==+2) {
1117  // pi+ + 3pi0
1118  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 4,
1119  "pi+","pi0","pi0","pi0");
1120  decayTable->Insert(mode);
1121  // 2pi+ + pi- + pi0
1122  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 4,
1123  "pi+","pi+","pi-","pi0");
1124  decayTable->Insert(mode);
1125  } else if (iIso3==-2) {
1126  // pi- + 3pi0
1127  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 4,
1128  "pi-","pi0","pi0","pi0");
1129  decayTable->Insert(mode);
1130  // 2pi- + pi+ + pi0
1131  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 4,
1132  "pi-","pi-","pi+","pi0");
1133  decayTable->Insert(mode);
1134  }
1135  return decayTable;
1136 }
1137 
1138 G4DecayTable* G4ExcitedMesonConstructor::Add2PiEtaMode(
1139  G4DecayTable* decayTable, const G4String& nameParent,
1140  G4double br, G4int , G4int iIso)
1141 {
1142  // f1-->eta + pi + pi mode
1143 
1144  if (iIso!=0) return decayTable;
1145 
1146  G4VDecayChannel* mode;
1147 
1148  // eta pi+ pi-
1149  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 3,
1150  "eta","pi+","pi-");
1151  decayTable->Insert(mode);
1152 
1153  // eta pi+ pi-
1154  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 3,
1155  "eta","pi0","pi0");
1156  decayTable->Insert(mode);
1157  return decayTable;
1158 }
1159 
1160 G4DecayTable* G4ExcitedMesonConstructor::Add2EtaMode(
1161  G4DecayTable* decayTable, const G4String& nameParent,
1162  G4double br, G4int , G4int iIso)
1163 {
1164  if (iIso!=0) return decayTable;
1165 
1166  G4VDecayChannel* mode;
1167 
1168  // eta eta
1169  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
1170  "eta","eta");
1171  decayTable->Insert(mode);
1172  return decayTable;
1173 }
1174 
1175 G4DecayTable* G4ExcitedMesonConstructor::Add2PiOmegaMode(
1176  G4DecayTable* decayTable, const G4String& nameParent,
1177  G4double br, G4int iIso3, G4int iIso)
1178 {
1179 
1180  G4VDecayChannel* mode;
1181  if (iIso==0) {
1182  // omega pi+ pi-
1183  mode = new G4PhaseSpaceDecayChannel(nameParent, br*2./3., 3,
1184  "omega","pi+","pi-");
1185  decayTable->Insert(mode);
1186 
1187  // omega pi+ pi-
1188  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 3,
1189  "omega","pi0","pi0");
1190  decayTable->Insert(mode);
1191  } else if (iIso==2) {
1192  if (iIso3==+2) {
1193  // omega pi+ pi0
1194  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 3,
1195  "omega","pi+","pi0");
1196  decayTable->Insert(mode);
1197  } else if (iIso3==0) {
1198  // omega pi+ pi-
1199  mode = new G4PhaseSpaceDecayChannel(nameParent, br/2., 3,
1200  "omega","pi-","pi+");
1201  decayTable->Insert(mode);
1202  // omega pi0 pi0
1203  mode = new G4PhaseSpaceDecayChannel(nameParent, br/2., 3,
1204  "omega","pi0","pi0");
1205  decayTable->Insert(mode);
1206  } else if (iIso3==-2) {
1207  // omega pi- pi0
1208  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 3,
1209  "omega","pi-","pi0");
1210  decayTable->Insert(mode);
1211  }
1212  }
1213  return decayTable;
1214 }
1215 
1216 
1217 
1218 G4DecayTable* G4ExcitedMesonConstructor::Add2PiRhoMode(
1219  G4DecayTable* decayTable, const G4String& nameParent,
1220  G4double br, G4int iIso3, G4int iIso)
1221 {
1222  G4VDecayChannel* mode;
1223 
1224  if (iIso==0) {
1225  // f1 --> rho0 + pi+ pi-
1226  // rho0 pi+ pi-
1227  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1228  "rho0","pi+","pi-");
1229  decayTable->Insert(mode);
1230  } else if (iIso==2) {
1231  if (iIso3==+2) {
1232  // rho+ pi0 pi0
1233  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1234  "rho+","pi0","pi0");
1235  decayTable->Insert(mode);
1236  // rho+ pi+ pi-
1237  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1238  "rho+","pi+","pi-");
1239  decayTable->Insert(mode);
1240  // rho0 pi+ pi0
1241  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1242  "rho0","pi+","pi0");
1243  decayTable->Insert(mode);
1244  // rho- pi+ pi+
1245  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1246  "rho-","pi+","pi+");
1247  decayTable->Insert(mode);
1248  } else if (iIso3==-2) {
1249  // rho- pi0 pi0
1250  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1251  "rho-","pi0","pi0");
1252  decayTable->Insert(mode);
1253  // rho- pi+ pi-
1254  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1255  "rho-","pi+","pi-");
1256  decayTable->Insert(mode);
1257  // rho0 pi- pi0
1258  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1259  "rho0","pi-","pi0");
1260  decayTable->Insert(mode);
1261  // rho+ pi- pi-
1262  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1263  "rho+","pi-","pi-");
1264  decayTable->Insert(mode);
1265  } else if (iIso3==0) {
1266  // rho+ pi- pi0
1267  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1268  "rho+","pi-","pi0");
1269  decayTable->Insert(mode);
1270  // rho0 pi+ pi-
1271  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1272  "rho0","pi+","pi-");
1273  decayTable->Insert(mode);
1274  // rho0 pi0 pi0
1275  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1276  "rho0","pi0","pi0");
1277  decayTable->Insert(mode);
1278  // rho- pi+ pi0
1279  mode = new G4PhaseSpaceDecayChannel(nameParent, br/5., 3,
1280  "rho-","pi+","pi-");
1281  decayTable->Insert(mode);
1282  }
1283  }
1284  return decayTable;
1285 }
1286 
1287 
1288 G4DecayTable* G4ExcitedMesonConstructor::AddKKStarMode(
1289  G4DecayTable* decayTable, const G4String& nameParent,
1290  G4double br, G4int iIso3, G4int )
1291 {
1292  G4VDecayChannel* mode;
1293 
1294  if (iIso3==0) {
1295  // X(I=0,J=1)-->K + Anti-K*, Anti_K + K* mode
1296  // K+ + K*-
1297  mode = new G4PhaseSpaceDecayChannel(nameParent, br/4., 2,
1298  "kaon+","k_star-");
1299  decayTable->Insert(mode);
1300 
1301  // K- + K*+
1302  mode = new G4PhaseSpaceDecayChannel(nameParent, br/4., 2,
1303  "kaon-","k_star0");
1304  decayTable->Insert(mode);
1305 
1306  // K0 + Anti_K*0
1307  mode = new G4PhaseSpaceDecayChannel(nameParent, br/4., 2,
1308  "kaon0","anti_k_star0");
1309  decayTable->Insert(mode);
1310 
1311  // Anti_K0 + K*0
1312  mode = new G4PhaseSpaceDecayChannel(nameParent, br/4., 2,
1313  "anti_kaon0","k_star0");
1314  decayTable->Insert(mode);
1315 
1316  } else if (iIso3==2) {
1317  // K+ + Anti_K*0
1318  mode = new G4PhaseSpaceDecayChannel(nameParent, br/2., 2,
1319  "kaon+","anti_k_star0");
1320  decayTable->Insert(mode);
1321 
1322  // K0 + K*+
1323  mode = new G4PhaseSpaceDecayChannel(nameParent, br/2., 2,
1324  "anti_kaon0","k_star+");
1325  decayTable->Insert(mode);
1326 
1327  } else if (iIso3==-2) {
1328  // K- + K*0
1329  mode = new G4PhaseSpaceDecayChannel(nameParent, br/2., 2,
1330  "kaon-","k_star0");
1331  decayTable->Insert(mode);
1332 
1333  // K0 + K*-
1334  mode = new G4PhaseSpaceDecayChannel(nameParent, br/2., 2,
1335  "kaon0","k_star-");
1336  decayTable->Insert(mode);
1337 
1338  }
1339 
1340  return decayTable;
1341 }
1342 
1343 G4DecayTable* G4ExcitedMesonConstructor::Add2KMode(
1344  G4DecayTable* decayTable, const G4String& nameParent,
1345  G4double br, G4int iIso3, G4int )
1346 {
1347  G4VDecayChannel* mode;
1348 
1349  if (iIso3==0) {
1350  // K+ + K-
1351  mode = new G4PhaseSpaceDecayChannel(nameParent, br/2., 2,
1352  "kaon+","kaon-");
1353  decayTable->Insert(mode);
1354 
1355  // K0 + Anti_K0
1356  mode = new G4PhaseSpaceDecayChannel(nameParent, br/2., 2,
1357  "kaon0","anti_kaon0");
1358  decayTable->Insert(mode);
1359  } else if (iIso3==+2) {
1360  // K+ + anti_K0
1361  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
1362  "kaon+","anti_kaon0");
1363  decayTable->Insert(mode);
1364  } else if (iIso3==-2) {
1365  // K- + K0
1366  mode = new G4PhaseSpaceDecayChannel(nameParent, br, 2,
1367  "kaon-","kaon0");
1368  decayTable->Insert(mode);
1369  }
1370 
1371  return decayTable;
1372 }
1373 
1374 G4DecayTable* G4ExcitedMesonConstructor::Add2KPiMode(
1375  G4DecayTable* decayTable, const G4String& nameParent,
1376  G4double br, G4int , G4int iIso)
1377 {
1378 
1379  // X(I=0)-->KKpi
1380  if (iIso!=0) return decayTable;
1381 
1382  G4VDecayChannel* mode;
1383 
1384  // K+ + K- + pi0
1385  mode = new G4PhaseSpaceDecayChannel(nameParent, br/6., 3,
1386  "kaon+","kaon-","pi0");
1387  decayTable->Insert(mode);
1388 
1389  // K0 + Anti_K0 + pi0
1390  mode = new G4PhaseSpaceDecayChannel(nameParent, br/6., 3,
1391  "kaon0","anti_kaon0","pi0");
1392  decayTable->Insert(mode);
1393 
1394  // K+ + anti_K0 + pi-
1395  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 3,
1396  "kaon+","anti_kaon0","pi-");
1397  decayTable->Insert(mode);
1398 
1399  // K- + K0 + pi+
1400  mode = new G4PhaseSpaceDecayChannel(nameParent, br/3., 3,
1401  "kaon-","kaon0","pi+");
1402  decayTable->Insert(mode);
1403 
1404 
1405  return decayTable;
1406 }
1407 
1408 // PDG2005
1409 // eta(1440) is renamed to eta(1475)
1410 // omega(1600) is renamed to omega(1650)
1411 //
1412 //
1413 
1414 
1415 const char* G4ExcitedMesonConstructor::name[G4ExcitedMesonConstructor::NMultiplets ][ G4ExcitedMesonConstructor::NMesonTypes ] =
1416 {
1417  { "b1(1235)", "h1(1170)", "h1(1380)", "k1(1270)", "k1(1270)" },
1418  { "a0(1450)", "f0(1370)", "", "k0_star(1430)", "k0_star(1430)" },
1419  { "a1(1260)", "f1(1285)", "f1(1420)", "k1(1400)", "k1(1400)" },
1420  { "a2(1320)", "f2(1270)","f2_prime(1525)","k2_star(1430)","k2_star(1430)"},
1421  {"pi2(1670)", "eta2(1645)", "eta2(1870)", "k2(1770)", "k2(1770)" },
1422  {"rho(1700)", "omega(1650)", "", "k_star(1680)", "k_star(1680)" },
1423  {"rho3(1690)","omega3(1670)","phi3(1850)", "k3_star(1780)", "k3_star(1780)" },
1424  { "pi(1300)", "eta(1295)", "eta(1475)", "k(1460)", "k(1460)" },
1425  {"rho(1450)","omega(1420)", "phi(1680)", "k_star(1410)", "k_star(1410)" },
1426  { "", "f2(1810)", "f2(2010)", "k2_star(1980)", "k2_star(1980)" }
1427 };
1428 
1429 const G4double G4ExcitedMesonConstructor::mass[G4ExcitedMesonConstructor::NMultiplets ][ G4ExcitedMesonConstructor::NMesonTypes ] =
1430 {
1431  { 1.2295*GeV, 1.170*GeV, 1.386*GeV, 1.272*GeV, 1.272*GeV },
1432  { 1.474*GeV, 1.350*GeV, 0.0, 1.430*GeV, 1.430*GeV },
1433  { 1.230*GeV,1.2819*GeV,1.4264*GeV, 1.403*GeV, 1.403*GeV },
1434  { 1.3183*GeV,1.2755*GeV, 1.525*GeV,1.4256*GeV, 1.4256*GeV },
1435  { 1.6722*GeV, 1.617*GeV, 1.842*GeV, 1.773*GeV, 1.773*GeV },
1436  { 1.720*GeV, 1.670*GeV, 0.0, 1.718*GeV, 1.718*GeV },
1437  { 1.6888*GeV, 1.667*GeV, 1.854*GeV, 1.776*GeV, 1.776*GeV },
1438  { 1.300*GeV, 1.294*GeV, 1.476*GeV, 1.460*GeV, 1.460*GeV },
1439  { 1.465*GeV, 1.425*GeV, 1.680*GeV, 1.421*GeV, 1.421*GeV },
1440  { 0.0, 1.815*GeV, 2.010*GeV, 1.973*GeV, 1.973*GeV }
1441 };
1442 
1443 const G4double G4ExcitedMesonConstructor::massKdiff[ G4ExcitedMesonConstructor::NMultiplets ] = {
1444  0.0*MeV, 0.0*MeV, 0.0*MeV, 6.8*MeV, 0.0*MeV,
1445  0.0*MeV, 0.0*MeV, 0.0*MeV, 0.0*MeV, 0.0*MeV
1446 };
1447 
1448 const G4double G4ExcitedMesonConstructor::widthKdiff[ G4ExcitedMesonConstructor::NMultiplets ] = {
1449  0.0*MeV, 0.0*MeV, 0.0*MeV, 10.5*MeV, 0.0*MeV,
1450  0.0*MeV, 0.0*MeV, 0.0*MeV, 0.0*MeV, 0.0*MeV
1451 };
1452 
1453 const G4double G4ExcitedMesonConstructor::width[G4ExcitedMesonConstructor::NMultiplets ][ G4ExcitedMesonConstructor::NMesonTypes ] =
1454 {
1455  { 142.0*MeV, 360.0*MeV, 91.0*MeV, 90.0*MeV, 90.0*MeV },
1456  { 265.0*MeV, 350.0*MeV, 0.0, 270.0*MeV, 270.0*MeV },
1457  { 420.0*MeV, 22.7*MeV, 54.9*MeV, 174.0*MeV, 174.0*MeV },
1458  { 107.0*MeV, 186.7*MeV, 73.0*MeV, 98.5*MeV, 98.5*MeV },
1459  { 260.0*MeV, 181.0*MeV, 225.0*MeV, 186.0*MeV, 186.0*MeV },
1460  { 250.0*MeV, 315.0*MeV, 0.0, 320.0*MeV, 320.0*MeV },
1461  { 161.0*MeV, 168.0*MeV, 87.0*MeV, 159.0*MeV, 159.0*MeV },
1462  { 400.0*MeV, 55.0*MeV, 85.0*MeV, 260.0*MeV, 260.0*MeV },
1463  { 400.0*MeV, 215.0*MeV, 150.0*MeV, 236.0*MeV, 236.0*MeV },
1464  { 0.0, 197.0*MeV, 200.0*MeV, 373.0*MeV, 373.0*MeV }
1465 };
1466 
1467 
1468 const G4int G4ExcitedMesonConstructor::iIsoSpin[] =
1469 {
1470 // Tpi TEta TEtaPrime TK TAntiK
1471  2, 0, 0, 1, 1
1472 };
1473 
1474 const G4int G4ExcitedMesonConstructor::iSpin[] =
1475 {
1476 //N 1 1 1 1 1 1 1 2 2 2
1477 // 1P1 3P0 3P1 3P2 1D2 3D1 3D3 1S0 3S1 3P2
1478  2, 0, 2, 4, 4, 2, 6, 0, 2, 4
1479 };
1480 
1481 const G4int G4ExcitedMesonConstructor::iParity[] =
1482 {
1483 //N 1 1 1 1 1 1 1 2 2 2
1484 // 1P1 3P0 3P1 3P2 1D2 3D1 3D3 1S0 3S1 3P2
1485  +1, +1, +1, +1, -1, -1, -1, -1, -1, +1
1486 };
1487 
1488 const G4int G4ExcitedMesonConstructor::iChargeConjugation[] =
1489 {
1490 //N 1 1 1 1 1 1 1 2 2 2
1491 // 1P1 3P0 3P1 3P2 1D2 3D1 3D3 1S0 3S1 3P2
1492  -1, +1, +1, +1, +1, -1, -1, +1, -1, +1
1493 };
1494 
1495 const G4int G4ExcitedMesonConstructor::iGParity[G4ExcitedMesonConstructor::NMultiplets ][ G4ExcitedMesonConstructor::NMesonTypes ]=
1496 {
1497  { +1, -1, -1, 0, 0},
1498  { -1, +1, 0, 0, 0},
1499  { -1, +1, +1, 0, 0},
1500  { -1, +1, +1, 0, 0},
1501  { -1, +1, +1, 0, 0},
1502  { +1, -1, 0, 0, 0},
1503  { +1, -1, -1, 0, 0},
1504  { -1, +1, +1, 0, 0},
1505  { +1, -1, -1, 0, 0},
1506  { 0, +1, +1, 0, 0}
1507 };
1508 
1509 
1510 const G4int G4ExcitedMesonConstructor::encodingOffset[]=
1511 { 10000, 10000, 20000, 0, 10000, 30000, 0, 100000,100000,100000};
1512 
1513 
1514 
1515 
1516 const G4double G4ExcitedMesonConstructor::bRatio[G4ExcitedMesonConstructor::NMultiplets][G4ExcitedMesonConstructor::NMesonTypes][G4ExcitedMesonConstructor::NumberOfDecayModes]=
1517 {
1518  // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1519  // "b1(1235)", "h1(1170)", "h1(1380)", "k1(1270)", "k1(1270)"
1520  {
1521  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 1.00, 0.00, 0.00, 0.00 },
1522  { 0.00, 0.00, 0.00, 0.90, 0.10, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1523  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 1.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1524  { 0.00, 0.47, 0.42, 0.11, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1525  { 0.00, 0.47, 0.42, 0.11, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }
1526  },
1527  // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1528  // "a0(1450)", "f0(1370)", "", "k0_star(1430)", "k0_star(1430)"
1529  {
1530  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.90, 0.00, 0.00, 0.00, 0.10, 0.00, 0.00, 0.00, 0.10, 0.00, 0.00, 0.00, 0.00, 0.00 },
1531  { 0.00, 0.00, 0.10, 0.00, 0.00, 0.00, 0.70, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.20, 0.00, 0.00, 0.00, 0.00, 0.00 },
1532  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1533  { 1.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1534  { 1.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }
1535  },
1536  // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1537  // "a1(1260)", "f1(1285)", "f1(1420)", "k1(1400)", "k1(1400)"
1538  {
1539  { 0.10, 0.00, 0.00, 0.90, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1540  { 0.00, 0.07, 0.00, 0.00, 0.00, 0.00, 0.20, 0.00, 0.54, 0.00, 0.10, 0.00, 0.00, 0.00, 0.09, 0.00, 0.00, 0.00, 0.00 },
1541  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.50, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.50, 0.00, 0.00, 0.00, 0.00 },
1542  { 0.00, 0.96, 0.03, 0.01, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1543  { 0.00, 0.96, 0.03, 0.01, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }
1544  },
1545  // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1546  //"a2(1320)", "f2(1270)","f2_prime(1525)","k2_star(1430)","k2_star(1430)"
1547  {
1548  { 0.00, 0.00, 0.00, 0.70, 0.00, 0.14, 0.00, 0.00, 0.00, 0.00, 0.00, 0.11, 0.00, 0.05, 0.00, 0.00, 0.00, 0.00, 0.00 },
1549  { 0.00, 0.00, 0.50, 0.00, 0.00, 0.00, 0.30, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.20, 0.00, 0.00, 0.00, 0.00, 0.00 },
1550  { 0.00, 0.00, 0.01, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.10, 0.89, 0.00, 0.00, 0.00, 0.00, 0.00 },
1551  { 0.50, 0.25, 0.09, 0.03, 0.13, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1552  { 0.50, 0.25, 0.09, 0.03, 0.13, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }
1553  },
1554  // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1555  // "pi2(1670)", "eta2(1645)", "eta2(1870)", "k2(1770)", "k2(1770)"
1556  {
1557  { 0.00, 0.00, 0.00, 0.30, 0.00, 0.00, 0.00, 0.04, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.56, 0.10, 0.00 },
1558  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.10, 0.00, 0.00, 0.00, 0.90 },
1559  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.50, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.50 },
1560  { 0.00, 0.00, 0.00, 0.00, 0.00, 1.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1561  { 0.00, 0.00, 0.00, 0.00, 0.00, 1.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }
1562  },
1563  // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1564  // "rho(1700)", "omega(1650)", "", "k_star(1680)", "k_star(1680)"
1565  {
1566  { 0.00, 0.00, 0.10, 0.00, 0.00, 0.20, 0.00, 0.00, 0.00, 0.00, 0.70, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1567  { 0.00, 0.00, 0.00, 0.50, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.50, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1568  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1569  { 0.00, 0.40, 0.30, 0.30, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1570  { 0.00, 0.40, 0.30, 0.30, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }
1571  },
1572  // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1573  // "rho3(1690)","omega3(1670)","phi3(1850)", "k3_star(1780)", "k3_star(1780)"
1574  {
1575  { 0.00, 0.00, 0.24, 0.00, 0.00, 0.00, 0.60, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.01, 0.04, 0.11, 0.00, 0.00, 0.00 },
1576  { 0.00, 0.00, 0.00, 0.50, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.50, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1577  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.40, 0.00, 0.00, 0.00, 0.00, 0.00, 0.60, 0.00, 0.00, 0.00, 0.00, 0.00 },
1578  { 0.19, 0.20, 0.31, 0.00, 0.00, 0.00, 0.30, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1579  { 0.19, 0.20, 0.31, 0.00, 0.00, 0.00, 0.30, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }
1580  },
1581  // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1582  // "pi(1300)", "eta(1295)", "eta(1475)", "k(1460)", "k(1460)"
1583  {
1584  { 0.00, 0.00, 0.00, 0.50, 0.50, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1585  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 1.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1586  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.20, 0.20, 0.00, 0.00, 0.00, 0.00, 0.00, 0.60, 0.00, 0.00, 0.00, 0.00 },
1587  { 0.00, 0.50, 0.50, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1588  { 0.00, 0.50, 0.50, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }
1589  },
1590  // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1591  // "rho(1450)","omega(1420)", "phi(1680)", "k_star(1410)", "k_star(1410)"
1592  {
1593  { 0.00, 0.00, 0.50, 0.00, 0.00, 0.00, 0.50, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1594  { 0.00, 0.00, 0.00, 1.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1595  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.80, 0.00, 0.00, 0.00, 0.00, 0.00, 0.10, 0.10, 0.00, 0.00, 0.00, 0.00 },
1596  { 0.30, 0.65, 0.05, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1597  { 0.30, 0.65, 0.05, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }
1598  },
1599  // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1600  // "", "f2(1810)", "f2(2010)", "k2_star(1980)", "k2_star(1980)"
1601  {
1602  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1603  { 0.00, 0.00, 0.50, 0.00, 0.00, 0.00, 0.00, 0.00, 0.30, 0.00, 0.00, 0.00, 0.00, 0.20, 0.00, 0.00, 0.00, 0.00, 0.00 },
1604  { 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 1.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1605  { 0.00, 0.00, 0.60, 0.40, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 },
1606  { 0.00, 0.00, 0.60, 0.40, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }
1607  }
1608 };
1609 
1610 
1611 
1612 
1613 
1614