ATLAS Offline Software
Loading...
Searching...
No Matches
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>

Inherits asg::AsgTool, and IJetDecorator.

Inherited by BJT::qgTagger, CP::JetQGTagger, CP::JetQGTaggerBDT, JSSTaggerUtils, JSSWTopTaggerANN, JSSWTopTaggerDNN, LundNetTagger, SmoothedTopTagger, and SmoothedWZTagger.

Public Member Functions

virtual StatusCode initialize () override
 Initialize the tool.
virtual StatusCode decorate (const xAOD::JetContainer &jets) const override
 Decorate jet collection with tagging info.
virtual void print () const
 Print the state of the tool.
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysInitialize () override
 Perform system initialization for an algorithm.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () 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.
Additional helper functions, not directly mimicking Athena
template<class T>
const T * getProperty (const std::string &name) const
 Get one of the tool's properties.
const std::string & msg_level_name () const __attribute__((deprecated))
 A deprecated function for getting the message level's name.
const std::string & getName (const void *ptr) const
 Get the name of an object that is / should be in the event store.
SG::sgkey_t getKey (const void *ptr) const
 Get the (hashed) key of an object that is in the event store.

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.
 ~JSSTaggerBase ()
 Default destructor - to be used in all derived classes.
virtual StatusCode tag (const xAOD::Jet &jet) const =0
 Decorate single jet with tagging info.
StatusCode getConfigReader ()
 Get configReader StatusCode.
StatusCode resetCuts (asg::AcceptData &acceptData) const
 Reset cuts.
bool passKinRange (const xAOD::Jet &jet) const
 Check if jet passes kinematic constraints.
StatusCode checkKinRange (const xAOD::Jet &jet, asg::AcceptData &acceptData) const
 Check and record if jet passes kinematic constraints.
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.
void decorateJSSRatios (const xAOD::JetContainer &jets) const
int findPV () const
 Find the PV (to be used for Ntrk)
int GetUnGroomTracks (const xAOD::Jet &jet, int indexPV) const
 Retrieve Ntrk variable from the ungroomed parent jet.
StatusCode getWeight (const xAOD::Jet &jet, bool passSel, asg::AcceptData &acceptData) const
 Get SF weight.
std::pair< double, double > getSF (const xAOD::Jet &jet, const std::string &truthLabelStr) const
 Get scale factor and efficiency.
std::string getTruthLabelStr (const xAOD::Jet &jet, asg::AcceptData &acceptData) const
void printCuts () const
 Print configured cuts.
void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
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.

Protected Attributes

TEnv m_configReader
 TEnv instance to read config files.
