|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
|
// -*- C++ -*-
//
// Package: FastSimulation/Particle
// Class : pdg_functions
//
// Implementation:
// [Notes on implementation]
//
// Original Author: Christopher Jones
// Created: Mon, 04 Mar 2019 19:53:06 GMT
//
// system include files
// user include files
#include "FastSimulation/Particle/interface/pdg_functions.h"
double pdg::mass(int pdgID, HepPDT::ParticleDataTable const* theTable) {
auto info = theTable->particle(HepPDT::ParticleID(pdgID));
if (info) {
return info->mass().value();
}
return kInvalidMass;
}
double pdg::cTau(int pdgID, HepPDT::ParticleDataTable const* theTable) {
auto info = theTable->particle(HepPDT::ParticleID(pdgID));
double ct = kInvalidCtau;
if (info) {
// The lifetime is 0. in the Pythia Particle Data Table !
// ct=tab->theTable()->particle(ParticleID(myId))->lifetime().value();
// Get it from the width (apparently Gamma/c!)
double w = info->totalWidth().value();
if (w != 0. && pdgID != 1000022) {
ct = 6.582119e-25 / w / 10.; // ctau in cm
} else {
// Temporary fix of a bug in the particle data table
unsigned amyId = abs(pdgID);
if (amyId != 22 && // photon
amyId != 11 && // e+/-
amyId != 10 && // nu_e
amyId != 12 && // nu_mu
amyId != 14 && // nu_tau
amyId != 1000022 && // Neutralino
amyId != 1000039 && // Gravitino
amyId != 2112 && // neutron/anti-neutron
amyId != 2212 && // proton/anti-proton
amyId != 101 && // Deutreron etc..
amyId != 102 && // Deutreron etc..
amyId != 103 && // Deutreron etc..
amyId != 104) { // Deutreron etc..
ct = 0.;
/* */
}
}
}
/*
std::cout << setw(20) << setprecision(18)
<< "myId/ctau/width = " << myId << " "
<< ct << " " << w << endl;
*/
return ct;
}
|