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

#include <SmoothedTopTagger.h>

Inheritance diagram for SmoothedTopTagger:
Collaboration diagram for SmoothedTopTagger:

Public Member Functions

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

Protected Types

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

Protected Member Functions

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

Protected Attributes

TEnv m_configReader
 TEnv instance to read config files. More...
 
asg::AcceptInfo m_acceptInfo
 Object that stores the results for a jet. More...
 
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTaggedKey {this, "TaggedName", "Tagged", "SG key for Tagged"}
 WriteDecorHandle keys for tagging bools. More...
 
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"}
 
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readD2Key {this, "D2Name", "D2", "SG key for D2"}
 
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decNtrk500Key {this, "Ntrk500NameOut", "ParentJetNTrkPt500", "SG key for Ntrk500 from ungroomed jet"}
 
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readNtrk500Key {this, "Ntrk500Name", "ParentJetNTrkPt500", "SG key for Ntrk500 from ungroomed jet"}
 
const int m_nWarnMax = 10
 Maximum number of warnings. More...
 
std::atomic< int > m_nWarnVar
 Warning counters. More...
 
bool m_ptGeV = false
 Flag to indicate units pT is defined in Set to false by default. More...
 
bool m_useMassCut = false
 Flag to indicate if mass window cut is used. More...
 
bool m_useScoreCut = false
 Flag to indicate if a discriminant score is used. More...
 
TAGCLASS m_tagClass {Unknown}
 
std::string m_containerName
 Configurable members. More...
 
std::string m_weightConfigPath
 Path to the SF configuration root file. More...
 
std::string m_configFile
 Configuration file name. More...
 
std::string m_calibArea
 Location where config files live on cvmfs. More...
 
std::string m_kerasConfigFileName
 Keras configurations for ML taggers. More...
 
std::string m_kerasConfigFilePath
 
std::string m_kerasConfigOutputName
 
std::string m_kerasCalibArea
 
std::string m_tmvaConfigFileName
 TMVA configurations for BDT taggers. More...
 
std::string m_tmvaConfigFilePath
 
std::string m_tmvaCalibArea
 
std::string m_wkpt
 Tagger information. More...
 
std::string m_tagType
 
float m_jetPtMin {}
 Kinematic bounds for the jet - the units are controlled by m_ptGeV. More...
 
float m_jetPtMax {}
 
float m_jetEtaMax {}
 
bool m_useECFG = false
 Flags controlling whether generalized ECF moments or L-series ratios are needed TODO: Implement the functionality controlled by these. More...
 
bool m_useLSeries = false
 
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau21WTAKey {this, "Tau21WTANameOut", "Tau21_wta", "SG key for Tau21_wta"}
 WriteDecorHandle keys for JSS moments. More...
 
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. More...
 
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"}
 
bool m_suppressOutputDependence {}
 
std::string m_strMassCutLow
 Strings for cut functions. More...
 
std::string m_strMassCutHigh
 
std::string m_strScoreCut
 
std::unique_ptr< TF1 > m_funcMassCutLow
 TF1 for cut functions. More...
 
std::unique_ptr< TF1 > m_funcMassCutHigh
 
std::unique_ptr< TF1 > m_funcScoreCut
 
std::string m_decorationName
 Decoration name. More...
 
bool m_calcSF {}
 Flag to calculate scale factor. More...
 
bool m_isMC {}
 
std::string m_weightDecorationName
 String for scale factor decoration names. More...
 
std::string m_weightFileName
 
std::string m_weightHistogramName
 
std::string m_efficiencyHistogramName
 
std::string m_weightFlavors
 
std::unique_ptr< TFile > m_weightConfig
 Histograms for scale factors. More...
 
std::map< std::string, std::unique_ptr< TH2D > > m_weightHistograms
 
std::map< std::string, std::unique_ptr< TH2D > > m_efficiencyHistograms
 
bool m_truthLabelUseTRUTH3 {}
 Truth label options. More...
 
std::string m_truthParticleContainerName
 
std::string m_truthBosonContainerName
 
std::string m_truthTopQuarkContainerName
 
std::string m_truthLabelName
 
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTruthLabelKey {this, "truthLabelName", "truthLabel", "SG key for truthLabel"}
 Truth label ReadDecorHandle key. More...
 
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutMLowKey {this, "CutMLowName", "Cut_mlow", "SG key for Cut_mlow"}
 WriteDecorHandle keys for cut values. More...
 
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. More...
 
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 > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

int m_numTaggerVars {}
 
std::vector< std::string > m_varCutExprs
 
std::vector< std::string > m_varCutNames
 
std::vector< std::unique_ptr< TF1 > > m_varCutFuncs
 
SG::WriteDecorHandleKey< xAOD::JetContainer > m_dec_mcut {this, "mcutName", "Cut_m", "SG key for Cut_m"}
 
SG::WriteDecorHandleKey< xAOD::JetContainer > m_dec_sphericitycut {this, "sphericitycutName", "Cut_Sphericity", "SG key for Cut_Sphericity"}
 
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassSphericityKey {this, "PassSphericityName", "PassSphericity", "SG key for PassSphericity"}
 
std::vector< std::string > m_recognisedCuts = {"Mass", "mass", "Sphericity", "sphericity"}
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 10 of file SmoothedTopTagger.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ TAGCLASS

enum JSSTaggerBase::TAGCLASS
protectedinherited

TAGTYPE enum.

Enumerator
Unknown 
WBoson 
ZBoson 
TopQuark 

Definition at line 100 of file JSSTaggerBase.h.

100 {Unknown, WBoson, ZBoson, TopQuark};

Constructor & Destructor Documentation

◆ SmoothedTopTagger()

SmoothedTopTagger::SmoothedTopTagger ( const std::string &  name)

Definition at line 7 of file SmoothedTopTagger.cxx.

7  :
8  JSSTaggerBase( name )
9 {
10  // minimum and maximum pT of jets to tag
11  declareProperty( "JetPtMin", m_jetPtMin = 350.0e3);
12  declareProperty( "JetPtMax", m_jetPtMax = 3000.0e3);
13 
14  // cut functions that describe the tagger cuts that will be made
15  declareProperty( "VarCutFuncs", m_varCutExprs={}, "") ;
16  // names of the variables that are used for cuts
17  declareProperty( "VarCutNames", m_varCutNames={}, "") ;
18 
19 
20 }

Member Function Documentation

◆ calculateJSSRatios()

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

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

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

This should be changed in the future

Create write decor handles

Create read decor handles

WTA N-subjettiness ratios

ECF ratios

Definition at line 399 of file JSSTaggerBase.cxx.

