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