ATLAS Offline Software
Classes | Public Member Functions | Protected Types | Protected Member Functions | Protected Attributes | Private Types | Private Member Functions | Private Attributes | List of all members
CP::JetQGTaggerBDT Class Reference

#include <JetQGTaggerBDT.h>

Inheritance diagram for CP::JetQGTaggerBDT:
Collaboration diagram for CP::JetQGTaggerBDT:

Classes

struct  Tagger
 Slot-specific TMVA tool and associated variables. More...
 

Public Member Functions

 JetQGTaggerBDT (const std::string &name)
 Constructor. More...
 
virtual StatusCode initialize () override
 Run once at the start of the job to setup everything. More...
 
virtual StatusCode tag (const xAOD::Jet &jet) const override
 IBoostedJetTagger interface. More...
 
virtual StatusCode decorate (const xAOD::JetContainer &jets) const override
 Decorate jet collection with tagging info. More...
 
virtual void print () const
 Print the state of the tool. More...
 
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & evtStore () const
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
 
virtual StatusCode sysInitialize () override
 Perform system initialization for an algorithm. More...
 
virtual StatusCode sysStart () override
 Handle START transition. More...
 
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles. More...
 
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles. More...
 
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T > &t)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc="none")
 Declare a new Gaudi property. More...
 
void updateVHKA (Gaudi::Details::PropertyBase &)
 
MsgStream & msg () const
 
MsgStream & msg (const MSG::Level lvl) const
 
bool msgLvl (const MSG::Level lvl) const
 
virtual StatusCode modify (xAOD::JetContainer &jets) const override final
 Concrete implementation of the function inherited from IJetModifier. More...
 

Protected Types

enum  TAGCLASS { Unknown, WBoson, ZBoson, TopQuark }
 TAGTYPE enum. More...
 

Protected Member Functions

StatusCode getConfigReader ()
 Get configReader StatusCode. More...
 
StatusCode resetCuts (asg::AcceptData &acceptData) const
 Reset cuts. More...
 
bool passKinRange (const xAOD::Jet &jet) const
 Check if jet passes kinematic constraints. More...
 
StatusCode checkKinRange (const xAOD::Jet &jet, asg::AcceptData &acceptData) const
 Check and record if jet passes kinematic constraints. More...
 
int calculateJSSRatios (const xAOD::Jet &jet) const
 Calculate JSS moment ratios in case they are not already saved TODO: Remove this once JSSMomentTools is modified to take const jets. More...
 
int findPV () const
 Find the PV (to be used for Ntrk) More...
 
StatusCode GetUnGroomTracks (const xAOD::Jet &jet, int indexPV) const
 Retrieve Ntrk variable from the ungroomed parent jet. More...
 
StatusCode getWeight (const xAOD::Jet &jet, bool passSel, asg::AcceptData &acceptData) const
 Get SF weight. More...
 
std::pair< double, double > getSF (const xAOD::Jet &jet, const std::string &truthLabelStr) const
 Get scale factor and efficiency. More...
 
std::string getTruthLabelStr (const xAOD::Jet &jet, asg::AcceptData &acceptData) const
 
void printCuts () const
 Print configured cuts. More...
 
void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
 
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed. More...
 

Protected Attributes

TEnv m_configReader
 TEnv instance to read config files. More...
 
asg::AcceptInfo m_acceptInfo
 Object that stores the results for a jet. More...
 
