AtlasPID.h File Reference

#include <vector>
#include <cmath>
#include <algorithm>
#include <array>
#include <cstdlib>

Include dependency graph for AtlasPID.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
class	DecodedPID
	Implementation of classification functions according to PDG2022. More...

Namespaces
	SUSY

Functions
template<class T >
bool	isQuark (const T &p)
	PDG rule 2: Quarks and leptons are numbered consecutively starting from 1 and 11 respectively; to dothis they are first ordered by family and within families by weak isospin. More...
template<>
bool	isQuark (const int &p)
template<>
bool	isQuark (const DecodedPID &p)
template<class T >
bool	isStrange (const T &p)
template<>
bool	isStrange (const int &p)
template<class T >
bool	isCharm (const T &p)
template<>
bool	isCharm (const int &p)
template<class T >
bool	isBottom (const T &p)
template<>
bool	isBottom (const int &p)
template<class T >
bool	isTop (const T &p)
template<>
bool	isTop (const int &p)
template<class T >
bool	isLepton (const T &p)
	APID: the fourth generation leptons are leptons. More...
template<>
bool	isLepton (const int &p)
template<>
bool	isLepton (const DecodedPID &p)
template<class T >
bool	isSMLepton (const T &p)
template<>
bool	isSMLepton (const int &p)
template<>
bool	isSMLepton (const DecodedPID &p)
template<class T >
bool	isChLepton (const T &p)
	APID: the fourth generation leptons are leptons. More...
template<>
bool	isChLepton (const int &p)
template<class T >
bool	isElectron (const T &p)
template<>
bool	isElectron (const int &p)
template<class T >
bool	isMuon (const T &p)
template<>
bool	isMuon (const int &p)
template<class T >
bool	isTau (const T &p)
template<>
bool	isTau (const int &p)
template<class T >
bool	isNeutrino (const T &p)
	APID: the fourth generation neutrinos are neutrinos. More...
template<>
bool	isNeutrino (const int &p)
template<class T >
bool	isSMNeutrino (const T &p)
template<>
bool	isSMNeutrino (const int &p)
template<class T >
bool	isGluon (const T &p)
template<>
bool	isGluon (const int &p)
template<class T >
bool	isPhoton (const T &p)
template<>
bool	isPhoton (const int &p)
template<class T >
bool	isZ (const T &p)
template<>
bool	isZ (const int &p)
template<class T >
bool	isW (const T &p)
template<>
bool	isW (const int &p)
template<class T >
bool	isDM (const T &p)
	PDG rule 11j: The nature of Dark Matter (DM) is not known, and therefore a definitive classificationis too early. More...
template<>
bool	isDM (const int &p)
template<class T >
bool	isTrajectory (const T &p)
	PDG rule 8: The pomeron and odderon trajectories and a generic reggeon trajectory of states in QCD areassigned codes 990, 9990, and 110 respectively. More...
template<>
bool	isTrajectory (const int &p)
template<class T >
bool	isHiggs (const T &p)
	APID: HIGGS boson is only one particle. More...
template<>
bool	isHiggs (const int &p)
template<class T >
bool	isResonance (const T &p)
template<class T >
bool	isGraviton (const T &p)
template<>
bool	isGraviton (const int &p)
template<class T >
bool	isLeptoQuark (const T &p)
template<>
bool	isLeptoQuark (const int &p)
template<class T >
bool	isSUSY (const T &p)
template<class T >
bool	isTechnicolor (const T &p)
template<class T >
bool	isExcited (const T &p)
template<class T >
bool	isKK (const T &p)
template<class T >
bool	isMonopole (const T &p)
template<class T >
bool	isHiddenValley (const T &p)
template<class T >
bool	isGenericMultichargedParticle (const T &p)
template<class T >
bool	isDiquark (const T &p)
template<class T >
bool	isHadron (const T &p)
template<class T >
bool	isMeson (const T &p)
template<class T >
bool	isBaryon (const T &p)
template<class T >
bool	isTetraquark (const T &p)
template<class T >
bool	isPentaquark (const T &p)
template<class T >
bool	isNucleus (const T &p)
template<class T >
bool	isBSM (const T &p)
template<class T >
bool	isValid (const T &p)
template<>
bool	isValid (const DecodedPID &p)
	Av: we implement here an ATLAS-sepcific convention: all particles which are 99xxxxx are fine. More...
template<class T >
bool	isTransportable (const T &p)
template<class T >
bool	isGenSpecific (const T &p)
template<class T >
bool	isGeantino (const T &p)
template<class T >
bool	isPythia8Specific (const T &p)
template<>
bool	isPythia8Specific (const DecodedPID &p)
template<>
bool	isPythia8Specific (const int &p)
template<>
bool	isGenSpecific (const int &p)
	Main Table for MC internal use 81–100,901–930,998-999,1901–1930,2901–2930, and 3901–3930. More...
template<>
bool	isGeantino (const int &p)
template<class T >
bool	isGlueball (const T &p)
	APID: Definition of Glueballs: SM glueballs 99X (X=1,3,5), 999Y (Y=3,5,7) More...
template<>
bool	isGlueball (const DecodedPID &p)
template<>
bool	isGlueball (const int &p)
template<>
bool	isSUSY (const DecodedPID &p)
	PDG rule 11d Fundamental supersymmetric particles are identified by adding a nonzero n to the particle number. More...
template<>
bool	isTechnicolor (const DecodedPID &p)
	PDG rule 11e Technicolor states have n = 3, with technifermions treated like ordinary fermions. More...
template<>
bool	isTechnicolor (const int &p)
template<>
bool	isExcited (const DecodedPID &p)
	PDG rule 11f Excited (composite) quarks and leptons are identified by setting n= 4 and nr= 0. More...
template<>
bool	isExcited (const int &p)
template<>
bool	isKK (const DecodedPID &p)
	PDG rule 11h A black hole in models with extra dimensions has code 5000040. More...
template<>
bool	isKK (const int &p)
template<>
bool	isMonopole (const DecodedPID &p)
	PDG rule 11i Magnetic monopoles and dyons are assumed to have one unit of Dirac monopole charge and a variable integer number nq1nq2 nq3 units of electric charge. More...
template<>
bool	isMonopole (const int &p)
template<>
bool	isHiddenValley (const DecodedPID &p)
	PDG rule 11k Hidden Valley particles have n = 4 and n_r = 9, and trailing numbers in agreement with their nearest-analog standard particles, as far as possible. More...
template<>
bool	isHiddenValley (const int &p)
template<>
bool	isGenericMultichargedParticle (const DecodedPID &p)
	In addition, there is a need to identify ”Q-ball” and similar very exotic (multi-charged) particles which may have large, non-integer charge. More...
template<>
bool	isGenericMultichargedParticle (const int &p)
template<>
bool	isDiquark (const DecodedPID &p)
	PDG rule 4 Diquarks have 4-digit numbers with nq1 >= nq2 and nq3 = 0 APID: the diquarks with fourth generation are not diquarks. More...
template<>
bool	isMeson (const DecodedPID &p)
	Table 43.1 PDG rule 5a: The numbers specifying the meson’s quark content conform to the convention nq1= 0 and nq2 >= nq3. More...
template<>
bool	isBaryon (const DecodedPID &p)
	Table 43.2. More...
template<>
bool	isPentaquark (const DecodedPID &p)
	PDG rule 15 The 9-digit penta-quark codes are±1nrnLnq1nq2nq3nq4nq5nJ, sorted such thatnq1≥nq2≥nq3≥nq4. More...
template<>
bool	isTetraquark (const DecodedPID &p)
	PDG rule 14 The 9-digit tetra-quark codes are±1nrnLnq1nq20nq3nq4nJ. More...
template<>
bool	isNucleus (const DecodedPID &p)
	PDG rule 16: Nuclear codes are given as 10-digit numbers±10LZZZAAAI. More...
template<>
bool	isBSM (const DecodedPID &p)
	APID: graviton and all Higgs extensions are BSM. More...
