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