ATLAS Offline Software
Public Member Functions | Protected Types | Protected Member Functions | Protected Attributes | Private Types | Private Member Functions | Private Attributes | List of all members
JSSTaggerBase Class Referenceabstract

#include <JSSTaggerBase.h>

Inheritance diagram for JSSTaggerBase:
Collaboration diagram for JSSTaggerBase:

Public Member Functions

virtual StatusCode initialize () override
 Initialize the tool. 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, V, H > &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

 JSSTaggerBase (const std::string &name)
 Default constructor - to be used in all derived classes. More...
 
 ~JSSTaggerBase ()
 Default destructor - to be used in all derived classes. More...
 
virtual StatusCode tag (const xAOD::Jet &jet) const =0
 Decorate single jet with tagging info. More...
 
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...
 
void decorateJSSRatios (const xAOD::JetContainer &jets) const
 
int findPV () const
 Find the PV (to be used for Ntrk) More...
 
int 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::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"}
 
SG::ReadDecorHandleKey< xAOD::JetContainerm_readNtrk500Key {this, "Ntrk500Name", "ParentJetNTrkPt500", "SG key for Ntrk500 from ungroomed jet"}
 
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 {Unknown}
 
std::string m_containerName
 Configurable members. More...
 
bool m_isSmallRJet = false
 
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_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_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_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"}
 
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

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

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 37 of file JSSTaggerBase.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
protected

TAGTYPE enum.

Enumerator
Unknown 
WBoson 
ZBoson 
TopQuark 

Definition at line 98 of file JSSTaggerBase.h.

Constructor & Destructor Documentation

◆ JSSTaggerBase()

JSSTaggerBase::JSSTaggerBase ( const std::string &  name)
protected

Default constructor - to be used in all derived classes.

Tagger configuration properties

Jet kinematics properties

Truth labeling properties

Keras properties

TMVA properties

Tagging scale factors

Definition at line 11 of file JSSTaggerBase.cxx.

11  :
13  m_calibArea(""),
14  m_jetPtMin(200000.),
15  m_jetPtMax(3000000.),
16  m_jetEtaMax(2.0),
17  m_strMassCutLow(""),
18  m_strMassCutHigh(""),
19  m_strScoreCut("")
20 {
21 
23  declareProperty( "ContainerName", m_containerName = "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets", "Name of jet container" );
24  declareProperty( "ConfigFile", m_configFile = "", "Name of config file" );
25  declareProperty( "CalibArea", m_calibArea = "", "Path to config file" );
26  declareProperty( "CalcSF", m_calcSF = false, "Flag to calculate the efficiency SF" );
27  declareProperty( "WorkingPoint", m_wkpt = "", "If specified, name of working point is added to the variable names" );
28  declareProperty( "TaggerType", m_tagType = "XXX", "Tagger type (e.g. SmoothedWZTagger, JSSWTopTaggerDNN, etc.)" );
29  declareProperty( "Decoration", m_decorationName = "XX", "Prefix for the variables decorated to xAOD::Jet" );
30  declareProperty( "IsMC", m_isMC = true, "Flag to identify data or MC" );
31 
33  declareProperty( "JetEtaMax", m_jetEtaMax = 2.0, "Eta cut to define fiducial phase space for the tagger");
34 
36  declareProperty( "UseTRUTH3", m_truthLabelUseTRUTH3 = true,
37  "Flag to use TRUTH3 containers. If false, TRUTH1 format is used." );
38  declareProperty( "TruthParticleContainerName", m_truthParticleContainerName = "TruthParticles",
39  "Name of truth-particle container (with UseTRUTH3=false). TruthParticles by default" );
40  declareProperty( "TruthBosonContainerName", m_truthBosonContainerName = "TruthBosonsWithDecayParticles",
41  "Name of truth-boson container (with UseTRUTH3=true). TruthBosonWithDecayParticles by default" );
42  declareProperty( "TruthTopQuarkContainerName", m_truthTopQuarkContainerName = "TruthTopQuarkWithDecayParticles",
43  "Name of truth-top container (with UseTRUTH3=true). TruthTopQuarkWithDecayParticles by default" );
44 
46  declareProperty( "CalibAreaKeras", m_kerasCalibArea = "BoostedJetTaggers/TopoclusterTopTagger/Boost2017/",
47  "Path to json file to configure ML-taggers (Keras)" );
48  declareProperty( "KerasConfigFile", m_kerasConfigFileName = "XXX",
49  "Name of json file to configure ML-taggers (Keras)" );
50  declareProperty( "KerasOutput", m_kerasConfigOutputName = "XXX",
51  "Name of output variable by the ML-tagger (Keras)" );
52 
54  declareProperty( "CalibAreaTMVA", m_tmvaCalibArea = "BoostedJetTaggers/JSSWTopTaggerBDT/Boost2017/",
55  "Path to xml file to configure ML-taggers (TMVA)" );
56  declareProperty( "TMVAConfigFile", m_tmvaConfigFileName = "XXX",
57  "Name of xml file to configure ML-taggers (TMVA)" );
58 
60  declareProperty( "WeightDecorationName", m_weightDecorationName = "SF",
61  "Name of SF variable decorated to xAOD::Jet" );
62  declareProperty( "WeightFile", m_weightFileName = "",
63  "Name of config ROOT file for SF calculation" );
64  declareProperty( "WeightHistogramName", m_weightHistogramName = "",
65  "Name of SF histograms in the ROOT file" );
66  declareProperty( "EfficiencyHistogramName", m_efficiencyHistogramName = "",
67  "Name of efficiency histograms in the ROOT file" );
68  declareProperty( "WeightFlavors", m_weightFlavors = "",
69  "List of jet flavours for which the SF is available. Divided by comma" );
70  declareProperty( "SuppressOutputDependence", m_suppressOutputDependence = false );
71 }

◆ ~JSSTaggerBase()

JSSTaggerBase::~JSSTaggerBase ( )
inlineprotected

Default destructor - to be used in all derived classes.

Definition at line 56 of file JSSTaggerBase.h.

56 {};

Member Function Documentation

◆ calculateJSSRatios()

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

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 398 of file JSSTaggerBase.cxx.

