MCTruthPartClassifier Namespace Reference

Classes
struct	ParticleDef
class	Info

Enumerations
enum	ParticleType : unsigned int {   Unknown = 0 , UnknownElectron = 1 , IsoElectron = 2 , NonIsoElectron = 3 ,   BkgElectron = 4 , UnknownMuon = 5 , IsoMuon = 6 , NonIsoMuon = 7 ,   BkgMuon = 8 , UnknownTau = 9 , IsoTau = 10 , NonIsoTau = 11 ,   BkgTau = 12 , UnknownPhoton = 13 , IsoPhoton = 14 , NonIsoPhoton = 15 ,   BkgPhoton = 16 , Hadron = 17 , Neutrino = 18 , NuclFrag = 19 ,   NonPrimary = 20 , GenParticle = 21 , SUSYParticle = 22 , OtherBSMParticle = 39 ,   BBbarMesonPart = 23 , BottomMesonPart = 24 , CCbarMesonPart = 25 , CharmedMesonPart = 26 ,   BottomBaryonPart = 27 , CharmedBaryonPart = 28 , StrangeBaryonPart = 29 , LightBaryonPart = 30 ,   StrangeMesonPart = 31 , LightMesonPart = 32 , BJet = 33 , CJet = 34 ,   LJet = 35 , GJet = 36 , TauJet = 37 , UnknownJet = 38 ,   Unknown = 0 , UnknownElectron = 1 , IsoElectron = 2 , NonIsoElectron = 3 ,   BkgElectron = 4 , UnknownMuon = 5 , IsoMuon = 6 , NonIsoMuon = 7 ,   BkgMuon = 8 , UnknownTau = 9 , IsoTau = 10 , NonIsoTau = 11 ,   BkgTau = 12 , UnknownPhoton = 13 , IsoPhoton = 14 , NonIsoPhoton = 15 ,   BkgPhoton = 16 , Hadron = 17 , Neutrino = 18 , NuclFrag = 19 ,   NonPrimary = 20 , GenParticle = 21 , SUSYParticle = 22 , OtherBSMParticle = 39 ,   BBbarMesonPart = 23 , BottomMesonPart = 24 , CCbarMesonPart = 25 , CharmedMesonPart = 26 ,   BottomBaryonPart = 27 , CharmedBaryonPart = 28 , StrangeBaryonPart = 29 , LightBaryonPart = 30 ,   StrangeMesonPart = 31 , LightMesonPart = 32 , BJet = 33 , CJet = 34 ,   LJet = 35 , GJet = 36 , TauJet = 37 , UnknownJet = 38 }
enum	ParticleOrigin : unsigned int {   NonDefined = 0 , SingleElec = 1 , SingleMuon = 2 , SinglePhot = 3 ,   SingleTau = 4 , PhotonConv = 5 , DalitzDec = 6 , ElMagProc = 7 ,   Mu = 8 , TauLep = 9 , top = 10 , QuarkWeakDec = 11 ,   WBoson = 12 , ZBoson = 13 , Higgs = 14 , HiggsMSSM = 15 ,   HeavyBoson = 16 , WBosonLRSM = 17 , NuREle = 18 , NuRMu = 19 ,   NuRTau = 20 , LQ = 21 , SUSY = 22 , OtherBSM = 46 ,   LightMeson = 23 , StrangeMeson = 24 , CharmedMeson = 25 , BottomMeson = 26 ,   CCbarMeson = 27 , JPsi = 28 , BBbarMeson = 29 , LightBaryon = 30 ,   StrangeBaryon = 31 , CharmedBaryon = 32 , BottomBaryon = 33 , PionDecay = 34 ,   KaonDecay = 35 , BremPhot = 36 , PromptPhot = 37 , UndrPhot = 38 ,   ISRPhot = 39 , FSRPhot = 40 , NucReact = 41 , PiZero = 42 ,   DiBoson = 43 , ZorHeavyBoson = 44 , MultiBoson = 47 , QCD = 45 ,   NonDefined = 0 , SingleElec = 1 , SingleMuon = 2 , SinglePhot = 3 ,   SingleTau = 4 , PhotonConv = 5 , DalitzDec = 6 , ElMagProc = 7 ,   Mu = 8 , TauLep = 9 , top = 10 , QuarkWeakDec = 11 ,   WBoson = 12 , ZBoson = 13 , Higgs = 14 , HiggsMSSM = 15 ,   HeavyBoson = 16 , WBosonLRSM = 17 , NuREle = 18 , NuRMu = 19 ,   NuRTau = 20 , LQ = 21 , SUSY = 22 , OtherBSM = 46 ,   LightMeson = 23 , StrangeMeson = 24 , CharmedMeson = 25 , BottomMeson = 26 ,   CCbarMeson = 27 , JPsi = 28 , BBbarMeson = 29 , LightBaryon = 30 ,   StrangeBaryon = 31 , CharmedBaryon = 32 , BottomBaryon = 33 , PionDecay = 34 ,   KaonDecay = 35 , BremPhot = 36 , PromptPhot = 37 , UndrPhot = 38 ,   ISRPhot = 39 , FSRPhot = 40 , NucReact = 41 , PiZero = 42 ,   DiBoson = 43 , ZorHeavyBoson = 44 , MultiBoson = 47 , QCD = 45 }
enum	ParticleOutCome : unsigned int {   NonDefinedOutCome = 0 , UnknownOutCome = 1 , UnConverted = 2 , Converted = 3 ,   NonInteract = 4 , NuclInteraction = 5 , ElectrMagInter = 6 , DecaytoElectron = 7 ,   DecaytoMuon = 8 , OneProng = 9 , ThreeProng = 10 , FiveProng = 11 ,   NonDefinedOutCome = 0 , UnknownOutCome = 1 , UnConverted = 2 , Converted = 3 ,   NonInteract = 4 , NuclInteraction = 5 , ElectrMagInter = 6 , DecaytoElectron = 7 ,   DecaytoMuon = 8 , OneProng = 9 , ThreeProng = 10 , FiveProng = 11 }
enum	MCTC_bits : unsigned int {   HadTau =0 , Tau , hadron , frombsm ,   uncat , isbsm , isgeant , stable ,   totalBits }
enum	ParticleType : unsigned int {   Unknown = 0 , UnknownElectron = 1 , IsoElectron = 2 , NonIsoElectron = 3 ,   BkgElectron = 4 , UnknownMuon = 5 , IsoMuon = 6 , NonIsoMuon = 7 ,   BkgMuon = 8 , UnknownTau = 9 , IsoTau = 10 , NonIsoTau = 11 ,   BkgTau = 12 , UnknownPhoton = 13 , IsoPhoton = 14 , NonIsoPhoton = 15 ,   BkgPhoton = 16 , Hadron = 17 , Neutrino = 18 , NuclFrag = 19 ,   NonPrimary = 20 , GenParticle = 21 , SUSYParticle = 22 , OtherBSMParticle = 39 ,   BBbarMesonPart = 23 , BottomMesonPart = 24 , CCbarMesonPart = 25 , CharmedMesonPart = 26 ,   BottomBaryonPart = 27 , CharmedBaryonPart = 28 , StrangeBaryonPart = 29 , LightBaryonPart = 30 ,   StrangeMesonPart = 31 , LightMesonPart = 32 , BJet = 33 , CJet = 34 ,   LJet = 35 , GJet = 36 , TauJet = 37 , UnknownJet = 38 ,   Unknown = 0 , UnknownElectron = 1 , IsoElectron = 2 , NonIsoElectron = 3 ,   BkgElectron = 4 , UnknownMuon = 5 , IsoMuon = 6 , NonIsoMuon = 7 ,   BkgMuon = 8 , UnknownTau = 9 , IsoTau = 10 , NonIsoTau = 11 ,   BkgTau = 12 , UnknownPhoton = 13 , IsoPhoton = 14 , NonIsoPhoton = 15 ,   BkgPhoton = 16 , Hadron = 17 , Neutrino = 18 , NuclFrag = 19 ,   NonPrimary = 20 , GenParticle = 21 , SUSYParticle = 22 , OtherBSMParticle = 39 ,   BBbarMesonPart = 23 , BottomMesonPart = 24 , CCbarMesonPart = 25 , CharmedMesonPart = 26 ,   BottomBaryonPart = 27 , CharmedBaryonPart = 28 , StrangeBaryonPart = 29 , LightBaryonPart = 30 ,   StrangeMesonPart = 31 , LightMesonPart = 32 , BJet = 33 , CJet = 34 ,   LJet = 35 , GJet = 36 , TauJet = 37 , UnknownJet = 38 }
enum	ParticleOrigin : unsigned int {   NonDefined = 0 , SingleElec = 1 , SingleMuon = 2 , SinglePhot = 3 ,   SingleTau = 4 , PhotonConv = 5 , DalitzDec = 6 , ElMagProc = 7 ,   Mu = 8 , TauLep = 9 , top = 10 , QuarkWeakDec = 11 ,   WBoson = 12 , ZBoson = 13 , Higgs = 14 , HiggsMSSM = 15 ,   HeavyBoson = 16 , WBosonLRSM = 17 , NuREle = 18 , NuRMu = 19 ,   NuRTau = 20 , LQ = 21 , SUSY = 22 , OtherBSM = 46 ,   LightMeson = 23 , StrangeMeson = 24 , CharmedMeson = 25 , BottomMeson = 26 ,   CCbarMeson = 27 , JPsi = 28 , BBbarMeson = 29 , LightBaryon = 30 ,   StrangeBaryon = 31 , CharmedBaryon = 32 , BottomBaryon = 33 , PionDecay = 34 ,   KaonDecay = 35 , BremPhot = 36 , PromptPhot = 37 , UndrPhot = 38 ,   ISRPhot = 39 , FSRPhot = 40 , NucReact = 41 , PiZero = 42 ,   DiBoson = 43 , ZorHeavyBoson = 44 , MultiBoson = 47 , QCD = 45 ,   NonDefined = 0 , SingleElec = 1 , SingleMuon = 2 , SinglePhot = 3 ,   SingleTau = 4 , PhotonConv = 5 , DalitzDec = 6 , ElMagProc = 7 ,   Mu = 8 , TauLep = 9 , top = 10 , QuarkWeakDec = 11 ,   WBoson = 12 , ZBoson = 13 , Higgs = 14 , HiggsMSSM = 15 ,   HeavyBoson = 16 , WBosonLRSM = 17 , NuREle = 18 , NuRMu = 19 ,   NuRTau = 20 , LQ = 21 , SUSY = 22 , OtherBSM = 46 ,   LightMeson = 23 , StrangeMeson = 24 , CharmedMeson = 25 , BottomMeson = 26 ,   CCbarMeson = 27 , JPsi = 28 , BBbarMeson = 29 , LightBaryon = 30 ,   StrangeBaryon = 31 , CharmedBaryon = 32 , BottomBaryon = 33 , PionDecay = 34 ,   KaonDecay = 35 , BremPhot = 36 , PromptPhot = 37 , UndrPhot = 38 ,   ISRPhot = 39 , FSRPhot = 40 , NucReact = 41 , PiZero = 42 ,   DiBoson = 43 , ZorHeavyBoson = 44 , MultiBoson = 47 , QCD = 45 }
enum	ParticleOutCome : unsigned int {   NonDefinedOutCome = 0 , UnknownOutCome = 1 , UnConverted = 2 , Converted = 3 ,   NonInteract = 4 , NuclInteraction = 5 , ElectrMagInter = 6 , DecaytoElectron = 7 ,   DecaytoMuon = 8 , OneProng = 9 , ThreeProng = 10 , FiveProng = 11 ,   NonDefinedOutCome = 0 , UnknownOutCome = 1 , UnConverted = 2 , Converted = 3 ,   NonInteract = 4 , NuclInteraction = 5 , ElectrMagInter = 6 , DecaytoElectron = 7 ,   DecaytoMuon = 8 , OneProng = 9 , ThreeProng = 10 , FiveProng = 11 }

