Gaudi::ParticleID Class Reference

Holds PDG + LHCb extension particle code, following the PDG particle numbering scheme (pdg.lbl.gov/2017/reviews/rpp2017-rev-monte-carlo-numbering.pdf). More...

#include </builds/gaudi/Gaudi/GaudiPartProp/include/Gaudi/ParticleID.h>

Public Types
enum	Location {   nj = 1 , nq3 , nq2 , nq1 ,   nl , nr , n , n8 ,   n9 , n10 }
	PDG ID digits (base 10) are: n nr nl nq1 ne2 nq3 nj. More...
enum	Quark {   down = 1 , up , strange , charm ,   bottom , top , bottom_prime , top_prime ,   first = down , last = top_prime }
	Quark PDG IDs. More...

Public Member Functions
	ParticleID (const int pid=0)
	Constructor with PDG code.
int	pid () const
	Retrieve the PDG ID.
constexpr unsigned int	abspid () const
	Absolute value of the PDG ID.
void	setPid (const int pid)
	Update the PDG ID.
bool	isValid () const
	Return if the PID is valid.
bool	isSM () const
	Return if the PID is from the standard model.
bool	isMeson () const
	Return if the PID is for a meson.
bool	isBaryon () const
	Return if the PID is for a baryon.
bool	isDiQuark () const
	Return if the PID is for a di-quark.
bool	isHadron () const
	Return if the PID is for a hadron.
bool	isLepton () const
	Return if the PID is for a lepton.
bool	isNucleus () const
	Return if the PID is for a nucleus.
bool	isQuark () const
	Return if the PID is for a bare quark.
bool	hasQuarks () const
	Return if the PID is a particle with quarks, but not a nucleus.
bool	hasQuark (const Quark &q) const
	Return if the PID is a particle containing a specified quark flavor.
bool	hasDown () const
	Return if the PID is a particle with a down quark.
bool	hasUp () const
	Return if the PID is a particle with an up quark.
bool	hasStrange () const
	Return if the PID is a particle with a down quark.
bool	hasCharm () const
	Return if the PID is a particle with a charm quark.
bool	hasBottom () const
	Return if the PID is a particle with a bottom quark.
bool	hasTop () const
	Return if the PID is a particle with a top quark.
bool	hasBottomPrime () const
	Return if the PID is a particle with a bottom' quark.
bool	hasTopPrime () const
	Return if the PID is a particle with a top' quark.
int	threeCharge () const
	Return three times the charge, in units of e+, valid for all particles.
int	jSpin () const
	Return 2J+1, where J is the total spin, valid for all particles.
int	sSpin () const
	Return 2S+1, where S is the spin, valid only for mesons.
int	lSpin () const
	Return 2L+1, where L is the orbital angular momentum, valid only for mesons.
int	Z () const
	Return the atomic number for a nucleus.
int	A () const
	Return the nucleon number for a nucleus.
int	nLambda () const
	Return the number of strange quarks for a nucleus.
int	fundamentalID () const
	Return the fundamental ID.
int	extraBits () const
	Return everything beyond the 7th PDG ID digit.
constexpr unsigned short	digit (const Location &loc) const
	Return the digit for a given PDG ID digit location.
bool	operator== (const ParticleID &o) const
	Equality operator.
bool	operator!= (const ParticleID &o) const
	Non-equality operator.
bool	operator< (const ParticleID &o) const
	Comparison operator.
std::ostream &	fillStream (std::ostream &s) const
	Fill a stream with the PID.
std::string	toString () const
	Return the PID stream representation as a string.

Static Public Member Functions
static std::ostream &	printLocation (const long l, std::ostream &s)
	Fill a stream with the PID digit enumeration.
static std::string	printLocation (const long l)
	Return the PID digit enumeration stream representation as a string.
static std::ostream &	printQuark (const long q, std::ostream &s)
	Fill a stream with the PID quark enumeration.
static std::string	printQuark (const long q)
	Return the PID quark enumeration stream representation as a string.

Private Attributes
int	m_pid { 0 }
	PDG ID.

Detailed Description

