Logo ROOT  

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