398  {
399 
400  int result = 0;
401 
411 
417 
419  float tau21_wta = -999.0;
420  float tau32_wta = -999.0;
421  float tau42_wta = -999.0;
422 
423  float tau1_wta = readTau1WTA(jet);
424  float tau2_wta = readTau2WTA(jet);
425  float tau3_wta = readTau3WTA(jet);
426  float tau4_wta = -999.0;
427  if(readTau4WTA.isAvailable()){
428  tau4_wta = readTau4WTA(jet);
429  }
430 
431  if ( tau1_wta > 1e-8 ) {
432  tau21_wta = tau2_wta / tau1_wta;
433  }
434  else result = 1;
435 
436  if ( tau2_wta > 1e-8 ) {
437  tau32_wta = tau3_wta / tau2_wta;
438  if(readTau4WTA.isAvailable()){
439  tau42_wta = tau4_wta / tau2_wta;
440  }
441  }
442  else result = 1;
443 
444  decTau21WTA(jet) = tau21_wta;
445  decTau32WTA(jet) = tau32_wta;
446  decTau42WTA(jet) = tau42_wta;
447 
449  float C2 = -999.0;
450  float D2 = -999.0;
451  float e3 = -999.0;
452 
453  // L-series for UFO top taggers
454  float L2 = -999.0;
455  float L3 = -999.0;
456 
457  if(!m_isSmallRJet){
461 
462  float ECF1 = readECF1(jet);
463  float ECF2 = readECF2(jet);
464  float ECF3 = readECF3(jet);
465 
466  if ( ECF2 > 1e-8 ) {
467  C2 = ECF3 * ECF1 / std::pow( ECF2, 2.0 );
468  D2 = ECF3 * std::pow( ECF1, 3.0 ) / std::pow( ECF2, 3.0 );
469  }
470  else result = 1;
471 
472  e3 = ECF3 / std::pow( ECF1, 3.0 );
473 
477 
478  if(readECFG331.isAvailable() && readECFG212.isAvailable()){
479  if(readECFG212(jet) > 1e-8){
480  L2 = readECFG331(jet) / std::pow(readECFG212(jet), 1.5);
481  }
482  else result = 1;
483  }
484 
485  if(readECFG331.isAvailable() && readECFG311.isAvailable()){
486  if(readECFG331(jet) > 1e-8){
487  L3 = readECFG311(jet) / std::pow(readECFG331(jet), 1./3.);
488  }
489  else result = 1;
490  }
491  }
492 
493  decC2(jet) = C2;
494  decD2(jet) = D2;
495  decE3(jet) = e3;
496 
497  static const SG::AuxElement::ConstAccessor<float> accL2("L2");
498  if(!accL2.isAvailable(jet)) decL2(jet) = L2;
499 
500  static const SG::AuxElement::ConstAccessor<float> accL3("L3");
501  if(!accL3.isAvailable(jet)) decL3(jet) = L3;
502 
503  // TODO: Add ECFG for ANN tagger whenever it is defined
504 
505  return result;
506 
507 }

◆ checkKinRange()

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

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 352 of file JSSTaggerBase.cxx.

352  {
353 
354  float scale = 1.0;
355  if ( m_ptGeV ) scale = 1.e3;
356 
357  bool passKinRange = true;
358 
361  if ( std::abs(jet.eta()) > m_jetEtaMax ) {
362  ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (|eta| < " << m_jetEtaMax << "). Jet eta = " << jet.eta() );
363  acceptData.setCutResult( "ValidEtaRange", false );
364  passKinRange = false;
365  }
366 
367  if ( jet.pt() < m_jetPtMin * scale ) {
368  ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (pT > " << m_jetPtMin * scale / 1.e3 << "). Jet pT = " << jet.pt() / 1.e3 << " GeV" );
369  acceptData.setCutResult( "ValidPtRangeLow", false );
370  passKinRange = false;
371  }
372 
373  if ( jet.pt() > m_jetPtMax * scale ) {
374  ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (pT < " << m_jetPtMax * scale / 1.e3 << "). Jet pT = " << jet.pt() / 1.e3 << " GeV" );
375  acceptData.setCutResult( "ValidPtRangeHigh", false );
376  passKinRange = false;
377  }
378 
384 
386  decValidPtRangeHigh(jet) = acceptData.getCutResult( "ValidPtRangeHigh" );
387  decValidPtRangeLow(jet) = acceptData.getCutResult( "ValidPtRangeLow" );
388  decValidEtaRange(jet) = acceptData.getCutResult( "ValidEtaRange" );
389  decValidKinRange(jet) = passKinRange;
390 
391  return StatusCode::SUCCESS;
392 
393 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T, V, H > &  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, V, H > &  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, V, H > &  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, V, H > &  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, V, H > &  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
overridevirtual

Decorate jet collection with tagging info.

Loop over jet collection and decorate each jet.

Implements IJetDecorator.

Reimplemented in JSSWTopTaggerANN, and JSSWTopTaggerDNN.

Definition at line 270 of file JSSTaggerBase.cxx.

270  {
271 
272  for ( const auto *jet : jets ) {
273  ATH_CHECK( tag(*jet) );
274  }
275 
276  return StatusCode::SUCCESS;
277 
278 }

◆ decorateJSSRatios()

void JSSTaggerBase::decorateJSSRatios ( const xAOD::JetContainer jets) const
protected

Create write decor handles

Create read decor handles

WTA N-subjettiness ratios

ECF ratios

Definition at line 510 of file JSSTaggerBase.cxx.

510  {
511 
514  bool decor_available = decTau21WTA.isAvailable();
515  if(decor_available) return;
516 
524 
530 
531  // Use pointers here so we can create only the ones we're configured for
532  std::unique_ptr<SG::ReadDecorHandle<xAOD::JetContainer, float>> readECF1;
533  std::unique_ptr<SG::ReadDecorHandle<xAOD::JetContainer, float>> readECF2;
534  std::unique_ptr<SG::ReadDecorHandle<xAOD::JetContainer, float>> readECF3;
535 
536  std::unique_ptr<SG::ReadDecorHandle<xAOD::JetContainer, float>> readECFG331;
537  std::unique_ptr<SG::ReadDecorHandle<xAOD::JetContainer, float>> readECFG311;
538  std::unique_ptr<SG::ReadDecorHandle<xAOD::JetContainer, float>> readECFG212;
539 
540  if(!m_isSmallRJet){
541  readECF1 = std::make_unique<SG::ReadDecorHandle<xAOD::JetContainer, float> >(m_readECF1Key);
542  readECF2 = std::make_unique<SG::ReadDecorHandle<xAOD::JetContainer, float> >(m_readECF2Key);
543  readECF3 = std::make_unique<SG::ReadDecorHandle<xAOD::JetContainer, float> >(m_readECF3Key);
544 
545  readECFG331 = std::make_unique<SG::ReadDecorHandle<xAOD::JetContainer, float> >(m_readECFG331Key);
546  readECFG311 = std::make_unique<SG::ReadDecorHandle<xAOD::JetContainer, float> >(m_readECFG311Key);
547  readECFG212 = std::make_unique<SG::ReadDecorHandle<xAOD::JetContainer, float> >(m_readECFG212Key);
548  }
549 
550  for(const xAOD::Jet* jet : jets){
551 
553  float tau21_wta = -999.0;
554  float tau32_wta = -999.0;
555  float tau42_wta = -999.0;
556 
557  float tau1_wta = readTau1WTA(*jet);
558  float tau2_wta = readTau2WTA(*jet);
559  float tau3_wta = readTau3WTA(*jet);
560  float tau4_wta = -999.0;
561  if(readTau4WTA.isAvailable()){
562  tau4_wta = readTau4WTA(*jet);
563  }
564 
565  if ( tau1_wta > 1e-8 ) {
566  tau21_wta = tau2_wta / tau1_wta;
567  }
568 
569  if ( tau2_wta > 1e-8 ) {
570  tau32_wta = tau3_wta / tau2_wta;
571  if(readTau4WTA.isAvailable()){
572  tau42_wta = tau4_wta / tau2_wta;
573  }
574  }
575 
576  decTau21WTA(*jet) = tau21_wta;
577  decTau32WTA(*jet) = tau32_wta;
578  decTau42WTA(*jet) = tau42_wta;
579 
580  float C2 = -999.0;
581  float D2 = -999.0;
582  float e3 = -999.0;
583 
584  float L2 = -999.0;
585  float L3 = -999.0;
586 
587  if(!m_isSmallRJet){
589  float ECF1 = (*readECF1)(*jet);
590  float ECF2 = (*readECF2)(*jet);
591  float ECF3 = (*readECF3)(*jet);
592 
593  if ( ECF2 > 1e-8 ) {
594  C2 = ECF3 * ECF1 / std::pow( ECF2, 2.0 );
595  D2 = ECF3 * std::pow( ECF1, 3.0 ) / std::pow( ECF2, 3.0 );
596  }
597 
598  e3 = ECF3 / std::pow( ECF1, 3.0 );
599 
600  // L-series for UFO top taggers
601  if((*readECFG331).isAvailable() && (*readECFG212).isAvailable()){
602  if((*readECFG212)(*jet) > 1e-8){
603  L2 = (*readECFG331)(*jet) / std::pow((*readECFG212)(*jet), 1.5);
604  }
605  }
606 
607  if((*readECFG331).isAvailable() && (*readECFG311).isAvailable()){
608  if((*readECFG331)(*jet) > 1e-8){
609  L3 = (*readECFG311)(*jet) / std::pow((*readECFG331)(*jet), 1./3.);
610  }
611  }
612  }
613 
614  decC2(*jet) = C2;
615  decD2(*jet) = D2;
616  decE3(*jet) = e3;
617 
618  static const SG::AuxElement::ConstAccessor<float> accL2("L2");
619  if(!accL2.isAvailable(*jet)) decL2(*jet) = L2;
620 
621  static const SG::AuxElement::ConstAccessor<float> accL3("L3");
622  if(!accL3.isAvailable(*jet)) decL3(*jet) = L3;
623 
624  // TODO: Add ECFG for ANN tagger whenever it is defined
625 
626  }
627 
628  return;
629 
630 }

◆ 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
protected

Find the PV (to be used for Ntrk)

Definition at line 633 of file JSSTaggerBase.cxx.

633  {
634 
635  int indexPV = -1;
636 
637  const xAOD::VertexContainer* vxCont = nullptr;
638  if ( evtStore()->retrieve( vxCont, "PrimaryVertices" ).isFailure() ) {
639  ATH_MSG_WARNING( "Unable to retrieve primary vertex container PrimaryVertices" );
640  }
641  else {
642  int vtx_counter = 0;
643  for ( const auto *vx : *vxCont ) {
644  if ( vx->vertexType()==xAOD::VxType::PriVtx ) {
645  indexPV = vtx_counter;
646  break;
647  }
648  vtx_counter++;
649  }
650  }
651 
652  return indexPV;
653 
654 }

◆ getConfigReader()

StatusCode JSSTaggerBase::getConfigReader ( )
protected

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 281 of file JSSTaggerBase.cxx.

281  {
282 
283  ATH_MSG_INFO( "Using config file : " << m_configFile );
284 
286  std::string configPath;
287 
288  if ( m_calibArea.compare("Local") == 0 ) {
290  }
291  else if ( m_calibArea.find("eos") != std::string::npos) {
292  configPath = PathResolverFindCalibFile((m_calibArea+"/"+m_configFile).c_str());
293  configPath = (m_calibArea+"/"+m_configFile).c_str();
294  }
295  else {
296  configPath = PathResolverFindCalibFile(("BoostedJetTaggers/"+m_calibArea+"/"+m_configFile).c_str());
297  }
298 
300  FileStat_t fStats;
301  int fSuccess = gSystem->GetPathInfo(configPath.c_str(), fStats);
302  if ( fSuccess ) {
303  ATH_MSG_ERROR( "Recommendations file " << m_configFile << " could not be found" );
304  return StatusCode::FAILURE;
305  }
306  else {
307  ATH_MSG_DEBUG( "Recommendations file was found : " << configPath );
308  }
309 
310  if ( m_configReader.ReadFile( configPath.c_str(), EEnvLevel(0) ) ) {
311  ATH_MSG_ERROR( "Error while reading config file : "<< configPath );
312  return StatusCode::FAILURE;
313  }
314 
315  return StatusCode::SUCCESS;
316 
317 }

◆ 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  }