399  {
400 
401  int result = 0;
402 
404  SG::WriteDecorHandle<xAOD::JetContainer, float> decTau21WTA(m_decTau21WTAKey);
405  SG::WriteDecorHandle<xAOD::JetContainer, float> decTau32WTA(m_decTau32WTAKey);
406  SG::WriteDecorHandle<xAOD::JetContainer, float> decTau42WTA(m_decTau42WTAKey);
407  SG::WriteDecorHandle<xAOD::JetContainer, float> decC2(m_decC2Key);
408  SG::WriteDecorHandle<xAOD::JetContainer, float> decD2(m_decD2Key);
409  SG::WriteDecorHandle<xAOD::JetContainer, float> decE3(m_decE3Key);
410  SG::WriteDecorHandle<xAOD::JetContainer, float> decL2(m_decL2Key);
411  SG::WriteDecorHandle<xAOD::JetContainer, float> decL3(m_decL3Key);
412 
414  SG::ReadDecorHandle<xAOD::JetContainer, float> readTau1WTA(m_readTau1WTAKey);
415  SG::ReadDecorHandle<xAOD::JetContainer, float> readTau2WTA(m_readTau2WTAKey);
416  SG::ReadDecorHandle<xAOD::JetContainer, float> readTau3WTA(m_readTau3WTAKey);
417  SG::ReadDecorHandle<xAOD::JetContainer, float> readTau4WTA(m_readTau4WTAKey);
418 
419  SG::ReadDecorHandle<xAOD::JetContainer, float> readECF1(m_readECF1Key);
420  SG::ReadDecorHandle<xAOD::JetContainer, float> readECF2(m_readECF2Key);
421  SG::ReadDecorHandle<xAOD::JetContainer, float> readECF3(m_readECF3Key);
422 
423  SG::ReadDecorHandle<xAOD::JetContainer, float> readECFG331(m_readECFG331Key);
424  SG::ReadDecorHandle<xAOD::JetContainer, float> readECFG311(m_readECFG311Key);
425  SG::ReadDecorHandle<xAOD::JetContainer, float> readECFG212(m_readECFG212Key);
426 
427 
429  float tau21_wta = -999.0;
430  float tau32_wta = -999.0;
431  float tau42_wta = -999.0;
432 
433  float tau1_wta = readTau1WTA(jet);
434  float tau2_wta = readTau2WTA(jet);
435  float tau3_wta = readTau3WTA(jet);
436  float tau4_wta = -999.0;
437  if(readTau4WTA.isAvailable()){
438  tau4_wta = readTau4WTA(jet);
439  }
440 
441  if ( tau1_wta > 1e-8 ) {
442  tau21_wta = tau2_wta / tau1_wta;
443  }
444  else result = 1;
445 
446  if ( tau2_wta > 1e-8 ) {
447  tau32_wta = tau3_wta / tau2_wta;
448  if(readTau4WTA.isAvailable()){
449  tau42_wta = tau4_wta / tau2_wta;
450  }
451  }
452  else result = 1;
453 
454  decTau21WTA(jet) = tau21_wta;
455  decTau32WTA(jet) = tau32_wta;
456  decTau42WTA(jet) = tau42_wta;
457 
459  float C2 = -999.0;
460  float D2 = -999.0;
461  float e3 = -999.0;
462 
463  float ECF1 = readECF1(jet);
464  float ECF2 = readECF2(jet);
465  float ECF3 = readECF3(jet);
466 
467  if ( ECF2 > 1e-8 ) {
468  C2 = ECF3 * ECF1 / std::pow( ECF2, 2.0 );
469  D2 = ECF3 * std::pow( ECF1, 3.0 ) / std::pow( ECF2, 3.0 );
470  }
471  else result = 1;
472 
473  e3 = ECF3 / std::pow( ECF1, 3.0 );
474 
475  decC2(jet) = C2;
476  decD2(jet) = D2;
477  decE3(jet) = e3;
478 
479  // L-series for UFO top taggers
480  float L2 = -999.0;
481  float L3 = -999.0;
482 
483  static const SG::AuxElement::ConstAccessor<float> accL2("L2");
484  if(!accL2.isAvailable(jet)){
485  if(readECFG331.isAvailable() && readECFG212.isAvailable()){
486  if(readECFG212(jet) > 1e-8){
487  L2 = readECFG331(jet) / pow(readECFG212(jet), (3.0/2.0));
488  }
489  else result = 1;
490  }
491  decL2(jet) = L2;
492  }
493 
494  static const SG::AuxElement::ConstAccessor<float> accL3("L3");
495  if(!accL3.isAvailable(jet)){
496  if(readECFG331.isAvailable() && readECFG311.isAvailable()){
497  if(readECFG331(jet) > 1e-8){
498  L3 = readECFG311(jet) / pow(readECFG331(jet), (1.0/3.0));
499  }
500  else result = 1;
501  }
502  decL3(jet) = L3;
503  }
504 
505  // TODO: Add ECFG for ANN tagger whenever it is defined
506 
507  return result;
508 
509 }

◆ checkKinRange()

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

Check and record if jet passes kinematic constraints.

Check each kinematic constraint Print warnings using counters

Create write decor handles

Decorate kinematic pass information

Definition at line 353 of file JSSTaggerBase.cxx.

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

◆ declareGaudiProperty() [1/4]

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

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

Definition at line 170 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [2/4]

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

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [3/4]

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

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

Definition at line 184 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [4/4]

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

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

Definition at line 199 of file AthCommonDataStore.h.

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

◆ declareProperty() [1/6]

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

Declare a new Gaudi property.

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

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

Definition at line 245 of file AthCommonDataStore.h.

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

◆ declareProperty() [2/6]

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

Declare a new Gaudi property.

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

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

Definition at line 221 of file AthCommonDataStore.h.

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

◆ declareProperty() [3/6]

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

Definition at line 259 of file AthCommonDataStore.h.

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

◆ declareProperty() [4/6]

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

Declare a new Gaudi property.

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

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

Definition at line 333 of file AthCommonDataStore.h.

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

◆ declareProperty() [5/6]

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

Declare a new Gaudi property.

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

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

Definition at line 352 of file AthCommonDataStore.h.

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

◆ declareProperty() [6/6]

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

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
147  return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
148  }

◆ decorate()

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

Decorate jet collection with tagging info.

Loop over jet collection and decorate each jet.

Implements IJetDecorator.

Reimplemented in JSSWTopTaggerANN, and JSSWTopTaggerDNN.

Definition at line 272 of file JSSTaggerBase.cxx.

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

◆ decorateJSSRatios()

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

Create write decor handles

Create read decor handles

WTA N-subjettiness ratios

ECF ratios

Definition at line 512 of file JSSTaggerBase.cxx.

512  {
513 
515  SG::WriteDecorHandle<xAOD::JetContainer, float> decTau21WTA(m_decTau21WTAKey);
516  bool decor_available = decTau21WTA.isAvailable();
517  if(decor_available) return;
518 
519  SG::WriteDecorHandle<xAOD::JetContainer, float> decTau32WTA(m_decTau32WTAKey);
520  SG::WriteDecorHandle<xAOD::JetContainer, float> decTau42WTA(m_decTau42WTAKey);
521  SG::WriteDecorHandle<xAOD::JetContainer, float> decC2(m_decC2Key);
522  SG::WriteDecorHandle<xAOD::JetContainer, float> decD2(m_decD2Key);
523  SG::WriteDecorHandle<xAOD::JetContainer, float> decE3(m_decE3Key);
524  SG::WriteDecorHandle<xAOD::JetContainer, float> decL2(m_decL2Key);
525  SG::WriteDecorHandle<xAOD::JetContainer, float> decL3(m_decL3Key);
526 
528  SG::ReadDecorHandle<xAOD::JetContainer, float> readTau1WTA(m_readTau1WTAKey);
529  SG::ReadDecorHandle<xAOD::JetContainer, float> readTau2WTA(m_readTau2WTAKey);
530  SG::ReadDecorHandle<xAOD::JetContainer, float> readTau3WTA(m_readTau3WTAKey);
531  SG::ReadDecorHandle<xAOD::JetContainer, float> readTau4WTA(m_readTau4WTAKey);
532 
533  SG::ReadDecorHandle<xAOD::JetContainer, float> readECF1(m_readECF1Key);
534  SG::ReadDecorHandle<xAOD::JetContainer, float> readECF2(m_readECF2Key);
535  SG::ReadDecorHandle<xAOD::JetContainer, float> readECF3(m_readECF3Key);
536 
537  SG::ReadDecorHandle<xAOD::JetContainer, float> readECFG331(m_readECFG331Key);
538  SG::ReadDecorHandle<xAOD::JetContainer, float> readECFG311(m_readECFG311Key);
539  SG::ReadDecorHandle<xAOD::JetContainer, float> readECFG212(m_readECFG212Key);
540 
541  for(const xAOD::Jet* jet : jets){
542 
544  float tau21_wta = -999.0;
545  float tau32_wta = -999.0;
546  float tau42_wta = -999.0;
547 
548  float tau1_wta = readTau1WTA(*jet);
549  float tau2_wta = readTau2WTA(*jet);
550  float tau3_wta = readTau3WTA(*jet);
551  float tau4_wta = -999.0;
552  if(readTau4WTA.isAvailable()){
553  tau4_wta = readTau4WTA(*jet);
554  }
555 
556  if ( tau1_wta > 1e-8 ) {
557  tau21_wta = tau2_wta / tau1_wta;
558  }
559 
560  if ( tau2_wta > 1e-8 ) {
561  tau32_wta = tau3_wta / tau2_wta;
562  if(readTau4WTA.isAvailable()){
563  tau42_wta = tau4_wta / tau2_wta;
564  }
565  }
566 
567  decTau21WTA(*jet) = tau21_wta;
568  decTau32WTA(*jet) = tau32_wta;
569  decTau42WTA(*jet) = tau42_wta;
570 
572  float C2 = -999.0;
573  float D2 = -999.0;
574  float e3 = -999.0;
575 
576  float ECF1 = readECF1(*jet);
577  float ECF2 = readECF2(*jet);
578  float ECF3 = readECF3(*jet);
579 
580  if ( ECF2 > 1e-8 ) {
581  C2 = ECF3 * ECF1 / std::pow( ECF2, 2.0 );
582  D2 = ECF3 * std::pow( ECF1, 3.0 ) / std::pow( ECF2, 3.0 );
583  }
584 
585  e3 = ECF3 / std::pow( ECF1, 3.0 );
586 
587  decC2(*jet) = C2;
588  decD2(*jet) = D2;
589  decE3(*jet) = e3;
590 
591  // L-series for UFO top taggers
592  float L2 = -999.0;
593  float L3 = -999.0;
594 
595  static const SG::AuxElement::ConstAccessor<float> accL2("L2");
596  if(!accL2.isAvailable(*jet)){
597  if(readECFG331.isAvailable() && readECFG212.isAvailable()){
598  if(readECFG212(*jet) > 1e-8){
599  L2 = readECFG331(*jet) / pow(readECFG212(*jet), (3.0/2.0));
600  }
601  }
602  decL2(*jet) = L2;
603  }
604 
605  static const SG::AuxElement::ConstAccessor<float> accL3("L3");
606  if(!accL3.isAvailable(*jet)){
607  if(readECFG331.isAvailable() && readECFG311.isAvailable()){
608  if(readECFG331(*jet) > 1e-8){
609  L3 = readECFG311(*jet) / pow(readECFG331(*jet), (1.0/3.0));
610  }
611  }
612  decL3(*jet) = L3;
613  }
614 
615  // TODO: Add ECFG for ANN tagger whenever it is defined
616 
617  }
618 
619  return;
620 
621 }

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

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

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

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

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