Holds PDG + LHCb extension particle code, following the PDG particle numbering scheme (pdg.lbl.gov/2017/reviews/rpp2017-rev-monte-carlo-numbering.pdf).

Specific conventions followed by Pythia 8 for beyond the standard model physics and color-octet quarkonia have been introduced.

Nuclei with the PDG 2017 convention (following the 2006 Monte Carlo nuclear code scheme) have the numbering +/- 10LZZZAAAI. Where AAA is A - the total baryon number, ZZZ is Z - the total number of protons, L is the total number of strange quarks, and I is the isomer number where I = 0 corresponds to the ground state. Backwards compatibility with the old heavy ion scheme has also been kept.

Date

19/02/2002

Author

Gloria Corti

Date

22/03/2018

Author

Philip Ilten

Definition at line 43 of file ParticleID.h.

Member Enumeration Documentation

Location

enum Gaudi::ParticleID::Location

PDG ID digits (base 10) are: n nr nl nq1 ne2 nq3 nj.

Enumerator
nj
nq3
nq2
nq1
nl
nr
n
n8
n9
n10

Definition at line 46 of file ParticleID.h.

{ nj = 1, nq3, nq2, nq1, nl, nr, n, n8, n9, n10 };

Quark

enum Gaudi::ParticleID::Quark

Quark PDG IDs.

Enumerator
down
up
strange
charm
bottom
top
bottom_prime
top_prime
first
last

Definition at line 48 of file ParticleID.h.

{ down = 1, up, strange, charm, bottom, top, bottom_prime, top_prime, first = down, last = top_prime };

Constructor & Destructor Documentation

ParticleID()

Gaudi::ParticleID::ParticleID ( const int pid = 0 )

inlineexplicit

Constructor with PDG code.

Definition at line 51 of file ParticleID.h.

{ setPid( pid ); }

Member Function Documentation

A()

int Gaudi::ParticleID::A

(

)

const

Return the nucleon number for a nucleus.

Definition at line 416 of file ParticleID.cpp.

                             {
  if ( extraBits() < 100 ) { return 0; }
  //
  DATA               data  = decode( *this );
  const unsigned int m_nq1 = data[Gaudi::ParticleID::nq1];
  const unsigned int m_nq2 = data[Gaudi::ParticleID::nq2];
  const unsigned int m_nq3 = data[Gaudi::ParticleID::nq3];
  //
  return m_nq1 * 100 + m_nq2 * 10 + m_nq3;
}

abspid()

unsigned int Gaudi::ParticleID::abspid ( ) const

inlineconstexpr

Absolute value of the PDG ID.

Definition at line 56 of file ParticleID.h.

{ return 0 > m_pid ? -m_pid : m_pid; }

digit()

unsigned short Gaudi::ParticleID::digit ( const Location & loc ) const

inlineconstexpr

Return the digit for a given PDG ID digit location.

Definition at line 125 of file ParticleID.h.

                                                                {
      constexpr std::uint32_t pows[10] = { 1, 10, 100, 1000, 10000, 100000, 1000000, 10000000, 100000000, 1000000000 };
      return ( abspid() / pows[loc - 1] ) % 10;
    }

extraBits()

int Gaudi::ParticleID::extraBits

(

)

const

Return everything beyond the 7th PDG ID digit.

Definition at line 158 of file ParticleID.cpp.

                                     {
  const unsigned int m_extra = abspid() / 10000000;
  return 45 == m_extra ? 100u : 46 == m_extra ? 100u : 47 == m_extra ? 100u : 49 == m_extra ? 100u : m_extra;
}

fillStream()

std::ostream & Gaudi::ParticleID::fillStream

(

std::ostream &

s

)

const

Fill a stream with the PID.

Definition at line 432 of file ParticleID.cpp.

                                                             {
  return s << "Gaudi.ParticleID(" << m_pid << ")";
}

fundamentalID()

int Gaudi::ParticleID::fundamentalID

(

)

const

Return the fundamental ID.

This is 0 for nuclei, mesons, baryons, and di-quarks. Otherwise, this is the first two digits of the PDG ID

