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