Add StoreName to extra input/output deps as needed.

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

◆ findPV()

int JSSTaggerBase::findPV ( ) const
protectedinherited

Find the PV (to be used for Ntrk)

Definition at line 624 of file JSSTaggerBase.cxx.

624  {
625 
626  int indexPV = -1;
627 
628  const xAOD::VertexContainer* vxCont = nullptr;
629  if ( evtStore()->retrieve( vxCont, "PrimaryVertices" ).isFailure() ) {
630  ATH_MSG_WARNING( "Unable to retrieve primary vertex container PrimaryVertices" );
631  }
632  else {
633  int vtx_counter = 0;
634  for ( const auto *vx : *vxCont ) {
635  if ( vx->vertexType()==xAOD::VxType::PriVtx ) {
636  indexPV = vtx_counter;
637  break;
638  }
639  vtx_counter++;
640  }
641  }
642 
643  return indexPV;
644 
645 }

◆ getConfigReader()

StatusCode JSSTaggerBase::getConfigReader ( )
protectedinherited

Get configReader StatusCode.

Get configReader TEnv.

Check for the existence of the configuration file

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

Definition at line 283 of file JSSTaggerBase.cxx.

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

◆ getKey()

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

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

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

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

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

Definition at line 119 of file AsgTool.cxx.

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

◆ getName()

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

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

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

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

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

Definition at line 106 of file AsgTool.cxx.

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

◆ getProperty()

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

Get one of the tool's properties.

◆ getSF()

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

Get scale factor and efficiency.

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

Definition at line 760 of file JSSTaggerBase.cxx.

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

◆ getTruthLabelStr()

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

Truth label string

Truth label value

Contained top tagger

Contained top

QCD

TCC W/Z 2-var tagger

Top

W/Z

QCD

Pass mass and D2

Fail mass, pass D2

Pass mass, fail D2

Fail mass and D2

W/Z tagger

Top

W/Z

QCD

Top

QCD

Definition at line 818 of file JSSTaggerBase.cxx.

818  {
820  std::string truthLabelStr;
821 
823  SG::ReadDecorHandle<xAOD::JetContainer, int> readTruthLabel(m_readTruthLabelKey);
824  LargeRJetTruthLabel::TypeEnum jetContainment = LargeRJetTruthLabel::intToEnum(readTruthLabel(jet));
825 
827  if ( m_weightHistograms.count("t_qqb") ) {
828 
830  if ( jetContainment==LargeRJetTruthLabel::tqqb ) {
831  truthLabelStr = "t_qqb";
832  }
834  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
835  truthLabelStr = "q";
836  }
837 
838  }
840  else if ( m_weightHistograms.count("V_qq_passMpassD2") ) {
841 
843  if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t ) {
844  truthLabelStr = "t_";
845  }
847  else if ( jetContainment==LargeRJetTruthLabel::Wqq || jetContainment==LargeRJetTruthLabel::Zqq || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
848  truthLabelStr = "V_qq_";
849  }
851  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
852  truthLabelStr = "q_";
853  }
854 
856  if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") && acceptData.getCutResult("PassD2") ) {
857  truthLabelStr += "passMpassD2";
858  }
860  else if ( !(acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh")) && acceptData.getCutResult("PassD2") ) {
861  truthLabelStr += "failMpassD2";
862  }
864  else if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") && !acceptData.getCutResult("PassD2") ) {
865  truthLabelStr += "passMfailD2";
866  }
868  else{
869  truthLabelStr += "failMfailD2";
870  }
871 
872  }
873 
875  else if ( m_weightHistograms.count("V_qq") ) {
876 
878  if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t ) {
879  truthLabelStr = "t";
880  }
882  else if ( jetContainment==LargeRJetTruthLabel::Wqq || jetContainment==LargeRJetTruthLabel::Zqq || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
883  truthLabelStr = "V_qq";
884  }
886  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
887  truthLabelStr = "q";
888  }
889  }
890 
891  // inclusive top tagger
892  else {
894  if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
895  truthLabelStr = "t";
896  }
898  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
899  truthLabelStr = "q";
900  }
901 
902  }
903 
904  return truthLabelStr;
905 }

◆ GetUnGroomTracks()

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

Retrieve Ntrk variable from the ungroomed parent jet.

Retrieve the Ntrk variable from the ungroomed parent jet.

Definition at line 648 of file JSSTaggerBase.cxx.

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

◆ getWeight()

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

Get SF weight.

Inefficiency SF is directly used

Efficiency SF

Calculate inefficiency SF

If inefficiency SF is not available, SF is always 1.0

Create write decor handles

Decorate values

Definition at line 686 of file JSSTaggerBase.cxx.

686  {
687  if ( !m_calcSF ) return StatusCode::SUCCESS;
688 
689  float weight = 1.0;
690  float effSF = 1.0;
691  float sigeffSF = 1.0;
692  float efficiency = 1.0;
693 
694  if ( m_isMC ) {
695 
696  std::string truthLabelStr = getTruthLabelStr( jet, acceptData );
697  std::tie(effSF, efficiency) = getSF( jet, truthLabelStr );
698 
699  // calculate signal efficiency SF
700  if ( m_weightHistograms.count("t_qqb") ) {
701  sigeffSF = getSF(jet, "t_qqb").first;
702  } else if ( m_weightHistograms.count("V_qq") ) {
703  sigeffSF = getSF(jet, "V_qq").first;
704  } else if ( m_weightHistograms.count("t") ){
705  sigeffSF = getSF(jet, "t").first;
706  } else {
707  sigeffSF = 1.0;
708  }
709 
711  if ( m_weightFlavors.find("fail") != std::string::npos ) {
712  weight = effSF;
713  }
714 
715  else {
716 
718  if ( passSel ) {
719  weight = effSF;
720  }
721 
723  else {
725  if ( m_efficiencyHistogramName.empty() ) {
726  weight = 1.0;
727  }
728  else if ( efficiency < 1.0 ) {
729  weight = ( 1. - effSF * efficiency ) / ( 1. - efficiency );
730  }
731  else {
732  weight = 1.0;
733  }
734  }
735  }
736 
737  }
738 
739  else {
740  weight = 1.0;
741  }
742 
744  SG::WriteDecorHandle<xAOD::JetContainer, float> decWeight(m_decWeightKey);
745  SG::WriteDecorHandle<xAOD::JetContainer, float> decEfficiency(m_decEfficiencyKey);
746  SG::WriteDecorHandle<xAOD::JetContainer, float> decEffSF(m_decEffSFKey);
747  SG::WriteDecorHandle<xAOD::JetContainer, float> decSigeffSF(m_decSigeffSFKey);
748 
750  decWeight(jet) = weight;
751  decEfficiency(jet) = efficiency;
752  decEffSF(jet) = effSF;
753  decSigeffSF(jet) = sigeffSF;
754 
755  return StatusCode::SUCCESS;
756 
757 }

◆ initialize()

StatusCode SmoothedTopTagger::initialize ( )
overridevirtual

Initialize the tool.

Get configReader

Call base class initialize

Loop over and print out the cuts that have been configured

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

Definition at line 22 of file SmoothedTopTagger.cxx.