Definition at line 166 of file ParticleID.cpp.

                                         {
  //
  DATA               data    = decode( *this );
  const unsigned int m_extra = data[0];
  const unsigned int m_nl    = data[Gaudi::ParticleID::nl];
  const unsigned int m_nq1   = data[Gaudi::ParticleID::nq1];
  const unsigned int m_nq2   = data[Gaudi::ParticleID::nq2];
  const unsigned int m_nq3   = data[Gaudi::ParticleID::nq3];
  const unsigned int m_nj    = data[Gaudi::ParticleID::nj];
  //
  return !m_extra && !( m_nl + m_nq1 + m_nq2 ) ? 10 * m_nq3 + m_nj : 0;
}

hasBottom()

bool Gaudi::ParticleID::hasBottom ( ) const

inline

Return if the PID is a particle with a bottom quark.

Definition at line 92 of file ParticleID.h.

{ return hasQuark( bottom ); }

hasBottomPrime()

bool Gaudi::ParticleID::hasBottomPrime ( ) const

inline

Return if the PID is a particle with a bottom' quark.

Definition at line 96 of file ParticleID.h.

{ return hasQuark( bottom_prime ); }

hasCharm()

bool Gaudi::ParticleID::hasCharm ( ) const

inline

Return if the PID is a particle with a charm quark.

Definition at line 90 of file ParticleID.h.

{ return hasQuark( charm ); }

hasDown()

bool Gaudi::ParticleID::hasDown ( ) const

inline

Return if the PID is a particle with a down quark.

Definition at line 84 of file ParticleID.h.

{ return hasQuark( down ); }

hasQuark()

bool Gaudi::ParticleID::hasQuark

(

const Quark &

q

)

const

Return if the PID is a particle containing a specified quark flavor.

Definition at line 257 of file ParticleID.cpp.

                                                     {
  if ( abspid() == q ) { return true; }
  if ( !hasQuarks() ) { return false; }
  //
  DATA               data  = decode( *this );
  const unsigned int m_nq1 = data[Gaudi::ParticleID::nq1];
  const unsigned int m_nq2 = data[Gaudi::ParticleID::nq2];
  const unsigned int m_nq3 = data[Gaudi::ParticleID::nq3];
  //
  return q == m_nq3 || q == m_nq2 || q == m_nq1;
}

hasQuarks()

bool Gaudi::ParticleID::hasQuarks

(

)

const

Return if the PID is a particle with quarks, but not a nucleus.

Definition at line 255 of file ParticleID.cpp.

{ return !extraBits() && !fundamentalID(); }

hasStrange()

bool Gaudi::ParticleID::hasStrange ( ) const

inline

Return if the PID is a particle with a down quark.

Definition at line 88 of file ParticleID.h.

{ return hasQuark( strange ); }

hasTop()

bool Gaudi::ParticleID::hasTop ( ) const

inline

Return if the PID is a particle with a top quark.

Definition at line 94 of file ParticleID.h.

{ return hasQuark( top ); }

hasTopPrime()

bool Gaudi::ParticleID::hasTopPrime ( ) const

inline

Return if the PID is a particle with a top' quark.

Definition at line 98 of file ParticleID.h.

{ return hasQuark( top_prime ); }

hasUp()

bool Gaudi::ParticleID::hasUp ( ) const

inline

Return if the PID is a particle with an up quark.

Definition at line 86 of file ParticleID.h.

{ return hasQuark( up ); }

isBaryon()

bool Gaudi::ParticleID::isBaryon

(

)

const

Return if the PID is for a baryon.

Definition at line 215 of file ParticleID.cpp.

                                     {
  //
  if ( !isSM() ) { return false; }
  //
  DATA               data  = decode( *this );
  const unsigned int m_nq1 = data[Gaudi::ParticleID::nq1];
  const unsigned int m_nq2 = data[Gaudi::ParticleID::nq2];
  const unsigned int m_nq3 = data[Gaudi::ParticleID::nq3];
  const unsigned int m_nj  = data[Gaudi::ParticleID::nj];
  //
  return m_nq1 != 0 && m_nq2 != 0 && m_nq3 != 0 && m_nj != 0;
}