template<>
bool	isSUSY (const int &p)
template<>
bool	isDiquark (const int &p)
template<>
bool	isMeson (const int &p)
template<>
bool	isBaryon (const int &p)
template<>
bool	isTetraquark (const int &p)
template<>
bool	isPentaquark (const int &p)
template<>
bool	isNucleus (const int &p)
template<>
bool	isBSM (const int &p)
template<>
bool	isHadron (const DecodedPID &p)
template<>
bool	isHadron (const int &p)
template<>
bool	isTransportable (const DecodedPID &p)
template<>
bool	isTransportable (const int &p)
template<>
bool	isValid (const int &p)
template<class T >
bool	hasQuark (const T &p, const int &q)
template<>
bool	hasQuark (const DecodedPID &p, const int &q)
template<>
bool	hasQuark (const int &p, const int &q)
template<class T >
int	leadingQuark (const T &p)
template<>
int	leadingQuark (const DecodedPID &p)
template<>
int	leadingQuark (const int &p)
template<class T >
bool	hasStrange (const T &p)
template<class T >
bool	hasCharm (const T &p)
template<class T >
bool	hasBottom (const T &p)
template<class T >
bool	hasTop (const T &p)
template<class T >
bool	isLightHadron (const T &p)
template<class T >
bool	isHeavyHadron (const T &p)
template<class T >
bool	isStrangeHadron (const T &p)
template<class T >
bool	isCharmHadron (const T &p)
template<class T >
bool	isBottomHadron (const T &p)
template<class T >
bool	isTopHadron (const T &p)
template<class T >
bool	isLightMeson (const T &p)
template<class T >
bool	isHeavyMeson (const T &p)
template<class T >
bool	isStrangeMeson (const T &p)
template<class T >
bool	isCharmMeson (const T &p)
template<class T >
bool	isBottomMeson (const T &p)
template<class T >
bool	isTopMeson (const T &p)
template<class T >
bool	isCCbarMeson (const T &p)
template<>
bool	isCCbarMeson (const DecodedPID &p)
template<>
bool	isCCbarMeson (const int &p)
template<class T >
bool	isBBbarMeson (const T &p)
template<>
bool	isBBbarMeson (const DecodedPID &p)
template<>
bool	isBBbarMeson (const int &p)
template<class T >
bool	isLightBaryon (const T &p)
template<class T >
bool	isHeavyBaryon (const T &p)
template<class T >
bool	isStrangeBaryon (const T &p)
template<class T >
bool	isCharmBaryon (const T &p)
template<class T >
bool	isBottomBaryon (const T &p)
template<class T >
bool	isTopBaryon (const T &p)
template<class T >
int	charge3 (const T &p)
template<class T >
double	fractionalCharge (const T &p)
template<class T >
double	charge (const T &p)
template<class T >
double	threeCharge (const T &p)
template<class T >
bool	isCharged (const T &p)
template<class T >
bool	isNeutral (const T &p)
template<>
int	charge3 (const DecodedPID &p)
template<>
int	charge3 (const int &p)
template<>
double	fractionalCharge (const DecodedPID &p)
template<>
double	fractionalCharge (const int &p)
template<class T >
bool	isEMInteracting (const T &p)
template<>
bool	isEMInteracting (const int &p)
template<class T >
bool	isParton (const T &p)
template<class T >
bool	SUSY::isRGlueball (const T &p)
	PDG rule 11g: Within several scenarios of new physics, it is possible to have colored particles sufficiently long-lived for color-singlet hadronic states to form around them. More...
template<>
bool	SUSY::isRGlueball (const DecodedPID &p)
template<>
bool	SUSY::isRGlueball (const int &p)
template<class T >
bool	SUSY::isRHadron (const T &p)
template<>
bool	SUSY::isRHadron (const DecodedPID &p)
template<>
bool	SUSY::isRHadron (const int &p)
template<class T >
bool	SUSY::isRMeson (const T &p)
template<>
bool	SUSY::isRMeson (const DecodedPID &p)
template<>
bool	SUSY::isRMeson (const int &p)
template<class T >
bool	SUSY::isRBaryon (const T &p)
template<>
bool	SUSY::isRBaryon (const DecodedPID &p)
template<>
bool	SUSY::isRBaryon (const int &p)
template<class T >
bool	SUSY::isSLepton (const T &p)
template<>
bool	SUSY::isSLepton (const DecodedPID &p)
template<>
bool	SUSY::isSLepton (const int &p)
template<class T >
bool	SUSY::isSBaryon (const T &p)
template<>
bool	SUSY::isSBaryon (const DecodedPID &p)
template<>
bool	SUSY::isSBaryon (const int &p)
template<class T >
bool	SUSY::isSMeson (const T &p)
template<>
bool	SUSY::isSMeson (const DecodedPID &p)
template<>
bool	SUSY::isSMeson (const int &p)
template<class T >
bool	SUSY::spin (const T &p)
template<>
bool	SUSY::spin (const int &p)
std::vector< int >	SUSY::containedQuarks (int p)
template<class T >
bool	SUSY::isSquark (const T &p)
template<>
bool	SUSY::isSquark (const DecodedPID &p)
template<>
bool	SUSY::isSquark (const int &p)
template<class T >
bool	SUSY::hasSquark (const T &p, const int &q)
template<>
bool	SUSY::hasSquark (const DecodedPID &p, const int &q)
template<>
bool	SUSY::hasSquark (const int &p, const int &q)
template<class T >
bool	isStrongInteracting (const T &p)
template<>
bool	isStrongInteracting (const int &p)

Function Documentation

charge()

template<class T >

double charge ( const T &  p )

inline

Definition at line 538 of file AtlasPID.h.

                                                   {
  if (isGenericMultichargedParticle(p)) // BSM multi-charged particles might have a fractional charge that's not a multiple of 1/3
    return fractionalCharge(p);
  else 
    return 1.0*charge3(p)/3.0; 
}

charge3() [1/3]

template<>

int charge3 ( const DecodedPID &  p )

inline

Codes 411nq1nq2 nq3 0 are then used when the magnetic and electrical charge sign agree and 412nq1nq2 nq3 0 when they disagree, with the overall sign of the particle set by the magnetic charge.

Definition at line 550 of file AtlasPID.h.

                                                   {
  auto ap = std::abs(p.pid());
  if (ap < TABLESIZE ) return p.pid() > 0 ? triple_charge.at(ap) : -triple_charge.at(ap);
  if (ap == K0) return 0;
  if (ap == GEANTINO0) return 0;
  if (ap == GEANTINOPLUS) return p.pid() > 0 ? 3 : -3;
  if (ap == MAVTOP) return p.pid() > 0 ? 2 : -2;
  size_t nq = 0;
  int sign = 1;
  int signmult = 1;
  int result=0;
  bool classified = false;
  if (!classified && isMeson(p)) { classified = true; nq = 2; if ((*(p.second.rbegin()+2)) == 2||(*(p.second.rbegin()+2)) == 4 ) { sign=-1;} signmult =-1; }
  if (!classified && isDiquark(p)) {return triple_charge.at(p(0))+triple_charge.at(p(1)); }
  if (!classified && isBaryon(p)) { classified = true; nq = 3; }
  if (!classified && isTetraquark(p)){ return triple_charge.at(p(3)) + triple_charge.at(p(4)) - triple_charge.at(p(6)) - triple_charge.at(p(7)); }
  if (!classified && isPentaquark(p)){ return triple_charge.at(p(3)) + triple_charge.at(p(4)) + triple_charge.at(p(5)) + triple_charge.at(p(6)) - triple_charge.at(p(7)); }
  if (!classified && isNucleus(p)) { classified = true; nq=0; result = 3*(p(3)*100 + p(4)*10 + p(5)) + (-1)*p(2);}
  if (!classified && isSUSY(p)) { nq = 0;
    classified = true;
    auto pp = p.shift(1); if (pp.ndigits() > 2) pp = pp.shift(1);
    return charge3(pp);
  }
  if (!classified && isMonopole(p)) {
    result = 3*(p(3)*100 + p(4)*10 + p(5));
    return ( (p.pid() > 0 && p(2) == 1) ||  (p.pid() < 0 && p(2) == 2) ) ? result : -result;
  }
  if (!classified && isGenericMultichargedParticle(p)) {
    double abs_charge = 0.0;
    if (p(0) == 1) abs_charge = p(3)*100. + p(4)*10. + p(5)*1 + p(6)*0.1; // multi-charged particle PDG ID is +/-100XXXY0, where the charge is XXX.Y
    if (p(0) == 2) abs_charge = (p(3)*10. + p(4))/(p(5)*10.0 + p(6)); // multi-charged particle PDG ID is +/-200XXYY0, where the charge is XX/YY
    int abs_threecharge = static_cast<int>(std::round(abs_charge * 3.)); // the multi-charged particles might have a fractional charge that's not a multiple of 1/3, in that case round to the closest multiple of 1/3 for charge3 and threecharge
    return p.pid() > 0 ? abs_threecharge : -1 * abs_threecharge;
  }
  for (auto r = p.second.rbegin() + 1; r != p.second.rbegin() + 1 + nq; ++r) {
    result += triple_charge.at(*r)*sign;
    sign*=signmult;
  }
  return p.pid() > 0 ? result : -result;
}