asg::AcceptInfo m_acceptInfo
 Object that stores the results for a jet.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTaggedKey {this, "TaggedName", "Tagged", "SG key for Tagged"}
 WriteDecorHandle keys for tagging bools.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeHighKey {this, "ValidPtRangeHighName", "ValidPtRangeHigh", "SG key for ValidPtRangeHigh"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeLowKey {this, "ValidPtRangeLowName", "ValidPtRangeLow", "SG key for ValidPtRangeLow"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidEtaRangeKey {this, "ValidEtaRangeName", "ValidEtaRange", "SG key for ValidEtaRange"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidKinRangeKey {this, "ValidKinRangeName", "ValidKinRange", "SG key for ValidKinRange"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidJetContentKey {this, "ValidJetContentName", "ValidJetContent", "SG key for ValidJetContent"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidEventContentKey {this, "ValidEventContentName", "ValidEventContent", "SG key for ValidEventContent"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassMassKey {this, "PassMassName", "PassMass", "SG key for PassMass"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassScoreKey {this, "PassScoreName", "PassScore", "SG key for PassScore"}
const int m_nWarnMax = 10
 Maximum number of warnings.
std::atomic< int > m_nWarnVar
 Warning counters.
bool m_ptGeV = false
 Flag to indicate units pT is defined in Set to false by default.
bool m_useMassCut = false
 Flag to indicate if mass window cut is used.
bool m_useScoreCut = false
 Flag to indicate if a discriminant score is used.
TAGCLASS m_tagClass {Unknown}
std::string m_containerName
 Configurable members.
bool m_isSmallRJet = false
std::string m_weightConfigPath
 Path to the SF configuration root file.
std::string m_configFile
 Configuration file name.
std::string m_calibArea
 Location where config files live on cvmfs.
std::string m_kerasConfigFileName
 Keras configurations for ML taggers.
std::string m_kerasConfigFilePath
std::string m_kerasConfigOutputName
std::string m_kerasCalibArea
std::string m_tmvaConfigFileName
 TMVA configurations for BDT taggers.
std::string m_tmvaConfigFilePath
std::string m_tmvaCalibArea
std::string m_wkpt
 Tagger information.
std::string m_tagType
float m_jetPtMin {}
 Kinematic bounds for the jet - the units are controlled by m_ptGeV.
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.
bool m_useLSeries = false
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau21WTAKey {this, "Tau21WTANameOut", "Tau21_wta", "SG key for Tau21_wta"}
 WriteDecorHandle keys for JSS moments.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau32WTAKey {this, "Tau32WTANameOut", "Tau32_wta", "SG key for Tau32_wta"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau42WTAKey {this, "Tau42WTANameOut", "Tau42_wta", "SG key for Tau42_wta"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decC2Key {this, "C2NameOut", "C2", "SG key for C2"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decD2Key {this, "D2NameOut", "D2", "SG key for D2"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decE3Key {this, "e3NameOut", "e3", "SG key for e3"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL2Key {this, "L2NameOut", "L2", "SG key for L2"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL3Key {this, "L3NameOut", "L3", "SG key for L3"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau1WTAKey {this, "Tau1WTAName", "Tau1_wta", "SG key for Tau1_wta"}
 ReadDecorHandle keys for JSS moments.
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau2WTAKey {this, "Tau2WTAName", "Tau2_wta", "SG key for Tau2_wta"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau3WTAKey {this, "Tau3WTAName", "Tau3_wta", "SG key for Tau3_wta"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau4WTAKey {this, "Tau4WTAName", "Tau4_wta", "SG key for Tau4_wta"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF1Key {this, "ECF1Name", "ECF1", "SG key for ECF1"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF2Key {this, "ECF2Name", "ECF2", "SG key for ECF2"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF3Key {this, "ECF3Name", "ECF3", "SG key for ECF3"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSplit12Key {this, "Split12Name", "Split12", "SG key for Split12"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSplit23Key {this, "Split23Name", "Split23", "SG key for Split23"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readQwKey {this, "QwName", "Qw", "SG key for Qw"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readThrustMajKey {this, "ThrustMajName", "ThrustMaj", "SG key for ThrustMaj"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSphericityKey {this, "SphericityName", "Sphericity", "SG key for Sphericity"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG331Key {this, "ECFG331Name", "ECFG_3_3_1", "SG key for ECFG_3_3_1"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG311Key {this, "ECFG311Name", "ECFG_3_1_1", "SG key for ECFG_3_1_1"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG212Key {this, "ECFG212Name", "ECFG_2_1_2", "SG key for ECFG_2_1_2"}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readParentKey {this, "ParentName", "Parent", "SG key for Parent"}
Gaudi::Property< bool > m_suppressInputDependence {this, "SuppressInputDependence", false}
Gaudi::Property< bool > m_suppressOutputDependence {this, "SuppressOutputDependence", false}
std::string m_strMassCutLow
 Strings for cut functions.
std::string m_strMassCutHigh
std::string m_strScoreCut
std::unique_ptr< TF1 > m_funcMassCutLow
 TF1 for cut functions.
std::unique_ptr< TF1 > m_funcMassCutHigh
std::unique_ptr< TF1 > m_funcScoreCut
std::string m_decorationName
 Decoration name.
bool m_calcSF {}
 Flag to calculate scale factor.
bool m_isMC {}
std::string m_weightDecorationName
 String for scale factor decoration names.
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.
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.
std::string m_truthParticleContainerName
std::string m_truthBosonContainerName
std::string m_truthTopQuarkContainerName
std::string m_truthLabelName
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTruthLabelKey {this, "truthLabelName", "truthLabel", "SG key for truthLabel"}
 Truth label ReadDecorHandle key.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutMLowKey {this, "CutMLowName", "Cut_mlow", "SG key for Cut_mlow"}
 WriteDecorHandle keys for cut values.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutMHighKey {this, "CutMHighName", "Cut_mhigh", "SG key for Cut_mhigh"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decScoreCutKey {this, "CutScoreName", "Cut_score", "SG key for Cut_score"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decScoreValueKey {this, "ScoreName", "Score", "SG key for Score"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decWeightKey {this, "weightName", "weight", "SG key for weight"}
 WriteDecorHandle keys for SF.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decEfficiencyKey {this, "efficiencyName", "efficiency", "SG key for efficiency"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decEffSFKey {this, "effSFName", "effSF", "SG key for effSF"}
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decSigeffSFKey {this, "sigeffSFName", "sigeffSF", "SG key for effSF"}

Private Types

typedef ServiceHandle< StoreGateSvc > StoreGateSvc_t

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default)
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared

Detailed Description

Definition at line 38 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 97 of file JSSTaggerBase.h.

97{Unknown, WBoson, ZBoson, TopQuark};

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 :
12 asg::AsgTool(name),
13 m_calibArea(""),
14 m_strMassCutLow(""),
15 m_strMassCutHigh(""),
16 m_strScoreCut("")
17{
18
20 declareProperty( "ContainerName", m_containerName = "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets", "Name of jet container" );
21 declareProperty( "ConfigFile", m_configFile = "", "Name of config file" );
22 declareProperty( "CalibArea", m_calibArea = "", "Path to config file" );
23 declareProperty( "CalcSF", m_calcSF = false, "Flag to calculate the efficiency SF" );
24 declareProperty( "WorkingPoint", m_wkpt = "", "If specified, name of working point is added to the variable names" );
25 declareProperty( "TaggerType", m_tagType = "XXX", "Tagger type (e.g. SmoothedWZTagger, JSSWTopTaggerDNN, etc.)" );
26 declareProperty( "Decoration", m_decorationName = "XX", "Prefix for the variables decorated to xAOD::Jet" );
27 declareProperty( "IsMC", m_isMC = true, "Flag to identify data or MC" );
28
30 declareProperty( "JetEtaMax", m_jetEtaMax = 2.0, "Eta cut to define fiducial phase space for the tagger");
31
33 declareProperty( "UseTRUTH3", m_truthLabelUseTRUTH3 = true,
34 "Flag to use TRUTH3 containers. If false, TRUTH1 format is used." );
35 declareProperty( "TruthParticleContainerName", m_truthParticleContainerName = "TruthParticles",
36 "Name of truth-particle container (with UseTRUTH3=false). TruthParticles by default" );
37 declareProperty( "TruthBosonContainerName", m_truthBosonContainerName = "TruthBosonsWithDecayParticles",
38 "Name of truth-boson container (with UseTRUTH3=true). TruthBosonWithDecayParticles by default" );
39 declareProperty( "TruthTopQuarkContainerName", m_truthTopQuarkContainerName = "TruthTopQuarkWithDecayParticles",
40 "Name of truth-top container (with UseTRUTH3=true). TruthTopQuarkWithDecayParticles by default" );
41
43 declareProperty( "CalibAreaKeras", m_kerasCalibArea = "BoostedJetTaggers/TopoclusterTopTagger/Boost2017/",
44 "Path to json file to configure ML-taggers (Keras)" );
45 declareProperty( "KerasConfigFile", m_kerasConfigFileName = "XXX",
46 "Name of json file to configure ML-taggers (Keras)" );
47 declareProperty( "KerasOutput", m_kerasConfigOutputName = "XXX",
48 "Name of output variable by the ML-tagger (Keras)" );
49
51 declareProperty( "CalibAreaTMVA", m_tmvaCalibArea = "BoostedJetTaggers/JSSWTopTaggerBDT/Boost2017/",
52 "Path to xml file to configure ML-taggers (TMVA)" );
53 declareProperty( "TMVAConfigFile", m_tmvaConfigFileName = "XXX",
54 "Name of xml file to configure ML-taggers (TMVA)" );
55
57 declareProperty( "WeightDecorationName", m_weightDecorationName = "SF",
58 "Name of SF variable decorated to xAOD::Jet" );
59 declareProperty( "WeightFile", m_weightFileName = "",
60 "Name of config ROOT file for SF calculation" );
61 declareProperty( "WeightHistogramName", m_weightHistogramName = "",
62 "Name of SF histograms in the ROOT file" );
63 declareProperty( "EfficiencyHistogramName", m_efficiencyHistogramName = "",
64 "Name of efficiency histograms in the ROOT file" );
65 declareProperty( "WeightFlavors", m_weightFlavors = "",
66 "List of jet flavours for which the SF is available. Divided by comma" );
67}
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition AthCommonDataStore.h:145
JSSTaggerBase::m_kerasCalibArea
std::string m_kerasCalibArea
Definition JSSTaggerBase.h:119
JSSTaggerBase::m_calcSF
bool m_calcSF
Flag to calculate scale factor.
Definition JSSTaggerBase.h:192
JSSTaggerBase::m_weightDecorationName
std::string m_weightDecorationName
String for scale factor decoration names.
Definition JSSTaggerBase.h:196
JSSTaggerBase::m_decorationName
std::string m_decorationName
Decoration name.
Definition JSSTaggerBase.h:189
JSSTaggerBase::m_kerasConfigFileName
std::string m_kerasConfigFileName
Keras configurations for ML taggers.
Definition JSSTaggerBase.h:116
JSSTaggerBase::m_wkpt
std::string m_wkpt
Tagger information.
Definition JSSTaggerBase.h:127
JSSTaggerBase::m_truthParticleContainerName
std::string m_truthParticleContainerName
Definition JSSTaggerBase.h:209
JSSTaggerBase::m_calibArea
std::string m_calibArea
Location where config files live on cvmfs.
Definition JSSTaggerBase.h:113
JSSTaggerBase::m_truthBosonContainerName
std::string m_truthBosonContainerName
Definition JSSTaggerBase.h:210
JSSTaggerBase::m_configFile
std::string m_configFile
Configuration file name.
Definition JSSTaggerBase.h:110
JSSTaggerBase::m_strMassCutHigh
std::string m_strMassCutHigh
Definition JSSTaggerBase.h:180
JSSTaggerBase::m_truthTopQuarkContainerName
std::string m_truthTopQuarkContainerName
Definition JSSTaggerBase.h:211
JSSTaggerBase::m_weightFileName
std::string m_weightFileName
Definition JSSTaggerBase.h:197
JSSTaggerBase::m_weightFlavors
std::string m_weightFlavors
Definition JSSTaggerBase.h:200
JSSTaggerBase::m_strMassCutLow
std::string m_strMassCutLow
Strings for cut functions.
Definition JSSTaggerBase.h:179
JSSTaggerBase::m_tagType
std::string m_tagType
Definition JSSTaggerBase.h:128
JSSTaggerBase::m_tmvaConfigFileName
std::string m_tmvaConfigFileName
TMVA configurations for BDT taggers.
Definition JSSTaggerBase.h:122
JSSTaggerBase::m_containerName
std::string m_containerName
Configurable members.
Definition JSSTaggerBase.h:103
JSSTaggerBase::m_strScoreCut
std::string m_strScoreCut
Definition JSSTaggerBase.h:181
JSSTaggerBase::m_efficiencyHistogramName
std::string m_efficiencyHistogramName
Definition JSSTaggerBase.h:199
JSSTaggerBase::m_weightHistogramName
std::string m_weightHistogramName
Definition JSSTaggerBase.h:198
JSSTaggerBase::m_truthLabelUseTRUTH3
bool m_truthLabelUseTRUTH3
Truth label options.
Definition JSSTaggerBase.h:208
JSSTaggerBase::m_kerasConfigOutputName
std::string m_kerasConfigOutputName
Definition JSSTaggerBase.h:118
JSSTaggerBase::m_tmvaCalibArea
std::string m_tmvaCalibArea
Definition JSSTaggerBase.h:124

◆ ~JSSTaggerBase()

JSSTaggerBase::~JSSTaggerBase ( )
inlineprotected

Default destructor - to be used in all derived classes.

Definition at line 57 of file JSSTaggerBase.h.

57{};

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

WTA N-subjettiness ratios

ECF ratios

Definition at line 410 of file JSSTaggerBase.cxx.

410 {
411
412 int result = 0;
413
415 SG::WriteDecorHandle<xAOD::JetContainer, float> decTau21WTA(m_decTau21WTAKey);
416 SG::WriteDecorHandle<xAOD::JetContainer, float> decTau32WTA(m_decTau32WTAKey);
417 SG::WriteDecorHandle<xAOD::JetContainer, float> decTau42WTA(m_decTau42WTAKey);
418 SG::WriteDecorHandle<xAOD::JetContainer, float> decC2(m_decC2Key);
419 SG::WriteDecorHandle<xAOD::JetContainer, float> decD2(m_decD2Key);
420 SG::WriteDecorHandle<xAOD::JetContainer, float> decE3(m_decE3Key);
421 SG::WriteDecorHandle<xAOD::JetContainer, float> decL2(m_decL2Key);
422 SG::WriteDecorHandle<xAOD::JetContainer, float> decL3(m_decL3Key);
423
425 float tau21_wta = -999.0;
426 float tau32_wta = -999.0;
427 float tau42_wta = -999.0;
428
429 if(!m_isSmallRJet){
430 SG::ReadDecorHandle<xAOD::JetContainer, float> readTau1WTA(m_readTau1WTAKey);
431 SG::ReadDecorHandle<xAOD::JetContainer, float> readTau2WTA(m_readTau2WTAKey);
432 SG::ReadDecorHandle<xAOD::JetContainer, float> readTau3WTA(m_readTau3WTAKey);
433 SG::ReadDecorHandle<xAOD::JetContainer, float> readTau4WTA(m_readTau4WTAKey);
434
435 float tau1_wta = readTau1WTA(jet);
436 float tau2_wta = readTau2WTA(jet);
437 float tau3_wta = readTau3WTA(jet);
438 float tau4_wta = readTau4WTA(jet);
439
440 if ( tau1_wta > 1e-8 ) {
441 tau21_wta = tau2_wta / tau1_wta;
442 }
443 else result = 1;
444
445 if ( tau2_wta > 1e-8 ) {
446 tau32_wta = tau3_wta / tau2_wta;
447 tau42_wta = tau4_wta / tau2_wta;
448 }
449 else result = 1;
450 }
451
452 decTau21WTA(jet) = tau21_wta;
453 decTau32WTA(jet) = tau32_wta;
454 decTau42WTA(jet) = tau42_wta;
455
457 float C2 = -999.0;
458 float D2 = -999.0;
459 float e3 = -999.0;
460
461 // L-series for UFO top taggers
462 float L2 = -999.0;
463 float L3 = -999.0;
464
465 if(!m_isSmallRJet){
466 SG::ReadDecorHandle<xAOD::JetContainer, float> readECF1(m_readECF1Key);
467 SG::ReadDecorHandle<xAOD::JetContainer, float> readECF2(m_readECF2Key);
468 SG::ReadDecorHandle<xAOD::JetContainer, float> readECF3(m_readECF3Key);
469
470 float ECF1 = readECF1(jet);
471 float ECF2 = readECF2(jet);
472 float ECF3 = readECF3(jet);
473
474 if ( ECF2 > 1e-8 ) {
475 C2 = ECF3 * ECF1 / std::pow( ECF2, 2.0 );
476 D2 = ECF3 * std::pow( ECF1, 3.0 ) / std::pow( ECF2, 3.0 );
477 }
478 else result = 1;
479
480 e3 = ECF3 / std::pow( ECF1, 3.0 );
481
482 SG::ReadDecorHandle<xAOD::JetContainer, float> readECFG331(m_readECFG331Key);
483 SG::ReadDecorHandle<xAOD::JetContainer, float> readECFG311(m_readECFG311Key);
484 SG::ReadDecorHandle<xAOD::JetContainer, float> readECFG212(m_readECFG212Key);
485
486 if(readECFG331.isAvailable() && readECFG212.isAvailable()){
487 if(readECFG212(jet) > 1e-8){
488 L2 = readECFG331(jet) / std::pow(readECFG212(jet), 1.5);
489 }
490 else result = 1;
491 }
492
493 if(readECFG331.isAvailable() && readECFG311.isAvailable()){
494 if(readECFG331(jet) > 1e-8){
495 L3 = readECFG311(jet) / std::pow(readECFG331(jet), 1./3.);
496 }
497 else result = 1;
498 }
499 }
500
501 decC2(jet) = C2;
502 decD2(jet) = D2;
503 decE3(jet) = e3;
504
505 static const SG::AuxElement::ConstAccessor<float> accL2("L2");
506 if(!accL2.isAvailable(jet)) decL2(jet) = L2;
507
508 static const SG::AuxElement::ConstAccessor<float> accL3("L3");
509 if(!accL3.isAvailable(jet)) decL3(jet) = L3;
510
511 // TODO: Add ECFG for ANN tagger whenever it is defined
512
513 return result;
514
515}
JSSTaggerBase::m_readECF2Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF2Key
Definition JSSTaggerBase.h:158
JSSTaggerBase::m_readECFG331Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG331Key
Definition JSSTaggerBase.h:169
JSSTaggerBase::m_decTau21WTAKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau21WTAKey
WriteDecorHandle keys for JSS moments.
Definition JSSTaggerBase.h:141
JSSTaggerBase::m_readECF1Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF1Key
Definition JSSTaggerBase.h:157
JSSTaggerBase::m_readTau4WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau4WTAKey
Definition JSSTaggerBase.h:155
JSSTaggerBase::m_readECFG311Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG311Key
Definition JSSTaggerBase.h:170
JSSTaggerBase::m_readTau3WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau3WTAKey
Definition JSSTaggerBase.h:154
JSSTaggerBase::m_readTau2WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau2WTAKey
Definition JSSTaggerBase.h:153
JSSTaggerBase::m_readECF3Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF3Key
Definition JSSTaggerBase.h:159
JSSTaggerBase::m_decC2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decC2Key
Definition JSSTaggerBase.h:145
JSSTaggerBase::m_readTau1WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau1WTAKey
ReadDecorHandle keys for JSS moments.
Definition JSSTaggerBase.h:152
JSSTaggerBase::m_readECFG212Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG212Key
Definition JSSTaggerBase.h:171
JSSTaggerBase::m_decL2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL2Key
Definition JSSTaggerBase.h:148
JSSTaggerBase::m_decE3Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decE3Key
Definition JSSTaggerBase.h:147
JSSTaggerBase::m_decD2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decD2Key
Definition JSSTaggerBase.h:146
JSSTaggerBase::m_decTau42WTAKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau42WTAKey
Definition JSSTaggerBase.h:143
JSSTaggerBase::m_decL3Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL3Key
Definition JSSTaggerBase.h:149
JSSTaggerBase::m_decTau32WTAKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau32WTAKey
Definition JSSTaggerBase.h:142
get_generator_info.result
result
Definition get_generator_info.py:21

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

364 {
365
366 float scale = 1.0;
367 if ( m_ptGeV ) scale = 1.e3;
368
369 bool passKinRange = true;
370
373 if ( std::abs(jet.eta()) > m_jetEtaMax ) {
374 ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (|eta| < " << m_jetEtaMax << "). Jet eta = " << jet.eta() );
375 acceptData.setCutResult( "ValidEtaRange", false );
376 passKinRange = false;
377 }
378
379 if ( jet.pt() < m_jetPtMin * scale ) {
380 ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (pT > " << m_jetPtMin * scale / 1.e3 << "). Jet pT = " << jet.pt() / 1.e3 << " GeV" );
381 acceptData.setCutResult( "ValidPtRangeLow", false );
382 passKinRange = false;
383 }
384
385 if ( jet.pt() > m_jetPtMax * scale ) {
386 ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (pT < " << m_jetPtMax * scale / 1.e3 << "). Jet pT = " << jet.pt() / 1.e3 << " GeV" );
387 acceptData.setCutResult( "ValidPtRangeHigh", false );
388 passKinRange = false;
389 }
390
392 SG::WriteDecorHandle<xAOD::JetContainer, char> decValidPtRangeHigh(m_decValidPtRangeHighKey);
393 SG::WriteDecorHandle<xAOD::JetContainer, char> decValidPtRangeLow(m_decValidPtRangeLowKey);
394 SG::WriteDecorHandle<xAOD::JetContainer, char> decValidEtaRange(m_decValidEtaRangeKey);
395 SG::WriteDecorHandle<xAOD::JetContainer, char> decValidKinRange(m_decValidKinRangeKey);
396
398 decValidPtRangeHigh(jet) = acceptData.getCutResult( "ValidPtRangeHigh" );
399 decValidPtRangeLow(jet) = acceptData.getCutResult( "ValidPtRangeLow" );
400 decValidEtaRange(jet) = acceptData.getCutResult( "ValidEtaRange" );
401 decValidKinRange(jet) = passKinRange;
402
403 return StatusCode::SUCCESS;
404
405}
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition AthMsgStreamMacros.h:28
JSSTaggerBase::m_decValidPtRangeHighKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeHighKey
Definition JSSTaggerBase.h:70
JSSTaggerBase::m_jetPtMin
float m_jetPtMin
Kinematic bounds for the jet - the units are controlled by m_ptGeV.
Definition JSSTaggerBase.h:131
JSSTaggerBase::m_decValidPtRangeLowKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeLowKey
Definition JSSTaggerBase.h:71
JSSTaggerBase::m_decValidKinRangeKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidKinRangeKey
Definition JSSTaggerBase.h:73
JSSTaggerBase::m_decValidEtaRangeKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidEtaRangeKey
Definition JSSTaggerBase.h:72
JSSTaggerBase::passKinRange
bool passKinRange(const xAOD::Jet &jet) const
Check if jet passes kinematic constraints.
Definition JSSTaggerBase.cxx:350
JSSTaggerBase::m_ptGeV
bool m_ptGeV
Flag to indicate units pT is defined in Set to false by default.
Definition JSSTaggerBase.h:88
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
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
xAOD::Jet_v1::pt
virtual double pt() const
The transverse momentum ( ) of the particle.
Definition Jet_v1.cxx:44
xAOD::Jet_v1::eta
virtual double eta() const
The pseudorapidity ( ) of the particle.
Definition Jet_v1.cxx:49
yodamerge_tmp.scale
scale
Definition yodamerge_tmp.py:138

◆ declareGaudiProperty()

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 }

◆ declareProperty()

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;
147 return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
148 }
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

◆ 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 BJT::qgTagger, JSSWTopTaggerANN, JSSWTopTaggerDNN, LundNetTagger, SmoothedTopTagger, and SmoothedWZTagger.

Definition at line 282 of file JSSTaggerBase.cxx.

282 {
283
284 for ( const auto *jet : jets ) {
285 ATH_CHECK( tag(*jet) );
286 }
287
288 return StatusCode::SUCCESS;
289
290}
ATH_CHECK
#define ATH_CHECK
Evaluate an expression and check for errors.
Definition AthCheckMacros.h:40
JSSTaggerBase::tag
virtual StatusCode tag(const xAOD::Jet &jet) const =0
Decorate single jet with tagging info.

◆ decorateJSSRatios()

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

Create write decor handles

WTA N-subjettiness ratios

ECF ratios

Definition at line 518 of file JSSTaggerBase.cxx.

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

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

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

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

641 {
642
643 int indexPV = -1;
644
645 const xAOD::VertexContainer* vxCont = nullptr;
646 if ( evtStore()->retrieve( vxCont, "PrimaryVertices" ).isFailure() ) {
647 ATH_MSG_WARNING( "Unable to retrieve primary vertex container PrimaryVertices" );
648 }
649 else {
650 int vtx_counter = 0;
651 for ( const auto *vx : *vxCont ) {
652 if ( vx->vertexType()==xAOD::VxType::PriVtx ) {
653 indexPV = vtx_counter;
654 break;
655 }
656 vtx_counter++;
657 }
658 }
659
660 return indexPV;
661
662}
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition AthMsgStreamMacros.h:32
AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
Definition AthCommonDataStore.h:85
xAOD::VxType::PriVtx
@ PriVtx
Primary vertex.
Definition TrackingPrimitives.h:572
xAOD::VertexContainer
VertexContainer_v1 VertexContainer
Definition of the current "Vertex container version".
Definition VertexContainer.h:14

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

293 {
294
295 ATH_MSG_INFO( "Using config file : " << m_configFile );
296
298 std::string configPath;
299
300 if ( m_calibArea.compare("Local") == 0 ) {
301 configPath = PathResolverFindCalibFile(m_configFile);
302 }
303 else if ( m_calibArea.find("eos") != std::string::npos) {
304 configPath = PathResolverFindCalibFile((m_calibArea+"/"+m_configFile).c_str());
305 configPath = (m_calibArea+"/"+m_configFile).c_str();
306 }
307 else {
308 configPath = PathResolverFindCalibFile(("BoostedJetTaggers/"+m_calibArea+"/"+m_configFile).c_str());
309 }
310
312 FileStat_t fStats;
313 int fSuccess = gSystem->GetPathInfo(configPath.c_str(), fStats);
314 if ( fSuccess ) {
315 ATH_MSG_ERROR( "Recommendations file " << m_configFile << " could not be found" );
316 return StatusCode::FAILURE;
317 }
318 else {
319 ATH_MSG_DEBUG( "Recommendations file was found : " << configPath );
320 }
321
322 if ( m_configReader.ReadFile( configPath.c_str(), EEnvLevel(0) ) ) {
323 ATH_MSG_ERROR( "Error while reading config file : "<< configPath );
324 return StatusCode::FAILURE;
325 }
326
327 return StatusCode::SUCCESS;
328
329}
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition AthMsgStreamMacros.h:33
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition AthMsgStreamMacros.h:31
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition AthMsgStreamMacros.h:29
PathResolverFindCalibFile
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
Definition PathResolver.cxx:325
JSSTaggerBase::m_configReader
TEnv m_configReader
TEnv instance to read config files.
Definition JSSTaggerBase.h:63

◆ 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 }
StateLessPT_NewConfig.proxy
proxy
Definition StateLessPT_NewConfig.py:415

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

777 {
778
779 if ( !passKinRange(jet) ) return std::make_pair( 1.0, 1.0 );
780
781
782 double logmOverPt = std::log(jet.m()/jet.pt());
783 if ( m_decorationName.find("SmoothZ") != std::string::npos ||
784 m_decorationName.find("SmoothInclusiveZ") != std::string::npos ) {
786 const double WtoZmassShift = 10803;
787 logmOverPt = std::log((jet.m()-WtoZmassShift)/jet.pt());
788 }
789
790 if ( logmOverPt > 0 ) logmOverPt = 0;
791
792 double SF = 1.0;
793 double eff = 1.0;
794
795 if ( m_weightHistograms.count(truthLabelStr.c_str()) ) {
796
797 int pt_mPt_bin = (m_weightHistograms.find(truthLabelStr.c_str())->second)->FindBin(jet.pt()*0.001, logmOverPt);
798 SF = (m_weightHistograms.find(truthLabelStr.c_str())->second)->GetBinContent(pt_mPt_bin);
799
800 if ( !m_efficiencyHistogramName.empty() ) {
801 eff = (m_efficiencyHistograms.find(truthLabelStr.c_str())->second)->GetBinContent(pt_mPt_bin);
802 }
803
804 }
805 else {
806 // set the efficiency for "Other" category to be the signal efficiency
807 std::string signal_truthLabel="";
808 if ( m_weightHistograms.count("t_qqb") ) {
809 signal_truthLabel="t_qqb";
810 }else if ( m_weightHistograms.count("V_qq") ){
811 signal_truthLabel="V_qq";
812 }else if ( m_weightHistograms.count("t") ){
813 signal_truthLabel="t";
814 }
815 if ( !signal_truthLabel.empty() && !m_efficiencyHistogramName.empty() ){
816 int pt_mPt_bin = (m_weightHistograms.find(signal_truthLabel.c_str())->second)->FindBin(jet.pt()*0.001, logmOverPt);
817 eff = (m_efficiencyHistograms.find(signal_truthLabel.c_str())->second)->GetBinContent(pt_mPt_bin);
818 }
819
820 return std::make_pair( 1.0, eff );
821 }
822
823 if ( SF < 1e-3 ) {
824 ATH_MSG_DEBUG( "(pt, m/pt) (" << jet.pt()/1.e3 << ", " << jet.m()/jet.pt() << ") is out of range for SF calculation. Returning 1.0" );
825 return std::make_pair( 1.0, 1.0 );
826 }
827 else {
828 return std::make_pair( SF, eff );
829 }
830
831}
JSSTaggerBase::m_efficiencyHistograms
std::map< std::string, std::unique_ptr< TH2D > > m_efficiencyHistograms
Definition JSSTaggerBase.h:205
JSSTaggerBase::m_weightHistograms
std::map< std::string, std::unique_ptr< TH2D > > m_weightHistograms
Definition JSSTaggerBase.h:204
xAOD::Jet_v1::m
virtual double m() const
The invariant mass of the particle.
Definition Jet_v1.cxx:59
TRT_PAI_gasdata::SF
const float SF[NF]
Cross sections for Fluor.
Definition TRT_PAI_gasdata.h:285
dqt_zlumi_alleff_HIST.eff
int eff
Definition dqt_zlumi_alleff_HIST.py:113

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

835 {
837 std::string truthLabelStr;
838
840 SG::ReadDecorHandle<xAOD::JetContainer, int> readTruthLabel(m_readTruthLabelKey);
841 LargeRJetTruthLabel::TypeEnum jetContainment = LargeRJetTruthLabel::intToEnum(readTruthLabel(jet));
842
844 if ( m_weightHistograms.count("t_qqb") ) {
845
847 if ( jetContainment==LargeRJetTruthLabel::tqqb ) {
848 truthLabelStr = "t_qqb";
849 }
851 else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
852 truthLabelStr = "q";
853 }
854
855 }
857 else if ( m_weightHistograms.count("V_qq_passMpassD2") ) {
858
860 if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t ) {
861 truthLabelStr = "t_";
862 }
864 else if ( jetContainment==LargeRJetTruthLabel::Wqq || jetContainment==LargeRJetTruthLabel::Zqq || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
865 truthLabelStr = "V_qq_";
866 }
868 else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
869 truthLabelStr = "q_";
870 }
871
873 if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") && acceptData.getCutResult("PassD2") ) {
874 truthLabelStr += "passMpassD2";
875 }
877 else if ( !(acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh")) && acceptData.getCutResult("PassD2") ) {
878 truthLabelStr += "failMpassD2";
879 }
881 else if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") && !acceptData.getCutResult("PassD2") ) {
882 truthLabelStr += "passMfailD2";
883 }
885 else{
886 truthLabelStr += "failMfailD2";
887 }
888
889 }
890
892 else if ( m_weightHistograms.count("V_qq") ) {
893
895 if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t ) {
896 truthLabelStr = "t";
897 }
899 else if ( jetContainment==LargeRJetTruthLabel::Wqq || jetContainment==LargeRJetTruthLabel::Zqq || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
900 truthLabelStr = "V_qq";
901 }
903 else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
904 truthLabelStr = "q";
905 }
906 }
907
908 // inclusive top tagger
909 else {
911 if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
912 truthLabelStr = "t";
913 }
915 else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
916 truthLabelStr = "q";
917 }
918
919 }
920
921 return truthLabelStr;
922}
JSSTaggerBase::m_readTruthLabelKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTruthLabelKey
Truth label ReadDecorHandle key.
Definition JSSTaggerBase.h:215
LargeRJetTruthLabel::intToEnum
TypeEnum intToEnum(const int type)
Definition LargeRJetLabelEnum.h:63

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

665 {
666
667 SG::ReadDecorHandle<xAOD::JetContainer, ElementLink<xAOD::JetContainer> > readParent(m_readParentKey);
668
669 const xAOD::Jet * ungroomedJet = nullptr;
670 if ( readParent.isAvailable() ) {
671 const ElementLink<xAOD::JetContainer>& linkToUngroomed = readParent(jet);
672 if ( linkToUngroomed.isValid() ) {
673 ungroomedJet = *linkToUngroomed;
674
675 static const SG::AuxElement::ConstAccessor< std::vector<int> >acc_Ntrk("NumTrkPt500");
676
677 if ( acc_Ntrk.isAvailable(*ungroomedJet) ) {
678
679 const std::vector<int>& NTrkPt500 = acc_Ntrk(*ungroomedJet);
680
681 int jet_ntrk = NTrkPt500.at(indexPV);
682 return jet_ntrk;
683
684 }
685 else {
686 ATH_MSG_ERROR("WARNING: Unable to retrieve Ntrk of the ungroomed parent jet. Please make sure this variable is in your derivations!!!");
687 }
688 }
689 else {
690 ATH_MSG_ERROR("WARNING: Unable to retrieve the parent ungroomed jet. Please make sure this variable is in your derivations!!!");
691 }
692 }
693 else {
694 ATH_MSG_ERROR("WARNING: Unable to retrieve the link to the parent ungroomed jet. Please make sure this variable is in your derivations!!!");
695 }
696
697 return -999;
698
699}
JSSTaggerBase::m_readParentKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readParentKey
Definition JSSTaggerBase.h:173
SG::AuxElement::ConstAccessor
SG::ConstAccessor< T, ALLOC > ConstAccessor
Definition AuxElement.h:570

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

703 {
704 if ( !m_calcSF ) return StatusCode::SUCCESS;
705
706 float weight = 1.0;
707 float effSF = 1.0;
708 float sigeffSF = 1.0;
709 float efficiency = 1.0;
710
711 if ( m_isMC ) {
712
713 std::string truthLabelStr = getTruthLabelStr( jet, acceptData );
714 std::tie(effSF, efficiency) = getSF( jet, truthLabelStr );
715
716 // calculate signal efficiency SF
717 if ( m_weightHistograms.count("t_qqb") ) {
718 sigeffSF = getSF(jet, "t_qqb").first;
719 } else if ( m_weightHistograms.count("V_qq") ) {
720 sigeffSF = getSF(jet, "V_qq").first;
721 } else if ( m_weightHistograms.count("t") ){
722 sigeffSF = getSF(jet, "t").first;
723 } else {
724 sigeffSF = 1.0;
725 }
726
728 if ( m_weightFlavors.find("fail") != std::string::npos ) {
729 weight = effSF;
730 }
731
732 else {
733
735 if ( passSel ) {
736 weight = effSF;
737 }
738
740 else {
742 if ( m_efficiencyHistogramName.empty() ) {
743 weight = 1.0;
744 }
745 else if ( efficiency < 1.0 ) {
746 weight = ( 1. - effSF * efficiency ) / ( 1. - efficiency );
747 }
748 else {
749 weight = 1.0;
750 }
751 }
752 }
753
754 }
755
756 else {
757 weight = 1.0;
758 }
759
761 SG::WriteDecorHandle<xAOD::JetContainer, float> decWeight(m_decWeightKey);
762 SG::WriteDecorHandle<xAOD::JetContainer, float> decEfficiency(m_decEfficiencyKey);
763 SG::WriteDecorHandle<xAOD::JetContainer, float> decEffSF(m_decEffSFKey);
764 SG::WriteDecorHandle<xAOD::JetContainer, float> decSigeffSF(m_decSigeffSFKey);
765
767 decWeight(jet) = weight;
768 decEfficiency(jet) = efficiency;
769 decEffSF(jet) = effSF;
770 decSigeffSF(jet) = sigeffSF;
771
772 return StatusCode::SUCCESS;
773
774}
JSSTaggerBase::getTruthLabelStr
std::string getTruthLabelStr(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Definition JSSTaggerBase.cxx:835
JSSTaggerBase::m_decWeightKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decWeightKey
WriteDecorHandle keys for SF.
Definition JSSTaggerBase.h:224
JSSTaggerBase::m_decEffSFKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decEffSFKey
Definition JSSTaggerBase.h:226
JSSTaggerBase::m_decEfficiencyKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decEfficiencyKey
Definition JSSTaggerBase.h:225
JSSTaggerBase::getSF
std::pair< double, double > getSF(const xAOD::Jet &jet, const std::string &truthLabelStr) const
Get scale factor and efficiency.
Definition JSSTaggerBase.cxx:777
JSSTaggerBase::m_decSigeffSFKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decSigeffSFKey
Definition JSSTaggerBase.h:227
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
dqt_zlumi_pandas.weight
int weight
Definition dqt_zlumi_pandas.py:190

◆ initialize()

StatusCode JSSTaggerBase::initialize ( void )
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 BJT::qgTagger, CP::JetQGTagger, CP::JetQGTaggerBDT, JSSTaggerUtils, JSSWTopTaggerANN, JSSWTopTaggerDNN, LundNetTagger, SmoothedTopTagger, and SmoothedWZTagger.

Definition at line 69 of file JSSTaggerBase.cxx.

69 {
70
72 if ( m_containerName.empty() ) {
73 ATH_MSG_ERROR( "ContainerName has not been set. Exiting" );
74 return StatusCode::FAILURE;
75 }
76 m_isSmallRJet = m_containerName.find("AntiKt4")!=std::string::npos;
77
79 m_nWarnVar = 0;
80
82 m_acceptInfo.addCut( "ValidPtRangeHigh" , "True if the jet is not too high pT" );
83 m_acceptInfo.addCut( "ValidPtRangeLow" , "True if the jet is not too low pT" );
84 m_acceptInfo.addCut( "ValidEtaRange" , "True if the jet is not too forward" );
85
86 m_acceptInfo.addCut( "ValidJetContent" , "True if the jet is alright technically (e.g. all attributes necessary for tag)" );
87 m_acceptInfo.addCut( "ValidEventContent" , "True if the event is alright technically (e.g. primary vertices)" );
88
90 ATH_MSG_INFO( "Decorators that will be attached to jet :" );
91
92 m_decTaggedKey = m_containerName + "." + m_decorationName + "_" + m_decTaggedKey.key();
93 m_decValidPtRangeHighKey = m_containerName + "." + m_decorationName + "_" + m_decValidPtRangeHighKey.key();
94 m_decValidPtRangeLowKey = m_containerName + "." + m_decorationName + "_" + m_decValidPtRangeLowKey.key();
95 m_decValidEtaRangeKey = m_containerName + "." + m_decorationName + "_" + m_decValidEtaRangeKey.key();
96 m_decValidKinRangeKey = m_containerName + "." + m_decorationName + "_" + m_decValidKinRangeKey.key();
97 m_decValidJetContentKey = m_containerName + "." + m_decorationName + "_" + m_decValidJetContentKey.key();
98 m_decValidEventContentKey = m_containerName + "." + m_decorationName + "_" + m_decValidEventContentKey.key();
99
100 m_decTau21WTAKey = m_containerName + "." + m_decTau21WTAKey.key();
101 m_decTau32WTAKey = m_containerName + "." + m_decTau32WTAKey.key();
102 m_decTau42WTAKey = m_containerName + "." + m_decTau42WTAKey.key();
103 m_decC2Key = m_containerName + "." + m_decC2Key.key();
104 m_decD2Key = m_containerName + "." + m_decD2Key.key();
105 m_decE3Key = m_containerName + "." + m_decE3Key.key();
106 m_decL2Key = m_containerName + "." + m_decL2Key.key();
107 m_decL3Key = m_containerName + "." + m_decL3Key.key();
108
109 if(!m_isSmallRJet){
110 m_readSplit12Key = m_containerName + "." + m_readSplit12Key.key();
111 m_readSplit23Key = m_containerName + "." + m_readSplit23Key.key();
112 m_readTau1WTAKey = m_containerName + "." + m_readTau1WTAKey.key();
113 m_readTau2WTAKey = m_containerName + "." + m_readTau2WTAKey.key();
114 m_readTau3WTAKey = m_containerName + "." + m_readTau3WTAKey.key();
115 m_readTau4WTAKey = m_containerName + "." + m_readTau4WTAKey.key();
116 m_readThrustMajKey = m_containerName + "." + m_readThrustMajKey.key();
117 m_readSphericityKey = m_containerName + "." + m_readSphericityKey.key();
118 m_readECF1Key = m_containerName + "." + m_readECF1Key.key();
119 m_readECF2Key = m_containerName + "." + m_readECF2Key.key();
120 m_readECF3Key = m_containerName + "." + m_readECF3Key.key();
121 m_readECFG331Key = m_containerName + "." + m_readECFG331Key.key();
122 m_readECFG311Key = m_containerName + "." + m_readECFG311Key.key();
123 m_readECFG212Key = m_containerName + "." + m_readECFG212Key.key();
124 m_readQwKey = m_containerName + "." + m_readQwKey.key();
125 m_readParentKey = m_containerName + "." + m_readParentKey.key();
126 }
127
128 ATH_CHECK( m_decTau21WTAKey.initialize() );
129 ATH_CHECK( m_decTau32WTAKey.initialize() );
130 ATH_CHECK( m_decTau42WTAKey.initialize() );
131 ATH_CHECK( m_decC2Key.initialize() );
132 ATH_CHECK( m_decD2Key.initialize() );
133 ATH_CHECK( m_decE3Key.initialize() );
134 ATH_CHECK( m_decL2Key.initialize() );
135 ATH_CHECK( m_decL3Key.initialize() );
136
137 ATH_CHECK( m_readTau1WTAKey.initialize(!m_isSmallRJet) );
138 ATH_CHECK( m_readTau2WTAKey.initialize(!m_isSmallRJet) );
139 ATH_CHECK( m_readTau3WTAKey.initialize(!m_isSmallRJet) );
140 ATH_CHECK( m_readTau4WTAKey.initialize(!m_isSmallRJet) );
141
142 ATH_CHECK( m_readSplit12Key.initialize(!m_isSmallRJet) );
143 ATH_CHECK( m_readSplit23Key.initialize(!m_isSmallRJet) );
144 ATH_CHECK( m_readThrustMajKey.initialize(!m_isSmallRJet) );
145 ATH_CHECK( m_readSphericityKey.initialize(!m_isSmallRJet) );
146 ATH_CHECK( m_readECF1Key.initialize(!m_isSmallRJet) );
147 ATH_CHECK( m_readECF2Key.initialize(!m_isSmallRJet) );
148 ATH_CHECK( m_readECF3Key.initialize(!m_isSmallRJet) );
149 ATH_CHECK( m_readECFG331Key.initialize(!m_isSmallRJet) );
150 ATH_CHECK( m_readECFG311Key.initialize(!m_isSmallRJet) );
151 ATH_CHECK( m_readECFG212Key.initialize(!m_isSmallRJet) );
152 ATH_CHECK( m_readQwKey.initialize(!m_isSmallRJet) );
153
154 ATH_CHECK( m_readParentKey.initialize(!m_isSmallRJet) );
155
156 ATH_CHECK( m_decTaggedKey.initialize() );
157 ATH_CHECK( m_decValidPtRangeHighKey.initialize() );
158 ATH_CHECK( m_decValidPtRangeLowKey.initialize() );
159 ATH_CHECK( m_decValidEtaRangeKey.initialize() );
160 ATH_CHECK( m_decValidKinRangeKey.initialize() );
161 ATH_CHECK( m_decValidJetContentKey.initialize() );
162 ATH_CHECK( m_decValidEventContentKey.initialize() );
163
164 ATH_MSG_INFO( " " << m_decTaggedKey.key() << " : pass tagging criteria" );
165 ATH_MSG_INFO( " " << m_decValidPtRangeHighKey.key() << " : pass upper pt range" );
166 ATH_MSG_INFO( " " << m_decValidPtRangeLowKey.key() << " : pass lower pt range" );
167 ATH_MSG_INFO( " " << m_decValidEtaRangeKey.key() << " : pass eta range" );
168 ATH_MSG_INFO( " " << m_decValidKinRangeKey.key() << " : pass kinematic range" );
169 ATH_MSG_INFO( " " << m_decValidJetContentKey.key() << " : has valid jet content" );
170 ATH_MSG_INFO( " " << m_decValidEventContentKey.key() << " : has valid event content" );
171
172
173 m_decPassMassKey = m_containerName + "." + m_decorationName + "_" + m_decPassMassKey.key();
174 m_decCutMLowKey = m_containerName + "." + m_decorationName + "_" + m_decCutMLowKey.key();
175 m_decCutMHighKey = m_containerName + "." + m_decorationName + "_" + m_decCutMHighKey.key();
176
177 ATH_CHECK( m_decPassMassKey.initialize() );
178 ATH_CHECK( m_decCutMLowKey.initialize() );
179 ATH_CHECK( m_decCutMHighKey.initialize() );
180
181 ATH_MSG_INFO( " " << m_decPassMassKey.key() << " : pass mass cut" );
182 ATH_MSG_INFO( " " << m_decCutMLowKey.key() << " : lower mass cut" );
183 ATH_MSG_INFO( " " << m_decCutMHighKey.key() << " : upper mass cut" );
184
185
186 m_decPassScoreKey = m_containerName + "." + m_decorationName + "_" + m_decPassScoreKey.key();
187 m_decScoreCutKey = m_containerName + "." + m_decorationName + "_" + m_decScoreCutKey.key();
188 m_decScoreValueKey = m_containerName + "." + m_decorationName + "_" + m_decScoreValueKey.key();
189
190 ATH_CHECK( m_decPassScoreKey.initialize() );
191 ATH_CHECK( m_decScoreCutKey.initialize() );
192 ATH_CHECK( m_decScoreValueKey.initialize() );
193
194 ATH_MSG_INFO( " " << m_decPassScoreKey.key() << " : pass MVA score cut" );
195 ATH_MSG_INFO( " " << m_decScoreCutKey.key() << " : MVA score cut" );
196 ATH_MSG_INFO( " " << m_decScoreValueKey.key() << " : evaluated MVA score" );
197
198
199 m_decWeightKey = m_containerName + "." + m_decorationName + "_" + m_weightDecorationName;
200 m_decEfficiencyKey = m_containerName + "." + m_decorationName + "_" + m_decEfficiencyKey.key();
201 m_decEffSFKey = m_containerName + "." + m_decorationName + "_" + m_decEffSFKey.key();
202 m_decSigeffSFKey = m_containerName + "." + m_decorationName + "_" + m_decSigeffSFKey.key();
203
204 ATH_CHECK( m_decWeightKey.initialize() );
205 ATH_CHECK( m_decEfficiencyKey.initialize() );
206 ATH_CHECK( m_decEffSFKey.initialize() );
207 ATH_CHECK( m_decSigeffSFKey.initialize() );
208
209 ATH_MSG_INFO( " " << m_decWeightKey.key() << " : tagging SF" );
210
211 m_readTruthLabelKey = m_containerName + "." + m_truthLabelName;
212 ATH_CHECK( m_readTruthLabelKey.initialize(!m_truthLabelName.empty()) ); // TODO Could also depend on m_isMC ?
213
214#ifndef XAOD_STANDALONE
215 if (m_suppressInputDependence) {
216 renounce(m_readQwKey);
217 renounce(m_readTau1WTAKey);
218 renounce(m_readTau2WTAKey);
219 renounce(m_readTau3WTAKey);
220 renounce(m_readTau4WTAKey);
221 renounce(m_readThrustMajKey);
222 renounce(m_readSphericityKey);
223 renounce(m_readECF1Key);
224 renounce(m_readECF2Key);
225 renounce(m_readECF3Key);
226 renounce(m_readSplit12Key);
227 renounce(m_readSplit23Key);
228 renounce(m_readECFG331Key);
229 renounce(m_readECFG311Key);
230 renounce(m_readECFG212Key);
231 renounce(m_readParentKey);
232 }
233 if (m_suppressOutputDependence) {
234 renounce(m_decValidPtRangeHighKey);
235 renounce(m_decValidPtRangeLowKey);
236 renounce(m_decValidEtaRangeKey);
237 renounce(m_decValidKinRangeKey);
238 renounce(m_decTau21WTAKey);
239 renounce(m_decTau32WTAKey);
240 renounce(m_decTau42WTAKey);
241 renounce(m_decC2Key);
242 renounce(m_decD2Key);
243 renounce(m_decE3Key);
244 renounce(m_decL2Key);
245 renounce(m_decL3Key);
246 }
247#endif
248
250 if ( m_calcSF ) {
251
253 m_weightConfig = std::make_unique<TFile>( m_weightConfigPath.c_str() );
254 if( !m_weightConfig ) {
255 ATH_MSG_ERROR( "SmoothedWZTagger: Error openning config file : " << m_weightConfigPath );
256 return StatusCode::FAILURE;
257 }
258
260 std::stringstream ss{m_weightFlavors};
261 std::string flavor;
262 while ( std::getline(ss, flavor, ',') ) {
263 m_weightHistograms.insert( std::make_pair( flavor, (TH2D*)m_weightConfig->Get((m_weightHistogramName+"_"+flavor).c_str()) ) );
264 if ( !m_efficiencyHistogramName.empty() ) {
265 m_efficiencyHistograms.insert( std::make_pair( flavor, (TH2D*)m_weightConfig->Get((m_efficiencyHistogramName+"_"+flavor).c_str()) ) );
266 }
267 ATH_MSG_INFO( "Tagging SF histogram for " << flavor << " is installed." );
268 }
269
270 }
271
273 if ( !m_strMassCutLow.empty() ) m_funcMassCutLow = std::make_unique<TF1>("strMassCutLow", m_strMassCutLow.c_str(), 0, 14000);
274 if ( !m_strMassCutHigh.empty() ) m_funcMassCutHigh = std::make_unique<TF1>("strMassCutHigh", m_strMassCutHigh.c_str(), 0, 14000);
275 if ( !m_strScoreCut.empty() ) m_funcScoreCut = std::make_unique<TF1>("strScoreCut", m_strScoreCut.c_str(), 0, 14000);
276
277 return StatusCode::SUCCESS;
278
279}
ss
static Double_t ss
Definition LArPhysWaveHECTool.cxx:37
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_funcMassCutLow
std::unique_ptr< TF1 > m_funcMassCutLow
TF1 for cut functions.
Definition JSSTaggerBase.h:184
JSSTaggerBase::m_decValidJetContentKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidJetContentKey
Definition JSSTaggerBase.h:74
JSSTaggerBase::m_decCutMHighKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutMHighKey
Definition JSSTaggerBase.h:219
JSSTaggerBase::m_weightConfig
std::unique_ptr< TFile > m_weightConfig
Histograms for scale factors.
Definition JSSTaggerBase.h:203
JSSTaggerBase::m_suppressOutputDependence
Gaudi::Property< bool > m_suppressOutputDependence
Definition JSSTaggerBase.h:176
JSSTaggerBase::m_truthLabelName
std::string m_truthLabelName
Definition JSSTaggerBase.h:212
JSSTaggerBase::m_suppressInputDependence
Gaudi::Property< bool > m_suppressInputDependence
Definition JSSTaggerBase.h:175
JSSTaggerBase::m_weightConfigPath
std::string m_weightConfigPath
Path to the SF configuration root file.
Definition JSSTaggerBase.h:107
JSSTaggerBase::m_readSphericityKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSphericityKey
Definition JSSTaggerBase.h:167
JSSTaggerBase::m_funcScoreCut
std::unique_ptr< TF1 > m_funcScoreCut
Definition JSSTaggerBase.h:186
JSSTaggerBase::m_readThrustMajKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readThrustMajKey
Definition JSSTaggerBase.h:166
JSSTaggerBase::m_decCutMLowKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutMLowKey
WriteDecorHandle keys for cut values.
Definition JSSTaggerBase.h:218
JSSTaggerBase::m_readQwKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readQwKey
Definition JSSTaggerBase.h:164
JSSTaggerBase::m_readSplit23Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSplit23Key
Definition JSSTaggerBase.h:162
JSSTaggerBase::m_decPassMassKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassMassKey
Definition JSSTaggerBase.h:77
JSSTaggerBase::m_nWarnVar
std::atomic< int > m_nWarnVar
Warning counters.
Definition JSSTaggerBase.h:84
JSSTaggerBase::m_decScoreValueKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decScoreValueKey
Definition JSSTaggerBase.h:221
JSSTaggerBase::m_decScoreCutKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decScoreCutKey
Definition JSSTaggerBase.h:220
JSSTaggerBase::m_decValidEventContentKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidEventContentKey
Definition JSSTaggerBase.h:75
JSSTaggerBase::m_funcMassCutHigh
std::unique_ptr< TF1 > m_funcMassCutHigh
Definition JSSTaggerBase.h:185
JSSTaggerBase::m_acceptInfo
asg::AcceptInfo m_acceptInfo
Object that stores the results for a jet.
Definition JSSTaggerBase.h:66
JSSTaggerBase::m_decTaggedKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTaggedKey
WriteDecorHandle keys for tagging bools.
Definition JSSTaggerBase.h:69
JSSTaggerBase::m_decPassScoreKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassScoreKey
Definition JSSTaggerBase.h:78
JSSTaggerBase::m_readSplit12Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSplit12Key
Definition JSSTaggerBase.h:161

◆ 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);};
IJetDecorator::decorate
virtual StatusCode decorate(const xAOD::JetContainer &jets) const =0
Decorate a jet collection without otherwise modifying it.

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ 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 }
AthCommonMsg< AlgTool >::msg
MsgStream & msg() const
Definition AthCommonMsg.h:24
MSG::name
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition MsgLevel.cxx:19

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

350 {
351
352 float scale = 1.0;
353 if ( m_ptGeV ) scale = 1.e3;
354
355 if ( jet.pt() < m_jetPtMin * scale ) return false;
356 if ( jet.pt() > m_jetPtMax * scale ) return false;
357 if ( std::abs( jet.eta() ) > m_jetEtaMax ) return false;
358
359 return true;
360
361}

◆ print()

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

Print the state of the tool.

Implements asg::IAsgTool.

Reimplemented in AsgHelloTool, HI::HIPileupTool, JetBottomUpSoftDrop, JetConstituentsRetriever, JetDumper, JetFinder, JetFromPseudojet, JetModifiedMassDrop, JetPileupLabelingTool, JetPruner, JetPseudojetRetriever, JetReclusterer, JetReclusteringTool, JetRecTool, JetRecursiveSoftDrop, JetSoftDrop, JetSplitter, JetSubStructureMomentToolsBase, JetToolRunner, JetTrimmer, JetTruthLabelingTool, and KtDeltaRTool.

Definition at line 131 of file AsgTool.cxx.

131 {
132
133 ATH_MSG_INFO( "AsgTool " << name() << " @ " << this );
134 return;
135 }

◆ printCuts()

void JSSTaggerBase::printCuts ( ) const
protected

Print configured cuts.

Definition at line 924 of file JSSTaggerBase.cxx.

924 {
925 ATH_MSG_INFO( "After tagging, you will have access to the following cuts as an asg::AcceptData : (<NCut>) <cut> : <description>)" );
926 int nCuts = m_acceptInfo.getNCuts();
927 for ( int iCut=0; iCut < nCuts; iCut++ ) {
928 std::string cut_string = "";
929 cut_string += " (";
930 cut_string += std::to_string(iCut);
931 cut_string += ") ";
932 cut_string += m_acceptInfo.getCutName(iCut).data();
933 cut_string += " : ";
934 cut_string += m_acceptInfo.getCutDescription(iCut).data();
935 ATH_MSG_INFO( cut_string );
936 }
937}

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

332 {
333
335 acceptData.clear();
336
338 acceptData.setCutResult( "ValidJetContent", true );
339 acceptData.setCutResult( "ValidEventContent", true );
340
341 acceptData.setCutResult( "ValidPtRangeHigh", true );
342 acceptData.setCutResult( "ValidPtRangeLow" , true );
343 acceptData.setCutResult( "ValidEtaRange" , true );
344
345 return StatusCode::SUCCESS;
346
347}
asg::AcceptData::clear
void clear()
Clear all bits.
Definition AcceptData.h:54

◆ 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 asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

◆ 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

Decorate single jet with tagging info.

Implemented in BJT::qgTagger, CP::JetQGTagger, CP::JetQGTaggerBDT, JSSTaggerUtils, JSSWTopTaggerANN, JSSWTopTaggerDNN, LundNetTagger, SmoothedTopTagger, and SmoothedWZTagger.

◆ 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 }
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition AthCommonDataStore.h:398

Member Data Documentation

◆ m_acceptInfo

asg::AcceptInfo JSSTaggerBase::m_acceptInfo
protected

Object that stores the results for a jet.

Definition at line 66 of file JSSTaggerBase.h.

◆ m_calcSF

bool JSSTaggerBase::m_calcSF {}
protected

Flag to calculate scale factor.

Definition at line 192 of file JSSTaggerBase.h.

192{};

◆ m_calibArea

std::string JSSTaggerBase::m_calibArea
protected

Location where config files live on cvmfs.

Definition at line 113 of file JSSTaggerBase.h.

◆ m_configFile

std::string JSSTaggerBase::m_configFile
protected

Configuration file name.

Definition at line 110 of file JSSTaggerBase.h.

◆ m_configReader

TEnv JSSTaggerBase::m_configReader
protected

TEnv instance to read config files.

Definition at line 63 of file JSSTaggerBase.h.

◆ m_containerName

std::string JSSTaggerBase::m_containerName
protected

Configurable members.

Jet container name

Definition at line 103 of file JSSTaggerBase.h.

◆ m_decC2Key

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

Definition at line 145 of file JSSTaggerBase.h.

145{this, "C2NameOut", "C2", "SG key for C2"};

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

219{this, "CutMHighName", "Cut_mhigh", "SG key for Cut_mhigh"};

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

218{this, "CutMLowName", "Cut_mlow", "SG key for Cut_mlow"};

◆ m_decD2Key

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

Definition at line 146 of file JSSTaggerBase.h.

146{this, "D2NameOut", "D2", "SG key for D2"};

◆ m_decE3Key

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

Definition at line 147 of file JSSTaggerBase.h.

147{this, "e3NameOut", "e3", "SG key for e3"};

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

225{this, "efficiencyName", "efficiency", "SG key for efficiency"};

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

226{this, "effSFName", "effSF", "SG key for effSF"};

◆ m_decL2Key

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

Definition at line 148 of file JSSTaggerBase.h.

148{this, "L2NameOut", "L2", "SG key for L2"};

◆ m_decL3Key

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

Definition at line 149 of file JSSTaggerBase.h.

149{this, "L3NameOut", "L3", "SG key for L3"};

◆ 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 77 of file JSSTaggerBase.h.

77{this, "PassMassName", "PassMass", "SG key for PassMass"};

◆ m_decPassScoreKey

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

Definition at line 78 of file JSSTaggerBase.h.

78{this, "PassScoreName", "PassScore", "SG key for PassScore"};

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

220{this, "CutScoreName", "Cut_score", "SG key for Cut_score"};

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

221{this, "ScoreName", "Score", "SG key for Score"};

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

227{this, "sigeffSFName", "sigeffSF", "SG key for effSF"};

◆ 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 69 of file JSSTaggerBase.h.

69{this, "TaggedName", "Tagged", "SG key for Tagged"};

◆ 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 141 of file JSSTaggerBase.h.

141{this, "Tau21WTANameOut", "Tau21_wta", "SG key for Tau21_wta"};

◆ m_decTau32WTAKey

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

Definition at line 142 of file JSSTaggerBase.h.

142{this, "Tau32WTANameOut", "Tau32_wta", "SG key for Tau32_wta"};

◆ m_decTau42WTAKey

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

Definition at line 143 of file JSSTaggerBase.h.

143{this, "Tau42WTANameOut", "Tau42_wta", "SG key for Tau42_wta"};

◆ m_decValidEtaRangeKey

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

Definition at line 72 of file JSSTaggerBase.h.

72{this, "ValidEtaRangeName", "ValidEtaRange", "SG key for ValidEtaRange"};

◆ m_decValidEventContentKey

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

Definition at line 75 of file JSSTaggerBase.h.

75{this, "ValidEventContentName", "ValidEventContent", "SG key for ValidEventContent"};

◆ m_decValidJetContentKey

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

Definition at line 74 of file JSSTaggerBase.h.

74{this, "ValidJetContentName", "ValidJetContent", "SG key for ValidJetContent"};

◆ m_decValidKinRangeKey

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

Definition at line 73 of file JSSTaggerBase.h.

73{this, "ValidKinRangeName", "ValidKinRange", "SG key for ValidKinRange"};

◆ m_decValidPtRangeHighKey

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

Definition at line 70 of file JSSTaggerBase.h.

70{this, "ValidPtRangeHighName", "ValidPtRangeHigh", "SG key for ValidPtRangeHigh"};

◆ m_decValidPtRangeLowKey

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

Definition at line 71 of file JSSTaggerBase.h.

71{this, "ValidPtRangeLowName", "ValidPtRangeLow", "SG key for ValidPtRangeLow"};

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

224{this, "weightName", "weight", "SG key for weight"};

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

193{};

◆ m_isSmallRJet

bool JSSTaggerBase::m_isSmallRJet = false
protected

Definition at line 104 of file JSSTaggerBase.h.

◆ m_jetEtaMax

float JSSTaggerBase::m_jetEtaMax {}
protected

Definition at line 133 of file JSSTaggerBase.h.

133{};

◆ m_jetPtMax

float JSSTaggerBase::m_jetPtMax {}
protected

Definition at line 132 of file JSSTaggerBase.h.

132{};

◆ m_jetPtMin

float JSSTaggerBase::m_jetPtMin {}
protected

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

Definition at line 131 of file JSSTaggerBase.h.

131{};

◆ m_kerasCalibArea

std::string JSSTaggerBase::m_kerasCalibArea
protected

Definition at line 119 of file JSSTaggerBase.h.

◆ m_kerasConfigFileName

std::string JSSTaggerBase::m_kerasConfigFileName
protected

Keras configurations for ML taggers.

Definition at line 116 of file JSSTaggerBase.h.

◆ m_kerasConfigFilePath

std::string JSSTaggerBase::m_kerasConfigFilePath
protected

Definition at line 117 of file JSSTaggerBase.h.

◆ m_kerasConfigOutputName

std::string JSSTaggerBase::m_kerasConfigOutputName
protected

Definition at line 118 of file JSSTaggerBase.h.

◆ m_nWarnMax

const int JSSTaggerBase::m_nWarnMax = 10
protected

Maximum number of warnings.

Definition at line 81 of file JSSTaggerBase.h.

◆ m_nWarnVar

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

Warning counters.

Definition at line 84 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 88 of file JSSTaggerBase.h.

◆ m_readECF1Key

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

Definition at line 157 of file JSSTaggerBase.h.

157{this, "ECF1Name", "ECF1", "SG key for ECF1"};

◆ m_readECF2Key

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

Definition at line 158 of file JSSTaggerBase.h.

158{this, "ECF2Name", "ECF2", "SG key for ECF2"};

◆ m_readECF3Key

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

Definition at line 159 of file JSSTaggerBase.h.

159{this, "ECF3Name", "ECF3", "SG key for ECF3"};

◆ 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 171 of file JSSTaggerBase.h.

171{this, "ECFG212Name", "ECFG_2_1_2", "SG key for ECFG_2_1_2"};

◆ 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 170 of file JSSTaggerBase.h.

170{this, "ECFG311Name", "ECFG_3_1_1", "SG key for ECFG_3_1_1"};

◆ 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 169 of file JSSTaggerBase.h.

169{this, "ECFG331Name", "ECFG_3_3_1", "SG key for ECFG_3_3_1"};

◆ m_readParentKey

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

Definition at line 173 of file JSSTaggerBase.h.

173{this, "ParentName", "Parent", "SG key for Parent"};

◆ m_readQwKey

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

Definition at line 164 of file JSSTaggerBase.h.

164{this, "QwName", "Qw", "SG key for Qw"};

◆ m_readSphericityKey

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

Definition at line 167 of file JSSTaggerBase.h.

167{this, "SphericityName", "Sphericity", "SG key for Sphericity"};

◆ m_readSplit12Key

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

Definition at line 161 of file JSSTaggerBase.h.

161{this, "Split12Name", "Split12", "SG key for Split12"};

◆ m_readSplit23Key

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

Definition at line 162 of file JSSTaggerBase.h.

162{this, "Split23Name", "Split23", "SG key for Split23"};

◆ 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 152 of file JSSTaggerBase.h.

152{this, "Tau1WTAName", "Tau1_wta", "SG key for Tau1_wta"};

◆ m_readTau2WTAKey

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

Definition at line 153 of file JSSTaggerBase.h.

153{this, "Tau2WTAName", "Tau2_wta", "SG key for Tau2_wta"};

◆ m_readTau3WTAKey

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

Definition at line 154 of file JSSTaggerBase.h.

154{this, "Tau3WTAName", "Tau3_wta", "SG key for Tau3_wta"};

◆ m_readTau4WTAKey

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

Definition at line 155 of file JSSTaggerBase.h.

155{this, "Tau4WTAName", "Tau4_wta", "SG key for Tau4_wta"};

◆ m_readThrustMajKey

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

Definition at line 166 of file JSSTaggerBase.h.

166{this, "ThrustMajName", "ThrustMaj", "SG key for ThrustMaj"};

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

215{this, "truthLabelName", "truthLabel", "SG key for truthLabel"};

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

Gaudi::Property<bool> JSSTaggerBase::m_suppressInputDependence {this, "SuppressInputDependence", false}
protected

Definition at line 175 of file JSSTaggerBase.h.

175{this, "SuppressInputDependence", false};

◆ m_suppressOutputDependence

Gaudi::Property<bool> JSSTaggerBase::m_suppressOutputDependence {this, "SuppressOutputDependence", false}
protected

Definition at line 176 of file JSSTaggerBase.h.

176{this, "SuppressOutputDependence", false};

◆ m_tagClass

TAGCLASS JSSTaggerBase::m_tagClass {Unknown}
protected

Definition at line 98 of file JSSTaggerBase.h.

98{Unknown};

◆ m_tagType

std::string JSSTaggerBase::m_tagType
protected

Definition at line 128 of file JSSTaggerBase.h.

◆ m_tmvaCalibArea

std::string JSSTaggerBase::m_tmvaCalibArea
protected

Definition at line 124 of file JSSTaggerBase.h.

◆ m_tmvaConfigFileName

std::string JSSTaggerBase::m_tmvaConfigFileName
protected

TMVA configurations for BDT taggers.

Definition at line 122 of file JSSTaggerBase.h.

◆ m_tmvaConfigFilePath

std::string JSSTaggerBase::m_tmvaConfigFilePath
protected

Definition at line 123 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.

208{};

◆ 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 137 of file JSSTaggerBase.h.

◆ m_useLSeries

bool JSSTaggerBase::m_useLSeries = false
protected

Definition at line 138 of file JSSTaggerBase.h.

◆ m_useMassCut

bool JSSTaggerBase::m_useMassCut = false
protected

Flag to indicate if mass window cut is used.

Definition at line 91 of file JSSTaggerBase.h.

◆ m_useScoreCut

bool JSSTaggerBase::m_useScoreCut = false
protected

Flag to indicate if a discriminant score is used.

Definition at line 94 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 107 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 127 of file JSSTaggerBase.h.

The documentation for this class was generated from the following files:
Generated on for ATLAS Offline Software by doxygen 1.14.0