isDiQuark()

bool Gaudi::ParticleID::isDiQuark

(

)

const

Return if the PID is for a di-quark.

Definition at line 228 of file ParticleID.cpp.

                                      {
  //
  if ( !isSM() ) { return false; }
  //
  DATA               data  = decode( *this );
  const unsigned int m_nq1 = data[Gaudi::ParticleID::nq1];
  const unsigned int m_nq2 = data[Gaudi::ParticleID::nq2];
  const unsigned int m_nq3 = data[Gaudi::ParticleID::nq3];
  const unsigned int m_nj  = data[Gaudi::ParticleID::nj];
  //
  return m_nq1 != 0 && m_nq2 != 0 && m_nq3 == 0 && m_nj != 0 && !( m_nj == 1 && m_nq1 == m_nq2 );
}

isHadron()

bool Gaudi::ParticleID::isHadron

(

)

const

Return if the PID is for a hadron.

Definition at line 241 of file ParticleID.cpp.

{ return isMeson() || isBaryon(); }

isLepton()

bool Gaudi::ParticleID::isLepton

(

)

const

Return if the PID is for a lepton.

Definition at line 243 of file ParticleID.cpp.

                                     {
  const unsigned int apid = abspid();
  return 11 <= apid && apid <= 18;
}

isMeson()

bool Gaudi::ParticleID::isMeson

(

)

const

Return if the PID is for a meson.

Definition at line 202 of file ParticleID.cpp.

                                    {
  //
  if ( !isSM() ) { return false; }
  //
  DATA               data  = decode( *this );
  const unsigned int m_nq1 = data[Gaudi::ParticleID::nq1];
  const unsigned int m_nq2 = data[Gaudi::ParticleID::nq2];
  const unsigned int m_nq3 = data[Gaudi::ParticleID::nq3];
  const unsigned int m_nj  = data[Gaudi::ParticleID::nj];
  //
  return m_nq1 == 0 && m_nq2 != 0 && m_nq3 != 0 && m_nj != 0 && !( pid() < 0 && m_nq2 == m_nq3 );
}

isNucleus()

bool Gaudi::ParticleID::isNucleus

(

)

const

Return if the PID is for a nucleus.

Definition at line 248 of file ParticleID.cpp.

{ return 100 <= extraBits() && Z() <= A(); }

isQuark()

bool Gaudi::ParticleID::isQuark

(

)

const

Return if the PID is for a bare quark.

Definition at line 250 of file ParticleID.cpp.

                                    {
  const unsigned int apid = abspid();
  return Gaudi::ParticleID::Quark::first <= apid && apid <= Gaudi::ParticleID::Quark::last;
}

isSM()

bool Gaudi::ParticleID::isSM

(

)

const

Return if the PID is from the standard model.

Definition at line 194 of file ParticleID.cpp.

                                 {
  DATA               data = decode( *this );
  const unsigned int m_n  = data[Gaudi::ParticleID::n];
  const unsigned int m_nr = data[Gaudi::ParticleID::nr];
  //
  return m_n < 1 || ( m_n == 9 && m_nr != 9 );
}

isValid()

bool Gaudi::ParticleID::isValid

(

)

const

Return if the PID is valid.

Definition at line 179 of file ParticleID.cpp.

                                    {
  // Nuclei, geantino, and intermediate.
  if ( extraBits() ) {
    return m_pid == 480000000 || m_pid == -990000000 || isNucleus();
  }
  // Standard model particles.
  else if ( abspid() < 100 || isHadron() || isDiQuark() ) {
    return true;
  }
  // Beyond the standard model particles.
  else {
    return !isSM();
  }
}

jSpin()

int Gaudi::ParticleID::jSpin

(

)

const

Return 2J+1, where J is the total spin, valid for all particles.