charge3() [2/3]

template<>

int charge3 ( const int &  p )

inline

Definition at line 592 of file AtlasPID.h.

                                           {
  int ap = std::abs(p);
  if (ap < TABLESIZE) return p > 0 ? triple_charge.at(ap):-triple_charge.at(ap);
  auto value_digits = DecodedPID(p);
  return charge3(value_digits);
}

charge3() [3/3]

template<class T >

int charge3 ( const T &  p )

inline

Definition at line 536 of file AtlasPID.h.

{return charge3(p->pdg_id());}

fractionalCharge() [1/3]

template<>

double fractionalCharge ( const DecodedPID &  p )

inline

Definition at line 599 of file AtlasPID.h.

                                                               {
  if(!isGenericMultichargedParticle(p)) return 1.0*charge3(p)/3.0; // this method is written for multi-charged particles, still make sure other cases are handled properly
  double abs_charge = 0;
  if (p(0) == 1) abs_charge = p(3)*100. + p(4)*10. + p(5)*1 + p(6)*0.1; // multi-charged particle PDG ID is +/-100XXXY0, where the charge is XXX.Y
  if (p(0) == 2) abs_charge = (p(3)*10. + p(4))/(p(5)*10.0 + p(6)); // multi-charged particle PDG ID is +/-200XXYY0, where the charge is XX/YY 
  return p.pid() > 0 ? abs_charge : -1 * abs_charge;
}

fractionalCharge() [2/3]

template<>

double fractionalCharge ( const int &  p )

inline

Definition at line 606 of file AtlasPID.h.

{auto value_digits = DecodedPID(p); return fractionalCharge(value_digits);}

fractionalCharge() [3/3]

template<class T >

double fractionalCharge ( const T &  p )

inline

Definition at line 537 of file AtlasPID.h.

{return fractionalCharge(p->pdg_id());}

hasBottom()

template<class T >

bool hasBottom ( const T &  p )

inline

Definition at line 502 of file AtlasPID.h.

{ return  hasQuark(p,BQUARK); }

hasCharm()

template<class T >

bool hasCharm ( const T &  p )

inline

Definition at line 501 of file AtlasPID.h.

{ return  hasQuark(p,CQUARK); }

hasQuark() [1/3]

template<>

bool hasQuark ( const DecodedPID &  p, const int &  q  )

inline

Definition at line 474 of file AtlasPID.h.

                                                                  {
  if (isQuark(p.pid())) { return (std::abs(p.pid()) == q );}
  if (isMeson(p)) {  return *(p.second.rbegin() + 1) == q ||*(p.second.rbegin()+2) ==q;}
  if (isDiquark(p)) {  auto i = std::find(p.second.rbegin() + 1,p.second.rbegin()+3,q); return (i!=p.second.rbegin()+3);}
  if (isBaryon(p)) { auto i = std::find(p.second.rbegin() + 1,p.second.rbegin()+4,q); return (i!=p.second.rbegin()+4);}
  if (isTetraquark(p)) { auto i = std::find(p.second.rbegin() + 1,p.second.rbegin()+5,q); return (i!=p.second.rbegin()+5);}
  if (isPentaquark(p)) { auto i = std::find(p.second.rbegin() + 1,p.second.rbegin()+6,q); return (i!=p.second.rbegin()+6);}
  if (isNucleus(p) && p.first != PROTON) { return q==3 && p(2) > 0;}
  return false;
}

hasQuark() [2/3]

template<>

bool hasQuark ( const int &  p, const int &  q  )

inline

Definition at line 485 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return hasQuark(value_digits, q);}

hasQuark() [3/3]

template<class T >

bool hasQuark ( const T &  p, const int &  q  )

inline

hasStrange()

template<class T >

bool hasStrange ( const T &  p )

inline

Definition at line 500 of file AtlasPID.h.

{ return  hasQuark(p,SQUARK); }

hasTop()

template<class T >

bool hasTop ( const T &  p )

inline

Definition at line 503 of file AtlasPID.h.

{ return  hasQuark(p,TQUARK); }

isBaryon() [1/3]

template<>

bool isBaryon ( const DecodedPID &  p )

inline

Table 43.2.

Definition at line 359 of file AtlasPID.h.

                                                    {
  if (p.ndigits() < 4 ) return false;
  if (p.max_digit(1,4) >= 6 ) return false;
  if (p.min_digit(1,4) == 0) return false;
  if (p.ndigits() == 4 && (p.last() == 2 || p.last() == 4|| p.last() == 6|| p.last() == 8) ) return true;
 
  if (p.ndigits() == 5 && p(0) == 1 &&  (p.last() == 2 || p.last() == 4) ) return true;
  if (p.ndigits() == 5 && p(0) == 3 &&  (p.last() == 2 || p.last() == 4) ) return true;
 
  if (p.ndigits() == 6 ) {
    if (p(0) == 1 && p(1) == 0 && p.last() == 2 ) return true;
    if (p(0) == 1 && p(1) == 1 && p.last() == 2 ) return true;
    if (p(0) == 1 && p(1) == 2 && p.last() == 4 ) return true;
 
    if (p(0) == 2 && p(1) == 0 && p.last() == 2 ) return true;
    if (p(0) == 2 && p(1) == 0 && p.last() == 4 ) return true;
    if (p(0) == 2 && p(1) == 1 && p.last() == 2 ) return true;
 
    if (p(0) == 1 && p(1) == 0 && p.last() == 4 ) return true;
    if (p(0) == 1 && p(1) == 0 && p.last() == 6 ) return true;
    if (p(0) == 2 && p(1) == 0 && p.last() == 6 ) return true;
    if (p(0) == 2 && p(1) == 0 && p.last() == 8 ) return true;
  }
 
  if (p.ndigits() == 5 ) {
    if (p(0) == 2 && p.last() == 2 ) return true;
    if (p(0) == 2 && p.last() == 4 ) return true;
    if (p(0) == 2 && p.last() == 6 ) return true;
    if (p(0) == 5 && p.last() == 2 ) return true;
    if (p(0) == 1 && p.last() == 6 ) return true;
    if (p(0) == 4 && p.last() == 2 ) return true;
  }
  return false;
}

isBaryon() [2/3]

template<>

bool isBaryon ( const int &  p )

inline

Definition at line 449 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isBaryon(value_digits);}

isBaryon() [3/3]

template<class T >

bool isBaryon ( const T &  p )

inline

Definition at line 220 of file AtlasPID.h.

{return isBaryon(p->pdg_id());}

isBBbarMeson() [1/3]

template<>

bool isBBbarMeson ( const DecodedPID &  p )

inline

Definition at line 524 of file AtlasPID.h.

{ return leadingQuark(p) == BQUARK && isMeson(p) && (*(p.second.rbegin()+2)) == BQUARK && (*(p.second.rbegin()+1)) == BQUARK; }

isBBbarMeson() [2/3]

template<>

bool isBBbarMeson ( const int &  p )

inline