22  {
23 
24  ATH_MSG_INFO( "Initializing SmoothedTopTagger tool" );
25 
26  if( ! m_configFile.empty() ) {
27 
29  ATH_CHECK( getConfigReader() );
30 
31  TString prefix = "";
32  if ( ! m_wkpt.empty() ) prefix = m_wkpt+".";
33 
34  // read the number of variables for tagger from file
35  m_numTaggerVars = std::stoi(m_configReader.GetValue( prefix+"NumVars", ""));
36 
37  ATH_MSG_DEBUG("Number of variables used by tagger is " << std::to_string(m_numTaggerVars));
38 
39  std::string varName, varCutExpr;
40  for (int i = 1; i <= m_numTaggerVars; i++) {
41  // read the cut name corresponding to this variable
42  varName = m_configReader.GetValue( prefix+"Var"+std::to_string(i), "");
43 
44  if (varName.empty()) {
45  ATH_MSG_ERROR("Config file does not specify Var" << std::to_string(i) << "!") ;
46  return StatusCode::FAILURE;
47  }
48 
49  if (std::find(m_recognisedCuts.begin(), m_recognisedCuts.end(), varName) == m_recognisedCuts.end()) {
50  ATH_MSG_ERROR("Unrecognised variable " << varName << " in config file!") ;
51  return StatusCode::FAILURE;
52  }
53 
54  m_varCutNames.push_back(varName);
55 
56  // read cut expression
57  varCutExpr = m_configReader.GetValue( prefix+m_varCutNames.back()+"Cut", "");
58 
59  if (varCutExpr.empty()) {
60  ATH_MSG_ERROR("Config file does not specify Var" << std::to_string(i) << " cut!") ;
61  return StatusCode::FAILURE;
62  }
63 
64  m_varCutExprs.push_back(varCutExpr);
65 
66  ATH_MSG_DEBUG("Read Var" << std::to_string(i) << " from file:");
67  ATH_MSG_DEBUG(" Variable: " << m_varCutNames.back());
68  ATH_MSG_DEBUG(" Cut: " << m_varCutExprs.back());
69  }
70 
71  // get the decoration name
72  m_decorationName = m_configReader.GetValue("DecorationName" ,"");
73  } else { // no config file
74 
75  // determine number of tagger variables from size of vectors
76  m_numTaggerVars = int(m_varCutNames.size());
77 
78  // make sure all of the tagger variables are recognised
79  for (const std::string& var: m_varCutNames) {
80  if (std::find(m_recognisedCuts.begin(), m_recognisedCuts.end(), var) == m_recognisedCuts.end()) {
81  ATH_MSG_ERROR("Unrecognised variable " << var << " provided!") ;
82  return StatusCode::FAILURE;
83  }
84  }
85  } // if config file
86 
87  // make sure cut names vector is not empty
88  if (m_varCutNames.empty()) {
89  ATH_MSG_ERROR( "Tagger variable names vector is empty." ) ;
90  return StatusCode::FAILURE;
91  }
92 
93  // make sure cut expressions vector is not empty
94  if (m_varCutExprs.empty()) {
95  ATH_MSG_ERROR( "Tagger variable cuts vector is empty." ) ;
96  return StatusCode::FAILURE;
97  }
98 
99  // make sure all vectors have the same length
100  if (int(m_varCutNames.size()) != int(m_varCutExprs.size())) {
101  ATH_MSG_ERROR( "Tagger variable names and cuts vectors don't have the same size." ) ;
102  return StatusCode::FAILURE;
103  }
104 
105  // make cut functions into TFI objects
106  for (int i = 0; i < m_numTaggerVars; i++) {
107  m_varCutFuncs.push_back(
108  std::make_unique<TF1>(m_varCutNames[i].c_str(), m_varCutExprs[i].c_str())
109  );
110  ATH_MSG_DEBUG("Configured " << m_varCutNames[i] << " : " << m_varCutExprs[i] << " cut TF1");
111  }
112 
113  ATH_MSG_INFO( "Smoothed top Tagger tool initialized with cuts:" );
114  for (int i = 0; i < m_numTaggerVars; i++) {
115  ATH_MSG_INFO( m_varCutNames[i]+" cut : "<< m_varCutExprs[i] );
116  }
117  ATH_MSG_INFO( "DecorationName : "<< m_decorationName );
118 
119  // add cuts for the output TAccept
120  // initialize decorators as decorationName+_decorator
121  ATH_MSG_INFO( "Additional decorators that will be attached to jet :" );
122 
123  if (std::find(m_varCutNames.begin(), m_varCutNames.end(), "Mass") != m_varCutNames.end() || std::find(m_varCutNames.begin(), m_varCutNames.end(), "mass") != m_varCutNames.end()) {
124  m_acceptInfo.addCut( "PassMass" , "mJet > mCut" );
125  // initialize decorators for passing cuts
126  // this uses m_decPassMassKey inherited from JSSTaggerBase
127  m_decPassMassKey = m_containerName + "." + m_decorationName + "_" + m_decPassMassKey.key();
128  ATH_CHECK( m_decPassMassKey.initialize() );
129  m_dec_mcut = m_containerName + "." + m_decorationName + "_" + m_dec_mcut.key();
130  ATH_CHECK( m_dec_mcut.initialize() );
131  ATH_MSG_INFO( " " << m_dec_mcut.key() << " : mass cut" );
132  ATH_MSG_INFO( " " << m_decPassMassKey.key() << " : pass mass cut");
133  }
134 
135  if (std::find(m_varCutNames.begin(), m_varCutNames.end(), "Sphericity") != m_varCutNames.end() || std::find(m_varCutNames.begin(), m_varCutNames.end(), "sphericity") != m_varCutNames.end()) {
136  m_acceptInfo.addCut( "PassSphericity" , "SphericityJet > SphericityCut" );
137  m_decPassSphericityKey = m_containerName + "." + m_decorationName + "_" + m_decPassSphericityKey.key();
138  ATH_CHECK( m_decPassSphericityKey.initialize() );
139  m_dec_sphericitycut = m_containerName + "." + m_decorationName + "_" + m_dec_sphericitycut.key();
140  ATH_CHECK( m_dec_sphericitycut.initialize() );
141  ATH_MSG_INFO( " " << m_dec_sphericitycut.key() << " : Sphericity cut" );
142  ATH_MSG_INFO( " " << m_decPassSphericityKey.key() << " : pass Sphericity cut" );
143  }
144 
146  ATH_CHECK( JSSTaggerBase::initialize() );
147 
149  printCuts();
150 
151  return StatusCode::SUCCESS;
152 } // end initialize()

◆ inputHandles()

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

Return this algorithm's input handles.

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

◆ modify()

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

Concrete implementation of the function inherited from IJetModifier.

Implements IJetModifier.

Definition at line 32 of file IJetDecorator.h.

32 {return decorate(jets);};

◆ msg() [1/2]

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

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

Definition at line 27 of file AthCommonMsg.h.

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

◆ msg_level_name()

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

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

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

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

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

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

Definition at line 101 of file AsgTool.cxx.

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

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

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

◆ outputHandles()

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

Return this algorithm's output handles.

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

◆ passKinRange()

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

Check if jet passes kinematic constraints.

Definition at line 339 of file JSSTaggerBase.cxx.

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

◆ print()

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

Print the state of the tool.

Implements asg::IAsgTool.

Reimplemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, 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
protectedinherited

Print configured cuts.

Definition at line 907 of file JSSTaggerBase.cxx.

907  {
908  ATH_MSG_INFO( "After tagging, you will have access to the following cuts as an asg::AcceptData : (<NCut>) <cut> : <description>)" );
909  int nCuts = m_acceptInfo.getNCuts();
910  for ( int iCut=0; iCut < nCuts; iCut++ ) {
911  std::string cut_string = "";
912  cut_string += " (";
913  cut_string += std::to_string(iCut);
914  cut_string += ") ";
915  cut_string += m_acceptInfo.getCutName(iCut).data();
916  cut_string += " : ";
917  cut_string += m_acceptInfo.getCutDescription(iCut).data();
918  ATH_MSG_INFO( cut_string );
919  }
920 }

◆ renounce()

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

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ resetCuts()

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

Reset cuts.

Reset the AcceptData cut results to false

Initialize common cuts to true by default

Definition at line 321 of file JSSTaggerBase.cxx.

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

◆ sysInitialize()

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

Perform system initialization for an algorithm.

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

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

◆ sysStart()

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

Handle START transition.

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

◆ tag()

StatusCode SmoothedTopTagger::tag ( const xAOD::Jet &  jet) const
overridevirtual

Decorate single jet with tagging info.

Create asg::AcceptData object

Reset the AcceptData cut results

