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