Definition at line 525 of file AtlasPID.h.

{ return isBBbarMeson(DecodedPID(p)); }

isBBbarMeson() [3/3]

template<class T >

bool isBBbarMeson ( const T &  p )

inline

Definition at line 523 of file AtlasPID.h.

{ return isBBbarMeson(p->pdg_id());}

isBottom() [1/2]

template<>

bool isBottom ( const int &  p )

inline

Definition at line 133 of file AtlasPID.h.

{ return std::abs(p) == 5;}

isBottom() [2/2]

template<class T >

bool isBottom ( const T &  p )

inline

Definition at line 132 of file AtlasPID.h.

{return isBottom(p->pdg_id());}

isBottomBaryon()

template<class T >

bool isBottomBaryon ( const T &  p )

inline

Definition at line 532 of file AtlasPID.h.

{ return  leadingQuark(p) == BQUARK && isBaryon(p); }

isBottomHadron()

template<class T >

bool isBottomHadron ( const T &  p )

inline

Definition at line 509 of file AtlasPID.h.

{ return  leadingQuark(p) == BQUARK && isHadron(p); }

isBottomMeson()

template<class T >

bool isBottomMeson ( const T &  p )

inline

Definition at line 516 of file AtlasPID.h.

{ return  leadingQuark(p) == BQUARK && isMeson(p); }

isBSM() [1/3]

template<>

bool isBSM ( const DecodedPID &  p )

inline

APID: graviton and all Higgs extensions are BSM.

Definition at line 431 of file AtlasPID.h.

                                                 {
  if (p.pid() == GRAVITON || std::abs(p.pid()) == MAVTOP || p.pid() == DARKPHOTON) return true;
  if (std::abs(p.pid()) > 16 && std::abs(p.pid()) < 19) return true;
  if (std::abs(p.pid()) > 31 && std::abs(p.pid()) < 38) return true;
  if (std::abs(p.pid()) > 39 && std::abs(p.pid()) < 81) return true;
  if (std::abs(p.pid()) > 6 && std::abs(p.pid()) < 9) return true;
  if (isSUSY(p)) return true;
  if (isTechnicolor(p)) return true;
  if (isExcited(p)) return true;
  if (isKK(p)) return true;
  if (isHiddenValley(p)) return true;
  if (isGenericMultichargedParticle(p)) return true;
  return false;
}

isBSM() [2/3]

template<>

bool isBSM ( const int &  p )

inline

Definition at line 453 of file AtlasPID.h.

                                          {
  if (p == GRAVITON || std::abs(p) == MAVTOP || p == DARKPHOTON) return true;
  if (std::abs(p) > 16 && std::abs(p) < 19) return true;
  if (std::abs(p) > 31 && std::abs(p) < 38) return true;
  if (std::abs(p) > 39 && std::abs(p) < 81) return true;
  if (std::abs(p) > 6 && std::abs(p) < 9) return true;
  auto value_digits = DecodedPID(p); return isBSM(value_digits);
}

isBSM() [3/3]

template<class T >

bool isBSM ( const T &  p )

inline

Definition at line 224 of file AtlasPID.h.

{return isBSM(p->pdg_id());}

isCCbarMeson() [1/3]

template<>

bool isCCbarMeson ( const DecodedPID &  p )

inline

Definition at line 520 of file AtlasPID.h.

{ return leadingQuark(p) == CQUARK && isMeson(p) && (*(p.second.rbegin()+2)) == CQUARK && (*(p.second.rbegin()+1)) == CQUARK; }

isCCbarMeson() [2/3]

template<>

bool isCCbarMeson ( const int &  p )

inline

Definition at line 521 of file AtlasPID.h.

{ return isCCbarMeson(DecodedPID(p)); }

isCCbarMeson() [3/3]

template<class T >

bool isCCbarMeson ( const T &  p )

inline

Definition at line 519 of file AtlasPID.h.

{ return isCCbarMeson(p->pdg_id());}

isCharged()

template<class T >

bool isCharged ( const T &  p )

inline

Definition at line 545 of file AtlasPID.h.

{ return charge3(p) != 0;}

isCharm() [1/2]

template<>

bool isCharm ( const int &  p )

inline

Definition at line 130 of file AtlasPID.h.

{ return std::abs(p) == 4;}

isCharm() [2/2]

template<class T >

bool isCharm ( const T &  p )

inline

Definition at line 129 of file AtlasPID.h.

{return isCharm(p->pdg_id());}

isCharmBaryon()

template<class T >

bool isCharmBaryon ( const T &  p )

inline

Definition at line 531 of file AtlasPID.h.

{ return  leadingQuark(p) == CQUARK && isBaryon(p); }

isCharmHadron()

template<class T >

bool isCharmHadron ( const T &  p )

inline

Definition at line 508 of file AtlasPID.h.

{ return  leadingQuark(p) == CQUARK && isHadron(p); }

isCharmMeson()

template<class T >

bool isCharmMeson ( const T &  p )

inline

Definition at line 515 of file AtlasPID.h.

{ return  leadingQuark(p) == CQUARK && isMeson(p); }

isChLepton() [1/2]

template<>

bool isChLepton ( const int &  p )

inline

Definition at line 149 of file AtlasPID.h.

{ auto sp = std::abs(p); return sp >= 11 && sp <= 18 && sp%2 == 1; }

isChLepton() [2/2]

template<class T >

bool isChLepton ( const T &  p )

inline

APID: the fourth generation leptons are leptons.

Definition at line 148 of file AtlasPID.h.

{return isChLepton(p->pdg_id());}

isDiquark() [1/3]

template<>

bool isDiquark ( const DecodedPID &  p )

inline

PDG rule 4 Diquarks have 4-digit numbers with nq1 >= nq2 and nq3 = 0 APID: the diquarks with fourth generation are not diquarks.

Definition at line 318 of file AtlasPID.h.

                                                     {
  if ( p.ndigits() == 4 &&  p(0) >= p(1) && p(2) == 0 &&  p.last() % 2 == 1
       && p.max_digit(1,3) <= 6
       ) return true;
  return false;
}

isDiquark() [2/3]

template<>

bool isDiquark ( const int &  p )

inline

Definition at line 447 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isDiquark(value_digits);}

isDiquark() [3/3]

template<class T >

bool isDiquark ( const T &  p )

inline

Definition at line 217 of file AtlasPID.h.

{return isDiquark(p->pdg_id());}

isDM() [1/2]

template<>

bool isDM ( const int &  p )

inline

Definition at line 190 of file AtlasPID.h.

{ auto sp = std::abs(p); return (sp >= 51 && sp <= 60) || sp == DARKPHOTON; }

isDM() [2/2]

template<class T >

bool isDM ( const T &  p )

inline

PDG rule 11j: The nature of Dark Matter (DM) is not known, and therefore a definitive classificationis too early.

Candidates within specific scenarios are classified therein, such as 1000022 for the lightest neutralino. Generic fundamental states can be given temporary codes in the range 51 - 60, with 51, 52 and 53 reserved for spin 0, 1/2 and 1 ones (this could also be an axion state). Generic mediators of s-channel DM pair creation of annihilationcan be given codes 54 and 55 for spin 0 or 1 ones. Separate antiparticles, with negativecodes, may or may not exist. More elaborate new scenarios should be constructed with n= 5 and nr = 9. APID: Only the 51-60 range is considered DM. The antiparticles are assumed to be existing.

Definition at line 189 of file AtlasPID.h.

{return isDM(p->pdg_id());}

isElectron() [1/2]

template<>

bool isElectron ( const int &  p )

inline

Definition at line 152 of file AtlasPID.h.

{ return std::abs(p) == ELECTRON;}

isElectron() [2/2]

template<class T >

bool isElectron ( const T &  p )

inline

Definition at line 151 of file AtlasPID.h.

{return isElectron(p->pdg_id());}

isEMInteracting() [1/2]

template<>

bool isEMInteracting ( const int &  p )

inline

Definition at line 609 of file AtlasPID.h.

{return (isPhoton(p) || isZ(p) || std::abs(charge(p))>std::numeric_limits<double>::epsilon() || isMonopole(p));}