◆ getProperty()

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

Get one of the tool's properties.

◆ getSF()

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

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 769 of file JSSTaggerBase.cxx.

769  {
770 
771  if ( !passKinRange(jet) ) return std::make_pair( 1.0, 1.0 );
772 
773 
774  double logmOverPt = std::log(jet.m()/jet.pt());
775  if ( m_decorationName.find("SmoothZ") != std::string::npos ||
776  m_decorationName.find("SmoothInclusiveZ") != std::string::npos ) {
778  const double WtoZmassShift = 10803;
779  logmOverPt = std::log((jet.m()-WtoZmassShift)/jet.pt());
780  }
781 
782  if ( logmOverPt > 0 ) logmOverPt = 0;
783 
784  double SF = 1.0;
785  double eff = 1.0;
786 
787  if ( m_weightHistograms.count(truthLabelStr.c_str()) ) {
788 
789  int pt_mPt_bin = (m_weightHistograms.find(truthLabelStr.c_str())->second)->FindBin(jet.pt()*0.001, logmOverPt);
790  SF = (m_weightHistograms.find(truthLabelStr.c_str())->second)->GetBinContent(pt_mPt_bin);
791 
792  if ( !m_efficiencyHistogramName.empty() ) {
793  eff = (m_efficiencyHistograms.find(truthLabelStr.c_str())->second)->GetBinContent(pt_mPt_bin);
794  }
795 
796  }
797  else {
798  // set the efficiency for "Other" category to be the signal efficiency
799  std::string signal_truthLabel="";
800  if ( m_weightHistograms.count("t_qqb") ) {
801  signal_truthLabel="t_qqb";
802  }else if ( m_weightHistograms.count("V_qq") ){
803  signal_truthLabel="V_qq";
804  }else if ( m_weightHistograms.count("t") ){
805  signal_truthLabel="t";
806  }
807  if ( !signal_truthLabel.empty() && !m_efficiencyHistogramName.empty() ){
808  int pt_mPt_bin = (m_weightHistograms.find(signal_truthLabel.c_str())->second)->FindBin(jet.pt()*0.001, logmOverPt);
809  eff = (m_efficiencyHistograms.find(signal_truthLabel.c_str())->second)->GetBinContent(pt_mPt_bin);
810  }
811 
812  return std::make_pair( 1.0, eff );
813  }
814 
815  if ( SF < 1e-3 ) {
816  ATH_MSG_DEBUG( "(pt, m/pt) (" << jet.pt()/1.e3 << ", " << jet.m()/jet.pt() << ") is out of range for SF calculation. Returning 1.0" );
817  return std::make_pair( 1.0, 1.0 );
818  }
819  else {
820  return std::make_pair( SF, eff );
821  }
822 
823 }

◆ getTruthLabelStr()

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

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 827 of file JSSTaggerBase.cxx.