Functions
ParticleOrigin	convHadronTypeToOrig (ParticleType pType, int motherPDG)
ParticleOrigin	defHadronType (int pdg)
ParticleType	defTypeOfHadron (int pdg)
ParticleType	defTypeOfElectron (ParticleOrigin EleOrig, bool isPrompt)
ParticleType	defTypeOfMuon (ParticleOrigin MuOrig, bool isPrompt)
ParticleType	defTypeOfTau (ParticleOrigin TauOrig)
ParticleType	defTypeOfPhoton (ParticleOrigin PhotOrig)
template<class T>
ParticleOrigin	defJetOrig (const T &allJetMothers)
template<class T>
std::tuple< unsigned int, T >	defOrigOfParticle (T thePart)
int	isPrompt (const unsigned int classify, bool allow_prompt_tau_decays=true)
template<class T>
ParticleOutCome	defOutComeOfElectron (T thePart)
template<class T>
ParticleOutCome	defOutComeOfMuon (T thePart)
template<class T>
ParticleOutCome	defOutComeOfTau (T thePart)
template<class T>
ParticleOutCome	defOutComeOfPhoton (T thePart)

Enumeration Type Documentation

MCTC_bits

enum MCTruthPartClassifier::MCTC_bits : unsigned int

Enumerator
HadTau
Tau
hadron
frombsm
uncat
isbsm
isgeant
stable
totalBits

