CMS 3D CMS Logo

List of all members | Static Public Member Functions
CustomPDGParser Class Reference

#include <CustomPDGParser.h>

Static Public Member Functions

static double s_charge (int pdg)
 
static std::vector< int > s_containedQuarks (int pdg)
 
static int s_containedQuarksCode (int pdg)
 
static bool s_isChargino (int pdg)
 
static bool s_isDphoton (int pdg)
 
static bool s_isMesonino (int pdg)
 
static bool s_isRBaryon (int pdg)
 
static bool s_isRGlueball (int pdg)
 
static bool s_isRHadron (int pdg)
 
static bool s_isRMeson (int pdg)
 
static bool s_isSbaryon (int pdg)
 
static bool s_issbottomHadron (int pdg)
 
static bool s_isSIMP (int pdg)
 
static bool s_isSLepton (int pdg)
 
static bool s_isstopHadron (int pdg)
 
static double s_spin (int pdg)
 

Detailed Description

Definition at line 6 of file CustomPDGParser.h.

Member Function Documentation

◆ s_charge()

double CustomPDGParser::s_charge ( int  pdg)
static

Definition at line 70 of file CustomPDGParser.cc.

70  {
71  float charge = 0, sign = 1;
72  int pdgAbs = abs(pdg);
73  if (pdg < 0)
74  sign = -1;
75 
76  if (s_isSLepton(pdg)) //Sleptons
77  {
78  if (pdgAbs % 2 == 0)
79  return 0;
80  else
81  return -sign;
82  }
83 
84  if (s_isDphoton(pdg)) {
85  return charge;
86  }
87  if (s_isChargino(pdg)) {
88  return sign;
89  }
90  if (s_isSIMP(pdg)) {
91  return 0;
92  }
93  if (s_isRMeson(pdg)) {
94  std::vector<int> quarks = s_containedQuarks(pdg);
95  if ((quarks[1] % 2 == 0 && quarks[0] % 2 == 1) || (quarks[1] % 2 == 1 && quarks[0] % 2 == 0))
96  charge = 1;
97  charge *= sign;
98  return charge;
99  }
100 
101  if (s_isRBaryon(pdg)) {
102  int baryon = s_containedQuarksCode(pdg);
103  for (int q = 1; q < 1000; q *= 10) {
104  if (baryon / q % 2 == 0)
105  charge += 2;
106  else
107  charge -= 1;
108  }
109  charge /= 3;
110  charge *= sign;
111  return charge;
112  }
113 
114  if (s_isMesonino(pdg)) {
115  int quark = s_containedQuarks(pdg)[0];
116  int squark = abs(pdg / 100 % 10);
117  if (squark % 2 == 0 && quark % 2 == 1)
118  charge = 1;
119  if (squark % 2 == 1 && quark % 2 == 0)
120  charge = 1;
121  charge *= sign;
122  if (s_issbottomHadron(pdg))
123  charge *= -1;
124  return charge;
125  }
126 
127  if (s_isSbaryon(pdg)) {
128  int baryon = s_containedQuarksCode(pdg) + 100 * (abs(pdg / 1000 % 10)); //Adding the squark back on
129  for (int q = 1; q < 1000; q *= 10) {
130  if (baryon / q % 2 == 0)
131  charge += 2;
132  else
133  charge -= 1;
134  }
135  charge /= 3;
136  charge *= sign;
137  if (s_issbottomHadron(pdg))
138  charge *= -1;
139  return charge;
140  }
141 
142  return 0;
143 }

References funct::abs(), ALCARECOTkAlJpsiMuMu_cff::charge, data-class-funcs::q, s_containedQuarks(), s_containedQuarksCode(), s_isChargino(), s_isDphoton(), s_isMesonino(), s_isRBaryon(), s_isRMeson(), s_isSbaryon(), s_issbottomHadron(), s_isSIMP(), s_isSLepton(), and Validation_hcalonly_cfi::sign.

Referenced by CustomParticleFactory::addCustomParticle().

◆ s_containedQuarks()

std::vector< int > CustomPDGParser::s_containedQuarks ( int  pdg)
static

Definition at line 156 of file CustomPDGParser.cc.

156  {
157  std::vector<int> quarks;
158  for (int i = s_containedQuarksCode(pdg); i > 0; i /= 10) {
159  quarks.push_back(i % 10);
160  }
161  return quarks;
162 }

References mps_fire::i, and s_containedQuarksCode().

Referenced by G4ProcessHelper::GetInclusiveCrossSection(), HadronicProcessHelper::inclusiveCrossSection(), and s_charge().

◆ s_containedQuarksCode()