Definition at line 313 of file ParticleID.cpp.

                                 {
  // Lookup table for spins. The Xu+ (43) and Xu0 (44) break the PDG
  // conventions and are EvtGen specific. The c-hadron (84), b-hadron
  // (85), t-hadron (86), b'-hadron (87), are LHCb specific but within
  // the acceptable 81 - 100 range. A double charged Higgs is inserted
  // as 9 and a gravintino is inserted as 10.
  static const unsigned int j[101] = { 0, 2, 2, 2, 2, 2, 2, 2, 2, 1, 4, // 1  - 10.
                                       2, 2, 2, 2, 2, 2, 2, 2, 0, 0,    // 11 - 20.
                                       3, 3, 3, 3, 1, 0, 0, 0, 0, 0,    // 21 - 30.
                                       0, 3, 3, 3, 1, 1, 1, 1, 5, 0,    // 31 - 40.
                                       3, 1, 0, 0, 1, 1, 0, 0, 0, 0,    // 41 - 50.
                                       1, 2, 3, 1, 3, 1, 1, 1, 1, 1,    // 51 - 60.
                                       0, 0, 0, 0, 0, 0, 0, 0, 0, 0,    // 61 - 70.
                                       0, 0, 0, 0, 0, 0, 0, 0, 0, 0,    // 71 - 80.
                                       0, 0, 0, 0, 0, 0, 0, 0, 0, 0,    // 81 - 90.
                                       0, 0, 3, 0, 0, 0, 0, 0, 0, 0 };  // 91 - 100.
  const unsigned int        fid    = fundamentalID();
  int                       spin   = 0;
  //
  DATA               data  = decode( *this );
  const unsigned int m_n   = data[Gaudi::ParticleID::n];
  const unsigned int m_nq3 = data[Gaudi::ParticleID::nq3];
  const unsigned int m_nj  = data[Gaudi::ParticleID::nj];
  //
  // Fundamental particles from the spin table.
  if ( fid ) {
    // Supersymmetric particles.
    if ( m_n == 1 || m_n == 2 || m_n == 3 ) {
      spin = j[fid] % 2 + 1;
    }
    // Standard particles.
    else {
      spin = j[fid];
    }
  }
  // Nuclei.
  else if ( isNucleus() ) {
    spin = m_nq3 % 2 + 1;
  }
  // Everything else.
  else {
    spin = m_nj;
  }
  //
  return spin;
}

lSpin()

int Gaudi::ParticleID::lSpin

(

)

const

Return 2L+1, where L is the orbital angular momentum, valid only for mesons.

Definition at line 382 of file ParticleID.cpp.

                                 {
  //
  if ( !isMeson() ) { return 0; }
  //
  DATA               data = decode( *this );
  const unsigned int m_n  = data[Gaudi::ParticleID::n];
  const unsigned int m_nl = data[Gaudi::ParticleID::nl];
  const unsigned int m_nj = data[Gaudi::ParticleID::nj];
  //
  if ( m_n == 9 ) {
    return 0;
  } else if ( m_nj == 1 ) {
    return m_nl;
  } else if ( m_nl == 0 ) {
    return ( m_nj - 3 ) / 2;
  } else if ( m_nl == 3 ) {
    return ( m_nj + 1 ) / 2;
  } else {
    return ( m_nj - 1 ) / 2;
  }
}

nLambda()

int Gaudi::ParticleID::nLambda

(

)

const

Return the number of strange quarks for a nucleus.

Definition at line 427 of file ParticleID.cpp.

                                   {
  const unsigned int m_extra = extraBits();
  return 100 <= m_extra ? m_extra - 100 : 0;
}

operator!=()

bool Gaudi::ParticleID::operator!= ( const ParticleID & o ) const

inline

Non-equality operator.

Definition at line 133 of file ParticleID.h.

{ return m_pid != o.m_pid; }

operator<()

bool Gaudi::ParticleID::operator< ( const ParticleID & o ) const

inline

Comparison operator.

Definition at line 135 of file ParticleID.h.

                                                {
      const unsigned int i1( abspid() ), i2( o.abspid() );
      return std::tie( i1, m_pid ) < std::tie( i2, o.m_pid );
    }

operator==()