Check basic kinematic selection

Calculate NSubjettiness and ECF ratios

Implements JSSTaggerBase.

Definition at line 154 of file SmoothedTopTagger.cxx.

154  {
155 
156  ATH_MSG_DEBUG( "Obtaining smoothed top result" );
157 
159  asg::AcceptData acceptData( &m_acceptInfo );
160 
162  ATH_CHECK( resetCuts( acceptData ) );
163 
165  ATH_CHECK( checkKinRange( jet, acceptData ) );
166 
167  // get the relevant attributes of the jet
168  // mass and pt - note that this will depend on the configuration of the calibration used
169  float jet_pt = jet.pt()/1000.0;
170  float jet_mass = jet.m()/1000.0;
171 
173  calculateJSSRatios(jet);
174 
175  // configure decorators from JSSTaggerBase class
176  SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidJetContent(m_decValidJetContentKey);
177  SG::WriteDecorHandle<xAOD::JetContainer, bool> decTagged(m_decTaggedKey);
178 
179  // initialize for use in other statements
180  bool passCuts = true;
181  float cut_var;
182  for (int i = 0; i < m_numTaggerVars; i++) {
183  // evaluate the cut value on this variable
184  cut_var = m_varCutFuncs[i]->Eval(jet_pt);
185 
186  // check which variable this cut corresponds to and make the
187  // selection
188  // when more taggers are implemented add the variables required
189  // into this if-else ladder
190  if (m_varCutNames[i] == "Mass" || m_varCutNames[i] == "mass") {
191  // decorators for jet after applying cuts
192  SG::WriteDecorHandle<xAOD::JetContainer, float> decMCut(m_dec_mcut);
193  SG::WriteDecorHandle<xAOD::JetContainer, bool> decPassMass(m_decPassMassKey);
194 
195  // decorate cut
196  decMCut(jet) = cut_var;
197 
198  // make cut and decorate on jet
199  if(jet_mass > cut_var) {
200  acceptData.setCutResult("PassMass",true);
201  }
202  decPassMass(jet) = acceptData.getCutResult("PassMass");
203  passCuts = passCuts && acceptData.getCutResult("PassMass");
204  }
205  else if (m_varCutNames[i] == "Sphericity" || m_varCutNames[i] == "sphericity") {
206  float sphericity = 0;
207 
208  // setup read/write handles for sphericity cut decorations & reading variables
209  SG::ReadDecorHandle<xAOD::JetContainer, float> readSphericity(m_readSphericityKey);
210  SG::WriteDecorHandle<xAOD::JetContainer, float> decSphericityCut(m_dec_sphericitycut);
211  SG::WriteDecorHandle<xAOD::JetContainer, bool> decPassSphericity(m_decPassSphericityKey);
212 
213  // decorate cut
214  decSphericityCut(jet) = cut_var;
215 
216  // read sphericity variable
217  if ( !readSphericity.isAvailable() ) {
218  ATH_MSG_VERBOSE( "The Sphericity variable is not available in your file" );
219  acceptData.setCutResult("ValidJetContent", false);
220  decValidJetContent(jet) = false;
221  } else {
222  // get sphericity only if it is decorated
223  sphericity = readSphericity(jet);
224  }
225 
226  // make cut and decorate results
227  if (sphericity > cut_var) {
228  acceptData.setCutResult("PassSphericity", true);
229  }
230  decPassSphericity(jet) = acceptData.getCutResult("PassSphericity");
231  // incorporate cut into complete set of selections
232  passCuts = passCuts && acceptData.getCutResult("PassSphericity");
233  }
234  }
235 
236  // decorate jets with tagging information and whether content is valid
237  decValidJetContent(jet) = acceptData.getCutResult("ValidJetContent");
238  decTagged(jet) = passCuts;
239 
240  return StatusCode::SUCCESS;
241 }

◆ updateVHKA()

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

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_acceptInfo

asg::AcceptInfo JSSTaggerBase::m_acceptInfo
protectedinherited

Object that stores the results for a jet.

Definition at line 65 of file JSSTaggerBase.h.

◆ m_calcSF

bool JSSTaggerBase::m_calcSF {}
protectedinherited

Flag to calculate scale factor.

Definition at line 193 of file JSSTaggerBase.h.

◆ m_calibArea

std::string JSSTaggerBase::m_calibArea
protectedinherited

Location where config files live on cvmfs.

Definition at line 115 of file JSSTaggerBase.h.

◆ m_configFile

std::string JSSTaggerBase::m_configFile
protectedinherited

Configuration file name.

Definition at line 112 of file JSSTaggerBase.h.

◆ m_configReader

TEnv JSSTaggerBase::m_configReader
protectedinherited

TEnv instance to read config files.

Definition at line 62 of file JSSTaggerBase.h.

◆ m_containerName

std::string JSSTaggerBase::m_containerName
protectedinherited

Configurable members.

Jet container name

Definition at line 106 of file JSSTaggerBase.h.

◆ m_dec_mcut

SG::WriteDecorHandleKey<xAOD::JetContainer> SmoothedTopTagger::m_dec_mcut {this, "mcutName", "Cut_m", "SG key for Cut_m"}
private

Definition at line 40 of file SmoothedTopTagger.h.

◆ m_dec_sphericitycut

SG::WriteDecorHandleKey<xAOD::JetContainer> SmoothedTopTagger::m_dec_sphericitycut {this, "sphericitycutName", "Cut_Sphericity", "SG key for Cut_Sphericity"}
private

Definition at line 41 of file SmoothedTopTagger.h.

◆ m_decC2Key

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

Definition at line 147 of file JSSTaggerBase.h.

◆ m_decCutMHighKey

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

Definition at line 220 of file JSSTaggerBase.h.

◆ m_decCutMLowKey

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

WriteDecorHandle keys for cut values.

Definition at line 219 of file JSSTaggerBase.h.

◆ m_decD2Key

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

Definition at line 148 of file JSSTaggerBase.h.

◆ m_decE3Key

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

Definition at line 149 of file JSSTaggerBase.h.

◆ m_decEfficiencyKey

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

Definition at line 226 of file JSSTaggerBase.h.

◆ m_decEffSFKey

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

Definition at line 227 of file JSSTaggerBase.h.

◆ m_decL2Key

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

Definition at line 150 of file JSSTaggerBase.h.

◆ m_decL3Key

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

Definition at line 151 of file JSSTaggerBase.h.

◆ m_decNtrk500Key

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

Definition at line 80 of file JSSTaggerBase.h.

◆ m_decorationName

std::string JSSTaggerBase::m_decorationName
protectedinherited

Decoration name.

Definition at line 190 of file JSSTaggerBase.h.

◆ m_decPassMassKey

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

Definition at line 76 of file JSSTaggerBase.h.

◆ m_decPassScoreKey

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

Definition at line 77 of file JSSTaggerBase.h.

◆ m_decPassSphericityKey

SG::WriteDecorHandleKey<xAOD::JetContainer> SmoothedTopTagger::m_decPassSphericityKey {this, "PassSphericityName", "PassSphericity", "SG key for PassSphericity"}
private

Definition at line 44 of file SmoothedTopTagger.h.

◆ m_decScoreCutKey

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

Definition at line 221 of file JSSTaggerBase.h.

◆ m_decScoreValueKey

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

Definition at line 222 of file JSSTaggerBase.h.

◆ m_decSigeffSFKey

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

Definition at line 228 of file JSSTaggerBase.h.

◆ m_decTaggedKey

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

WriteDecorHandle keys for tagging bools.

Definition at line 68 of file JSSTaggerBase.h.

◆ m_decTau21WTAKey

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

WriteDecorHandle keys for JSS moments.

Definition at line 143 of file JSSTaggerBase.h.

◆ m_decTau32WTAKey

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

Definition at line 144 of file JSSTaggerBase.h.

◆ m_decTau42WTAKey

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

Definition at line 145 of file JSSTaggerBase.h.

◆ m_decValidEtaRangeKey

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

Definition at line 71 of file JSSTaggerBase.h.

◆ m_decValidEventContentKey

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

Definition at line 74 of file JSSTaggerBase.h.

◆ m_decValidJetContentKey

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

Definition at line 73 of file JSSTaggerBase.h.

◆ m_decValidKinRangeKey

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

Definition at line 72 of file JSSTaggerBase.h.

◆ m_decValidPtRangeHighKey

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

Definition at line 69 of file JSSTaggerBase.h.

◆ m_decValidPtRangeLowKey

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

Definition at line 70 of file JSSTaggerBase.h.

◆ m_decWeightKey

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