827  {
829  std::string truthLabelStr;
830 
833  LargeRJetTruthLabel::TypeEnum jetContainment = LargeRJetTruthLabel::intToEnum(readTruthLabel(jet));
834 
836  if ( m_weightHistograms.count("t_qqb") ) {
837 
839  if ( jetContainment==LargeRJetTruthLabel::tqqb ) {
840  truthLabelStr = "t_qqb";
841  }
843  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
844  truthLabelStr = "q";
845  }
846 
847  }
849  else if ( m_weightHistograms.count("V_qq_passMpassD2") ) {
850 
852  if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t ) {
853  truthLabelStr = "t_";
854  }
856  else if ( jetContainment==LargeRJetTruthLabel::Wqq || jetContainment==LargeRJetTruthLabel::Zqq || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
857  truthLabelStr = "V_qq_";
858  }
860  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
861  truthLabelStr = "q_";
862  }
863 
865  if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") && acceptData.getCutResult("PassD2") ) {
866  truthLabelStr += "passMpassD2";
867  }
869  else if ( !(acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh")) && acceptData.getCutResult("PassD2") ) {
870  truthLabelStr += "failMpassD2";
871  }
873  else if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") && !acceptData.getCutResult("PassD2") ) {
874  truthLabelStr += "passMfailD2";
875  }
877  else{
878  truthLabelStr += "failMfailD2";
879  }
880 
881  }
882 
884  else if ( m_weightHistograms.count("V_qq") ) {
885 
887  if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t ) {
888  truthLabelStr = "t";
889  }
891  else if ( jetContainment==LargeRJetTruthLabel::Wqq || jetContainment==LargeRJetTruthLabel::Zqq || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
892  truthLabelStr = "V_qq";
893  }
895  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
896  truthLabelStr = "q";
897  }
898  }
899 
900  // inclusive top tagger
901  else {
903  if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
904  truthLabelStr = "t";
905  }
907  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
908  truthLabelStr = "q";
909  }
910 
911  }
912 
913  return truthLabelStr;
914 }

◆ GetUnGroomTracks()

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

Retrieve Ntrk variable from the ungroomed parent jet.

Retrieve the Ntrk variable from the ungroomed parent jet.

Definition at line 657 of file JSSTaggerBase.cxx.

657  {
658 
660 
661  const xAOD::Jet * ungroomedJet = nullptr;
662  if ( readParent.isAvailable() ) {
663  const ElementLink<xAOD::JetContainer>& linkToUngroomed = readParent(jet);
664  if ( linkToUngroomed.isValid() ) {
665  ungroomedJet = *linkToUngroomed;
666 
667  static const SG::AuxElement::ConstAccessor< std::vector<int> >acc_Ntrk("NumTrkPt500");
668 
669  if ( acc_Ntrk.isAvailable(*ungroomedJet) ) {
670 
671  const std::vector<int>& NTrkPt500 = acc_Ntrk(*ungroomedJet);
672 
673  int jet_ntrk = NTrkPt500.at(indexPV);
674  return jet_ntrk;
675 
676  }
677  else {
678  ATH_MSG_ERROR("WARNING: Unable to retrieve Ntrk of the ungroomed parent jet. Please make sure this variable is in your derivations!!!");
679  }
680  }
681  else {
682  ATH_MSG_ERROR("WARNING: Unable to retrieve the parent ungroomed jet. Please make sure this variable is in your derivations!!!");
683  }
684  }
685  else {
686  ATH_MSG_ERROR("WARNING: Unable to retrieve the link to the parent ungroomed jet. Please make sure this variable is in your derivations!!!");
687  }
688 
689  return -999;
690 
691 }

◆ getWeight()

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

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 695 of file JSSTaggerBase.cxx.

695  {
696  if ( !m_calcSF ) return StatusCode::SUCCESS;
697 
698  float weight = 1.0;
699  float effSF = 1.0;
700  float sigeffSF = 1.0;
701  float efficiency = 1.0;
702 
703  if ( m_isMC ) {
704 
705  std::string truthLabelStr = getTruthLabelStr( jet, acceptData );
706  std::tie(effSF, efficiency) = getSF( jet, truthLabelStr );
707 
708  // calculate signal efficiency SF
709  if ( m_weightHistograms.count("t_qqb") ) {
710  sigeffSF = getSF(jet, "t_qqb").first;
711  } else if ( m_weightHistograms.count("V_qq") ) {
712  sigeffSF = getSF(jet, "V_qq").first;
713  } else if ( m_weightHistograms.count("t") ){
714  sigeffSF = getSF(jet, "t").first;
715  } else {
716  sigeffSF = 1.0;
717  }
718 
720  if ( m_weightFlavors.find("fail") != std::string::npos ) {
721  weight = effSF;
722  }
723 
724  else {
725 
727  if ( passSel ) {
728  weight = effSF;
729  }
730 
732  else {
734  if ( m_efficiencyHistogramName.empty() ) {
735  weight = 1.0;
736  }
737  else if ( efficiency < 1.0 ) {
738  weight = ( 1. - effSF * efficiency ) / ( 1. - efficiency );
739  }
740  else {
741  weight = 1.0;
742  }
743  }
744  }
745 
746  }
747 
748  else {
749  weight = 1.0;
750  }
751 
757 
759  decWeight(jet) = weight;
760  decEfficiency(jet) = efficiency;
761  decEffSF(jet) = effSF;
762  decSigeffSF(jet) = sigeffSF;
763 
764  return StatusCode::SUCCESS;
765 
766 }

◆ initialize()

StatusCode JSSTaggerBase::initialize ( )
overridevirtual

Initialize the tool.

Make sure jet container is set

Initialize warning counters

Define common tagger states

Initialize decorators

Initialize SFs if they are needed

Get weight config file

Install histograms for tagging SFs

Initialize the TF1 functions

Reimplemented from asg::AsgTool.

Reimplemented in CP::JetQGTagger, CP::JetQGTaggerBDT, JSSWTopTaggerANN, JSSTaggerUtils, JSSWTopTaggerDNN, SmoothedTopTagger, and SmoothedWZTagger.

Definition at line 73 of file JSSTaggerBase.cxx.

