Changes needed for MC production of J/psi from B with radiative decays
[u/mrichter/AliRoot.git] / PYTHIA6 / AliDecayerPythia.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 /* $Id$ */
17
18 // Implementation of AliDecayer using Pythia
19 // Method forwarding to the AliPythia instance.
20 // Muonic and electronic decay modes can be forced
21 // for heavy flavor hadrons.
22 // Author: andreas.morsch@cern.ch
23
24
25 #include "AliDecayerPythia.h"
26 #include "AliPythia.h"
27 #include "AliLog.h"
28 #include <TPDGCode.h>
29 #include <TLorentzVector.h>
30 #include <TClonesArray.h>
31 #include <TParticle.h>
32
33 ClassImp(AliDecayerPythia)
34
35 #ifndef WIN32
36 # define py1ent py1ent_
37 # define opendecaytable opendecaytable_
38 # define closedecaytable closedecaytable_
39 # define type_of_call
40 #else
41 # define lu1ent PY1ENT
42 # define opendecaytable OPENDECAYTABLE
43 # define closedecaytable CLOSEDECAYTABLE
44 # define type_of_call _stdcall
45 #endif
46
47 extern "C" void type_of_call 
48           py1ent(Int_t&, Int_t&, Double_t&, Double_t&, Double_t&);
49
50 extern "C" void type_of_call 
51           opendecaytable(Int_t&);
52
53 extern "C" void type_of_call 
54           closedecaytable(Int_t&);
55
56 Bool_t AliDecayerPythia::fgInit = kFALSE;
57
58
59 AliDecayerPythia::AliDecayerPythia():
60     fPythia(AliPythia::Instance()),
61     fDecay(kAll),
62     fHeavyFlavour(kTRUE),
63     fLongLived(kFALSE),
64     fPatchOmegaDalitz(0)
65 {
66 // Default Constructor
67     for (Int_t i=0; i< 501; i++) fBraPart[i]= 1.;
68     ReadDecayTable();
69 }
70
71 AliDecayerPythia::AliDecayerPythia(const AliDecayerPythia &decayer):
72     AliDecayer(decayer),
73     fPythia(0),
74     fDecay(kAll),
75     fHeavyFlavour(kTRUE),
76     fLongLived(kFALSE),
77     fPatchOmegaDalitz(0)
78 {
79     // Copy Constructor
80     decayer.Copy(*this);
81     for (Int_t i = 0; i < 501; i++) fBraPart[i] = 0.;
82 }
83
84 void AliDecayerPythia::Init()
85 {
86     
87 // Initialisation
88 //
89     if (!fgInit) {
90         fgInit = kTRUE;
91         fPythia->SetDecayTable();
92     }
93
94 // Switch on heavy flavor decays
95     
96     Int_t kc, i, j;
97     Int_t heavy[14] = {411, 421, 431, 4122, 4132, 4232, 4332, 511, 521, 531, 5122, 5132, 5232, 5332};
98     fPythia->ResetDecayTable();
99     for (j=0; j < 14; j++) {
100         kc=fPythia->Pycomp(heavy[j]);
101         if (fDecay == kNoDecayHeavy) {
102           fPythia->SetMDCY(kc,1,0);
103         } else {
104           fPythia->SetMDCY(kc,1,1);
105           for (i=fPythia->GetMDCY(kc,2); 
106                i<fPythia->GetMDCY(kc,2)+fPythia->GetMDCY(kc,3); 
107                i++) {
108             fPythia->SetMDME(i,1,1);
109           }
110         }
111     }
112 //...Switch off decay of pi0, K0S, Lambda, Sigma+-, Xi0-, Omega-.
113
114 /*
115     if (fDecay != kNeutralPion) {
116         fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 0);
117     } else {
118         fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1);
119     }
120 */
121
122     fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1);
123     Int_t isw = 0;
124     if (fLongLived) isw = 1;
125     
126     fPythia->SetMDCY(fPythia->Pycomp(310) ,1, isw);
127     fPythia->SetMDCY(fPythia->Pycomp(3122),1, isw);
128     fPythia->SetMDCY(fPythia->Pycomp(3112),1, isw);
129 //    fPythia->SetMDCY(fPythia->Pycomp(3212),1, isw); // Sigma0 decays elem.
130     fPythia->SetMDCY(fPythia->Pycomp(3222),1, isw);
131     fPythia->SetMDCY(fPythia->Pycomp(3312),1, isw);
132     fPythia->SetMDCY(fPythia->Pycomp(3322),1, isw);
133     fPythia->SetMDCY(fPythia->Pycomp(3334),1, isw);
134     
135 // .. Force decay channels
136     ForceDecay();
137 }
138
139 void AliDecayerPythia::Decay(Int_t idpart, TLorentzVector* p)
140 {
141 //  Decay a particle
142 // 
143     Float_t energy = p->Energy();
144     Float_t theta  = p->Theta();
145     Float_t phi    = p->Phi();
146     if (!fPatchOmegaDalitz) {
147         Lu1Ent(0, idpart, energy, theta, phi);
148     } else {
149         fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 0);
150         Lu1Ent(0, idpart, energy, theta, phi);
151         fPythia->DalitzDecays();
152         fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1);
153         fPythia->Pyexec();
154     }
155     fPythia->GetPrimaries();
156 }
157
158 Int_t AliDecayerPythia::ImportParticles(TClonesArray *particles)
159 {
160 // Import the decay products
161 //
162   const Float_t kconvT=0.001/2.999792458e8; // mm/c to seconds conversion
163   const Float_t kconvL=1./10; // mm to cm conversion
164   int np = fPythia->ImportParticles(particles, "All");
165   // pythia assigns decay time in mm/c, convert to seconds
166   for (int ip=1;ip<np;ip++) {
167     TParticle* prod = (TParticle*)particles->At(ip);
168     if (!prod) continue;
169     prod->SetProductionVertex(prod->Vx()*kconvL,prod->Vy()*kconvL,prod->Vz()*kconvL,kconvT*prod->T());
170   }
171   return np;
172 }
173
174
175 void AliDecayerPythia::ForceDecay()
176 {
177 // Force a particle decay mode
178 // Switch heavy flavour production off if requested
179     if (!fHeavyFlavour) SwitchOffHeavyFlavour();
180 //
181     Decay_t decay=fDecay;
182     fPythia->SetMSTJ(21,2);
183
184 //
185 // select mode    
186     Int_t products[2];
187     Int_t mult[2];
188     Int_t products1[3];
189     Int_t mult1[3];
190     
191     switch (decay) 
192     {
193     case kHardMuons:
194         products1[0] =     13;
195         products1[1] =    443;
196         products1[2] = 100443;
197         mult1[0] = 1;
198         mult1[1] = 1;
199         mult1[2] = 1;
200         ForceParticleDecay(  511, products1, mult1, 3); 
201         ForceParticleDecay(  521, products1, mult1, 3); 
202         ForceParticleDecay(  531, products1, mult1, 3); 
203         ForceParticleDecay( 5122, products1, mult1, 3); 
204         ForceParticleDecay( 5132, products1, mult1, 3); 
205         ForceParticleDecay( 5232, products1, mult1, 3); 
206         ForceParticleDecay( 5332, products1, mult1, 3); 
207         ForceParticleDecay( 100443, 443, 1);  // Psi'  -> J/Psi X    
208         ForceParticleDecay(    443,  13, 2);  // J/Psi -> mu+ mu-
209         ForceParticleDecay(  411,13,1); // D+/-     
210         ForceParticleDecay(  421,13,1); // D0     
211         ForceParticleDecay(  431,13,1); // D_s     
212         ForceParticleDecay( 4122,13,1); // Lambda_c    
213         ForceParticleDecay( 4132,13,1); // Xsi_c     
214         ForceParticleDecay( 4232,13,1); // Sigma_c 
215         ForceParticleDecay( 4332,13,1); // Omega_c     
216         break;
217    case kChiToJpsiGammaToMuonMuon:
218         products[0] =    443;
219         products[1] =     22;
220         mult[0] = 1;
221         mult[1] = 1;
222         ForceParticleDecay( 20443, products, mult, 2);   // Chi_1c  -> J/Psi  Gamma
223         ForceParticleDecay(   445, products, mult, 2);   // Chi_2c  -> J/Psi  Gamma
224         ForceParticleDecay(   443, 13, 2);               // J/Psi -> mu+ mu-
225         break;
226     case kChiToJpsiGammaToElectronElectron:
227         products[0] =    443;
228         products[1] =     22;
229         mult[0] = 1;
230         mult[1] = 1;
231         ForceParticleDecay( 20443, products, mult, 2);   // Chi_1c  -> J/Psi  Gamma
232         ForceParticleDecay(   445, products, mult, 2);   // Chi_2c  -> J/Psi  Gamma
233         ForceParticleDecay(   443, 11, 2);               // J/Psi -> e+ e-
234         break;
235
236     case kBSemiMuonic:
237         ForceParticleDecay(  511,13,1); // B0     
238         ForceParticleDecay(  521,13,1); // B+/-     
239         ForceParticleDecay(  531,13,1); // B_s     
240         ForceParticleDecay( 5122,13,1); // Lambda_b    
241         ForceParticleDecay( 5132,13,1); // Xsi_b    
242         ForceParticleDecay( 5232,13,1); // Sigma_b    
243         ForceParticleDecay( 5332,13,1); // Omega_b    
244         break;
245     case kSemiMuonic:
246         ForceParticleDecay(  411,13,1); // D+/-     
247         ForceParticleDecay(  421,13,1); // D0     
248         ForceParticleDecay(  431,13,1); // D_s     
249         ForceParticleDecay( 4122,13,1); // Lambda_c    
250         ForceParticleDecay( 4132,13,1); // Xsi_c     
251         ForceParticleDecay( 4232,13,1); // Sigma_c 
252         ForceParticleDecay( 4332,13,1); // Omega_c     
253         ForceParticleDecay(  511,13,1); // B0     
254         ForceParticleDecay(  521,13,1); // B+/-     
255         ForceParticleDecay(  531,13,1); // B_s     
256         ForceParticleDecay( 5122,13,1); // Lambda_b    
257         ForceParticleDecay( 5132,13,1); // Xsi_b    
258         ForceParticleDecay( 5232,13,1); // Sigma_b    
259         ForceParticleDecay( 5332,13,1); // Omega_b    
260         break;
261     case kDiMuon:
262         ForceParticleDecay(  113,13,2); // rho
263         ForceParticleDecay(  221,13,2); // eta
264         ForceParticleDecay(  223,13,2); // omega
265         ForceParticleDecay(  333,13,2); // phi
266         ForceParticleDecay(  443,13,2); // J/Psi
267         ForceParticleDecay(100443,13,2);// Psi'
268         ForceParticleDecay(  553,13,2); // Upsilon
269         ForceParticleDecay(100553,13,2);// Upsilon'
270         ForceParticleDecay(200553,13,2);// Upsilon''
271         break;
272     case kJpsiDiMuon:
273         ForceParticleDecay(  443,13,2); // J/Psi
274         break;
275     case kBSemiElectronic:
276         ForceParticleDecay(  511,11,1); // B0     
277         ForceParticleDecay(  521,11,1); // B+/-     
278         ForceParticleDecay(  531,11,1); // B_s     
279         ForceParticleDecay( 5122,11,1); // Lambda_b     
280         ForceParticleDecay( 5132,11,1); // Xsi_b    
281         ForceParticleDecay( 5232,11,1); // Sigma_b    
282         ForceParticleDecay( 5332,11,1); // Omega_b    
283         break;
284     case kSemiElectronic:
285         ForceParticleDecay(  411,11,1); // D+/-     
286         ForceParticleDecay(  421,11,1); // D0     
287         ForceParticleDecay(  431,11,1); // D_s     
288         ForceParticleDecay( 4122,11,1); // Lambda_c     
289         ForceParticleDecay( 4132,11,1); // Xsi_c     
290         ForceParticleDecay( 4232,11,1); // Sigma_c 
291         ForceParticleDecay( 4332,11,1); // Omega_c     
292         ForceParticleDecay(  511,11,1); // B0     
293         ForceParticleDecay(  521,11,1); // B+/-     
294         ForceParticleDecay(  531,11,1); // B_s     
295         ForceParticleDecay( 5122,11,1); // Lambda_b     
296         ForceParticleDecay( 5132,11,1); // Xsi_b    
297         ForceParticleDecay( 5232,11,1); // Sigma_b    
298         ForceParticleDecay( 5332,11,1); // Omega_b    
299         break;
300     case kDiElectron:
301         ForceParticleDecay(  113,11,2); // rho
302         ForceParticleDecay(  333,11,2); // phi
303         ForceParticleDecay(  221,11,2); // eta
304         ForceParticleDecay(  223,11,2); // omega
305         ForceParticleDecay(  443,11,2); // J/Psi
306         ForceParticleDecay(100443,11,2);// Psi'
307         ForceParticleDecay(  553,11,2); // Upsilon
308         ForceParticleDecay(100553,11,2);// Upsilon'
309         ForceParticleDecay(200553,11,2);// Upsilon''
310         break;
311     case kPsiPrimeJpsiDiElectron:
312         products[0] =    443;
313         products[1] =    211;
314         mult[0] = 1;
315         mult[1] = 2;
316         ForceParticleDecay( 100443, products, mult, 2, 1);
317         ForceParticleDecay(    443,11,2);             
318         break;  
319     case kBJpsiDiMuon:
320
321         products[0] =    443;
322         products[1] = 100443;
323         mult[0] = 1;
324         mult[1] = 1;
325
326         ForceParticleDecay(  511, products, mult, 2); // B0   -> J/Psi (Psi') X   
327         ForceParticleDecay(  521, products, mult, 2); // B+/- -> J/Psi (Psi') X     
328         ForceParticleDecay(  531, products, mult, 2); // B_s  -> J/Psi (Psi') X     
329         ForceParticleDecay( 5122, products, mult, 2); // Lambda_b -> J/Psi (Psi') X 
330         ForceParticleDecay( 100443, 443, 1);          // Psi'  -> J/Psi X    
331         ForceParticleDecay(    443,13,2);             // J/Psi -> mu+ mu-   
332         break;
333     case kBPsiPrimeDiMuon:
334         ForceParticleDecay(  511,100443,1); // B0     
335         ForceParticleDecay(  521,100443,1); // B+/-     
336         ForceParticleDecay(  531,100443,1); // B_s     
337         ForceParticleDecay( 5122,100443,1); // Lambda_b 
338         ForceParticleDecay(100443,13,2);    // Psi'   
339         break;
340     case kBJpsiDiElectron:
341         ForceParticleDecay(  511,443,1); // B0     
342         ForceParticleDecay(  521,443,1); // B+/-     
343         ForceParticleDecay(  531,443,1); // B_s     
344         ForceParticleDecay( 5122,443,1); // Lambda_b
345         ForceParticleDecay(  443,11,2);  // J/Psi    
346         break;
347     case kBJpsi:
348         ForceParticleDecay(  511,443,1); // B0     
349         ForceParticleDecay(  521,443,1); // B+/-     
350         ForceParticleDecay(  531,443,1); // B_s     
351         ForceParticleDecay( 5122,443,1); // Lambda_b
352         break;
353     case kBJpsiUndecayed:
354         ForceParticleDecay(  511,443,1); // B0     
355         ForceParticleDecay(  521,443,1); // B+/-     
356         ForceParticleDecay(  531,443,1); // B_s     
357         ForceParticleDecay( 5122,443,1); // Lambda_b
358         fPythia->SetMDCY(fPythia->Pycomp(443),1,0); // switch-off J/psi 
359         break;
360     case kBPsiPrimeDiElectron:
361         ForceParticleDecay(  511,100443,1); // B0     
362         ForceParticleDecay(  521,100443,1); // B+/-     
363         ForceParticleDecay(  531,100443,1); // B_s     
364         ForceParticleDecay( 5122,100443,1); // Lambda_b 
365         ForceParticleDecay(100443,11,2);   // Psi'   
366         break;
367     case kPiToMu:
368         ForceParticleDecay(211,13,1); // pi->mu     
369         break;
370     case kKaToMu:
371         ForceParticleDecay(321,13,1); // K->mu     
372         break;
373     case kAllMuonic:
374         ForceParticleDecay(211,13,1); // pi->mu
375         ForceParticleDecay(321,13,1); // K->mu    
376         break;
377     case kWToMuon:
378         ForceParticleDecay(  24, 13,1); // W -> mu
379         break;
380     case kWToCharm:
381         ForceParticleDecay(   24, 4,1); // W -> c
382         break;
383     case kWToCharmToMuon:
384         ForceParticleDecay(   24, 4,1); // W -> c
385         ForceParticleDecay(  411,13,1); // D+/- -> mu
386         ForceParticleDecay(  421,13,1); // D0  -> mu
387         ForceParticleDecay(  431,13,1); // D_s  -> mu
388         ForceParticleDecay( 4122,13,1); // Lambda_c
389         ForceParticleDecay( 4132,13,1); // Xsi_c
390         ForceParticleDecay( 4232,13,1); // Sigma_c
391         ForceParticleDecay( 4332,13,1); // Omega_c
392         break;
393     case kZDiMuon:
394         ForceParticleDecay(  23, 13,2); // Z -> mu+ mu-
395         break;
396     case kZDiElectron:
397         ForceParticleDecay(  23, 11,2); // Z -> e+ e-
398         break;
399     case kHadronicD:
400         ForceHadronicD(1);
401         break;
402     case kHadronicDWithout4Bodies:
403         ForceHadronicD(0);
404         break;
405     case kPhiKK:
406         ForceParticleDecay(333,321,2); // Phi->K+K-
407         break;
408     case kOmega:
409         ForceOmega();
410         break;
411     case kLambda:
412         ForceLambda();
413         break;
414     case kAll:
415         break;
416     case kNoDecay:
417         fPythia->SetMSTJ(21,0);
418         break;
419     case kNoDecayHeavy:
420     case kNeutralPion:
421         break;
422     case kNoDecayBeauty:                                                    
423         SwitchOffBDecay(); 
424         break;
425     case kElectronEM:
426         ForceParticleDecay(  111,11,1); // pi^0     
427         ForceParticleDecay(  221,11,1); // eta     
428         ForceParticleDecay(  113,11,1); // rho     
429         ForceParticleDecay(  223,11,1); // omega     
430         ForceParticleDecay(  331,11,1); // etaprime     
431         ForceParticleDecay(  333,11,1); // phi     
432         break;
433     case kDiElectronEM:
434         ForceParticleDecay(  111,11,2); // pi^0     
435         ForceParticleDecay(  221,11,2); // eta     
436         ForceParticleDecay(  113,11,2); // rho     
437         ForceParticleDecay(  223,11,2); // omega     
438         ForceParticleDecay(  331,11,2); // etaprime     
439         ForceParticleDecay(  333,11,2); // phi     
440         break;
441     case kGammaEM:
442         ForceParticleDecay(  111,22,1); // pi^0     
443         ForceParticleDecay(  221,22,1); // eta     
444         ForceParticleDecay(  113,22,1); // rho     
445         ForceParticleDecay(  223,22,1); // omega     
446         ForceParticleDecay(  331,22,1); // etaprime     
447         ForceParticleDecay(  333,22,1); // phi     
448         break;
449     case kBeautyUpgrade:
450         ForceBeautyUpgrade();
451       break;
452     }
453 }
454
455 void  AliDecayerPythia::SwitchOffHeavyFlavour()
456 {
457     // Switch off heavy flavour production
458     //
459     // Maximum number of quark flavours used in pdf 
460     fPythia->SetMSTP(58, 3);
461     // Maximum number of flavors that can be used in showers
462     fPythia->SetMSTJ(45, 3);    
463     // Switch off g->QQbar splitting in decay table
464     for (Int_t i = 156; i <= 160; i++) fPythia->SetMDME(i, 1, 0);
465 }
466
467 void  AliDecayerPythia::ForceBeautyUpgrade()
468 {
469  //
470  // Force dedicated decay channels of signals ineresting
471  // for the ITS upgrade (Lb, Lc, Xi_c, B)
472  //
473
474    ForceParticleDecay( 5122, 4122, 1);
475    ForceParticleDecay( 4232, 3312, 1);
476    const Int_t prod[2]={421,211};
477    Int_t mult[2]={1,1};
478    ForceParticleDecay(521,prod,mult,2,1);
479    ForceHadronicD(1);
480 }
481
482 void  AliDecayerPythia::Lu1Ent(Int_t flag, Int_t idpart, 
483                       Double_t mom, Double_t theta, Double_t phi)
484 {
485 // Wrap of Pythia lu1ent subroutine
486 //  printf("%d %d %f %f %f\n",flag, idpart, mom, theta, phi);
487   py1ent(flag, idpart, mom, theta, phi);
488   
489 }
490
491
492
493 Int_t AliDecayerPythia::CountProducts(Int_t channel, Int_t particle)
494 {
495 // Count number of decay products
496     Int_t np=0;
497     for (Int_t i=1; i<=5; i++) {
498         if (TMath::Abs(fPythia->GetKFDP(channel,i)) == particle) np++;
499     }
500     return np;
501 }
502
503
504 void AliDecayerPythia::ForceHadronicD(Int_t optUse4Bodies)
505 {
506 //
507 // Force golden D decay modes
508 //
509   const Int_t kNHadrons = 5;
510   Int_t channel;
511   Int_t hadron[kNHadrons] = {411,  421, 431, 4112, 4122};
512
513   // for D+ -> K0* (-> K- pi+) pi+
514   Int_t iKstar0    =  313;
515   Int_t iKstarbar0 = -313;
516   Int_t products[2] = {kKPlus, kPiMinus}, mult[2] = {1, 1};
517   ForceParticleDecay(iKstar0, products, mult, 2);
518   // for Ds -> Phi pi+
519   Int_t iPhi=333;
520   ForceParticleDecay(iPhi,kKPlus,2); // Phi->K+K-
521   // for D0 -> rho0 pi+ k-
522   Int_t iRho0=113;
523   ForceParticleDecay(iRho0,kPiPlus,2); // Rho0->pi+pi-
524   // for Lambda_c -> Delta++ K-
525   Int_t iDeltaPP = 2224;
526   Int_t productsD[2] = {kProton, kPiPlus}, multD[2] = {1, 1};
527   ForceParticleDecay(iDeltaPP, productsD, multD, 2);
528   // for Lambda_c -> Lambda(1520) pi+ -> p K- pi+ 
529   Int_t iLambda1520 = 3124;
530   Int_t productsL[2] = {kProton, kKMinus}, multL[2] = {1, 1};
531   ForceParticleDecay(iLambda1520, productsL, multL, 2);
532   // for Lambda_c -> Lambda pi+
533   Int_t iLambda=3122;
534   //for Lambda_c -> antiK0 p
535   Int_t iK0bar=-311;
536
537
538   Int_t decayP1[kNHadrons][4] = 
539     { 
540       {kKMinus, kPiPlus,    kPiPlus, 0},
541       {kKMinus, kPiPlus,    0      , 0},
542       {kKPlus , iKstarbar0, 0      , 0},
543       {-1     , -1        , -1     , -1},
544       {kProton, iKstarbar0, 0      , 0}
545     };
546   Int_t decayP2[kNHadrons][4] = 
547     { 
548       {iKstarbar0, kPiPlus, 0      , 0},
549       {kKMinus   , kPiPlus, kPiPlus, kPiMinus},
550       {iPhi      , kPiPlus, 0      , 0},
551       {-1        , -1     , -1     , -1},
552       {iDeltaPP  , kKMinus, 0      , 0}
553     };
554   Int_t decayP3[kNHadrons][4] = 
555     { 
556       {kPiPlus   , iPhi   , 0       , 0},
557       {kKMinus   , kPiPlus, iRho0   , 0 },
558       {-1        , -1     , -1      , -1},
559       {-1        , -1     , -1      , -1},
560       {kProton   , kKMinus, kPiPlus , 0}
561     };
562   // for Lambda_c -> Lambda_1520 pi+ -> p K- pi+, D0-> K*0 pi+ pi- -> K3pi
563     Int_t decayP4[kNHadrons][4] =
564     {
565       {iKstarbar0  , kKPlus  , 0       ,  0},
566       {iKstarbar0  , kPiPlus , kPiMinus,  0},
567       {-1          , -1      , -1      , -1},
568       {-1          , -1      , -1      , -1},
569       {iLambda1520 , kPiPlus ,  0      ,  0}
570     };
571   // for Lambda_c -> Lambda pi+ 
572     Int_t decayP5[kNHadrons][4] =
573     {
574       {-1        , -1     , -1      , -1},
575       {-1        , -1     , -1      , -1},
576       {-1        , -1     , -1      , -1},
577       {-1        , -1     , -1      , -1},
578       {iLambda   , kPiPlus,  0      ,  0}
579     };
580
581     // for Lambda_c -> K0bar p
582     Int_t decayP6[kNHadrons][4] =
583     {
584       {-1        , -1     , -1      , -1},
585       {-1        , -1     , -1      , -1},
586       {-1        , -1     , -1      , -1},
587       {-1        , -1     , -1      , -1},
588       {kProton    , iK0bar,  0      ,  0}
589     };
590
591   if(optUse4Bodies==0){
592     for(Int_t iDau=0;iDau<4;iDau++){
593       decayP2[1][iDau]=-1;
594       decayP3[1][iDau]=-1;
595       decayP4[1][iDau]=-1; 
596     }
597   }  
598   
599   for (Int_t ihadron = 0; ihadron < kNHadrons; ihadron++)
600     {
601       Int_t kc = fPythia->Pycomp(hadron[ihadron]);      
602       Int_t ifirst = fPythia->GetMDCY(kc,2);
603       Int_t ilast  = ifirst + fPythia->GetMDCY(kc,3)-1;
604       Double_t norm  = 0.;
605       for (channel=ifirst; channel<=ilast;channel++) norm+=fPythia->GetBRAT(channel);
606       if (norm < 1.-1.e-12 || norm > 1.+1.e-12) {
607         char pName[16];
608         fPythia->Pyname(hadron[ihadron],pName);
609         AliWarning(Form("Total branching ratio of %s (PDG code = %d) not equal to 1 (= %f)",pName,hadron[ihadron],norm));
610       }
611       fBraPart[kc] = norm;
612       fPythia->SetMDCY(kc,1,1); 
613
614       for (channel = ifirst; channel <= ilast; channel++) {
615         if ((
616              fPythia->GetKFDP(channel,1) == decayP1[ihadron][0] &&
617              fPythia->GetKFDP(channel,2) == decayP1[ihadron][1] &&
618              fPythia->GetKFDP(channel,3) == decayP1[ihadron][2] &&
619              fPythia->GetKFDP(channel,4) == decayP1[ihadron][3] &&
620              fPythia->GetKFDP(channel,5) == 0
621        ) || (
622              fPythia->GetKFDP(channel,1) == decayP2[ihadron][0] &&
623              fPythia->GetKFDP(channel,2) == decayP2[ihadron][1] &&
624              fPythia->GetKFDP(channel,3) == decayP2[ihadron][2] &&
625              fPythia->GetKFDP(channel,4) == decayP2[ihadron][3] &&
626              fPythia->GetKFDP(channel,5) == 0
627         ) || (
628               fPythia->GetKFDP(channel,1) == decayP3[ihadron][0] &&
629               fPythia->GetKFDP(channel,2) == decayP3[ihadron][1] &&
630               fPythia->GetKFDP(channel,3) == decayP3[ihadron][2] &&
631               fPythia->GetKFDP(channel,4) == decayP3[ihadron][3] &&
632               fPythia->GetKFDP(channel,5) == 0
633         ) || (
634               fPythia->GetKFDP(channel,1) == decayP4[ihadron][0] &&
635               fPythia->GetKFDP(channel,2) == decayP4[ihadron][1] &&
636               fPythia->GetKFDP(channel,3) == decayP4[ihadron][2] &&
637               fPythia->GetKFDP(channel,4) == decayP4[ihadron][3] &&
638               fPythia->GetKFDP(channel,5) == 0
639         ) || (
640               fPythia->GetKFDP(channel,1) == decayP5[ihadron][0] &&
641               fPythia->GetKFDP(channel,2) == decayP5[ihadron][1] &&
642               fPythia->GetKFDP(channel,3) == decayP5[ihadron][2] &&
643               fPythia->GetKFDP(channel,4) == decayP5[ihadron][3] &&
644               fPythia->GetKFDP(channel,5) == 0
645         ) || (
646               fPythia->GetKFDP(channel,1) == decayP6[ihadron][0] &&
647               fPythia->GetKFDP(channel,2) == decayP6[ihadron][1] &&
648               fPythia->GetKFDP(channel,3) == decayP6[ihadron][2] &&
649               fPythia->GetKFDP(channel,4) == decayP6[ihadron][3] &&
650               fPythia->GetKFDP(channel,5) == 0
651
652               ))
653           
654           {
655             fPythia->SetMDME(channel,1,1);
656           } else {
657             fPythia->SetMDME(channel,1,0);
658             fBraPart[kc] -= fPythia->GetBRAT(channel);
659           } // selected channel ?
660       } // decay channels
661       if (norm > 0) fBraPart[kc] /= norm;
662     } // hadrons
663 }
664
665
666
667 void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t product, Int_t mult)
668 {
669 //
670 //  Force decay of particle into products with multiplicity mult
671
672     Int_t kc=fPythia->Pycomp(particle);
673     fPythia->SetMDCY(kc,1,1);
674     Int_t ifirst=fPythia->GetMDCY(kc,2);
675     Int_t ilast=ifirst+fPythia->GetMDCY(kc,3)-1;
676     Double_t norm = 0.;
677     for (Int_t channel=ifirst; channel<=ilast;channel++) norm+=fPythia->GetBRAT(channel);
678     if (norm < 1.-1.e-12 || norm > 1.+1.e-12) {
679         char pName[16];
680         fPythia->Pyname(particle,pName);
681         AliWarning(Form("Total branching ratio of %s (PDG code = %d) not equal to 1 (= %f)",pName,particle,norm));
682     }
683     fBraPart[kc] = norm;
684 //
685 //  Loop over decay channels
686     for (Int_t channel=ifirst; channel<=ilast;channel++) {
687         if (CountProducts(channel,product) >= mult) {
688             fPythia->SetMDME(channel,1,1);
689         } else {
690             fPythia->SetMDME(channel,1,0);
691             fBraPart[kc] -= fPythia->GetBRAT(channel);
692         }
693     }
694     if (norm > 0.) fBraPart[kc] /= norm;
695 }
696
697 void AliDecayerPythia::ForceParticleDecay(Int_t particle, const Int_t* products, Int_t* mult, Int_t npart, Bool_t flag)
698 {
699 //
700 //  Force decay of particle into products with multiplicity mult
701
702     Int_t kc=fPythia->Pycomp(particle);
703     fPythia->SetMDCY(kc,1,1);
704     Int_t ifirst=fPythia->GetMDCY(kc,2);
705     Int_t ilast=ifirst+fPythia->GetMDCY(kc,3)-1;
706     Double_t norm = 0.;
707     for (Int_t channel=ifirst; channel<=ilast;channel++) norm+=fPythia->GetBRAT(channel);
708     if (norm < 1.-1.e-12 || norm > 1.+1.e-12) {
709         char pName[16];
710         fPythia->Pyname(particle,pName);
711         AliWarning(Form("Total branching ratio of %s (PDG code = %d) not equal to 1 (= %f)",pName,particle,norm));
712     }
713     fBraPart[kc] = norm;
714 //
715 //  Loop over decay channels
716     for (Int_t channel = ifirst; channel <= ilast; channel++) {
717         Int_t nprod = 0;
718         for (Int_t i = 0; i < npart; i++) {
719             nprod += (CountProducts(channel, products[i]) >= mult[i]);
720         }
721         if ((nprod && !flag) || ((nprod == npart) && flag)) {
722             fPythia->SetMDME(channel,1,1);
723         } else {                // 
724             fPythia->SetMDME(channel,1,0);
725             fBraPart[kc] -= fPythia->GetBRAT(channel);
726         }
727     }
728     if (norm > 0.) fBraPart[kc] /= norm;
729 }
730
731 void AliDecayerPythia::DefineParticles()
732 {
733 //
734 //
735 // Some particles cloned for rare decays     
736 //
737 //  phi-> mu+mu- and phi -> e+e-
738 //  clone the original phi
739
740     Float_t mass;
741     Float_t tlife;
742     Int_t kc, nkc, i;
743
744     kc  = fPythia->Pycomp(333);
745     nkc = 41;
746     
747     for (i=1;i<=3;i++) {
748          fPythia->SetKCHG(nkc,i,fPythia->GetKCHG(kc,i));
749     }
750     
751     for (i=1;i<=4;i++) {
752          fPythia->SetPMAS(nkc,i,fPythia->GetPMAS(kc,i));
753     }
754 //    fPythia->SetCHAF(nkc, fPythia->GetCHAF(kc));
755     fBraPart[kc]=1;
756 //    
757 //  decay
758     fPythia-> SetMDCY(nkc,1,1);
759     fPythia-> SetMDCY(nkc,2,993);
760     fPythia-> SetMDCY(nkc,3,2);
761 //
762 //  phi-> e+e-
763     fPythia->SetMDME(993,1,1);
764     fPythia->SetMDME(993,2,0);
765     fPythia->SetBRAT(993,2.99e-4);
766     fPythia->SetKFDP(993,1,+11);
767     fPythia->SetKFDP(993,2,-11);
768     fPythia->SetKFDP(993,3,0);
769     fPythia->SetKFDP(993,4,0);
770     fPythia->SetKFDP(993,5,0);
771 //
772 //  phi-> mu+mu-
773     fPythia->SetMDME(994,1,1);
774     fPythia->SetMDME(994,2,0);
775     fPythia->SetBRAT(994,2.5e-4);
776     fPythia->SetKFDP(994,1,+13);
777     fPythia->SetKFDP(994,2,-13);
778     fPythia->SetKFDP(994,3,0);
779     fPythia->SetKFDP(994,4,0);
780     fPythia->SetKFDP(994,5,0);
781 //
782 //          Vector mesons
783 //
784 // phi clone for dilepton decay-channel
785     kc =   fPythia->Pycomp(41);     
786     mass = fPythia->GetPMAS(kc,1);
787     tlife= fPythia->GetPMAS(kc,4);
788 //    gMC->Gspart(113,"Phi",3,mass,0,tlife);
789 }
790
791 void  AliDecayerPythia::ForceOmega()
792 {
793     // Force Omega -> Lambda K- Decay
794     Int_t kc=fPythia->Pycomp(3334);
795     fPythia->SetMDCY(kc,1,1);   
796     Int_t ifirst = fPythia->GetMDCY(kc,2);
797     Int_t ilast  = ifirst + fPythia->GetMDCY(kc,3)-1;
798     for (Int_t channel = ifirst; channel <= ilast; channel++) {
799         if (
800             fPythia->GetKFDP(channel,1) == kLambda0 &&
801             fPythia->GetKFDP(channel,2) == kKMinus  &&
802             fPythia->GetKFDP(channel,3) == 0
803             )
804         {
805             fPythia->SetMDME(channel,1,1);
806         } else {
807             fPythia->SetMDME(channel,1,0);
808         } // selected channel ?
809     } // decay channels
810 }
811
812 void  AliDecayerPythia::ForceLambda()
813 {
814     // Force Lambda -> p pi-
815     Int_t kc=fPythia->Pycomp(3122);
816     fPythia->SetMDCY(kc,1,1);   
817     Int_t ifirst = fPythia->GetMDCY(kc,2);
818     Int_t ilast  = ifirst + fPythia->GetMDCY(kc,3)-1;
819     for (Int_t channel = ifirst; channel <= ilast; channel++) {
820         if (
821             fPythia->GetKFDP(channel,1) == kProton   &&
822             fPythia->GetKFDP(channel,2) == kPiMinus  &&
823             fPythia->GetKFDP(channel,3) == 0
824             )
825         {
826             fPythia->SetMDME(channel,1,1);
827         } else {
828             fPythia->SetMDME(channel,1,0);
829         } // selected channel ?
830     } // decay channels
831 }
832
833
834 void AliDecayerPythia::SwitchOffBDecay()
835 {
836 // Switch off B-decays
837   Int_t heavyB[]={511,521,531,5122,5132,5232,5332};
838   for(int i=0;i<4;i++)
839     {
840       fPythia->SetMDCY(fPythia->Pycomp(heavyB[i]),1,0);
841     }
842 }
843
844 Float_t  AliDecayerPythia::GetPartialBranchingRatio(Int_t kf)
845 {
846 // Get branching ratio
847     Int_t kc=fPythia->Pycomp(TMath::Abs(kf));
848     return fBraPart[kc];
849 }
850
851 Float_t  AliDecayerPythia::GetLifetime(Int_t kf)
852 {
853 // Get branching ratio
854     Int_t kc=fPythia->Pycomp(TMath::Abs(kf));
855     return fPythia->GetPMAS(kc,4)*3.3333e-12;
856 }
857
858 void AliDecayerPythia::WriteDecayTable()
859 {
860 //
861 // Write the decay table
862     fPythia->Pyupda(1,15);
863 }
864
865 void AliDecayerPythia::ReadDecayTable()
866 {
867 //
868 // Read the decay table
869     Int_t lun = 15;
870     opendecaytable(lun);
871     fPythia->Pyupda(3,lun);
872     closedecaytable(lun);
873     
874 }
875
876 void AliDecayerPythia::Copy(TObject &) const
877 {
878   //
879   // Copy *this onto AliDecayerPythia -- not implemented
880   //
881   Fatal("Copy","Not implemented!\n");
882 }
883
884
885 /*
886
887                               Particle/parton data table
888
889         KF     KC    particle        antiparticle      chg  col  anti        mass       width       w-cut     lifetime decay
890            IDC on/off ME   Br.rat.    decay products
891
892          1      1    d               dbar               -1    1    1      0.33000     0.00000     0.00000   0.00000E+00    0
893              1    1  102    0.000000    g               d                                                               
894              2    1  102    0.000000    gamma           d                                                               
895              3    1  102    0.000000    Z0              d                                                               
896              4    1  102    0.000000    W-              u                                                               
897              5    1  102    0.000000    W-              c                                                               
898              6    1  102    0.000000    W-              t                                                               
899              7   -1  102    0.000000    W-              t'                                                              
900              8    1  102    0.000000    h0              d                                                               
901
902          2      2    u               ubar                2    1    1      0.33000     0.00000     0.00000   0.00000E+00    0
903              9    1  102    0.000000    g               u                                                               
904             10    1  102    0.000000    gamma           u                                                               
905             11    1  102    0.000000    Z0              u                                                               
906             12    1  102    0.000000    W+              d                                                               
907             13    1  102    0.000000    W+              s                                                               
908             14    1  102    0.000000    W+              b                                                               
909             15   -1  102    0.000000    W+              b'                                                              
910             16    1  102    0.000000    h0              u                                                               
911
912          3      3    s               sbar               -1    1    1      0.50000     0.00000     0.00000   0.00000E+00    0
913             17    1  102    0.000000    g               s                                                               
914             18    1  102    0.000000    gamma           s                                                               
915             19    1  102    0.000000    Z0              s                                                               
916             20    1  102    0.000000    W-              u                                                               
917             21    1  102    0.000000    W-              c                                                               
918             22    1  102    0.000000    W-              t                                                               
919             23   -1  102    0.000000    W-              t'                                                              
920             24    1  102    0.000000    h0              s                                                               
921
922          4      4    c               cbar                2    1    1      1.20000     0.00000     0.00000   0.00000E+00    0
923             25    1  102    0.000000    g               c                                                               
924             26    1  102    0.000000    gamma           c                                                               
925             27    1  102    0.000000    Z0              c                                                               
926             28    1  102    0.000000    W+              d                                                               
927             29    1  102    0.000000    W+              s                                                               
928             30    1  102    0.000000    W+              b                                                               
929             31   -1  102    0.000000    W+              b'                                                              
930             32    1  102    0.000000    h0              c                                                               
931
932          5      5    b               bbar               -1    1    1      4.80000     0.00000     0.00000   0.00000E+00    0
933             33    1  102    0.000000    g               b                                                               
934             34    1  102    0.000000    gamma           b                                                               
935             35    1  102    0.000000    Z0              b                                                               
936             36    1  102    0.000000    W-              u                                                               
937             37    1  102    0.000000    W-              c                                                               
938             38    1  102    0.000000    W-              t                                                               
939             39   -1  102    0.000000    W-              t'                                                              
940             40    1  102    0.000000    h0              b                                                               
941
942          6      6    t               tbar                2    1    1    175.00000     1.40244    14.02444   0.00000E+00    1
943             41    1  102    0.000000    g               t                                                               
944             42    1  102    0.000000    gamma           t                                                               
945             43    1  102    0.000000    Z0              t                                                               
946             44    1    0    0.000030    W+              d                                                               
947             45    1    0    0.001765    W+              s                                                               
948             46    1    0    0.998205    W+              b                                                               
949             47   -1    0    0.000000    W+              b'                                                              
950             48    1  102    0.000000    h0              t                                                               
951             49   -1    0    0.000000    H+              b                                                               
952             50   -1   53    0.000000    ~chi_10         ~t_1                                                            
953             51   -1   53    0.000000    ~chi_20         ~t_1                                                            
954             52   -1   53    0.000000    ~chi_30         ~t_1                                                            
955             53   -1   53    0.000000    ~chi_40         ~t_1                                                            
956             54   -1   53    0.000000    ~g              ~t_1                                                            
957             55   -1   53    0.000000    ~gravitino      ~t_1                                                            
958
959          7      7    b'              b'bar              -1    1    1    400.00000     0.00000     0.00000   0.00000E+00    1
960             56    1  102    0.000000    g               b'                                                              
961             57    1  102    0.000000    gamma           b'                                                              
962             58    1  102    0.000000    Z0              b'                                                              
963             59    1    0    0.000000    W-              u                                                               
964             60    1    0    0.000000    W-              c                                                               
965             61    1    0    0.000000    W-              t                                                               
966             62    1    0    0.000000    W-              t'                                                              
967             63    1  102    0.000000    h0              b'                                                              
968             64   -1    0    0.000000    H-              c                                                               
969             65   -1    0    0.000000    H-              t                                                               
970
971          8      8    t'              t'bar               2    1    1    400.00000     0.00000     0.00000   0.00000E+00    1
972             66    1  102    0.000000    g               t'                                                              
973             67    1  102    0.000000    gamma           t'                                                              
974             68    1  102    0.000000    Z0              t'                                                              
975             69    1    0    0.000000    W+              d                                                               
976             70    1    0    0.000000    W+              s                                                               
977             71    1    0    0.000000    W+              b                                                               
978             72    1    0    0.000000    W+              b'                                                              
979             73    1  102    0.000000    h0              t'                                                              
980             74   -1    0    0.000000    H+              b                                                               
981             75   -1    0    0.000000    H+              b'                                                              
982
983         11     11    e-              e+                 -3    0    1      0.00051     0.00000     0.00000   0.00000E+00    0
984             76    1  102    0.000000    gamma           e-                                                              
985             77    1  102    0.000000    Z0              e-                                                              
986             78    1  102    0.000000    W-              nu_e                                                            
987             79    1  102    0.000000    h0              e-                                                              
988
989         12     12    nu_e            nu_ebar             0    0    1      0.00000     0.00000     0.00000   0.00000E+00    0
990             80    1  102    0.000000    Z0              nu_e                                                            
991             81    1  102    0.000000    W+              e-                                                              
992
993         13     13    mu-             mu+                -3    0    1      0.10566     0.00000     0.00000   6.58654E+05    0
994             82    1   42    1.000000    nu_ebar         e-              nu_mu                                           
995             83    1  102    0.000000    gamma           mu-                                                             
996             84    1  102    0.000000    Z0              mu-                                                             
997             85    1  102    0.000000    W-              nu_mu                                                           
998             86    1  102    0.000000    h0              mu-                                                             
999
1000         14     14    nu_mu           nu_mubar            0    0    1      0.00000     0.00000     0.00000   0.00000E+00    0
1001             87    1  102    0.000000    Z0              nu_mu                                                           
1002             88    1  102    0.000000    W+              mu-                                                             
1003
1004         15     15    tau-            tau+               -3    0    1      1.77700     0.00000     0.00000   8.72000E-02    1
1005             89    1   42    0.178300    nu_ebar         e-              nu_tau                                          
1006             90    1   42    0.173500    nu_mubar        mu-             nu_tau                                          
1007             91    1    0    0.113100    nu_tau          pi-                                                             
1008             92    1    0    0.249400    nu_tau          rho-                                                            
1009             93    1   41    0.003000    nu_tau          pi-             pi0                                             
1010             94    1   41    0.090000    nu_tau          rho-            pi0                                             
1011             95    1   41    0.002700    nu_tau          pi-             pi0             pi0                             
1012             96    1   41    0.010000    nu_tau          rho-            pi0             pi0                             
1013             97    1   41    0.001400    nu_tau          pi-             pi0             pi0             pi0             
1014             98    1   41    0.001200    nu_tau          rho-            pi0             pi0             pi0             
1015             99    1   41    0.000250    nu_tau          pi-             K_S0                                            
1016            100    1   41    0.000250    nu_tau          pi-             K_L0                                            
1017            101    1    0    0.007100    nu_tau          K-                                                              
1018            102    1    0    0.012000    nu_tau          K*-                                                             
1019            103    1   41    0.000400    nu_tau          K-              pi0                                             
1020            104    1   41    0.000750    nu_tau          K*-             pi0                                             
1021            105    1   41    0.000060    nu_tau          K*-             pi0             pi0                             
1022            106    1   41    0.000780    nu_tau          K-              K_S0                                            
1023            107    1   41    0.000780    nu_tau          K-              K_L0                                            
1024            108    1   41    0.003400    nu_tau          K-              K+              pi-                             
1025            109    1   41    0.080000    nu_tau          pi-             rho0                                            
1026            110    1   41    0.011000    nu_tau          pi-             pi+             pi-                             
1027            111    1   41    0.019100    nu_tau          pi-             omega                                           
1028            112    1   41    0.000060    nu_tau          pi-             eta                                             
1029            113    1   41    0.005000    nu_tau          rho-            rho0                                            
1030            114    1   41    0.013300    nu_tau          pi-             rho0            pi0                             
1031            115    1   41    0.006700    nu_tau          rho-            pi+             pi-                             
1032            116    1   41    0.000500    nu_tau          pi-             pi+             pi-             pi0             
1033            117    1   41    0.003500    nu_tau          rho-            omega                                           
1034            118    1   41    0.000600    nu_tau          pi-             omega           pi0                             
1035            119    1   41    0.001500    nu_tau          rho-            eta                                             
1036            120    1   41    0.000210    nu_tau          pi-             eta             pi0                             
1037            121    1   41    0.000200    nu_tau          rho-            rho0            pi0                             
1038            122    1   41    0.000750    nu_tau          pi-             rho0            rho0                            
1039            123    1   41    0.000100    nu_tau          pi-             eta             eta                             
1040            124    1   41    0.000200    nu_tau          pi-             rho0            pi0             pi0             
1041            125    1   41    0.001100    nu_tau          rho-            rho0            pi0             pi0             
1042            126    1   41    0.000200    nu_tau          pi-             rho+            rho-                            
1043            127    1   41    0.000200    nu_tau          pi-             rho+            pi-             pi0             
1044            128    1   41    0.000200    nu_tau          pi-             rho-            pi+             pi0             
1045            129    1   41    0.000220    nu_tau          pi-             rho0            rho0            pi0             
1046            130    1   41    0.000400    nu_tau          K*-             pi0             pi0                             
1047            131    1   41    0.000100    nu_tau          K-              pi0             pi0             pi0             
1048            132    1   41    0.002050    nu_tau          pi-             K_S0            pi0                             
1049            133    1   41    0.002050    nu_tau          pi-             K_L0            pi0                             
1050            134    1   41    0.000690    nu_tau          K-              K_S0            pi0                             
1051            135    1   41    0.000690    nu_tau          K-              K_L0            pi0                             
1052            136    1   41    0.000250    nu_tau          pi-             K_S0            K_S0                            
1053            137    1   41    0.000510    nu_tau          pi-             K_S0            K_L0                            
1054            138    1   41    0.000250    nu_tau          pi-             K_L0            K_L0                            
1055            139    1  102    0.000000    gamma           tau-                                                            
1056            140    1  102    0.000000    Z0              tau-                                                            
1057            141    1  102    0.000000    W-              nu_tau                                                          
1058            142    1  102    0.000000    h0              tau-                                                            
1059
1060         16     16    nu_tau          nu_taubar           0    0    1      0.00000     0.00000     0.00000   0.00000E+00    0
1061            143    1  102    0.000000    Z0              nu_tau                                                          
1062            144    1  102    0.000000    W+              tau-                                                            
1063
1064         17     17    tau'-           tau'+              -3    0    1    400.00000     0.00000     0.00000   0.00000E+00    1
1065            145    1  102    0.000000    gamma           tau'-                                                           
1066            146    1  102    0.000000    Z0              tau'-                                                           
1067            147    1    0    0.000000    W-              nu'_tau                                                         
1068            148    1  102    0.000000    h0              tau'-                                                           
1069            149   -1    0    0.000000    H-              nu'_tau                                                         
1070
1071         18     18    nu'_tau         nu'_taubar          0    0    1      0.00000     0.00000     0.00000   0.00000E+00    0
1072            150    1  102    0.000000    Z0              nu'_tau                                                         
1073            151    1    0    0.000000    W+              tau'-                                                           
1074            152   -1    0    0.000000    H+              tau'-                                                           
1075
1076         21     21    g                                   0    2    0      0.00000     0.00000     0.00000   0.00000E+00    0
1077            153    0  102    0.000000    d               dbar                                                            
1078            154    0  102    0.000000    u               ubar                                                            
1079            155    0  102    0.000000    s               sbar                                                            
1080            156    1  102    0.000000    c               cbar                                                            
1081            157    0  102    0.000000    b               bbar                                                            
1082            158    0  102    0.000000    t               tbar                                                            
1083            159    0  102    0.000000    b'              b'bar                                                           
1084            160    0  102    0.000000    t'              t'bar                                                           
1085            161    1  102    0.000000    g               g                                                               
1086
1087         22     22    gamma                               0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1088            162    0  102    0.000000    d               dbar                                                            
1089            163    0  102    0.000000    u               ubar                                                            
1090            164    0  102    0.000000    s               sbar                                                            
1091            165    1  102    0.000000    c               cbar                                                            
1092            166    0  102    0.000000    b               bbar                                                            
1093            167    0  102    0.000000    t               tbar                                                            
1094            168    0  102    0.000000    b'              b'bar                                                           
1095            169    0  102    0.000000    t'              t'bar                                                           
1096            170    0  102    0.000000    e-              e+                                                              
1097            171    0  102    0.000000    mu-             mu+                                                             
1098            172    0  102    0.000000    tau-            tau+                                                            
1099            173    0  102    0.000000    tau'-           tau'+                                                           
1100
1101         23     23    Z0                                  0    0    0     91.18700     2.47872    24.78720   0.00000E+00    1
1102            174    1   32    0.153998    d               dbar                                                            
1103            175    1   32    0.119422    u               ubar                                                            
1104            176    1   32    0.153988    s               sbar                                                            
1105            177    1   32    0.119322    c               cbar                                                            
1106            178    1   32    0.152275    b               bbar                                                            
1107            179    1   32    0.000000    t               tbar                                                            
1108            180   -1   32    0.000000    b'              b'bar                                                           
1109            181   -1   32    0.000000    t'              t'bar                                                           
1110            182    1    0    0.033568    e-              e+                                                              
1111            183    1    0    0.066789    nu_e            nu_ebar                                                         
1112            184    1    0    0.033567    mu-             mu+                                                             
1113            185    1    0    0.066789    nu_mu           nu_mubar                                                        
1114            186    1    0    0.033492    tau-            tau+                                                            
1115            187    1    0    0.066789    nu_tau          nu_taubar                                                       
1116            188   -1    0    0.000000    tau'-           tau'+                                                           
1117            189   -1    0    0.000000    nu'_tau         nu'_taubar                                                      
1118
1119         24     24    W+              W-                  3    0    1     80.33000     2.06856    20.68560   0.00000E+00    1
1120            190    1   32    0.321379    dbar            u                                                               
1121            191    1   32    0.016498    dbar            c                                                               
1122            192    1   32    0.000000    dbar            t                                                               
1123            193   -1   32    0.000000    dbar            t'                                                              
1124            194    1   32    0.016502    sbar            u                                                               
1125            195    1   32    0.320685    sbar            c                                                               
1126            196    1   32    0.000000    sbar            t                                                               
1127            197   -1   32    0.000000    sbar            t'                                                              
1128            198    1   32    0.000010    bbar            u                                                               
1129            199    1   32    0.000591    bbar            c                                                               
1130            200    1   32    0.000000    bbar            t                                                               
1131            201   -1   32    0.000000    bbar            t'                                                              
1132            202   -1   32    0.000000    b'bar           u                                                               
1133            203   -1   32    0.000000    b'bar           c                                                               
1134            204   -1   32    0.000000    b'bar           t                                                               
1135            205   -1   32    0.000000    b'bar           t'                                                              
1136            206    1    0    0.108138    e+              nu_e                                                            
1137            207    1    0    0.108138    mu+             nu_mu                                                           
1138            208    1    0    0.108059    tau+            nu_tau                                                          
1139            209   -1    0    0.000000    tau'+           nu'_tau                                                         
1140
1141         25     25    h0                                  0    0    0     80.00000     0.00240     0.02402   0.00000E+00    1
1142            210    1   32    0.000001    d               dbar                                                            
1143            211    1   32    0.000000    u               ubar                                                            
1144            212    1   32    0.000378    s               sbar                                                            
1145            213    1   32    0.054441    c               cbar                                                            
1146            214    1   32    0.853399    b               bbar                                                            
1147            215    1   32    0.000000    t               tbar                                                            
1148            216   -1   32    0.000000    b'              b'bar                                                           
1149            217   -1   32    0.000000    t'              t'bar                                                           
1150            218    1    0    0.000000    e-              e+                                                              
1151            219    1    0    0.000241    mu-             mu+                                                             
1152            220    1    0    0.067867    tau-            tau+                                                            
1153            221   -1    0    0.000000    tau'-           tau'+                                                           
1154            222    1    0    0.022178    g               g                                                               
1155            223    1    0    0.000867    gamma           gamma                                                           
1156            224    1    0    0.000000    gamma           Z0                                                              
1157            225    1    0    0.000134    Z0              Z0                                                              
1158            226    1    0    0.000494    W+              W-                                                              
1159            227   -1   53    0.000000    ~chi_10         ~chi_10                                                         
1160            228   -1   53    0.000000    ~chi_20         ~chi_10                                                         
1161            229   -1   53    0.000000    ~chi_20         ~chi_20                                                         
1162            230   -1   53    0.000000    ~chi_30         ~chi_10                                                         
1163            231   -1   53    0.000000    ~chi_30         ~chi_20                                                         
1164            232   -1   53    0.000000    ~chi_30         ~chi_30                                                         
1165            233   -1   53    0.000000    ~chi_40         ~chi_10                                                         
1166            234   -1   53    0.000000    ~chi_40         ~chi_20                                                         
1167            235   -1   53    0.000000    ~chi_40         ~chi_30                                                         
1168            236   -1   53    0.000000    ~chi_40         ~chi_40                                                         
1169            237   -1   53    0.000000    ~chi_1+         ~chi_1-                                                         
1170            238   -1   53    0.000000    ~chi_1+         ~chi_2-                                                         
1171            239   -1   53    0.000000    ~chi_2+         ~chi_1-                                                         
1172            240   -1   53    0.000000    ~chi_2+         ~chi_2-                                                         
1173            241   -1   53    0.000000    ~d_L            ~d_Lbar                                                         
1174            242   -1   53    0.000000    ~d_R            ~d_Rbar                                                         
1175            243   -1   53    0.000000    ~d_L            ~d_Rbar                                                         
1176            244   -1   53    0.000000    ~d_Lbar         ~d_R                                                            
1177            245   -1   53    0.000000    ~u_L            ~u_Lbar                                                         
1178            246   -1   53    0.000000    ~u_R            ~u_Rbar                                                         
1179            247   -1   53    0.000000    ~u_L            ~u_Rbar                                                         
1180            248   -1   53    0.000000    ~u_Lbar         ~u_R                                                            
1181            249   -1   53    0.000000    ~s_L            ~s_Lbar                                                         
1182            250   -1   53    0.000000    ~s_R            ~s_Rbar                                                         
1183            251   -1   53    0.000000    ~s_L            ~s_Rbar                                                         
1184            252   -1   53    0.000000    ~s_Lbar         ~s_R                                                            
1185            253   -1   53    0.000000    ~c_L            ~c_Lbar                                                         
1186            254   -1   53    0.000000    ~c_R            ~c_Rbar                                                         
1187            255   -1   53    0.000000    ~c_L            ~c_Rbar                                                         
1188            256   -1   53    0.000000    ~c_Lbar         ~c_R                                                            
1189            257   -1   53    0.000000    ~b_1            ~b_1bar                                                         
1190            258   -1   53    0.000000    ~b_2            ~b_2bar                                                         
1191            259   -1   53    0.000000    ~b_1            ~b_2bar                                                         
1192            260   -1   53    0.000000    ~b_1bar         ~b_2                                                            
1193            261   -1   53    0.000000    ~t_1            ~t_1bar                                                         
1194            262   -1   53    0.000000    ~t_2            ~t_2bar                                                         
1195            263   -1   53    0.000000    ~t_1            ~t_2bar                                                         
1196            264   -1   53    0.000000    ~t_1bar         ~t_2                                                            
1197            265   -1   53    0.000000    ~e_L-           ~e_L+                                                           
1198            266   -1   53    0.000000    ~e_R-           ~e_R+                                                           
1199            267   -1   53    0.000000    ~e_L-           ~e_R+                                                           
1200            268   -1   53    0.000000    ~e_L+           ~e_R-                                                           
1201            269   -1   53    0.000000    ~nu_eL          ~nu_eLbar                                                       
1202            270   -1   53    0.000000    ~nu_eR          ~nu_eRbar                                                       
1203            271   -1   53    0.000000    ~nu_eL          ~nu_eRbar                                                       
1204            272   -1   53    0.000000    ~nu_eLbar       ~nu_eR                                                          
1205            273   -1   53    0.000000    ~mu_L-          ~mu_L+                                                          
1206            274   -1   53    0.000000    ~mu_R-          ~mu_R+                                                          
1207            275   -1   53    0.000000    ~mu_L-          ~mu_R+                                                          
1208            276   -1   53    0.000000    ~mu_L+          ~mu_R-                                                          
1209            277   -1   53    0.000000    ~nu_muL         ~nu_muLbar                                                      
1210            278   -1   53    0.000000    ~nu_muR         ~nu_muRbar                                                      
1211            279   -1   53    0.000000    ~nu_muL         ~nu_muRbar                                                      
1212            280   -1   53    0.000000    ~nu_muLbar      ~nu_muR                                                         
1213            281   -1   53    0.000000    ~tau_1-         ~tau_1+                                                         
1214            282   -1   53    0.000000    ~tau_2-         ~tau_2+                                                         
1215            283   -1   53    0.000000    ~tau_1-         ~tau_2+                                                         
1216            284   -1   53    0.000000    ~tau_1+         ~tau_2-                                                         
1217            285   -1   53    0.000000    ~nu_tauL        ~nu_tauLbar                                                     
1218            286   -1   53    0.000000    ~nu_tauR        ~nu_tauRbar                                                     
1219            287   -1   53    0.000000    ~nu_tauL        ~nu_tauRbar                                                     
1220            288   -1   53    0.000000    ~nu_tauLbar     ~nu_tauR                                                        
1221
1222         28     28    reggeon                             0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1223
1224         29     29    pomeron                             0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1225
1226         32     32    Z'0                                 0    0    0    500.00000    14.54208   145.42084   0.00000E+00    1
1227            289    1   32    0.145842    d               dbar                                                            
1228            290    1   32    0.113282    u               ubar                                                            
1229            291    1   32    0.145842    s               sbar                                                            
1230            292    1   32    0.113278    c               cbar                                                            
1231            293    1   32    0.145788    b               bbar                                                            
1232            294    1   32    0.049004    t               tbar                                                            
1233            295   -1   32    0.000000    b'              b'bar                                                           
1234            296   -1   32    0.000000    t'              t'bar                                                           
1235            297    1    0    0.032021    e-              e+                                                              
1236            298    1    0    0.063634    nu_e            nu_ebar                                                         
1237            299    1    0    0.032021    mu-             mu+                                                             
1238            300    1    0    0.063634    nu_mu           nu_mubar                                                        
1239            301    1    0    0.032018    tau-            tau+                                                            
1240            302    1    0    0.063634    nu_tau          nu_taubar                                                       
1241            303   -1    0    0.000000    tau'-           tau'+                                                           
1242            304   -1    0    0.000000    nu'_tau         nu'_taubar                                                      
1243            305   -1    0    0.000000    W+              W-                                                              
1244            306   -1    0    0.000000    H+              H-                                                              
1245            307   -1    0    0.000000    Z0              gamma                                                           
1246            308   -1    0    0.000000    Z0              h0                                                              
1247            309   -1    0    0.000000    h0              A0                                                              
1248            310   -1    0    0.000000    H0              A0                                                              
1249
1250         33     33    Z"0                                 0    0    0    900.00000     0.00000     0.00000   0.00000E+00    0
1251
1252         34     34    W'+             W'-                 3    0    1    500.00000    16.66312   166.63122   0.00000E+00    1
1253            311    1   32    0.251235    dbar            u                                                               
1254            312    1   32    0.012901    dbar            c                                                               
1255            313    1   32    0.000006    dbar            t                                                               
1256            314   -1   32    0.000000    dbar            t'                                                              
1257            315    1   32    0.012901    sbar            u                                                               
1258            316    1   32    0.250776    sbar            c                                                               
1259            317    1   32    0.000380    sbar            t                                                               
1260            318   -1   32    0.000000    sbar            t'                                                              
1261            319    1   32    0.000008    bbar            u                                                               
1262            320    1   32    0.000465    bbar            c                                                               
1263            321    1   32    0.215427    bbar            t                                                               
1264            322   -1   32    0.000000    bbar            t'                                                              
1265            323   -1   32    0.000000    b'bar           u                                                               
1266            324   -1   32    0.000000    b'bar           c                                                               
1267            325   -1   32    0.000000    b'bar           t                                                               
1268            326   -1   32    0.000000    b'bar           t'                                                              
1269            327    1    0    0.085301    e+              nu_e                                                            
1270            328    1    0    0.085301    mu+             nu_mu                                                           
1271            329    1    0    0.085299    tau+            nu_tau                                                          
1272            330   -1    0    0.000000    tau'+           nu'_tau                                                         
1273            331   -1    0    0.000000    W+              Z0                                                              
1274            332   -1    0    0.000000    W+              gamma                                                           
1275            333   -1    0    0.000000    W+              h0                                                              
1276
1277         35     35    H0                                  0    0    0    300.00000     8.42840    84.28402   0.00000E+00    1
1278            334    1   32    0.000000    d               dbar                                                            
1279            335    1   32    0.000000    u               ubar                                                            
1280            336    1   32    0.000000    s               sbar                                                            
1281            337    1   32    0.000048    c               cbar                                                            
1282            338    1   32    0.000768    b               bbar                                                            
1283            339    1   32    0.000000    t               tbar                                                            
1284            340   -1   32    0.000000    b'              b'bar                                                           
1285            341   -1   32    0.000000    t'              t'bar                                                           
1286            342    1    0    0.000000    e-              e+                                                              
1287            343    1    0    0.000000    mu-             mu+                                                             
1288            344    1    0    0.000074    tau-            tau+                                                            
1289            345   -1    0    0.000000    tau'-           tau'+                                                           
1290            346    1    0    0.000422    g               g                                                               
1291            347    1    0    0.000015    gamma           gamma                                                           
1292            348    1    0    0.000061    gamma           Z0                                                              
1293            349    1    0    0.306171    Z0              Z0                                                              
1294            350    1    0    0.688641    W+              W-                                                              
1295            351    1    0    0.000000    Z0              h0                                                              
1296            352    1    0    0.003799    h0              h0                                                              
1297            353    1    0    0.000000    A0              A0                                                              
1298            354   -1   53    0.000000    ~chi_10         ~chi_10                                                         
1299            355   -1   53    0.000000    ~chi_20         ~chi_10                                                         
1300            356   -1   53    0.000000    ~chi_20         ~chi_20                                                         
1301            357   -1   53    0.000000    ~chi_30         ~chi_10                                                         
1302            358   -1   53    0.000000    ~chi_30         ~chi_20                                                         
1303            359   -1   53    0.000000    ~chi_30         ~chi_30                                                         
1304            360   -1   53    0.000000    ~chi_40         ~chi_10                                                         
1305            361   -1   53    0.000000    ~chi_40         ~chi_20                                                         
1306            362   -1   53    0.000000    ~chi_40         ~chi_30                                                         
1307            363   -1   53    0.000000    ~chi_40         ~chi_40                                                         
1308            364   -1   53    0.000000    ~chi_1+         ~chi_1-                                                         
1309            365   -1   53    0.000000    ~chi_1+         ~chi_2-                                                         
1310            366   -1   53    0.000000    ~chi_2+         ~chi_1-                                                         
1311            367   -1   53    0.000000    ~chi_2+         ~chi_2-                                                         
1312            368   -1   53    0.000000    ~d_L            ~d_Lbar                                                         
1313            369   -1   53    0.000000    ~d_R            ~d_Rbar                                                         
1314            370   -1   53    0.000000    ~d_L            ~d_Rbar                                                         
1315            371   -1   53    0.000000    ~d_Lbar         ~d_R                                                            
1316            372   -1   53    0.000000    ~u_L            ~u_Lbar                                                         
1317            373   -1   53    0.000000    ~u_R            ~u_Rbar                                                         
1318            374   -1   53    0.000000    ~u_L            ~u_Rbar                                                         
1319            375   -1   53    0.000000    ~u_Lbar         ~u_R                                                            
1320            376   -1   53    0.000000    ~s_L            ~s_Lbar                                                         
1321            377   -1   53    0.000000    ~s_R            ~s_Rbar                                                         
1322            378   -1   53    0.000000    ~s_L            ~s_Rbar                                                         
1323            379   -1   53    0.000000    ~s_Lbar         ~s_R                                                            
1324            380   -1   53    0.000000    ~c_L            ~c_Lbar                                                         
1325            381   -1   53    0.000000    ~c_R            ~c_Rbar                                                         
1326            382   -1   53    0.000000    ~c_L            ~c_Rbar                                                         
1327            383   -1   53    0.000000    ~c_Lbar         ~c_R                                                            
1328            384   -1   53    0.000000    ~b_1            ~b_1bar                                                         
1329            385   -1   53    0.000000    ~b_2            ~b_2bar                                                         
1330            386   -1   53    0.000000    ~b_1            ~b_2bar                                                         
1331            387   -1   53    0.000000    ~b_1bar         ~b_2                                                            
1332            388   -1   53    0.000000    ~t_1            ~t_1bar                                                         
1333            389   -1   53    0.000000    ~t_2            ~t_2bar                                                         
1334            390   -1   53    0.000000    ~t_1            ~t_2bar                                                         
1335            391   -1   53    0.000000    ~t_1bar         ~t_2                                                            
1336            392   -1   53    0.000000    ~e_L-           ~e_L+                                                           
1337            393   -1   53    0.000000    ~e_R-           ~e_R+                                                           
1338            394   -1   53    0.000000    ~e_L-           ~e_R+                                                           
1339            395   -1   53    0.000000    ~e_L+           ~e_R-                                                           
1340            396   -1   53    0.000000    ~nu_eL          ~nu_eLbar                                                       
1341            397   -1   53    0.000000    ~nu_eR          ~nu_eRbar                                                       
1342            398   -1   53    0.000000    ~nu_eL          ~nu_eRbar                                                       
1343            399   -1   53    0.000000    ~nu_eLbar       ~nu_eR                                                          
1344            400   -1   53    0.000000    ~mu_L-          ~mu_L+                                                          
1345            401   -1   53    0.000000    ~mu_R-          ~mu_R+                                                          
1346            402   -1   53    0.000000    ~mu_L-          ~mu_R+                                                          
1347            403   -1   53    0.000000    ~mu_L+          ~mu_R-                                                          
1348            404   -1   53    0.000000    ~nu_muL         ~nu_muLbar                                                      
1349            405   -1   53    0.000000    ~nu_muR         ~nu_muRbar                                                      
1350            406   -1   53    0.000000    ~nu_muL         ~nu_muRbar                                                      
1351            407   -1   53    0.000000    ~nu_muLbar      ~nu_muR                                                         
1352            408   -1   53    0.000000    ~tau_1-         ~tau_1+                                                         
1353            409   -1   53    0.000000    ~tau_2-         ~tau_2+                                                         
1354            410   -1   53    0.000000    ~tau_1-         ~tau_2+                                                         
1355            411   -1   53    0.000000    ~tau_1+         ~tau_2-                                                         
1356            412   -1   53    0.000000    ~nu_tauL        ~nu_tauLbar                                                     
1357            413   -1   53    0.000000    ~nu_tauR        ~nu_tauRbar                                                     
1358            414   -1   53    0.000000    ~nu_tauL        ~nu_tauRbar                                                     
1359            415   -1   53    0.000000    ~nu_tauLbar     ~nu_tauR                                                        
1360
1361         36     36    A0                                  0    0    0    300.00000     4.91995    49.19946   0.00000E+00    1
1362            416    1   32    0.000000    d               dbar                                                            
1363            417    1   32    0.000000    u               ubar                                                            
1364            418    1   32    0.000001    s               sbar                                                            
1365            419    1   32    0.000082    c               cbar                                                            
1366            420    1   32    0.001318    b               bbar                                                            
1367            421    1   32    0.000000    t               tbar                                                            
1368            422   -1   32    0.000000    b'              b'bar                                                           
1369            423   -1   32    0.000000    t'              t'bar                                                           
1370            424    1    0    0.000000    e-              e+                                                              
1371            425    1    0    0.000000    mu-             mu+                                                             
1372            426    1    0    0.000126    tau-            tau+                                                            
1373            427   -1    0    0.000000    tau'-           tau'+                                                           
1374            428    1    0    0.002164    g               g                                                               
1375            429    1    0    0.000010    gamma           gamma                                                           
1376            430    1    0    0.000002    gamma           Z0                                                              
1377            431    1    0    0.000000    Z0              Z0                                                              
1378            432    1    0    0.000000    W+              W-                                                              
1379            433    1    0    0.996296    Z0              h0                                                              
1380            434   -1   53    0.000000    ~chi_10         ~chi_10                                                         
1381            435   -1   53    0.000000    ~chi_20         ~chi_10                                                         
1382            436   -1   53    0.000000    ~chi_20         ~chi_20                                                         
1383            437   -1   53    0.000000    ~chi_30         ~chi_10                                                         
1384            438   -1   53    0.000000    ~chi_30         ~chi_20                                                         
1385            439   -1   53    0.000000    ~chi_30         ~chi_30                                                         
1386            440   -1   53    0.000000    ~chi_40         ~chi_10                                                         
1387            441   -1   53    0.000000    ~chi_40         ~chi_20                                                         
1388            442   -1   53    0.000000    ~chi_40         ~chi_30                                                         
1389            443   -1   53    0.000000    ~chi_40         ~chi_40                                                         
1390            444   -1   53    0.000000    ~chi_1+         ~chi_1-                                                         
1391            445   -1   53    0.000000    ~chi_1+         ~chi_2-                                                         
1392            446   -1   53    0.000000    ~chi_2+         ~chi_1-                                                         
1393            447   -1   53    0.000000    ~chi_2+         ~chi_2-                                                         
1394            448   -1   53    0.000000    ~d_L            ~d_Lbar                                                         
1395            449   -1   53    0.000000    ~d_R            ~d_Rbar                                                         
1396            450   -1   53    0.000000    ~d_L            ~d_Rbar                                                         
1397            451   -1   53    0.000000    ~d_Lbar         ~d_R                                                            
1398            452   -1   53    0.000000    ~u_L            ~u_Lbar                                                         
1399            453   -1   53    0.000000    ~u_R            ~u_Rbar                                                         
1400            454   -1   53    0.000000    ~u_L            ~u_Rbar                                                         
1401            455   -1   53    0.000000    ~u_Lbar         ~u_R                                                            
1402            456   -1   53    0.000000    ~s_L            ~s_Lbar                                                         
1403            457   -1   53    0.000000    ~s_R            ~s_Rbar                                                         
1404            458   -1   53    0.000000    ~s_L            ~s_Rbar                                                         
1405            459   -1   53    0.000000    ~s_Lbar         ~s_R                                                            
1406            460   -1   53    0.000000    ~c_L            ~c_Lbar                                                         
1407            461   -1   53    0.000000    ~c_R            ~c_Rbar                                                         
1408            462   -1   53    0.000000    ~c_L            ~c_Rbar                                                         
1409            463   -1   53    0.000000    ~c_Lbar         ~c_R                                                            
1410            464   -1   53    0.000000    ~b_1            ~b_1bar                                                         
1411            465   -1   53    0.000000    ~b_2            ~b_2bar                                                         
1412            466   -1   53    0.000000    ~b_1            ~b_2bar                                                         
1413            467   -1   53    0.000000    ~b_1bar         ~b_2                                                            
1414            468   -1   53    0.000000    ~t_1            ~t_1bar                                                         
1415            469   -1   53    0.000000    ~t_2            ~t_2bar                                                         
1416            470   -1   53    0.000000    ~t_1            ~t_2bar                                                         
1417            471   -1   53    0.000000    ~t_1bar         ~t_2                                                            
1418            472   -1   53    0.000000    ~e_L-           ~e_L+                                                           
1419            473   -1   53    0.000000    ~e_R-           ~e_R+                                                           
1420            474   -1   53    0.000000    ~e_L-           ~e_R+                                                           
1421            475   -1   53    0.000000    ~e_L+           ~e_R-                                                           
1422            476   -1   53    0.000000    ~nu_eL          ~nu_eLbar                                                       
1423            477   -1   53    0.000000    ~nu_eR          ~nu_eRbar                                                       
1424            478   -1   53    0.000000    ~nu_eL          ~nu_eRbar                                                       
1425            479   -1   53    0.000000    ~nu_eLbar       ~nu_eR                                                          
1426            480   -1   53    0.000000    ~mu_L-          ~mu_L+                                                          
1427            481   -1   53    0.000000    ~mu_R-          ~mu_R+                                                          
1428            482   -1   53    0.000000    ~mu_L-          ~mu_R+                                                          
1429            483   -1   53    0.000000    ~mu_L+          ~mu_R-                                                          
1430            484   -1   53    0.000000    ~nu_muL         ~nu_muLbar                                                      
1431            485   -1   53    0.000000    ~nu_muR         ~nu_muRbar                                                      
1432            486   -1   53    0.000000    ~nu_muL         ~nu_muRbar                                                      
1433            487   -1   53    0.000000    ~nu_muLbar      ~nu_muR                                                         
1434            488   -1   53    0.000000    ~tau_1-         ~tau_1+                                                         
1435            489   -1   53    0.000000    ~tau_2-         ~tau_2+                                                         
1436            490   -1   53    0.000000    ~tau_1-         ~tau_2+                                                         
1437            491   -1   53    0.000000    ~tau_1+         ~tau_2-                                                         
1438            492   -1   53    0.000000    ~nu_tauL        ~nu_tauLbar                                                     
1439            493   -1   53    0.000000    ~nu_tauR        ~nu_tauRbar                                                     
1440            494   -1   53    0.000000    ~nu_tauL        ~nu_tauRbar                                                     
1441            495   -1   53    0.000000    ~nu_tauLbar     ~nu_tauR                                                        
1442
1443         37     37    H+              H-                  3    0    1    300.00000     5.76067    57.60673   0.00000E+00    1
1444            496    1   32    0.000000    dbar            u                                                               
1445            497    1   32    0.000015    sbar            c                                                               
1446            498    1   32    0.067644    bbar            t                                                               
1447            499   -1   32    0.000000    b'bar           t'                                                              
1448            500    1    0    0.000000    e+              nu_e                                                            
1449            501    1    0    0.000010    mu+             nu_mu                                                           
1450            502    1    0    0.002701    tau+            nu_tau                                                          
1451            503   -1    0    0.000000    tau'+           nu'_tau                                                         
1452            504    1    0    0.929631    W+              h0                                                              
1453            505   -1   53    0.000000    ~chi_10         ~chi_1+                                                         
1454            506   -1   53    0.000000    ~chi_10         ~chi_2+                                                         
1455            507   -1   53    0.000000    ~chi_20         ~chi_1+                                                         
1456            508   -1   53    0.000000    ~chi_20         ~chi_2+                                                         
1457            509   -1   53    0.000000    ~chi_30         ~chi_1+                                                         
1458            510   -1   53    0.000000    ~chi_30         ~chi_2+                                                         
1459            511   -1   53    0.000000    ~chi_40         ~chi_1+                                                         
1460            512   -1   53    0.000000    ~chi_40         ~chi_2+                                                         
1461            513   -1   53    0.000000    ~t_1            ~b_1bar                                                         
1462            514   -1   53    0.000000    ~t_2            ~b_1bar                                                         
1463            515   -1   53    0.000000    ~t_1            ~b_2bar                                                         
1464            516   -1   53    0.000000    ~t_2            ~b_2bar                                                         
1465            517   -1   53    0.000000    ~d_Lbar         ~u_L                                                            
1466            518   -1   53    0.000000    ~s_Lbar         ~c_L                                                            
1467            519   -1   53    0.000000    ~e_L+           ~nu_eL                                                          
1468            520   -1   53    0.000000    ~mu_L+          ~nu_muL                                                         
1469            521   -1   53    0.000000    ~tau_1+         ~nu_tauL                                                        
1470            522   -1   53    0.000000    ~tau_2+         ~nu_tauL                                                        
1471
1472         38     38    eta_tech0                           0    2    0    350.00000     0.09572     0.95720   0.00000E+00    1
1473            523    1   32    0.442959    b               bbar                                                            
1474            524    1   32    0.000000    t               tbar                                                            
1475            525    1   32    0.557041    g               g                                                               
1476
1477         39     39    LQ_ue           LQ_uebar           -1    1    1    200.00000     0.39162     3.91621   0.00000E+00    1
1478            526    1    0    1.000000    u               e-                                                              
1479
1480         40     40    R0              Rbar0               0    0    1   5000.00000   417.32877  4173.28769   0.00000E+00    1
1481            527    1   32    0.215122    d               sbar                                                            
1482            528    1   32    0.215122    u               cbar                                                            
1483            529    1   32    0.215122    s               bbar                                                            
1484            530    1   32    0.214727    c               tbar                                                            
1485            531   -1   32    0.000000    b               b'bar                                                           
1486            532   -1   32    0.000000    t               t'bar                                                           
1487            533    1    0    0.069953    e-              mu+                                                             
1488            534    1    0    0.069953    mu-             tau+                                                            
1489            535   -1    0    0.000000    tau-            tau'+                                                           
1490
1491         51     51    pi_tech0                            0    0    0    110.00000     0.02899     0.28994   0.00000E+00    1
1492            536    1   32    0.017504    s               sbar                                                            
1493            537    1   32    0.053796    c               cbar                                                            
1494            538    1   32    0.857596    b               bbar                                                            
1495            539    1   32    0.000000    t               tbar                                                            
1496            540    1    0    0.000000    e-              e+                                                              
1497            541    1    0    0.000251    mu-             mu+                                                             
1498            542    1    0    0.070854    tau-            tau+                                                            
1499            543    1   32    0.000000    g               g                                                               
1500
1501         52     52    pi_tech+        pi_tech-            3    0    1    110.00000     0.01070     0.10704   0.00000E+00    1
1502            544    1   32    0.042758    c               sbar                                                            
1503            545    1   32    0.909078    c               bbar                                                            
1504            546    1   32    0.000000    W+              b               bbar                                            
1505            547    1    0    0.000000    e+              nu_e                                                            
1506            548    1    0    0.000170    mu+             nu_mu                                                           
1507            549    1    0    0.047994    tau+            nu_tau                                                          
1508
1509         53     53    pi'_tech0                           0    0    0    110.00000     0.04547     0.45469   0.00000E+00    1
1510            550    1   32    0.011162    s               sbar                                                            
1511            551    1   32    0.034304    c               cbar                                                            
1512            552    1   32    0.546865    b               bbar                                                            
1513            553    1   32    0.000000    t               tbar                                                            
1514            554    1    0    0.000000    e-              e+                                                              
1515            555    1    0    0.000160    mu-             mu+                                                             
1516            556    1    0    0.045181    tau-            tau+                                                            
1517            557    1   32    0.362328    g               g                                                               
1518
1519         54     54    rho_tech0                           0    0    0    210.00000     0.87415     8.74152   0.00000E+00    1
1520            558    1    0    0.144048    W+              W-                                                              
1521            559    1    0    0.352384    W+              pi_tech-                                                        
1522            560    1    0    0.352384    pi_tech+        W-                                                              
1523            561    1    0    0.000000    pi_tech+        pi_tech-                                                        
1524            562    1    0    0.081586    gamma           pi_tech0                                                        
1525            563    1    0    0.029378    gamma           pi'_tech0                                                       
1526            564    1    0    0.001501    Z0              pi_tech0                                                        
1527            565    1    0    0.000721    Z0              pi'_tech0                                                       
1528            566    1   32    0.004490    d               dbar                                                            
1529            567    1   32    0.006482    u               ubar                                                            
1530            568    1   32    0.004490    s               sbar                                                            
1531            569    1   32    0.006482    c               cbar                                                            
1532            570    1   32    0.004485    b               bbar                                                            
1533            571    1   32    0.000000    t               tbar                                                            
1534            572   -1   32    0.000000    b'              b'bar                                                           
1535            573   -1   32    0.000000    t'              t'bar                                                           
1536            574    1    0    0.002889    e-              e+                                                              
1537            575    1    0    0.000967    nu_e            nu_ebar                                                         
1538            576    1    0    0.002889    mu-             mu+                                                             
1539            577    1    0    0.000967    nu_mu           nu_mubar                                                        
1540            578    1    0    0.002889    tau-            tau+                                                            
1541            579    1    0    0.000967    nu_tau          nu_taubar                                                       
1542            580   -1    0    0.000000    tau'-           tau'+                                                           
1543            581   -1    0    0.000000    nu'_tau         nu'_taubar                                                      
1544
1545         55     55    rho_tech+       rho_tech-           3    0    1    210.00000     0.62673     6.26729   0.00000E+00    1
1546            582    1    0    0.143941    W+              Z0                                                              
1547            583    1    0    0.491500    W+              pi_tech0                                                        
1548            584    1    0    0.194259    pi_tech+        Z0                                                              
1549            585    1    0    0.000000    pi_tech+        pi_tech0                                                        
1550            586    1    0    0.113795    pi_tech+        gamma                                                           
1551            587    1    0    0.008460    W+              pi'_tech0                                                       
1552            588    1   32    0.014790    dbar            u                                                               
1553            589    1   32    0.000759    dbar            c                                                               
1554            590    1   32    0.000000    dbar            t                                                               
1555            591   -1   32    0.000000    dbar            t'                                                              
1556            592    1   32    0.000759    sbar            u                                                               
1557            593    1   32    0.014762    sbar            c                                                               
1558            594    1   32    0.000003    sbar            t                                                               
1559            595   -1   32    0.000000    sbar            t'                                                              
1560            596    1   32    0.000000    bbar            u                                                               
1561            597    1   32    0.000027    bbar            c                                                               
1562            598    1   32    0.001934    bbar            t                                                               
1563            599   -1   32    0.000000    bbar            t'                                                              
1564            600   -1   32    0.000000    b'bar           u                                                               
1565            601   -1   32    0.000000    b'bar           c                                                               
1566            602   -1   32    0.000000    b'bar           t                                                               
1567            603   -1   32    0.000000    b'bar           t'                                                              
1568            604    1    0    0.005003    e+              nu_e                                                            
1569            605    1    0    0.005003    mu+             nu_mu                                                           
1570            606    1    0    0.005002    tau+            nu_tau                                                          
1571            607   -1    0    0.000000    tau'+           nu'_tau                                                         
1572
1573         56     56    omega_tech                          0    0    0    210.00000     0.19204     1.92039   0.00000E+00    1
1574            608    1    0    0.133696    gamma           pi_tech0                                                        
1575            609    1    0    0.003283    Z0              pi_tech0                                                        
1576            610    1    0    0.371467    gamma           pi'_tech0                                                       
1577            611    1    0    0.006835    Z0              pi'_tech0                                                       
1578            612    1    0    0.031199    W+              pi_tech-                                                        
1579            613    1    0    0.031199    pi_tech+        W-                                                              
1580            614    1    0    0.001639    W+              W-                                                              
1581            615    1    0    0.000000    pi_tech+        pi_tech-                                                        
1582            616    1   32    0.047205    d               dbar                                                            
1583            617    1   32    0.073708    u               ubar                                                            
1584            618    1   32    0.047205    s               sbar                                                            
1585            619    1   32    0.073705    c               cbar                                                            
1586            620    1   32    0.047161    b               bbar                                                            
1587            621    1   32    0.000000    t               tbar                                                            
1588            622   -1   32    0.000000    b'              b'bar                                                           
1589            623   -1   32    0.000000    t'              t'bar                                                           
1590            624    1    0    0.034740    e-              e+                                                              
1591            625    1    0    0.009160    nu_e            nu_ebar                                                         
1592            626    1    0    0.034740    mu-             mu+                                                             
1593            627    1    0    0.009160    nu_mu           nu_mubar                                                        
1594            628    1    0    0.034738    tau-            tau+                                                            
1595            629    1    0    0.009160    nu_tau          nu_taubar                                                       
1596            630   -1    0    0.000000    tau'-           tau'+                                                           
1597            631   -1    0    0.000000    nu'_tau         nu'_taubar                                                      
1598
1599         61     61    H_L++           H_L--               6    0    1    200.00000     0.88161     8.81606   0.00000E+00    1
1600            632    1    0    0.090264    e+              e+                                                              
1601            633    1    0    0.001805    e+              mu+                                                             
1602            634    1    0    0.001805    e+              tau+                                                            
1603            635    1    0    0.090264    mu+             mu+                                                             
1604            636    1    0    0.001805    mu+             tau+                                                            
1605            637    1    0    0.812250    tau+            tau+                                                            
1606            638    1    0    0.001806    W+              W+                                                              
1607
1608         62     62    H_R++           H_R--               6    0    1    200.00000     0.88001     8.80013   0.00000E+00    1
1609            639    1    0    0.090428    e+              e+                                                              
1610            640    1    0    0.001809    e+              mu+                                                             
1611            641    1    0    0.001808    e+              tau+                                                            
1612            642    1    0    0.090428    mu+             mu+                                                             
1613            643    1    0    0.001808    mu+             tau+                                                            
1614            644    1    0    0.813720    tau+            tau+                                                            
1615            645    1    0    0.000000    W_R+            W_R+                                                            
1616
1617         63     63    W_R+            W_R-                3    0    1    750.00000    19.32815   193.28147   0.00000E+00    1
1618            646    1   32    0.325914    dbar            u                                                               
1619            647    1   32    0.016735    dbar            c                                                               
1620            648    1   32    0.000009    dbar            t                                                               
1621            649    1   32    0.016735    sbar            u                                                               
1622            650    1   32    0.325320    sbar            c                                                               
1623            651    1   32    0.000554    sbar            t                                                               
1624            652    1   32    0.000010    bbar            u                                                               
1625            653    1   32    0.000603    bbar            c                                                               
1626            654    1   32    0.314119    bbar            t                                                               
1627            655    1    0    0.000000    e+              nu_Re                                                           
1628            656    1    0    0.000000    mu+             nu_Rmu                                                          
1629            657    1    0    0.000000    tau+            nu_Rtau                                                         
1630
1631         64     64    nu_Re           nu_Rebar            0    0    1    750.00000     0.00000     0.00000   0.00000E+00    0
1632
1633         65     65    nu_Rmu          nu_Rmubar           0    0    1    750.00000     0.00000     0.00000   0.00000E+00    0
1634
1635         66     66    nu_Rtau         nu_Rtaubar          0    0    1    750.00000     0.00000     0.00000   0.00000E+00    0
1636
1637         81     81    specflav                            0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1638
1639         82     82    rndmflav        rndmflavbar         0    0    1      0.00000     0.00000     0.00000   0.00000E+00    0
1640
1641         83     83    phasespa                            0    0    0      1.00000     0.00000     0.00000   0.00000E+00    1
1642            658    1   12    1.000000    rndmflav        rndmflavbar                                                     
1643
1644         84     84    c-hadron        c-hadronbar         2    0    1      2.00000     0.00000     0.00000   1.00000E-01    1
1645            659    1   42    0.080000    e+              nu_e            s               specflav                        
1646            660    1   42    0.080000    mu+             nu_mu           s               specflav                        
1647            661    1   11    0.760000    u               dbar            s               specflav                        
1648            662    1   11    0.080000    u               sbar            s               specflav                        
1649
1650         85     85    b-hadron        b-hadronbar        -1    0    1      5.00000     0.00000     0.00000   3.87000E-01    1
1651            663    1   42    0.105000    nu_ebar         e-              c               specflav                        
1652            664    1   42    0.105000    nu_mubar        mu-             c               specflav                        
1653            665    1   42    0.040000    nu_taubar       tau-            c               specflav                        
1654            666    1   42    0.500000    ubar            d               c               specflav                        
1655            667    1   42    0.080000    ubar            c               d               specflav                        
1656            668    1   42    0.140000    cbar            s               c               specflav                        
1657            669    1   42    0.010000    cbar            c               s               specflav                        
1658            670    1   42    0.015000    ubar            d               u               specflav                        
1659            671    1   42    0.005000    cbar            s               u               specflav                        
1660
1661         91     91    cluster                             0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1662
1663         92     92    string                              0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1664
1665         93     93    indep.                              0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1666
1667         94     94    CMshower                            0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1668
1669         95     95    SPHEaxis                            0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1670
1671         96     96    THRUaxis                            0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1672
1673         97     97    CLUSjet                             0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1674
1675         98     98    CELLjet                             0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1676
1677         99     99    table                               0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1678
1679        110    101    rho_diff0                           0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1680
1681        111    102    pi0                                 0    0    0      0.13498     0.00000     0.00000   3.00000E-05    0
1682            672    1    0    0.988000    gamma           gamma                                                           
1683            673    1    2    0.012000    gamma           e-              e+                                              
1684
1685        113    103    rho0                                0    0    0      0.76850     0.15100     0.40000   0.00000E+00    1
1686            674    1    3    0.998739    pi+             pi-                                                             
1687            675    1    0    0.000790    pi0             gamma                                                           
1688            676    1    0    0.000380    eta             gamma                                                           
1689            677    1    0    0.000046    mu-             mu+                                                             
1690            678    1    0    0.000045    e-              e+                                                              
1691
1692        115    104    a_20                                0    0    0      1.31800     0.10700     0.25000   0.00000E+00    1
1693            679    1    0    0.347250    rho+            pi-                                                             
1694            680    1    0    0.347250    rho-            pi+                                                             
1695            681    1    0    0.144000    eta             pi0                                                             
1696            682    1    0    0.104000    omega           pi+             pi-                                             
1697            683    1    0    0.024500    K+              K-                                                              
1698            684    1    0    0.012250    K_L0            K_L0                                                            
1699            685    1    0    0.012250    K_S0            K_S0                                                            
1700            686    1    0    0.002800    pi0             gamma                                                           
1701            687    1    0    0.005700    eta'            pi0                                                             
1702
1703        130    105    K_L0                                0    0    0      0.49767     0.00000     0.00000   1.55000E+04    0
1704            688    1    0    0.211200    pi0             pi0             pi0                                             
1705            689    1    0    0.125600    pi+             pi-             pi0                                             
1706            690    1   42    0.193900    nu_ebar         e-              pi+                                             
1707            691    1   42    0.193900    nu_e            e+              pi-                                             
1708            692    1   42    0.135900    nu_mubar        mu-             pi+                                             
1709            693    1   42    0.135900    nu_mu           mu+             pi-                                             
1710            694    1    0    0.002000    pi+             pi-                                                             
1711            695    1    0    0.001000    pi0             pi0                                                             
1712            696    1    0    0.000600    gamma           gamma                                                           
1713
1714        210    106    pi_diffr+       pi_diffr-           3    0    1      0.00000     0.00000     0.00000   0.00000E+00    0
1715
1716        211    107    pi+             pi-                 3    0    1      0.13957     0.00000     0.00000   7.80450E+03    0
1717            697    1    0    0.999877    mu+             nu_mu                                                           
1718            698    1    0    0.000123    e+              nu_e                                                            
1719
1720        213    108    rho+            rho-                3    0    1      0.76690     0.14900     0.40000   0.00000E+00    1
1721            699    1    3    0.999550    pi+             pi0                                                             
1722            700    1    0    0.000450    pi+             gamma                                                           
1723
1724        215    109    a_2+            a_2-                3    0    1      1.31800     0.10700     0.25000   0.00000E+00    1
1725            701    1    0    0.347250    rho+            pi0                                                             
1726            702    1    0    0.347250    rho0            pi+                                                             
1727            703    1    0    0.144000    eta             pi+                                                             
1728            704    1    0    0.104000    omega           pi+             pi0                                             
1729            705    1    0    0.049000    K+              Kbar0                                                           
1730            706    1    0    0.002800    pi+             gamma                                                           
1731            707    1    0    0.005700    eta'            pi+                                                             
1732
1733        220    110    omega_di                            0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1734
1735        221    111    eta                                 0    0    0      0.54745     0.00000     0.00000   0.00000E+00    1
1736            708    1    0    0.392300    gamma           gamma                                                           
1737            709    1    0    0.321000    pi0             pi0             pi0                                             
1738            710    1    0    0.231700    pi+             pi-             pi0                                             
1739            711    1    0    0.047800    gamma           pi+             pi-                                             
1740            712    1    2    0.004900    gamma           e-              e+                                              
1741            713    1    0    0.001300    pi+             pi-             e-              e+                              
1742            714    1    0    0.000300    gamma           mu-             mu+                                             
1743            715    1    0    0.000700    pi0             gamma           gamma                                           
1744
1745        223    112    omega                               0    0    0      0.78194     0.00843     0.10000   0.00000E+00    1
1746            716    1    1    0.890000    pi+             pi-             pi0                                             
1747            717    1    0    0.086930    gamma           pi0                                                             
1748            718    1    3    0.022100    pi+             pi-                                                             
1749            719    1    0    0.000830    eta             gamma                                                           
1750            720    1    0    0.000070    pi0             pi0             gamma                                           
1751            721    1    0    0.000070    e-              e+                                                              
1752
1753        225    113    f_2                                 0    0    0      1.27500     0.18500     0.17000   0.00000E+00    1
1754            722    1    0    0.564000    pi+             pi-                                                             
1755            723    1    0    0.282000    pi0             pi0                                                             
1756            724    1    0    0.072000    pi+             pi-             pi0             pi0                             
1757            725    1    0    0.028000    pi+             pi-             pi+             pi-                             
1758            726    1    0    0.023000    K+              K-                                                              
1759            727    1    0    0.011500    K_L0            K_L0                                                            
1760            728    1    0    0.011500    K_S0            K_S0                                                            
1761            729    1    0    0.005000    eta             eta                                                             
1762            730    1    0    0.003000    pi0             pi0             pi0             pi0                             
1763
1764        310    114    K_S0                                0    0    0      0.49767     0.00000     0.00000   2.67620E+01    1
1765            731    1    0    0.686100    pi+             pi-                                                             
1766            732    1    0    0.313900    pi0             pi0                                                             
1767
1768        311    115    K0              Kbar0               0    0    1      0.49767     0.00000     0.00000   0.00000E+00    1
1769            733    1    0    0.500000    K_L0                                                                            
1770            734    1    0    0.500000    K_S0                                                                            
1771
1772        313    116    K*0             K*bar0              0    0    1      0.89610     0.05050     0.20000   0.00000E+00    1
1773            735    1    3    0.665000    K+              pi-                                                             
1774            736    1    3    0.333000    K0              pi0                                                             
1775            737    1    0    0.002000    K0              gamma                                                           
1776
1777        315    117    K*_20           K*_2bar0            0    0    1      1.43200     0.10900     0.12000   0.00000E+00    1
1778            738    1    0    0.333000    K+              pi-                                                             
1779            739    1    0    0.166000    K0              pi0                                                             
1780            740    1    0    0.168000    K*+             pi-                                                             
1781            741    1    0    0.084000    K*0             pi0                                                             
1782            742    1    0    0.087000    K*+             pi-             pi0                                             
1783            743    1    0    0.043000    K*0             pi+             pi-                                             
1784            744    1    0    0.059000    K+              rho-                                                            
1785            745    1    0    0.029000    K0              rho0                                                            
1786            746    1    0    0.029000    K0              omega                                                           
1787            747    1    0    0.002000    K0              eta                                                             
1788
1789        321    118    K+              K-                  3    0    1      0.49360     0.00000     0.00000   3.70900E+03    0
1790            748    1    0    0.635200    mu+             nu_mu                                                           
1791            749    1    0    0.211600    pi+             pi0                                                             
1792            750    1    0    0.055900    pi+             pi+             pi-                                             
1793            751    1    0    0.017300    pi+             pi0             pi0                                             
1794            752    1   42    0.048200    nu_e            e+              pi0                                             
1795            753    1   42    0.031800    nu_mu           mu+             pi0                                             
1796
1797        323    119    K*+             K*-                 3    0    1      0.89160     0.04980     0.20000   0.00000E+00    1
1798            754    1    3    0.666000    K0              pi+                                                             
1799            755    1    3    0.333000    K+              pi0                                                             
1800            756    1    0    0.001000    K+              gamma                                                           
1801
1802        325    120    K*_2+           K*_2-               3    0    1      1.42500     0.09800     0.12000   0.00000E+00    1
1803            757    1    0    0.332000    K0              pi+                                                             
1804            758    1    0    0.166000    K+              pi0                                                             
1805            759    1    0    0.168000    K*0             pi+                                                             
1806            760    1    0    0.084000    K*+             pi0                                                             
1807            761    1    0    0.086000    K*0             pi+             pi0                                             
1808            762    1    0    0.043000    K*+             pi+             pi-                                             
1809            763    1    0    0.059000    K0              rho+                                                            
1810            764    1    0    0.029000    K+              rho0                                                            
1811            765    1    0    0.029000    K+              omega                                                           
1812            766    1    0    0.002000    K+              eta                                                             
1813            767    1    0    0.002000    K+              gamma                                                           
1814
1815        330    121    phi_diff                            0    0    0      0.00000     0.00000     0.00000   0.00000E+00    0
1816
1817        331    122    eta'                                0    0    0      0.95777     0.00020     0.00200   0.00000E+00    1
1818            768    1    0    0.437000    pi+             pi-             eta                                             
1819            769    1    0    0.208000    pi0             pi0             eta                                             
1820            770    1    0    0.302000    gamma           rho0                                                            
1821            771    1    0    0.030200    gamma           omega                                                           
1822            772    1    0    0.021200    gamma           gamma                                                           
1823            773    1    0    0.001600    pi0             pi0             pi0                                             
1824
1825        333    123    phi                                 0    0    0      1.01940     0.00443     0.01500   0.00000E+00    1
1826            774    1    3    0.489470    K+              K-                                                              
1827            775    1    3    0.340000    K_L0            K_S0                                                            
1828            776    1    0    0.043000    rho-            pi+                                                             
1829            777    1    0    0.043000    rho0            pi0                                                             
1830            778    1    0    0.043000    rho+            pi-                                                             
1831            779    1    1    0.027000    pi+             pi-             pi0                                             
1832            780    1    0    0.012600    gamma           eta                                                             
1833            781    1    0    0.001300    pi0             gamma                                                           
1834            782    1    0    0.000300    e-              e+                                                              
1835            783    1    0    0.000250    mu-             mu+                                                             
1836            784    1    0    0.000080    pi+             pi-                                                             
1837
1838        335    124    f'_2                                0    0    0      1.52500     0.07600     0.20000   0.00000E+00    1
1839            785    1    0    0.444000    K+              K-                                                              
1840            786    1    0    0.222000    K_L0            K_L0                                                            
1841            787    1    0    0.222000    K_S0            K_S0                                                            
1842            788    1    0    0.104000    eta             eta                                                             
1843            789    1    0    0.004000    pi+             pi-                                                             
1844            790    1    0    0.004000    pi0             pi0                                                             
1845
1846        411    125    D+              D-                  3    0    1      1.86930     0.00000     0.00000   3.17000E-01    1
1847            791    0   42    0.070000    e+              nu_e            Kbar0                                           
1848            792    0   42    0.065000    e+              nu_e            K*bar0                                          
1849            793    0   42    0.005000    e+              nu_e            Kbar0           pi0                             
1850            794    0   42    0.005000    e+              nu_e            K-              pi+                             
1851            795    0   42    0.011000    e+              nu_e            K*bar0          pi0                             
1852            796    0   42    0.011000    e+              nu_e            K*-             pi+                             
1853            797    0   42    0.001000    e+              nu_e            pi0                                             
1854            798    0   42    0.001000    e+              nu_e            eta                                             
1855            799    0   42    0.001000    e+              nu_e            eta'                                            
1856
1857            800    0   42    0.001000    e+              nu_e            rho0                                            
1858            801    0   42    0.001000    e+              nu_e            omega                                           
1859            802    1   42    0.070000    mu+             nu_mu           Kbar0                                           
1860            803    1   42    0.065000    mu+             nu_mu           K*bar0                                          
1861            804    1   42    0.005000    mu+             nu_mu           Kbar0           pi0                             
1862            805    1   42    0.005000    mu+             nu_mu           K-              pi+                             
1863            806    1   42    0.011000    mu+             nu_mu           K*bar0          pi0                             
1864            807    1   42    0.011000    mu+             nu_mu           K*-             pi+                             
1865            808    1   42    0.001000    mu+             nu_mu           pi0                                             
1866            809    1   42    0.001000    mu+             nu_mu           eta                                             
1867            810    1   42    0.001000    mu+             nu_mu           eta'                                            
1868            811    1   42    0.001000    mu+             nu_mu           rho0                                            
1869            812    1   42    0.001000    mu+             nu_mu           omega                                           
1870            813    0    0    0.026000    Kbar0           pi+                                                             
1871            814    0    0    0.019000    K*bar0          pi+                                                             
1872            815    0    0    0.066000    Kbar0           rho+                                                            
1873            816    0    0    0.041000    K*bar0          rho+                                                            
1874            817    0    0    0.045000    K*_1bar0        pi+                                                             
1875            818    0    0    0.076000    Kbar0           a_1+                                                            
1876            819    0    0    0.007300    Kbar0           K+                                                              
1877            820    0    0    0.004700    K*bar0          K+                                                              
1878            821    0    0    0.004700    Kbar0           K*+                                                             
1879            822    0    0    0.026000    K*bar0          K*+                                                             
1880            823    0    0    0.001000    pi0             pi+                                                             
1881            824    0    0    0.000600    pi0             rho+                                                            
1882            825    0    0    0.006600    eta             pi+                                                             
1883            826    0    0    0.005000    eta             rho+                                                            
1884            827    0    0    0.003000    eta'            pi+