int CustomPDGParser::s_containedQuarksCode ( int  pdg)
static

Definition at line 164 of file CustomPDGParser.cc.

164  {
165  int pdgAbs = abs(pdg);
166  if (s_isRBaryon(pdg))
167  return pdgAbs / 10 % 1000;
168 
169  if (s_isRMeson(pdg))
170  return pdgAbs / 10 % 100;
171 
172  if (s_isMesonino(pdg))
173  return pdgAbs / 10 % 1000 % 10;
174 
175  if (s_isSbaryon(pdg))
176  return pdgAbs / 10 % 1000 % 100;
177 
178  return 0;
179 }

References funct::abs(), s_isMesonino(), s_isRBaryon(), s_isRMeson(), and s_isSbaryon().

Referenced by s_charge(), and s_containedQuarks().

◆ s_isChargino()

bool CustomPDGParser::s_isChargino ( int  pdg)
static

Definition at line 60 of file CustomPDGParser.cc.

60  {
61  int pdgAbs = abs(pdg);
62  return (pdgAbs == 1000024);
63 }

References funct::abs().

Referenced by s_charge().

◆ s_isDphoton()

bool CustomPDGParser::s_isDphoton ( int  pdg)
static

Definition at line 40 of file CustomPDGParser.cc.

40  {
41  int pdgAbs = abs(pdg);
42  return (pdgAbs == 1072000);
43 }

References funct::abs().

Referenced by CustomParticleFactory::addCustomParticle(), and s_charge().

◆ s_isMesonino()

bool CustomPDGParser::s_isMesonino ( int  pdg)
static

Definition at line 50 of file CustomPDGParser.cc.

50  {
51  int pdgAbs = abs(pdg);
52  return ((pdgAbs % 10000 / 100 == 6) || (pdgAbs % 10000 / 100 == 5));
53 }

References funct::abs().

Referenced by CustomParticleFactory::addCustomParticle(), G4ProcessHelper::GetFinalState(), G4ProcessHelper::GetInclusiveCrossSection(), FullModelHadronicProcess::PostStepDoIt(), s_charge(), and s_containedQuarksCode().

◆ s_isRBaryon()

bool CustomPDGParser::s_isRBaryon ( int  pdg)
static

Definition at line 30 of file CustomPDGParser.cc.

30  {
31  int pdgAbs = abs(pdg);
32  return (pdgAbs % 100000 / 10000 == 9);
33 }

References funct::abs().

Referenced by G4ProcessHelper::GetFinalState(), G4ProcessHelper::GetInclusiveCrossSection(), G4ProcessHelper::ReactionGivesBaryon(), s_charge(), and s_containedQuarksCode().

◆ s_isRGlueball()

bool CustomPDGParser::s_isRGlueball ( int  pdg)
static

Definition at line 35 of file CustomPDGParser.cc.

35  {
36  int pdgAbs = abs(pdg);
37  return (pdgAbs % 100000 / 10 == 99);
38 }

References funct::abs().

Referenced by G4ProcessHelper::GetInclusiveCrossSection(), HadronicProcessHelper::inclusiveCrossSection(), and s_isRHadron().

◆ s_isRHadron()

bool CustomPDGParser::s_isRHadron ( int  pdg)
static

Definition at line 10 of file CustomPDGParser.cc.

10  {
11  int pdgAbs = abs(pdg);
12  return ((pdgAbs % 100000 / 10000 == 9) || (pdgAbs % 10000 / 1000 == 9) || s_isRGlueball(pdg));
13 }

References funct::abs(), and s_isRGlueball().

Referenced by CustomParticleFactory::addCustomParticle(), CustomPhysicsListSS::ConstructProcess(), CustomPhysicsList::ConstructProcess(), and CustomParticleFactory::getMassTable().

◆ s_isRMeson()

bool CustomPDGParser::s_isRMeson ( int  pdg)
static

Definition at line 45 of file CustomPDGParser.cc.

45  {
46  int pdgAbs = abs(pdg);
47  return (pdgAbs % 10000 / 1000 == 9);
48 }

References funct::abs().

Referenced by G4ProcessHelper::GetFinalState(), G4ProcessHelper::GetInclusiveCrossSection(), FullModelHadronicProcess::PostStepDoIt(), s_charge(), and s_containedQuarksCode().

◆ s_isSbaryon()

bool CustomPDGParser::s_isSbaryon ( int  pdg)
static

Definition at line 55 of file CustomPDGParser.cc.

55  {
56  int pdgAbs = abs(pdg);
57  return ((pdgAbs % 10000 / 1000 == 6) || (pdgAbs % 10000 / 1000 == 5));
58 }

References funct::abs().