73  {
74 
76  if ( m_containerName.empty() ) {
77  ATH_MSG_ERROR( "ContainerName has not been set. Exiting" );
78  return StatusCode::FAILURE;
79  }
80  m_isSmallRJet = m_containerName.find("AntiKt4")!=std::string::npos;
81 
83  m_nWarnVar = 0;
84 
86  m_acceptInfo.addCut( "ValidPtRangeHigh" , "True if the jet is not too high pT" );
87  m_acceptInfo.addCut( "ValidPtRangeLow" , "True if the jet is not too low pT" );
88  m_acceptInfo.addCut( "ValidEtaRange" , "True if the jet is not too forward" );
89 
90  m_acceptInfo.addCut( "ValidJetContent" , "True if the jet is alright technically (e.g. all attributes necessary for tag)" );
91  m_acceptInfo.addCut( "ValidEventContent" , "True if the event is alright technically (e.g. primary vertices)" );
92 
94  ATH_MSG_INFO( "Decorators that will be attached to jet :" );
95 
103 
112 
119  m_readQwKey = m_containerName + "." + m_readQwKey.key();
122  if(!m_isSmallRJet){
129  }
130 
139 
140  ATH_CHECK( m_readTau1WTAKey.initialize() );
141  ATH_CHECK( m_readTau2WTAKey.initialize() );
142  ATH_CHECK( m_readTau3WTAKey.initialize() );
143  ATH_CHECK( m_readTau4WTAKey.initialize() );
144  ATH_CHECK( m_readSplit12Key.initialize() );
145  ATH_CHECK( m_readSplit23Key.initialize() );
146  ATH_CHECK( m_readQwKey.initialize() );
147  ATH_CHECK( m_readThrustMajKey.initialize() );
148  ATH_CHECK( m_readSphericityKey.initialize() );
149  ATH_CHECK( m_readECF1Key.initialize(!m_isSmallRJet) );
150  ATH_CHECK( m_readECF2Key.initialize(!m_isSmallRJet) );
151  ATH_CHECK( m_readECF3Key.initialize(!m_isSmallRJet) );
152  ATH_CHECK( m_readECFG331Key.initialize(!m_isSmallRJet) );
153  ATH_CHECK( m_readECFG311Key.initialize(!m_isSmallRJet) );
154  ATH_CHECK( m_readECFG212Key.initialize(!m_isSmallRJet) );
155 
157  ATH_CHECK( m_readParentKey.initialize() );
158 
166 
167  ATH_MSG_INFO( " " << m_decTaggedKey.key() << " : pass tagging criteria" );
168  ATH_MSG_INFO( " " << m_decValidPtRangeHighKey.key() << " : pass upper pt range" );
169  ATH_MSG_INFO( " " << m_decValidPtRangeLowKey.key() << " : pass lower pt range" );
170  ATH_MSG_INFO( " " << m_decValidEtaRangeKey.key() << " : pass eta range" );
171  ATH_MSG_INFO( " " << m_decValidKinRangeKey.key() << " : pass kinematic range" );
172  ATH_MSG_INFO( " " << m_decValidJetContentKey.key() << " : has valid jet content" );
173  ATH_MSG_INFO( " " << m_decValidEventContentKey.key() << " : has valid event content" );
174 
175 
179 
183 
184  ATH_MSG_INFO( " " << m_decPassMassKey.key() << " : pass mass cut" );
185  ATH_MSG_INFO( " " << m_decCutMLowKey.key() << " : lower mass cut" );
186  ATH_MSG_INFO( " " << m_decCutMHighKey.key() << " : upper mass cut" );
187 
188 
192 
196 
197  ATH_MSG_INFO( " " << m_decPassScoreKey.key() << " : pass MVA score cut" );
198  ATH_MSG_INFO( " " << m_decScoreCutKey.key() << " : MVA score cut" );
199  ATH_MSG_INFO( " " << m_decScoreValueKey.key() << " : evaluated MVA score" );
200 
201 
206 
211 
212  ATH_MSG_INFO( " " << m_decWeightKey.key() << " : tagging SF" );
213 
215  ATH_CHECK( m_readTruthLabelKey.initialize() );
216 
218  ATH_CHECK( m_readNtrk500Key.initialize() );
219 
220 #ifndef XAOD_STANDALONE
234  }
235 #endif
236 
238  if ( m_calcSF ) {
239 
241  m_weightConfig = std::make_unique<TFile>( m_weightConfigPath.c_str() );
242  if( !m_weightConfig ) {
243  ATH_MSG_ERROR( "SmoothedWZTagger: Error openning config file : " << m_weightConfigPath );
244  return StatusCode::FAILURE;
245  }
246 
248  std::stringstream ss{m_weightFlavors};
249  std::string flavor;
250  while ( std::getline(ss, flavor, ',') ) {
251  m_weightHistograms.insert( std::make_pair( flavor, (TH2D*)m_weightConfig->Get((m_weightHistogramName+"_"+flavor).c_str()) ) );
252  if ( !m_efficiencyHistogramName.empty() ) {
253  m_efficiencyHistograms.insert( std::make_pair( flavor, (TH2D*)m_weightConfig->Get((m_efficiencyHistogramName+"_"+flavor).c_str()) ) );
254  }
255  ATH_MSG_INFO( "Tagging SF histogram for " << flavor << " is installed." );
256  }
257 
258  }
259 
261  if ( !m_strMassCutLow.empty() ) m_funcMassCutLow = std::make_unique<TF1>("strMassCutLow", m_strMassCutLow.c_str(), 0, 14000);
262  if ( !m_strMassCutHigh.empty() ) m_funcMassCutHigh = std::make_unique<TF1>("strMassCutHigh", m_strMassCutHigh.c_str(), 0, 14000);
263  if ( !m_strScoreCut.empty() ) m_funcScoreCut = std::make_unique<TF1>("strScoreCut", m_strScoreCut.c_str(), 0, 14000);
264 
265  return StatusCode::SUCCESS;
266 
267 }

◆ 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.

◆ 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
protected

Check if jet passes kinematic constraints.

Definition at line 338 of file JSSTaggerBase.cxx.

338  {
339 
340  float scale = 1.0;
341  if ( m_ptGeV ) scale = 1.e3;
342 
343  if ( jet.pt() < m_jetPtMin * scale ) return false;
344  if ( jet.pt() > m_jetPtMax * scale ) return false;
345  if ( std::abs( jet.eta() ) > m_jetEtaMax ) return false;
346 
347  return true;
348 
349 }

◆ print()

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

◆ printCuts()

void JSSTaggerBase::printCuts ( ) const
protected

Print configured cuts.

Definition at line 916 of file JSSTaggerBase.cxx.

916  {
917  ATH_MSG_INFO( "After tagging, you will have access to the following cuts as an asg::AcceptData : (<NCut>) <cut> : <description>)" );
918  int nCuts = m_acceptInfo.getNCuts();
919  for ( int iCut=0; iCut < nCuts; iCut++ ) {
920  std::string cut_string = "";
921  cut_string += " (";
922  cut_string += std::to_string(iCut);
923  cut_string += ") ";
924  cut_string += m_acceptInfo.getCutName(iCut).data();
925  cut_string += " : ";
926  cut_string += m_acceptInfo.getCutDescription(iCut).data();
927  ATH_MSG_INFO( cut_string );
928  }
929 }

◆ 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
protected

Reset cuts.

Reset the AcceptData cut results to false

Initialize common cuts to true by default

Definition at line 320 of file JSSTaggerBase.cxx.

320  {
321 
323  acceptData.clear();
324 
326  acceptData.setCutResult( "ValidJetContent", true );
327  acceptData.setCutResult( "ValidEventContent", true );
328 
329  acceptData.setCutResult( "ValidPtRangeHigh", true );
330  acceptData.setCutResult( "ValidPtRangeLow" , true );
331  acceptData.setCutResult( "ValidEtaRange" , true );
332 
333  return StatusCode::SUCCESS;
334 
335 }

◆ 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()

virtual StatusCode JSSTaggerBase::tag ( const xAOD::Jet jet) const
protectedpure virtual

◆ 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

◆ m_acceptInfo

asg::AcceptInfo JSSTaggerBase::m_acceptInfo
protected

Object that stores the results for a jet.

Definition at line 65 of file JSSTaggerBase.h.

◆ m_calcSF

bool JSSTaggerBase::m_calcSF {}
protected

Flag to calculate scale factor.

Definition at line 192 of file JSSTaggerBase.h.

◆ m_calibArea

std::string JSSTaggerBase::m_calibArea
protected

Location where config files live on cvmfs.

Definition at line 114 of file JSSTaggerBase.h.

◆ m_configFile

std::string JSSTaggerBase::m_configFile
protected

Configuration file name.

Definition at line 111 of file JSSTaggerBase.h.

◆ m_configReader

TEnv JSSTaggerBase::m_configReader
protected

TEnv instance to read config files.

Definition at line 62 of file JSSTaggerBase.h.

◆ m_containerName

std::string JSSTaggerBase::m_containerName
protected

Configurable members.

Jet container name

Definition at line 104 of file JSSTaggerBase.h.

◆ m_decC2Key

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

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"}
protected

Definition at line 219 of file JSSTaggerBase.h.

◆ m_decCutMLowKey

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

WriteDecorHandle keys for cut values.

Definition at line 218 of file JSSTaggerBase.h.

◆ m_decD2Key

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

Definition at line 147 of file JSSTaggerBase.h.

