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