bool Gaudi::ParticleID::operator== ( const ParticleID & o ) const

inline

Equality operator.

Definition at line 131 of file ParticleID.h.

{ return m_pid == o.m_pid; }

pid()

int Gaudi::ParticleID::pid ( ) const

inline

Retrieve the PDG ID.

Definition at line 54 of file ParticleID.h.

{ return m_pid; }

printLocation() [1/2]

std::string Gaudi::ParticleID::printLocation ( const long l )

static

Return the PID digit enumeration stream representation as a string.

Definition at line 471 of file ParticleID.cpp.

                                                       {
  std::ostringstream s;
  printLocation( l, s );
  return s.str();
}

printLocation() [2/2]

std::ostream & Gaudi::ParticleID::printLocation ( const long l, std::ostream & s )

static

Fill a stream with the PID digit enumeration.

Definition at line 443 of file ParticleID.cpp.

                                                                        {
  switch ( l ) {
  case Gaudi::ParticleID::nj:
    return s << "Gaudi.ParticleID.nj";
  case Gaudi::ParticleID::nq3:
    return s << "Gaudi.ParticleID.nq3";
  case Gaudi::ParticleID::nq2:
    return s << "Gaudi.ParticleID.nq2";
  case Gaudi::ParticleID::nq1:
    return s << "Gaudi.ParticleID.nq1";
  case Gaudi::ParticleID::nl:
    return s << "Gaudi.ParticleID.nl";
  case Gaudi::ParticleID::nr:
    return s << "Gaudi.ParticleID.nr";
  case Gaudi::ParticleID::n:
    return s << "Gaudi.ParticleID.n";
  case Gaudi::ParticleID::n8:
    return s << "Gaudi.ParticleID.n8";
  case Gaudi::ParticleID::n9:
    return s << "Gaudi.ParticleID.n9";
  case Gaudi::ParticleID::n10:
    return s << "Gaudi.ParticleID.n10";
  default:;
  }
  return s << "Gaudi.ParticleID.Location(" << l << ")";
}

printQuark() [1/2]

std::string Gaudi::ParticleID::printQuark ( const long q )

static

Return the PID quark enumeration stream representation as a string.

Definition at line 502 of file ParticleID.cpp.

                                                    {
  std::ostringstream s;
  printQuark( q, s );
  return s.str();
}

printQuark() [2/2]

std::ostream & Gaudi::ParticleID::printQuark ( const long q, std::ostream & s )

static

Fill a stream with the PID quark enumeration.

Definition at line 478 of file ParticleID.cpp.

                                                                     {
  switch ( q ) {
  case Gaudi::ParticleID::down:
    return s << "Gaudi.ParticleID.down";
  case Gaudi::ParticleID::up:
    return s << "Gaudi.ParticleID.up";
  case Gaudi::ParticleID::strange:
    return s << "Gaudi.ParticleID.strange";
  case Gaudi::ParticleID::charm:
    return s << "Gaudi.ParticleID.charm";
  case Gaudi::ParticleID::bottom:
    return s << "Gaudi.ParticleID.bottom";
  case Gaudi::ParticleID::top:
    return s << "Gaudi.ParticleID.top";
  case Gaudi::ParticleID::bottom_prime:
    return s << "Gaudi.ParticleID.bottom_prime";
  case Gaudi::ParticleID::top_prime:
    return s << "Gaudi.ParticleID.top_prime";
  default:;
  }
  return s << "Gaudi.ParticleID.Quark(" << q << ")";
}

setPid()

void Gaudi::ParticleID::setPid ( const int pid )

inline

Update the PDG ID.

Definition at line 58 of file ParticleID.h.

{ m_pid = pid; }

sSpin()

int Gaudi::ParticleID::sSpin

(

)

const

Return 2S+1, where S is the spin, valid only for mesons.