◆ m_decE3Key

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

Definition at line 148 of file JSSTaggerBase.h.

◆ m_decEfficiencyKey

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

Definition at line 225 of file JSSTaggerBase.h.

◆ m_decEffSFKey

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

Definition at line 226 of file JSSTaggerBase.h.

◆ m_decL2Key

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

Definition at line 149 of file JSSTaggerBase.h.

◆ m_decL3Key

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

Definition at line 150 of file JSSTaggerBase.h.

◆ m_decorationName

std::string JSSTaggerBase::m_decorationName
protected

Decoration name.

Definition at line 189 of file JSSTaggerBase.h.

◆ m_decPassMassKey

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

Definition at line 76 of file JSSTaggerBase.h.

◆ m_decPassScoreKey

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

Definition at line 77 of file JSSTaggerBase.h.

◆ m_decScoreCutKey

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

Definition at line 220 of file JSSTaggerBase.h.

◆ m_decScoreValueKey

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

Definition at line 221 of file JSSTaggerBase.h.

◆ m_decSigeffSFKey

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

Definition at line 227 of file JSSTaggerBase.h.

◆ m_decTaggedKey

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

WriteDecorHandle keys for tagging bools.

Definition at line 68 of file JSSTaggerBase.h.

◆ m_decTau21WTAKey

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

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"}
protected

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"}
protected

Definition at line 144 of file JSSTaggerBase.h.

◆ m_decValidEtaRangeKey

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

Definition at line 71 of file JSSTaggerBase.h.

◆ m_decValidEventContentKey

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

Definition at line 74 of file JSSTaggerBase.h.

◆ m_decValidJetContentKey

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

Definition at line 73 of file JSSTaggerBase.h.

◆ m_decValidKinRangeKey

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

Definition at line 72 of file JSSTaggerBase.h.

◆ m_decValidPtRangeHighKey

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

Definition at line 69 of file JSSTaggerBase.h.

◆ m_decValidPtRangeLowKey

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

Definition at line 70 of file JSSTaggerBase.h.

◆ m_decWeightKey

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

WriteDecorHandle keys for SF.

Definition at line 224 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
protected

Definition at line 199 of file JSSTaggerBase.h.

◆ m_efficiencyHistograms

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

Definition at line 205 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
protected

Definition at line 185 of file JSSTaggerBase.h.

◆ m_funcMassCutLow

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

TF1 for cut functions.

Definition at line 184 of file JSSTaggerBase.h.

◆ m_funcScoreCut

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

Definition at line 186 of file JSSTaggerBase.h.

◆ m_isMC

bool JSSTaggerBase::m_isMC {}
protected

Definition at line 193 of file JSSTaggerBase.h.

◆ m_isSmallRJet

bool JSSTaggerBase::m_isSmallRJet = false
protected

Definition at line 105 of file JSSTaggerBase.h.

◆ m_jetEtaMax

float JSSTaggerBase::m_jetEtaMax {}
protected

Definition at line 134 of file JSSTaggerBase.h.

◆ m_jetPtMax

float JSSTaggerBase::m_jetPtMax {}
protected

Definition at line 133 of file JSSTaggerBase.h.

◆ m_jetPtMin

float JSSTaggerBase::m_jetPtMin {}
protected

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
protected

Definition at line 120 of file JSSTaggerBase.h.

◆ m_kerasConfigFileName

std::string JSSTaggerBase::m_kerasConfigFileName
protected

Keras configurations for ML taggers.

Definition at line 117 of file JSSTaggerBase.h.

◆ m_kerasConfigFilePath

std::string JSSTaggerBase::m_kerasConfigFilePath
protected

Definition at line 118 of file JSSTaggerBase.h.

◆ m_kerasConfigOutputName

std::string JSSTaggerBase::m_kerasConfigOutputName
protected

Definition at line 119 of file JSSTaggerBase.h.

◆ m_nWarnMax

const int JSSTaggerBase::m_nWarnMax = 10
protected

Maximum number of warnings.

Definition at line 82 of file JSSTaggerBase.h.

◆ m_nWarnVar

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

Warning counters.

Definition at line 85 of file JSSTaggerBase.h.

◆ m_ptGeV

bool JSSTaggerBase::m_ptGeV = false
protected

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

Definition at line 89 of file JSSTaggerBase.h.

◆ m_readECF1Key

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

Definition at line 158 of file JSSTaggerBase.h.

◆ m_readECF2Key

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

Definition at line 159 of file JSSTaggerBase.h.

◆ m_readECF3Key

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

Definition at line 160 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"}
protected

Definition at line 172 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"}
protected

Definition at line 171 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"}
protected

Definition at line 170 of file JSSTaggerBase.h.

◆ m_readNtrk500Key

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

Definition at line 79 of file JSSTaggerBase.h.

◆ m_readParentKey

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

Definition at line 174 of file JSSTaggerBase.h.

◆ m_readQwKey

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

Definition at line 165 of file JSSTaggerBase.h.

◆ m_readSphericityKey

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

Definition at line 168 of file JSSTaggerBase.h.

◆ m_readSplit12Key

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

Definition at line 162 of file JSSTaggerBase.h.

◆ m_readSplit23Key

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

Definition at line 163 of file JSSTaggerBase.h.

◆ m_readTau1WTAKey

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

ReadDecorHandle keys for JSS moments.

Definition at line 153 of file JSSTaggerBase.h.

◆ m_readTau2WTAKey

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

Definition at line 154 of file JSSTaggerBase.h.

◆ m_readTau3WTAKey

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

Definition at line 155 of file JSSTaggerBase.h.

◆ m_readTau4WTAKey

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

Definition at line 156 of file JSSTaggerBase.h.

◆ m_readThrustMajKey

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

Definition at line 167 of file JSSTaggerBase.h.

◆ m_readTruthLabelKey

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

Truth label ReadDecorHandle key.

Definition at line 215 of file JSSTaggerBase.h.

◆ m_strMassCutHigh

std::string JSSTaggerBase::m_strMassCutHigh
protected

Definition at line 180 of file JSSTaggerBase.h.

◆ m_strMassCutLow

std::string JSSTaggerBase::m_strMassCutLow
protected

Strings for cut functions.

Definition at line 179 of file JSSTaggerBase.h.

◆ m_strScoreCut

std::string JSSTaggerBase::m_strScoreCut
protected

Definition at line 181 of file JSSTaggerBase.h.

◆ m_suppressOutputDependence

bool JSSTaggerBase::m_suppressOutputDependence {}
protected

Definition at line 176 of file JSSTaggerBase.h.

◆ m_tagClass

TAGCLASS JSSTaggerBase::m_tagClass {Unknown}
protected

Definition at line 99 of file JSSTaggerBase.h.

◆ m_tagType

std::string JSSTaggerBase::m_tagType
protected

Definition at line 129 of file JSSTaggerBase.h.

◆ m_tmvaCalibArea

std::string JSSTaggerBase::m_tmvaCalibArea
protected

Definition at line 125 of file JSSTaggerBase.h.

◆ m_tmvaConfigFileName

std::string JSSTaggerBase::m_tmvaConfigFileName
protected

TMVA configurations for BDT taggers.

Definition at line 123 of file JSSTaggerBase.h.

◆ m_tmvaConfigFilePath

std::string JSSTaggerBase::m_tmvaConfigFilePath
protected

Definition at line 124 of file JSSTaggerBase.h.

◆ m_truthBosonContainerName