SG::ReadHandleKey< xAOD::JetContainerm_jetContainer_key {this, "JetContainer", "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets", "SG key for input jet container"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decTaggedKey {this, "TaggedName", "Tagged", "SG key for Tagged"}
 WriteDecorHandle keys for tagging bools. More...
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decValidPtRangeHighKey {this, "ValidPtRangeHighName", "ValidPtRangeHigh", "SG key for ValidPtRangeHigh"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decValidPtRangeLowKey {this, "ValidPtRangeLowName", "ValidPtRangeLow", "SG key for ValidPtRangeLow"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decValidEtaRangeKey {this, "ValidEtaRangeName", "ValidEtaRange", "SG key for ValidEtaRange"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decValidKinRangeKey {this, "ValidKinRangeName", "ValidKinRange", "SG key for ValidKinRange"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decValidJetContentKey {this, "ValidJetContentName", "ValidJetContent", "SG key for ValidJetContent"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decValidEventContentKey {this, "ValidEventContentName", "ValidEventContent", "SG key for ValidEventContent"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decPassMassKey {this, "PassMassName", "PassMass", "SG key for PassMass"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decPassScoreKey {this, "PassScoreName", "PassScore", "SG key for PassScore"}
 
const int m_nWarnMax = 10
 Maximum number of warnings. More...
 
std::atomic< int > m_nWarnVar
 Warning counters. More...
 
bool m_ptGeV = false
 Flag to indicate units pT is defined in Set to false by default. More...
 
bool m_useMassCut = false
 Flag to indicate if mass window cut is used. More...
 
bool m_useScoreCut = false
 Flag to indicate if a discriminant score is used. More...
 
TAGCLASS m_tagClass
 
std::string m_containerName
 Configurable members. More...
 
std::string m_weightConfigPath
 Path to the SF configuration root file. More...
 
std::string m_configFile
 Configuration file name. More...
 
std::string m_calibArea
 Location where config files live on cvmfs. More...
 
std::string m_kerasConfigFileName
 Keras configurations for ML taggers. More...
 
std::string m_kerasConfigFilePath
 
std::string m_kerasConfigOutputName
 
std::string m_kerasCalibArea
 
std::string m_tmvaConfigFileName
 TMVA configurations for BDT taggers. More...
 
std::string m_tmvaConfigFilePath
 
std::string m_tmvaCalibArea
 
std::string m_wkpt
 Tagger information. More...
 
std::string m_tagType
 
float m_jetPtMin
 Kinematic bounds for the jet - the units are controlled by m_ptGeV. More...
 
float m_jetPtMax
 
float m_jetEtaMax
 
bool m_useECFG = false
 Flags controlling whether generalized ECF moments or L-series ratios are needed TODO: Implement the functionality controlled by these. More...
 
bool m_useLSeries = false
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decTau21WTAKey {this, "Tau21WTANameOut", "Tau21_wta", "SG key for Tau21_wta"}
 WriteDecorHandle keys for JSS moments. More...
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decTau32WTAKey {this, "Tau32WTANameOut", "Tau32_wta", "SG key for Tau32_wta"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decTau42WTAKey {this, "Tau42WTANameOut", "Tau42_wta", "SG key for Tau42_wta"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decC2Key {this, "C2NameOut", "C2", "SG key for C2"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decD2Key {this, "D2NameOut", "D2", "SG key for D2"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decE3Key {this, "e3NameOut", "e3", "SG key for e3"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decL2Key {this, "L2NameOut", "L2", "SG key for L2"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decL3Key {this, "L3NameOut", "L3", "SG key for L3"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readTau1WTAKey {this, "Tau1WTAName", "Tau1_wta", "SG key for Tau1_wta"}
 ReadDecorHandle keys for JSS moments. More...
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readTau2WTAKey {this, "Tau2WTAName", "Tau2_wta", "SG key for Tau2_wta"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readTau3WTAKey {this, "Tau3WTAName", "Tau3_wta", "SG key for Tau3_wta"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readTau4WTAKey {this, "Tau4WTAName", "Tau4_wta", "SG key for Tau4_wta"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readTau21WTAKey {this, "Tau21WTAName", "Tau21_wta", "SG key for Tau21_wta"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readTau32WTAKey {this, "Tau32WTAName", "Tau32_wta", "SG key for Tau32_wta"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readTau42WTAKey {this, "Tau42WTAName", "Tau42_wta", "SG key for Tau42_wta"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readECF1Key {this, "ECF1Name", "ECF1", "SG key for ECF1"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readECF2Key {this, "ECF2Name", "ECF2", "SG key for ECF2"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readECF3Key {this, "ECF3Name", "ECF3", "SG key for ECF3"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readC2Key {this, "C2Name", "C2", "SG key for C2"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readD2Key {this, "D2Name", "D2", "SG key for D2"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readE3Key {this, "e3Name", "e3", "SG key for e3"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readSplit12Key {this, "Split12Name", "Split12", "SG key for Split12"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readSplit23Key {this, "Split23Name", "Split23", "SG key for Split23"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readQwKey {this, "QwName", "Qw", "SG key for Qw"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readThrustMajKey {this, "ThrustMajName", "ThrustMaj", "SG key for ThrustMaj"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readSphericityKey {this, "SphericityName", "Sphericity", "SG key for Sphericity"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readL2Key {this, "L2Name", "L2", "SG key for L2"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readL3Key {this, "L3Name", "L3", "SG key for L3"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readECFG331Key {this, "ECFG331Name", "ECFG_3_3_1", "SG key for ECFG_3_3_1"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readECFG311Key {this, "ECFG311Name", "ECFG_3_1_1", "SG key for ECFG_3_1_1"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readECFG212Key {this, "ECFG212Name", "ECFG_2_1_2", "SG key for ECFG_2_1_2"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readParentKey {this, "ParentName", "Parent", "SG key for Parent"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decNtrk500Key {this, "Ntrk500NameOut", "ParentJetNTrkPt500", "SG key for Ntrk500 from ungroomed jet"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readNtrk500Key {this, "Ntrk500Name", "ParentJetNTrkPt500", "SG key for Ntrk500 from ungroomed jet"}
 
bool m_suppressOutputDependence
 
std::string m_strMassCutLow
 Strings for cut functions. More...
 
std::string m_strMassCutHigh
 
std::string m_strScoreCut
 
std::unique_ptr< TF1 > m_funcMassCutLow
 TF1 for cut functions. More...
 
std::unique_ptr< TF1 > m_funcMassCutHigh
 
std::unique_ptr< TF1 > m_funcScoreCut
 
std::string m_decorationName
 Decoration name. More...
 
bool m_calcSF
 Flag to calculate scale factor. More...
 
bool m_isMC
 
std::string m_weightDecorationName
 String for scale factor decoration names. More...
 
std::string m_weightFileName
 
std::string m_weightHistogramName
 
std::string m_efficiencyHistogramName
 
std::string m_weightFlavors
 
std::unique_ptr< TFile > m_weightConfig
 Histograms for scale factors. More...
 
std::map< std::string, std::unique_ptr< TH2D > > m_weightHistograms
 
std::map< std::string, std::unique_ptr< TH2D > > m_efficiencyHistograms
 
bool m_truthLabelUseTRUTH3
 Truth label options. More...
 
std::string m_truthParticleContainerName
 
std::string m_truthBosonContainerName
 
std::string m_truthTopQuarkContainerName
 
std::string m_truthLabelName
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readTruthLabelKey {this, "truthLabelName", "truthLabel", "SG key for truthLabel"}
 Truth label ReadDecorHandle key. More...
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decCutMLowKey {this, "CutMLowName", "Cut_mlow", "SG key for Cut_mlow"}
 WriteDecorHandle keys for cut values. More...
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decCutMHighKey {this, "CutMHighName", "Cut_mhigh", "SG key for Cut_mhigh"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decScoreCutKey {this, "CutScoreName", "Cut_score", "SG key for Cut_score"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decScoreValueKey {this, "ScoreName", "Score", "SG key for Score"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decWeightKey {this, "weightName", "weight", "SG key for weight"}
 WriteDecorHandle keys for SF. More...
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decEfficiencyKey {this, "efficiencyName", "efficiency", "SG key for efficiency"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decEffSFKey {this, "effSFName", "effSF", "SG key for effSF"}
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decSigeffSFKey {this, "sigeffSFName", "sigeffSF", "SG key for effSF"}
 

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t
 

Private Member Functions

float getScore (const xAOD::Jet &jet, asg::AcceptData &acceptData) const
 Retrieve BDT score. More...
 
bool getJetProperties (const xAOD::Jet &jet, asg::AcceptData &acceptData) const
 Update the jet substructure variables for each jet to use in BDT. More...
 
bool getPrecomputedVariables (const xAOD::Jet &jet, asg::AcceptData &acceptData) const
 
bool calculateVariables (const xAOD::Jet &jet, asg::AcceptData &acceptData) const
 
bool isCorrectNumberOfTracks (int expectedNTracks, int nTracksFromGhostTracks) const
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

SG::SlotSpecificObj< Tagger > m_bdtTagger ATLAS_THREAD_SAFE
 
std::string m_BDTmethod
 
asg::AnaToolHandle< InDet::IInDetTrackSelectionToolm_trkSelectionTool
 
int m_mode
 
SG::WriteDecorHandleKey< xAOD::JetContainerm_decScoreKey {this, "QGTaggerBDTScore", "QGTaggerBDTScore", "SG key for QGTaggerBDTScore"}
 WriteDecorHandle keys. More...
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 22 of file JetQGTaggerBDT.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ TAGCLASS

enum JSSTaggerBase::TAGCLASS
protectedinherited

TAGTYPE enum.

Enumerator
Unknown 
WBoson 
ZBoson 
TopQuark 

Definition at line 99 of file JSSTaggerBase.h.

Constructor & Destructor Documentation

◆ JetQGTaggerBDT()

CP::JetQGTaggerBDT::JetQGTaggerBDT ( const std::string &  name)

Constructor.

Jet kinematics

Replace base class default value for JetEtaMax

Tagger configuration

Overwrite base class default

0 uses the tracks. 1 uses variables from the jets (default)

Definition at line 15 of file JetQGTaggerBDT.cxx.

15  :
17  m_BDTmethod("BDT_method"),
18  m_trkSelectionTool(name+"_trackselectiontool", this)
19  {
20 
22  declareProperty( "JetPtMin", m_jetPtMin = 20000.0);
23  declareProperty( "JetPtMax", m_jetPtMax = 1500000.0);
24  m_jetEtaMax = 2.5;
25 
27  m_calibArea = "BoostedJetTaggers/JetQGTaggerBDT/Oct18/";
28  declareProperty( "UseJetVars", m_mode = 1);
29 
30  m_tmvaConfigFileName="TMVAClassification_BDTQGTagger_Oct18_BDT.weights.xml";
31 
32  }

Member Function Documentation

◆ calculateJSSRatios()

int JSSTaggerBase::calculateJSSRatios ( const xAOD::Jet jet) const
protectedinherited

Calculate JSS moment ratios in case they are not already saved TODO: Remove this once JSSMomentTools is modified to take const jets.

Calculate JSS moment ratios in case they are not already saved These are calculated by hand here because JetSubStructureMomentTools does not operate on const jets.

This should be changed in the future

Create write decor handles

Create read decor handles

WTA N-subjettiness ratios

ECF ratios

Definition at line 416 of file JSSTaggerBase.cxx.

416  {
417 
418  int result = 0;
419 
429 
435 
439 
443 
446 
447 
449  float tau21_wta = -999.0;
450  float tau32_wta = -999.0;
451  float tau42_wta = -999.0;
452 
453  float tau1_wta = readTau1WTA(jet);
454  float tau2_wta = readTau2WTA(jet);
455  float tau3_wta = readTau3WTA(jet);
456  float tau4_wta = -999.0;
457  if(readTau4WTA.isAvailable()){
458  tau4_wta = readTau4WTA(jet);
459  }
460 
461  if ( tau1_wta > 1e-8 ) {
462  tau21_wta = tau2_wta / tau1_wta;
463  }
464  else result = 1;
465 
466  if ( tau2_wta > 1e-8 ) {
467  tau32_wta = tau3_wta / tau2_wta;
468  if(readTau4WTA.isAvailable()){
469  tau42_wta = tau4_wta / tau2_wta;
470  }
471  }
472  else result = 1;
473 
474  decTau21WTA(jet) = tau21_wta;
475  decTau32WTA(jet) = tau32_wta;
476  decTau42WTA(jet) = tau42_wta;
477 
479  float C2 = -999.0;
480  float D2 = -999.0;
481  float e3 = -999.0;
482 
483  float ECF1 = readECF1(jet);
484  float ECF2 = readECF2(jet);
485  float ECF3 = readECF3(jet);
486 
487  if ( ECF2 > 1e-8 ) {
488  C2 = ECF3 * ECF1 / std::pow( ECF2, 2.0 );
489  D2 = ECF3 * std::pow( ECF1, 3.0 ) / std::pow( ECF2, 3.0 );
490  }
491  else result = 1;
492 
493  e3 = ECF3 / std::pow( ECF1, 3.0 );
494 
495  decC2(jet) = C2;
496  decD2(jet) = D2;
497  decE3(jet) = e3;
498 
499  // L-series for UFO top taggers
500  float L2 = -999.0;
501  float L3 = -999.0;
502 
503  if(!readL2.isAvailable()){
504  if(readECFG331.isAvailable() && readECFG212.isAvailable()){
505  if(readECFG212(jet) > 1e-8){
506  L2 = readECFG331(jet) / pow(readECFG212(jet), (3.0/2.0));
507  }
508  else result = 1;
509  }
510  decL2(jet) = L2;
511  }
512 
513  if(!readL3.isAvailable()){
514  if(readECFG331.isAvailable() && readECFG311.isAvailable()){
515  if(readECFG331(jet) > 1e-8){
516  L3 = readECFG311(jet) / pow(readECFG331(jet), (1.0/3.0));
517  }
518  else result = 1;
519  }
520  decL3(jet) = L3;
521  }
522 
523  // TODO: Add ECFG for ANN tagger whenever it is defined
524 
525  return result;
526 
527 }

◆ calculateVariables()

bool CP::JetQGTaggerBDT::calculateVariables ( const xAOD::Jet jet,
asg::AcceptData acceptData 
) const
private

calculate q/g tagging variables from GhostTracks associated to jet some derivations apply slimming to these tracks, which would lead to wrong values. so we compare the number of GhostTracks to NumTrkPt500 (i.e. nTracks) if they are "close enough" we can proceed

take the first vertex in the list that is a primary vertex

NTracks

if NumTrkPt500 is not available, I can't confirm that the number of GhostTracks is correct (i.e. unslimmed)

TrackWidth

if TrackWidthPt500 is not available, we can maybe calculate it from tracks

TrackC1

track selection

calculate TrackC1 (and TrackWidth if necessary)

Definition at line 230 of file JetQGTaggerBDT.cxx.

230  {
235 
236  bool validVars = true;
237  bool isValid = true;
238  const xAOD::Vertex* primvertex {nullptr};
239 
240  const xAOD::VertexContainer* vxCont = nullptr;
241  if ( evtStore()->retrieve( vxCont, "PrimaryVertices" ).isFailure() ) {
242  if ( m_nWarnVar++ < m_nWarnMax ) ATH_MSG_WARNING( "Unable to retrieve primary vertex container PrimaryVertices" );
243  else ATH_MSG_DEBUG( "Unable to retrieve primary vertex container PrimaryVertices" );
244  acceptData.setCutResult("ValidEventContent", false);
245  isValid = false;
246  }
247  else if ( vxCont->empty() ) {
248  if ( m_nWarnVar++ < m_nWarnMax ) ATH_MSG_WARNING( "Event has no primary vertices!" );
249  ATH_MSG_DEBUG( "Event has no primary vertices!" );
250  acceptData.setCutResult("ValidEventContent", false);
251  isValid = false;
252  }
253  else {
254  for ( const auto *vx : *vxCont ) {
256  if ( vx->vertexType()==xAOD::VxType::PriVtx ) {
257  primvertex = vx;
258  break;
259  }
260  }
261  }
262  if ( !primvertex ) isValid = false;
263 
264  if ( !isValid ) {
265  validVars = false;
266  return validVars;
267  }
268 
270  std::vector<int> nTrkVec;
271  if(jet.getAttribute(xAOD::JetAttribute::NumTrkPt500, nTrkVec)){
272  ATH_MSG_DEBUG(nTrkVec.size());
273  m_bdtTagger->ntracks = (float) nTrkVec[primvertex->index()];
274  }
275  else
277  validVars = false;
278 
280  bool undefTrackWidth = false;
281  std::vector<float> trkWidthVec;
282  if(jet.getAttribute(xAOD::JetAttribute::TrackWidthPt500, trkWidthVec)){
283  ATH_MSG_DEBUG(trkWidthVec.size());
284  m_bdtTagger->trackwidth = trkWidthVec[primvertex->index()];
285  }
286  else
288  undefTrackWidth = true;
289  float weightedwidth = 0.;
290 
292  float beta = 0.2;
293  float weightedwidth2 = 0.;
294  float sumPt = 0.;
295 
296  std::vector<const xAOD::TrackParticle*> trackParttmp;
297  if(!jet.getAssociatedObjects("GhostTrack",trackParttmp)){
298  ATH_MSG_ERROR("This jet has no associated objects");
299  validVars = false;
300  }
302  for(unsigned i=trackParttmp.size();i>0; i--){
303  if(!trackParttmp[i-1]){
304  trackParttmp.erase(trackParttmp.begin()+i-1);
305  continue;
306  }
307  const xAOD::TrackParticle* trk = static_cast<const xAOD::TrackParticle*>(trackParttmp[i-1]);
308  bool accept = (trk->pt()>500 &&
310  );
311  if (!accept){
312  trackParttmp.erase(trackParttmp.begin()+i-1);
313  }
314  }
315 
316  if(! isCorrectNumberOfTracks(m_bdtTagger->ntracks,trackParttmp.size())){
317  ATH_MSG_ERROR("Number of ghost associated tracks wrong!");
318  validVars = false;
319  }
320 
322  for(unsigned i=0; i<trackParttmp.size(); i++){
323  double ipt = trackParttmp.at(i)->pt();
324  double ieta = trackParttmp.at(i)->eta();
325  double iphi = trackParttmp.at(i)->phi();
326  sumPt += ipt;
327  if(undefTrackWidth){
328  double deta_i = trackParttmp.at(i)->eta() - jet.eta();
329  double dphi_i = TVector2::Phi_mpi_pi(trackParttmp.at(i)->phi() - jet.phi());
330  double dR_i = sqrt( deta_i*deta_i + dphi_i*dphi_i );
331  weightedwidth += ipt * dR_i;
332  }
333 
334  for(unsigned j=i+1; j<trackParttmp.size(); j++){
335  double deta = ieta - trackParttmp.at(j)->eta();
336  double dphi = TVector2::Phi_mpi_pi(iphi - trackParttmp.at(j)->phi());
337  double dR = sqrt( deta*deta + dphi*dphi );
338  weightedwidth2 += ipt * trackParttmp.at(j)->pt() * pow(dR,beta);
339  }
340  }
341 
342  if(undefTrackWidth) {
343  m_bdtTagger->trackwidth = sumPt>0 ? weightedwidth/sumPt : -0.1;
344  }
345  m_bdtTagger->trackC1 = sumPt>0 ? weightedwidth2/(sumPt*sumPt) : -0.1;
346 
347  return validVars;
348  }

◆ checkKinRange()

StatusCode JSSTaggerBase::checkKinRange ( const xAOD::Jet jet,
asg::AcceptData acceptData 
) const
protectedinherited

Check and record if jet passes kinematic constraints.

Check each kinematic constraint Print warnings using counters

Create write decor handles

Decorate kinematic pass information

Definition at line 370 of file JSSTaggerBase.cxx.

370  {
371 
372  float scale = 1.0;
373  if ( m_ptGeV ) scale = 1.e3;
374 
375  bool passKinRange = true;
376 
379  if ( std::abs(jet.eta()) > m_jetEtaMax ) {
380  ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (|eta| < " << m_jetEtaMax << "). Jet eta = " << jet.eta() );
381  acceptData.setCutResult( "ValidEtaRange", false );
382  passKinRange = false;
383  }
384 
385  if ( jet.pt() < m_jetPtMin * scale ) {
386  ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (pT > " << m_jetPtMin * scale / 1.e3 << "). Jet pT = " << jet.pt() / 1.e3 << " GeV" );
387  acceptData.setCutResult( "ValidPtRangeLow", false );
388  passKinRange = false;
389  }
390 
391  if ( jet.pt() > m_jetPtMax * scale ) {
392  ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (pT < " << m_jetPtMax * scale / 1.e3 << "). Jet pT = " << jet.pt() / 1.e3 << " GeV" );
393  acceptData.setCutResult( "ValidPtRangeHigh", false );
394  passKinRange = false;
395  }
396 
402 
404  decValidPtRangeHigh(jet) = acceptData.getCutResult( "ValidPtRangeHigh" );
405  decValidPtRangeLow(jet) = acceptData.getCutResult( "ValidPtRangeLow" );
406  decValidEtaRange(jet) = acceptData.getCutResult( "ValidEtaRange" );
407  decValidKinRange(jet) = passKinRange;
408 
409  return StatusCode::SUCCESS;
410 
411 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleKeyArrayType  
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleKeyType  
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleType  
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  t,
const SG::NotHandleType  
)
inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleBase hndl,
const std::string &  doc,
const SG::VarHandleType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKey hndl,
const std::string &  doc,
const SG::VarHandleKeyType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType  
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc,
const SG::NotHandleType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

337  {
338  return PBASE::declareProperty(name, property, doc);
339  }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc = "none" 
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

355  {
356  typedef typename SG::HandleClassifier<T>::type htype;
357  return declareProperty (name, property, doc, htype());
358  }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > &  t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
148  }

◆ decorate()

StatusCode JSSTaggerBase::decorate ( const xAOD::JetContainer jets) const
overridevirtualinherited

Decorate jet collection with tagging info.

Loop over jet collection and decorate each jet.

Implements IJetDecorator.

Reimplemented in JSSWTopTaggerANN.

Definition at line 289 of file JSSTaggerBase.cxx.

289  {
290 
291  for ( const auto *jet : jets ) {
292  ATH_CHECK( tag(*jet) );
293  }
294 
295  return StatusCode::SUCCESS;
296 
297 }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ findPV()

int JSSTaggerBase::findPV ( ) const
protectedinherited

Find the PV (to be used for Ntrk)

Definition at line 529 of file JSSTaggerBase.cxx.

529  {
530 
531  int indexPV = -1;
532 
533  const xAOD::VertexContainer* vxCont = nullptr;
534  if ( evtStore()->retrieve( vxCont, "PrimaryVertices" ).isFailure() ) {
535  ATH_MSG_WARNING( "Unable to retrieve primary vertex container PrimaryVertices" );
536  }
537  else {
538  int vtx_counter = 0;
539  for ( const auto *vx : *vxCont ) {
540  if ( vx->vertexType()==xAOD::VxType::PriVtx ) {
541  indexPV = vtx_counter;
542  break;
543  }
544  vtx_counter++;
545  }
546  }
547 
548  return indexPV;
549 
550 }

◆ getConfigReader()

StatusCode JSSTaggerBase::getConfigReader ( )
protectedinherited

Get configReader StatusCode.

Get configReader TEnv.

Check for the existence of the configuration file

https://root.cern.ch/root/roottalk/roottalk02/5332.html

Definition at line 300 of file JSSTaggerBase.cxx.

300  {
301 
302  ATH_MSG_INFO( "Using config file : " << m_configFile );
303 
305  std::string configPath;
306 
307  if ( m_calibArea.compare("Local") == 0 ) {
308  configPath = PathResolverFindCalibFile(("$WorkDir_DIR/data/BoostedJetTaggers/"+m_configFile).c_str());
309  }
310  else if ( m_calibArea.find("eos") != std::string::npos) {
311  configPath = PathResolverFindCalibFile((m_calibArea+"/"+m_configFile).c_str());
312  }
313  else {
314  configPath = PathResolverFindCalibFile(("BoostedJetTaggers/"+m_calibArea+"/"+m_configFile).c_str());
315  }
316 
318  FileStat_t fStats;
319  int fSuccess = gSystem->GetPathInfo(configPath.c_str(), fStats);
320  if ( fSuccess ) {
321  ATH_MSG_ERROR( "Recommendations file " << m_configFile << " could not be found" );
322  return StatusCode::FAILURE;
323  }
324  else {
325  ATH_MSG_DEBUG( "Recommendations file was found : " << configPath );
326  }
327 
328  if ( m_configReader.ReadFile( configPath.c_str(), EEnvLevel(0) ) ) {
329  ATH_MSG_ERROR( "Error while reading config file : "<< configPath );
330  return StatusCode::FAILURE;
331  }
332 
333  return StatusCode::SUCCESS;
334 
335 }

◆ getJetProperties()

bool CP::JetQGTaggerBDT::getJetProperties ( const xAOD::Jet jet,
asg::AcceptData acceptData 
) const
private

Update the jet substructure variables for each jet to use in BDT.

Definition at line 166 of file JetQGTaggerBDT.cxx.

166  {
167  /* Update the jet substructure variables for this jet */
168 
169  m_bdtTagger->pt = jet.pt()/1000.0;
170  m_bdtTagger->eta = jet.eta();
171 
172  ATH_MSG_DEBUG( TString::Format("pT: %g, eta: %g", m_bdtTagger->pt, m_bdtTagger->eta) );
173 
174  m_bdtTagger->ntracks = -1.;
175  m_bdtTagger->trackwidth = -1.;
176  m_bdtTagger->trackC1 = -1.;
177 
178  bool validVars = true;
179 
180  if ( m_mode == 1 ) {
181  validVars = getPrecomputedVariables( jet, acceptData );
182  }
183  else if( m_mode == 0 ) {
184  validVars = calculateVariables( jet, acceptData );
185  }
186 
187  if ( !validVars ) {
188  ATH_MSG_ERROR( "Can't determine QG tagging variables! Try different mode." );
189  }
190 
191  return validVars;
192 
193  }

◆ getKey()

SG::sgkey_t asg::AsgTool::getKey ( const void *  ptr) const
inherited

Get the (hashed) key of an object that is in the event store.

This is a bit of a special one. StoreGateSvc and xAOD::TEvent both provide ways for getting the SG::sgkey_t key for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.

In order to allow tools to efficiently perform this operation, they can use this helper function.

See also
asg::AsgTool::getName
Parameters
ptrThe bare pointer to the object that the event store should know about
Returns
The hashed key of the object in the store. If not found, an invalid (zero) key.

Definition at line 119 of file AsgTool.cxx.

119  {
120 
121 #ifdef XAOD_STANDALONE
122  // In case we use @c xAOD::TEvent, we have a direct function call
123  // for this.
124  return evtStore()->event()->getKey( ptr );
125 #else
126  const SG::DataProxy* proxy = evtStore()->proxy( ptr );
127  return ( proxy == nullptr ? 0 : proxy->sgkey() );
128 #endif // XAOD_STANDALONE
129  }

◆ getName()

const std::string & asg::AsgTool::getName ( const void *  ptr) const
inherited

Get the name of an object that is / should be in the event store.

This is a bit of a special one. StoreGateSvc and xAOD::TEvent both provide ways for getting the std::string name for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.

In order to allow tools to efficiently perform this operation, they can use this helper function.

See also
asg::AsgTool::getKey
Parameters
ptrThe bare pointer to the object that the event store should know about
Returns
The string name of the object in the store. If not found, an empty string.

Definition at line 106 of file AsgTool.cxx.

106  {
107 
108 #ifdef XAOD_STANDALONE
109  // In case we use @c xAOD::TEvent, we have a direct function call
110  // for this.
111  return evtStore()->event()->getName( ptr );
112 #else
113  const SG::DataProxy* proxy = evtStore()->proxy( ptr );
114  static const std::string dummy = "";
115  return ( proxy == nullptr ? dummy : proxy->name() );
116 #endif // XAOD_STANDALONE
117  }

◆ getPrecomputedVariables()

bool CP::JetQGTaggerBDT::getPrecomputedVariables ( const xAOD::Jet jet,
asg::AcceptData acceptData 
) const
private

Definition at line 195 of file JetQGTaggerBDT.cxx.

195  {
196 
197  bool validVars = true;
198 
199  int ntrk = -1;
200  float trkWidth = -1.;
201  float trkC1 = -1.;
202 
203  if ( !jet.getAttribute<int>("DFCommonJets_QGTagger_NTracks", ntrk) ) {
204  if ( m_nWarnVar++ < m_nWarnMax ) ATH_MSG_WARNING( "Unable to retrieve DFCommonJets_QGTagger_NTracks" );
205  else ATH_MSG_DEBUG( "Unable to retrieve DFCommonJets_QGTagger_NTracks" );
206  acceptData.setCutResult("ValidEventContent", false);
207  validVars = false;
208  }
209  if ( !jet.getAttribute<float>("DFCommonJets_QGTagger_TracksWidth", trkWidth) ) {
210  if ( m_nWarnVar++ < m_nWarnMax )ATH_MSG_WARNING( "Unable to retrieve DFCommonJets_QGTagger_TracksWidth" );
211  else ATH_MSG_DEBUG( "Unable to retrieve DFCommonJets_QGTagger_TracksWidth" );
212  acceptData.setCutResult("ValidEventContent", false);
213  validVars = false;
214  }
215  if ( !jet.getAttribute<float>("DFCommonJets_QGTagger_TracksC1", trkC1) ) {
216  if ( m_nWarnVar++ < m_nWarnMax ) ATH_MSG_WARNING( "Unable to retrieve DFCommonJets_QGTagger_TracksC1" );
217  else ATH_MSG_DEBUG( "Unable to retrieve DFCommonJets_QGTagger_TracksC1" );
218  acceptData.setCutResult("ValidEventContent", false);
219  validVars = false;
220  }
221 
222  m_bdtTagger->ntracks = (float) ntrk;
223  m_bdtTagger->trackwidth = trkWidth;
224  m_bdtTagger->trackC1 = trkC1;
225 
226  return validVars;
227 
228  }

◆ getProperty()

template<class T >
const T* asg::AsgTool::getProperty ( const std::string &  name) const
inherited

Get one of the tool's properties.

◆ getScore()

float CP::JetQGTaggerBDT::getScore ( const xAOD::Jet jet,
asg::AcceptData acceptData 
) const
private

Retrieve BDT score.

Load the new values of the variables for this jet

evaluate bdt

Definition at line 150 of file JetQGTaggerBDT.cxx.

150  {
151 
153  bool validVars = getJetProperties( jet, acceptData );
154 
156  float bdt_score(-666.);
157  if ( !validVars ) {
158  ATH_MSG_WARNING( "One (or more) tagger input variable has an undefined value (NaN), setting score to -666" );
159  return bdt_score;
160  }
161  bdt_score = m_bdtTagger->tmva->EvaluateMVA( m_BDTmethod.c_str() );
162 
163  return bdt_score;
164  }

◆ getSF()

std::pair< double, double > JSSTaggerBase::getSF ( const xAOD::Jet jet,
const std::string &  truthLabelStr 
) const
protectedinherited

Get scale factor and efficiency.

To apply W-tagging efficiency SF to Z-tagger, jet mass is shifted by 10.803 GeV

Definition at line 670 of file JSSTaggerBase.cxx.

670  {
671 
672  if ( !passKinRange(jet) ) return std::make_pair( 1.0, 1.0 );
673 
674 
675  double logmOverPt = std::log(jet.m()/jet.pt());
676  if ( m_decorationName.find("SmoothZ") != std::string::npos ||
677  m_decorationName.find("SmoothInclusiveZ") != std::string::npos ) {
679  const double WtoZmassShift = 10803;
680  logmOverPt = std::log((jet.m()-WtoZmassShift)/jet.pt());
681  }
682 
683  if ( logmOverPt > 0 ) logmOverPt = 0;
684 
685  double SF = 1.0;
686  double eff = 1.0;
687 
688  if ( m_weightHistograms.count(truthLabelStr.c_str()) ) {
689 
690  int pt_mPt_bin = (m_weightHistograms.find(truthLabelStr.c_str())->second)->FindBin(jet.pt()*0.001, logmOverPt);
691  SF = (m_weightHistograms.find(truthLabelStr.c_str())->second)->GetBinContent(pt_mPt_bin);
692 
693  if ( !m_efficiencyHistogramName.empty() ) {
694  eff = (m_efficiencyHistograms.find(truthLabelStr.c_str())->second)->GetBinContent(pt_mPt_bin);
695  }
696 
697  }
698  else {
699  // set the efficiency for "Other" category to be the signal efficiency
700  std::string signal_truthLabel="";
701  if ( m_weightHistograms.count("t_qqb") ) {
702  signal_truthLabel="t_qqb";
703  }else if ( m_weightHistograms.count("V_qq") ){
704  signal_truthLabel="V_qq";
705  }else if ( m_weightHistograms.count("t") ){
706  signal_truthLabel="t";
707  }
708  if ( !signal_truthLabel.empty() && !m_efficiencyHistogramName.empty() ){
709  int pt_mPt_bin = (m_weightHistograms.find(signal_truthLabel.c_str())->second)->FindBin(jet.pt()*0.001, logmOverPt);
710  eff = (m_efficiencyHistograms.find(signal_truthLabel.c_str())->second)->GetBinContent(pt_mPt_bin);
711  }
712 
713  return std::make_pair( 1.0, eff );
714  }
715 
716  if ( SF < 1e-3 ) {
717  ATH_MSG_DEBUG( "(pt, m/pt) (" << jet.pt()/1.e3 << ", " << jet.m()/jet.pt() << ") is out of range for SF calculation. Returning 1.0" );
718  return std::make_pair( 1.0, 1.0 );
719  }
720  else {
721  return std::make_pair( SF, eff );
722  }
723 
724 }

◆ getTruthLabelStr()

std::string JSSTaggerBase::getTruthLabelStr ( const xAOD::Jet jet,
asg::AcceptData acceptData 
) const
protectedinherited

Truth label string

Truth label value

Contained top tagger

Contained top

QCD

TCC W/Z 2-var tagger

Top

W/Z

QCD

Pass mass and D2

Fail mass, pass D2

Pass mass, fail D2

Fail mass and D2

W/Z tagger

Top

W/Z

QCD

Top

QCD

Definition at line 728 of file JSSTaggerBase.cxx.

728  {
730  std::string truthLabelStr;
731 
734  LargeRJetTruthLabel::TypeEnum jetContainment = LargeRJetTruthLabel::intToEnum(readTruthLabel(jet));
735 
737  if ( m_weightHistograms.count("t_qqb") ) {
738 
740  if ( jetContainment==LargeRJetTruthLabel::tqqb ) {
741  truthLabelStr = "t_qqb";
742  }
744  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
745  truthLabelStr = "q";
746  }
747 
748  }
750  else if ( m_weightHistograms.count("V_qq_passMpassD2") ) {
751 
753  if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t ) {
754  truthLabelStr = "t_";
755  }
757  else if ( jetContainment==LargeRJetTruthLabel::Wqq || jetContainment==LargeRJetTruthLabel::Zqq || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
758  truthLabelStr = "V_qq_";
759  }
761  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
762  truthLabelStr = "q_";
763  }
764 
766  if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") && acceptData.getCutResult("PassD2") ) {
767  truthLabelStr += "passMpassD2";
768  }
770  else if ( !(acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh")) && acceptData.getCutResult("PassD2") ) {
771  truthLabelStr += "failMpassD2";
772  }
774  else if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") && !acceptData.getCutResult("PassD2") ) {
775  truthLabelStr += "passMfailD2";
776  }
778  else{
779  truthLabelStr += "failMfailD2";
780  }
781 
782  }
783 
785  else if ( m_weightHistograms.count("V_qq") ) {
786 
788  if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t ) {
789  truthLabelStr = "t";
790  }
792  else if ( jetContainment==LargeRJetTruthLabel::Wqq || jetContainment==LargeRJetTruthLabel::Zqq || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
793  truthLabelStr = "V_qq";
794  }
796  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
797  truthLabelStr = "q";
798  }
799  }
800 
801  // inclusive top tagger
802  else {
804  if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
805  truthLabelStr = "t";
806  }
808  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
809  truthLabelStr = "q";
810  }
811 
812  }
813 
814  return truthLabelStr;
815 }

◆ GetUnGroomTracks()

StatusCode JSSTaggerBase::GetUnGroomTracks ( const xAOD::Jet jet,
int  indexPV 
) const
protectedinherited

Retrieve Ntrk variable from the ungroomed parent jet.

Retrieve the Ntrk variable from the ungroomed parent jet.

Definition at line 553 of file JSSTaggerBase.cxx.

553  {
554 
556 
558 
559  const xAOD::Jet * ungroomedJet = nullptr;
560  if ( readParent.isAvailable() ) {
561  const ElementLink<xAOD::JetContainer>& linkToUngroomed = readParent(jet);
562  if ( linkToUngroomed.isValid() ) {
563  ungroomedJet = *linkToUngroomed;
564 
565  static const SG::AuxElement::ConstAccessor< std::vector<int> >acc_Ntrk("NumTrkPt500");
566 
567  if ( acc_Ntrk.isAvailable(*ungroomedJet) ) {
568 
569  const std::vector<int>& NTrkPt500 = acc_Ntrk(*ungroomedJet);
570 
571  int jet_ntrk = NTrkPt500.at(indexPV);
572  decNtrk500(jet) = jet_ntrk;
573 
574  }
575  else {
576  ATH_MSG_ERROR("WARNING: Unable to retrieve Ntrk of the ungroomed parent jet. Please make sure this variable is in your derivations!!!");
577  return StatusCode::FAILURE;
578  }
579  }
580  else {
581  ATH_MSG_ERROR("WARNING: Unable to retrieve the parent ungroomed jet. Please make sure this variable is in your derivations!!!");
582  return StatusCode::FAILURE;
583  }
584  }
585  else {
586  ATH_MSG_ERROR("WARNING: Unable to retrieve the link to the parent ungroomed jet. Please make sure this variable is in your derivations!!!");
587  return StatusCode::FAILURE;
588  }
589 
590  return StatusCode::SUCCESS;
591 
592 }

◆ getWeight()

StatusCode JSSTaggerBase::getWeight ( const xAOD::Jet jet,
bool  passSel,
asg::AcceptData acceptData 
) const
protectedinherited

Get SF weight.

Inefficiency SF is directly used

Efficiency SF

Calculate inefficiency SF

If inefficiency SF is not available, SF is always 1.0

Create write decor handles

Decorate values

Definition at line 596 of file JSSTaggerBase.cxx.

596  {
597  if ( !m_calcSF ) return StatusCode::SUCCESS;
598 
599  float weight = 1.0;
600  float effSF = 1.0;
601  float sigeffSF = 1.0;
602  float efficiency = 1.0;
603 
604  if ( m_isMC ) {
605 
606  std::string truthLabelStr = getTruthLabelStr( jet, acceptData );
607  std::tie(effSF, efficiency) = getSF( jet, truthLabelStr );
608 
609  // calculate signal efficiency SF
610  if ( m_weightHistograms.count("t_qqb") ) {
611  sigeffSF = getSF(jet, "t_qqb").first;
612  } else if ( m_weightHistograms.count("V_qq") ) {
613  sigeffSF = getSF(jet, "V_qq").first;
614  } else if ( m_weightHistograms.count("t") ){
615  sigeffSF = getSF(jet, "t").first;
616  } else {
617  sigeffSF = 1.0;
618  }
619 
621  if ( m_weightFlavors.find("fail") != std::string::npos ) {
622  weight = effSF;
623  }
624 
625  else {
626 
628  if ( passSel ) {
629  weight = effSF;
630  }
631 
633  else {
635  if ( m_efficiencyHistogramName.empty() ) {
636  weight = 1.0;
637  }
638  else if ( efficiency < 1.0 ) {
639  weight = ( 1. - effSF * efficiency ) / ( 1. - efficiency );
640  }
641  else {
642  weight = 1.0;
643  }
644  }
645  }
646 
647  }
648 
649  else {
650  weight = 1.0;
651  }
652 
658 
660  decWeight(jet) = weight;
661  decEfficiency(jet) = efficiency;
662  decEffSF(jet) = effSF;
663  decSigeffSF(jet) = sigeffSF;
664 
665  return StatusCode::SUCCESS;
666 
667 }

◆ initialize()

StatusCode CP::JetQGTaggerBDT::initialize ( )
overridevirtual

Run once at the start of the job to setup everything.

Get configReader

if the calibArea is specified to be "Local" then it looks in the same place as the top level configs

convert the JSON config file name to the full path

get the config file from CVMFS necessary because xml files are too large to house on the data space

Make sure score cut string is not empty

set up InDet selection tool

read json file for DNN weights

– Initialize TMVA for BDTs for each slot

configure the bdt

Initialize WriteDecorHandle keys

Call base class initialize

Reimplemented from JSSTaggerBase.

Definition at line 34 of file JetQGTaggerBDT.cxx.

34  {
35 
36  ATH_MSG_INFO( "Initializing JetQGTaggerBDT tool" );
37 
38  if( ! m_configFile.empty() ) {
39 
42 
43  // get the name/path of the JSON config
44  m_tmvaConfigFileName = m_configReader.GetValue("TMVAConfigFile" ,"");
45 
46  }
47 
49  if( m_calibArea.empty() ){
50  ATH_MSG_ERROR( "You need to specify where the calibArea is as either being Local or on CVMFS" );
51  return StatusCode::FAILURE;
52  }
53  else if(m_calibArea.compare("Local")==0){
54  std::string localCalibArea = "BoostedJetTaggers/share/JetQGTaggerBDT/";
55  ATH_MSG_INFO( "Using Local calibArea " << localCalibArea );
58  }
59  else{
60  ATH_MSG_INFO( "Using CVMFS calibArea" );
64  }
65 
67  m_strScoreCut = m_configReader.GetValue("ScoreCut" ,"");
68  if(m_strScoreCut.empty()){
69  ATH_MSG_ERROR( "Score cut function is empty!" );
70  return StatusCode::FAILURE;
71  }
72  else ATH_MSG_INFO( "scoreCut: " << m_strScoreCut );
73 
76  ANA_CHECK( m_trkSelectionTool.setProperty( "CutLevel", "Loose" ) );
78 
80  ATH_MSG_INFO( "BDT Tagger configured with: " << m_tmvaConfigFilePath );
81 
83  TMVA::Tools::Instance();
84  for (auto& tagger : m_bdtTagger) {
85  tagger.tmva = std::make_unique<TMVA::Reader>( "!Color:!Silent" );
86 
87  tagger.tmva->AddVariable( "NTracks", &tagger.ntracks );
88  tagger.tmva->AddVariable( "TrackWidth", &tagger.trackwidth );
89  tagger.tmva->AddVariable( "JetPt", &tagger.pt );
90  tagger.tmva->AddVariable( "JetEta", &tagger.eta );
91  tagger.tmva->AddVariable( "TrackC1", &tagger.trackC1 );
92 
94  tagger.tmva->BookMVA( m_BDTmethod.c_str(), m_tmvaConfigFilePath.c_str() );
95  }
96 
99 
100  ATH_CHECK( m_decScoreKey.initialize() );
101 
104 
105  return StatusCode::SUCCESS;
106 
107  }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const
overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ isCorrectNumberOfTracks()

bool CP::JetQGTaggerBDT::isCorrectNumberOfTracks ( int  expectedNTracks,
int  nTracksFromGhostTracks 
) const
private

some derivations do not store all tracks associated to the jet. In this case the calculation of the tagging variables will be wrong. The requirements are fairly loose, because a few tracks may get lost in the derivation production. But it will fail quickly if the too many tracks were slimmed away.

Definition at line 350 of file JetQGTaggerBDT.cxx.

350  {
355  if(nTracksFromGhostTracks == 0){
356  if(expectedNTracks == 0)
357  return true;
358  return abs(expectedNTracks-nTracksFromGhostTracks) < 3;
359  }else if(expectedNTracks/nTracksFromGhostTracks < 0.5 && abs(expectedNTracks-nTracksFromGhostTracks) > 5){
360  return false;
361  }
362  return true;
363  }

◆ modify()

virtual StatusCode IJetDecorator::modify ( xAOD::JetContainer jets) const
inlinefinaloverridevirtualinherited

Concrete implementation of the function inherited from IJetModifier.

Implements IJetModifier.

Definition at line 32 of file IJetDecorator.h.

32 {return decorate(jets);};

◆ msg() [1/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24  {
25  return this->msgStream();
26  }

◆ msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level  lvl) const
inlineinherited

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

◆ msg_level_name()

const std::string & asg::AsgTool::msg_level_name ( ) const
inherited

A deprecated function for getting the message level's name.

Instead of using this, weirdly named function, user code should get the string name of the current minimum message level (in case they really need it...), with:

MSG::name( msg().level() )

This function's name doesn't follow the ATLAS coding rules, and as such will be removed in the not too distant future.

Returns
The string name of the current minimum message level that's printed

Definition at line 101 of file AsgTool.cxx.

101  {
102 
103  return MSG::name( msg().level() );
104  }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level  lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30  {
31  return this->msgLevel(lvl);
32  }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const
overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ passKinRange()

bool JSSTaggerBase::passKinRange ( const xAOD::Jet jet) const
protectedinherited

Check if jet passes kinematic constraints.

Definition at line 356 of file JSSTaggerBase.cxx.

356  {
357 
358  float scale = 1.0;
359  if ( m_ptGeV ) scale = 1.e3;
360 
361  if ( jet.pt() < m_jetPtMin * scale ) return false;
362  if ( jet.pt() > m_jetPtMax * scale ) return false;
363  if ( std::abs( jet.eta() ) > m_jetEtaMax ) return false;
364 
365  return true;
366 
367 }

◆ print()

void asg::AsgTool::print ( ) const
virtualinherited

◆ printCuts()

void JSSTaggerBase::printCuts ( ) const
protectedinherited

Print configured cuts.

Definition at line 817 of file JSSTaggerBase.cxx.

817  {
818  ATH_MSG_INFO( "After tagging, you will have access to the following cuts as an asg::AcceptData : (<NCut>) <cut> : <description>)" );
819  int nCuts = m_acceptInfo.getNCuts();
820  for ( int iCut=0; iCut < nCuts; iCut++ ) {
821  std::string cut_string = "";
822  cut_string += " (";
823  cut_string += std::to_string(iCut);
824  cut_string += ") ";
825  cut_string += m_acceptInfo.getCutName(iCut).data();
826  cut_string += " : ";
827  cut_string += m_acceptInfo.getCutDescription(iCut).data();
828  ATH_MSG_INFO( cut_string );
829  }
830 }

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T &  h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381  {
382  h.renounce();
383  PBASE::renounce (h);
384  }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ resetCuts()

StatusCode JSSTaggerBase::resetCuts ( asg::AcceptData acceptData) const
protectedinherited

Reset cuts.

Reset the AcceptData cut results to false

Initialize common cuts to true by default

Definition at line 338 of file JSSTaggerBase.cxx.

338  {
339 
341  acceptData.clear();
342 
344  acceptData.setCutResult( "ValidJetContent", true );
345  acceptData.setCutResult( "ValidEventContent", true );
346 
347  acceptData.setCutResult( "ValidPtRangeHigh", true );
348  acceptData.setCutResult( "ValidPtRangeLow" , true );
349  acceptData.setCutResult( "ValidEtaRange" , true );
350 
351  return StatusCode::SUCCESS;
352 
353 }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )
overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ tag()

StatusCode CP::JetQGTaggerBDT::tag ( const xAOD::Jet jet) const
overridevirtual

IBoostedJetTagger interface.

Create asg::AcceptData object

Reset the AcceptData cut results

Check basic kinematic selection

Create WriteDecorHandles

TODO: Is this actually needed?

Get BDT score

Get cut from cut function

Decorate score to jet

Check score against cut and decorate result

Implements JSSTaggerBase.

Definition at line 109 of file JetQGTaggerBDT.cxx.

109  {
110 
111  ATH_MSG_DEBUG( "Obtaining BDT QG result" );
112 
114  asg::AcceptData acceptData( &m_acceptInfo );
115 
117  ATH_CHECK( resetCuts( acceptData ) );
118 
120  ATH_CHECK( checkKinRange( jet, acceptData ) );
121 
125 
127  if ( !passKinRange(jet) ) {
128  decTagged(jet) = false;
129  return StatusCode::SUCCESS;
130  }
131 
133  float jet_score = getScore( jet, acceptData );
134  ATH_MSG_DEBUG(TString::Format("jet score %g",jet_score) );
135 
137  float cut = m_funcScoreCut->Eval(jet.pt()/1000.);
138 
140  decScore(jet) = jet_score;
141 
143  if ( jet_score < cut ) acceptData.setCutResult("QuarkJetTag", true);
144  decTagged(jet) = acceptData.getCutResult( "QuarkJetTag" );
145 
146  return StatusCode::SUCCESS;
147 
148  }

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase &  )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308  {
309  // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310  // << " size: " << m_vhka.size() << endmsg;
311  for (auto &a : m_vhka) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

Member Data Documentation

◆ ATLAS_THREAD_SAFE

SG::SlotSpecificObj<Tagger> m_bdtTagger CP::JetQGTaggerBDT::ATLAS_THREAD_SAFE
mutableprivate

Definition at line 61 of file JetQGTaggerBDT.h.

◆ m_acceptInfo

asg::AcceptInfo JSSTaggerBase::m_acceptInfo
protectedinherited

Object that stores the results for a jet.

Definition at line 65 of file JSSTaggerBase.h.

◆ m_BDTmethod

std::string CP::JetQGTaggerBDT::m_BDTmethod
private

Definition at line 62 of file JetQGTaggerBDT.h.

◆ m_calcSF

bool JSSTaggerBase::m_calcSF
protectedinherited

Flag to calculate scale factor.

Definition at line 205 of file JSSTaggerBase.h.

◆ m_calibArea

std::string JSSTaggerBase::m_calibArea
protectedinherited

Location where config files live on cvmfs.

Definition at line 114 of file JSSTaggerBase.h.

◆ m_configFile

std::string JSSTaggerBase::m_configFile
protectedinherited

Configuration file name.

Definition at line 111 of file JSSTaggerBase.h.

◆ m_configReader

TEnv JSSTaggerBase::m_configReader
protectedinherited

TEnv instance to read config files.

Definition at line 62 of file JSSTaggerBase.h.

◆ m_containerName

std::string JSSTaggerBase::m_containerName
protectedinherited

Configurable members.

Jet container name

Definition at line 105 of file JSSTaggerBase.h.

◆ m_decC2Key

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decC2Key {this, "C2NameOut", "C2", "SG key for C2"}
protectedinherited

Definition at line 146 of file JSSTaggerBase.h.

◆ m_decCutMHighKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decCutMHighKey {this, "CutMHighName", "Cut_mhigh", "SG key for Cut_mhigh"}
protectedinherited

Definition at line 232 of file JSSTaggerBase.h.

◆ m_decCutMLowKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decCutMLowKey {this, "CutMLowName", "Cut_mlow", "SG key for Cut_mlow"}
protectedinherited

WriteDecorHandle keys for cut values.

Definition at line 231 of file JSSTaggerBase.h.

◆ m_decD2Key

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decD2Key {this, "D2NameOut", "D2", "SG key for D2"}
protectedinherited

Definition at line 147 of file JSSTaggerBase.h.

◆ m_decE3Key

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decE3Key {this, "e3NameOut", "e3", "SG key for e3"}
protectedinherited

Definition at line 148 of file JSSTaggerBase.h.

◆ m_decEfficiencyKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decEfficiencyKey {this, "efficiencyName", "efficiency", "SG key for efficiency"}
protectedinherited

Definition at line 238 of file JSSTaggerBase.h.

◆ m_decEffSFKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decEffSFKey {this, "effSFName", "effSF", "SG key for effSF"}
protectedinherited

Definition at line 239 of file JSSTaggerBase.h.

◆ m_decL2Key

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decL2Key {this, "L2NameOut", "L2", "SG key for L2"}
protectedinherited

Definition at line 149 of file JSSTaggerBase.h.

◆ m_decL3Key

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decL3Key {this, "L3NameOut", "L3", "SG key for L3"}
protectedinherited

Definition at line 150 of file JSSTaggerBase.h.

◆ m_decNtrk500Key

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decNtrk500Key {this, "Ntrk500NameOut", "ParentJetNTrkPt500", "SG key for Ntrk500 from ungroomed jet"}
protectedinherited

Definition at line 186 of file JSSTaggerBase.h.

◆ m_decorationName

std::string JSSTaggerBase::m_decorationName
protectedinherited

Decoration name.

Definition at line 202 of file JSSTaggerBase.h.

◆ m_decPassMassKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decPassMassKey {this, "PassMassName", "PassMass", "SG key for PassMass"}
protectedinherited

Definition at line 79 of file JSSTaggerBase.h.

◆ m_decPassScoreKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decPassScoreKey {this, "PassScoreName", "PassScore", "SG key for PassScore"}
protectedinherited

Definition at line 80 of file JSSTaggerBase.h.

◆ m_decScoreCutKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decScoreCutKey {this, "CutScoreName", "Cut_score", "SG key for Cut_score"}
protectedinherited

Definition at line 233 of file JSSTaggerBase.h.

◆ m_decScoreKey

SG::WriteDecorHandleKey<xAOD::JetContainer> CP::JetQGTaggerBDT::m_decScoreKey {this, "QGTaggerBDTScore", "QGTaggerBDTScore", "SG key for QGTaggerBDTScore"}
private

WriteDecorHandle keys.

Definition at line 69 of file JetQGTaggerBDT.h.

◆ m_decScoreValueKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decScoreValueKey {this, "ScoreName", "Score", "SG key for Score"}
protectedinherited

Definition at line 234 of file JSSTaggerBase.h.

◆ m_decSigeffSFKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decSigeffSFKey {this, "sigeffSFName", "sigeffSF", "SG key for effSF"}
protectedinherited

Definition at line 240 of file JSSTaggerBase.h.

◆ m_decTaggedKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decTaggedKey {this, "TaggedName", "Tagged", "SG key for Tagged"}
protectedinherited

WriteDecorHandle keys for tagging bools.

Definition at line 71 of file JSSTaggerBase.h.

◆ m_decTau21WTAKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decTau21WTAKey {this, "Tau21WTANameOut", "Tau21_wta", "SG key for Tau21_wta"}
protectedinherited

WriteDecorHandle keys for JSS moments.

Definition at line 142 of file JSSTaggerBase.h.

◆ m_decTau32WTAKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decTau32WTAKey {this, "Tau32WTANameOut", "Tau32_wta", "SG key for Tau32_wta"}
protectedinherited

Definition at line 143 of file JSSTaggerBase.h.

◆ m_decTau42WTAKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decTau42WTAKey {this, "Tau42WTANameOut", "Tau42_wta", "SG key for Tau42_wta"}
protectedinherited

Definition at line 144 of file JSSTaggerBase.h.

◆ m_decValidEtaRangeKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decValidEtaRangeKey {this, "ValidEtaRangeName", "ValidEtaRange", "SG key for ValidEtaRange"}
protectedinherited

Definition at line 74 of file JSSTaggerBase.h.

◆ m_decValidEventContentKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decValidEventContentKey {this, "ValidEventContentName", "ValidEventContent", "SG key for ValidEventContent"}
protectedinherited

Definition at line 77 of file JSSTaggerBase.h.

◆ m_decValidJetContentKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decValidJetContentKey {this, "ValidJetContentName", "ValidJetContent", "SG key for ValidJetContent"}
protectedinherited

Definition at line 76 of file JSSTaggerBase.h.

◆ m_decValidKinRangeKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decValidKinRangeKey {this, "ValidKinRangeName", "ValidKinRange", "SG key for ValidKinRange"}
protectedinherited

Definition at line 75 of file JSSTaggerBase.h.

◆ m_decValidPtRangeHighKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decValidPtRangeHighKey {this, "ValidPtRangeHighName", "ValidPtRangeHigh", "SG key for ValidPtRangeHigh"}
protectedinherited

Definition at line 72 of file JSSTaggerBase.h.

◆ m_decValidPtRangeLowKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decValidPtRangeLowKey {this, "ValidPtRangeLowName", "ValidPtRangeLow", "SG key for ValidPtRangeLow"}
protectedinherited

Definition at line 73 of file JSSTaggerBase.h.

◆ m_decWeightKey

SG::WriteDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_decWeightKey {this, "weightName", "weight", "SG key for weight"}
protectedinherited

WriteDecorHandle keys for SF.

Definition at line 237 of file JSSTaggerBase.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_efficiencyHistogramName

std::string JSSTaggerBase::m_efficiencyHistogramName
protectedinherited

Definition at line 212 of file JSSTaggerBase.h.

◆ m_efficiencyHistograms

std::map<std::string, std::unique_ptr<TH2D> > JSSTaggerBase::m_efficiencyHistograms
protectedinherited

Definition at line 218 of file JSSTaggerBase.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_funcMassCutHigh

std::unique_ptr<TF1> JSSTaggerBase::m_funcMassCutHigh
protectedinherited

Definition at line 198 of file JSSTaggerBase.h.

◆ m_funcMassCutLow

std::unique_ptr<TF1> JSSTaggerBase::m_funcMassCutLow
protectedinherited

TF1 for cut functions.

Definition at line 197 of file JSSTaggerBase.h.

◆ m_funcScoreCut

std::unique_ptr<TF1> JSSTaggerBase::m_funcScoreCut
protectedinherited

Definition at line 199 of file JSSTaggerBase.h.

◆ m_isMC

bool JSSTaggerBase::m_isMC
protectedinherited

Definition at line 206 of file JSSTaggerBase.h.

◆ m_jetContainer_key

SG::ReadHandleKey<xAOD::JetContainer> JSSTaggerBase::m_jetContainer_key {this, "JetContainer", "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets", "SG key for input jet container"}
protectedinherited

Definition at line 68 of file JSSTaggerBase.h.

◆ m_jetEtaMax

float JSSTaggerBase::m_jetEtaMax
protectedinherited

Definition at line 134 of file JSSTaggerBase.h.

◆ m_jetPtMax

float JSSTaggerBase::m_jetPtMax
protectedinherited

Definition at line 133 of file JSSTaggerBase.h.

◆ m_jetPtMin

float JSSTaggerBase::m_jetPtMin
protectedinherited

Kinematic bounds for the jet - the units are controlled by m_ptGeV.

Definition at line 132 of file JSSTaggerBase.h.

◆ m_kerasCalibArea

std::string JSSTaggerBase::m_kerasCalibArea
protectedinherited

Definition at line 120 of file JSSTaggerBase.h.

◆ m_kerasConfigFileName

std::string JSSTaggerBase::m_kerasConfigFileName
protectedinherited

Keras configurations for ML taggers.

Definition at line 117 of file JSSTaggerBase.h.

◆ m_kerasConfigFilePath

std::string JSSTaggerBase::m_kerasConfigFilePath
protectedinherited

Definition at line 118 of file JSSTaggerBase.h.

◆ m_kerasConfigOutputName

std::string JSSTaggerBase::m_kerasConfigOutputName
protectedinherited

Definition at line 119 of file JSSTaggerBase.h.

◆ m_mode

int CP::JetQGTaggerBDT::m_mode
private

Definition at line 66 of file JetQGTaggerBDT.h.

◆ m_nWarnMax

const int JSSTaggerBase::m_nWarnMax = 10
protectedinherited

Maximum number of warnings.

Definition at line 83 of file JSSTaggerBase.h.

◆ m_nWarnVar

std::atomic<int> JSSTaggerBase::m_nWarnVar
mutableprotectedinherited

Warning counters.

Definition at line 86 of file JSSTaggerBase.h.

◆ m_ptGeV

bool JSSTaggerBase::m_ptGeV = false
protectedinherited

Flag to indicate units pT is defined in Set to false by default.

Definition at line 90 of file JSSTaggerBase.h.

◆ m_readC2Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readC2Key {this, "C2Name", "C2", "SG key for C2"}
protectedinherited

Definition at line 166 of file JSSTaggerBase.h.

◆ m_readD2Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readD2Key {this, "D2Name", "D2", "SG key for D2"}
protectedinherited

Definition at line 167 of file JSSTaggerBase.h.

◆ m_readE3Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readE3Key {this, "e3Name", "e3", "SG key for e3"}
protectedinherited

Definition at line 168 of file JSSTaggerBase.h.

◆ m_readECF1Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readECF1Key {this, "ECF1Name", "ECF1", "SG key for ECF1"}
protectedinherited

Definition at line 162 of file JSSTaggerBase.h.

◆ m_readECF2Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readECF2Key {this, "ECF2Name", "ECF2", "SG key for ECF2"}
protectedinherited

Definition at line 163 of file JSSTaggerBase.h.

◆ m_readECF3Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readECF3Key {this, "ECF3Name", "ECF3", "SG key for ECF3"}
protectedinherited

Definition at line 164 of file JSSTaggerBase.h.

◆ m_readECFG212Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readECFG212Key {this, "ECFG212Name", "ECFG_2_1_2", "SG key for ECFG_2_1_2"}
protectedinherited

Definition at line 183 of file JSSTaggerBase.h.

◆ m_readECFG311Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readECFG311Key {this, "ECFG311Name", "ECFG_3_1_1", "SG key for ECFG_3_1_1"}
protectedinherited

Definition at line 182 of file JSSTaggerBase.h.

◆ m_readECFG331Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readECFG331Key {this, "ECFG331Name", "ECFG_3_3_1", "SG key for ECFG_3_3_1"}
protectedinherited

Definition at line 181 of file JSSTaggerBase.h.

◆ m_readL2Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readL2Key {this, "L2Name", "L2", "SG key for L2"}
protectedinherited

Definition at line 178 of file JSSTaggerBase.h.

◆ m_readL3Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readL3Key {this, "L3Name", "L3", "SG key for L3"}
protectedinherited

Definition at line 179 of file JSSTaggerBase.h.

◆ m_readNtrk500Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readNtrk500Key {this, "Ntrk500Name", "ParentJetNTrkPt500", "SG key for Ntrk500 from ungroomed jet"}
protectedinherited

Definition at line 187 of file JSSTaggerBase.h.

◆ m_readParentKey

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readParentKey {this, "ParentName", "Parent", "SG key for Parent"}
protectedinherited

Definition at line 185 of file JSSTaggerBase.h.

◆ m_readQwKey

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readQwKey {this, "QwName", "Qw", "SG key for Qw"}
protectedinherited

Definition at line 173 of file JSSTaggerBase.h.

◆ m_readSphericityKey

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readSphericityKey {this, "SphericityName", "Sphericity", "SG key for Sphericity"}
protectedinherited

Definition at line 176 of file JSSTaggerBase.h.

◆ m_readSplit12Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readSplit12Key {this, "Split12Name", "Split12", "SG key for Split12"}
protectedinherited

Definition at line 170 of file JSSTaggerBase.h.

◆ m_readSplit23Key

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readSplit23Key {this, "Split23Name", "Split23", "SG key for Split23"}
protectedinherited

Definition at line 171 of file JSSTaggerBase.h.

◆ m_readTau1WTAKey

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readTau1WTAKey {this, "Tau1WTAName", "Tau1_wta", "SG key for Tau1_wta"}
protectedinherited

ReadDecorHandle keys for JSS moments.

Definition at line 153 of file JSSTaggerBase.h.

◆ m_readTau21WTAKey

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readTau21WTAKey {this, "Tau21WTAName", "Tau21_wta", "SG key for Tau21_wta"}
protectedinherited

Definition at line 158 of file JSSTaggerBase.h.

◆ m_readTau2WTAKey

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readTau2WTAKey {this, "Tau2WTAName", "Tau2_wta", "SG key for Tau2_wta"}
protectedinherited

Definition at line 154 of file JSSTaggerBase.h.

◆ m_readTau32WTAKey

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readTau32WTAKey {this, "Tau32WTAName", "Tau32_wta", "SG key for Tau32_wta"}
protectedinherited

Definition at line 159 of file JSSTaggerBase.h.

◆ m_readTau3WTAKey

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readTau3WTAKey {this, "Tau3WTAName", "Tau3_wta", "SG key for Tau3_wta"}
protectedinherited

Definition at line 155 of file JSSTaggerBase.h.

◆ m_readTau42WTAKey

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readTau42WTAKey {this, "Tau42WTAName", "Tau42_wta", "SG key for Tau42_wta"}
protectedinherited

Definition at line 160 of file JSSTaggerBase.h.

◆ m_readTau4WTAKey

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readTau4WTAKey {this, "Tau4WTAName", "Tau4_wta", "SG key for Tau4_wta"}
protectedinherited

Definition at line 156 of file JSSTaggerBase.h.

◆ m_readThrustMajKey

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readThrustMajKey {this, "ThrustMajName", "ThrustMaj", "SG key for ThrustMaj"}
protectedinherited

Definition at line 175 of file JSSTaggerBase.h.

◆ m_readTruthLabelKey

SG::ReadDecorHandleKey<xAOD::JetContainer> JSSTaggerBase::m_readTruthLabelKey {this, "truthLabelName", "truthLabel", "SG key for truthLabel"}
protectedinherited

Truth label ReadDecorHandle key.

Definition at line 228 of file JSSTaggerBase.h.

◆ m_strMassCutHigh

std::string JSSTaggerBase::m_strMassCutHigh
protectedinherited

Definition at line 193 of file JSSTaggerBase.h.

◆ m_strMassCutLow

std::string JSSTaggerBase::m_strMassCutLow
protectedinherited

Strings for cut functions.

Definition at line 192 of file JSSTaggerBase.h.

◆ m_strScoreCut

std::string JSSTaggerBase::m_strScoreCut
protectedinherited

Definition at line 194 of file JSSTaggerBase.h.

◆ m_suppressOutputDependence

bool JSSTaggerBase::m_suppressOutputDependence
protectedinherited

Definition at line 189 of file JSSTaggerBase.h.

◆ m_tagClass

TAGCLASS JSSTaggerBase::m_tagClass
protectedinherited

Definition at line 100 of file JSSTaggerBase.h.

◆ m_tagType

std::string JSSTaggerBase::m_tagType
protectedinherited

Definition at line 129 of file JSSTaggerBase.h.

◆ m_tmvaCalibArea

std::string JSSTaggerBase::m_tmvaCalibArea
protectedinherited

Definition at line 125 of file JSSTaggerBase.h.

◆ m_tmvaConfigFileName

std::string JSSTaggerBase::m_tmvaConfigFileName
protectedinherited

TMVA configurations for BDT taggers.

Definition at line 123 of file JSSTaggerBase.h.

◆ m_tmvaConfigFilePath

std::string JSSTaggerBase::m_tmvaConfigFilePath
protectedinherited

Definition at line 124 of file JSSTaggerBase.h.

◆ m_trkSelectionTool

asg::AnaToolHandle<InDet::IInDetTrackSelectionTool> CP::JetQGTaggerBDT::m_trkSelectionTool
private

Definition at line 64 of file JetQGTaggerBDT.h.

◆ m_truthBosonContainerName

std::string JSSTaggerBase::m_truthBosonContainerName
protectedinherited

Definition at line 223 of file JSSTaggerBase.h.

◆ m_truthLabelName

std::string JSSTaggerBase::m_truthLabelName
protectedinherited

Definition at line 225 of file JSSTaggerBase.h.

◆ m_truthLabelUseTRUTH3

bool JSSTaggerBase::m_truthLabelUseTRUTH3
protectedinherited

Truth label options.

Definition at line 221 of file JSSTaggerBase.h.

◆ m_truthParticleContainerName

std::string JSSTaggerBase::m_truthParticleContainerName
protectedinherited

Definition at line 222 of file JSSTaggerBase.h.

◆ m_truthTopQuarkContainerName

std::string JSSTaggerBase::m_truthTopQuarkContainerName
protectedinherited

Definition at line 224 of file JSSTaggerBase.h.

◆ m_useECFG

bool JSSTaggerBase::m_useECFG = false
protectedinherited

Flags controlling whether generalized ECF moments or L-series ratios are needed TODO: Implement the functionality controlled by these.

Definition at line 138 of file JSSTaggerBase.h.

◆ m_useLSeries

bool JSSTaggerBase::m_useLSeries = false
protectedinherited

Definition at line 139 of file JSSTaggerBase.h.

◆ m_useMassCut

bool JSSTaggerBase::m_useMassCut = false
protectedinherited

Flag to indicate if mass window cut is used.

Definition at line 93 of file JSSTaggerBase.h.

◆ m_useScoreCut

bool JSSTaggerBase::m_useScoreCut = false
protectedinherited

Flag to indicate if a discriminant score is used.

Definition at line 96 of file JSSTaggerBase.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_weightConfig

std::unique_ptr<TFile> JSSTaggerBase::m_weightConfig
protectedinherited

Histograms for scale factors.

Definition at line 216 of file JSSTaggerBase.h.

◆ m_weightConfigPath

std::string JSSTaggerBase::m_weightConfigPath
protectedinherited

Path to the SF configuration root file.

Definition at line 108 of file JSSTaggerBase.h.

◆ m_weightDecorationName

std::string JSSTaggerBase::m_weightDecorationName
protectedinherited

String for scale factor decoration names.

Definition at line 209 of file JSSTaggerBase.h.

◆ m_weightFileName

std::string JSSTaggerBase::m_weightFileName
protectedinherited

Definition at line 210 of file JSSTaggerBase.h.

◆ m_weightFlavors

std::string JSSTaggerBase::m_weightFlavors
protectedinherited

Definition at line 213 of file JSSTaggerBase.h.

◆ m_weightHistogramName

std::string JSSTaggerBase::m_weightHistogramName
protectedinherited

Definition at line 211 of file JSSTaggerBase.h.

◆ m_weightHistograms

std::map<std::string, std::unique_ptr<TH2D> > JSSTaggerBase::m_weightHistograms
protectedinherited

Definition at line 217 of file JSSTaggerBase.h.

◆ m_wkpt

std::string JSSTaggerBase::m_wkpt
protectedinherited

Tagger information.

Definition at line 128 of file JSSTaggerBase.h.


The documentation for this class was generated from the following files:
JSSTaggerBase::m_readParentKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readParentKey
Definition: JSSTaggerBase.h:185
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
LargeRJetTruthLabel::TypeEnum
TypeEnum
Definition: LargeRJetLabelEnum.h:14
CP::JetQGTaggerBDT::m_BDTmethod
std::string m_BDTmethod
Definition: JetQGTaggerBDT.h:62
xAOD::TrackParticle_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition: TrackParticle_v1.cxx:73
python.SystemOfUnits.second
int second
Definition: SystemOfUnits.py:120
JSSTaggerBase::m_readTau4WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau4WTAKey
Definition: JSSTaggerBase.h:156
Trk::L2
@ L2
Definition: AlignModuleList.h:32
JSSTaggerBase::m_readECF1Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF1Key
Definition: JSSTaggerBase.h:162
JSSTaggerBase::getTruthLabelStr
std::string getTruthLabelStr(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Definition: JSSTaggerBase.cxx:728
get_generator_info.result
result
Definition: get_generator_info.py:21
StateLessPT_NewConfig.proxy
proxy
Definition: StateLessPT_NewConfig.py:392
python.PerfMonSerializer.p
def p
Definition: PerfMonSerializer.py:743
JSSTaggerBase::m_funcScoreCut
std::unique_ptr< TF1 > m_funcScoreCut
Definition: JSSTaggerBase.h:199
JSSTaggerBase::ZBoson
@ ZBoson
Definition: JSSTaggerBase.h:99
JSSTaggerBase::resetCuts
StatusCode resetCuts(asg::AcceptData &acceptData) const
Reset cuts.
Definition: JSSTaggerBase.cxx:338
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
JSSTaggerBase::m_readECFG311Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG311Key
Definition: JSSTaggerBase.h:182
JSSTaggerBase::m_jetEtaMax
float m_jetEtaMax
Definition: JSSTaggerBase.h:134
JSSTaggerBase::m_nWarnVar
std::atomic< int > m_nWarnVar
Warning counters.
Definition: JSSTaggerBase.h:86
asg::AnaToolHandle::retrieve
StatusCode retrieve()
initialize the tool
JSSTaggerBase::m_decValidEtaRangeKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidEtaRangeKey
Definition: JSSTaggerBase.h:74
JSSTaggerBase::m_efficiencyHistograms
std::map< std::string, std::unique_ptr< TH2D > > m_efficiencyHistograms
Definition: JSSTaggerBase.h:218
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
JSSTaggerBase::m_configFile
std::string m_configFile
Configuration file name.
Definition: JSSTaggerBase.h:111
CP::JetQGTaggerBDT::m_decScoreKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decScoreKey
WriteDecorHandle keys.
Definition: JetQGTaggerBDT.h:69
JSSTaggerBase::m_configReader
TEnv m_configReader
TEnv instance to read config files.
Definition: JSSTaggerBase.h:62
RootHelpers::FindBin
Int_t FindBin(const TAxis *axis, const double x)
Definition: RootHelpers.cxx:14
StateLessPT_NewConfig.Format
Format
Definition: StateLessPT_NewConfig.py:146
JSSTaggerBase::TopQuark
@ TopQuark
Definition: JSSTaggerBase.h:99
CutsMETMaker::accept
StatusCode accept(const xAOD::Muon *mu)
Definition: CutsMETMaker.cxx:18
JSSTaggerBase::m_nWarnMax
const int m_nWarnMax
Maximum number of warnings.
Definition: JSSTaggerBase.h:83
conifer::pow
constexpr int pow(int x)
Definition: conifer.h:20
JSSTaggerBase::tag
virtual StatusCode tag(const xAOD::Jet &jet) const =0
Decorate single jet with tagging info.
TRT_PAI_gasdata::SF
const float SF[NF]
Cross sections for Fluor.
Definition: TRT_PAI_gasdata.h:285
LargeRJetTruthLabel::other_From_t
@ other_From_t
Definition: LargeRJetLabelEnum.h:22
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
JSSTaggerBase::m_readTau2WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau2WTAKey
Definition: JSSTaggerBase.h:154
JSSTaggerBase::m_jetPtMin
float m_jetPtMin
Kinematic bounds for the jet - the units are controlled by m_ptGeV.
Definition: JSSTaggerBase.h:132
JSSTaggerBase::m_decL3Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL3Key
Definition: JSSTaggerBase.h:150
ANA_CHECK
#define ANA_CHECK(EXP)
check whether the given expression was successful
Definition: Control/AthToolSupport/AsgMessaging/AsgMessaging/MessageCheck.h:324
SG::ConstAccessor
Helper class to provide constant type-safe access to aux data.
Definition: ConstAccessor.h:54
JSSTaggerBase::m_decC2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decC2Key
Definition: JSSTaggerBase.h:146
JSSTaggerBase::m_decorationName
std::string m_decorationName
Decoration name.
Definition: JSSTaggerBase.h:202
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
JSSTaggerBase::m_readECFG331Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG331Key
Definition: JSSTaggerBase.h:181
LargeRJetTruthLabel::Zqq
@ Zqq
Definition: LargeRJetLabelEnum.h:20
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
isValid
bool isValid(const T &p)
Definition: AtlasPID.h:214
ASG_MAKE_ANA_TOOL
#define ASG_MAKE_ANA_TOOL(handle, type)
create the tool in the given tool handle
Definition: AnaToolHandle.h:690
yodamerge_tmp.scale
scale
Definition: yodamerge_tmp.py:138
asg::AnaToolHandle::setProperty
StatusCode setProperty(const std::string &property, const T2 &value)
set the given property of the tool.
Phi_mpi_pi
__HOSTDEV__ double Phi_mpi_pi(double)
Definition: GeoRegion.cxx:7
LargeRJetTruthLabel::qcd
@ qcd
Definition: LargeRJetLabelEnum.h:25
JSSTaggerBase::m_decValidPtRangeHighKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeHighKey
Definition: JSSTaggerBase.h:72
asg::AcceptInfo::getCutName
const std::string & getCutName(unsigned int cutPosition) const
Get the name of a cut, based on the cut position (slow, avoid usage)
Definition: AcceptInfo.cxx:36
JSSTaggerBase::m_jetPtMax
float m_jetPtMax
Definition: JSSTaggerBase.h:133
JSSTaggerBase::m_calibArea
std::string m_calibArea
Location where config files live on cvmfs.
Definition: JSSTaggerBase.h:114
JSSTaggerBase::m_decWeightKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decWeightKey
WriteDecorHandle keys for SF.
Definition: JSSTaggerBase.h:237
JSSTaggerBase::m_readTruthLabelKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTruthLabelKey
Truth label ReadDecorHandle key.
Definition: JSSTaggerBase.h:228
LargeRJetTruthLabel::notruth
@ notruth
Definition: LargeRJetLabelEnum.h:24
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
JSSTaggerBase::m_decSigeffSFKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decSigeffSFKey
Definition: JSSTaggerBase.h:240
dqt_zlumi_pandas.weight
int weight
Definition: dqt_zlumi_pandas.py:200
JSSTaggerBase::m_decEfficiencyKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decEfficiencyKey
Definition: JSSTaggerBase.h:238
AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:85
JSSTaggerBase::m_decE3Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decE3Key
Definition: JSSTaggerBase.h:148
JSSTaggerBase::m_ptGeV
bool m_ptGeV
Flag to indicate units pT is defined in Set to false by default.
Definition: JSSTaggerBase.h:90
efficiency
void efficiency(std::vector< double > &bins, std::vector< double > &values, const std::vector< std::string > &files, const std::string &histname, const std::string &tplotname, const std::string &label="")
Definition: dependence.cxx:128
JSSTaggerBase::m_decTau32WTAKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau32WTAKey
Definition: JSSTaggerBase.h:143
jet
Definition: JetCalibTools_PlotJESFactors.cxx:23
AthCommonDataStore
Definition: AthCommonDataStore.h:52
LargeRJetTruthLabel::Wqq
@ Wqq
Definition: LargeRJetLabelEnum.h:17
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
CheckAppliedSFs.e3
e3
Definition: CheckAppliedSFs.py:264
SG::ReadDecorHandle
Handle class for reading a decoration on an object.
Definition: StoreGate/StoreGate/ReadDecorHandle.h:94
JSSTaggerBase::getSF
std::pair< double, double > getSF(const xAOD::Jet &jet, const std::string &truthLabelStr) const
Get scale factor and efficiency.
Definition: JSSTaggerBase.cxx:670
lumiFormat.i
int i
Definition: lumiFormat.py:92
Trk::L3
@ L3
Definition: AlignModuleList.h:32
JSSTaggerBase::m_readECFG212Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG212Key
Definition: JSSTaggerBase.h:183
JSSTaggerBase::m_tmvaConfigFileName
std::string m_tmvaConfigFileName
TMVA configurations for BDT taggers.
Definition: JSSTaggerBase.h:123
asg::AcceptInfo::getNCuts
unsigned int getNCuts() const
Get the number of cuts defined.
Definition: AcceptInfo.h:46
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
JSSTaggerBase::checkKinRange
StatusCode checkKinRange(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Check and record if jet passes kinematic constraints.
Definition: JSSTaggerBase.cxx:370
CP::JetQGTaggerBDT::getJetProperties
bool getJetProperties(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Update the jet substructure variables for each jet to use in BDT.
Definition: JetQGTaggerBDT.cxx:166
SG::WriteDecorHandle
Handle class for adding a decoration to an object.
Definition: StoreGate/StoreGate/WriteDecorHandle.h:99
xAOD::VxType::PriVtx
@ PriVtx
Primary vertex.
Definition: TrackingPrimitives.h:571
JSSTaggerBase::Unknown
@ Unknown
Definition: JSSTaggerBase.h:99
BindingsTest.cut
cut
This script demonstrates how to call a C++ class from Python Also how to use PyROOT is shown.
Definition: BindingsTest.py:13
JSSTaggerBase::m_readL3Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readL3Key
Definition: JSSTaggerBase.h:179
CP::JetQGTaggerBDT::getScore
float getScore(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Retrieve BDT score.
Definition: JetQGTaggerBDT.cxx:150
JSSTaggerBase::m_efficiencyHistogramName
std::string m_efficiencyHistogramName
Definition: JSSTaggerBase.h:212
JSSTaggerBase::m_strScoreCut
std::string m_strScoreCut
Definition: JSSTaggerBase.h:194
JSSTaggerBase::m_decTau42WTAKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau42WTAKey
Definition: JSSTaggerBase.h:144
python.xAODType.dummy
dummy
Definition: xAODType.py:4
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
CP::JetQGTaggerBDT::m_trkSelectionTool
asg::AnaToolHandle< InDet::IInDetTrackSelectionTool > m_trkSelectionTool
Definition: JetQGTaggerBDT.h:64
MSG::name
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition: MsgLevel.cxx:19
JSSTaggerBase::m_decEffSFKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decEffSFKey
Definition: JSSTaggerBase.h:239
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
JSSTaggerBase::m_decNtrk500Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decNtrk500Key
Definition: JSSTaggerBase.h:186
JSSTaggerBase::m_calcSF
bool m_calcSF
Flag to calculate scale factor.
Definition: JSSTaggerBase.h:205
DataVector
Derived DataVector<T>.
Definition: DataVector.h:581
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
JSSTaggerBase::m_tmvaConfigFilePath
std::string m_tmvaConfigFilePath
Definition: JSSTaggerBase.h:124
JSSTaggerBase::m_decD2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decD2Key
Definition: JSSTaggerBase.h:147
CP::JetQGTaggerBDT::m_mode
int m_mode
Definition: JetQGTaggerBDT.h:66
LargeRJetTruthLabel::Wqq_From_t
@ Wqq_From_t
Definition: LargeRJetLabelEnum.h:21
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:192
JSSTaggerBase::initialize
virtual StatusCode initialize() override
Initialize the tool.
Definition: JSSTaggerBase.cxx:73
JSSTaggerBase::m_acceptInfo
asg::AcceptInfo m_acceptInfo
Object that stores the results for a jet.
Definition: JSSTaggerBase.h:65
ActsTrk::to_string
std::string to_string(const DetectorType &type)
Definition: GeometryDefs.h:34
JSSTaggerBase::m_decTau21WTAKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau21WTAKey
WriteDecorHandle keys for JSS moments.
Definition: JSSTaggerBase.h:142
InDet::IInDetTrackSelectionTool::accept
virtual asg::AcceptData accept(const xAOD::IParticle *p) const =0
Get the decision using a generic IParticle pointer.
PathResolverFindCalibFile
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
Definition: PathResolver.cxx:431
JSSTaggerBase::m_isMC
bool m_isMC
Definition: JSSTaggerBase.h:206
JSSTaggerBase::WBoson
@ WBoson
Definition: JSSTaggerBase.h:99
xAOD::Jet_v1
Class describing a jet.
Definition: Jet_v1.h:57
JSSTaggerBase::m_decValidKinRangeKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidKinRangeKey
Definition: JSSTaggerBase.h:75
JSSTaggerBase::m_readTau1WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau1WTAKey
ReadDecorHandle keys for JSS moments.
Definition: JSSTaggerBase.h:153
JSSTaggerBase::m_weightHistograms
std::map< std::string, std::unique_ptr< TH2D > > m_weightHistograms
Definition: JSSTaggerBase.h:217
DiTauMassTools::MaxHistStrategyV2::e
e
Definition: PhysicsAnalysis/TauID/DiTauMassTools/DiTauMassTools/HelperFunctions.h:26
a
TList * a
Definition: liststreamerinfos.cxx:10
asg::AcceptData::setCutResult
void setCutResult(const std::string &cutName, bool cutResult)
Set the result of a cut, based on the cut name (safer)
Definition: AcceptData.h:134
h
xAOD::Vertex_v1
Class describing a Vertex.
Definition: Vertex_v1.h:42
IJetDecorator::decorate
virtual StatusCode decorate(const xAOD::JetContainer &jets) const =0
Decorate a jet collection without otherwise modifying it.
JSSTaggerBase::passKinRange
bool passKinRange(const xAOD::Jet &jet) const
Check if jet passes kinematic constraints.
Definition: JSSTaggerBase.cxx:356
asg::AcceptData::clear
void clear()
Clear all bits.
Definition: AcceptData.h:54
CP::JetQGTaggerBDT::getPrecomputedVariables
bool getPrecomputedVariables(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Definition: JetQGTaggerBDT.cxx:195
asg::AcceptData::getCutResult
bool getCutResult(const std::string &cutName) const
Get the result of a cut, based on the cut name (safer)
Definition: AcceptData.h:98
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
JSSTaggerBase::m_decValidPtRangeLowKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeLowKey
Definition: JSSTaggerBase.h:73
JSSTaggerBase::m_readTau3WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau3WTAKey
Definition: JSSTaggerBase.h:155
JSSTaggerBase::m_readECF2Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF2Key
Definition: JSSTaggerBase.h:163
JSSTaggerBase::m_readECF3Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF3Key
Definition: JSSTaggerBase.h:164
AthCommonMsg< AlgTool >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
JSSTaggerBase::getConfigReader
StatusCode getConfigReader()
Get configReader StatusCode.
Definition: JSSTaggerBase.cxx:300
LargeRJetTruthLabel::intToEnum
TypeEnum intToEnum(const int type)
Definition: LargeRJetLabelEnum.h:57
xAOD::JetAttribute::TrackWidthPt500
@ TrackWidthPt500
Definition: JetAttributes.h:110
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
InDet::InDetTrackSelectionTool
Implementation of the track selector tool.
Definition: InDetTrackSelectionTool.h:51
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:616
JSSTaggerBase::JSSTaggerBase
JSSTaggerBase(const std::string &name)
Default constructor - to be used in all derived classes.
Definition: JSSTaggerBase.cxx:11
declareProperty
#define declareProperty(n, p, h)
Definition: BaseFakeBkgTool.cxx:15
CP::JetQGTaggerBDT::calculateVariables
bool calculateVariables(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Definition: JetQGTaggerBDT.cxx:230
defineDB.jets
list jets
Definition: JetTagCalibration/share/defineDB.py:24
dqt_zlumi_alleff_HIST.eff
int eff
Definition: dqt_zlumi_alleff_HIST.py:113
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:790
LargeRJetTruthLabel::tqqb
@ tqqb
Definition: LargeRJetLabelEnum.h:16
xAOD::TrackParticle_v1
Class describing a TrackParticle.
Definition: TrackParticle_v1.h:43
JSSTaggerBase::m_decL2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL2Key
Definition: JSSTaggerBase.h:149
JSSTaggerBase::m_weightFlavors
std::string m_weightFlavors
Definition: JSSTaggerBase.h:213
MuonParameters::beta
@ beta
Definition: MuonParamDefs.h:144
CP::JetQGTaggerBDT::isCorrectNumberOfTracks
bool isCorrectNumberOfTracks(int expectedNTracks, int nTracksFromGhostTracks) const
Definition: JetQGTaggerBDT.cxx:350
asg::AcceptData
Definition: AcceptData.h:30
xAOD::JetAttribute::NumTrkPt500
@ NumTrkPt500
Definition: JetAttributes.h:106
SG::DataProxy
Definition: DataProxy.h:44
JSSTaggerBase::m_readL2Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readL2Key
Definition: JSSTaggerBase.h:178
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
readCCLHist.float
float
Definition: readCCLHist.py:83
DataVector::empty
bool empty() const noexcept
Returns true if the collection is empty.
fitman.k
k
Definition: fitman.py:528
asg::AcceptInfo::getCutDescription
const std::string & getCutDescription(const std::string &cutName) const
Get the description of a cut, based on the cut name.
Definition: AcceptInfo.cxx:119
JSSTaggerBase::m_decTaggedKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTaggedKey
WriteDecorHandle keys for tagging bools.
Definition: JSSTaggerBase.h:71
JSSTaggerBase::m_containerName
std::string m_containerName
Configurable members.
Definition: JSSTaggerBase.h:105