Definition at line 361 of file ParticleID.cpp.

                                 {
  //
  if ( !isMeson() ) { return 0; }
  //
  DATA               data = decode( *this );
  const unsigned int m_n  = data[Gaudi::ParticleID::n];
  const unsigned int m_nl = data[Gaudi::ParticleID::nl];
  const unsigned int m_nj = data[Gaudi::ParticleID::nj];
  //
  if ( m_n == 9 ) {
    return 0;
  } else if ( m_nj == 1 ) {
    return m_nl;
  } else if ( m_nl == 1 ) {
    return 0;
  } else {
    return 1;
  }
}

threeCharge()

int Gaudi::ParticleID::threeCharge

(

)

const

Return three times the charge, in units of e+, valid for all particles.

Definition at line 269 of file ParticleID.cpp.

                                       {
  // Lookup table for charges. The Xu+ (43) and Xu0 (44) break the PDG
  // conventions and are EvtGen specific. The c-hadron (84), b-hadron
  // (85), t-hadron (86), b'-hadron (87), are LHCb specific but within
  // the acceptable 81 - 100 range. A double charged Higgs is inserted
  // as 9 and a gravintino as 10.
  static const int   q[101] = { 0,  -1, 2,  -1, 2,  -1, 2,  -1, 2, 6, 0, // 1  - 10.
                                -3, 0,  -3, 0,  -3, 0,  -3, 0,  0, 0,    // 11 - 20.
                                0,  0,  0,  3,  0,  0,  0,  0,  0, 0,    // 21 - 30.
                                0,  0,  0,  3,  0,  0,  3,  0,  0, 0,    // 31 - 40.
                                0,  -1, 0,  3,  0,  0,  0,  0,  0, 0,    // 41 - 50.
                                0,  0,  0,  0,  0,  0,  0,  0,  0, 0,    // 51 - 60.
                                0,  0,  0,  0,  0,  0,  0,  0,  0, 0,    // 61 - 70.
                                0,  0,  0,  0,  0,  0,  0,  0,  0, 0,    // 71 - 80.
                                0,  0,  0,  2,  -1, 2,  -1, 0,  0, 0,    // 81 - 90.
                                0,  0,  0,  0,  0,  0,  0,  0,  0, 0 };  // 91 - 100.
  const unsigned int fid    = fundamentalID();
  int                charge = 0;
  //
  DATA               data  = decode( *this );
  const unsigned int m_nq1 = data[Gaudi::ParticleID::nq1];
  const unsigned int m_nq2 = data[Gaudi::ParticleID::nq2];
  const unsigned int m_nq3 = data[Gaudi::ParticleID::nq3];
 
  // Nuclei.
  if ( isNucleus() ) {
    charge = 3 * Z();
  }
  // Fundamental particles from the charge table.
  else if ( fid ) {
    charge = q[fid];
  }
  // Mesons.
  else if ( !m_nq1 ) {
    charge = abs( q[m_nq2] - q[m_nq3] );
  }
  // Everything else.
  else {
    charge = q[m_nq3] + q[m_nq2] + q[m_nq1];
  }
 
  return ( pid() < 0 ? -1 : 1 ) * charge;
}

toString()

std::string Gaudi::ParticleID::toString

(

)

const

Return the PID stream representation as a string.

Definition at line 436 of file ParticleID.cpp.

                                          {
  std::ostringstream s;
  fillStream( s );
  return s.str();
}

Z()

int Gaudi::ParticleID::Z

(

)

const

Return the atomic number for a nucleus.

Definition at line 404 of file ParticleID.cpp.

                             {
  if ( extraBits() < 100 ) { return 0; }
  //
  DATA data = decode( *this );
  //
  const unsigned int m_n  = data[Gaudi::ParticleID::n];
  const unsigned int m_nr = data[Gaudi::ParticleID::nr];
  const unsigned int m_nl = data[Gaudi::ParticleID::nl];
  //
  return m_n * 100 + m_nr * 10 + m_nl;
}

Member Data Documentation

m_pid

int Gaudi::ParticleID::m_pid { 0 }

private

PDG ID.

Definition at line 155 of file ParticleID.h.

{ 0 };

---

The documentation for this class was generated from the following files:

• GaudiPartProp/include/Gaudi/ParticleID.h
• GaudiPartProp/src/ParticleID.cpp