isEMInteracting() [2/2]

template<class T >

bool isEMInteracting ( const T &  p )

inline

Definition at line 608 of file AtlasPID.h.

{return isEMInteracting(p->pdg_id());}

isExcited() [1/3]

template<>

bool isExcited ( const DecodedPID &  p )

inline

PDG rule 11f Excited (composite) quarks and leptons are identified by setting n= 4 and nr= 0.

Definition at line 276 of file AtlasPID.h.

{return (p.ndigits() == 7 && (p(0) == 4 && p(1) == 0 ) && (isLepton(p.shift(2))||isQuark(p.shift(2))) );}

isExcited() [2/3]

template<>

bool isExcited ( const int &  p )

inline

Definition at line 277 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isExcited(value_digits);}

isExcited() [3/3]

template<class T >

bool isExcited ( const T &  p )

inline

Definition at line 212 of file AtlasPID.h.

{return isExcited(p->pdg_id());}

isGeantino() [1/2]

template<>

bool isGeantino ( const int &  p )

inline

Definition at line 245 of file AtlasPID.h.

{ return (std::abs(p) ==  GEANTINO0 || std::abs(p) ==  GEANTINOPLUS);}

isGeantino() [2/2]

template<class T >

bool isGeantino ( const T &  p )

inline

Definition at line 229 of file AtlasPID.h.

{return isGeantino(p->pdg_id());}

isGenericMultichargedParticle() [1/3]

template<>

bool isGenericMultichargedParticle ( const DecodedPID &  p )

inline

In addition, there is a need to identify ”Q-ball” and similar very exotic (multi-charged) particles which may have large, non-integer charge.

These particles are assigned the ad-hoc numbering +/-100XXXY0, where the charge is XXX.Y. or +/-200XXYY0, where the charge is XX/YY. The case of +/-200XXYY0 is legacy, see https://gitlab.cern.ch/atlas/athena/-/merge_requests/25862 Note that no other quantum numbers besides the charge are considered for these generic multi-charged particles (e.g. isSUSY() is false for them). Such a model was used in previous Run-1 (1301.5272,1504.04188) and Run-2 (1812.03673,2303.13613) ATLAS searches.

Definition at line 312 of file AtlasPID.h.

{return (p.ndigits() == 8 && (p(0) == 1 || p(0) == 2) && p(1) == 0 && p(2) == 0 && p(7) == 0);}

isGenericMultichargedParticle() [2/3]

template<>

bool isGenericMultichargedParticle ( const int &  p )

inline

Definition at line 313 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isGenericMultichargedParticle(value_digits);}

isGenericMultichargedParticle() [3/3]

template<class T >

bool isGenericMultichargedParticle ( const T &  p )

inline

Definition at line 216 of file AtlasPID.h.

{return isGenericMultichargedParticle(p->pdg_id());}

isGenSpecific() [1/2]

template<>

bool isGenSpecific ( const int &  p )

inline

Main Table for MC internal use 81–100,901–930,998-999,1901–1930,2901–2930, and 3901–3930.

Definition at line 235 of file AtlasPID.h.

                                                  {
  if (p >= 81 && p <= 100) return true;
  if (p >= 901 && p <= 930) return true;
  if (p >= 998 && p <= 999) return true;
  if (p >= 1901 && p <= 1930) return true;
  if (p >= 2901 && p <= 2930) return true;
  if (p >= 3901 && p <= 3930) return true;
  return false;
}

isGenSpecific() [2/2]

template<class T >

bool isGenSpecific ( const T &  p )

inline

Definition at line 228 of file AtlasPID.h.

{return isGenSpecific(p->pdg_id());}

isGlueball() [1/3]

template<>

bool isGlueball ( const DecodedPID &  p )

inline

Definition at line 249 of file AtlasPID.h.

                                                       {
  if (p.ndigits() > 4) return false; // APID avoid classifying R-Glueballs as SM Glueballs
  size_t ng = 0;
  for (size_t i = 1; i + 1 < p.ndigits(); ++i) {
    if (p(i) == COMPOSITEGLUON) ng++;
  }
  return  (*(p.second.rbegin()+2)) == COMPOSITEGLUON && (*(p.second.rbegin()+1)) == COMPOSITEGLUON && ng > 0;
}

isGlueball() [2/3]

template<>

bool isGlueball ( const int &  p )

inline

Definition at line 257 of file AtlasPID.h.

{  auto value_digits = DecodedPID(p); return isGlueball(value_digits); }

isGlueball() [3/3]

template<class T >

bool isGlueball ( const T &  p )

inline

APID: Definition of Glueballs: SM glueballs 99X (X=1,3,5), 999Y (Y=3,5,7)

Definition at line 248 of file AtlasPID.h.

{ return isGlueball(p->pdg_id()); }

isGluon() [1/2]

template<>

bool isGluon ( const int &  p )

inline

Definition at line 168 of file AtlasPID.h.

{ return p == GLUON; }

isGluon() [2/2]

template<class T >

bool isGluon ( const T &  p )

inline

Definition at line 167 of file AtlasPID.h.

{return isGluon(p->pdg_id());}

isGraviton() [1/2]

template<>

bool isGraviton ( const int &  p )

inline

Definition at line 205 of file AtlasPID.h.

{ return p == GRAVITON; }

isGraviton() [2/2]

template<class T >

bool isGraviton ( const T &  p )

inline

Definition at line 204 of file AtlasPID.h.

{return isGraviton(p->pdg_id());}

isHadron() [1/3]

template<>

bool isHadron ( const DecodedPID &  p )

inline

Definition at line 462 of file AtlasPID.h.

{ return isMeson(p)||isBaryon(p)||isTetraquark(p)||isPentaquark(p);}

isHadron() [2/3]

template<>

bool isHadron ( const int &  p )

inline

Definition at line 463 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isHadron(value_digits);}

isHadron() [3/3]

template<class T >

bool isHadron ( const T &  p )

inline

Definition at line 218 of file AtlasPID.h.

{return isHadron(p->pdg_id());}

isHeavyBaryon()

template<class T >

bool isHeavyBaryon ( const T &  p )

inline

Definition at line 529 of file AtlasPID.h.

{  auto lq = leadingQuark(p); return  (lq == CQUARK || lq == BQUARK || lq == TQUARK) && isBaryon(p); }

isHeavyHadron()

template<class T >

bool isHeavyHadron ( const T &  p )

inline

Definition at line 506 of file AtlasPID.h.

{  auto lq = leadingQuark(p); return  (lq == CQUARK || lq == BQUARK || lq == TQUARK ) && isHadron(p); }

isHeavyMeson()

template<class T >

bool isHeavyMeson ( const T &  p )

inline

Definition at line 513 of file AtlasPID.h.

{ auto lq = leadingQuark(p); return  (lq == CQUARK || lq == BQUARK || lq == TQUARK) && isMeson(p); }

isHiddenValley() [1/3]

template<>

bool isHiddenValley ( const DecodedPID &  p )

inline

PDG rule 11k Hidden Valley particles have n = 4 and n_r = 9, and trailing numbers in agreement with their nearest-analog standard particles, as far as possible.

Thus 4900021 is the gauge boson gv of a confining gauge field, 490000 nqv and 490001 nlv fundamental constituents charged or not under this, 4900022 is the γv of a non-confining field, and 4900 nqv1 nqv2 nJ a Hidden Valley meson.

Definition at line 303 of file AtlasPID.h.

{return (p.ndigits() == 7 &&  p(0) == 4 && p(1) == 9 && isValid(p.shift(2)));}

isHiddenValley() [2/3]

template<>

bool isHiddenValley ( const int &  p )

inline

Definition at line 304 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isHiddenValley(value_digits);}

isHiddenValley() [3/3]

template<class T >

bool isHiddenValley ( const T &  p )

inline

Definition at line 215 of file AtlasPID.h.

{return isHiddenValley(p->pdg_id());}

isHiggs() [1/2]

template<>

bool isHiggs ( const int &  p )

inline

Definition at line 200 of file AtlasPID.h.