Definition at line 148 of file TruthClassifiers.h.

: unsigned int { HadTau=0, Tau, hadron, frombsm, uncat, isbsm, isgeant, stable, totalBits };

ParticleOrigin [1/2]

enum MCTruthPartClassifier::ParticleOrigin : unsigned int

Enumerator
NonDefined
SingleElec
SingleMuon
SinglePhot
SingleTau
PhotonConv
DalitzDec
ElMagProc
Mu
TauLep
top
QuarkWeakDec
WBoson
ZBoson
Higgs
HiggsMSSM
HeavyBoson
WBosonLRSM
NuREle
NuRMu
NuRTau
LQ
SUSY
OtherBSM
LightMeson
StrangeMeson
CharmedMeson
BottomMeson
CCbarMeson
JPsi
BBbarMeson
LightBaryon
StrangeBaryon
CharmedBaryon
BottomBaryon
PionDecay
KaonDecay
BremPhot
PromptPhot
UndrPhot
ISRPhot
FSRPhot
NucReact
PiZero
DiBoson
ZorHeavyBoson
MultiBoson
QCD
NonDefined
SingleElec
SingleMuon
SinglePhot
SingleTau
PhotonConv
DalitzDec
ElMagProc
Mu
TauLep
top
QuarkWeakDec
WBoson
ZBoson
Higgs
HiggsMSSM
HeavyBoson
WBosonLRSM
NuREle
NuRMu
NuRTau
LQ
SUSY
OtherBSM
LightMeson
StrangeMeson
CharmedMeson
BottomMeson
CCbarMeson
JPsi
BBbarMeson
LightBaryon
StrangeBaryon
CharmedBaryon
BottomBaryon
PionDecay
KaonDecay
BremPhot
PromptPhot
UndrPhot
ISRPhot
FSRPhot
NucReact
PiZero
DiBoson
ZorHeavyBoson
MultiBoson
QCD

Definition at line 52 of file TruthClassifiers.h.

                                                                             {
 
  if (EleOrig == NonDefined)
    return UnknownElectron;
 
  if (EleOrig == WBoson || EleOrig == ZBoson || EleOrig == top || EleOrig == SingleElec || EleOrig == Higgs ||
      EleOrig == HiggsMSSM || EleOrig == HeavyBoson || EleOrig == WBosonLRSM || EleOrig == NuREle || EleOrig == NuRMu ||
      EleOrig == NuRTau || EleOrig == LQ || EleOrig == SUSY || EleOrig == DiBoson || EleOrig == ZorHeavyBoson ||
      EleOrig == OtherBSM || EleOrig == MultiBoson || isPrompt) {
    return IsoElectron;
  }
  if (EleOrig == JPsi || EleOrig == BottomMeson || EleOrig == CharmedMeson || EleOrig == BottomBaryon ||
      EleOrig == CharmedBaryon || EleOrig == TauLep || EleOrig == Mu || EleOrig == QuarkWeakDec) {
    return NonIsoElectron;
  }
  return BkgElectron;
}
 
 
inline ParticleType defTypeOfMuon(ParticleOrigin MuOrig, bool isPrompt) {
 
  if (MuOrig == NonDefined) return UnknownMuon;
 
  if (MuOrig == WBoson || MuOrig == ZBoson || MuOrig == top || MuOrig == SingleMuon || MuOrig == Higgs ||
      MuOrig == HiggsMSSM || MuOrig == HeavyBoson || MuOrig == WBosonLRSM || MuOrig == NuREle || MuOrig == NuRMu ||
      MuOrig == NuRTau || MuOrig == LQ || MuOrig == SUSY || MuOrig == DiBoson || MuOrig == ZorHeavyBoson ||
      MuOrig == OtherBSM || MuOrig == MultiBoson || isPrompt) {
    return IsoMuon;
  }
  if (MuOrig == JPsi || MuOrig == BottomMeson || MuOrig == CharmedMeson || MuOrig == BottomBaryon ||
      MuOrig == CharmedBaryon || MuOrig == TauLep || MuOrig == QuarkWeakDec) {
    return NonIsoMuon;
  }
  //  if (MuOrig == Pion  || MuOrig == Kaon ) return  DecayMuon;
  return BkgMuon;
}
 
