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