CbmRoot
Loading...
Searching...
No Matches
PStdData.cxx
Go to the documentation of this file.
1/* Copyright (C) 1999-2019 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
2 SPDX-License-Identifier: GPL-3.0-only
3 Authors: M.A. Kagarlis, Florian Uhlig [committer] */
4
6// Pluto Particle Data (The standard version)
7//
8// This class contains a permanent data base of particles
9// and their properties, functions to load additional
10// temporary particles and decay modes, as well as functions
11// to build decay widths and branching ratios and sample masses
12// of arbitrary hadronic resonances.
13// The PID convention is consistent with GEANT, with up to
14// 999 particles supported.
15//
16// Author: M.A. Kagarlis
17// Written: 31/01/99
18// Revised: 24/05/2004 R.H.
19// Revised particle tables: 09/06/2004 R.H.
20// Bug fixes: 13/10/2004 R.H.
21// mu+mu-: 02/02/2005 R.H.
22// N1535: 25/01/2007 R.H.
23//
24// The static arrays habe been copied from the original PData class
25// to a "filler" to have more flexibility (IF 28.4.2007)
26//
27// Ref 1: Particle Data Group, Review of Particle Properties,
28// Phys. Rev. D54 (1996) 1 (and earlier editions)
30
31// local arrays
32#include "PStdData.h"
33
34#include "PDataBase.h"
35#include "TArrayD.h"
36#include "TArrayI.h"
37
38#include <cmath>
39
40PStdData& fStdData()
41{
42 static PStdData* ans = new PStdData();
43 return *ans;
44}
45
46PStdData* makeStdData() { return &fStdData(); }
47
48PStdData::PStdData()
49{
50 // Info("PStdData()","(%s)", PRINT_CTOR);
51
52 disable = 0;
53}
54
55PStdData::~PStdData() { cout << "Killing standard particle data table" << endl; }
56
57
58Bool_t PStdData::fillDataBase(void)
59{
60
61 if (disable) return kTRUE;
62
63 PPosition = new int[maxnumpar];
64 resetPosition();
65
66 PDataBase* base = makeDataBase();
67 //make particle header
68 base->MakeParamInt("setx", "Particle set index");
69 base->MakeParamInt("snpart", "Number of particles");
70 base->MakeParamInt("slink", "Particle link");
71
72 Int_t pidkey = base->MakeParamInt("pid", "Pluto particle ID");
73
74 base->SetFastKey(pidkey, 100000);
75 base->MakeParamDouble("width", "Particle static width [GeV]");
76 base->MakeParamDouble("mass", "Particle pole mass [GeV]");
77 base->MakeParamInt("meson", "Meson number");
78 base->MakeParamInt("baryon", "Baryon number");
79 base->MakeParamInt("lepton", "Lepton number");
80 base->MakeParamInt("charge", "Charge");
81 base->MakeParamInt("spin", "Spin");
82 base->MakeParamInt("parity", "Parity");
83 base->MakeParamInt("ispin", "Isospin");
84 base->MakeParamInt("pythiakf", "Pythia KF code");
85
86 base->MakeParamInt("pnmodes", "Number of decay modes");
87 base->MakeParamInt("link", "Decay link");
88 base->MakeParamInt("d1", "Decay product 1");
89 base->MakeParamInt("d2", "Decay product 2");
90 base->MakeParamInt("d3", "Decay product 3");
91 base->MakeParamInt("d4", "Decay product 4");
92 base->MakeParamInt("d5", "Decay product 5");
93 base->MakeParamInt("d6", "Decay product 6");
94 base->MakeParamInt("d7", "Decay product 7");
95 base->MakeParamInt("ppid", "Parent Pluto ID"); //parent pid
96 Int_t idxkey = base->MakeParamInt("didx", "Decay index");
97 base->MakeParamDouble("br", "Branching ratio");
98 base->MakeParamDouble("brorig", "Original branching ratio (for norm.)");
99 base->SetFastKey(idxkey, 1000 * 32);
100
101 base->MakeParamInt("widx", "Width flag"); //partial/total width switched on by default
102 base->MakeParamInt("tdepth", "total depth");
103 base->MakeParamInt("hdepth", "hadronic depth");
104 base->MakeParamDouble("ethreshold", "Energy threshold");
105 base->MakeParamDouble("scfactor", "Self consistency factor");
106 base->MakeParamInt("sccount", "Self consistency count (max tries)");
107
108
109 base->MakeParamTObj("mesh", "Mesh object");
110 base->MakeParamTObj("tf1", "TF1 object");
111 base->MakeParamTObj("model", "Basic PModel");
112
113 base->MakeParamInt("maxmesh", "Number of Mesh points");
114 base->MakeParamDouble("lmass", "Lower mass");
115 base->MakeParamDouble("umass", "Upper mass");
116
117 base->MakeParamInt("brflag", "=1: use integral normalization");
118
119 //alias
120 base->MakeParamInt("nalias", "Number of aliases");
121 base->MakeParamInt("lalias", "Alias link");
122
123 //this is for the decay manager
124 base->MakeParamTObj("decaychannel", "Decay channel for decay manager");
125
126 //...for the decay_all option
127 base->MakeParamTObj("stackchannel", "Stack of channels for the decayall option in PReaction");
128
129 Int_t skey = -1;
130 if ((skey = base->AddEntry("std_set")) < 0) return kFALSE;
131
132
133 int pkey = -1, *ii;
134 Double_t* dd;
135 for (int i = 0; i < maxnumpar; i++) {
136 // if ((pkey=base->AddEntry(PStdData::PName[i]))<0) return kFALSE;
137 // Int_t dkey = base->AddListEntry(PStdData::PName[i],"pnmodes", "link",PMDescription [j]);
138 if ((pkey = base->AddListEntry("std_set", "snpart", "slink", PStdData::PName[i])) < 0) return kFALSE;
139 ii = new int(i); //never destructed, but called only once!
140 if (!base->SetParamInt(pkey, "pid", ii)) return kFALSE;
141 if (!base->SetParamDouble(pkey, "width", &(PStdData::PWidth[i]))) return kFALSE;
142 if (!base->SetParamDouble(pkey, "mass", &(PStdData::PMass[i]))) return kFALSE;
143 if (PMeson[i])
144 if (!base->SetParamInt(pkey, "meson", &(PStdData::PMeson[i]))) return kFALSE;
145 if (PBaryon[i])
146 if (!base->SetParamInt(pkey, "baryon", &(PStdData::PBaryon[i]))) return kFALSE;
147 if (PLepton[i])
148 if (!base->SetParamInt(pkey, "lepton", &(PStdData::PLepton[i]))) return kFALSE;
149 if (!base->SetParamInt(pkey, "charge", &(PStdData::PCharge[i]))) return kFALSE;
150 if (!base->SetParamInt(pkey, "spin", &(PStdData::PJ[i]))) return kFALSE;
151 if (!base->SetParamInt(pkey, "parity", &(PStdData::PParity[i]))) return kFALSE;
152 if (!base->SetParamInt(pkey, "ispin", &(PStdData::PI[i]))) return kFALSE;
153 if (!base->SetParamInt(pkey, "pythiakf", &(PStdData::Pkf[i]))) return kFALSE;
154 ii = new int(0); //never destructed, but called only once!
155 //0 means on, -1 means off
156 if (!base->SetParamInt(pkey, "widx", ii)) return kFALSE;
157 ii = new int(0);
158 if (!base->SetParamInt(pkey, "tdepth", ii)) return kFALSE;
159 ii = new int(0);
160 if (!base->SetParamInt(pkey, "hdepth", ii)) return kFALSE;
161 dd = new Double_t(PStdData::PMass[i]
162 - 2 * PStdData::PWidth[i]); //BUGBUG look to ethreshold later -> only 1st guess!!!
163 if (!base->SetParamDouble(pkey, "ethreshold", dd)) return kFALSE;
164
165 //Adding Fireballs!
166 size_t buf_size = 100;
167 char* name = new char[buf_size];
168 snprintf(name, buf_size - 1, "Fireball: %s", PStdData::PName[i]);
169 if ((pkey = base->AddListEntry("std_set", "snpart", "slink", name)) < 0) return kFALSE;
170 ii = new int(i + 500); //never destructed, but called only once!
171 if (!base->SetParamInt(pkey, "pid", ii)) return kFALSE;
172 if (!base->SetParamDouble(pkey, "mass", &(PStdData::PMass[i]))) return kFALSE;
173 if (!base->SetParamDouble(pkey, "width", &(PStdData::PWidth[i]))) return kFALSE;
174 }
175
176 for (int i = 0; i < maxnumpar; i++) {
177 // Int_t pkey=base->getEntry(PStdData::PName[i]);
178 //now the decays
179 int jmin = 0, jmax = PStdData::maxnummodes;
180 jmin = PPosition[i];
181 jmax = jmin + PStdData::PNModes[i];
182 //debug info
183 // cout << "decay of " << PStdData::PName[i] << endl;
184
185 for (int j = jmin; j < jmax; ++j) {
186
187 TString s = PMode[j]; // retrieve decay mode string
188
189 Int_t dkey = base->AddListEntry(PStdData::PName[i], "pnmodes", "link", PMDescription[j]);
190
191 ii = new int(i);
192 base->SetParamInt(dkey, "ppid", ii); //set parent id
193
194 ii = new int(j + 1); //Shifted by 1 because 0 is not good
195 base->SetParamInt(dkey, "didx", ii); //decay mode index
196
197 dd = new Double_t(0.); //BUGBUG look to ethreshold
198 if (!base->SetParamDouble(dkey, "ethreshold", dd)) return kFALSE;
199
200 ii = new int(0); //never destructed, but called only once!
201 if (!base->SetParamInt(dkey, "widx", ii)) return kFALSE;
202
203 if (!base->SetParamDouble(dkey, "br", &(PStdData::PBR[j]))) return kFALSE;
204
205 if (!base->SetParamDouble(dkey, "brorig", new Double_t(PStdData::PBR[j]))) return kFALSE;
206
207 dd = new Double_t(0.);
208 if (!base->SetParamDouble(dkey, "ethreshold", dd)) return kFALSE;
209
210 dd = new Double_t(0.); //partial width resetted later
211 if (!base->SetParamDouble(dkey, "width", dd)) return kFALSE;
212
213 dd = new Double_t(1.); //sc factor=1.
214 if (!base->SetParamDouble(dkey, "scfactor", dd)) return kFALSE;
215
216 Int_t len = s.Length();
217 Int_t res = len % 3;
218 Int_t np = len / 3 + (res > 0);
219 Int_t pid, k = 0;
220 for (int iPart = 0; iPart < np; ++iPart) { // loop over product particles
221 Int_t m = (!iPart && res) ? res : 3; // number of digits in current pid
222 pid = 0; // reset id
223 for (int jj = 1; jj <= m; ++jj)
224 pid += (*(s(k + jj - 1, 1).Data()) - 48) * int(pow(10., m - jj));
225 //cout << s << endl;
226 //cout << pid << endl;
227 Int_t* pkey1 = new int(base->GetEntryInt("pid", pid));
228 if (*pkey1 < 0) { cout << "Error: processing decay: do not find pid " << pid << endl; }
229 if (iPart == 0) base->SetParamInt(dkey, "d1", pkey1);
230 if (iPart == 1) base->SetParamInt(dkey, "d2", pkey1);
231 if (iPart == 2) base->SetParamInt(dkey, "d3", pkey1);
232 if (iPart == 3) base->SetParamInt(dkey, "d4", pkey1);
233 if (iPart == 4) base->SetParamInt(dkey, "d5", pkey1);
234 if (iPart == 5) base->SetParamInt(dkey, "d6", pkey1);
235 if (iPart == 6) base->SetParamInt(dkey, "d7", pkey1);
236 if (iPart > 6)
237 cout << "PStdData::fillDataBase: More then 7 decay products not "
238 "supported:"
239 << PMDescription[j] << endl;
240 k += m;
241
242 } //end decay products
243
244 } //end decay mode loop
245 }
246
247
248 disable = 1;
249 return kTRUE;
250}
251
252
253//Please do not touch the tables below
254
255#define nmax 999 // maximum number of supported particles
256#define mnpar 70 // number of particles stored permanently
257#define mnmodes 220 // number of decay modes stored permanently
258
259
260// particle naming convention for use with the PParticle constructor
261const char* PStdData::NAME[mnpar] = {
262 "dummy", "g", "e+", "e-", "nu", "mu+", "mu-", "pi0", "pi+",
263 "pi-", "K0L", "K+", "K-", "n", "p", "anti_p", "K0S", "eta",
264 "Lambda", "Sigma+", "Sigma0", "Sigma-", "Xi0", "Xi-", "Omega", "anti_n", "anti_Lambda",
265 "anti_Sigma-", "anti_Sigma0", "anti_Sigma+", "anti_Xi0", "anti_Xi+", "anti_Omega+", "File", "D0", "D++",
266 "D+", "D-", "NP11+", "ND13+", "NS11+", "rho0", "rho+", "rho-", "BOZO",
267 "d", "t", "alpha", "BOZO2", "He3", //BOZO->BOZO2
268 "dimuon", "dilepton", "w", "eta'", "sigma", "phi", "DP330", "DP33++", "DP33+",
269 "DP33-", "DS310", "DS31++", "DS31+", "DS31-", "NP110", "ND130", "NS110", "J/Psi",
270 "Psi'", "pn"};
271
272// particle masses (GeV/c**2)
273//const double MASS[mnpar]={
274double PStdData::MASS[mnpar] = {
275 /* 0: dummy */ 0.0,
276 /* 1: Photon */ 0.0,
277 /* 2: Positron */ 0.51099906e-3,
278 /* 3: Electron */ 0.51099906e-3,
279 /* 4: Neutrino */ 0.0,
280 /* 5: mu+ */ 0.105658389,
281 /* 6: mu- */ 0.105658389,
282 /* 7: pi0 */ 0.1349764,
283 /* 8: pi+ */ 0.13956995,
284 /* 9: pi- */ 0.13956995,
285 /*10: K0 long */ 0.497672,
286 /*11: K+ */ 0.493677,
287 /*12: K- */ 0.493677,
288 /*13: Neutron */ 0.93956563,
289 /*14: Proton */ 0.93827231,
290 /*15: Antiproton*/ 0.93827231,
291 /*16: K0 short */ 0.497672,
292 /*17: Eta */ 0.54745,
293 /*18: Lambda */ 1.115684,
294 /*19: Sigma+ */ 1.18937,
295 /*20: Sigma0 */ 1.19255,
296 /*21: Sigma- */ 1.197436,
297 /*22: Xi0 */ 1.3149,
298 /*23: Xi- */ 1.32132,
299 /*24: Omega */ 1.67245,
300 /*25: Antineutrn*/ 0.93956563,
301 /*26: Antilambda*/ 1.115684,
302 /*27: Antisigma-*/ 1.18937, //1.197436,
303 /*28: Antisigma0*/ 1.19255,
304 /*29: Antisigma+*/ 1.197436, //1.18937,
305 /*30: Antixi0 */ 1.3149,
306 /*31: Antixi+ */ 1.32132,
307 /*32: Antiomega+*/ 1.67245,
308 /*33: File */ 0.0,
309 /*34: Delta0 */ 1.232,
310 /*35: Delta++ */ 1.232,
311 /*36: Delta+ */ 1.232,
312 /*37: Delta- */ 1.232,
313 /*38: NP11+ */ 1.44,
314 /*39: ND13+ */ 1.520,
315 /*40: NS11+ */ 1.535,
316 /*41: rho0 */ 0.7699,
317 /*42: rho+ */ 0.7699,
318 /*43: rho- */ 0.7699,
319 /*44: NULL */ 0.0,
320 /*45: Deuteron */ 1.875613,
321 /*46: Tritium */ 2.80925,
322 /*47: Alpha */ 3.727417,
323 /*48: NULL */ 0.0,
324 /*49: He3 */ 2.80923,
325 /*50: dimuon */ 0.21131678,
326 /*51: dilepton */ 0.001022,
327 /*52: omega */ 0.78194,
328 /*53: eta' */ 0.9577,
329 /*54: sigma */ 0.6,
330 /*55: phi */ 1.019413,
331 /*56: Delta0*P33*/ 1.6,
332 /*57: Delta++ *P33*/ 1.6,
333 /*58: Delta+*P33*/ 1.6,
334 /*59: Delta- *P33 */ 1.6,
335 /*60: Delta0*S31*/ 1.62,
336 /*61: Delta++ *S31*/ 1.62,
337 /*62: Delta+*S31*/ 1.62,
338 /*63: Delta- *S31 */ 1.62,
339 /*64: NP110 */ 1.44,
340 /*65: ND130 */ 1.520,
341 /*66: NS110 */ 1.535,
342 /*67: J/Psi */ 3.09688,
343 /*68: Psi' */ 3.68596,
344 /*69: pn */ 2.65
345
346};
347
348// particle widths (GeV/c**2)
349//const double WIDTH[mnpar]={
350double PStdData::WIDTH[mnpar] = {
351 /* 0: dummy */ 0.0,
352 /* 1: Photon */ 0.0,
353 /* 2: Positron */ 0.0,
354 /* 3: Electron */ 0.0,
355 /* 4: Neutrino */ 0.0,
356 /* 5: mu+ */ hbar / 2.19703e-6,
357 /* 6: mu- */ PStdData::hbar / 2.19703e-6,
358 /* 7: pi0 */ PStdData::hbar / 8.4e-17,
359 /* 8: pi+ */ PStdData::hbar / 2.6033e-8,
360 /* 9: pi- */ PStdData::hbar / 2.6033e-8,
361 /*10: K0 long */ PStdData::hbar / 5.17e-8,
362 /*11: K+ */ PStdData::hbar / 1.2386e-8,
363 /*12: K- */ PStdData::hbar / 1.23861e-8,
364 /*13: Neutron */ PStdData::hbar / 887.,
365 /*14: Proton */ 0.0,
366 /*15: Antiproton*/ 0.0,
367 /*16: K0 short */ PStdData::hbar / 8.927e-11,
368 /*17: eta */ 1.29e-6,
369 /*18: Lambda */ PStdData::hbar / 2.632e-10,
370 /*19: Sigma+ */ PStdData::hbar / 7.99e-11,
371 /*20: Sigma0 */ PStdData::hbar / 7.4e-20,
372 /*21: Sigma- */ PStdData::hbar / 1.479e-10,
373 /*22: Xi0 */ PStdData::hbar / 2.9e-10,
374 /*23: Xi- */ PStdData::hbar / 1.639e-10,
375 /*24: Omega- */ PStdData::hbar / 8.22e-11,
376 /*25: Antineutrn*/ PStdData::hbar / 887.,
377 /*26: Antilambda*/ PStdData::hbar / 2.632e-10,
378 /*27: Antisigma-*/ PStdData::hbar / 7.99e-11,
379 /*28: Antisigma0*/ PStdData::hbar / 7.4e-20,
380 /*29: Antisigma+*/ PStdData::hbar / 1.479e-10,
381 /*30: Antixi0 */ PStdData::hbar / 2.9e-10,
382 /*31: Antixi+ */ PStdData::hbar / 1.639e-10,
383 /*32: Antiomega+*/ PStdData::hbar / 8.22e-11,
384 /*33: File */ 0.0,
385 /*34: Delta0 */ 0.12,
386 /*35: Delta++ */ 0.12,
387 /*36: Delta+ */ 0.12,
388 /*37: Delta- */ 0.12,
389 /*38: NP11+ */ 0.35,
390 /*39: ND13+ */ 0.12,
391 /*40: NS11+ */ 0.15,
392 /*41: rho0 */ 0.1507,
393 /*42: rho+ */ 0.1507,
394 /*43: rho- */ 0.1507,
395 /*44: NULL */ 0.0,
396 /*45: Deuteron */ 0.0,
397 /*46: Tritium */ 0.0,
398 /*47: Alpha */ 0.0,
399 /*48: NULL */ 0.0,
400 /*49: He3 */ 0.0,
401 /*50: dimuon */ 0.0,
402 /*51: dilepton */ 0.0,
403 /*52: omega */ 0.00843,
404 /*53: eta' */ 0.000201,
405 /*54: sigma */ 0.5,
406 /*55: phi */ 0.00443,
407 /*56: Delta0*P33*/ 0.35,
408 /*57: Delta++ *P33*/ 0.35,
409 /*58: Delta+*P33*/ 0.35,
410 /*59: Delta- *P33 */ 0.35,
411 /*60: Delta0*S31*/ 0.15,
412 /*61: Delta++ *S31*/ 0.15,
413 /*62: Delta+*S31*/ 0.15,
414 /*63: Delta- *S31 */ 0.15,
415 /*64: NP110 */ 0.35,
416 /*65: ND130 */ 0.12,
417 /*66: NS110 */ 0.15,
418 /*67: J/Psi */ 0.000087,
419 /*68: Psi' */ 0.000277,
420 /*69: pn */ 0.5};
421
422// Pythia6 KF code
423const int PStdData::PYTHIAKF[mnpar] = {
424 /* 0: dummy */ 0, /* 1: Photon */ 22,
425 /* 2: Positron */ -11, /* 3: Electron */ 11,
426 /* 4: Neutrino */ 12, /* 5: mu+ */ -13,
427 /* 6: mu- */ 13, /* 7: pi0 */ 111,
428 /* 8: pi+ */ 211, /* 9: pi- */ -211,
429 /*10: K0 long */ 130, /*11: K+ */ 321,
430 /*12: K- */ -321, /*13: Neutron */ 2112,
431 /*14: Proton */ 2212, /*15: Antiproton*/ -2212,
432 /*16: K0 short */ 310, /*17: eta */ 221,
433 /*18: Lambda */ 3122, /*19: Sigma+ */ 3222,
434 /*20: Sigma0 */ 3212, /*21: Sigma- */ 3112,
435 /*22: Xi0 */ 3322, /*23: Xi- */ 3312,
436 /*24: Omega- */ 3334, /*25: Antineutrn*/ -2112,
437 /*26: Antilambda*/ -3122, /*27: Antisigma-*/ -3112,
438 /*28: Antisigma0*/ -3212, /*29: Antisigma+*/ -3222,
439 /*30: Antixi0 */ -3322, /*31: Antixi+ */ -3312,
440 /*32: Antiomega+*/ -3334, /*33: File */ 0,
441 /*34: Delta0 */ 2114, /*35: Delta++ */ 2224,
442 /*36: Delta+ */ 2214, /*37: Delta- */ 1114,
443 /*38: NP11+ */ 0, /*39: ND13+ */ 0,
444 /*40: NS11+ */ 0, /*41: rho0 */ 113,
445 /*42: rho+ */ 213, /*43: rho- */ -213,
446 /*44: NULL */ 0, /*45: Deuteron */ 0,
447 /*46: Tritium */ 0, /*47: Alpha */ 0,
448 /*48: NULL */ 0, /*49: He3 */ 0,
449 /*50: dimuon */ 0, /*51: dilepton */ 0,
450 /*52: omega */ 223, /*53: eta' */ 331,
451 /*54: sigma */ 0, /*55: phi */ 333,
452 /*56: Delta0*P33*/ 0, /*57: Delta++ *P33*/ 0,
453 /*58: Delta+*P33*/ 0, /*59: Delta- *P33 */ 0,
454 /*60: Delta0*S31*/ 0, /*61: Delta++ *S31*/ 0,
455 /*62: Delta+*S31*/ 0, /*63: Delta- *S31 */ 0,
456 /*64: NP110 */ 0, /*65: ND130 */ 0,
457 /*66: NS110 */ 0, /*67: J/Psi */ 443,
458 /*68: Psi' */ 0, /*69: pn */ 0};
459
460// "is Meson" flag
461const int PStdData::MESON[mnpar] = {
462 /* 0: dummy */ 0, /* 1: Photon */ 0,
463 /* 2: Positron */ 0, /* 3: Electron */ 0,
464 /* 4: Neutrino */ 0, /* 5: mu+ */ 0,
465 /* 6: mu- */ 0, /* 7: pi0 */ 1,
466 /* 8: pi+ */ 1, /* 9: pi- */ 1,
467 /*10: K0 long */ 1, /*11: K+ */ 1,
468 /*12: K- */ 1, /*13: Neutron */ 0,
469 /*14: Proton */ 0, /*15: Antiproton*/ 0,
470 /*16: K0 short */ 1, /*17: eta */ 1,
471 /*18: Lambda */ 0, /*19: Sigma+ */ 0,
472 /*20: Sigma0 */ 0, /*21: Sigma- */ 0,
473 /*22: Xi0 */ 0, /*23: Xi- */ 0,
474 /*24: Omega- */ 0, /*25: Antineutrn*/ 0,
475 /*26: Antilambda*/ 0, /*27: Antisigma-*/ 0,
476 /*28: Antisigma0*/ 0, /*29: Antisigma+*/ 0,
477 /*30: Antixi0 */ 0, /*31: Antixi+ */ 0,
478 /*32: Antiomega+*/ 0, /*33: File */ 0,
479 /*34: Delta0 */ 0, /*35: Delta++ */ 0,
480 /*36: Delta+ */ 0, /*37: Delta- */ 0,
481 /*38: NP11+ */ 0, /*39: ND13+ */ 0,
482 /*40: NS11+ */ 0, /*41: rho0 */ 1,
483 /*42: rho+ */ 1, /*43: rho- */ 1,
484 /*44: NULL */ 0, /*45: Deuteron */ 0,
485 /*46: Tritium */ 0, /*47: Alpha */ 0,
486 /*48: NULL */ 0, /*49: He3 */ 0,
487 /*50: dimuon */ 0, /*51: dilepton */ 0,
488 /*52: omega */ 1, /*53: eta' */ 1,
489 /*54: sigma */ 1, /*55: phi */ 1,
490 /*56: Delta0*P33*/ 0, /*57: Delta++ *P33*/ 0,
491 /*58: Delta+*P33*/ 0, /*59: Delta- *P33 */ 0,
492 /*60: Delta0*S31*/ 0, /*61: Delta++ *S31*/ 0,
493 /*62: Delta+*S31*/ 0, /*63: Delta- *S31 */ 0,
494 /*64: NP110 */ 0, /*65: ND130 */ 0,
495 /*66: NS110 */ 0, /*67: J/Psi */ 1,
496 /*68: Psi' */ 1, /*69: pn */ 0};
497
498// Baryon number
499const int PStdData::BARYON[mnpar] = {
500 /* 0: dummy */ 0, /* 1: Photon */ 0,
501 /* 2: Positron */ 0, /* 3: Electron */ 0,
502 /* 4: Neutrino */ 0, /* 5: mu+ */ 0,
503 /* 6: mu- */ 0, /* 7: pi0 */ 0,
504 /* 8: pi+ */ 0, /* 9: pi- */ 0,
505 /*10: K0 long */ 0, /*11: K+ */ 0,
506 /*12: K- */ 0, /*13: Neutron */ 1,
507 /*14: Proton */ 1, /*15: Antiproton*/ -1,
508 /*16: K0 short */ 0, /*17: eta */ 0,
509 /*18: Lambda */ 1, /*19: Sigma+ */ 1,
510 /*20: Sigma0 */ 1, /*21: Sigma- */ 1,
511 /*22: Xi0 */ 1, /*23: Xi- */ 1,
512 /*24: Omega- */ 1, /*25: Antineutrn*/ -1,
513 /*26: Antilambda*/ -1, /*27: Antisigma-*/ -1,
514 /*28: Antisigma0*/ -1, /*29: Antisigma+*/ -1,
515 /*30: Antixi0 */ -1, /*31: Antixi+ */ -1,
516 /*32: Antiomega+*/ -1, /*33: File */ 0,
517 /*34: Delta0 */ 1, /*35: Delta++ */ 1,
518 /*36: Delta+ */ 1, /*37: Delta- */ 1,
519 /*38: NP11+ */ 1, /*39: ND13+ */ 1,
520 /*40: NS11+ */ 1, /*41: rho0 */ 0,
521 /*42: rho+ */ 0, /*43: rho- */ 0,
522 /*44: NULL */ 0, /*45: Deuteron */ 2,
523 /*46: Tritium */ 3, /*47: Alpha */ 4,
524 /*48: NULL */ 0, /*49: He3 */ 3,
525 /*50: dimuon */ 0, /*51: dilepton */ 0,
526 /*52: omega */ 0, /*53: eta' */ 0,
527 /*54: sigma */ 0, /*55: phi */ 0,
528 /*56: Delta0*P33*/ 1, /*57: Delta++ *P33*/ 1,
529 /*58: Delta+*P33*/ 1, /*59: Delta- *P33 */ 1,
530 /*60: Delta0*S31*/ 1, /*61: Delta++ *S31*/ 1,
531 /*62: Delta+*S31*/ 1, /*63: Delta- *S31 */ 1,
532 /*64: NP110 */ 1, /*65: ND130 */ 1,
533 /*66: NS110 */ 1, /*67: J/Psi */ 0,
534 /*68: Psi' */ 0, /*69: pn */ 2};
535
536// Lepton number
537const int PStdData::LEPTON[mnpar] = {
538 /* 0: dummy */ 0, /* 1: Photon */ 0,
539 /* 2: Positron */ -1, /* 3: Electron */ 1,
540 /* 4: Neutrino */ 1, /* 5: mu+ */ -1,
541 /* 6: mu- */ 1, /* 7: pi0 */ 0,
542 /* 8: pi+ */ 0, /* 9: pi- */ 0,
543 /*10: K0 long */ 0, /*11: K+ */ 0,
544 /*12: K- */ 0, /*13: Neutron */ 0,
545 /*14: Proton */ 0, /*15: Antiproton*/ 0,
546 /*16: K0 short */ 0, /*17: eta */ 0,
547 /*18: Lambda */ 0, /*19: Sigma+ */ 0,
548 /*20: Sigma0 */ 0, /*21: Sigma- */ 0,
549 /*22: Xi0 */ 0, /*23: Xi- */ 0,
550 /*24: Omega- */ 0, /*25: Antineutrn*/ 0,
551 /*26: Antilambda*/ 0, /*27: Antisigma-*/ 0,
552 /*28: Antisigma0*/ 0, /*29: Antisigma+*/ 0,
553 /*30: Antixi0 */ 0, /*31: Antixi+ */ 0,
554 /*32: Antiomega+*/ 0, /*33: File */ 0,
555 /*34: Delta0 */ 0, /*35: Delta++ */ 0,
556 /*36: Delta+ */ 0, /*37: Delta- */ 0,
557 /*38: NP11+ */ 0, /*39: ND13+ */ 0,
558 /*40: NS11+ */ 0, /*41: rho0 */ 0,
559 /*42: rho+ */ 0, /*43: rho- */ 0,
560 /*44: NULL */ 0, /*45: Deuteron */ 0,
561 /*46: Tritium */ 0, /*47: Alpha */ 0,
562 /*48: NULL */ 0, /*49: He3 */ 0,
563 /*50: dimuon */ 0, /*51: dilepton */ 0,
564 /*52: omega */ 0, /*53: eta' */ 0,
565 /*54: sigma */ 0, /*55: phi */ 0,
566 /*56: Delta0*P33*/ 0, /*57: Delta++ *P33*/ 0,
567 /*58: Delta+*P33*/ 0, /*59: Delta- *P33 */ 0,
568 /*60: Delta0*S31*/ 0, /*61: Delta++ *S31*/ 0,
569 /*62: Delta+*S31*/ 0, /*63: Delta- *S31 */ 0,
570 /*64: NP110 */ 0, /*65: ND130 */ 0,
571 /*66: NS110 */ 0, /*67: J/Psi */ 0,
572 /*68: Psi' */ 0, /*69: pn */ 0};
573
574// Particle charge
575const int PStdData::CHARGE[mnpar] = {
576 /* 0: dummy */ 0, /* 1: Photon */ 0,
577 /* 2: Positron */ 1, /* 3: Electron */ -1,
578 /* 4: Neutrino */ 0, /* 5: mu+ */ 1,
579 /* 6: mu- */ -1, /* 7: pi0 */ 0,
580 /* 8: pi+ */ 1, /* 9: pi- */ -1,
581 /*10: K0 long */ 0, /*11: K+ */ 1,
582 /*12: K- */ -1, /*13: Neutron */ 0,
583 /*14: Proton */ 1, /*15: Antiproton*/ -1,
584 /*16: K0 short */ 0, /*17: eta */ 0,
585 /*18: Lambda */ 0, /*19: Sigma+ */ 1,
586 /*20: Sigma0 */ 0, /*21: Sigma- */ -1,
587 /*22: Xi0 */ 0, /*23: Xi- */ -1,
588 /*24: Omega- */ -1, /*25: Antineutrn*/ 0,
589 /*26: Antilambda*/ 0, /*27: Antisigma-*/ -1,
590 /*28: Antisigma0*/ 0, /*29: Antisigma+*/ 1,
591 /*30: Antixi0 */ 0, /*31: Antixi+ */ 1,
592 /*32: Antiomega+*/ 1, /*33: File */ 0,
593 /*34: Delta0 */ 0, /*35: Delta++ */ 2,
594 /*36: Delta+ */ 1, /*37: Delta- */ -1,
595 /*38: NP11+ */ 1, /*39: ND13+ */ 1,
596 /*40: NS11+ */ 1, /*41: rho0 */ 0,
597 /*42: rho+ */ 1, /*43: rho- */ -1,
598 /*44: NULL */ 0, /*45: Deuteron */ 1,
599 /*46: Tritium */ 1, /*47: Alpha */ 2,
600 /*48: NULL */ 0, /*49: He3 */ 2,
601 /*50: dimuon */ 0, /*51: dilepton */ 0,
602 /*52: omega */ 0, /*53: eta' */ 0,
603 /*54: sigma */ 0, /*55: phi */ 0,
604 /*56: Delta0*P33*/ 0, /*57: Delta++ *P33*/ 2,
605 /*58: Delta+*P33*/ 1, /*59: Delta- *P33 */ -1,
606 /*60: Delta0*S31*/ 0, /*61: Delta++ *S31*/ 2,
607 /*62: Delta+*S31*/ 1, /*63: Delta- *S31 */ -1,
608 /*64: NP110 */ 0, /*65: ND130 */ 0,
609 /*66: NS110 */ 0, /*67: J/Psi */ 0,
610 /*68: Psi' */ 0, /*69: pn */ 1};
611
612// 2 x angular momentum
613const int PStdData::SPIN[mnpar] = {
614 /* 0: dummy */ 0, /* 1: Photon */ 2,
615 /* 2: Positron */ 1, /* 3: Electron */ 1,
616 /* 4: Neutrino */ 1, /* 5: mu+ */ 1,
617 /* 6: mu- */ 1, /* 7: pi0 */ 0,
618 /* 8: pi+ */ 0, /* 9: pi- */ 0,
619 /*10: K0 long */ 0, /*11: K+ */ 0,
620 /*12: K- */ 0, /*13: Neutron */ 1,
621 /*14: Proton */ 1, /*15: Antiproton*/ 1,
622 /*16: K0 short */ 0, /*17: eta */ 0,
623 /*18: Lambda */ 1, /*19: Sigma+ */ 1,
624 /*20: Sigma0 */ 1, /*21: Sigma- */ 1,
625 /*22: Xi0 */ 1, /*23: Xi- */ 1,
626 /*24: Omega- */ 3, /*25: Antineutrn*/ 1,
627 /*26: Antilambda*/ 1, /*27: Antisigma-*/ 1,
628 /*28: Antisigma0*/ 1, /*29: Antisigma+*/ 1,
629 /*30: Antixi0 */ 1, /*31: Antixi+ */ 1,
630 /*32: Antiomega+*/ 3, /*33: File */ 0,
631 /*34: Delta0 */ 3, /*35: Delta++ */ 3,
632 /*36: Delta+ */ 3, /*37: Delta- */ 3,
633 /*38: NP11+ */ 1, /*39: ND13+ */ 3,
634 /*40: NS11+ */ 1, /*41: rho0 */ 2,
635 /*42: rho+ */ 2, /*43: rho- */ 2,
636 /*44: NULL */ 0, /*45: Deuteron */ 2,
637 /*46: Tritium */ 1, /*47: Alpha */ 0,
638 /*48: NULL */ 0, /*49: He3 */ 1,
639 /*50: dimuon */ 2, /*51: dilepton */ 2,
640 /*52: omega */ 2, /*53: eta' */ 0,
641 /*54: sigma */ 0, /*55: phi */ 2,
642 /*56: Delta0*P33*/ 3, /*57: Delta++ *P33*/ 3,
643 /*58: Delta+*P33*/ 3, /*59: Delta- *P33 */ 3,
644 /*60: Delta0*S31*/ 1, /*61: Delta++ *S31*/ 1,
645 /*62: Delta+*S31*/ 1, /*63: Delta- *S31 */ 1,
646 /*64: NP110 */ 1, /*65: ND130 */ 3,
647 /*66: NS110 */ 1, /*67: J/Psi */ 2,
648 /*68: Psi' */ 2, /*69: pn */ 0};
649
650// Parity (+/-1, 0 if irrelevant)
651const int PStdData::PARITY[mnpar] = {
652 /* 0: dummy */ 0, /* 1: Photon */ -1,
653 /* 2: Positron */ 0, /* 3: Electron */ 0,
654 /* 4: Neutrino */ 0, /* 5: mu+ */ 0,
655 /* 6: mu- */ 0, /* 7: pi0 */ -1,
656 /* 8: pi+ */ -1, /* 9: pi- */ -1,
657 /*10: K0 long */ -1, /*11: K+ */ -1,
658 /*12: K- */ -1, /*13: Neutron */ 1,
659 /*14: Proton */ 1, /*15: Antiproton*/ -1,
660 /*16: K0 short */ -1, /*17: eta */ -1,
661 /*18: Lambda */ 1, /*19: Sigma+ */ 1,
662 /*20: Sigma0 */ 1, /*21: Sigma- */ 1,
663 /*22: Xi0 */ 1, /*23: Xi- */ 1,
664 /*24: Omega- */ 1, /*25: Antineutrn*/ -1,
665 /*26: Antilambda*/ -1, /*27: Antisigma-*/ -1,
666 /*28: Antisigma0*/ -1, /*29: Antisigma+*/ -1,
667 /*30: Antixi0 */ -1, /*31: Antixi+ */ -1,
668 /*32: Antiomega+*/ -1, /*33: File */ 0,
669 /*34: Delta0 */ 1, /*35: Delta++ */ 1,
670 /*36: Delta+ */ 1, /*37: Delta- */ 1,
671 /*38: NP11+ */ 1, /*39: ND13+ */ -1,
672 /*40: NS11+ */ -1, /*41: rho0 */ -1,
673 /*42: rho+ */ -1, /*43: rho- */ -1,
674 /*44: NULL */ 0, /*45: Deuteron */ 0,
675 /*46: Tritium */ 0, /*47: Alpha */ 0,
676 /*48: NULL */ 0, /*49: He3 */ 0,
677 /*50: dimuon */ -1, /*51: dilepton */ -1,
678 /*52: omega */ -1, /*53: eta' */ -1,
679 /*54: sigma */ 1, /*55: phi */ -1,
680 /*56: Delta0*P33*/ 1, /*57: Delta++ *P33*/ 1,
681 /*58: Delta+*P33*/ 1, /*59: Delta- *P33 */ 1,
682 /*60: Delta0*S31*/ -1, /*61: Delta++ *S31*/ -1,
683 /*62: Delta+*S31*/ -1, /*63: Delta- *S31 */ -1,
684 /*64: NP110 */ 1, /*65: ND130 */ -1,
685 /*66: NS110 */ -1, /*67: J/Psi */ -1,
686 /*68: Psi' */ -1, /*69: pn */ 1};
687
688// 2 x isospin (also 0 if irrelevant)
689const int PStdData::ISPIN[mnpar] = {
690 /* 0: dummy */ 0, /* 1: Photon */ 0,
691 /* 2: Positron */ 0, /* 3: Electron */ 0,
692 /* 4: Neutrino */ 0, /* 5: mu+ */ 0,
693 /* 6: mu- */ 0, /* 7: pi0 */ 2,
694 /* 8: pi+ */ 2, /* 9: pi- */ 2,
695 /*10: K0 long */ 1, /*11: K+ */ 1,
696 /*12: K- */ 1, /*13: Neutron */ 1,
697 /*14: Proton */ 1, /*15: Antiproton*/ 1,
698 /*16: K0 short */ 1, /*17: eta */ 0,
699 /*18: Lambda */ 0, /*19: Sigma+ */ 2,
700 /*20: Sigma0 */ 2, /*21: Sigma- */ 2,
701 /*22: Xi0 */ 1, /*23: Xi- */ 1,
702 /*24: Omega- */ 0, /*25: Antineutrn*/ 0,
703 /*26: Antilambda*/ 0, /*27: Antisigma-*/ 2,
704 /*28: Antisigma0*/ 2, /*29: Antisigma+*/ 2,
705 /*30: Antixi0 */ 1, /*31: Antixi+ */ 1,
706 /*32: Antiomega+*/ 0, /*33: File */ 0,
707 /*34: Delta0 */ 3, /*35: Delta++ */ 3,
708 /*36: Delta+ */ 3, /*37: Delta- */ 3,
709 /*38: NP11+ */ 1, /*39: ND13+ */ 1,
710 /*40: NS11+ */ 1, /*41: rho0 */ 2,
711 /*42: rho+ */ 2, /*43: rho- */ 2,
712 /*44: NULL */ 0, /*45: Deuteron */ 0,
713 /*46: Tritium */ 0, /*47: Alpha */ 0,
714 /*48: NULL */ 0, /*49: He3 */ 0,
715 /*50: dimuon */ 0, /*51: dilepton */ 0,
716 /*52: omega */ 0, /*53: eta' */ 0,
717 /*54: sigma */ 0, /*55: phi */ 0,
718 /*56: Delta0*P33*/ 3, /*57: Delta++ *P33*/ 3,
719 /*58: Delta+*P33*/ 3, /*59: Delta- *P33 */ 3,
720 /*60: Delta0*S31*/ 3, /*61: Delta++ *S31*/ 3,
721 /*62: Delta+*S31*/ 3, /*63: Delta- *S31 */ 3,
722 /*64: NP110 */ 1, /*65: ND130 */ 1,
723 /*66: NS110 */ 1, /*67: J/Psi */ 0,
724 /*68: Psi' */ 0, /*69: pn */ 0};
725
726// Number of decay modes per particle
727const int PStdData::NMODES[mnpar] = {
728 /* 0: dummy */ 0, /* 1: Photon */ 0,
729 /* 2: Positron */ 0, /* 3: Electron */ 0,
730 /* 4: Neutrino */ 0, /* 5: mu+ */ 1,
731 /* 6: mu- */ 1, /* 7: pi0 */ 2,
732 /* 8: pi+ */ 1, /* 9: pi- */ 1,
733 /*10: K0 long */ 6, /*11: K+ */ 6,
734 /*12: K- */ 6, /*13: Neutron */ 1,
735 /*14: Proton */ 0, /*15: Antiproton*/ 0,
736 /*16: K0 short */ 3, /*17: eta */ 8,
737 /*18: Lambda */ 3, /*19: Sigma+ */ 0,
738 /*20: Sigma0 */ 0, /*21: Sigma- */ 0,
739 /*22: Xi0 */ 0, /*23: Xi- */ 0,
740 /*24: Omega- */ 0, /*25: Antineutrn*/ 1,
741 /*26: Antilambda*/ 0, /*27: Antisigma-*/ 0,
742 /*28: Antisigma0*/ 0, /*29: Antisigma+*/ 0,
743 /*30: Antixi0 */ 0, /*31: Antixi+ */ 0,
744 /*32: Antiomega+*/ 0, /*33: File */ 0,
745 /*34: Delta0 */ 4, /*35: Delta++ */ 1,
746 /*36: Delta+ */ 4, /*37: Delta- */ 1,
747 /*38: NP11+ */ 10, /*39: ND13+ */ 10,
748 /*40: NS11+ */ 13, /*41: rho0 */ 3,
749 /*42: rho+ */ 1, /*43: rho- */ 1,
750 /*44: NULL */ 0, /*45: Deuteron */ 0,
751 /*46: Tritium */ 0, /*47: Alpha */ 0,
752 /*48: NULL */ 0, /*49: He3 */ 0,
753 /*50: dimuon */ 1, /*51: dilepton */ 1,
754 /*52: omega */ 7, /*53: eta' */ 7,
755 /*54: sigma */ 3, /*55: phi */ 8,
756 /*56: Delta0*P33*/ 10, /*57: Delta++ *P33*/ 5,
757 /*58: Delta+*P33*/ 10, /*59: Delta- *P33 */ 5,
758 /*60: Delta0*S31*/ 8, /*61: Delta++ *S31*/ 4,
759 /*62: Delta+*S31*/ 8, /*63: Delta- *S31 */ 4,
760 /*64: NP110 */ 10, /*65: ND130 */ 10,
761 /*66: NS110 */ 13, /*67: J/Psi */ 8,
762 /*68: Psi' */ 8, /*69: pn */ 2};
763
764// Static branching ratios:
765// different charge states are included explicitly
766// with the appropriate isospin factors.
767// Note: The static branching ratios for decay modes of
768// the type 1) hadron -> hadron + hadron and
769// 2) hadron -> N + pi + pi are averages of PDG data, and
770// are only included for completeness. They are actually
771// calculated as functions of mass explicitly (see Width()
772// and Width1()). Static BR's are only used for modes that
773// are not entirely hadronic, or involve more than two
774// decay products.
775//const double BRR[mnmodes]={
776double PStdData::BRR[mnmodes] = {
777 // mu+, 1 channel
778 1., // id=5 mu+ --> e+ + neutrino + neutrino
779 // mu-, 1 channel
780 1., // id=6 mu- --> e- + neutrino + neutrino
781 // pi0, 2 channels
782 0.988, // id=7 pi0 --> photon + photon
783 0.012, // id=7 pi0 --> photon + dilepton (Dalitz)
784 // pi+, 1 channel
785 1., // id=8 pi+ --> mu+ + neutrino
786 // pi-, 1 channel
787 1., // id=9 pi- --> mu- + neutrino
788 // K0L, 6 channels
789 0.211, // id=10 K0 long --> pi0 + pi0 + pi0
790 0.126, // id=10 K0 long --> pi+ + pi- + pi0
791 0.136, // id=10 K0 long --> pi+ + mu- + neutrino
792 0.136, // id=10 K0 long --> pi- + mu+ + neutrino
793 0.194, // id=10 K0 long --> pi+ + e- + neutrino
794 0.194, // id=10 K0 long --> pi- + e+ + neutrino
795 // K+, 6 channels
796 0.635, // id=11 K+ --> mu+ + neutrino
797 0.212, // id=11 K+ --> pi+ + pi0
798 0.056, // id=11 K+ --> pi+ + pi+ + pi-
799 0.017, // id=11 K+ --> pi+ + pi0 + pi0
800 0.032, // id=11 K+ --> pi0 + mu+ + neutrino
801 0.048, // id=11 K+ --> pi0 + e+ + neutrino
802 // K-, 6 channels
803 0.635, // id=12 K- --> mu- + neutrino
804 0.212, // id=12 K- --> pi- + pi0
805 0.056, // id=12 K- --> pi- + pi- + pi+
806 0.017, // id=12 K- --> pi- + pi0 + pi0
807 0.032, // id=12 K- --> pi0 + mu- + neutrino
808 0.048, // id=12 K- --> pi0 + e- + neutrino
809 // n, 1 channel
810 1., // id=13 n --> p + e- + neutrino
811 // K0S, 3 channels
812 0.6851, // id=16 K0 short --> pi+ + pi-
813 0.3129, // id=16 K0 short --> pi0 + pi0
814 0.0018, // id=16 K0 short --> pi+ + pi- + photon
815 // eta, 8 channels
816 0.394, // id=17 eta --> photon + photon
817 0.325, // id=17 eta --> pi0 + pi0 + pi0
818 0.226, // id=17 eta --> pi+ + pi- + pi0
819 0.0468, // id=17 eta --> pi+ + pi- + photon
820 0.006, // id=17 eta --> photon + dilepton (Dalitz)
821 0.00031, // id=17 eta --> photon + dimuon
822 6.0e-6, // id=17 eta --> e+ + e-
823 5.8e-6, // id=17 eta --> mu+ + mu-
824 // Lambda, 3 channels
825 0.639, // id=18 Lambda --> p + pi-
826 0.358, // id=18 Lambda --> n + pi0
827 0.0018, // id=18 Lambda --> n + photon
828 // anti_n, 1 channel
829 1., // id=25 anti_n --> anti_p + e+ + neutrino
830 // D0, 4 channels
831 0.99 / 3., // id=34 Delta0 --> p + pi-
832 0.99 * 2. / 3., // id=34 Delta0 --> n + pi0
833 0.0055, // id=34 Delta0 --> n + photon
834 0.0055 / 137., // id=34 Delta0 --> n + dilepton
835 // D++, 1 channel
836 1., // id=35 Delta++ --> p + pi+
837 // D+, 4 channels
838 0.99 * 2. / 3, // id=36 Delta+ --> p + pi0
839 0.99 / 3., // id=36 Delta+ --> n + pi+
840 0.0055, // id=36 Delta+ --> p + photon
841 0.0055 / 137., // id=36 Delta+ --> p + dilepton
842 // D-, 1 channel
843 1., // id=36 Delta- --> n + pi-
844 // NP11+, 10 channels
845 0.65 / 3., // id=38 N*(1440)+ --> p + pi0
846 0.65 * 2. / 3., // id=38 N*(1440)+ --> n + pi+
847 0.25 / 2., // id=38 N*(1440)+ --> Delta++ + pi-
848 0.25 / 3., // id=38 N*(1440)+ --> Delta+ + pi0
849 0.25 / 6., // id=38 N*(1440)+ --> Delta0 + pi+
850 0.024585 / 3., // id=38 N*(1440)+ --> p + rho0
851 0.024585 * 2 / 3., // id=38 N*(1440)+ --> n + rho+
852 0.075 / 3., // id=38 N*(1440)+ --> p + pi0 + pi0
853 0.075 * 2. / 3., // id=38 N*(1440)+ --> p + pi+ + pi-
854 0.000415, // id=38 N*(1440)+ --> p + photon
855 // ND13+, 10 channels
856 0.55 / 3., // id=39 N*(1520)+ --> p + pi0
857 0.55 * 2. / 3., // id=39 N*(1520)+ --> n + pi+
858 0.20 / 2., // id=39 N*(1520)+ --> Delta++ + pi-
859 0.20 / 3., // id=39 N*(1520)+ --> Delta+ + pi0
860 0.20 / 6., // id=39 N*(1520)+ --> Delta0 + pi+
861 0.20 / 3., // id=39 N*(1520)+ --> p + rho0
862 0.20 * 2 / 3., // id=39 N*(1520)+ --> n + rho+
863 0.0449 / 3., // id=39 N*(1520)+ --> p + pi0 + pi0
864 0.0449 * 2. / 3., // id=39 N*(1520)+ --> p + pi+ + pi-
865 0.0051, // id=39 N*(1520)+ --> p + photon
866 // NS11+, 13 channels
867 0.46 / 3., // id=40 N*(1535)+ --> p + pi0
868 0.46 * 2. / 3., // id=40 N*(1535)+ --> n + pi+
869 0.3875, // id=40 N*(1535)+ --> p + eta
870 0.01 / 2., // id=40 N*(1535)+ --> Delta++ + pi-
871 0.01 / 3., // id=40 N*(1535)+ --> Delta+ + pi0
872 0.01 / 6., // id=40 N*(1535)+ --> Delta0 + pi+
873 0.04 / 3., // id=40 N*(1535)+ --> p + rho0
874 0.04 * 2 / 3., // id=40 N*(1535)+ --> n + rho+
875 0.03 / 3., // id=40 N*(1535)+ --> p + pi0 + pi0
876 0.03 * 2. / 3., // id=40 N*(1535)+ --> p + pi+ + pi-
877 0.07 / 3., // id=40 N*(1535)+ --> N*(1440)+ + pi0
878 0.07 * 2. / 3., // id=40 N*(1535)+ --> N*(1440)0 + pi+
879 0.0025, // id=40 N*(1535)+ --> p + photon
880 // rho0, 3 channels
881 0.9999091, // id=41 rho0 --> pi+ + pi-
882 4.67e-5, // id=41 rho0 --> e+ + e-
883 4.55e-5, // id=41 rho0 --> mu+ + mu-
884 // rho+, 1 channel
885 1., // id=42 rho+ --> pi+ + pi0
886 // rho-, 1 channel
887 1., // id=43 rho- --> pi- + pi0
888 // dimuon, 1 channel
889 1., // id=50 dimuon --> mu+ + mu-
890 // dilepton, 1 channel
891 1., // id=51 dilepton --> e+ + e-
892 // w, 7 channels
893 0.888, // id=52 omega --> pi+ + pi- + pi0
894 0.085, // id=52 omega --> pi0 + photon
895 0.017, // id=52 omega --> pi+ + pi-
896 5.9e-4, // id=52 omega --> pi0 + dilepton (Dalitz)
897 9.6e-5, // id=52 omega --> pi0 + dimuon
898 7.07e-5, // id=52 omega --> e+ + e-
899 8.0e-5, // id=52 omega --> mu+ + mu-
900 // eta', 7 channels
901 0.443, // id=53 eta' --> pi- + pi+ + eta
902 0.295, // id=53 eta' --> rho0 + photon
903 0.209, // id=53 eta' --> pi0 + pi0 + eta
904 0.0303, // id=53 eta' --> omega + photon
905 0.0212, // id=53 eta' --> photon + photon
906 0.00156, // id=53 eta' --> pi0 + pi0 + pi0
907 0.000104, // id=53 eta' --> photon + dimuon
908 // sigma, 3 channels
909 0.99999, // id=54 sigma --> pi+ + pi-
910 0.000005, // id=54 sigma --> e+ + e-
911 0.000005, // id=54 sigma --> mu+ + mu-
912 // phi, 8 channels
913 0.492, // id=55 phi --> K+ + K-
914 0.338, // id=55 phi --> K0L + K0S
915 0.155, // id=55 phi --> pi+ + pi- + pi0
916 0.01297, // id=55 phi --> eta + photon
917 0.00126, // id=55 phi --> pi0 + photon
918 0.000296, // id=55 phi --> e+ + e-
919 0.000287, // id=55 phi --> mu+ + mu-
920 0.000115, // id=55 phi --> dilepton + eta
921 // DP330, 10 channels
922 .175 / 3., // id=56 Delta(1600)0 --> p + pi-
923 .175 * 2. / 3., // id=56 Delta(1600)0 --> n + pi0
924 .55 * 8. / 15., // id=56 Delta(1600)0 --> Delta+ + pi-
925 .55 / 15., // id=56 Delta(1600)0 --> Delta0 + pi0
926 .55 * 2. / 5., // id=56 Delta(1600)0 --> Delta- + pi+
927 .225 / 3., // id=56 Delta(1600)0 --> p + rho-
928 .225 * 2. / 3., // id=56 Delta(1600)0 --> n + rho0
929 .0499 / 3., // id=56 Delta(1600)0 --> N(1440)+ + rho-
930 .0499 * 2. / 3., // id=56 Delta(1600)0 --> N(1440)0 + rho0
931 .0001, // id=56 Delta(1600)0 --> n + photon
932 // DP33++, 5 channels
933 .175, // id=57 Delta(1600)++ --> p + pi+
934 .55 * 3. / 5., // id=57 Delta(1600)++ --> Delta++ + pi0
935 .55 * 2. / 5., // id=57 Delta(1600)++ --> Delta+ + pi+
936 .225, // id=57 Delta(1600)++ --> p + rho+
937 .05, // id=57 Delta(1600)++ --> N(1440)+ + rho+
938 // DP33+, 10 channels
939 .175 * 2. / 3., // id=58 Delta(1600)+ --> p + pi0
940 .175 / 3., // id=58 Delta(1600)+ --> n + pi+
941 .55 * 2. / 5., // id=58 Delta(1600)+ --> Delta++ + pi-
942 .55 / 15., // id=58 Delta(1600)+ --> Delta+ + pi0
943 .55 * 8. / 15., // id=58 Delta(1600)+ --> Delta0 + pi+
944 .225 * 2. / 3., // id=58 Delta(1600)+ --> p + rho0
945 .225 / 3., // id=58 Delta(1600)+ --> n + rho+
946 .0499 * 2. / 3., // id=58 Delta(1600)+ --> N(1440)+ + rho0
947 .0499 / 3., // id=58 Delta(1600)+ --> N(1440)0 + rho+
948 .0001, // id=58 Delta(1600)+ --> p + photon
949 // DP33-, 5 channels
950 .175, // id=59 Delta(1600)- --> n + pi-
951 .55 * 2. / 5., // id=59 Delta(1600)- --> Delta0 + pi-
952 .55 * 3. / 5., // id=59 Delta(1600)- --> Delta- + pi0
953 .225, // id=59 Delta(1600)- --> n + rho-
954 .05, // id=59 Delta(1600)- --> N(1440)0 + rho-
955 // DS310, 8 channels
956 .25 / 3., // id=60 Delta(1620)0 --> p + pi-
957 .25 * 2. / 3., // id=60 Delta(1620)0 --> n + pi0
958 .5897 * 8. / 15., // id=60 Delta(1620)0 --> Delta+ + pi-
959 .5897 / 15., // id=60 Delta(1620)0 --> Delta0 + pi0
960 .5897 * 2. / 5., // id=60 Delta(1620)0 --> Delta- + pi+
961 .16 / 3., // id=60 Delta(1620)0 --> p + rho-
962 .16 * 2. / 3., // id=60 Delta(1620)0 --> n + rho0
963 .0003, // id=60 Delta(1620)0 --> n + photon
964 // DS31++, 4 channels
965 .25, // id=61 Delta(1620)++ --> p + pi+
966 .5897 * 3. / 5., // id=61 Delta(1620)++ --> Delta++ + pi0
967 .5897 * 2. / 5., // id=61 Delta(1620)++ --> Delta+ + pi+
968 .16, // id=61 Delta(1620)++ --> p + rho+
969 // DS31+, 8 channels
970 .25 * 2. / 3., // id=62 Delta(1620)+ --> p + pi0
971 .25 / 3., // id=62 Delta(1620)+ --> n + pi+
972 .5897 * 2. / 5., // id=62 Delta(1620)+ --> Delta++ + pi-
973 .5897 / 15., // id=62 Delta(1620)+ --> Delta+ + pi0
974 .5897 * 8. / 15., // id=62 Delta(1620)+ --> Delta0 + pi+
975 .16 * 2. / 3., // id=62 Delta(1620)+ --> p + rho0
976 .16 / 3., // id=62 Delta(1620)+ --> n + rho+
977 .0003, // id=62 Delta(1620)+ --> p + photon
978 // DS31-, 4 channels
979 .25, // id=63 Delta(1620)- --> n + pi-
980 .5897 * 2. / 5., // id=63 Delta(1620)- --> Delta0 + pi-
981 .5897 * 3. / 5., // id=63 Delta(1620)- --> Delta- + pi0
982 .16, // id=63 Delta(1620)- --> n + rho-
983 // NP110, 10 channels
984 0.65 / 3., // id=64 N*(1440)0 --> p + pi-
985 0.65 * 2. / 3., // id=64 N*(1440)0 --> n + pi0
986 0.25 / 2., // id=64 N*(1440)0 --> Delta+ + pi-
987 0.25 / 3., // id=64 N*(1440)0 --> Delta0 + pi0
988 0.25 / 6., // id=64 N*(1440)0 --> Delta- + pi+
989 0.024585 / 3., // id=64 N*(1440)0 --> p + rho-
990 0.024585 * 2 / 3., // id=64 N*(1440)0 --> n + rho0
991 0.075 * 2. / 3., // id=64 N*(1440)0 --> n + pi+ + pi-
992 0.075 / 3., // id=64 N*(1440)0 --> n + pi0 + pi0
993 0.000415, // id=64 N*(1440)0 --> n + photon
994 // ND130, 10 channels
995 0.55 * 2. / 3., // id=65 N*(1520)0 --> p + pi-
996 0.55 / 3., // id=65 N*(1520)0 --> n + pi0
997 0.20 / 6., // id=65 N*(1520)0 --> Delta+ + pi-
998 0.20 / 3., // id=65 N*(1520)0 --> Delta0 + pi0
999 0.20 / 2., // id=65 N*(1520)0 --> Delta- + pi+
1000 0.20 * 2. / 3., // id=65 N*(1520)0 --> p + rho-
1001 0.20 / 3., // id=65 N*(1520)0 --> n + rho0
1002 0.0449 * 2. / 3., // id=65 N*(1520)0 --> n + pi+ + pi-
1003 0.0449 / 3., // id=65 N*(1520)0 --> n + pi0 + pi0
1004 0.0051, // id=65 N*(1520)0 --> n + photon
1005 // NS110, 13 channels
1006 0.46 / 3., // id=66 N*(1535)0 --> p + pi-
1007 0.46 * 2. / 3., // id=66 N*(1535)0 --> n + pi0
1008 0.3875, // id=66 N*(1535)0 --> n + eta
1009 0.01 / 2., // id=66 N*(1535)0 --> Delta+ + pi-
1010 0.01 / 3., // id=66 N*(1535)0 --> Delta0 + pi0
1011 0.01 / 6., // id=66 N*(1535)0 --> Delta- + pi+
1012 0.04 / 3., // id=66 N*(1535)0 --> p + rho-
1013 0.04 * 2 / 3., // id=66 N*(1535)0 --> n + rho0
1014 0.03 * 2. / 3., // id=66 N*(1535)0 --> n + pi+ + pi-
1015 0.03 / 3., // id=66 N*(1535)0 --> n + pi0 + pi0
1016 0.07 / 3., // id=66 N*(1535)0 --> N*(1440)+ + pi-
1017 0.07 * 2. / 3., // id=66 N*(1535)0 --> N*(1440)0 + pi0
1018 0.0025, // id=66 N*(1535)0 --> n + photon
1019 // J/Psi, 8 channels
1020 0.0602, // id=67 J/Psi --> e+ + e-
1021 0.0601, // id=67 J/Psi --> mu+ + mu-
1022 0.0088, // id=67 J/Psi --> dilepton + photon
1023 0.0337, // id=67 J/Psi --> 2pi+ + 2pi- + pi0
1024 0.0290, // id=67 J/Psi --> 3pi+ + 3pi- + pi0
1025 0.0150, // id=67 J/Psi --> pi+ + pi- + pi0
1026 0.0120, // id=67 J/Psi --> K+ + K- + pi+ + pi- + pi0
1027 0.7812, // id=67 J/Psi --> junk
1028 // Psi' = Psi(2S), 8 channels
1029 0.0073, // id=68 Psi' --> e+e-
1030 0.0070, // id=68 Psi' --> mu+mu-
1031 0.31, // id=68 Psi' --> J/Psi pi+ pi-
1032 0.182, // id=68 Psi' --> J/Psi pi0 pi0
1033 0.027, // id=68 Psi' --> J/Psi eta
1034 0.0035, // id=68 Psi' --> 2pi+ + 2pi- + pi0
1035 0.0030, // id=68 Psi' --> 3pi+ + 3pi- + pi0
1036 0.4554, // id=68 Psi' --> junk
1037 // pn, 2 channels
1038 0.999, // id=69 pn --> dilepton + p + n
1039 0.001 // id=69 pn --> photon + p + n
1040};
1041
1042// Decay-mode nomenclature:
1043// id1 + id2*1000 + .. + idn*1000^(n-1), or id of non-virtual-photon if Dalitz
1044const char* PStdData::MODE[mnmodes] = {
1045 // mu+, 1 channel
1046 "4004002", // id=5 mu+ --> e+ + neutrino + neutrino
1047 // mu-, 1 channel
1048 "4004003", // id=6 mu- --> e- + neutrino + neutrino
1049 // pi0, 2 channels
1050 "1001", // id=7 pi0 --> photon + photon
1051 "1051", // id=7 pi0 --> dilepton + photon (Dalitz)
1052 // pi+, 1 channel
1053 "4005", // id=8 pi+ --> mu+ + neutrino
1054 // pi-, 1 channel
1055 "4006", // id=9 pi- --> mu- + neutrino
1056 // K0L, 6 channels
1057 "7007007", // id=10 K0 long --> pi0 + pi0 + pi0
1058 "7009008", // id=10 K0 long --> pi+ + pi- + pi0
1059 "4006008", // id=10 K0 long --> pi+ + mu- + neutrino
1060 "4005009", // id=10 K0 long --> pi- + mu+ + neutrino
1061 "4003008", // id=10 K0 long --> pi+ + e- + neutrino
1062 "4002009", // id=10 K0 long --> pi- + e+ + neutrino
1063 // K+, 6 channels
1064 "4005", // id=11 K+ --> mu+ + neutrino
1065 "8007", // id=11 K+ --> pi+ + pi0
1066 "9008008", // id=11 K+ --> pi+ + pi+ + pi-
1067 "7007008", // id=11 K+ --> pi+ + pi0 + pi0
1068 "4005007", // id=11 K+ --> pi0 + mu+ + neutrino
1069 "4002007", // id=11 K+ --> pi0 + e+ + neutrino
1070 // K-, 6 channels
1071 "4006", // id=12 K- --> mu- + neutrino
1072 "7009", // id=12 K- --> pi- + pi0
1073 "8009009", // id=12 K- --> pi- + pi- + pi+
1074 "7007009", // id=12 K- --> pi- + pi0 + pi0
1075 "4006007", // id=12 K- --> pi0 + mu- + neutrino
1076 "4003007", // id=12 K- --> pi0 + e- + neutrino
1077 // n, 1 channel
1078 "4003014", // id=13 n --> p + e- + neutrino
1079 // K0S, 3 channels
1080 "9008", // id=16 K0 short --> pi+ + pi-
1081 "7007", // id=16 K0 short --> pi0 + pi0
1082 "1009008", // id=16 K0 short --> pi+ + pi- + photon
1083 // eta, 8 channels
1084 "1001", // id=17 eta --> photon + photon
1085 "7007007", // id=17 eta --> pi0 + pi0 + pi0
1086 "8009007", // id=17 eta --> pi+ + pi- + pi0
1087 "8009001", // id=17 eta --> pi+ + pi- + photon
1088 "1051", // id=17 eta --> dilepton + photon (Dalitz)
1089 "1050", // id=17 eta --> photon + dimuon
1090 "3002", // id=17 eta --> e+ + e-
1091 "5006", // id=17 eta --> mu+ + mu-
1092 // Lambda, 3 channels
1093 "9014", // id=18 Lambda --> p + pi-
1094 "7013", // id=18 Lambda --> n + pi0
1095 "1013", // id=18 Lambda --> n + photon
1096 // anti_n, 1 channel
1097 "4002015", // id=25 anti_n --> anti_p + e+ + neutrino
1098 // D0, 4 channels
1099 "9014", // id=34 Delta0 --> p + pi-
1100 "7013", // id=34 Delta0 --> n + pi0
1101 "1013", // id=34 Delta0 --> n + photon
1102 "13051", // id=34 Delta0 --> dilepton + n (Dalitz)
1103 // D++, 1 channel
1104 "8014", // id=35 Delta++ --> p + pi+
1105 // D+, 4 channels
1106 "7014", // id=36 Delta+ --> p + pi0
1107 "8013", // id=36 Delta+ --> n + pi+
1108 "1014", // id=36 Delta+ --> p + photon
1109 "14051", // id=36 Delta+ --> dilepton + p (Dalitz)
1110 // D-, 1 channel
1111 "9013", // id=36 Delta- --> n + pi-
1112 // NP11+, 10 channels
1113 "7014", // id=38 N*(1440)+ --> p + pi0
1114 "8013", // id=38 N*(1440)+ --> n + pi+
1115 "9035", // id=38 N*(1440)+ --> Delta++ + pi-
1116 "7036", // id=38 N*(1440)+ --> Delta+ + pi0
1117 "8034", // id=38 N*(1440)+ --> Delta0 + pi+
1118 "41014", // id=38 N*(1440)+ --> p + rho0
1119 "42013", // id=38 N*(1440)+ --> n + rho+
1120 "7007014", // id=38 N*(1440)+ --> p + pi0 + pi0
1121 "9008014", // id=38 N*(1440)+ --> p + pi+ + pi-
1122 "1014", // id=38 N*(1440)+ --> p + photon
1123 // ND13+, 10 channels
1124 "7014", // id=39 N*(1520)+ --> p + pi0
1125 "8013", // id=39 N*(1520)+ --> n + pi+
1126 "9035", // id=39 N*(1520)+ --> Delta++ + pi-
1127 "7036", // id=39 N*(1520)+ --> Delta+ + pi0
1128 "8034", // id=39 N*(1520)+ --> Delta0 + pi+
1129 "41014", // id=39 N*(1520)+ --> p + rho0
1130 "42013", // id=39 N*(1520)+ --> n + rho+
1131 "7007014", // id=39 N*(1520)+ --> p + pi0 + pi0
1132 "9008014", // id=39 N*(1520)+ --> p + pi+ + pi-
1133 "1014", // id=39 N*(1520)+ --> p + photon
1134 // NS11+, 13 channels
1135 "7014", // id=40 N*(1535)+ --> p + pi0
1136 "8013", // id=40 N*(1535)+ --> n + pi+
1137 "17014", // id=40 N*(1535)+ --> p + eta
1138 "9035", // id=40 N*(1535)+ --> Delta++ + pi-
1139 "7036", // id=40 N*(1535)+ --> Delta+ + pi0
1140 "8034", // id=40 N*(1535)+ --> Delta0 + pi+
1141 "41014", // id=40 N*(1535)+ --> p + rho0
1142 "42013", // id=40 N*(1535)+ --> n + rho+
1143 "7007014", // id=40 N*(1535)+ --> p + pi0 + pi0
1144 "9008014", // id=40 N*(1535)+ --> p + pi+ + pi-
1145 "7038", // id=40 N*(1535)+ --> N*(1440)+ + pi0
1146 "8064", // id=40 N*(1535)+ --> N*(1440)0 + pi+
1147 "1014", // id=40 N*(1535)+ --> p + photon
1148 // rho0, 3 channels
1149 "9008", // id=41 rho0 --> pi+ + pi-
1150 "3002", // id=41 rho0 --> e+ + e-
1151 "5006", // id=41 rho0 --> mu+ + mu-
1152 // rho+, 1 channel
1153 "7008", // id=42 rho+ --> pi+ + pi0
1154 // rho-, 1 channel
1155 "7009", // id=43 rho- --> pi- + pi0
1156 // dimuon, 1 channel
1157 "6005", // id=50 dimuon --> mu+ + mu-
1158 // dilepton, 1 channel
1159 "3002", // id=51 dilepton --> e+ + e-
1160 // w, 7 channels
1161 "7009008", // id=52 omega --> pi+ + pi- + pi0
1162 "1007", // id=52 omega --> pi0 + photon
1163 "9008", // id=52 omega --> pi+ + pi-
1164 "7051", // id=52 omega --> dilepton + pi0 (Dalitz)
1165 "7050", // id=52 omega --> dimuon + pi0
1166 "3002", // id=52 omega --> e+ + e-
1167 "5006", // id=52 omega --> mu+ + mu-
1168 // eta', 7 channels
1169 "9008017", // id=53 eta' --> eta + pi- + pi+
1170 "1041", // id=53 eta' --> rho0 + photon
1171 "7007017", // id=53 eta' --> eta + pi0 + pi0
1172 "1052", // id=53 eta' --> omega + photon
1173 "1001", // id=53 eta' --> photon + photon
1174 "7007007", // id=53 eta' --> pi0 + pi0 + pi0
1175 "1050", // id=53 eta' --> dimuon + photon
1176 // sigma, 3 channels
1177 "9008", // id=54 sigma --> pi+ + pi-
1178 "3002", // id=54 sigma --> e+ + e-
1179 "5006", // id=54 sigma --> mu+ + mu-
1180 // phi, 8 channels
1181 "12011", // id=55 phi --> K+ + K-
1182 "16010", // id=55 phi --> K0L + K0S
1183 "7009008", // id=55 phi --> pi+ + pi- + pi0
1184 "1017", // id=55 phi --> eta + photon
1185 "1007", // id=55 phi --> pi0 + photon
1186 "3002", // id=55 phi --> e+ + e-
1187 "5006", // id=55 phi --> mu+ + mu-
1188 "17051", // id=55 phi --> dilepton + eta
1189 // DP330, 10 channels
1190 "9014", // id=56 Delta(1600)0 --> p + pi-
1191 "7013", // id=56 Delta(1600)0 --> n + pi0
1192 "9036", // id=56 Delta(1600)0 --> Delta+ + pi-
1193 "7034", // id=56 Delta(1600)0 --> Delta0 + pi0
1194 "8037", // id=56 Delta(1600)0 --> Delta- + pi+
1195 "43014", // id=56 Delta(1600)0 --> p + rho-
1196 "41013", // id=56 Delta(1600)0 --> n + rho0
1197 "43038", // id=56 Delta(1600)0 --> N(1440)+ + rho-
1198 "41064", // id=56 Delta(1600)0 --> N(1440)0 + rho0
1199 "1013", // id=56 Delta(1600)0 --> n + photon
1200 // DP33++, 5 channels
1201 "8014", // id=57 Delta(1600)++ --> p + pi+
1202 "7035", // id=57 Delta(1600)++ --> Delta++ + pi0
1203 "8036", // id=57 Delta(1600)++ --> Delta+ + pi+
1204 "42014", // id=57 Delta(1600)++ --> p + rho+
1205 "42038", // id=57 Delta(1600)++ --> N(1440)+ + rho+
1206 // DP33+, 10 channels
1207 "7014", // id=58 Delta(1600)+ --> p + pi0
1208 "8013", // id=58 Delta(1600)+ --> n + pi+
1209 "9035", // id=58 Delta(1600)+ --> Delta++ + pi-
1210 "7036", // id=58 Delta(1600)+ --> Delta+ + pi0
1211 "8034", // id=58 Delta(1600)+ --> Delta0 + pi+
1212 "41014", // id=58 Delta(1600)+ --> p + rho0
1213 "42013", // id=58 Delta(1600)+ --> n + rho+
1214 "41038", // id=58 Delta(1600)+ --> N(1440)+ + rho0
1215 "42064", // id=58 Delta(1600)+ --> N(1440)0 + rho+
1216 "1014", // id=58 Delta(1600)+ --> p + photon
1217 // DP33-, 5 channels
1218 "9013", // id=59 Delta(1600)- --> n + pi-
1219 "9034", // id=59 Delta(1600)- --> Delta0 + pi-
1220 "7037", // id=59 Delta(1600)- --> Delta- + pi0
1221 "43013", // id=59 Delta(1600)- --> n + rho-
1222 "43064", // id=59 Delta(1600)- --> N(1440)0 + rho-
1223 // DS310, 8 channels
1224 "9014", // id=60 Delta(1620)0 --> p + pi-
1225 "7013", // id=60 Delta(1620)0 --> n + pi0
1226 "9036", // id=60 Delta(1620)0 --> Delta+ + pi-
1227 "7034", // id=60 Delta(1620)0 --> Delta0 + pi0
1228 "8037", // id=60 Delta(1620)0 --> Delta- + pi+
1229 "43014", // id=60 Delta(1620)0 --> p + rho-
1230 "41013", // id=60 Delta(1620)0 --> n + rho0
1231 "1013", // id=60 Delta(1620)0 --> n + photon
1232 // DS31++, 4 channels
1233 "8014", // id=61 Delta(1620)++ --> p + pi+
1234 "7035", // id=61 Delta(1620)++ --> Delta++ + pi0
1235 "8036", // id=61 Delta(1620)++ --> Delta+ + pi+
1236 "42014", // id=61 Delta(1620)++ --> p + rho+
1237 // DS31+, 8 channels
1238 "7014", // id=62 Delta(1620)+ --> p + pi0
1239 "8013", // id=62 Delta(1620)+ --> n + pi+
1240 "9035", // id=62 Delta(1620)+ --> Delta++ + pi-
1241 "7036", // id=62 Delta(1620)+ --> Delta+ + pi0
1242 "8034", // id=62 Delta(1620)+ --> Delta0 + pi+
1243 "41014", // id=62 Delta(1620)+ --> p + rho0
1244 "42013", // id=62 Delta(1620)+ --> n + rho+
1245 "1014", // id=62 Delta(1620)+ --> p + photon
1246 // DS31-, 4 channels
1247 "9013", // id=63 Delta(1620)- --> n + pi-
1248 "9034", // id=63 Delta(1620)- --> Delta0 + pi-
1249 "7037", // id=63 Delta(1620)- --> Delta- + pi0
1250 "43013", // id=63 Delta(1620)- --> n + rho-
1251 // NP110, 10 channels
1252 "9014", // id=64 N*(1440)0 --> p + pi-
1253 "7013", // id=64 N*(1440)0 --> n + pi0
1254 "9036", // id=64 N*(1440)0 --> Delta+ + pi-
1255 "7034", // id=64 N*(1440)0 --> Delta0 + pi0
1256 "8037", // id=64 N*(1440)0 --> Delta- + pi+
1257 "43014", // id=64 N*(1440)0 --> p + rho-
1258 "41013", // id=64 N*(1440)0 --> n + rho0
1259 "9008013", // id=64 N*(1440)0 --> n + pi+ + pi-
1260 "7007013", // id=64 N*(1440)0 --> n + pi0 + pi0
1261 "1013", // id=64 N*(1440)0 --> n + photon
1262 // ND130, 10 channels
1263 "9014", // id=65 N*(1520)0 --> p + pi-
1264 "7013", // id=65 N*(1520)0 --> n + pi0
1265 "9036", // id=65 N*(1520)0 --> Delta+ + pi-
1266 "7034", // id=65 N*(1520)0 --> Delta0 + pi0
1267 "8037", // id=65 N*(1520)0 --> Delta- + pi+
1268 "43014", // id=65 N*(1520)0 --> p + rho-
1269 "41013", // id=65 N*(1520)0 --> n + rho0
1270 "9008013", // id=65 N*(1520)0 --> n + pi+ + pi-
1271 "7007013", // id=65 N*(1520)0 --> n + pi0 + pi0
1272 "1013", // id=65 N*(1520)0 --> n + photon
1273 // NS110, 13 channels
1274 "9014", // id=66 N*(1535)0 --> p + pi-
1275 "7013", // id=66 N*(1535)0 --> n + pi0
1276 "17013", // id=66 N*(1535)0 --> n + eta
1277 "9036", // id=66 N*(1535)0 --> Delta+ + pi-
1278 "7034", // id=66 N*(1535)0 --> Delta0 + pi0
1279 "8037", // id=66 N*(1535)0 --> Delta- + pi+
1280 "43014", // id=66 N*(1535)0 --> p + rho-
1281 "41013", // id=66 N*(1535)0 --> n + rho0
1282 "9008013", // id=66 N*(1535)0 --> n + pi+ + pi-
1283 "7007013", // id=66 N*(1535)0 --> n + pi0 + pi0
1284 "9038", // id=66 N*(1535)0 --> N*(1440)+ + pi-
1285 "7064", // id=66 N*(1535)0 --> N*(1440)0 + pi0
1286 "1013", // id=66 N*(1535)0 --> n + photon
1287 // J/Psi, 8 channels
1288 "2003", // J/Psi --> e+ + e-
1289 "5006", // J/Psi --> mu+ + mu-
1290 "1051", // J/Psi --> dilepton + photon
1291 "8009008009007", // J/Psi --> 2pi+ + 2pi- + pi0
1292 "8009008009008009007", // J/Psi --> 3pi+ + 3pi- + pi0
1293 "8009007", // J/Psi --> pi+ + pi- + pi0
1294 "12011008009007", // J/Psi --> K+ + K- + pi+ + pi- + pi0
1295 "4004", // J/Psi --> 2 neutrinos (this is a placeholder)
1296 // Psi' = Psi(2S), 8 channels
1297 "2003", // id=68 Psi' --> e+e-
1298 "5006", // id=68 Psi' --> mu+mu-
1299 "67008009", // id=68 Psi' --> J/Psi pi+ pi-
1300 "67007007", // id=68 Psi' --> J/Psi pi0 pi0
1301 "67017", // id=68 Psi' --> J/Psi eta
1302 "8009008009007", // id=68 Psi' --> 2pi+ + 2pi- + pi0
1303 "8009008009008009007", // id=68 Psi' --> 3pi+ + 3pi- + pi0
1304 "4004", // id=68 Psi' --> junk (=2 neutrinos as a placeholder)
1305 // pn, 2 channels
1306 "51013014", // id=69 pn --> dilepton + p + n
1307 "1013014" // id=69 pn --> photon + p + n
1308};
1309
1310// Decay-mode text description
1311const char* PStdData::DESCRIPTION[mnmodes] = {"mu+ --> e+ + neutrino + neutrino",
1312 "mu- --> e- + neutrino + neutrino",
1313 "pi0 --> photon + photon",
1314 "pi0 --> dilepton + photon (Dalitz)",
1315 "pi+ --> mu+ + neutrino",
1316 "pi- --> mu- + neutrino",
1317 "K0 long --> pi0 + pi0 + pi0",
1318 "K0 long --> pi+ + pi- + pi0",
1319 "K0 long --> pi+ + mu- + neutrino",
1320 "K0 long --> pi- + mu+ + neutrino",
1321 "K0 long --> pi+ + e- + neutrino",
1322 "K0 long --> pi- + e+ + neutrino",
1323 "K+ --> mu+ + neutrino",
1324 "K+ --> pi+ + pi0",
1325 "K+ --> pi+ + pi+ + pi-",
1326 "K+ --> pi+ + pi0 + pi0",
1327 "K+ --> pi0 + mu+ + neutrino",
1328 "K+ --> pi0 + e+ + neutrino",
1329 "K- --> mu- + neutrino",
1330 "K- --> pi- + pi0",
1331 "K- --> pi- + pi- + pi+",
1332 "K- --> pi- + pi0 + pi0",
1333 "K- --> pi0 + mu- + neutrino",
1334 "K- --> pi0 + e- + neutrino",
1335 "n --> p + e- + neutrino",
1336 "K0 short --> pi+ + pi-",
1337 "K0 short --> pi0 + pi0",
1338 "K0 short --> pi+ + pi- + photon",
1339 "eta --> photon + photon",
1340 "eta --> pi0 + pi0 + pi0",
1341 "eta --> pi+ + pi- + pi0",
1342 "eta --> pi+ + pi- + photon",
1343 "eta --> dilepton + photon (Dalitz)",
1344 "eta --> dimuon + photon",
1345 "eta --> e+ + e-",
1346 "eta --> mu+ + mu-",
1347 "Lambda --> p + pi-",
1348 "Lambda --> n + pi0",
1349 "Lambda --> n + photon",
1350 "anti_n --> anti_p + e+ + neutrino",
1351 "Delta0 --> p + pi-",
1352 "Delta0 --> n + pi0",
1353 "Delta0 --> n + photon",
1354 "Delta0 --> dilepton + n (Dalitz)",
1355 "Delta++ --> p + pi+",
1356 "Delta+ --> p + pi0",
1357 "Delta+ --> n + pi+",
1358 "Delta+ --> p + photon",
1359 "Delta+ --> dilepton + p (Dalitz)",
1360 "Delta- --> n + pi-",
1361 "N*(1440)+ --> p + pi0",
1362 "N*(1440)+ --> n + pi+",
1363 "N*(1440)+ --> Delta++ + pi-",
1364 "N*(1440)+ --> Delta+ + pi0",
1365 "N*(1440)+ --> Delta0 + pi+",
1366 "N*(1440)+ --> p + rho0",
1367 "N*(1440)+ --> n + rho+",
1368 "N*(1440)+ --> p + pi0 + pi0",
1369 "N*(1440)+ --> p + pi+ + pi-",
1370 "N*(1440)+ --> p + photon",
1371 "N*(1520)+ --> p + pi0",
1372 "N*(1520)+ --> n + pi+",
1373 "N*(1520)+ --> Delta++ + pi-",
1374 "N*(1520)+ --> Delta+ + pi0",
1375 "N*(1520)+ --> Delta0 + pi+",
1376 "N*(1520)+ --> p + rho0",
1377 "N*(1520)+ --> n + rho+",
1378 "N*(1520)+ --> p + pi0 + pi0",
1379 "N*(1520)+ --> p + pi+ + pi-",
1380 "N*(1520)+ --> p + photon",
1381 "N*(1535)+ --> p + pi0",
1382 "N*(1535)+ --> n + pi+",
1383 "N*(1535)+ --> p + eta",
1384 "N*(1535)+ --> Delta++ + pi-",
1385 "N*(1535)+ --> Delta+ + pi0",
1386 "N*(1535)+ --> Delta0 + pi+",
1387 "N*(1535)+ --> p + rho0",
1388 "N*(1535)+ --> n + rho+",
1389 "N*(1535)+ --> p + pi0 + pi0",
1390 "N*(1535)+ --> p + pi+ + pi-",
1391 "N*(1535)+ --> N*(1440)+ + pi0",
1392 "N*(1535)+ --> N*(1440)0 + pi+",
1393 "N*(1535)+ --> p + photon",
1394 "rho0 --> pi+ + pi-",
1395 "rho0 --> e+ + e-",
1396 "rho0 --> mu+ + mu-",
1397 "rho+ --> pi+ + pi0",
1398 "rho- --> pi- + pi0",
1399 "dimuon --> mu+ + mu-",
1400 "dilepton --> e+ + e-",
1401 "omega --> pi+ + pi- + pi0",
1402 "omega --> pi0 + photon",
1403 "omega --> pi+ + pi-",
1404 "omega --> dilepton + pi0 (Dalitz)",
1405 "omega --> dimuon + pi0",
1406 "omega --> e+ + e-",
1407 "omega --> mu+ + mu-",
1408 "eta' --> eta + pi- + pi+ ",
1409 "eta' --> rho0 + photon",
1410 "eta' --> eta + pi0 + pi0",
1411 "eta' --> omega + photon",
1412 "eta' --> photon + photon",
1413 "eta' --> pi0 + pi0 + pi0",
1414 "eta' --> dimuon + photon",
1415 "sigma --> pi+ + pi-",
1416 "sigma --> e+ + e-",
1417 "sigma --> mu+ + mu-",
1418 "phi --> K+ + K-",
1419 "phi --> K0L + K0S",
1420 "phi --> pi+ + pi- + pi0",
1421 "phi --> eta + photon",
1422 "phi --> pi0 + photon",
1423 "phi --> e+ + e-",
1424 "phi --> mu+ + mu-",
1425 "phi --> eta + dilepton",
1426 "Delta(1600)0 --> p + pi-",
1427 "Delta(1600)0 --> n + pi0",
1428 "Delta(1600)0 --> Delta+ + pi-",
1429 "Delta(1600)0 --> Delta0 + pi0",
1430 "Delta(1600)0 --> Delta- + pi+",
1431 "Delta(1600)0 --> p + rho-",
1432 "Delta(1600)0 --> n + rho0",
1433 "Delta(1600)0 --> N(1440)+ + rho-",
1434 "Delta(1600)0 --> N(1440)0 + rho0",
1435 "Delta(1600)0 --> n + photon",
1436 "Delta(1600)++ --> p + pi+",
1437 "Delta(1600)++ --> Delta++ + pi0",
1438 "Delta(1600)++ --> Delta+ + pi+",
1439 "Delta(1600)++ --> p + rho+",
1440 "Delta(1600)++ --> N(1440)+ + rho+",
1441 "Delta(1600)+ --> p + pi0",
1442 "Delta(1600)+ --> n + pi+",
1443 "Delta(1600)+ --> Delta++ + pi-",
1444 "Delta(1600)+ --> Delta+ + pi0",
1445 "Delta(1600)+ --> Delta0 + pi+",
1446 "Delta(1600)+ --> p + rho0",
1447 "Delta(1600)+ --> n + rho+",
1448 "Delta(1600)+ --> N(1440)+ + rho0",
1449 "Delta(1600)+ --> N(1440)0 + rho+",
1450 "Delta(1600)+ --> p + photon",
1451 "Delta(1600)- --> n + pi-",
1452 "Delta(1600)- --> Delta0 + pi-",
1453 "Delta(1600)- --> Delta- + pi0",
1454 "Delta(1600)- --> n + rho-",
1455 "Delta(1600)- --> N(1440)0 + rho-",
1456 "Delta(1620)0 --> p + pi-",
1457 "Delta(1620)0 --> n + pi0",
1458 "Delta(1620)0 --> Delta+ + pi-",
1459 "Delta(1620)0 --> Delta0 + pi0",
1460 "Delta(1620)0 --> Delta- + pi+",
1461 "Delta(1620)0 --> p + rho-",
1462 "Delta(1620)0 --> n + rho0",
1463 "Delta(1620)0 --> n + photon",
1464 "Delta(1620)++ --> p + pi+",
1465 "Delta(1620)++ --> Delta++ + pi0",
1466 "Delta(1620)++ --> Delta+ + pi+",
1467 "Delta(1620)++ --> p + rho+",
1468 "Delta(1620)+ --> p + pi0",
1469 "Delta(1620)+ --> n + pi+",
1470 "Delta(1620)+ --> Delta++ + pi-",
1471 "Delta(1620)+ --> Delta+ + pi0",
1472 "Delta(1620)+ --> Delta0 + pi+",
1473 "Delta(1620)+ --> p + rho0",
1474 "Delta(1620)+ --> n + rho+",
1475 "Delta(1620)+ --> p + photon",
1476 "Delta(1620)- --> n + pi-",
1477 "Delta(1620)- --> Delta0 + pi-",
1478 "Delta(1620)- --> Delta- + pi0",
1479 "Delta(1620)- --> n + rho-",
1480 "N*(1440)0 --> p + pi-",
1481 "N*(1440)0 --> n + pi0",
1482 "N*(1440)0 --> Delta+ + pi-",
1483 "N*(1440)0 --> Delta0 + pi0",
1484 "N*(1440)0 --> Delta- + pi+",
1485 "N*(1440)0 --> p + rho-",
1486 "N*(1440)0 --> n + rho0",
1487 "N*(1440)0 --> n + pi+ + pi-",
1488 "N*(1440)0 --> n + pi0 + pi0",
1489 "N*(1440)0 --> n + photon",
1490 "N*(1520)0 --> p + pi-",
1491 "N*(1520)0 --> n + pi0",
1492 "N*(1520)0 --> Delta+ + pi-",
1493 "N*(1520)0 --> Delta0 + pi0",
1494 "N*(1520)0 --> Delta- + pi+",
1495 "N*(1520)0 --> p + rho-",
1496 "N*(1520)0 --> n + rho0",
1497 "N*(1520)0 --> n + pi+ + pi-",
1498 "N*(1520)0 --> n + pi0 + pi0",
1499 "N*(1520)0 --> n + photon",
1500 "N*(1535)0 --> p + pi-",
1501 "N*(1535)0 --> n + pi0",
1502 "N*(1535)0 --> n + eta",
1503 "N*(1535)0 --> Delta+ + pi-",
1504 "N*(1535)0 --> Delta0 + pi0",
1505 "N*(1535)0 --> Delta- + pi+",
1506 "N*(1535)0 --> p + rho-",
1507 "N*(1535)0 --> n + rho0",
1508 "N*(1535)0 --> n + pi+ + pi-",
1509 "N*(1535)0 --> n + pi0 + pi0",
1510 "N*(1535)0 --> N*(1440)+ + pi-",
1511 "N*(1535)0 --> N*(1440)0 + pi0",
1512 "N*(1535)0 --> n + photon",
1513 "J/Psi --> e+ + e-",
1514 "J/Psi --> mu+ + mu-",
1515 "J/Psi --> dilepton + photon",
1516 "J/Psi --> 2pi+ + 2pi- + pi0",
1517 "J/Psi --> 3pi+ + 3pi- + pi0",
1518 "J/Psi --> pi+ + pi- + pi0",
1519 "J/Psi --> K+ + K- + pi+ + pi- + pi0",
1520 "J/Psi --> junk",
1521 "Psi' --> e+ + e-",
1522 "Psi' --> mu+ + mu-",
1523 "Psi' --> J/Psi + pi+ + pi-",
1524 "Psi' --> J/Psi + pi0 + pi0",
1525 "Psi' --> J/Psi + eta",
1526 "Psi' --> 2pi+ + 2pi- + pi0",
1527 "Psi' --> 3pi+ + 3pi- + pi0",
1528 "Psi' --> junk",
1529 "pn --> dilepton + p + n",
1530 "pn --> photon + p + n"};
1531
1532
1533// static members
1534int PStdData::maxnumpar = mnpar; // number of particles
1535int PStdData::maxnummodes = mnmodes; // number of decay modes
1536
1537char** PStdData::PName = (char**) NAME; // particle names
1538double* PStdData::PMass = (double*) MASS; // masses
1539double* PStdData::PWidth = (double*) WIDTH; // widths
1540int* PStdData::PMeson = (int*) MESON; // meson (1) or not (0)
1541int* PStdData::PBaryon = (int*) BARYON; // baryon number
1542int* PStdData::PLepton = (int*) LEPTON; // lepton number
1543int* PStdData::PCharge = (int*) CHARGE; // charge
1544int* PStdData::PJ = (int*) SPIN; // 2 x J
1545int* PStdData::PParity = (int*) PARITY; // Parity
1546int* PStdData::PI = (int*) ISPIN; // 2 x I
1547int* PStdData::PNModes = (int*) NMODES; // number of decay modes
1548int* PStdData::Pkf = (int*) PYTHIAKF; // Pythia6 kf code
1549double* PStdData::PBR = (double*) BRR; // branching ratio
1550char** PStdData::PMode = (char**) MODE; // decay mode coded by product ids
1551char** PStdData::PMDescription = (char**) DESCRIPTION; // text description of decay mode
1552
1553
1554//const long double PStdData::hbar=6.582122e-25; // units of (GeV s)
1555const double PStdData::hbar = 6.582122e-25; // units of (GeV s)
1556
1557
1558ClassImp(PStdData)
ClassImp
ClassImp(CbmConverterManager)
makeDataBase
PDataBase * makeDataBase()
Definition PDataBase.cxx:22
mnmodes
#define mnmodes
Definition PStdData.cxx:257
makeStdData
PStdData * makeStdData()
Definition PStdData.cxx:46
mnpar
#define mnpar
Definition PStdData.cxx:256
fStdData
PStdData & fStdData()
Definition PStdData.cxx:40
Int_t
int Int_t
Definition RootTypesDef.h:16
Bool_t
bool Bool_t
Definition RootTypesDef.h:15
Data
Class to handle QA-objects in the online reconstruction.
Definition QaData.h:27
PDataBase::MakeParamDouble
Int_t MakeParamDouble(const char *paramname, const char *descr)
Definition PDataBase.cxx:105
PDataBase::AddEntry
Bool_t AddEntry(Int_t key, const char *name)
Definition PDataBase.cxx:508
PDataBase::MakeParamTObj
Int_t MakeParamTObj(const char *paramname, const char *descr)
Definition PDataBase.cxx:155
PDataBase::MakeParamInt
Int_t MakeParamInt(const char *paramname, const char *descr)
Definition PDataBase.cxx:140
PDataBase::SetParamDouble
Bool_t SetParamDouble(Int_t key, const char *paramname, Double_t *result)
Definition PDataBase.cxx:370
PDataBase::SetFastKey
void SetFastKey(Int_t pkey, Int_t maxkey)
Definition PDataBase.cxx:75
PDataBase::AddListEntry
Int_t AddListEntry(const char *name, const char *count, const char *link, const char *newname)
Definition PDataBase.cxx:533
PDataBase::GetEntryInt
Int_t GetEntryInt(const char *paramname, Int_t value)
Definition PDataBase.cxx:464
PDataBase::SetParamInt
Bool_t SetParamInt(Int_t key, const char *paramname, Int_t *result)
Definition PDataBase.cxx:406
PStdData::MESON
static const int MESON[]
Definition PStdData.h:461
PStdData::DESCRIPTION
static const char * DESCRIPTION[]
Definition PStdData.h:1311
PStdData::PJ
static int * PJ
Definition PStdData.h:62
PStdData::PMDescription
static char ** PMDescription
Definition PStdData.h:65
PStdData::hbar
static const double hbar
Definition PStdData.h:91
PStdData::WIDTH
static double WIDTH[]
Definition PStdData.h:350
PStdData::maxnumpar
static int maxnumpar
Definition PStdData.h:62
PStdData::PMeson
static int * PMeson
Definition PStdData.h:62
PStdData::Pkf
static int * Pkf
Definition PStdData.h:62
PStdData::PLepton
static int * PLepton
Definition PStdData.h:62
PStdData::disable
int disable
Definition PStdData.h:42
PStdData::resetPosition
void resetPosition()
Definition PStdData.h:44
PStdData::NMODES
static const int NMODES[]
Definition PStdData.h:727
PStdData::CHARGE
static const int CHARGE[]
Definition PStdData.h:575
PStdData::fillDataBase
Bool_t fillDataBase(void)
Definition PStdData.cxx:58
PStdData::PWidth
static double * PWidth
Definition PStdData.h:64
PStdData::PParity
static int * PParity
Definition PStdData.h:62
PStdData::PName
static char ** PName
Definition PStdData.h:65
PStdData::PCharge
static int * PCharge
Definition PStdData.h:62
PStdData::SPIN
static const int SPIN[]
Definition PStdData.h:613
PStdData::BARYON
static const int BARYON[]
Definition PStdData.h:499
PStdData::PNModes
static int * PNModes
Definition PStdData.h:62
PStdData::PARITY
static const int PARITY[]
Definition PStdData.h:651
PStdData::BRR
static double BRR[]
Definition PStdData.h:776
PStdData::NAME
static const char * NAME[]
Definition PStdData.h:261
PStdData::PYTHIAKF
static const int PYTHIAKF[]
Definition PStdData.h:423
PStdData::PBaryon
static int * PBaryon
Definition PStdData.h:62
PStdData::PBR
static double * PBR
Definition PStdData.h:64
PStdData::MODE
static const char * MODE[]
Definition PStdData.h:1044
PStdData::PMode
static char ** PMode
Definition PStdData.h:65
PStdData::LEPTON
static const int LEPTON[]
Definition PStdData.h:537
PStdData::PI
static int * PI
Definition PStdData.h:62
PStdData::PPosition
int * PPosition
Definition PStdData.h:59
PStdData::ISPIN
static const int ISPIN[]
Definition PStdData.h:689
PStdData::maxnummodes
static int maxnummodes
Definition PStdData.h:62
PStdData::PMass
static double * PMass
Definition PStdData.h:64
PStdData::MASS
static double MASS[]
Definition PStdData.h:274