std::string JSSTaggerBase::m_truthBosonContainerName
protected

Definition at line 210 of file JSSTaggerBase.h.

◆ m_truthLabelName

std::string JSSTaggerBase::m_truthLabelName
protected

Definition at line 212 of file JSSTaggerBase.h.

◆ m_truthLabelUseTRUTH3

bool JSSTaggerBase::m_truthLabelUseTRUTH3 {}
protected

Truth label options.

Definition at line 208 of file JSSTaggerBase.h.

◆ m_truthParticleContainerName

std::string JSSTaggerBase::m_truthParticleContainerName
protected

Definition at line 209 of file JSSTaggerBase.h.

◆ m_truthTopQuarkContainerName

std::string JSSTaggerBase::m_truthTopQuarkContainerName
protected

Definition at line 211 of file JSSTaggerBase.h.

◆ m_useECFG

bool JSSTaggerBase::m_useECFG = false
protected

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
protected

Definition at line 139 of file JSSTaggerBase.h.

◆ m_useMassCut

bool JSSTaggerBase::m_useMassCut = false
protected

Flag to indicate if mass window cut is used.

Definition at line 92 of file JSSTaggerBase.h.

◆ m_useScoreCut

bool JSSTaggerBase::m_useScoreCut = false
protected

Flag to indicate if a discriminant score is used.

Definition at line 95 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
protected

Histograms for scale factors.

Definition at line 203 of file JSSTaggerBase.h.

◆ m_weightConfigPath

std::string JSSTaggerBase::m_weightConfigPath
protected

Path to the SF configuration root file.

Definition at line 108 of file JSSTaggerBase.h.

◆ m_weightDecorationName

std::string JSSTaggerBase::m_weightDecorationName
protected

String for scale factor decoration names.

Definition at line 196 of file JSSTaggerBase.h.

◆ m_weightFileName

std::string JSSTaggerBase::m_weightFileName
protected

Definition at line 197 of file JSSTaggerBase.h.

◆ m_weightFlavors

std::string JSSTaggerBase::m_weightFlavors
protected

Definition at line 200 of file JSSTaggerBase.h.

◆ m_weightHistogramName

std::string JSSTaggerBase::m_weightHistogramName
protected

Definition at line 198 of file JSSTaggerBase.h.

◆ m_weightHistograms

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

Definition at line 204 of file JSSTaggerBase.h.

◆ m_wkpt

