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