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