{ return p == HIGGSBOSON; }

isHiggs() [2/2]

template<class T >

bool isHiggs ( const T &  p )

inline

APID: HIGGS boson is only one particle.

Definition at line 199 of file AtlasPID.h.

{return isHiggs(p->pdg_id());}

isKK() [1/3]

template<>

bool isKK ( const DecodedPID &  p )

inline

PDG rule 11h A black hole in models with extra dimensions has code 5000040.

Kaluza-Klein excitations in models with extra dimensions have n = 5 or n = 6, to distinguish excitations of left-or right-handed fermions or, in case of mixing, the lighter or heavier state (cf. 11d). The non zero nr digit gives the radial excitation number, in scenarios where the level spacing allows these to be distinguished. Should the model also contain supersymmetry, excited SUSY states would be denoted by a nn_r > 0, with n = 1 or 2 as usual. Should some colored states be long-lived enough that hadrons would form around them, the coding strategy of 11g applies, with the initial two nnr digits preserved in the combined code.

Definition at line 286 of file AtlasPID.h.

{return (p.ndigits() == 7 && (p(0) == 5 || p(0) == 6 ) );}

isKK() [2/3]

template<>

bool isKK ( const int &  p )

inline

Definition at line 287 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isKK(value_digits);}

isKK() [3/3]

template<class T >

bool isKK ( const T &  p )

inline

Definition at line 213 of file AtlasPID.h.

{return isKK(p->pdg_id());}

isLepton() [1/3]

template<>

bool isLepton ( const DecodedPID &  p )

inline

Definition at line 141 of file AtlasPID.h.

{ return isLepton(p.pid()); }

isLepton() [2/3]

template<>

bool isLepton ( const int &  p )

inline

Definition at line 140 of file AtlasPID.h.

{ auto sp = std::abs(p); return sp >= 11 && sp <= 18; }

isLepton() [3/3]

template<class T >

bool isLepton ( const T &  p )

inline

APID: the fourth generation leptons are leptons.

Definition at line 139 of file AtlasPID.h.

{return isLepton(p->pdg_id());}

isLeptoQuark() [1/2]

template<>

bool isLeptoQuark ( const int &  p )

inline

Definition at line 208 of file AtlasPID.h.

{ return std::abs(p) == LEPTOQUARK; }

isLeptoQuark() [2/2]

template<class T >

bool isLeptoQuark ( const T &  p )

inline

Definition at line 207 of file AtlasPID.h.

{return isLeptoQuark(p->pdg_id());}

isLightBaryon()

template<class T >

bool isLightBaryon ( const T &  p )

inline

Definition at line 528 of file AtlasPID.h.

{ auto lq = leadingQuark(p); return  (lq == DQUARK || lq == UQUARK||lq == SQUARK) && isBaryon(p); }

isLightHadron()

template<class T >

bool isLightHadron ( const T &  p )

inline

Definition at line 505 of file AtlasPID.h.

{ auto lq = leadingQuark(p); return  (lq == DQUARK || lq == UQUARK||lq == SQUARK) && isHadron(p); }

isLightMeson()

template<class T >

bool isLightMeson ( const T &  p )

inline

Definition at line 512 of file AtlasPID.h.

{ auto lq = leadingQuark(p); return  (lq == DQUARK || lq == UQUARK||lq == SQUARK) && isMeson(p); }

isMeson() [1/3]

template<>

bool isMeson ( const DecodedPID &  p )

inline

Table 43.1 PDG rule 5a: The numbers specifying the meson’s quark content conform to the convention nq1= 0 and nq2 >= nq3.

The special case K0L is the sole exception to this rule. PDG rule 5C: The special numbers 310 and 130 are given to the K0S and K0L respectively. APID: The special code K0 is used when a generator uses K0S/K0L

Definition at line 333 of file AtlasPID.h.

                                                   {
  if (p.ndigits() < 3 ) return false;
  if (std::abs(p.pid()) == K0S) return true;
  if (std::abs(p.pid()) == K0L) return true;
  if (std::abs(p.pid()) == K0) return true;
  if (p.last() % 2 != 1 ) return false;
  if (p.max_digit(1,3) >= 6 ) return false;
  if (p.max_digit(1,3) == 0 ) return false;
 
  if (p.ndigits() == 3 && p(0) == p(1) && p.pid() < 0 ) return false;
  if (p.ndigits() == 5 && p(2) == p(3) && p.pid() < 0 ) return false;
  if (p.ndigits() == 7 && p(4) == p(5) && p.pid() < 0 ) return false;
 
 
  if (p.ndigits() == 3 && p(0) >= p(1) && p(1) != 0  ) return true;
  if (p.ndigits() == 5 && p(2) >= p(3) && p(3) != 0 && p(0) == 1 && p(1) == 0) return true;
  if (p.ndigits() == 5 && p(2) >= p(3) && p(3) != 0 && p(0) == 2 && p(1) == 0 && p.last() > 1 ) return true;
  if (p.ndigits() == 5 && p(2) >= p(3) && p(3) != 0 && p(0) == 3 && p(1) == 0 && p.last() > 1 ) return true;
 
  if (p.ndigits() == 6 && p(3) >= p(4) && p(4) != 0 && p.last() % 2 == 1  ) return true;
 
  if (p.ndigits() == 7 && p(4) >= p(5) && p(5) != 0) return true;
 
  return false;
}

isMeson() [2/3]

template<>

bool isMeson ( const int &  p )

inline

Definition at line 448 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isMeson(value_digits);}

isMeson() [3/3]

template<class T >

bool isMeson ( const T &  p )

inline

Definition at line 219 of file AtlasPID.h.

{return isMeson(p->pdg_id());}

isMonopole() [1/3]

template<>

bool isMonopole ( const DecodedPID &  p )

inline

PDG rule 11i Magnetic monopoles and dyons are assumed to have one unit of Dirac monopole charge and a variable integer number nq1nq2 nq3 units of electric charge.

Codes 411nq1nq2 nq3 0 are then used when the magnetic and electrical charge sign agree and 412nq1nq2 nq3 0 when they disagree, with the overall sign of the particle set by the magnetic charge. For now no spin information is provided.

Definition at line 295 of file AtlasPID.h.

{return (p.ndigits() == 7 && p(0) == 4 && p(1) == 1  && (p(2) == 1 || p(2) == 2 ) && p(6) == 0);}

isMonopole() [2/3]

template<>

bool isMonopole ( const int &  p )

inline

Definition at line 296 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isMonopole(value_digits);}

isMonopole() [3/3]

template<class T >

bool isMonopole ( const T &  p )

inline

Definition at line 214 of file AtlasPID.h.

{return isMonopole(p->pdg_id());}

isMuon() [1/2]

template<>

bool isMuon ( const int &  p )

inline

Definition at line 155 of file AtlasPID.h.

{ return std::abs(p) == MUON;}

isMuon() [2/2]

template<class T >

bool isMuon ( const T &  p )

inline

Definition at line 154 of file AtlasPID.h.

{return isMuon(p->pdg_id());}

isNeutral()

template<class T >

bool isNeutral ( const T &  p )

inline

Definition at line 546 of file AtlasPID.h.

{ return charge3(p) == 0;}

isNeutrino() [1/2]

template<>

bool isNeutrino ( const int &  p )

inline

Definition at line 162 of file AtlasPID.h.

{ auto sp = std::abs(p); return sp == NU_E || sp == NU_MU || sp == NU_TAU || sp == 18;  }

isNeutrino() [2/2]

template<class T >

bool isNeutrino ( const T &  p )

inline

APID: the fourth generation neutrinos are neutrinos.

Definition at line 161 of file AtlasPID.h.

{return isNeutrino(p->pdg_id());}

isNucleus() [1/3]

template<>

bool isNucleus ( const DecodedPID &  p )

inline

PDG rule 16: Nuclear codes are given as 10-digit numbers±10LZZZAAAI.