WriteDecorHandle keys for SF.

Definition at line 225 of file JSSTaggerBase.h.

◆ m_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_efficiencyHistogramName

std::string JSSTaggerBase::m_efficiencyHistogramName
protectedinherited

Definition at line 200 of file JSSTaggerBase.h.

◆ m_efficiencyHistograms

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

Definition at line 206 of file JSSTaggerBase.h.

◆ m_evtStore

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

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_funcMassCutHigh

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

Definition at line 186 of file JSSTaggerBase.h.

◆ m_funcMassCutLow

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

TF1 for cut functions.

Definition at line 185 of file JSSTaggerBase.h.

◆ m_funcScoreCut

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

Definition at line 187 of file JSSTaggerBase.h.

◆ m_isMC

bool JSSTaggerBase::m_isMC {}
protectedinherited

Definition at line 194 of file JSSTaggerBase.h.

◆ m_jetEtaMax

float JSSTaggerBase::m_jetEtaMax {}
protectedinherited

Definition at line 135 of file JSSTaggerBase.h.

◆ m_jetPtMax

float JSSTaggerBase::m_jetPtMax {}
protectedinherited

Definition at line 134 of file JSSTaggerBase.h.

◆ m_jetPtMin

float JSSTaggerBase::m_jetPtMin {}
protectedinherited

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

Definition at line 133 of file JSSTaggerBase.h.

◆ m_kerasCalibArea

std::string JSSTaggerBase::m_kerasCalibArea
protectedinherited

Definition at line 121 of file JSSTaggerBase.h.

◆ m_kerasConfigFileName

std::string JSSTaggerBase::m_kerasConfigFileName
protectedinherited

Keras configurations for ML taggers.

Definition at line 118 of file JSSTaggerBase.h.

◆ m_kerasConfigFilePath

std::string JSSTaggerBase::m_kerasConfigFilePath
protectedinherited

Definition at line 119 of file JSSTaggerBase.h.

◆ m_kerasConfigOutputName

std::string JSSTaggerBase::m_kerasConfigOutputName
protectedinherited

Definition at line 120 of file JSSTaggerBase.h.

◆ m_numTaggerVars

int SmoothedTopTagger::m_numTaggerVars {}
private

Definition at line 28 of file SmoothedTopTagger.h.

◆ m_nWarnMax

const int JSSTaggerBase::m_nWarnMax = 10
protectedinherited

Maximum number of warnings.

Definition at line 84 of file JSSTaggerBase.h.

◆ m_nWarnVar

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

Warning counters.

Definition at line 87 of file JSSTaggerBase.h.

◆ m_ptGeV

bool JSSTaggerBase::m_ptGeV = false
protectedinherited

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

Definition at line 91 of file JSSTaggerBase.h.

◆ m_readD2Key

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

Definition at line 79 of file JSSTaggerBase.h.

◆ m_readECF1Key

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

Definition at line 159 of file JSSTaggerBase.h.

◆ m_readECF2Key

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

Definition at line 160 of file JSSTaggerBase.h.

◆ m_readECF3Key

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

Definition at line 161 of file JSSTaggerBase.h.

◆ m_readECFG212Key

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

Definition at line 173 of file JSSTaggerBase.h.

◆ m_readECFG311Key

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

Definition at line 172 of file JSSTaggerBase.h.

◆ m_readECFG331Key

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

Definition at line 171 of file JSSTaggerBase.h.

◆ m_readNtrk500Key

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

Definition at line 81 of file JSSTaggerBase.h.

◆ m_readParentKey

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

Definition at line 175 of file JSSTaggerBase.h.

◆ m_readQwKey

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

Definition at line 166 of file JSSTaggerBase.h.

◆ m_readSphericityKey

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

Definition at line 169 of file JSSTaggerBase.h.

◆ m_readSplit12Key

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

Definition at line 163 of file JSSTaggerBase.h.

◆ m_readSplit23Key

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

Definition at line 164 of file JSSTaggerBase.h.

◆ m_readTau1WTAKey

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

ReadDecorHandle keys for JSS moments.

Definition at line 154 of file JSSTaggerBase.h.

◆ m_readTau2WTAKey

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

Definition at line 155 of file JSSTaggerBase.h.

◆ m_readTau3WTAKey

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

Definition at line 156 of file JSSTaggerBase.h.

◆ m_readTau4WTAKey

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

Definition at line 157 of file JSSTaggerBase.h.

◆ m_readThrustMajKey

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

Definition at line 168 of file JSSTaggerBase.h.

◆ m_readTruthLabelKey

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

Truth label ReadDecorHandle key.

Definition at line 216 of file JSSTaggerBase.h.

◆ m_recognisedCuts

std::vector<std::string> SmoothedTopTagger::m_recognisedCuts = {"Mass", "mass", "Sphericity", "sphericity"}
private

Definition at line 51 of file SmoothedTopTagger.h.

◆ m_strMassCutHigh

std::string JSSTaggerBase::m_strMassCutHigh
protectedinherited

Definition at line 181 of file JSSTaggerBase.h.

◆ m_strMassCutLow

std::string JSSTaggerBase::m_strMassCutLow
protectedinherited

Strings for cut functions.

Definition at line 180 of file JSSTaggerBase.h.

◆ m_strScoreCut

std::string JSSTaggerBase::m_strScoreCut
protectedinherited

Definition at line 182 of file JSSTaggerBase.h.

◆ m_suppressOutputDependence

bool JSSTaggerBase::m_suppressOutputDependence {}
protectedinherited

Definition at line 177 of file JSSTaggerBase.h.

◆ m_tagClass

TAGCLASS JSSTaggerBase::m_tagClass {Unknown}
protectedinherited

Definition at line 101 of file JSSTaggerBase.h.

◆ m_tagType

std::string JSSTaggerBase::m_tagType
protectedinherited

Definition at line 130 of file JSSTaggerBase.h.

◆ m_tmvaCalibArea

std::string JSSTaggerBase::m_tmvaCalibArea
protectedinherited

Definition at line 126 of file JSSTaggerBase.h.

◆ m_tmvaConfigFileName

std::string JSSTaggerBase::m_tmvaConfigFileName
protectedinherited

TMVA configurations for BDT taggers.

Definition at line 124 of file JSSTaggerBase.h.

◆ m_tmvaConfigFilePath

std::string JSSTaggerBase::m_tmvaConfigFilePath
protectedinherited

Definition at line 125 of file JSSTaggerBase.h.

◆ m_truthBosonContainerName

std::string JSSTaggerBase::m_truthBosonContainerName
protectedinherited

Definition at line 211 of file JSSTaggerBase.h.

◆ m_truthLabelName

std::string JSSTaggerBase::m_truthLabelName
protectedinherited

Definition at line 213 of file JSSTaggerBase.h.

◆ m_truthLabelUseTRUTH3

bool JSSTaggerBase::m_truthLabelUseTRUTH3 {}
protectedinherited

Truth label options.

Definition at line 209 of file JSSTaggerBase.h.

◆ m_truthParticleContainerName

std::string JSSTaggerBase::m_truthParticleContainerName
protectedinherited

Definition at line 210 of file JSSTaggerBase.h.

◆ m_truthTopQuarkContainerName

std::string JSSTaggerBase::m_truthTopQuarkContainerName
protectedinherited

Definition at line 212 of file JSSTaggerBase.h.

◆ m_useECFG

bool JSSTaggerBase::m_useECFG = false
protectedinherited

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

Definition at line 139 of file JSSTaggerBase.h.

◆ m_useLSeries

bool JSSTaggerBase::m_useLSeries = false
protectedinherited

Definition at line 140 of file JSSTaggerBase.h.

◆ m_useMassCut

bool JSSTaggerBase::m_useMassCut = false
protectedinherited

Flag to indicate if mass window cut is used.

Definition at line 94 of file JSSTaggerBase.h.

◆ m_useScoreCut

bool JSSTaggerBase::m_useScoreCut = false
protectedinherited

Flag to indicate if a discriminant score is used.

Definition at line 97 of file JSSTaggerBase.h.

◆ m_varCutExprs

std::vector<std::string> SmoothedTopTagger::m_varCutExprs
private

Definition at line 31 of file SmoothedTopTagger.h.

◆ m_varCutFuncs

std::vector<std::unique_ptr<TF1> > SmoothedTopTagger::m_varCutFuncs
private

Definition at line 37 of file SmoothedTopTagger.h.

◆ m_varCutNames

std::vector<std::string> SmoothedTopTagger::m_varCutNames
private

Definition at line 34 of file SmoothedTopTagger.h.

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.

