Logo ROOT   6.30.04
Reference Guide
 All Namespaces Files Pages
TAttParticle.cxx
Go to the documentation of this file.
1 // @(#)root/eg:$Id$
2 // Author: Ola Nordmann 29/09/95
3 
4 /*************************************************************************
5  * Copyright (C) 1995-2000, Rene Brun and Fons Rademakers. *
6  * All rights reserved. *
7  * *
8  * For the licensing terms see $ROOTSYS/LICENSE. *
9  * For the list of contributors see $ROOTSYS/README/CREDITS. *
10  *************************************************************************/
11 
12 /** \class TAttParticle
13  \ingroup eg
14 
15 Particle definition, partly based on GEANT3 particle definition
16 */
17 
18 #include "TAttParticle.h"
19 #include "THashList.h"
20 #include "TMath.h"
21 #include "TRandom.h"
22 
23 ClassImp(TAttParticle);
24 
25 THashList *TAttParticle::fgList = new THashList;
26 
27 ////////////////////////////////////////////////////////////////////////////////
28 ///
29 /// Particle definition default constructor
30 ///
31 
32 TAttParticle::TAttParticle()
33 {
34  //do nothing just set some dummy values
35  fPDGMass = 0.0;
36  fPDGStable = kTRUE;
37  fPDGDecayWidth = 0.0;
38  fPDGCharge = 0.0;
39  fParticleType = "";
40  fMCnumberOfPDG = 0;
41  fEnergyCut = 1.e-5;
42  fEnergyLimit = 1.e4;
43  fGranularity = 90;
44 }
45 
46 ////////////////////////////////////////////////////////////////////////////////
47 ///
48 /// Particle definition normal constructor. If the particle is set to be
49 /// stable, the decay width parameter does have no meaning and can be set to
50 /// any value. The parameters granularity, LowerCutOff and HighCutOff are
51 /// used for the construction of the mean free path look up tables. The
52 /// granularity will be the number of logwise energy points for which the
53 /// mean free path will be calculated.
54 ///
55 
56 TAttParticle::TAttParticle(const char *name, const char *title,
57  Double_t Mass, Bool_t Stable,
58  Double_t DecayWidth, Double_t Charge, const char *Type,
59  Int_t MCnumber, Int_t granularity, Double_t LowerCutOff,
60  Double_t HighCutOff) : TNamed(name,title)
61 {
62  fPDGMass = Mass;
63  fPDGStable = Stable;
64  fPDGDecayWidth = DecayWidth;
65  fPDGCharge = Charge;
66  fParticleType = Type;
67  fMCnumberOfPDG = MCnumber;
68  fEnergyCut = LowerCutOff;
69  fEnergyLimit = HighCutOff;
70  fGranularity = granularity;
71 
72  fgList->Add(this);
73 }
74 
75 ////////////////////////////////////////////////////////////////////////////////
76 ///
77 /// Particle destructor
78 ///
79 
80 TAttParticle::~TAttParticle()
81 {
82 }
83 
84 ////////////////////////////////////////////////////////////////////////////////
85 ///
86 /// Converts the ISAJET Particle number into the PDG MC number
87 ///
88 
89 Int_t TAttParticle::ConvertISAtoPDG(Int_t isaNumber)
90 {
91  switch (isaNumber) {
92  case 1 : return 2; // UP .30000E+00 .67
93  case -1 : return -2; // UB .30000E+00 -.67
94  case 2 : return 1; // DN .30000E+00 -.33
95  case -2 : return -1; // DB .30000E+00 .33
96  case 3 : return 3; // ST .50000E+00 -.33
97  case -3 : return -3; // SB .50000E+00 .33
98  case 4 : return 4; // CH .16000E+01 .67
99  case -4 : return -4; // CB .16000E+01 -.67
100  case 5 : return 5; // BT .49000E+01 -.33
101  case -5 : return -5; // BB .49000E+01 .33
102  case 6 : return 7; // TP .30000E+02 .67
103  case -6 : return -7; // TB .30000E+02 -.67
104  case 9 : return 21; // GL 0. 0.00
105  case 80 : return 24; // W+ SIN2W=.23 1.00
106  case -80 : return -24; // W- SIN2W=.23 -1.00
107  case 90 : return 23; // Z0 SIN2W=.23 0.00
108  case 230 : return 311; // K0 .49767E+00 0.00
109  case -230 : return -311; // AK0 .49767E+00 0.00
110  case 330 : return 331; // ETAP .95760E+00 0.00
111  case 340 : return 0; // F- .20300E+01 -1.00
112  case -340 : return 0; // F+ .20300E+01 1.00
113  case 440 : return 441; // ETAC .29760E+01 0.00
114  case 111 : return 113; // RHO0 .77000E+00 0.00
115  case 121 : return 213; // RHO+ .77000E+00 1.00
116  case -121 : return -213; // RHO- .77000E+00 -1.00
117  case 221 : return 223; // OMEG .78260E+00 0.00
118  case 131 : return 323; // K*+ .88810E+00 1.00
119  case -131 : return -323; // K*- .88810E+00 -1.00
120  case 231 : return 313; // K*0 .89220E+00 0.00
121  case -231 : return -313; // AK*0 .89220E+00 0.00
122  case 331 : return 333; // PHI .10196E+01 0.00
123  case 141 : return -423; // AD*0 .20060E+01 0.00
124  case -141 : return 423; // D*0 .20060E+01 0.00
125  case 241 : return -413; // D*- .20086E+01 -1.00
126  case -241 : return 413; // D*+ .20086E+01 1.00
127  case 341 : return 0; // F*- .21400E+01 -1.00
128  case -341 : return 0; // F*+ .21400E+01 1.00
129  case 441 : return 443; // JPSI .30970E+01 0.00
130  case 1140 : return 4222; // SC++ .24300E+01 2.00
131  case -1140 : return -4222; // ASC-- .24300E+01 -2.00
132  case 1240 : return 4212; // SC+ .24300E+01 1.00
133  case -1240 : return -4212; // ASC- .24300E+01 -1.00
134  case 2140 : return 4122; // LC+ .22600E+01 1.00
135  case -2140 : return -4122; // ALC- .22600E+01 -1.00
136  case 2240 : return 4112; // SC0 .24300E+01 0.00
137  case -2240 : return -4112; // ASC0 .24300E+01 0.00
138  case 1340 : return 0; // USC. .25000E+01 1.00
139  case -1340 : return 0; // AUSC. .25000E+01 -1.00
140  case 3140 : return 0; // SUC. .24000E+01 1.00
141  case -3140 : return 0; // ASUC. .24000E+01 -1.00
142  case 2340 : return 0; // DSC. .25000E+01 0.00
143  case -2340 : return 0; // ADSC. .25000E+01 0.00
144  case 3240 : return 0; // SDC. .24000E+01 0.00
145  case -3240 : return 0; // ASDC. .24000E+01 0.00
146  case 3340 : return 0; // SSC. .26000E+01 0.00
147  case -3340 : return 0; // ASSC. .26000E+01 0.00
148  case 1440 : return 0; // UCC. .35500E+01 2.00
149  case -1440 : return 0; // AUCC. .35500E+01 -2.00
150  case 2440 : return 0; // DCC. .35500E+01 1.00
151  case -2440 : return 0; // ADCC. .35500E+01 -1.00
152  case 3440 : return 0; // SCC. .37000E+01 1.00
153  case -3440 : return 0; // ASCC. .37000E+01 -1.00
154  case 1111 : return 2224; // DL++ .12320E+01 2.00
155  case -1111 : return -2224; // ADL-- .12320E+01 -2.00
156  case 1121 : return 2214; // DL+ .12320E+01 1.00
157  case -1121 : return -2214; // ADL- .12320E+01 -1.00
158  case 1221 : return 2114; // DL0 .12320E+01 0.00
159  case -1221 : return -2114; // ADL0 .12320E+01 0.00
160  case 2221 : return 1114; // DL- .12320E+01 -1.00
161  case -2221 : return -1114; // ADL+ .12320E+01 1.00
162  case 1131 : return 3224; // S*+ .13823E+01 1.00
163  case -1131 : return -3224; // AS*- .13823E+01 -1.00
164  case 1231 : return 3214; // S*0 .13820E+01 0.00
165  case -1231 : return -3214; // AS*0 .13820E+01 0.00
166  case 2231 : return 3114; // S*- .13875E+01 -1.00
167  case -2231 : return -3114; // AS*+ .13875E+01 1.00
168  case 1331 : return 3324; // XI*0 .15318E+01 0.00
169  case -1331 : return -3324; // AXI*0 .15318E+01 0.00
170  case 2331 : return 3314; // XI*- .15350E+01 -1.00
171  case -2331 : return -3314; // AXI*+ .15350E+01 1.00
172  case 3331 : return 3334; // OM- .16722E+01 -1.00
173  case -3331 : return -3334; // AOM+ .16722E+01 1.00
174  case 1141 : return 0; // UUC* .26300E+01 2.00
175  case -1141 : return 0; // AUUC* .26300E+01 -2.00
176  case 1241 : return 0; // UDC* .26300E+01 1.00
177  case -1241 : return 0; // AUDC* .26300E+01 -1.00
178  case 2241 : return 0; // DDC* .26300E+01 0.00
179  case -2241 : return 0; // ADDC* .26300E+01 0.00
180  case 1341 : return 0; // USC* .27000E+01 1.00
181  case -1341 : return 0; // AUSC* .27000E+01 -1.00
182  case 2341 : return 0; // DSC* .27000E+01 0.00
183  case -2341 : return 0; // ADSC* .27000E+01 0.00
184  case 3341 : return 0; // SSC* .28000E+01 0.00
185  case -3341 : return 0; // ASSC* .28000E+01 0.00
186  case 1441 : return 0; // UCC* .37500E+01 2.00
187  case -1441 : return 0; // AUCC* .37500E+01 -2.00
188  case 2441 : return 0; // DCC* .37500E+01 1.00
189  case -2441 : return 0; // ADCC* .37500E+01 -1.00
190  case 3441 : return 0; // SCC* .39000E+01 1.00
191  case -3441 : return 0; // ASCC* .39000E+01 -1.00
192  case 4441 : return 0; // CCC* .48000E+01 2.00
193  case -4441 : return 0; // ACCC* .48000E+01 -2.00
194  case 10 : return 22; // Photon
195  case 12 : return 11; // Electron
196  case -12 : return -11; // Positron
197  case 14 : return 13; // Muon-
198  case -14 : return -13; // Muon+
199  case 16 : return 15; // Tau-
200  case -16 : return -15; // Tau+
201  case 11 : return 12; // Neutrino e
202  case -11 : return -12; // Anti Neutrino e
203  case 13 : return 14; // Neutrino Muon
204  case -13 : return -14; // Anti Neutrino Muon
205  case 15 : return 16; // Neutrino Tau
206  case -15 : return -16; // Anti Neutrino Tau
207  case 110 : return 111; // Pion0
208  case 120 : return 211; // Pion+
209  case -120 : return -211; // Pion-
210  case 220 : return 221; // Eta
211  case 130 : return 321; // Kaon+
212  case -130 : return -321; // Kaon-
213  case -20 : return 130; // Kaon Long
214  case 20 : return 310; // Kaon Short
215  case -240 : return -411; // D+
216  case 240 : return 411; // D-
217  case -140 : return 421; // D0
218  case 140 : return -421; // D0 bar
219  case 1120 : return 2212; // Proton
220  case -1120 : return -2212; // Anti Proton
221  case 1220 : return 2112; // Neutron
222  case -1220 : return -2112; // Anti Neutron
223  case 2130 : return 3122; // Lambda
224  case -2130 : return -3122; // Lambda bar
225  case 1130 : return 3222; // Sigma+
226  case -1130 : return -3222; // Sigma bar -
227  case 1230 : return 3212; // Sigma0
228  case -1230 : return -3212; // Sigma bar 0
229  case 2230 : return 3112; // Sigma-
230  case -2230 : return -3112; // Sigma bar +
231  case 1330 : return 3322; // Xi0
232  case -1330 : return -3322; // Xi bar 0
233  case 2330 : return 3312; // Xi-
234  case -2330 : return -3312; // Xi bar +
235  default : return 0; // isajet or pdg number does not exist
236  }
237 }
238 
239 ////////////////////////////////////////////////////////////////////////////////
240 ///
241 /// Defines particles according to the Particle Data Group
242 ///
243 /// For questions regarding distribution or content of the MC particle
244 /// codes, contact
245 /// Gary Wagman (GSWagman@LBL.BITNET, LBL::GSWagman, or GSWagman@LBL.GOV).
246 /// (510)486-6610
247 ///
248 
249 void TAttParticle::DefinePDG()
250 {
251  if (!fgList->IsEmpty()) return;
252 
253  new TAttParticle("down","Q001",0.005,kTRUE, .0,-0.333333333333333,"Quark", 1);
254  new TAttParticle("down bar","Q001",
255  0.005,kTRUE, .0,
256  0.333333333333333,"Quark", -1);
257  new TAttParticle("up","Q002",
258  0.003,kTRUE, .0,
259  0.666666666666666,"Quark", 2);
260  new TAttParticle("up bar","Q002",
261  0.003,kTRUE, .0,
262  -0.666666666666666,"Quark", -2);
263  new TAttParticle("strange","Q003",
264  0.1,kTRUE, .0,
265  -0.333333333333333,"Quark", 3);
266  new TAttParticle("strange bar","Q003",
267  0.1,kTRUE, .0,
268  0.333333333333333,"Quark", -3);
269  new TAttParticle("charm","Q004",
270  1.4,kTRUE, .0,
271  0.666666666666666,"Quark", 4);
272  new TAttParticle("charm bar","Q004",
273  1.4,kTRUE, .0,
274  -0.666666666666666,"Quark", -4);
275  new TAttParticle("bottom","Q005",
276  4.4,kTRUE, .0,
277  -0.333333333333333,"Quark", 5);
278  new TAttParticle("bottom bar","Q005",
279  4.4,kTRUE, .0,
280  0.333333333333333,"Quark", -5);
281  new TAttParticle("top","Q007",
282  173.8,kTRUE, .0,
283  0.666666666666666,"Quark", 6);
284  new TAttParticle("top bar","Q007",
285  173.8,kTRUE, .0,
286  -0.666666666666666,"Quark", -6);
287  new TAttParticle("gluon","G021",
288  .0,kTRUE, .0,
289  0.0,"Gauge Boson", 21);
290 // Entry point of the pdg table conversion
291  new TAttParticle("Searches0","S054",
292  169.0,kTRUE, .0,
293  0.,"Meson", 7);
294  new TAttParticle("e-","S003",
295  5.10999E-04,kTRUE, .0,
296  -1.,"Lepton", 11);
297  new TAttParticle("e+","S003",
298  5.10999E-04,kTRUE, .0,
299  1.,"Lepton", -11);
300  new TAttParticle("nu(e)","S001",
301  .0,kTRUE, .0,
302  0.0,"Lepton", 12);
303  new TAttParticle("nu(e) bar","S001",
304  .0,kTRUE, .0,
305  0.0,"Lepton", -12);
306  new TAttParticle("mu-","S004",
307  .1056583,kFALSE, 2.99591E-19,
308  -1.,"Lepton", 13);
309  new TAttParticle("mu+","S004",
310  .1056583,kFALSE, 2.99591E-19,
311  1.,"Lepton", -13);
312  new TAttParticle("nu(mu)","S002",
313  .0,kTRUE, .0,
314  0.0,"Lepton", 14);
315  new TAttParticle("nu(mu) bar","S002",
316  .0,kTRUE, .0,
317  0.0,"Lepton", -14);
318  new TAttParticle("tau-","S035",
319  1.7771,kFALSE, 2.22700E-12,
320  -1.,"Lepton", 15);
321  new TAttParticle("tau+","S035",
322  1.7771,kFALSE, 2.22700E-12,
323  1.,"Lepton", -15);
324  new TAttParticle("nu(tau)","S036",
325  .0,kTRUE, .0,
326  0.0,"Lepton", 16);
327  new TAttParticle("nu(tau) bar","S036",
328  .0,kTRUE, .0,
329  0.0,"Lepton", -16);
330  new TAttParticle("gamma","S000",
331  .0,kTRUE, .0,
332  0.0,"Gauge Boson", 22);
333  new TAttParticle("Z0","S044",
334  91.18699,kFALSE, 2.49,
335  0.0,"Gauge Boson", 23);
336  new TAttParticle("W+","S043",
337  80.41,kFALSE, 2.06,
338  +1.,"Gauge Boson", 24);
339  new TAttParticle("W-","S043",
340  80.41,kFALSE, 2.06,
341  -1.,"Gauge Boson", -24);
342  new TAttParticle("pi0","S009",
343  .1349764,kFALSE, 7.80000E-09,
344  0.0,"Meson", 111);
345  new TAttParticle("rho(770)0","M009",
346  .7699,kFALSE, .1511999,
347  0.0,"Meson", 113);
348  new TAttParticle("a(2)(1320)0","M012",
349  1.3181,kFALSE, .107,
350  0.0,"Meson", 115);
351  new TAttParticle("rho(3)(1690)0","M015",
352  1.691,kFALSE, .160,
353  0.0,"Meson", 117);
354  new TAttParticle("K(L)0","S013",
355  .4976719,kFALSE, 1.27400E-17,
356  0.0,"Meson", 130);
357  new TAttParticle("pi+","S008",
358  .1395699,kFALSE, 2.52860E-17,
359  1.,"Meson", 211);
360  new TAttParticle("pi-","S008",
361  .1395699,kFALSE, 2.52860E-17,
362  -1.,"Meson", -211);
363  new TAttParticle("rho(770)+","M009",
364  .7699,kFALSE, .1507,
365  1.,"Meson", 213);
366  new TAttParticle("rho(770)-","M009",
367  .7699,kFALSE, .1507,
368  -1.,"Meson", -213);
369  new TAttParticle("a(2)(1320)+","M012",
370  1.3181,kFALSE, .107,
371  1.,"Meson", 215);
372  new TAttParticle("a(2)(1320)-","M012",
373  1.3181,kFALSE, .107,
374  -1.,"Meson", -215);
375  new TAttParticle("rho(3)(1690)+","M015",
376  1.691,kFALSE, .160,
377  1.,"Meson", 217);
378  new TAttParticle("rho(3)(1690)-","M015",
379  1.691,kFALSE, .160,
380  -1.,"Meson", -217);
381  new TAttParticle("eta0","S014",
382  .54730,kFALSE, 1.20000E-06,
383  0.0,"Meson", 221);
384  new TAttParticle("omega(782)0","M001",
385  .78194,kFALSE, 8.43000E-03,
386  0.0,"Meson", 223);
387  new TAttParticle("f(2)(1270)0","M005",
388  1.275,kFALSE, .1855,
389  0.0,"Meson", 225);
390  new TAttParticle("omega(3)(1670)0","M045",
391  1.667,kFALSE, .168,
392  0.0,"Meson", 227);
393  new TAttParticle("f(4)(2050)0","M016",
394  2.044,kFALSE, .208,
395  0.0,"Meson", 229);
396  new TAttParticle("K(S)0","S012",
397  .497672,kFALSE, 7.37400E-15,
398  0.0,"Meson", 310);
399  new TAttParticle("K0","S011",
400  .497672,kFALSE, .0,
401  0.0,"Meson", 311);
402  new TAttParticle("K0 bar","S011",
403  .497672,kFALSE, .0,
404  0.0,"Meson", -311);
405  new TAttParticle("K*(892)0","M018",
406  .89610,kFALSE, 5.05000E-02,
407  0.0,"Meson", 313);
408  new TAttParticle("K*(892)0 bar","M018",
409  .89610,kFALSE, 5.05000E-02,
410  0.0,"Meson", -313);
411  new TAttParticle("K(2)*(1430)0","M022",
412  1.4324,kFALSE, .1089999,
413  0.0,"Meson", 315);
414  new TAttParticle("K(2)*(1430)0 bar","M022",
415  1.4324,kFALSE, .1089999,
416  0.0,"Meson", -315);
417  new TAttParticle("K(3)*(1780)0","M060",
418  1.776,kFALSE, .159,
419  0.0,"Meson", 317);
420  new TAttParticle("K(3)*(1780)0 bar","M060",
421  1.776,kFALSE, .159,
422  0.0,"Meson", -317);
423  new TAttParticle("K(4)*(2045)0","M035",
424  2.045,kFALSE, .198,
425  0.0,"Meson", 319);
426  new TAttParticle("K(4)*(2045)0 bar","M035",
427  2.045,kFALSE, .198,
428  0.0,"Meson", -319);
429  new TAttParticle("K+","S010",
430  .493677,kFALSE, 5.32100E-17,
431  1.,"Meson", 321);
432  new TAttParticle("K-","S010",
433  .493677,kFALSE, 5.32100E-17,
434  -1.,"Meson", -321);
435  new TAttParticle("K*(892)+","M018",
436  .8915899,kFALSE, 5.08000E-02,
437  1.,"Meson", 323);
438  new TAttParticle("K*(892)-","M018",
439  .8915899,kFALSE, 5.08000E-02,
440  -1.,"Meson", -323);
441  new TAttParticle("K(2)*(1430)+","M022",
442  1.4256,kFALSE, 9.85000E-02,
443  1.,"Meson", 325);
444  new TAttParticle("K(2)*(1430)-","M022",
445  1.4256,kFALSE, 9.85000E-02,
446  -1.,"Meson", -325);
447  new TAttParticle("K(3)*(1780)+","M060",
448  1.776,kFALSE, .159,
449  1.,"Meson", 327);
450  new TAttParticle("K(3)*(1780)-","M060",
451  1.776,kFALSE, .159,
452  -1.,"Meson", -327);
453  new TAttParticle("K(4)*(2045)+","M035",
454  2.045,kFALSE, .198,
455  1.,"Meson", 329);
456  new TAttParticle("K(4)*(2045)-","M035",
457  2.045,kFALSE, .198,
458  -1.,"Meson", -329);
459  new TAttParticle("eta'(958)0","M002",
460  .95778,kFALSE, 2.03000E-04,
461  0.0,"Meson", 331);
462  new TAttParticle("phi(1020)0","M004",
463  1.01941,kFALSE, 4.43000E-03,
464  0.0,"Meson", 333);
465  new TAttParticle("f(2)'(1525)0","M013",
466  1.525,kFALSE, 7.60000E-02,
467  0.0,"Meson", 335);
468  new TAttParticle("phi(3)(1850)0","M054",
469  1.854,kFALSE, 8.70000E-02,
470  0.0,"Meson", 337);
471  new TAttParticle("D+","S031",
472  1.8693,kFALSE, 6.23000E-13,
473  1.,"Meson", 411);
474  new TAttParticle("D-","S031",
475  1.8693,kFALSE, 6.23000E-13,
476  -1.,"Meson", -411);
477  new TAttParticle("D*(2010)+","M062",
478  2.01,kTRUE, .0,
479  1.,"Meson", 413);
480  new TAttParticle("D*(2010)-","M062",
481  2.01,kTRUE, .0,
482  -1.,"Meson", -413);
483  new TAttParticle("D(2)*(2460)+","M150",
484  2.4589,kFALSE, 2.30000E-02,
485  1.,"Meson", 415);
486  new TAttParticle("D(2)*(2460)-","M150",
487  2.4589,kFALSE, 2.30000E-02,
488  -1.,"Meson", -415);
489  new TAttParticle("D0","S032",
490  1.8646,kFALSE, 1.58600E-12,
491  0.0,"Meson", 421);
492  new TAttParticle("D*(2007)0","M061",
493  2.0067,kTRUE, .0,
494  0.0,"Meson", 423);
495  new TAttParticle("D(2)*(2460)0","M119",
496  2.4589,kFALSE, 2.30000E-02,
497  0.0,"Meson", 425);
498  new TAttParticle("D(s)+","S034",
499  1.9685,kFALSE, 1.41000E-12,
500  1.,"Meson", 431);
501  new TAttParticle("D(s)-","S034",
502  1.9685,kFALSE, 1.41000E-12,
503  -1.,"Meson", -431);
504  new TAttParticle("D(s)*+","S074",
505  2.1124,kTRUE, .0,
506  1.,"Meson", 433);
507  new TAttParticle("D(s)*-","S074",
508  2.1124,kTRUE, .0,
509  -1.,"Meson", -433);
510  new TAttParticle("eta(c)(1S)0","M026",
511  2.9798,kFALSE, 1.32000E-02,
512  0.0,"Meson", 441);
513  new TAttParticle("J/psi(1S)0","M070",
514  3.09688,kFALSE, 8.70000E-05,
515  0.0,"Meson", 443);
516  new TAttParticle("chi(c2)(1P)0","M057",
517  3.55617,kFALSE, 2.00000E-03,
518  0.0,"Meson", 445);
519  new TAttParticle("B0","S049",
520  5.2792,kFALSE, 4.39000E-13,
521  0.0,"Meson", 511);
522  new TAttParticle("B*0","S085",
523  5.3249,kTRUE, .0,
524  0.0,"Meson", 513);
525  new TAttParticle("B+","S049",
526  5.2789,kFALSE, 4.28000E-13,
527  1.,"Meson", 521);
528  new TAttParticle("B-","S049",
529  5.2789,kFALSE, 4.28000E-13,
530  -1.,"Meson", -521);
531  new TAttParticle("B*+","S085",
532  5.3249,kTRUE, .0,
533  1.,"Meson", 523);
534  new TAttParticle("B*-","S085",
535  5.3249,kTRUE, .0,
536  -1.,"Meson", -523);
537  new TAttParticle("B(s)0","S086",
538  5.3693,kFALSE, 4.90000E-13,
539  0.0,"Meson", 531);
540  new TAttParticle("chi(b0)(1P)0","M076",
541  9.8598,kTRUE, .0,
542  0.0,"Meson", 551);
543  new TAttParticle("Upsilon(1S)0","M049",
544  9.46037,kFALSE, 5.25000E-05,
545  0.0,"Meson", 553);
546  new TAttParticle("chi(b2)(1P)0","M078",
547  9.9132,kTRUE, .0,
548  0.0,"Meson", 555);
549  new TAttParticle("Delta(1620)-","B082",
550  1.62,kFALSE, .15,
551  -1.,"Baryon", 1112);
552  new TAttParticle("Delta(1620)+ bar","B082",
553  1.62,kFALSE, .15,
554  +1.,"Baryon", -1112);
555  new TAttParticle("Delta(1232)-","B033",
556  1.232,kFALSE, .1199999,
557  -1.,"Baryon", 1114);
558  new TAttParticle("Delta(1232)+ bar","B033",
559  1.232,kFALSE, .1199999,
560  +1.,"Baryon", -1114);
561  new TAttParticle("Delta(1905)-","B011",
562  1.905,kFALSE, .3499999,
563  -1.,"Baryon", 1116);
564  new TAttParticle("Delta(1905)+ bar","B011",
565  1.905,kFALSE, .3499999,
566  +1.,"Baryon", -1116);
567  new TAttParticle("Delta(1950)-","B083",
568  1.95,kFALSE, .3,
569  -1.,"Baryon", 1118);
570  new TAttParticle("Delta(1950)+ bar","B083",
571  1.95,kFALSE, .3,
572  +1.,"Baryon", -1118);
573  new TAttParticle("Delta(1620)0","B082",
574  1.62,kFALSE, .15,
575  0.0,"Baryon", 1212);
576  new TAttParticle("Delta(1620)0 bar","B082",
577  1.62,kFALSE, .15,
578  0.0,"Baryon", -1212);
579  new TAttParticle("N(1520)0","B062",
580  1.52,kFALSE, .1199999,
581  0.0,"Baryon", 1214);
582  new TAttParticle("N(1520)0 bar","B062",
583  1.52,kFALSE, .1199999,
584  0.0,"Baryon", -1214);
585  new TAttParticle("Delta(1905)0","B011",
586  1.905,kFALSE, .3499999,
587  0.0,"Baryon", 1216);
588  new TAttParticle("Delta(1905)0 bar","B011",
589  1.905,kFALSE, .3499999,
590  0.0,"Baryon", -1216);
591  new TAttParticle("N(2190)0","B071",
592  2.19,kFALSE, .4499999,
593  0.0,"Baryon", 1218);
594  new TAttParticle("N(2190)0 bar","B071",
595  2.19,kFALSE, .4499999,
596  0.0,"Baryon", -1218);
597  new TAttParticle("n","S017",
598  .9395656,kFALSE, 7.42100E-28,
599  0.0,"Baryon", 2112);
600  new TAttParticle("n bar","S017",
601  .9395656,kFALSE, 7.42100E-28,
602  0.0,"Baryon", -2112);
603  new TAttParticle("Delta(1232)0","B033",
604  1.232,kFALSE, .1199999,
605  0.0,"Baryon", 2114);
606  new TAttParticle("Delta(1232)0 bar","B033",
607  1.232,kFALSE, .1199999,
608  0.0,"Baryon", -2114);
609  new TAttParticle("N(1675)0","B064",
610  1.675,kFALSE, .15,
611  0.0,"Baryon", 2116);
612  new TAttParticle("N(1675)0 bar","B064",
613  1.675,kFALSE, .15,
614  0.0,"Baryon", -2116);
615  new TAttParticle("Delta(1950)0","B083",
616  1.95,kFALSE, .3,
617  0.0,"Baryon", 2118);
618  new TAttParticle("Delta(1950)0 bar","B083",
619  1.95,kFALSE, .3,
620  0.0,"Baryon", -2118);
621  new TAttParticle("Delta(1620)+","B082",
622  1.62,kFALSE, .15,
623  +1.,"Baryon", 2122);
624  new TAttParticle("Delta(1620)- bar","B082",
625  1.62,kFALSE, .15,
626  -1.,"Baryon", -2122);
627  new TAttParticle("N(1520)+","B062",
628  1.52,kFALSE, .1199999,
629  +1.,"Baryon", 2124);
630  new TAttParticle("N(1520)- bar","B062",
631  1.52,kFALSE, .1199999,
632  -1.,"Baryon", -2124);
633  new TAttParticle("Delta(1905)+","B011",
634  1.905,kFALSE, .3499999,
635  +1.,"Baryon", 2126);
636  new TAttParticle("Delta(1905)- bar","B011",
637  1.905,kFALSE, .3499999,
638  -1.,"Baryon", -2126);
639  new TAttParticle("N(2190)+","B071",
640  2.19,kFALSE, .4499999,
641  +1.,"Baryon", 2128);
642  new TAttParticle("N(2190)- bar","B071",
643  2.19,kFALSE, .4499999,
644  -1.,"Baryon", -2128);
645  new TAttParticle("p","S016",
646  .9382722,kTRUE, .0,
647  +1.,"Baryon", 2212);
648  new TAttParticle("p bar","S016",
649  .9382722,kTRUE, .0,
650  -1.,"Baryon", -2212);
651  new TAttParticle("Delta(1232)+","B033",
652  1.232,kFALSE, .1199999,
653  +1.,"Baryon", 2214);
654  new TAttParticle("Delta(1232)- bar","B033",
655  1.232,kFALSE, .1199999,
656  -1.,"Baryon", -2214);
657  new TAttParticle("N(1675)+","B064",
658  1.675,kFALSE, .15,
659  +1.,"Baryon", 2216);
660  new TAttParticle("N(1675)- bar","B064",
661  1.675,kFALSE, .15,
662  -1.,"Baryon", -2216);
663  new TAttParticle("Delta(1950)+","B083",
664  1.95,kFALSE, .3,
665  +1.,"Baryon", 2218);
666  new TAttParticle("Delta(1950)- bar","B083",
667  1.95,kFALSE, .3,
668  -1.,"Baryon", -2218);
669  new TAttParticle("Delta(1620)++","B082",
670  1.62,kFALSE, .15,
671  +2.,"Baryon", 2222);
672  new TAttParticle("Delta(1620)-- bar","B082",
673  1.62,kFALSE, .15,
674  -2.,"Baryon", -2222);
675  new TAttParticle("Delta(1232)++","B033",
676  1.232,kFALSE, .1199999,
677  +2.,"Baryon", 2224);
678  new TAttParticle("Delta(1232)-- bar","B033",
679  1.232,kFALSE, .1199999,
680  -2.,"Baryon", -2224);
681  new TAttParticle("Delta(1905)++","B011",
682  1.905,kFALSE, .3499999,
683  +2.,"Baryon", 2226);
684  new TAttParticle("Delta(1905)-- bar","B011",
685  1.905,kFALSE, .3499999,
686  -2.,"Baryon", -2226);
687  new TAttParticle("Delta(1950)++","B083",
688  1.95,kFALSE, .3,
689  +2.,"Baryon", 2228);
690  new TAttParticle("Delta(1950)-- bar","B083",
691  1.95,kFALSE, .3,
692  -2.,"Baryon", -2228);
693  new TAttParticle("Sigma-","S020",
694  1.19744,kFALSE, 4.45000E-15,
695  -1.,"Baryon", 3112);
696  new TAttParticle("Sigma+ bar","S020",
697  1.19744,kFALSE, 4.45000E-15,
698  +1.,"Baryon", -3112);
699  new TAttParticle("Sigma(1385)-","B043",
700  1.3872,kFALSE, 3.94000E-02,
701  -1.,"Baryon", 3114);
702  new TAttParticle("Sigma(1385)+ bar","B043",
703  1.3872,kFALSE, 3.94000E-02,
704  +1.,"Baryon", -3114);
705  new TAttParticle("Sigma(1775)-","B045",
706  1.775,kFALSE, .1199999,
707  -1.,"Baryon", 3116);
708  new TAttParticle("Sigma(1775)+ bar","B045",
709  1.775,kFALSE, .1199999,
710  +1.,"Baryon", -3116);
711  new TAttParticle("Sigma(2030)-","B047",
712  2.03,kFALSE, .18,
713  -1.,"Baryon", 3118);
714  new TAttParticle("Sigma(2030)+ bar","B047",
715  2.03,kFALSE, .18,
716  +1.,"Baryon", -3118);
717  new TAttParticle("Lambda0","S018",
718  1.11568,kFALSE, 2.50100E-15,
719  0.0,"Baryon", 3122);
720  new TAttParticle("Lambda0 bar","S018",
721  1.11568,kFALSE, 2.50100E-15,
722  0.0,"Baryon", -3122);
723  new TAttParticle("Lambda(1520)0","B038",
724  1.5195,kFALSE, 1.56000E-02,
725  0.0,"Baryon", 3124);
726  new TAttParticle("Lambda(1520)0 bar","B038",
727  1.5195,kFALSE, 1.56000E-02,
728  0.0,"Baryon", -3124);
729  new TAttParticle("Lambda(1820)0","B039",
730  1.82,kFALSE, 8.00000E-02,
731  0.0,"Baryon", 3126);
732  new TAttParticle("Lambda(1820)0 bar","B039",
733  1.82,kFALSE, 8.00000E-02,
734  0.0,"Baryon", -3126);
735  new TAttParticle("Lambda(2100)0","B041",
736  2.1,kFALSE, .2,
737  0.0,"Baryon", 3128);
738  new TAttParticle("Lambda(2100)0 bar","B041",
739  2.1,kFALSE, .2,
740  0.0,"Baryon", -3128);
741  new TAttParticle("Sigma0","S021",
742  1.19255,kFALSE, 8.90000E-06,
743  0.0,"Baryon", 3212);
744  new TAttParticle("Sigma0 bar","S021",
745  1.19255,kFALSE, 8.90000E-06,
746  0.0,"Baryon", -3212);
747  new TAttParticle("Sigma(1385)0","B043",
748  1.3837,kFALSE, 3.60000E-02,
749  0.0,"Baryon", 3214);
750  new TAttParticle("Sigma(1385)0 bar","B043",
751  1.3837,kFALSE, 3.60000E-02,
752  0.0,"Baryon", -3214);
753  new TAttParticle("Sigma(1775)0","B045",
754  1.775,kFALSE, .1199999,
755  0.0,"Baryon", 3216);
756  new TAttParticle("Sigma(1775)0 bar","B045",
757  1.775,kFALSE, .1199999,
758  0.0,"Baryon", -3216);
759  new TAttParticle("Sigma(2030)0","B047",
760  2.03,kFALSE, .18,
761  0.0,"Baryon", 3218);
762  new TAttParticle("Sigma(2030)0 bar","B047",
763  2.03,kFALSE, .18,
764  0.0,"Baryon", -3218);
765  new TAttParticle("Sigma+","S019",
766  1.18937,kFALSE, 8.24000E-15,
767  +1.,"Baryon", 3222);
768  new TAttParticle("Sigma- bar","S019",
769  1.18937,kFALSE, 8.24000E-15,
770  -1.,"Baryon", -3222);
771  new TAttParticle("Sigma(1385)+","B043",
772  1.3828,kFALSE, 3.58000E-02,
773  +1.,"Baryon", 3224);
774  new TAttParticle("Sigma(1385)- bar","B043",
775  1.3828,kFALSE, 3.58000E-02,
776  -1.,"Baryon", -3224);
777  new TAttParticle("Sigma(1775)+","B045",
778  1.775,kFALSE, .1199999,
779  +1.,"Baryon", 3226);
780  new TAttParticle("Sigma(1775)- bar","B045",
781  1.775,kFALSE, .1199999,
782  -1.,"Baryon", -3226);
783  new TAttParticle("Sigma(2030)+","B047",
784  2.03,kFALSE, .18,
785  +1.,"Baryon", 3228);
786  new TAttParticle("Sigma(2030)- bar","B047",
787  2.03,kFALSE, .18,
788  -1.,"Baryon", -3228);
789  new TAttParticle("Xi-","S022",
790  1.32132,kFALSE, 4.02000E-15,
791  -1.,"Baryon", 3312);
792  new TAttParticle("Xi+ bar","S022",
793  1.32132,kFALSE, 4.02000E-15,
794  +1.,"Baryon", -3312);
795  new TAttParticle("Xi(1530)-","B049",
796  1.535,kFALSE, 9.90000E-03,
797  -1.,"Baryon", 3314);
798  new TAttParticle("Xi(1530)+ bar","B049",
799  1.535,kFALSE, 9.90000E-03,
800  +1.,"Baryon", -3314);
801  new TAttParticle("Xi0","S023",
802  1.3149,kFALSE, 2.27000E-15,
803  0.0,"Baryon", 3322);
804  new TAttParticle("Xi0 bar","S023",
805  1.3149,kFALSE, 2.27000E-15,
806  0.0,"Baryon", -3322);
807  new TAttParticle("Xi(1530)0","B049",
808  1.5318,kFALSE, 9.10000E-03,
809  0.0,"Baryon", 3324);
810  new TAttParticle("Xi(1530)0 bar","B049",
811  1.5318,kFALSE, 9.10000E-03,
812  0.0,"Baryon", -3324);
813  new TAttParticle("Omega-","S024",
814  1.67245,kFALSE, 8.01000E-15,
815  -1.,"Baryon", 3334);
816  new TAttParticle("Omega+ bar","S024",
817  1.67245,kFALSE, 8.01000E-15,
818  +1.,"Baryon", -3334);
819  new TAttParticle("Sigma(c)(2455)0","B104",
820  2.4524,kTRUE, .0,
821  0.0,"Baryon", 4112);
822  new TAttParticle("Sigma(c)(2455)0 bar","B104",
823  2.4524,kTRUE, .0,
824  0.0,"Baryon", -4112);
825  new TAttParticle("Lambda(c)+","S033",
826  2.2849,kFALSE, 3.29000E-12,
827  +1.,"Baryon", 4122);
828  new TAttParticle("Lambda(c)- bar","S033",
829  2.2849,kFALSE, 3.29000E-12,
830  -1.,"Baryon", -4122);
831  new TAttParticle("Sigma(c)(2455)+","B104",
832  2.4538,kTRUE, .0,
833  +1.,"Baryon", 4212);
834  new TAttParticle("Sigma(c)(2455)- bar","B104",
835  2.4538,kTRUE, .0,
836  -1.,"Baryon", -4212);
837  new TAttParticle("Sigma(c)(2455)++","B104",
838  2.4531,kTRUE, .0,
839  +2.,"Baryon", 4222);
840  new TAttParticle("Sigma(c)(2455)-- bar","B104",
841  2.4531,kTRUE, .0,
842  -2.,"Baryon", -4222);
843  new TAttParticle("Xi(c)0","S048",
844  2.4703,kFALSE, 6.70000E-12,
845  0.0,"Baryon", 4312);
846  new TAttParticle("Xi(c)0 bar","S048",
847  2.4703,kFALSE, 6.70000E-12,
848  0.0,"Baryon", -4312);
849  new TAttParticle("Xi(c)+","S045",
850  2.4651,kFALSE, 1.86000E-12,
851  +1.,"Baryon", 4322);
852  new TAttParticle("Xi(c)- bar","S045",
853  2.4651,kFALSE, 1.86000E-12,
854  -1.,"Baryon", -4322);
855  new TAttParticle("Lambda(b)0","S040",
856  5.64,kFALSE, 6.20000E-13,
857  0.0,"Baryon", 5122);
858  new TAttParticle("Lambda(b)0 bar","S040",
859  5.64,kFALSE, 6.20000E-13,
860  0.0,"Baryon", -5122);
861  new TAttParticle("a(0)(980)0","M036",
862  .9824,kTRUE, .0,
863  0.0,"Meson", 10111);
864  new TAttParticle("b(1)(1235)0","M011",
865  1.231,kFALSE, .142,
866  0.0,"Meson", 10113);
867  new TAttParticle("pi(2)(1670)0","M034",
868  1.67,kFALSE, .2399999,
869  0.0,"Meson", 10115);
870  new TAttParticle("a(0)(980)+","M036",
871  .9834,kTRUE, .0,
872  1.,"Meson", 10211);
873  new TAttParticle("a(0)(980)-","M036",
874  .9834,kTRUE, .0,
875  -1.,"Meson", -10211);
876  new TAttParticle("b(1)(1235)+","M011",
877  1.2295,kFALSE, .142,
878  1.,"Meson", 10213);
879  new TAttParticle("b(1)(1235)-","M011",
880  1.2295,kFALSE, .142,
881  -1.,"Meson", -10213);
882  new TAttParticle("pi(2)(1670)+","M034",
883  1.67,kFALSE, .2399999,
884  1.,"Meson", 10215);
885  new TAttParticle("pi(2)(1670)-","M034",
886  1.67,kFALSE, .2399999,
887  -1.,"Meson", -10215);
888  new TAttParticle("f(0)(980)0","M003",
889  .98,kTRUE, .0,
890  0.0,"Meson", 10221);
891  new TAttParticle("h(1)(1170)0","M030",
892  1.17,kFALSE, .36,
893  0.0,"Meson", 10223);
894  new TAttParticle("K(0)*(1430)0","M019",
895  1.429,kFALSE, .287,
896  0.0,"Meson", 10311);
897  new TAttParticle("K(0)*(1430)0 bar","M019",
898  1.429,kFALSE, .287,
899  0.0,"Meson", -10311);
900  new TAttParticle("K(1)(1270)0","M028",
901  1.272,kFALSE, 9.00000E-02,
902  0.0,"Meson", 10313);
903  new TAttParticle("K(1)(1270)0 bar","M028",
904  1.272,kFALSE, 9.00000E-02,
905  0.0,"Meson", -10313);
906  new TAttParticle("K(2)(1770)0","M023",
907  1.773,kFALSE, .186,
908  0.0,"Meson", 10315);
909  new TAttParticle("K(2)(1770)0 bar","M023",
910  1.773,kFALSE, .186,
911  0.0,"Meson", -10315);
912  new TAttParticle("K(0)*(1430)+","M019",
913  1.429,kFALSE, .287,
914  1.,"Meson", 10321);
915  new TAttParticle("K(0)*(1430)-","M019",
916  1.429,kFALSE, .287,
917  -1.,"Meson", -10321);
918  new TAttParticle("K(1)(1270)+","M028",
919  1.272,kFALSE, 9.00000E-02,
920  1.,"Meson", 10323);
921  new TAttParticle("K(1)(1270)-","M028",
922  1.272,kFALSE, 9.00000E-02,
923  -1.,"Meson", -10323);
924  new TAttParticle("K(2)(1770)+","M023",
925  1.773,kFALSE, .186,
926  1.,"Meson", 10325);
927  new TAttParticle("K(2)(1770)-","M023",
928  1.773,kFALSE, .186,
929  -1.,"Meson", -10325);
930  new TAttParticle("phi(1680)0","M067",
931  1.68,kFALSE, .15,
932  0.0,"Meson", 10333);
933  new TAttParticle("D(1)(2420)0","M097",
934  2.4228,kFALSE, 1.80000E-02,
935  0.0,"Meson", 10423);
936  new TAttParticle("D(s1)(2536)+","M121",
937  2.53535,kTRUE, .0,
938  1.,"Meson", 10433);
939  new TAttParticle("D(s1)(2536)-","M121",
940  2.53535,kTRUE, .0,
941  -1.,"Meson", -10433);
942  new TAttParticle("chi(c0)(1P)0","M056",
943  3.4151,kFALSE, 1.40000E-02,
944  0.0,"Meson", 10441);
945  new TAttParticle("chi(c1)(1P)0","M055",
946  3.51053,kFALSE, 8.80000E-04,
947  0.0,"Meson", 10443);
948  new TAttParticle("chi(b0)(2P)0","M079",
949  10.23209,kTRUE, .0,
950  0.0,"Meson", 10551);
951  new TAttParticle("chi(b1)(1P)0","M077",
952  9.8919,kTRUE, .0,
953  0.0,"Meson", 10553);
954  new TAttParticle("chi(b2)(2P)0","M081",
955  10.2685,kTRUE, .0,
956  0.0,"Meson", 10555);
957  new TAttParticle("Delta(1900)-","B030",
958  1.9,kFALSE, .2,
959  -1.,"Baryon", 11112);
960  new TAttParticle("Delta(1900)+ bar","B030",
961  1.9,kFALSE, .2,
962  +1.,"Baryon", -11112);
963  new TAttParticle("Delta(1700)-","B010",
964  1.7,kFALSE, .3,
965  -1.,"Baryon", 11114);
966  new TAttParticle("Delta(1700)+ bar","B010",
967  1.7,kFALSE, .3,
968  +1.,"Baryon", -11114);
969  new TAttParticle("Delta(1930)-","B013",
970  1.93,kFALSE, .3499999,
971  -1.,"Baryon", 11116);
972  new TAttParticle("Delta(1930)+ bar","B013",
973  1.93,kFALSE, .3499999,
974  +1.,"Baryon", -11116);
975  new TAttParticle("Delta(1900)0","B030",
976  1.9,kFALSE, .2,
977  0.0,"Baryon", 11212);
978  new TAttParticle("Delta(1900)0 bar","B030",
979  1.9,kFALSE, .2,
980  0.0,"Baryon", -11212);
981  new TAttParticle("Delta(1930)0","B013",
982  1.93,kFALSE, .3499999,
983  0.0,"Baryon", 11216);
984  new TAttParticle("Delta(1930)0 bar","B013",
985  1.93,kFALSE, .3499999,
986  0.0,"Baryon", -11216);
987  new TAttParticle("N(1440)0","B061",
988  1.44,kFALSE, .3499999,
989  0.0,"Baryon", 12112);
990  new TAttParticle("N(1440)0 bar","B061",
991  1.44,kFALSE, .3499999,
992  0.0,"Baryon", -12112);
993  new TAttParticle("Delta(1700)0","B010",
994  1.7,kFALSE, .3,
995  0.0,"Baryon", 12114);
996  new TAttParticle("Delta(1700)0 bar","B010",
997  1.7,kFALSE, .3,
998  0.0,"Baryon", -12114);
999  new TAttParticle("N(1680)0","B065",
1000  1.68,kFALSE, .1299999,
1001  0.0,"Baryon", 12116);
1002  new TAttParticle("N(1680)0 bar","B065",
1003  1.68,kFALSE, .1299999,
1004  0.0,"Baryon", -12116);
1005  new TAttParticle("Delta(1900)+","B030",
1006  1.9,kFALSE, .2,
1007  +1.,"Baryon", 12122);
1008  new TAttParticle("Delta(1900)- bar","B030",
1009  1.9,kFALSE, .2,
1010  -1.,"Baryon", -12122);
1011  new TAttParticle("Delta(1930)+","B013",
1012  1.93,kFALSE, .3499999,
1013  +1.,"Baryon", 12126);
1014  new TAttParticle("Delta(1930)- bar","B013",
1015  1.93,kFALSE, .3499999,
1016  -1.,"Baryon", -12126);
1017  new TAttParticle("N(1440)+","B061",
1018  1.44,kFALSE, .3499999,
1019  +1.,"Baryon", 12212);
1020  new TAttParticle("N(1440)- bar","B061",
1021  1.44,kFALSE, .3499999,
1022  -1.,"Baryon", -12212);
1023  new TAttParticle("Delta(1700)+","B010",
1024  1.7,kFALSE, .3,
1025  +1.,"Baryon", 12214);
1026  new TAttParticle("Delta(1700)- bar","B010",
1027  1.7,kFALSE, .3,
1028  -1.,"Baryon", -12214);
1029  new TAttParticle("N(1680)+","B065",
1030  1.68,kFALSE, .1299999,
1031  +1.,"Baryon", 12216);
1032  new TAttParticle("N(1680)- bar","B065",
1033  1.68,kFALSE, .1299999,
1034  -1.,"Baryon", -12216);
1035  new TAttParticle("Delta(1900)++","B030",
1036  1.9,kFALSE, .2,
1037  +2.,"Baryon", 12222);
1038  new TAttParticle("Delta(1900)-- bar","B030",
1039  1.9,kFALSE, .2,
1040  -2.,"Baryon", -12222);
1041  new TAttParticle("Delta(1700)++","B010",
1042  1.7,kFALSE, .3,
1043  +2.,"Baryon", 12224);
1044  new TAttParticle("Delta(1700)-- bar","B010",
1045  1.7,kFALSE, .3,
1046  -2.,"Baryon", -12224);
1047  new TAttParticle("Delta(1930)++","B013",
1048  1.93,kFALSE, .3499999,
1049  +2.,"Baryon", 12226);
1050  new TAttParticle("Delta(1930)-- bar","B013",
1051  1.93,kFALSE, .3499999,
1052  -2.,"Baryon", -12226);
1053  new TAttParticle("Sigma(1660)-","B079",
1054  1.66,kFALSE, .1,
1055  -1.,"Baryon", 13112);
1056  new TAttParticle("Sigma(1660)+ bar","B079",
1057  1.66,kFALSE, .1,
1058  +1.,"Baryon", -13112);
1059  new TAttParticle("Sigma(1670)-","B051",
1060  1.67,kFALSE, 6.00000E-02,
1061  -1.,"Baryon", 13114);
1062  new TAttParticle("Sigma(1670)+ bar","B051",
1063  1.67,kFALSE, 6.00000E-02,
1064  +1.,"Baryon", -13114);
1065  new TAttParticle("Sigma(1915)-","B046",
1066  1.915,kFALSE, .1199999,
1067  -1.,"Baryon", 13116);
1068  new TAttParticle("Sigma(1915)+ bar","B046",
1069  1.915,kFALSE, .1199999,
1070  +1.,"Baryon", -13116);
1071  new TAttParticle("Lambda(1405)0","B037",
1072  1.407,kFALSE, 5.00000E-02,
1073  0.0,"Baryon", 13122);
1074  new TAttParticle("Lambda(1405)0 bar","B037",
1075  1.407,kFALSE, 5.00000E-02,
1076  0.0,"Baryon", -13122);
1077  new TAttParticle("Lambda(1690)0","B055",
1078  1.69,kFALSE, 6.00000E-02,
1079  0.0,"Baryon", 13124);
1080  new TAttParticle("Lambda(1690)0 bar","B055",
1081  1.69,kFALSE, 6.00000E-02,
1082  0.0,"Baryon", -13124);
1083  new TAttParticle("Lambda(1830)0","B056",
1084  1.83,kFALSE, 9.50000E-02,
1085  0.0,"Baryon", 13126);
1086  new TAttParticle("Lambda(1830)0 bar","B056",
1087  1.83,kFALSE, 9.50000E-02,
1088  0.0,"Baryon", -13126);
1089  new TAttParticle("Sigma(1660)0","B079",
1090  1.66,kFALSE, .1,
1091  0.0,"Baryon", 13212);
1092  new TAttParticle("Sigma(1660)0 bar","B079",
1093  1.66,kFALSE, .1,
1094  0.0,"Baryon", -13212);
1095  new TAttParticle("Sigma(1670)0","B051",
1096  1.67,kFALSE, 6.00000E-02,
1097  0.0,"Baryon", 13214);
1098  new TAttParticle("Sigma(1670)0 bar","B051",
1099  1.67,kFALSE, 6.00000E-02,
1100  0.0,"Baryon", -13214);
1101  new TAttParticle("Sigma(1915)0","B046",
1102  1.915,kFALSE, .1199999,
1103  0.0,"Baryon", 13216);
1104  new TAttParticle("Sigma(1915)0 bar","B046",
1105  1.915,kFALSE, .1199999,
1106  0.0,"Baryon", -13216);
1107  new TAttParticle("Sigma(1660)+","B079",
1108  1.66,kFALSE, .1,
1109  +1.,"Baryon", 13222);
1110  new TAttParticle("Sigma(1660)- bar","B079",
1111  1.66,kFALSE, .1,
1112  -1.,"Baryon", -13222);
1113  new TAttParticle("Sigma(1670)+","B051",
1114  1.67,kFALSE, 6.00000E-02,
1115  +1.,"Baryon", 13224);
1116  new TAttParticle("Sigma(1670)- bar","B051",
1117  1.67,kFALSE, 6.00000E-02,
1118  -1.,"Baryon", -13224);
1119  new TAttParticle("Sigma(1915)+","B046",
1120  1.915,kFALSE, .1199999,
1121  +1.,"Baryon", 13226);
1122  new TAttParticle("Sigma(1915)- bar","B046",
1123  1.915,kFALSE, .1199999,
1124  -1.,"Baryon", -13226);
1125  new TAttParticle("Xi(1820)-","B050",
1126  1.823,kFALSE, 2.40000E-02,
1127  -1.,"Baryon", 13314);
1128  new TAttParticle("Xi(1820)+ bar","B050",
1129  1.823,kFALSE, 2.40000E-02,
1130  +1.,"Baryon", -13314);
1131  new TAttParticle("Xi(1820)0","B050",
1132  1.823,kFALSE, 2.40000E-02,
1133  0.0,"Baryon", 13324);
1134  new TAttParticle("Xi(1820)0 bar","B050",
1135  1.823,kFALSE, 2.40000E-02,
1136  0.0,"Baryon", -13324);
1137  new TAttParticle("pi(1300)0","M058",
1138  1.3,kTRUE, .4,
1139  0.0,"Meson", 20111);
1140  new TAttParticle("a(1)(1260)0","M010",
1141  1.23,kTRUE, .4,
1142  0.0,"Meson", 20113);
1143  new TAttParticle("pi(1300)+","M058",
1144  1.3,kTRUE, .4,
1145  1.,"Meson", 20211);
1146  new TAttParticle("pi(1300)-","M058",
1147  1.3,kTRUE, .4,
1148  -1.,"Meson", -20211);
1149  new TAttParticle("a(1)(1260)+","M010",
1150  1.23,kTRUE, .4,
1151  1.,"Meson", 20213);
1152  new TAttParticle("a(1)(1260)-","M010",
1153  1.23,kTRUE, .4,
1154  -1.,"Meson", -20213);
1155  new TAttParticle("eta(1295)0","M037",
1156  1.297,kFALSE, 5.30000E-02,
1157  0.0,"Meson", 20221);
1158  new TAttParticle("f(1)(1285)0","M008",
1159  1.282,kFALSE, 2.40000E-02,
1160  0.0,"Meson", 20223);
1161  new TAttParticle("f(2)(2010)0","M106",
1162  2.01,kFALSE, .2,
1163  0.0,"Meson", 20225);
1164  new TAttParticle("K(1)(1400)0","M064",
1165  1.402,kFALSE, .1739999,
1166  0.0,"Meson", 20313);
1167  new TAttParticle("K(1)(1400)0 bar","M064",
1168  1.402,kFALSE, .1739999,
1169  0.0,"Meson", -20313);
1170  new TAttParticle("K(2)(1820)0","M146",
1171  1.816,kFALSE, .2759999,
1172  0.0,"Meson", 20315);
1173  new TAttParticle("K(2)(1820)0 bar","M146",
1174  1.816,kFALSE, .2759999,
1175  0.0,"Meson", -20315);
1176  new TAttParticle("K(1)(1400)+","M064",
1177  1.402,kFALSE, .1739999,
1178  1.,"Meson", 20323);
1179  new TAttParticle("K(1)(1400)-","M064",
1180  1.402,kFALSE, .1739999,
1181  -1.,"Meson", -20323);
1182  new TAttParticle("K(2)(1820)+","M146",
1183  1.816,kFALSE, .2759999,
1184  1.,"Meson", 20325);
1185  new TAttParticle("K(2)(1820)-","M146",
1186  1.816,kFALSE, .2759999,
1187  -1.,"Meson", -20325);
1188  new TAttParticle("psi(2S)0","M071",
1189  3.686,kFALSE, 2.77000E-04,
1190  0.0,"Meson", 20443);
1191  new TAttParticle("Upsilon(2S)0","M052",
1192  10.0233,kFALSE, 4.40000E-05,
1193  0.0,"Meson", 20553);
1194  new TAttParticle("Delta(1910)-","B012",
1195  1.91,kFALSE, .25,
1196  -1.,"Baryon", 21112);
1197  new TAttParticle("Delta(1910)+ bar","B012",
1198  1.91,kFALSE, .25,
1199  +1.,"Baryon", -21112);
1200  new TAttParticle("Delta(1920)-","B117",
1201  1.92,kFALSE, .2,
1202  -1.,"Baryon", 21114);
1203  new TAttParticle("Delta(1920)+ bar","B117",
1204  1.92,kFALSE, .2,
1205  +1.,"Baryon", -21114);
1206  new TAttParticle("Delta(1910)0","B012",
1207  1.91,kFALSE, .25,
1208  0.0,"Baryon", 21212);
1209  new TAttParticle("Delta(1910)0 bar","B012",
1210  1.91,kFALSE, .25,
1211  0.0,"Baryon", -21212);
1212  new TAttParticle("N(1700)0","B018",
1213  1.7,kFALSE, .1,
1214  0.0,"Baryon", 21214);
1215  new TAttParticle("N(1700)0 bar","B018",
1216  1.7,kFALSE, .1,
1217  0.0,"Baryon", -21214);
1218  new TAttParticle("N(1535)0","B063",
1219  1.535,kFALSE, .15,
1220  0.0,"Baryon", 22112);
1221  new TAttParticle("N(1535)0 bar","B063",
1222  1.535,kFALSE, .15,
1223  0.0,"Baryon", -22112);
1224  new TAttParticle("Delta(1920)0","B117",
1225  1.92,kFALSE, .2,
1226  0.0,"Baryon", 22114);
1227  new TAttParticle("Delta(1920)0 bar","B117",
1228  1.92,kFALSE, .2,
1229  0.0,"Baryon", -22114);
1230  new TAttParticle("Delta(1910)+","B012",
1231  1.91,kFALSE, .25,
1232  +1.,"Baryon", 22122);
1233  new TAttParticle("Delta(1910)- bar","B012",
1234  1.91,kFALSE, .25,
1235  -1.,"Baryon", -22122);
1236  new TAttParticle("N(1700)+","B018",
1237  1.7,kFALSE, .1,
1238  +1.,"Baryon", 22124);
1239  new TAttParticle("N(1700)- bar","B018",
1240  1.7,kFALSE, .1,
1241  -1.,"Baryon", -22124);
1242  new TAttParticle("N(1535)+","B063",
1243  1.535,kFALSE, .15,
1244  +1.,"Baryon", 22212);
1245  new TAttParticle("N(1535)- bar","B063",
1246  1.535,kFALSE, .15,
1247  -1.,"Baryon", -22212);
1248  new TAttParticle("Delta(1920)+","B117",
1249  1.92,kFALSE, .2,
1250  +1.,"Baryon", 22214);
1251  new TAttParticle("Delta(1920)- bar","B117",
1252  1.92,kFALSE, .2,
1253  -1.,"Baryon", -22214);
1254  new TAttParticle("Delta(1910)++","B012",
1255  1.91,kFALSE, .25,
1256  +2.,"Baryon", 22222);
1257  new TAttParticle("Delta(1910)-- bar","B012",
1258  1.91,kFALSE, .25,
1259  -2.,"Baryon", -22222);
1260  new TAttParticle("Delta(1920)++","B117",
1261  1.92,kFALSE, .2,
1262  +2.,"Baryon", 22224);
1263  new TAttParticle("Delta(1920)-- bar","B117",
1264  1.92,kFALSE, .2,
1265  -2.,"Baryon", -22224);
1266  new TAttParticle("Sigma(1750)-","B057",
1267  1.75,kFALSE, 9.00000E-02,
1268  -1.,"Baryon", 23112);
1269  new TAttParticle("Sigma(1750)+ bar","B057",
1270  1.75,kFALSE, 9.00000E-02,
1271  +1.,"Baryon", -23112);
1272  new TAttParticle("Sigma(1940)-","B098",
1273  1.94,kFALSE, .2199999,
1274  -1.,"Baryon", 23114);
1275  new TAttParticle("Sigma(1940)+ bar","B098",
1276  1.94,kFALSE, .2199999,
1277  +1.,"Baryon", -23114);
1278  new TAttParticle("Lambda(1600)0","B101",
1279  1.6,kFALSE, .15,
1280  0.0,"Baryon", 23122);
1281  new TAttParticle("Lambda(1600)0 bar","B101",
1282  1.6,kFALSE, .15,
1283  0.0,"Baryon", -23122);
1284  new TAttParticle("Lambda(1890)0","B060",
1285  1.89,kFALSE, .1,
1286  0.0,"Baryon", 23124);
1287  new TAttParticle("Lambda(1890)0 bar","B060",
1288  1.89,kFALSE, .1,
1289  0.0,"Baryon", -23124);
1290  new TAttParticle("Lambda(2110)0","B035",
1291  2.11,kFALSE, .2,
1292  0.0,"Baryon", 23126);
1293  new TAttParticle("Lambda(2110)0 bar","B035",
1294  2.11,kFALSE, .2,
1295  0.0,"Baryon", -23126);
1296  new TAttParticle("Sigma(1750)0","B057",
1297  1.75,kFALSE, 9.00000E-02,
1298  0.0,"Baryon", 23212);
1299  new TAttParticle("Sigma(1750)0 bar","B057",
1300  1.75,kFALSE, 9.00000E-02,
1301  0.0,"Baryon", -23212);
1302  new TAttParticle("Sigma(1940)0","B098",
1303  1.94,kFALSE, .2199999,
1304  0.0,"Baryon", 23214);
1305  new TAttParticle("Sigma(1940)0 bar","B098",
1306  1.94,kFALSE, .2199999,
1307  0.0,"Baryon", -23214);
1308  new TAttParticle("Sigma(1750)+","B057",
1309  1.75,kFALSE, 9.00000E-02,
1310  +1.,"Baryon", 23222);
1311  new TAttParticle("Sigma(1750)- bar","B057",
1312  1.75,kFALSE, 9.00000E-02,
1313  -1.,"Baryon", -23222);
1314  new TAttParticle("Sigma(1940)+","B098",
1315  1.94,kFALSE, .2199999,
1316  +1.,"Baryon", 23224);
1317  new TAttParticle("Sigma(1940)- bar","B098",
1318  1.94,kFALSE, .2199999,
1319  -1.,"Baryon", -23224);
1320  new TAttParticle("rho(1700)0","M065",
1321  1.7,kFALSE, .24,
1322  0.0,"Meson", 30113);
1323  new TAttParticle("rho(1700)+","M065",
1324  1.7,kFALSE, .24,
1325  1.,"Meson", 30213);
1326  new TAttParticle("rho(1700)-","M065",
1327  1.7,kFALSE, .24,
1328  -1.,"Meson", -30213);
1329  new TAttParticle("f(1)(1420)0","M006",
1330  1.4268,kFALSE, 5.20000E-02,
1331  0.0,"Meson", 30223);
1332  new TAttParticle("f(2)(2300)0","M107",
1333  2.297,kFALSE, .15,
1334  0.0,"Meson", 30225);
1335  new TAttParticle("K*(1410)0","M094",
1336  1.412,kFALSE, .2269999,
1337  0.0,"Meson", 30313);
1338  new TAttParticle("K*(1410)0 bar","M094",
1339  1.412,kFALSE, .2269999,
1340  0.0,"Meson", -30313);
1341  new TAttParticle("K*(1410)+","M094",
1342  1.412,kFALSE, .2269999,
1343  1.,"Meson", 30323);
1344  new TAttParticle("K*(1410)-","M094",
1345  1.412,kFALSE, .2269999,
1346  -1.,"Meson", -30323);
1347  new TAttParticle("psi(3770)0","M053",
1348  3.7699,kFALSE, 2.36000E-02,
1349  0.0,"Meson", 30443);
1350  new TAttParticle("Upsilon(3S)0","M048",
1351  10.35529,kFALSE, 2.63000E-05,
1352  0.0,"Meson", 30553);
1353  new TAttParticle("Delta(1600)-","B019",
1354  1.6,kFALSE, .3499999,
1355  -1.,"Baryon", 31114);
1356  new TAttParticle("Delta(1600)+ bar","B019",
1357  1.6,kFALSE, .3499999,
1358  +1.,"Baryon", -31114);
1359  new TAttParticle("N(1720)0","B015",
1360  1.72,kFALSE, .15,
1361  0.0,"Baryon", 31214);
1362  new TAttParticle("N(1720)0 bar","B015",
1363  1.72,kFALSE, .15,
1364  0.0,"Baryon", -31214);
1365  new TAttParticle("N(1650)0","B066",
1366  1.65,kFALSE, .15,
1367  0.0,"Baryon", 32112);
1368  new TAttParticle("N(1650)0 bar","B066",
1369  1.65,kFALSE, .15,
1370  0.0,"Baryon", -32112);
1371  new TAttParticle("Delta(1600)0","B019",
1372  1.6,kFALSE, .3499999,
1373  0.0,"Baryon", 32114);
1374  new TAttParticle("Delta(1600)0 bar","B019",
1375  1.6,kFALSE, .3499999,
1376  0.0,"Baryon", -32114);
1377  new TAttParticle("N(1720)+","B015",
1378  1.72,kFALSE, .15,
1379  +1.,"Baryon", 32124);
1380  new TAttParticle("N(1720)- bar","B015",
1381  1.72,kFALSE, .15,
1382  -1.,"Baryon", -32124);
1383  new TAttParticle("N(1650)+","B066",
1384  1.65,kFALSE, .15,
1385  +1.,"Baryon", 32212);
1386  new TAttParticle("N(1650)- bar","B066",
1387  1.65,kFALSE, .15,
1388  -1.,"Baryon", -32212);
1389  new TAttParticle("Delta(1600)+","B019",
1390  1.6,kFALSE, .3499999,
1391  +1.,"Baryon", 32214);
1392  new TAttParticle("Delta(1600)- bar","B019",
1393  1.6,kFALSE, .3499999,
1394  -1.,"Baryon", -32214);
1395  new TAttParticle("Delta(1600)++","B019",
1396  1.6,kFALSE, .3499999,
1397  +2.,"Baryon", 32224);
1398  new TAttParticle("Delta(1600)-- bar","B019",
1399  1.6,kFALSE, .3499999,
1400  -2.,"Baryon", -32224);
1401  new TAttParticle("Lambda(1670)0","B040",
1402  1.67,kFALSE, 3.50000E-02,
1403  0.0,"Baryon", 33122);
1404  new TAttParticle("Lambda(1670)0 bar","B040",
1405  1.67,kFALSE, 3.50000E-02,
1406  0.0,"Baryon", -33122);
1407  new TAttParticle("rho(1450)0","M105",
1408  1.465,kFALSE, .31,
1409  0.0,"Meson", 40113);
1410  new TAttParticle("rho(1450)+","M105",
1411  1.465,kFALSE, .31,
1412  1.,"Meson", 40213);
1413  new TAttParticle("rho(1450)-","M105",
1414  1.465,kFALSE, .31,
1415  -1.,"Meson", -40213);
1416  new TAttParticle("eta(1440)0","M027",
1417  1.42,kFALSE, 6.00000E-02,
1418  0.0,"Meson", 40221);
1419  new TAttParticle("f(1)(1510)0","M084",
1420  1.512,kFALSE, 3.50000E-02,
1421  0.0,"Meson", 40223);
1422  new TAttParticle("f(2)(2340)0","M108",
1423  2.34,kFALSE, .3199999,
1424  0.0,"Meson", 40225);
1425  new TAttParticle("K*(1680)0","M095",
1426  1.714,kFALSE, .3199999,
1427  0.0,"Meson", 40313);
1428  new TAttParticle("K*(1680)0 bar","M095",
1429  1.714,kFALSE, .3199999,
1430  0.0,"Meson", -40313);
1431  new TAttParticle("K*(1680)+","M095",
1432  1.714,kFALSE, .3199999,
1433  1.,"Meson", 40323);
1434  new TAttParticle("K*(1680)-","M095",
1435  1.714,kFALSE, .3199999,
1436  -1.,"Meson", -40323);
1437  new TAttParticle("psi(4040)0","M072",
1438  4.04,kFALSE, 5.20000E-02,
1439  0.0,"Meson", 40443);
1440  new TAttParticle("Upsilon(4S)0","M047",
1441  10.57999,kFALSE, 2.38000E-02,
1442  0.0,"Meson", 40553);
1443  new TAttParticle("N(1710)0","B014",
1444  1.71,kFALSE, .1,
1445  0.0,"Baryon", 42112);
1446  new TAttParticle("N(1710)0 bar","B014",
1447  1.71,kFALSE, .1,
1448  0.0,"Baryon", -42112);
1449  new TAttParticle("N(1710)+","B014",
1450  1.71,kFALSE, .1,
1451  +1.,"Baryon", 42212);
1452  new TAttParticle("N(1710)- bar","B014",
1453  1.71,kFALSE, .1,
1454  -1.,"Baryon", -42212);
1455  new TAttParticle("Lambda(1800)0","B036",
1456  1.8,kFALSE, .3,
1457  0.0,"Baryon", 43122);
1458  new TAttParticle("Lambda(1800)0 bar","B036",
1459  1.8,kFALSE, .3,
1460  0.0,"Baryon", -43122);
1461  new TAttParticle("f(0)(1590)0","M096",
1462  1.581,kFALSE, .18,
1463  0.0,"Meson", 50221);
1464  new TAttParticle("omega(1420)0","M125",
1465  1.419,kFALSE, .17,
1466  0.0,"Meson", 50223);
1467  new TAttParticle("psi(4160)0","M025",
1468  4.159,kFALSE, 7.80000E-02,
1469  0.0,"Meson", 50443);
1470  new TAttParticle("Upsilon(10860)0","M092",
1471  10.86499,kFALSE, .1099999,
1472  0.0,"Meson", 50553);
1473  new TAttParticle("Lambda(1810)0","B077",
1474  1.81,kFALSE, .15,
1475  0.0,"Baryon", 53122);
1476  new TAttParticle("Lambda(1810)0 bar","B077",
1477  1.81,kFALSE, .15,
1478  0.0,"Baryon", -53122);
1479  new TAttParticle("f(J)(1710)0","M068",
1480  1.709,kFALSE, .14,
1481  0.0,"Meson", 60221);
1482  new TAttParticle("omega(1600)0","M126",
1483  1.662,kFALSE, .28,
1484  0.0,"Meson", 60223);
1485  new TAttParticle("psi(4415)0","M073",
1486  4.415,kFALSE, 4.30000E-02,
1487  0.0,"Meson", 60443);
1488  new TAttParticle("Upsilon(11020)0","M093",
1489  11.019,kFALSE, 7.90000E-02,
1490  0.0,"Meson", 60553);
1491  new TAttParticle("chi(b1)(2P)0","M080",
1492  10.2552,kTRUE, .0,
1493  0.0,"Meson", 70553);
1494 // End of the entry point of the pdg table conversion
1495  new TAttParticle("Rootino","",
1496  0.0,kTRUE,
1497  1.e38,0.0,"Artificial",0);
1498 }
1499 
1500 ////////////////////////////////////////////////////////////////////////////////
1501 ///
1502 /// Get a pointer to the particle object according to the name given
1503 ///
1504 
1505 TAttParticle* TAttParticle::GetParticle(const char *name)
1506 {
1507  TAttParticle *def = (TAttParticle *)fgList->FindObject(name);
1508  if (!def) {
1509  fgList->Error("GetParticle","No match for %s exists !",name);
1510  }
1511  return def;
1512 }
1513 
1514 ////////////////////////////////////////////////////////////////////////////////
1515 ///
1516 /// Get a pointer to the particle object according to the MC code number
1517 ///
1518 
1519 TAttParticle* TAttParticle::GetParticle(Int_t mcnumber)
1520 {
1521  TIter next(fgList);
1522  TAttParticle *par;
1523  while ((par = (TAttParticle *)next())) {
1524  if (par->GetMCNumber() == mcnumber) return par;
1525  }
1526  fgList->Error("GetParticle","No match for %d exists !",mcnumber);
1527  return 0;
1528 }
1529 
1530 ////////////////////////////////////////////////////////////////////////////////
1531 ///
1532 /// Print the entire information of this kind of particle
1533 ///
1534 
1535 void TAttParticle::Print(Option_t *) const
1536 {
1537  Printf("\nParticle: %-15s ",
1538  this->GetName());
1539  if (!fPDGStable) {
1540  Printf("Mass: %8f DecayWidth: %8f Charge : %8f",
1541  fPDGMass, fPDGDecayWidth, fPDGCharge);
1542  }
1543  else {
1544  Printf("Mass: %8f DecayWidth: Stable Charge : %8f",
1545  fPDGMass, fPDGCharge);
1546  }
1547  Printf(" ");
1548 }
1549 
1550 ////////////////////////////////////////////////////////////////////////////////
1551 ///
1552 /// Samples a mass according to the Breit-Wigner resonance distribution
1553 ///
1554 
1555 Double_t TAttParticle::SampleMass() const
1556 {
1557  if ( fPDGStable || fPDGDecayWidth == 0.0 )
1558  return fPDGMass;
1559  else {
1560  return (fPDGMass+
1561  0.5*fPDGDecayWidth*
1562  TMath::Tan((2.0*gRandom->Rndm()-1.0)*TMath::Pi()*0.5));
1563  }
1564 }
1565 
1566 ////////////////////////////////////////////////////////////////////////////////
1567 ///
1568 /// Samples a mass in the interval:
1569 ///
1570 /// fPDGMass-widthcut*fPDGDecayWidtht - fPDGMass+widthcut*fPDGDecayWidth
1571 ///
1572 /// according to the Breit-Wigner resonance distribution
1573 ///
1574 
1575 Double_t TAttParticle::SampleMass(Double_t widthcut) const
1576 {
1577  if ( fPDGStable || fPDGDecayWidth == 0.0 )
1578  return fPDGMass;
1579  else {
1580  return (fPDGMass+
1581  0.5*fPDGDecayWidth*
1582  TMath::Tan((2.0*gRandom->Rndm()-1.0)*TMath::ATan(2.0*widthcut)));
1583  }
1584 }