Logo ROOT  

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/** \class TAttParticle
13 \ingroup eg
14
15Particle 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
23ClassImp(TAttParticle);
24
25THashList *TAttParticle::fgList = new THashList;
26
27////////////////////////////////////////////////////////////////////////////////
28///
29/// Particle definition default constructor
30///
31
32TAttParticle::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
56TAttParticle::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
80TAttParticle::~TAttParticle()
81{
82}
83
84////////////////////////////////////////////////////////////////////////////////
85///
86/// Converts the ISAJET Particle number into the PDG MC number
87///
88
89Int_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
249void 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
1505TAttParticle* 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
1519TAttParticle* 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
1535void 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
1555Double_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
1575Double_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}
Int_t
int Int_t
Definition: RtypesCore.h:41
kFALSE
const Bool_t kFALSE
Definition: RtypesCore.h:88
Bool_t
bool Bool_t
Definition: RtypesCore.h:59
Double_t
double Double_t
Definition: RtypesCore.h:55
kTRUE
const Bool_t kTRUE
Definition: RtypesCore.h:87
Option_t
const char Option_t
Definition: RtypesCore.h:62
ClassImp
#define ClassImp(name)
Definition: Rtypes.h:365
name
char name[80]
Definition: TGX11.cxx:109
gRandom
R__EXTERN TRandom * gRandom
Definition: TRandom.h:62
Printf
void Printf(const char *fmt,...)
TAttParticle
Particle definition, partly based on GEANT3 particle definition.
Definition: TAttParticle.h:30
TAttParticle::GetParticle
static TAttParticle * GetParticle(const char *name)
Get a pointer to the particle object according to the name given.
Definition: TAttParticle.cxx:1505
TAttParticle::GetMCNumber
virtual Int_t GetMCNumber() const
Definition: TAttParticle.h:60
TAttParticle::TAttParticle
TAttParticle()
Particle definition default constructor.
Definition: TAttParticle.cxx:32
TAttParticle::fEnergyCut
Double_t fEnergyCut
Definition: TAttParticle.h:39
TAttParticle::ConvertISAtoPDG
static Int_t ConvertISAtoPDG(Int_t isaNumber)
Converts the ISAJET Particle number into the PDG MC number.
Definition: TAttParticle.cxx:89
TAttParticle::fPDGMass
Double_t fPDGMass
Definition: TAttParticle.h:32
TAttParticle::SampleMass
virtual Double_t SampleMass() const
Samples a mass according to the Breit-Wigner resonance distribution.
Definition: TAttParticle.cxx:1555
TAttParticle::DefinePDG
static void DefinePDG()
Defines particles according to the Particle Data Group.
Definition: TAttParticle.cxx:249
TAttParticle::fPDGDecayWidth
Double_t fPDGDecayWidth
Definition: TAttParticle.h:34
TAttParticle::fgList
static THashList * fgList
Definition: TAttParticle.h:51
TAttParticle::~TAttParticle
virtual ~TAttParticle()
Particle destructor.
Definition: TAttParticle.cxx:80
TAttParticle::Print
virtual void Print(Option_t *option="") const
Print the entire information of this kind of particle.
Definition: TAttParticle.cxx:1535
TAttParticle::fPDGCharge
Double_t fPDGCharge
Definition: TAttParticle.h:35
TAttParticle::fPDGStable
Bool_t fPDGStable
Definition: TAttParticle.h:33
TAttParticle::fMCnumberOfPDG
Int_t fMCnumberOfPDG
Definition: TAttParticle.h:37
TAttParticle::fParticleType
TString fParticleType
Definition: TAttParticle.h:36
TAttParticle::fEnergyLimit
Double_t fEnergyLimit
Definition: TAttParticle.h:40
TAttParticle::fGranularity
Double_t fGranularity
Definition: TAttParticle.h:41
TCollection::IsEmpty
virtual Bool_t IsEmpty() const
Definition: TCollection.h:186
THashList
THashList implements a hybrid collection class consisting of a hash table and a list to store TObject...
Definition: THashList.h:34
THashList::FindObject
TObject * FindObject(const char *name) const
Find object using its name.
Definition: THashList.cxx:262
TList::Add
virtual void Add(TObject *obj)
Definition: TList.h:87
TNamed
The TNamed class is the base class for all named ROOT classes.
Definition: TNamed.h:29
TNamed::GetName
virtual const char * GetName() const
Returns name of object.
Definition: TNamed.h:47
TObject::Error
virtual void Error(const char *method, const char *msgfmt,...) const
Issue error message.
Definition: TObject.cxx:880
TRandom::Rndm
virtual Double_t Rndm()
Machine independent random number generator.
Definition: TRandom.cxx:541
ROOT::Math::IntegrationOneDim::Type
Type
enumeration specifying the integration types.
Definition: AllIntegrationTypes.h:45
TMath::ATan
Double_t ATan(Double_t)
Definition: TMath.h:665
TMath::Pi
constexpr Double_t Pi()
Definition: TMath.h:38
TMath::Tan
Double_t Tan(Double_t)
Definition: TMath.h:635