For a (hyper)nucleus consistingofnpprotons,nnneutrons andnΛΛ’s,A=np+nn+nΛgives the total baryon number, Z=np the total charge andL=nΛthe total number of strange quarks.Igives the isomerlevel, withI= 0corresponding to the ground state andI >0to excitations, see [2], wherestates denotedm,n,p,qtranslate toI= 1–4. As examples, the deuteron is 1000010020 and 235U is

1. To avoid ambiguities, nuclear codes should not be applied to a singlehadron, like p,n or Λ0, where quark-contents-based codes already exist.

Definition at line 426 of file AtlasPID.h.

                                                     {
  if (std::abs(p.pid()) == PROTON) return true;
  return (p.ndigits() == 10 &&  p(0) == 1 &&  p(1) == 0 );
}

isNucleus() [2/3]

template<>

bool isNucleus ( const int &  p )

inline

Definition at line 452 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isNucleus(value_digits);}

isNucleus() [3/3]

template<class T >

bool isNucleus ( const T &  p )

inline

Definition at line 223 of file AtlasPID.h.

{return isNucleus(p->pdg_id());}

isParton()

template<class T >

bool isParton ( const T &  p )

inline

Definition at line 611 of file AtlasPID.h.

{ return isQuark(p)||isGluon(p);}

isPentaquark() [1/3]

template<>

bool isPentaquark ( const DecodedPID &  p )

inline

PDG rule 15 The 9-digit penta-quark codes are±1nrnLnq1nq2nq3nq4nq5nJ, sorted such thatnq1≥nq2≥nq3≥nq4.

In the particle the first four are quarks and the fifth an antiquark while t heopposite holds in the antiparticle, which is given with a negative sign. Thenr,nL, andnJnumbers have the same meaning as for ordinary hadrons.

Definition at line 399 of file AtlasPID.h.

                                                        {
  return (p.ndigits() == 9 && p(0) == 1 &&
          p.max_digit(1,6) <= 6  && p.min_digit(1,6) > 0 &&
          ( p(3) >= p(4) && p(4) >= p(5) && p(5) >= p(6)) );
}

isPentaquark() [2/3]

template<>

bool isPentaquark ( const int &  p )

inline

Definition at line 451 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isPentaquark(value_digits);}

isPentaquark() [3/3]

template<class T >

bool isPentaquark ( const T &  p )

inline

Definition at line 222 of file AtlasPID.h.

{return isPentaquark(p->pdg_id());}

isPhoton() [1/2]

template<>

bool isPhoton ( const int &  p )

inline

Definition at line 171 of file AtlasPID.h.

{ return p == PHOTON; }

isPhoton() [2/2]

template<class T >

bool isPhoton ( const T &  p )

inline

Definition at line 170 of file AtlasPID.h.

{return isPhoton(p->pdg_id());}

isPythia8Specific() [1/3]

template<>

bool isPythia8Specific ( const DecodedPID &  p )

inline

Definition at line 231 of file AtlasPID.h.

{ return (p.ndigits() == 7 && p(0) == 9 && p(1) == 9);}

isPythia8Specific() [2/3]

template<>

bool isPythia8Specific ( const int &  p )

inline

Definition at line 232 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isPythia8Specific(value_digits);}

isPythia8Specific() [3/3]

template<class T >

bool isPythia8Specific ( const T &  p )

inline

Definition at line 230 of file AtlasPID.h.

{return isPythia8Specific(p->pdg_id());}

isQuark() [1/3]

template<>

bool isQuark ( const DecodedPID &  p )

inline

Definition at line 124 of file AtlasPID.h.

{ return isQuark(p.pid()); }

isQuark() [2/3]

template<>

bool isQuark ( const int &  p )

inline

Definition at line 123 of file AtlasPID.h.

{ return p != 0 && (std::abs(p) <= 8 || std::abs(p) == MAVTOP);}

isQuark() [3/3]

template<class T >

bool isQuark ( const T &  p )

inline

PDG rule 2: Quarks and leptons are numbered consecutively starting from 1 and 11 respectively; to dothis they are first ordered by family and within families by weak isospin.

APID: the fourth generation quarks are quarks.

Definition at line 122 of file AtlasPID.h.

{return isQuark(p->pdg_id());}

isResonance()

template<class T >

bool isResonance ( const T &  p )

inline

Definition at line 202 of file AtlasPID.h.

{ return isZ(p) || isW(p) || isHiggs(p) || isTop(p); } // APID: not including t' (pdg_id=8)

isSMLepton() [1/3]

template<>

bool isSMLepton ( const DecodedPID &  p )

inline

Definition at line 145 of file AtlasPID.h.

{ return isSMLepton(p.pid()); }

isSMLepton() [2/3]

template<>

bool isSMLepton ( const int &  p )

inline

Definition at line 144 of file AtlasPID.h.

{ auto sp = std::abs(p); return sp >= 11 && sp <= 16; }

isSMLepton() [3/3]

template<class T >

bool isSMLepton ( const T &  p )

inline

Definition at line 143 of file AtlasPID.h.

{return isSMLepton(p->pdg_id());}

isSMNeutrino() [1/2]

template<>

bool isSMNeutrino ( const int &  p )

inline

Definition at line 165 of file AtlasPID.h.

{ auto sp = std::abs(p); return sp == NU_E || sp == NU_MU || sp == NU_TAU;  }

isSMNeutrino() [2/2]

template<class T >

bool isSMNeutrino ( const T &  p )

inline

Definition at line 164 of file AtlasPID.h.

{return isSMNeutrino(p->pdg_id());}

isStrange() [1/2]

template<>

bool isStrange ( const int &  p )

inline

Definition at line 127 of file AtlasPID.h.

{ return std::abs(p) == 3;}

isStrange() [2/2]

template<class T >

bool isStrange ( const T &  p )

inline

Definition at line 126 of file AtlasPID.h.

{return isStrange(p->pdg_id());}

isStrangeBaryon()

template<class T >

bool isStrangeBaryon ( const T &  p )

inline

Definition at line 530 of file AtlasPID.h.

{ return  leadingQuark(p) == SQUARK && isBaryon(p); }

isStrangeHadron()

template<class T >

bool isStrangeHadron ( const T &  p )

inline

Definition at line 507 of file AtlasPID.h.

{ return  leadingQuark(p) == SQUARK && isHadron(p); }

isStrangeMeson()

template<class T >

bool isStrangeMeson ( const T &  p )

inline

Definition at line 514 of file AtlasPID.h.

{ return  leadingQuark(p) == SQUARK && isMeson(p); }

isStrongInteracting() [1/2]

template<>

bool isStrongInteracting ( const int &  p )

inline

Definition at line 698 of file AtlasPID.h.

{ return (isGluon(p) || isQuark(p) || isDiquark(p) || isGlueball(p) || isLeptoQuark(p) || isHadron(p) || SUSY::isRHadron(p));} // APID: Glueballs and R-Hadrons are also strong-interacting

isStrongInteracting() [2/2]

template<class T >

bool isStrongInteracting ( const T &  p )

inline

Definition at line 697 of file AtlasPID.h.

{return isStrongInteracting(p->pdg_id());}

isSUSY() [1/3]

template<>

bool isSUSY ( const DecodedPID &  p )

inline

PDG rule 11d Fundamental supersymmetric particles are identified by adding a nonzero n to the particle number.

The superpartner of a boson or a left-handed fermion has n = 1 while the superpartner of a right-handed fermion has n = 2. When mixing occurs, such as between the winos and charged Higgsinos to give charginos, or between left and right sfermions, the lighter physical state is given the smaller basis state number.

Definition at line 264 of file AtlasPID.h.

{return (p.ndigits() == 7 && (p(0) == 1 || p(0) == 2 ) && ((isValid(p.shift(2)) && !isGenSpecific(p.shift(2).pid())) || isGlueball(p.shift(2))) );}

isSUSY() [2/3]

template<>

bool isSUSY ( const int &  p )

inline

Definition at line 446 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isSUSY(value_digits);}

isSUSY() [3/3]

template<class T >

bool isSUSY ( const T &  p )

inline

Definition at line 210 of file AtlasPID.h.

{return isSUSY(p->pdg_id());}

isTau() [1/2]

template<>

bool isTau ( const int &  p )

inline

