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