inline ParticleType defTypeOfTau(ParticleOrigin TauOrig) {
  if (TauOrig == NonDefined) return UnknownTau;
 
  if (TauOrig == WBoson || TauOrig == ZBoson || TauOrig == top || TauOrig == SingleMuon || TauOrig == Higgs ||
      TauOrig == HiggsMSSM || TauOrig == HeavyBoson || TauOrig == WBosonLRSM || TauOrig == NuREle || TauOrig == NuRMu ||
      TauOrig == NuRTau || TauOrig == SUSY || TauOrig == DiBoson || TauOrig == ZorHeavyBoson || TauOrig == OtherBSM ||
      TauOrig == MultiBoson)

ParticleOrigin [2/2]

enum MCTruthPartClassifier::ParticleOrigin : unsigned int

Enumerator
NonDefined
SingleElec
SingleMuon
SinglePhot
SingleTau
PhotonConv
DalitzDec
ElMagProc
Mu
TauLep
top
QuarkWeakDec
WBoson
ZBoson
Higgs
HiggsMSSM
HeavyBoson
WBosonLRSM
NuREle
NuRMu
NuRTau
LQ
SUSY
OtherBSM
LightMeson
StrangeMeson
CharmedMeson
BottomMeson
CCbarMeson
JPsi
BBbarMeson
LightBaryon
StrangeBaryon
CharmedBaryon
BottomBaryon
PionDecay
KaonDecay
BremPhot
PromptPhot
UndrPhot
ISRPhot
FSRPhot
NucReact
PiZero
DiBoson
ZorHeavyBoson
MultiBoson
QCD
NonDefined
SingleElec
SingleMuon
SinglePhot
SingleTau
PhotonConv
DalitzDec
ElMagProc
Mu
TauLep
top
QuarkWeakDec
WBoson
ZBoson
Higgs
HiggsMSSM
HeavyBoson
WBosonLRSM
NuREle
NuRMu
NuRTau
LQ
SUSY
OtherBSM
LightMeson
StrangeMeson
CharmedMeson
BottomMeson
CCbarMeson
JPsi
BBbarMeson
LightBaryon
StrangeBaryon
CharmedBaryon
BottomBaryon
PionDecay
KaonDecay
BremPhot
PromptPhot
UndrPhot
ISRPhot
FSRPhot
NucReact
PiZero
DiBoson
ZorHeavyBoson
MultiBoson
QCD

Definition at line 52 of file MCTruthClassifierDefs.h.

ParticleOutCome [1/2]

enum MCTruthPartClassifier::ParticleOutCome : unsigned int

Enumerator
NonDefinedOutCome
UnknownOutCome
UnConverted
Converted
NonInteract
NuclInteraction
ElectrMagInter
DecaytoElectron
DecaytoMuon
OneProng
ThreeProng
FiveProng
NonDefinedOutCome
UnknownOutCome
UnConverted
Converted
NonInteract
NuclInteraction
ElectrMagInter
DecaytoElectron
DecaytoMuon
OneProng
ThreeProng
FiveProng

Definition at line 106 of file TruthClassifiers.h.

{
  if (PhotOrig == NonDefined) return UnknownPhoton;
 
  if (PhotOrig == WBoson || PhotOrig == ZBoson || PhotOrig == SinglePhot || PhotOrig == Higgs ||

ParticleOutCome [2/2]

enum MCTruthPartClassifier::ParticleOutCome : unsigned int

Enumerator
NonDefinedOutCome
UnknownOutCome
UnConverted
Converted
NonInteract
NuclInteraction
ElectrMagInter
DecaytoElectron
DecaytoMuon
OneProng
ThreeProng
FiveProng
NonDefinedOutCome
UnknownOutCome
UnConverted
Converted
NonInteract
NuclInteraction
ElectrMagInter
DecaytoElectron
DecaytoMuon
OneProng
ThreeProng
FiveProng

Definition at line 106 of file MCTruthClassifierDefs.h.

                                {
#include "TruthUtils/TruthClasses.h"
}
 
#endif // MCTRUTHCLASSIFIER_MCTRUTHCLASSIFIERDEFS_H

ParticleType [1/2]

enum MCTruthPartClassifier::ParticleType : unsigned int

Enumerator
Unknown
UnknownElectron
IsoElectron
NonIsoElectron
BkgElectron
UnknownMuon
IsoMuon
NonIsoMuon
BkgMuon
UnknownTau
IsoTau
NonIsoTau
BkgTau
UnknownPhoton
IsoPhoton
NonIsoPhoton
BkgPhoton
Hadron
Neutrino
NuclFrag
NonPrimary
GenParticle
SUSYParticle
OtherBSMParticle
BBbarMesonPart
BottomMesonPart
CCbarMesonPart
CharmedMesonPart
BottomBaryonPart
CharmedBaryonPart
StrangeBaryonPart
LightBaryonPart
StrangeMesonPart
LightMesonPart
BJet
CJet
LJet
GJet
TauJet
UnknownJet
Unknown
UnknownElectron
IsoElectron
NonIsoElectron
BkgElectron
UnknownMuon
IsoMuon
NonIsoMuon
BkgMuon
UnknownTau
IsoTau
NonIsoTau
BkgTau
UnknownPhoton
IsoPhoton
NonIsoPhoton
BkgPhoton
Hadron
Neutrino
NuclFrag
NonPrimary
GenParticle
SUSYParticle
OtherBSMParticle
BBbarMesonPart
BottomMesonPart
CCbarMesonPart
CharmedMesonPart
BottomBaryonPart
CharmedBaryonPart
StrangeBaryonPart
LightBaryonPart
StrangeMesonPart
LightMesonPart
BJet
CJet
LJet
GJet
TauJet
UnknownJet

Definition at line 9 of file TruthClassifiers.h.

                                {
#include "TruthUtils/TruthClasses.h"
 
inline ParticleOrigin convHadronTypeToOrig(ParticleType pType, int motherPDG)
{
  if (pType == CCbarMesonPart && abs(motherPDG) == MC::JPSI) return JPsi;
  if (pType == BBbarMesonPart) return BBbarMeson;
  if (pType == BottomMesonPart) return BottomMeson;
  if (pType == BottomBaryonPart) return BottomBaryon;
  if (pType == CCbarMesonPart) return CCbarMeson;
  if (pType == CharmedMesonPart) return CharmedMeson;
  if (pType == CharmedBaryonPart) return CharmedBaryon;
  if (pType == StrangeBaryonPart) return StrangeBaryon;
  if (pType == StrangeMesonPart) return StrangeMeson;
  if (pType == LightBaryonPart) return LightBaryon;
  if (pType == LightMesonPart) return LightMeson;
  return NonDefined;
}
 
inline ParticleOrigin defHadronType(int pdg) {
  if (abs(pdg) == MC::JPSI) return JPsi;
  if (MC::isBBbarMeson(pdg)) return BBbarMeson;
  if (MC::isCCbarMeson(pdg)) return CCbarMeson;
  if (MC::isBottomMeson(pdg)) return BottomMeson;
  if (MC::isCharmMeson(pdg)) return CharmedMeson;
  if (MC::isBottomBaryon(pdg)) return BottomBaryon;
  if (MC::isCharmBaryon(pdg)) return CharmedBaryon;
  if (MC::isStrangeBaryon(pdg)) return StrangeBaryon;
  if (MC::isLightBaryon(pdg)) return LightBaryon;
  if (MC::isStrangeMeson(pdg)) return StrangeMeson;
  if (MC::isLightMeson(pdg)) return LightMeson;
  return NonDefined;
}
 
inline ParticleType defTypeOfHadron(int pdg) {
  if (MC::isBBbarMeson(pdg)) return BBbarMesonPart;
  if (MC::isCCbarMeson(pdg)) return CCbarMesonPart;
  if (MC::isBottomMeson(pdg)) return BottomMesonPart;
  if (MC::isCharmMeson(pdg)) return CharmedMesonPart;

ParticleType [2/2]

enum MCTruthPartClassifier::ParticleType : unsigned int

Enumerator
Unknown
UnknownElectron
IsoElectron
NonIsoElectron
BkgElectron
UnknownMuon
IsoMuon
NonIsoMuon
BkgMuon
UnknownTau
IsoTau
NonIsoTau
BkgTau
UnknownPhoton
IsoPhoton
NonIsoPhoton
BkgPhoton
Hadron
Neutrino
NuclFrag
NonPrimary
GenParticle
SUSYParticle
OtherBSMParticle
BBbarMesonPart
BottomMesonPart
CCbarMesonPart
CharmedMesonPart
BottomBaryonPart
CharmedBaryonPart
StrangeBaryonPart
LightBaryonPart
StrangeMesonPart
LightMesonPart
BJet
CJet
LJet
GJet
TauJet
UnknownJet
Unknown
UnknownElectron
IsoElectron
NonIsoElectron
BkgElectron
UnknownMuon
IsoMuon
NonIsoMuon
BkgMuon
UnknownTau
IsoTau
NonIsoTau
BkgTau
UnknownPhoton
IsoPhoton
NonIsoPhoton
BkgPhoton
Hadron
Neutrino
NuclFrag
NonPrimary
GenParticle
SUSYParticle
OtherBSMParticle
BBbarMesonPart
BottomMesonPart
CCbarMesonPart
CharmedMesonPart
BottomBaryonPart
CharmedBaryonPart
StrangeBaryonPart
LightBaryonPart
StrangeMesonPart
LightMesonPart
BJet
CJet
LJet
GJet
TauJet
UnknownJet

Definition at line 9 of file MCTruthClassifierDefs.h.

    :     MCTruthClassifierDefs.h
PACKAGE:  atlasoff/PhysicsAnalysis/MCTruthClassifier
 
AUTHORS:  O. Fedin
CREATED:  Oct 2007
 
For more information, please see the twiki page:
https://twiki.cern.ch/twiki/bin/view/AtlasProtected/MCTruthClassifier
********************************************************************
Updated Feb 2024 by Andrii Verbytskyi <andrii.verbytskyi@mpp.mpg.de>
********************************************************************/
#include <vector>
#include <string>
namespace MCTruthPartClassifier {
#include "TruthUtils/TruthClasses.h"
}
 
#endif // MCTRUTHCLASSIFIER_MCTRUTHCLASSIFIERDEFS_H

Function Documentation

convHadronTypeToOrig()

ParticleOrigin MCTruthPartClassifier::convHadronTypeToOrig ( ParticleType pType, int motherPDG )

inline

Definition at line 15 of file TruthClassifiers.h.

{
  if (pType == CCbarMesonPart && abs(motherPDG) == MC::JPSI) return JPsi;
  if (pType == BBbarMesonPart) return BBbarMeson;
  if (pType == BottomMesonPart) return BottomMeson;
  if (pType == BottomBaryonPart) return BottomBaryon;
  if (pType == CCbarMesonPart) return CCbarMeson;
  if (pType == CharmedMesonPart) return CharmedMeson;
  if (pType == CharmedBaryonPart) return CharmedBaryon;
  if (pType == StrangeBaryonPart) return StrangeBaryon;
  if (pType == StrangeMesonPart) return StrangeMeson;
  if (pType == LightBaryonPart) return LightBaryon;
  if (pType == LightMesonPart) return LightMeson;
  return NonDefined;
}

defHadronType()

ParticleOrigin MCTruthPartClassifier::defHadronType ( int pdg )

inline

Definition at line 31 of file TruthClassifiers.h.

                                             {
  if (abs(pdg) == MC::JPSI) return JPsi;
  if (MC::isBBbarMeson(pdg)) return BBbarMeson;
  if (MC::isCCbarMeson(pdg)) return CCbarMeson;
  if (MC::isBottomMeson(pdg)) return BottomMeson;
  if (MC::isCharmMeson(pdg)) return CharmedMeson;
  if (MC::isBottomBaryon(pdg)) return BottomBaryon;
  if (MC::isCharmBaryon(pdg)) return CharmedBaryon;
  if (MC::isStrangeBaryon(pdg)) return StrangeBaryon;
  if (MC::isLightBaryon(pdg)) return LightBaryon;
  if (MC::isStrangeMeson(pdg)) return StrangeMeson;
  if (MC::isLightMeson(pdg)) return LightMeson;
  return NonDefined;
}

defJetOrig()

template<class T>

ParticleOrigin MCTruthPartClassifier::defJetOrig

(

const T &

allJetMothers

)

Definition at line 131 of file TruthClassifiers.h.

                                                                     {
  ParticleOrigin partOrig = NonDefined;
  for (const auto& it: allJetMothers) {
    int pdg = abs(it->pdg_id());
    if (MC::isTop(pdg)) partOrig = top;
    if (MC::isZ(pdg))   partOrig = ZBoson;
    if (MC::isW(pdg) && !(partOrig == top)) partOrig = WBoson;
    if ((MC::isQuark(pdg) || MC::isGluon(pdg)) && partOrig != top && partOrig != ZBoson && partOrig != WBoson) partOrig = QCD;
    if (MC::isHiggs(pdg)) return Higgs;
    if (pdg == 35 || pdg == 36 || pdg == 37) return HiggsMSSM;
    if (pdg == 32 || pdg == 33 || pdg == 34) return HeavyBoson;
    if (pdg == 42) return LQ;
    if (MC::isSUSY(pdg)) return SUSY;
    if (MC::isBSM(pdg)) return OtherBSM;
  }
  return partOrig;
}

defOrigOfParticle()

template<class T>

std::tuple< unsigned int, T > MCTruthPartClassifier::defOrigOfParticle

(

T

thePart

)

Definition at line 150 of file TruthClassifiers.h.

                                                                          {
 
  T parent_hadron_ptr = nullptr;
 
  bool uncat = 0, fromHad = 0, fromTau = 0;
  bool isPhysical = MC::isPhysical(thePart);
  bool isGeant = HepMC::is_simulation_particle(thePart);
  bool isBSM = MC::isBSM(thePart);
  bool fromBSM = isBSM; // just to initialise
 
  auto prodVtx = thePart->production_vertex();
  if (isPhysical && prodVtx && !isGeant) {
    fromHad = MC::isFromHadron(thePart, parent_hadron_ptr, fromTau, fromBSM);
  }
  else  uncat = 1;
 
  std::bitset<MCTC_bits::totalBits> classifier;
  classifier[MCTC_bits::stable]  = isPhysical;
  classifier[MCTC_bits::isgeant] = isGeant;
  classifier[MCTC_bits::isbsm]   = isBSM;
  classifier[MCTC_bits::uncat]   = uncat;
  classifier[MCTC_bits::frombsm] = fromBSM;
  classifier[MCTC_bits::hadron]  = fromHad;
  classifier[MCTC_bits::Tau]     = fromTau;
  classifier[MCTC_bits::HadTau]  = fromHad && fromTau;
  unsigned int outputvalue = static_cast<unsigned int>(classifier.to_ulong());
 
  return std::make_tuple(outputvalue,parent_hadron_ptr);
}

defOutComeOfElectron()

template<class T>

ParticleOutCome MCTruthPartClassifier::defOutComeOfElectron

(

T

thePart

)

Definition at line 188 of file TruthClassifiers.h.

                                                                   {
  ParticleOutCome PartOutCome = UnknownOutCome;
  auto EndVert = MC::findSimulatedEndVertex(thePart);
  if (EndVert == nullptr) return NonInteract;
 
  int ElecOutNumOfNucFr(0);
  int ElecOutNumOfElec(0);
  int NumOfHadr(0);
  auto outgoing = EndVert->particles_out();
  int NumOfElecDaug = outgoing.size();
  for (const auto& p: outgoing) {
    if (!p) continue;
    int EndDaugType = p->pdg_id();
    if (MC::isElectron(EndDaugType)) ElecOutNumOfElec++;
    if (MC::isHadron(p) && !MC::isBeam(p)) NumOfHadr++;
    if (EndDaugType > 1000000000 || EndDaugType == 0 || abs(EndDaugType) == 2212 || abs(EndDaugType) == 2112) ElecOutNumOfNucFr++;
  }
 
  if (ElecOutNumOfNucFr != 0 || NumOfHadr != 0)  PartOutCome = NuclInteraction;
  if (ElecOutNumOfElec == 1 && NumOfElecDaug == 1) PartOutCome = ElectrMagInter;
 
  return PartOutCome;
}

defOutComeOfMuon()

template<class T>

ParticleOutCome MCTruthPartClassifier::defOutComeOfMuon

(

T

thePart

)

Definition at line 212 of file TruthClassifiers.h.

                                                               {
  ParticleOutCome PartOutCome = UnknownOutCome;
  auto EndVert = MC::findSimulatedEndVertex(thePart);
  if (EndVert == nullptr) return NonInteract;
  int MuOutNumOfNucFr(0);
  int NumOfHadr(0);
  int NumOfEleNeutr(0);
  int NumOfMuonNeutr(0);
  int NumOfElec(0);
  auto outgoing = EndVert->particles_out();
  int NumOfMuDaug = outgoing.size();
  for (const auto& p: outgoing) {
    if (!p) continue;
    int EndDaugType = p->pdg_id();
    if (MC::isElectron(EndDaugType)) NumOfElec++;
    if (abs(EndDaugType) == 12) NumOfEleNeutr++;
    if (abs(EndDaugType) == 14) NumOfMuonNeutr++;
    if (MC::isHadron(p) && !MC::isBeam(p)) NumOfHadr++;
    if (EndDaugType > 1000000000 || EndDaugType == 0 || abs(EndDaugType) == 2212 || abs(EndDaugType) == 2112) MuOutNumOfNucFr++;
  }
 
  if (MuOutNumOfNucFr != 0 || NumOfHadr != 0) PartOutCome = NuclInteraction;
  if (NumOfMuDaug == 3 && NumOfElec == 1 && NumOfEleNeutr == 1 && NumOfMuonNeutr == 1) PartOutCome = DecaytoElectron;
 
  return PartOutCome;
}

defOutComeOfPhoton()

template<class T>

ParticleOutCome MCTruthPartClassifier::defOutComeOfPhoton

(

T

thePart

)

Definition at line 266 of file TruthClassifiers.h.

                                                                 {
  ParticleOutCome PartOutCome = UnknownOutCome;
  auto EndVert = MC::findSimulatedEndVertex(thePart);
  if (EndVert == nullptr) return UnConverted;
 
  int PhtOutNumOfNucFr(0);
  int PhtOutNumOfEl(0);
  int PhtOutNumOfPos(0);
  int PhtOutNumOfHadr(0);
 
  auto outgoing = EndVert->particles_out();
  int NumOfPhtDaug = outgoing.size();
  for (const auto& p: outgoing) {
    if (!p) continue;
    int EndDaugType = p->pdg_id();
    if (EndDaugType > 1000000000 || EndDaugType == 0 || abs(EndDaugType) == 2212 || abs(EndDaugType) == 2112) PhtOutNumOfNucFr++;
    if (EndDaugType == 11)  PhtOutNumOfEl++;
    if (EndDaugType == -11) PhtOutNumOfPos++;
    if (MC::isHadron(p) && !MC::isBeam(p)) PhtOutNumOfHadr++;
  }
 
  if (PhtOutNumOfEl == 1 && PhtOutNumOfPos == 1 && NumOfPhtDaug == 2) PartOutCome = Converted;
  if ((NumOfPhtDaug > 1 && PhtOutNumOfNucFr != 0) || PhtOutNumOfHadr > 0) PartOutCome = NuclInteraction;
 
  return PartOutCome;
}

defOutComeOfTau()

template<class T>

ParticleOutCome MCTruthPartClassifier::defOutComeOfTau

(

T

thePart

)

Definition at line 238 of file TruthClassifiers.h.

                                                              {
  ParticleOutCome PartOutCome = UnknownOutCome;
  auto EndVert = MC::findSimulatedEndVertex(thePart);
  if (EndVert == nullptr) return NonInteract;
  int NumOfTauDaug = EndVert->nOutgoingParticles();
  auto tauFinalStatePart = MC::findFinalStateParticles(EndVert);
  auto PD = DecayProducts(tauFinalStatePart);
  int NumOfElec = PD.apd(11);
  int NumOfMuon = PD.apd(13);
  int NumOfElecNeut = PD.apd(12);
  int NumOfMuonNeut = PD.apd(14);
  int NumOfPhot = PD.apd(22);
  int NumOfPi = PD.apd(211);
  int NumOfKaon = PD.apd(321);
  int NumOfNucFr = PD.apd(0) + PD.apd(1000000000, std::numeric_limits<int>::max());
 
  if (NumOfNucFr != 0) PartOutCome = NuclInteraction;
  if ((NumOfTauDaug == 3 && NumOfElec == 1 && NumOfElecNeut == 1) || (NumOfTauDaug == (3 + NumOfPhot) && NumOfElecNeut == 1)) PartOutCome = DecaytoElectron;
  if ((NumOfTauDaug == 3 && NumOfMuon == 1 && NumOfMuonNeut == 1) || (NumOfTauDaug == (3 + NumOfPhot) && NumOfMuonNeut == 1)) PartOutCome = DecaytoMuon;
 
  if (NumOfPi == 1 || NumOfKaon == 1) PartOutCome = OneProng;
  if (NumOfPi + NumOfKaon == 3) PartOutCome = ThreeProng;
  if (NumOfPi + NumOfKaon == 5) PartOutCome = FiveProng;
 
 
  return PartOutCome;
}

defTypeOfElectron()

ParticleType MCTruthPartClassifier::defTypeOfElectron ( ParticleOrigin EleOrig, bool isPrompt )

inline

Definition at line 61 of file TruthClassifiers.h.

                                                                             {
 
  if (EleOrig == NonDefined)
    return UnknownElectron;
 
  if (EleOrig == WBoson || EleOrig == ZBoson || EleOrig == top || EleOrig == SingleElec || EleOrig == Higgs ||
      EleOrig == HiggsMSSM || EleOrig == HeavyBoson || EleOrig == WBosonLRSM || EleOrig == NuREle || EleOrig == NuRMu ||
      EleOrig == NuRTau || EleOrig == LQ || EleOrig == SUSY || EleOrig == DiBoson || EleOrig == ZorHeavyBoson ||
      EleOrig == OtherBSM || EleOrig == MultiBoson || isPrompt) {
    return IsoElectron;
  }
  if (EleOrig == JPsi || EleOrig == BottomMeson || EleOrig == CharmedMeson || EleOrig == BottomBaryon ||
      EleOrig == CharmedBaryon || EleOrig == TauLep || EleOrig == Mu || EleOrig == QuarkWeakDec) {
    return NonIsoElectron;
  }
  return BkgElectron;
}

defTypeOfHadron()

ParticleType MCTruthPartClassifier::defTypeOfHadron ( int pdg )

inline

Definition at line 46 of file TruthClassifiers.h.

                                             {
  if (MC::isBBbarMeson(pdg)) return BBbarMesonPart;
  if (MC::isCCbarMeson(pdg)) return CCbarMesonPart;
  if (MC::isBottomMeson(pdg)) return BottomMesonPart;
  if (MC::isCharmMeson(pdg)) return CharmedMesonPart;
  if (MC::isBottomBaryon(pdg)) return BottomBaryonPart;
  if (MC::isCharmBaryon(pdg)) return CharmedBaryonPart;
  if (MC::isStrangeBaryon(pdg)) return StrangeBaryonPart;
  if (MC::isLightBaryon(pdg)) return LightBaryonPart;
  if (MC::isStrangeMeson(pdg)) return StrangeMesonPart;
  if (MC::isLightMeson(pdg)) return LightMesonPart;
  return Unknown;
}

defTypeOfMuon()

ParticleType MCTruthPartClassifier::defTypeOfMuon ( ParticleOrigin MuOrig, bool isPrompt )

inline

Definition at line 80 of file TruthClassifiers.h.

                                                                        {
 
  if (MuOrig == NonDefined) return UnknownMuon;
 
  if (MuOrig == WBoson || MuOrig == ZBoson || MuOrig == top || MuOrig == SingleMuon || MuOrig == Higgs ||
      MuOrig == HiggsMSSM || MuOrig == HeavyBoson || MuOrig == WBosonLRSM || MuOrig == NuREle || MuOrig == NuRMu ||
      MuOrig == NuRTau || MuOrig == LQ || MuOrig == SUSY || MuOrig == DiBoson || MuOrig == ZorHeavyBoson ||
      MuOrig == OtherBSM || MuOrig == MultiBoson || isPrompt) {
    return IsoMuon;
  }
  if (MuOrig == JPsi || MuOrig == BottomMeson || MuOrig == CharmedMeson || MuOrig == BottomBaryon ||
      MuOrig == CharmedBaryon || MuOrig == TauLep || MuOrig == QuarkWeakDec) {
    return NonIsoMuon;
  }
  //  if (MuOrig == Pion  || MuOrig == Kaon ) return  DecayMuon;
  return BkgMuon;
}

defTypeOfPhoton()

ParticleType MCTruthPartClassifier::defTypeOfPhoton ( ParticleOrigin PhotOrig )

inline

Definition at line 115 of file TruthClassifiers.h.

{
  if (PhotOrig == NonDefined) return UnknownPhoton;
 
  if (PhotOrig == WBoson || PhotOrig == ZBoson || PhotOrig == SinglePhot || PhotOrig == Higgs ||
      PhotOrig == HiggsMSSM || PhotOrig == HeavyBoson || PhotOrig == PromptPhot || PhotOrig == SUSY ||
      PhotOrig == OtherBSM)
    return IsoPhoton;
 
  if (PhotOrig == ISRPhot || PhotOrig == FSRPhot || PhotOrig == TauLep || PhotOrig == Mu || PhotOrig == NuREle ||
      PhotOrig == NuRMu || PhotOrig == NuRTau)
    return NonIsoPhoton;
 
  return BkgPhoton;
}

defTypeOfTau()

ParticleType MCTruthPartClassifier::defTypeOfTau ( ParticleOrigin TauOrig )

inline

Definition at line 98 of file TruthClassifiers.h.

                                                         {
  if (TauOrig == NonDefined) return UnknownTau;
 
  if (TauOrig == WBoson || TauOrig == ZBoson || TauOrig == top || TauOrig == SingleMuon || TauOrig == Higgs ||
      TauOrig == HiggsMSSM || TauOrig == HeavyBoson || TauOrig == WBosonLRSM || TauOrig == NuREle || TauOrig == NuRMu ||
      TauOrig == NuRTau || TauOrig == SUSY || TauOrig == DiBoson || TauOrig == ZorHeavyBoson || TauOrig == OtherBSM ||
      TauOrig == MultiBoson)
    return IsoTau;
 
  if (TauOrig == JPsi || TauOrig == BottomMeson || TauOrig == CharmedMeson || TauOrig == BottomBaryon ||
      TauOrig == CharmedBaryon || TauOrig == QuarkWeakDec)
    return NonIsoTau;
 
  return BkgTau;
}

isPrompt()

int MCTruthPartClassifier::isPrompt ( const unsigned int classify, bool allow_prompt_tau_decays = true )

inline

Definition at line 180 of file TruthClassifiers.h.

                                                                                    {
    std::bitset<MCTC_bits::totalBits> res(classify);
    if (res.test(MCTC_bits::uncat)) return -1;
    bool fromPromptTau = res.test(MCTC_bits::Tau) && !res.test(MCTC_bits::HadTau);
    if (fromPromptTau) return int(allow_prompt_tau_decays);
    return !res.test(MCTC_bits::hadron);
  }