CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
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_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_isSLepton (int pdg)
 
static bool s_isstopHadron (int pdg)
 
static double s_spin (int pdg)
 

Detailed Description

Definition at line 7 of file CustomPDGParser.h.

Member Function Documentation

double CustomPDGParser::s_charge ( int  pdg)
static

Definition at line 77 of file CustomPDGParser.cc.

References funct::abs(), DeDxDiscriminatorTools::charge(), lumiQueryAPI::q, s_containedQuarks(), s_containedQuarksCode(), s_isChargino(), s_isMesonino(), s_isRBaryon(), s_isRMeson(), s_isSbaryon(), s_issbottomHadron(), and s_isSLepton().

Referenced by CustomParticleFactory::addCustomParticle().

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

Definition at line 153 of file CustomPDGParser.cc.

References i, and s_containedQuarksCode().

Referenced by s_charge().

154 {
155  std::vector<int> quarks;
156  for(int i=s_containedQuarksCode(pdg); i > 0; i /= 10)
157  {
158  quarks.push_back(i % 10);
159  }
160  return quarks;
161 }
static int s_containedQuarksCode(int pdg)
int i
Definition: DBlmapReader.cc:9
int CustomPDGParser::s_containedQuarksCode ( int  pdg)
static

Definition at line 163 of file CustomPDGParser.cc.

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

Referenced by s_charge(), and s_containedQuarks().

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 
179 return 0;
180 }
static bool s_isSbaryon(int pdg)
static bool s_isRMeson(int pdg)
static bool s_isMesonino(int pdg)
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
static bool s_isRBaryon(int pdg)
bool CustomPDGParser::s_isChargino ( int  pdg)
static

Definition at line 70 of file CustomPDGParser.cc.

References funct::abs().

Referenced by s_charge().

71 {
72  int pdgAbs = abs(pdg);
73  return (pdgAbs == 1000024);
74 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
bool CustomPDGParser::s_isMesonino ( int  pdg)
static

Definition at line 55 of file CustomPDGParser.cc.

References funct::abs().

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

56 {
57  int pdgAbs=abs(pdg);
58  return ((pdgAbs % 10000 / 100 == 6 ) || (pdgAbs % 10000 / 100 == 5));
59 
60 
61 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
bool CustomPDGParser::s_isRBaryon ( int  pdg)
static

Definition at line 34 of file CustomPDGParser.cc.

References funct::abs().

Referenced by s_charge(), and s_containedQuarksCode().

35 {
36  int pdgAbs=abs(pdg);
37  return (pdgAbs % 100000 / 10000 == 9);
38 
39 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
bool CustomPDGParser::s_isRGlueball ( int  pdg)
static

Definition at line 41 of file CustomPDGParser.cc.

References funct::abs().

Referenced by s_isRHadron().

42 {
43  int pdgAbs=abs(pdg);
44  return (pdgAbs % 100000 / 10 == 99);
45 
46 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
bool CustomPDGParser::s_isRHadron ( int  pdg)
static

Definition at line 10 of file CustomPDGParser.cc.

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

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

11 {
12  int pdgAbs=abs(pdg);
13  return ( (pdgAbs % 100000 / 10000 == 9) || (pdgAbs % 10000 / 1000 == 9) || s_isRGlueball(pdg) );
14 }
static bool s_isRGlueball(int pdg)
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
bool CustomPDGParser::s_isRMeson ( int  pdg)
static

Definition at line 48 of file CustomPDGParser.cc.

References funct::abs().

Referenced by s_charge(), and s_containedQuarksCode().

49 {
50  int pdgAbs=abs(pdg);
51  return (pdgAbs % 10000 / 1000 == 9);
52 
53 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
bool CustomPDGParser::s_isSbaryon ( int  pdg)
static

Definition at line 63 of file CustomPDGParser.cc.

References funct::abs().

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

64 {
65  int pdgAbs=abs(pdg);
66  return ((pdgAbs % 10000 / 1000 == 6) || (pdgAbs % 10000 / 1000 == 5));
67 
68 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
bool CustomPDGParser::s_issbottomHadron ( int  pdg)
static

Definition at line 22 of file CustomPDGParser.cc.

References funct::abs().

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

23 {
24  int pdgAbs=abs(pdg);
25  return ( (pdgAbs % 10000 / 1000 == 5) || (pdgAbs % 10000 / 100 == 5) );
26 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
bool CustomPDGParser::s_isSLepton ( int  pdg)
static

Definition at line 28 of file CustomPDGParser.cc.

References funct::abs().

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

29 {
30  int pdgAbs=abs(pdg);
31  return (pdgAbs / 100 % 10000 == 0 && pdgAbs / 10 % 10 == 1);
32 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
bool CustomPDGParser::s_isstopHadron ( int  pdg)
static

Definition at line 16 of file CustomPDGParser.cc.

References funct::abs().

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

17 {
18  int pdgAbs=abs(pdg);
19  return ( (pdgAbs % 10000 / 1000 == 6) || (pdgAbs % 1000 / 100 == 6) );
20 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
double CustomPDGParser::s_spin ( int  pdg)
static

Definition at line 141 of file CustomPDGParser.cc.

References funct::abs().

Referenced by CustomParticleFactory::addCustomParticle().

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