◆ m_weightConfig

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

Histograms for scale factors.

Definition at line 204 of file JSSTaggerBase.h.

◆ m_weightConfigPath

std::string JSSTaggerBase::m_weightConfigPath
protectedinherited

Path to the SF configuration root file.

Definition at line 109 of file JSSTaggerBase.h.

◆ m_weightDecorationName

std::string JSSTaggerBase::m_weightDecorationName
protectedinherited

String for scale factor decoration names.

Definition at line 197 of file JSSTaggerBase.h.

◆ m_weightFileName

std::string JSSTaggerBase::m_weightFileName
protectedinherited

Definition at line 198 of file JSSTaggerBase.h.

◆ m_weightFlavors

std::string JSSTaggerBase::m_weightFlavors
protectedinherited

Definition at line 201 of file JSSTaggerBase.h.

◆ m_weightHistogramName

std::string JSSTaggerBase::m_weightHistogramName
protectedinherited

Definition at line 199 of file JSSTaggerBase.h.

◆ m_weightHistograms

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

Definition at line 205 of file JSSTaggerBase.h.

◆ m_wkpt

std::string JSSTaggerBase::m_wkpt
protectedinherited

Tagger information.

Definition at line 129 of file JSSTaggerBase.h.

The documentation for this class was generated from the following files:
JSSTaggerBase::m_readParentKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readParentKey
Definition: JSSTaggerBase.h:175
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
LargeRJetTruthLabel::TypeEnum
TypeEnum
Definition: LargeRJetLabelEnum.h:14
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
python.SystemOfUnits.second
int second
Definition: SystemOfUnits.py:120
JSSTaggerBase::m_readTau4WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau4WTAKey
Definition: JSSTaggerBase.h:157
beamspotnt.var
var
Definition: bin/beamspotnt.py:1394
Trk::L2
@ L2
Definition: AlignModuleList.h:32
JSSTaggerBase::m_readECF1Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF1Key
Definition: JSSTaggerBase.h:159
JSSTaggerBase::getTruthLabelStr
std::string getTruthLabelStr(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Definition: JSSTaggerBase.cxx:818
get_generator_info.result
result
Definition: get_generator_info.py:21
StateLessPT_NewConfig.proxy
proxy
Definition: StateLessPT_NewConfig.py:392
JSSTaggerBase::ZBoson
@ ZBoson
Definition: JSSTaggerBase.h:100
JSSTaggerBase::resetCuts
StatusCode resetCuts(asg::AcceptData &acceptData) const
Reset cuts.
Definition: JSSTaggerBase.cxx:321
SmoothedTopTagger::m_varCutFuncs
std::vector< std::unique_ptr< TF1 > > m_varCutFuncs
Definition: SmoothedTopTagger.h:37
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
find
std::string find(const std::string &s)
return a remapped string
Definition: hcg.cxx:135
CaloCellPos2Ntuple.int
int
Definition: CaloCellPos2Ntuple.py:24
JSSTaggerBase::m_readECFG311Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG311Key
Definition: JSSTaggerBase.h:172
JSSTaggerBase::m_jetEtaMax
float m_jetEtaMax
Definition: JSSTaggerBase.h:135
SmoothedTopTagger::m_decPassSphericityKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassSphericityKey
Definition: SmoothedTopTagger.h:44
JSSTaggerBase::m_decValidEtaRangeKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidEtaRangeKey
Definition: JSSTaggerBase.h:71
JSSTaggerBase::m_efficiencyHistograms
std::map< std::string, std::unique_ptr< TH2D > > m_efficiencyHistograms
Definition: JSSTaggerBase.h:206
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
JSSTaggerBase::m_configFile
std::string m_configFile
Configuration file name.
Definition: JSSTaggerBase.h:112
JSSTaggerBase::m_configReader
TEnv m_configReader
TEnv instance to read config files.
Definition: JSSTaggerBase.h:62
RootHelpers::FindBin
Int_t FindBin(const TAxis *axis, const double x)
Definition: RootHelpers.cxx:14
JSSTaggerBase::TopQuark
@ TopQuark
Definition: JSSTaggerBase.h:100
JSSTaggerBase::tag
virtual StatusCode tag(const xAOD::Jet &jet) const =0
Decorate single jet with tagging info.
SmoothedTopTagger::m_dec_sphericitycut
SG::WriteDecorHandleKey< xAOD::JetContainer > m_dec_sphericitycut
Definition: SmoothedTopTagger.h:41
TRT_PAI_gasdata::SF
const float SF[NF]
Cross sections for Fluor.
Definition: TRT_PAI_gasdata.h:285
LargeRJetTruthLabel::other_From_t
@ other_From_t
Definition: LargeRJetLabelEnum.h:22
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
JSSTaggerBase::m_readTau2WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau2WTAKey
Definition: JSSTaggerBase.h:155
JSSTaggerBase::m_jetPtMin
float m_jetPtMin
Kinematic bounds for the jet - the units are controlled by m_ptGeV.
Definition: JSSTaggerBase.h:133
JSSTaggerBase::m_decL3Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL3Key
Definition: JSSTaggerBase.h:151
SG::ConstAccessor
Helper class to provide constant type-safe access to aux data.
Definition: ConstAccessor.h:55
JSSTaggerBase::m_decC2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decC2Key
Definition: JSSTaggerBase.h:147
JSSTaggerBase::m_decorationName
std::string m_decorationName
Decoration name.
Definition: JSSTaggerBase.h:190
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
JSSTaggerBase::m_readECFG331Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG331Key
Definition: JSSTaggerBase.h:171
LargeRJetTruthLabel::Zqq
@ Zqq
Definition: LargeRJetLabelEnum.h:20
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
dbg::ptr
void * ptr(T *p)
Definition: SGImplSvc.cxx:74
yodamerge_tmp.scale
scale
Definition: yodamerge_tmp.py:138
LargeRJetTruthLabel::qcd
@ qcd
Definition: LargeRJetLabelEnum.h:25
JSSTaggerBase::m_decValidPtRangeHighKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeHighKey
Definition: JSSTaggerBase.h:69
JSSTaggerBase::m_decValidJetContentKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidJetContentKey
Definition: JSSTaggerBase.h:73
asg::AcceptInfo::getCutName
const std::string & getCutName(unsigned int cutPosition) const
Get the name of a cut, based on the cut position (slow, avoid usage)
Definition: AcceptInfo.cxx:36
JSSTaggerBase::m_jetPtMax
float m_jetPtMax
Definition: JSSTaggerBase.h:134
JSSTaggerBase::m_calibArea
std::string m_calibArea
Location where config files live on cvmfs.
Definition: JSSTaggerBase.h:115
JSSTaggerBase::m_decWeightKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decWeightKey
WriteDecorHandle keys for SF.
Definition: JSSTaggerBase.h:225
JSSTaggerBase::m_readTruthLabelKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTruthLabelKey
Truth label ReadDecorHandle key.
Definition: JSSTaggerBase.h:216
LargeRJetTruthLabel::notruth
@ notruth
Definition: LargeRJetLabelEnum.h:24
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
JSSTaggerBase::m_decSigeffSFKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decSigeffSFKey
Definition: JSSTaggerBase.h:228
dqt_zlumi_pandas.weight
int weight
Definition: dqt_zlumi_pandas.py:189
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
JSSTaggerBase::m_decEfficiencyKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decEfficiencyKey
Definition: JSSTaggerBase.h:226
AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:85
JSSTaggerBase::m_decE3Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decE3Key
Definition: JSSTaggerBase.h:149
JSSTaggerBase::m_ptGeV
bool m_ptGeV
Flag to indicate units pT is defined in Set to false by default.
Definition: JSSTaggerBase.h:91
efficiency
void efficiency(std::vector< double > &bins, std::vector< double > &values, const std::vector< std::string > &files, const std::string &histname, const std::string &tplotname, const std::string &label="")
Definition: dependence.cxx:128
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
JSSTaggerBase::m_decTau32WTAKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau32WTAKey
Definition: JSSTaggerBase.h:144
jet
Definition: JetCalibTools_PlotJESFactors.cxx:23
AthCommonDataStore
Definition: AthCommonDataStore.h:52
LargeRJetTruthLabel::Wqq
@ Wqq
Definition: LargeRJetLabelEnum.h:17
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
CheckAppliedSFs.e3
e3
Definition: CheckAppliedSFs.py:264
SG::ReadDecorHandle
Handle class for reading a decoration on an object.
Definition: StoreGate/StoreGate/ReadDecorHandle.h:94
JSSTaggerBase::getSF
std::pair< double, double > getSF(const xAOD::Jet &jet, const std::string &truthLabelStr) const
Get scale factor and efficiency.
Definition: JSSTaggerBase.cxx:760
lumiFormat.i
int i
Definition: lumiFormat.py:85
Trk::L3
@ L3
Definition: AlignModuleList.h:32
JSSTaggerBase::m_readECFG212Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG212Key
Definition: JSSTaggerBase.h:173
PixelAthClusterMonAlgCfg.varName
string varName
end cluster ToT and charge
Definition: PixelAthClusterMonAlgCfg.py:125
asg::AcceptInfo::getNCuts
unsigned int getNCuts() const
Get the number of cuts defined.
Definition: AcceptInfo.h:46
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
JSSTaggerBase::checkKinRange
StatusCode checkKinRange(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Check and record if jet passes kinematic constraints.
Definition: JSSTaggerBase.cxx:353
SG::WriteDecorHandle
Handle class for adding a decoration to an object.
Definition: StoreGate/StoreGate/WriteDecorHandle.h:100
xAOD::VxType::PriVtx
@ PriVtx
Primary vertex.
Definition: TrackingPrimitives.h:571
JSSTaggerBase::Unknown
@ Unknown
Definition: JSSTaggerBase.h:100
checkCorrelInHIST.prefix
dictionary prefix
Definition: checkCorrelInHIST.py:391
JSSTaggerBase::m_efficiencyHistogramName
std::string m_efficiencyHistogramName
Definition: JSSTaggerBase.h:200
JSSTaggerBase::m_decTau42WTAKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau42WTAKey
Definition: JSSTaggerBase.h:145
JSSTaggerBase::calculateJSSRatios
int calculateJSSRatios(const xAOD::Jet &jet) const
Calculate JSS moment ratios in case they are not already saved TODO: Remove this once JSSMomentTools ...
Definition: JSSTaggerBase.cxx:399
python.xAODType.dummy
dummy
Definition: xAODType.py:4
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
MSG::name
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition: MsgLevel.cxx:19
JSSTaggerBase::m_decEffSFKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decEffSFKey
Definition: JSSTaggerBase.h:227
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
JSSTaggerBase::m_wkpt
std::string m_wkpt
Tagger information.
Definition: JSSTaggerBase.h:129
JSSTaggerBase::m_calcSF
bool m_calcSF
Flag to calculate scale factor.
Definition: JSSTaggerBase.h:193
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
JSSTaggerBase::m_decD2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decD2Key
Definition: JSSTaggerBase.h:148
LargeRJetTruthLabel::Wqq_From_t
@ Wqq_From_t
Definition: LargeRJetLabelEnum.h:21
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
JSSTaggerBase::initialize
virtual StatusCode initialize() override
Initialize the tool.
Definition: JSSTaggerBase.cxx:73
JSSTaggerBase::m_acceptInfo
asg::AcceptInfo m_acceptInfo
Object that stores the results for a jet.
Definition: JSSTaggerBase.h:65
ActsTrk::to_string
std::string to_string(const DetectorType &type)
Definition: GeometryDefs.h:34
JSSTaggerBase::m_decTau21WTAKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau21WTAKey
WriteDecorHandle keys for JSS moments.
Definition: JSSTaggerBase.h:143
SmoothedTopTagger::m_dec_mcut
SG::WriteDecorHandleKey< xAOD::JetContainer > m_dec_mcut
Definition: SmoothedTopTagger.h:40
PathResolverFindCalibFile
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
Definition: PathResolver.cxx:431
JSSTaggerBase::printCuts
void printCuts() const
Print configured cuts.
Definition: JSSTaggerBase.cxx:907
JSSTaggerBase::m_isMC
bool m_isMC
Definition: JSSTaggerBase.h:194
JSSTaggerBase::WBoson
@ WBoson
Definition: JSSTaggerBase.h:100
xAOD::Jet_v1
Class describing a jet.
Definition: Jet_v1.h:57
JSSTaggerBase::m_decValidKinRangeKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidKinRangeKey
Definition: JSSTaggerBase.h:72
JSSTaggerBase::m_readTau1WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau1WTAKey
ReadDecorHandle keys for JSS moments.
Definition: JSSTaggerBase.h:154
JSSTaggerBase::m_weightHistograms
std::map< std::string, std::unique_ptr< TH2D > > m_weightHistograms
Definition: JSSTaggerBase.h:205
a
TList * a
Definition: liststreamerinfos.cxx:10
asg::AcceptData::setCutResult
void setCutResult(const std::string &cutName, bool cutResult)
Set the result of a cut, based on the cut name (safer)
Definition: AcceptData.h:134
h
IJetDecorator::decorate
virtual StatusCode decorate(const xAOD::JetContainer &jets) const =0
Decorate a jet collection without otherwise modifying it.
JSSTaggerBase::passKinRange
bool passKinRange(const xAOD::Jet &jet) const
Check if jet passes kinematic constraints.
Definition: JSSTaggerBase.cxx:339
asg::AcceptData::clear
void clear()
Clear all bits.
Definition: AcceptData.h:54
asg::AcceptData::getCutResult
bool getCutResult(const std::string &cutName) const
Get the result of a cut, based on the cut name (safer)
Definition: AcceptData.h:98
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
JSSTaggerBase::m_decValidPtRangeLowKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeLowKey
Definition: JSSTaggerBase.h:70
SmoothedTopTagger::m_numTaggerVars
int m_numTaggerVars
Definition: SmoothedTopTagger.h:28
JSSTaggerBase::m_readSphericityKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSphericityKey
Definition: JSSTaggerBase.h:169
JSSTaggerBase::m_readTau3WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau3WTAKey
Definition: JSSTaggerBase.h:156
JSSTaggerBase::m_readECF2Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF2Key
Definition: JSSTaggerBase.h:160
JSSTaggerBase::m_readECF3Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF3Key
Definition: JSSTaggerBase.h:161
AthCommonMsg< AlgTool >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
JSSTaggerBase::getConfigReader
StatusCode getConfigReader()
Get configReader StatusCode.
Definition: JSSTaggerBase.cxx:283
LargeRJetTruthLabel::intToEnum
TypeEnum intToEnum(const int type)
Definition: LargeRJetLabelEnum.h:57
JSSTaggerBase::m_decPassMassKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassMassKey
Definition: JSSTaggerBase.h:76
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
SmoothedTopTagger::m_recognisedCuts
std::vector< std::string > m_recognisedCuts
Definition: SmoothedTopTagger.h:51
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
JSSTaggerBase::JSSTaggerBase
JSSTaggerBase(const std::string &name)
Default constructor - to be used in all derived classes.
Definition: JSSTaggerBase.cxx:11
defineDB.jets
list jets
Definition: JetTagCalibration/share/defineDB.py:24
dqt_zlumi_alleff_HIST.eff
int eff
Definition: dqt_zlumi_alleff_HIST.py:113
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
LargeRJetTruthLabel::tqqb
@ tqqb
Definition: LargeRJetLabelEnum.h:16
JSSTaggerBase::m_decL2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL2Key
Definition: JSSTaggerBase.h:150
JSSTaggerBase::m_weightFlavors
std::string m_weightFlavors
Definition: JSSTaggerBase.h:201
SmoothedTopTagger::m_varCutExprs
std::vector< std::string > m_varCutExprs
Definition: SmoothedTopTagger.h:31
asg::AcceptData
Definition: AcceptData.h:30
SG::DataProxy
Definition: DataProxy.h:45
pow
constexpr int pow(int base, int exp) noexcept
Definition: ap_fixedTest.cxx:15
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
SmoothedTopTagger::m_varCutNames
std::vector< std::string > m_varCutNames
Definition: SmoothedTopTagger.h:34
fitman.k
k
Definition: fitman.py:528
asg::AcceptInfo::getCutDescription
const std::string & getCutDescription(const std::string &cutName) const
Get the description of a cut, based on the cut name.
Definition: AcceptInfo.cxx:119
JSSTaggerBase::m_decTaggedKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTaggedKey
WriteDecorHandle keys for tagging bools.
Definition: JSSTaggerBase.h:68
asg::AcceptInfo::addCut
int addCut(const std::string &cutName, const std::string &cutDescription)
Add a cut; returning the cut position.
Definition: AcceptInfo.h:53
JSSTaggerBase::m_containerName
std::string m_containerName
Configurable members.
Definition: JSSTaggerBase.h:106
Generated on Sun Dec 22 2024 21:44:31 for ATLAS Offline Software by   doxygen 1.8.18