Referenced by CustomParticleFactory::addCustomParticle(), G4ProcessHelper::GetFinalState(), G4ProcessHelper::GetInclusiveCrossSection(), G4ProcessHelper::ReactionGivesBaryon(), s_charge(), and s_containedQuarksCode().

◆ s_issbottomHadron()

bool CustomPDGParser::s_issbottomHadron ( int  pdg)
static

Definition at line 20 of file CustomPDGParser.cc.

20  {
21  int pdgAbs = abs(pdg);
22  return ((pdgAbs % 10000 / 1000 == 5) || (pdgAbs % 10000 / 100 == 5));
23 }

References funct::abs().

Referenced by CustomParticleFactory::addCustomParticle(), and s_charge().

◆ s_isSIMP()

bool CustomPDGParser::s_isSIMP ( int  pdg)
static

Definition at line 65 of file CustomPDGParser.cc.

65  {
66  int pdgAbs = abs(pdg);
67  return (pdgAbs == 9000006);
68 }

References funct::abs().

Referenced by CustomParticleFactory::addCustomParticle(), CustomParticleFactory::getMassTable(), and s_charge().

◆ s_isSLepton()

bool CustomPDGParser::s_isSLepton ( int  pdg)
static

Definition at line 25 of file CustomPDGParser.cc.

25  {
26  int pdgAbs = abs(pdg);
27  return (pdgAbs / 100 % 10000 == 0 && pdgAbs / 10 % 10 == 1);
28 }

References funct::abs().

Referenced by CustomParticleFactory::addCustomParticle(), and s_charge().

◆ s_isstopHadron()

bool CustomPDGParser::s_isstopHadron ( int  pdg)
static

Definition at line 15 of file CustomPDGParser.cc.

15  {
16  int pdgAbs = abs(pdg);
17  return ((pdgAbs % 10000 / 1000 == 6) || (pdgAbs % 1000 / 100 == 6));
18 }

References funct::abs().

Referenced by CustomParticleFactory::addCustomParticle(), and CustomParticleFactory::getMassTable().

◆ s_spin()

double CustomPDGParser::s_spin ( int  pdg)
static

Definition at line 145 of file CustomPDGParser.cc.

145  {
146  // The PDG numbering is described in the Review of Particle Physics:
147  // "3. In composite quark systems (diquarks, mesons, and baryons) ...
148  // the rightmost digit nJ = 2J + 1 gives the system's spin."
149  // Since this does not apply to SUSY / exotic particles,
150  // if the spin is important for the simulation
151  // it should be hard-coded based on PDG ID in this function.
152  int pdgAbs = abs(pdg);
153  return pdgAbs % 10;
154 }

References funct::abs().

Referenced by CustomParticleFactory::addCustomParticle().

CustomPDGParser::s_containedQuarksCode
static int s_containedQuarksCode(int pdg)
Definition: CustomPDGParser.cc:164
mps_fire.i
i
Definition: mps_fire.py:355
CustomPDGParser::s_isRMeson
static bool s_isRMeson(int pdg)
Definition: CustomPDGParser.cc:45
CustomPDGParser::s_isDphoton
static bool s_isDphoton(int pdg)
Definition: CustomPDGParser.cc:40
CustomPDGParser::s_issbottomHadron
static bool s_issbottomHadron(int pdg)
Definition: CustomPDGParser.cc:20
data-class-funcs.q
q
Definition: data-class-funcs.py:169
Validation_hcalonly_cfi.sign
sign
Definition: Validation_hcalonly_cfi.py:32
CustomPDGParser::s_isRBaryon
static bool s_isRBaryon(int pdg)
Definition: CustomPDGParser.cc:30
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
CustomPDGParser::s_isSIMP
static bool s_isSIMP(int pdg)
Definition: CustomPDGParser.cc:65
CustomPDGParser::s_isChargino
static bool s_isChargino(int pdg)
Definition: CustomPDGParser.cc:60
CustomPDGParser::s_isRGlueball
static bool s_isRGlueball(int pdg)
Definition: CustomPDGParser.cc:35
CustomPDGParser::s_isMesonino
static bool s_isMesonino(int pdg)
Definition: CustomPDGParser.cc:50
pdg
Definition: pdg_functions.h:28
CustomPDGParser::s_isSbaryon
static bool s_isSbaryon(int pdg)
Definition: CustomPDGParser.cc:55
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
CustomPDGParser::s_isSLepton
static bool s_isSLepton(int pdg)
Definition: CustomPDGParser.cc:25
CustomPDGParser::s_containedQuarks
static std::vector< int > s_containedQuarks(int pdg)
Definition: CustomPDGParser.cc:156