std::string JSSTaggerBase::m_wkpt
protected

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:174
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
LargeRJetTruthLabel::TypeEnum
TypeEnum
Definition: LargeRJetLabelEnum.h:14
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
JSSTaggerBase::m_readTau4WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau4WTAKey
Definition: JSSTaggerBase.h:156
JSSTaggerBase::m_suppressOutputDependence
bool m_suppressOutputDependence
Definition: JSSTaggerBase.h:176
JSSTaggerBase::m_readSplit23Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSplit23Key
Definition: JSSTaggerBase.h:163
asg::AsgTool
Base class for the dual-use tool implementation classes.
Definition: AsgTool.h:47
Trk::L2
@ L2
Definition: AlignModuleList.h:32
JSSTaggerBase::m_readECF1Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF1Key
Definition: JSSTaggerBase.h:158
JSSTaggerBase::getTruthLabelStr
std::string getTruthLabelStr(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Definition: JSSTaggerBase.cxx:827
JSSTaggerBase::m_decScoreCutKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decScoreCutKey
Definition: JSSTaggerBase.h:220
JSSTaggerBase::m_strMassCutLow
std::string m_strMassCutLow
Strings for cut functions.
Definition: JSSTaggerBase.h:179
get_generator_info.result
result
Definition: get_generator_info.py:21
StateLessPT_NewConfig.proxy
proxy
Definition: StateLessPT_NewConfig.py:407
JSSTaggerBase::m_funcScoreCut
std::unique_ptr< TF1 > m_funcScoreCut
Definition: JSSTaggerBase.h:186
JSSTaggerBase::ZBoson
@ ZBoson
Definition: JSSTaggerBase.h:98
PowhegControl_ttHplus_NLO.ss
ss
Definition: PowhegControl_ttHplus_NLO.py:83
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
JSSTaggerBase::m_readNtrk500Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readNtrk500Key
Definition: JSSTaggerBase.h:79
JSSTaggerBase::m_readECFG311Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG311Key
Definition: JSSTaggerBase.h:171
JSSTaggerBase::m_jetEtaMax
float m_jetEtaMax
Definition: JSSTaggerBase.h:134
JSSTaggerBase::m_nWarnVar
std::atomic< int > m_nWarnVar
Warning counters.
Definition: JSSTaggerBase.h:85
JSSTaggerBase::m_decValidEtaRangeKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidEtaRangeKey
Definition: JSSTaggerBase.h:71
JSSTaggerBase::m_kerasCalibArea
std::string m_kerasCalibArea
Definition: JSSTaggerBase.h:120
JSSTaggerBase::m_efficiencyHistograms
std::map< std::string, std::unique_ptr< TH2D > > m_efficiencyHistograms
Definition: JSSTaggerBase.h:205
JSSTaggerBase::m_configFile
std::string m_configFile
Configuration file name.
Definition: JSSTaggerBase.h:111
JSSTaggerBase::m_configReader
TEnv m_configReader
TEnv instance to read config files.
Definition: JSSTaggerBase.h:62
JSSTaggerBase::m_tmvaCalibArea
std::string m_tmvaCalibArea
Definition: JSSTaggerBase.h:125
RootHelpers::FindBin
Int_t FindBin(const TAxis *axis, const double x)
Definition: RootHelpers.cxx:14
JSSTaggerBase::TopQuark
@ TopQuark
Definition: JSSTaggerBase.h:98
JSSTaggerBase::m_weightFileName
std::string m_weightFileName
Definition: JSSTaggerBase.h:197
JSSTaggerBase::tag
virtual StatusCode tag(const xAOD::Jet &jet) const =0
Decorate single jet with tagging info.
AthCommonDataStore< AthCommonMsg< AlgTool > >::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 > renounce(T &h)
Definition: AthCommonDataStore.h:380
JSSTaggerBase::m_decPassScoreKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassScoreKey
Definition: JSSTaggerBase.h:77
python.SystemOfUnits.second
float second
Definition: SystemOfUnits.py:135
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
defineDB.jets
jets
Definition: JetTagCalibration/share/defineDB.py:24
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_readSplit12Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSplit12Key
Definition: JSSTaggerBase.h:162
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
SG::ConstAccessor
Helper class to provide constant type-safe access to aux data.
Definition: ConstAccessor.h:55
JSSTaggerBase::m_decCutMLowKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutMLowKey
WriteDecorHandle keys for cut values.
Definition: JSSTaggerBase.h:218
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:189
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
JSSTaggerBase::m_readECFG331Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG331Key
Definition: JSSTaggerBase.h:170
LargeRJetTruthLabel::Zqq
@ Zqq
Definition: LargeRJetLabelEnum.h:20
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
SG::VarHandleKey::key
const std::string & key() const
Return the StoreGate ID for the referenced object.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:141
dbg::ptr
void * ptr(T *p)
Definition: SGImplSvc.cxx:74
yodamerge_tmp.scale
scale
Definition: yodamerge_tmp.py:138
LargeRJetTruthLabel::qcd
@ qcd
Definition: LargeRJetLabelEnum.h:25
JSSTaggerBase::m_decValidPtRangeHighKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeHighKey
Definition: JSSTaggerBase.h:69
JSSTaggerBase::m_funcMassCutHigh
std::unique_ptr< TF1 > m_funcMassCutHigh
Definition: JSSTaggerBase.h:185
JSSTaggerBase::m_weightHistogramName
std::string m_weightHistogramName
Definition: JSSTaggerBase.h:198
JSSTaggerBase::m_decValidJetContentKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidJetContentKey
Definition: JSSTaggerBase.h:73
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:224
JSSTaggerBase::m_readTruthLabelKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTruthLabelKey
Truth label ReadDecorHandle key.
Definition: JSSTaggerBase.h:215
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:227
dqt_zlumi_pandas.weight
int weight
Definition: dqt_zlumi_pandas.py:190
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
JSSTaggerBase::m_truthLabelName
std::string m_truthLabelName
Definition: JSSTaggerBase.h:212
JSSTaggerBase::m_decEfficiencyKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decEfficiencyKey
Definition: JSSTaggerBase.h:225
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
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_decCutMHighKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutMHighKey
Definition: JSSTaggerBase.h:219
JSSTaggerBase::m_ptGeV
bool m_ptGeV
Flag to indicate units pT is defined in Set to false by default.
Definition: JSSTaggerBase.h:89
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
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:209
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
JSSTaggerBase::m_readQwKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readQwKey
Definition: JSSTaggerBase.h:165
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:769
JSSTaggerBase::m_weightConfig
std::unique_ptr< TFile > m_weightConfig
Histograms for scale factors.
Definition: JSSTaggerBase.h:203
JSSTaggerBase::m_kerasConfigOutputName
std::string m_kerasConfigOutputName
Definition: JSSTaggerBase.h:119
Trk::L3
@ L3
Definition: AlignModuleList.h:32
JSSTaggerBase::m_readECFG212Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG212Key
Definition: JSSTaggerBase.h:172
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
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition: AthCommonDataStore.h:145
SG::WriteDecorHandle
Handle class for adding a decoration to an object.
Definition: StoreGate/StoreGate/WriteDecorHandle.h:100
xAOD::VxType::PriVtx
@ PriVtx
Primary vertex.
Definition: TrackingPrimitives.h:572
CalibDbCompareRT.dummy
dummy
Definition: CalibDbCompareRT.py:59
JSSTaggerBase::Unknown
@ Unknown
Definition: JSSTaggerBase.h:98
JSSTaggerBase::m_efficiencyHistogramName
std::string m_efficiencyHistogramName
Definition: JSSTaggerBase.h:199
JSSTaggerBase::m_strScoreCut
std::string m_strScoreCut
Definition: JSSTaggerBase.h:181
JSSTaggerBase::m_decTau42WTAKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau42WTAKey
Definition: JSSTaggerBase.h:144
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
JSSTaggerBase::m_truthBosonContainerName
std::string m_truthBosonContainerName
Definition: JSSTaggerBase.h:210
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:226
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
JSSTaggerBase::m_wkpt
std::string m_wkpt
Tagger information.
Definition: JSSTaggerBase.h:128
JSSTaggerBase::m_calcSF
bool m_calcSF
Flag to calculate scale factor.
Definition: JSSTaggerBase.h:192
JSSTaggerBase::m_decValidEventContentKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidEventContentKey
Definition: JSSTaggerBase.h:74
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
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_decD2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decD2Key
Definition: JSSTaggerBase.h:147
JSSTaggerBase::m_truthTopQuarkContainerName
std::string m_truthTopQuarkContainerName
Definition: JSSTaggerBase.h:211
LargeRJetTruthLabel::Wqq_From_t
@ Wqq_From_t
Definition: LargeRJetLabelEnum.h:21
JSSTaggerBase::m_readThrustMajKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readThrustMajKey
Definition: JSSTaggerBase.h:167
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
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
JSSTaggerBase::m_truthLabelUseTRUTH3
bool m_truthLabelUseTRUTH3
Truth label options.
Definition: JSSTaggerBase.h:208
JSSTaggerBase::m_decScoreValueKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decScoreValueKey
Definition: JSSTaggerBase.h:221
PathResolverFindCalibFile
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
Definition: PathResolver.cxx:283
JSSTaggerBase::m_isMC
bool m_isMC
Definition: JSSTaggerBase.h:193
JSSTaggerBase::m_weightDecorationName
std::string m_weightDecorationName
String for scale factor decoration names.
Definition: JSSTaggerBase.h:196
JSSTaggerBase::WBoson
@ WBoson
Definition: JSSTaggerBase.h:98
xAOD::Jet_v1
Class describing a jet.
Definition: Jet_v1.h:57
JSSTaggerBase::m_decValidKinRangeKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidKinRangeKey
Definition: JSSTaggerBase.h:72
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:204
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
SG::WriteDecorHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
h
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:338
asg::AcceptData::clear
void clear()
Clear all bits.
Definition: AcceptData.h:54
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:70
JSSTaggerBase::m_truthParticleContainerName
std::string m_truthParticleContainerName
Definition: JSSTaggerBase.h:209
JSSTaggerBase::m_readSphericityKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSphericityKey
Definition: JSSTaggerBase.h:168
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:159
JSSTaggerBase::m_readECF3Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF3Key
Definition: JSSTaggerBase.h:160
AthCommonMsg< AlgTool >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
LargeRJetTruthLabel::intToEnum
TypeEnum intToEnum(const int type)
Definition: LargeRJetLabelEnum.h:57
JSSTaggerBase::m_decPassMassKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassMassKey
Definition: JSSTaggerBase.h:76
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
JSSTaggerBase::m_isSmallRJet
bool m_isSmallRJet
Definition: JSSTaggerBase.h:105
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:629
JSSTaggerBase::m_tagType
std::string m_tagType
Definition: JSSTaggerBase.h:129
JSSTaggerBase::m_funcMassCutLow
std::unique_ptr< TF1 > m_funcMassCutLow
TF1 for cut functions.
Definition: JSSTaggerBase.h:184
dqt_zlumi_alleff_HIST.eff
int eff
Definition: dqt_zlumi_alleff_HIST.py:113
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:801
LargeRJetTruthLabel::tqqb
@ tqqb
Definition: LargeRJetLabelEnum.h:16
JSSTaggerBase::m_strMassCutHigh
std::string m_strMassCutHigh
Definition: JSSTaggerBase.h:180
JSSTaggerBase::m_kerasConfigFileName
std::string m_kerasConfigFileName
Keras configurations for ML taggers.
Definition: JSSTaggerBase.h:117
JSSTaggerBase::m_decL2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL2Key
Definition: JSSTaggerBase.h:149
JSSTaggerBase::m_weightConfigPath
std::string m_weightConfigPath
Path to the SF configuration root file.
Definition: JSSTaggerBase.h:108
JSSTaggerBase::m_weightFlavors
std::string m_weightFlavors
Definition: JSSTaggerBase.h:200
SG::DataProxy
Definition: DataProxy.h:45
pow
constexpr int pow(int base, int exp) noexcept
Definition: ap_fixedTest.cxx:15
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:68
asg::AcceptInfo::addCut
int addCut(const std::string &cutName, const std::string &cutDescription)
Add a cut; returning the cut position.
Definition: AcceptInfo.h:53
JSSTaggerBase::m_containerName
std::string m_containerName
Configurable members.
Definition: JSSTaggerBase.h:104