Definition at line 158 of file AtlasPID.h.

{ return std::abs(p) == TAU;}

isTau() [2/2]

template<class T >

bool isTau ( const T &  p )

inline

Definition at line 157 of file AtlasPID.h.

{return isTau(p->pdg_id());}

isTechnicolor() [1/3]

template<>

bool isTechnicolor ( const DecodedPID &  p )

inline

PDG rule 11e Technicolor states have n = 3, with technifermions treated like ordinary fermions.

States which are ordinary color singlets have n_r = 0. Color octets have n_r = 1. If a state has non-trivial quantum numbers under the topcolor groups SU(3)1×SU(3)2, the quantum numbers are specified by tech, ij, where i and j are 1 or 2. nLis then 2i+j. The coloron V8, is a heavy gluon color octet and thus is 3100021

Definition at line 271 of file AtlasPID.h.

{return (p.ndigits() == 7 &&  p(0) == 3 && (p(1) == 0 || p(0) == 1) && isValid(p.shift(2)) && !isGenSpecific(p.shift(2).pid()));}

isTechnicolor() [2/3]

template<>

bool isTechnicolor ( const int &  p )

inline

Definition at line 272 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isTechnicolor(value_digits);}

isTechnicolor() [3/3]

template<class T >

bool isTechnicolor ( const T &  p )

inline

Definition at line 211 of file AtlasPID.h.

{return isTechnicolor(p->pdg_id());}

isTetraquark() [1/3]

template<>

bool isTetraquark ( const DecodedPID &  p )

inline

PDG rule 14 The 9-digit tetra-quark codes are±1nrnLnq1nq20nq3nq4nJ.

For the particleq1q2is a diquarkand ̄q3 ̄q4an antidiquark, sorted such thatnq1≥nq2,nq3≥nq4,nq1≥nq3, andnq2≥nq4ifnq1=nq3. For the antiparticle, given with a negative sign, ̄q1 ̄q2is an antidiquark andq3q4a diquark, with the same sorting except that eithernq1> nq3ornq2> nq4(so thatflavour-diagonal states are particles). Thenr,nL, andnJnumbers have the same meaningas for ordinary hadrons.

Definition at line 410 of file AtlasPID.h.

                                                        {
  return (p.ndigits() == 9 && p(0) == 1 && p(5) == 0 &&
          p.max_digit(1,3) <= 6  && p.min_digit(1,3) > 0 &&
          p.max_digit(1+3,3+3) <= 6  && p.min_digit(1+3,3+3) > 0 &&
          ( p(3) >= p(4)  && p(6) >= p(7) ) &&  ( ( p(3) > p(6) ) || ( p(3) == p(6) && (p(4) >= p(7))))
          );
}

isTetraquark() [2/3]

template<>

bool isTetraquark ( const int &  p )

inline

Definition at line 450 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isTetraquark(value_digits);}

isTetraquark() [3/3]

template<class T >

bool isTetraquark ( const T &  p )

inline

Definition at line 221 of file AtlasPID.h.

{return isTetraquark(p->pdg_id());}

isTop() [1/2]

template<>

bool isTop ( const int &  p )

inline

Definition at line 136 of file AtlasPID.h.

{ return std::abs(p) == 6;}

isTop() [2/2]

template<class T >

bool isTop ( const T &  p )

inline

Definition at line 135 of file AtlasPID.h.

{return isTop(p->pdg_id());}

isTopBaryon()

template<class T >

bool isTopBaryon ( const T &  p )

inline

Definition at line 533 of file AtlasPID.h.

{ return  leadingQuark(p) == TQUARK && isBaryon(p); }

isTopHadron()

template<class T >

bool isTopHadron ( const T &  p )

inline

Definition at line 510 of file AtlasPID.h.

{ return  leadingQuark(p) == TQUARK && isHadron(p); }

isTopMeson()

template<class T >

bool isTopMeson ( const T &  p )

inline

Definition at line 517 of file AtlasPID.h.

{ return  leadingQuark(p) == TQUARK && isMeson(p); }

isTrajectory() [1/2]

template<>

bool isTrajectory ( const int &  p )

inline

Definition at line 196 of file AtlasPID.h.

{ return std::abs(p) == POMERON || std::abs(p) == ODDERON || std::abs(p) == REGGEON; }

isTrajectory() [2/2]

template<class T >

bool isTrajectory ( const T &  p )

inline

PDG rule 8: The pomeron and odderon trajectories and a generic reggeon trajectory of states in QCD areassigned codes 990, 9990, and 110 respectively.

Definition at line 195 of file AtlasPID.h.

{return isTrajectory(p->pdg_id());}

isTransportable() [1/3]

template<>

bool isTransportable ( const DecodedPID &  p )

inline

Definition at line 464 of file AtlasPID.h.

{ return isPhoton(p.pid()) || isGeantino(p.pid()) || isHadron(p) || isLepton(p.pid()) || p.pid() == DARKPHOTON;}

isTransportable() [2/3]

template<>

bool isTransportable ( const int &  p )

inline

Definition at line 465 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return isTransportable(value_digits);}

isTransportable() [3/3]

template<class T >

bool isTransportable ( const T &  p )

inline

Definition at line 227 of file AtlasPID.h.

{return isTransportable(p->pdg_id());}

isValid() [1/3]

template<>

bool isValid ( const DecodedPID &  p )

inline

Av: we implement here an ATLAS-sepcific convention: all particles which are 99xxxxx are fine.

Definition at line 467 of file AtlasPID.h.

{ return isHadron(p) || isTrajectory(p.pid()) || isDiquark(p) || isBSM(p) || isNucleus(p) || (std::abs(p.pid()) < 42) || isGenSpecific(p.pid()) || isGeantino(p.pid()) || isPythia8Specific(p) || isGlueball(p); }

isValid() [2/3]

template<>

bool isValid ( const int &  p )

inline

Definition at line 468 of file AtlasPID.h.

                                            { if (!p) return false; if (std::abs(p) < 42) return true;
  if (isGenSpecific(p)) return true;
  auto value_digits = DecodedPID(p); return isValid(value_digits);
}

isValid() [3/3]

template<class T >

bool isValid ( const T &  p )

inline

Definition at line 225 of file AtlasPID.h.

{return isValid(p->pdg_id());}

isW() [1/2]

template<>

bool isW ( const int &  p )

inline

Definition at line 177 of file AtlasPID.h.

{ return std::abs(p) == WPLUSBOSON; }

isW() [2/2]

template<class T >

bool isW ( const T &  p )

inline

Definition at line 176 of file AtlasPID.h.

{return isW(p->pdg_id());}

isZ() [1/2]

template<>

bool isZ ( const int &  p )

inline

Definition at line 174 of file AtlasPID.h.

{ return p == Z0BOSON; }

isZ() [2/2]

template<class T >

bool isZ ( const T &  p )

inline

Definition at line 173 of file AtlasPID.h.

{return isZ(p->pdg_id());}

leadingQuark() [1/3]

template<>

int leadingQuark ( const DecodedPID &  p )

inline

Definition at line 488 of file AtlasPID.h.

                                                       {
  if (isQuark(p.pid())) { return std::abs(p.pid());}
  if (isMeson(p)) {return p.max_digit(1,3);}
  if (isDiquark(p)) {return p.max_digit(1,3);}
  if (isBaryon(p)) {return p.max_digit(1,4); }
  if (isTetraquark(p)) { return p.max_digit(1,5);}
  if (isPentaquark(p)) { return p.max_digit(1,6);}
  return 0;
}

leadingQuark() [2/3]

template<>

int leadingQuark ( const int &  p )

inline

Definition at line 498 of file AtlasPID.h.

{ auto value_digits = DecodedPID(p); return leadingQuark(value_digits);}

leadingQuark() [3/3]

template<class T >

int leadingQuark ( const T &  p )

inline

Definition at line 487 of file AtlasPID.h.

{return leadingQuark(p->pdg_id());}

threeCharge()

template<class T >

double threeCharge ( const T &  p )

inline

Definition at line 544 of file AtlasPID.h.

{ return charge3(p);}