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

#include <JSSWTopTaggerANN.h>

Inheritance diagram for JSSWTopTaggerANN:
Collaboration diagram for JSSWTopTaggerANN:

Public Member Functions

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

Protected Types

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

Protected Member Functions

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"}
 
const int m_nWarnMax = 10
 Maximum number of warnings. More...
 
std::atomic< int > m_nWarnVar
 Warning counters. More...
 
bool m_ptGeV = false
 Flag to indicate units pT is defined in Set to false by default. More...
 
bool m_useMassCut = false
 Flag to indicate if mass window cut is used. More...
 
bool m_useScoreCut = false
 Flag to indicate if a discriminant score is used. More...
 
TAGCLASS m_tagClass {Unknown}
 
std::string m_containerName
 Configurable members. More...
 
bool m_isSmallRJet = false
 
std::string m_weightConfigPath
 Path to the SF configuration root file. More...
 
std::string m_configFile
 Configuration file name. More...
 
std::string m_calibArea
 Location where config files live on cvmfs. More...
 
std::string m_kerasConfigFileName
 Keras configurations for ML taggers. More...
 
std::string m_kerasConfigFilePath
 
std::string m_kerasConfigOutputName
 
std::string m_kerasCalibArea
 
std::string m_tmvaConfigFileName
 TMVA configurations for BDT taggers. More...
 
std::string m_tmvaConfigFilePath
 
std::string m_tmvaCalibArea
 
std::string m_wkpt
 Tagger information. More...
 
std::string m_tagType
 
float m_jetPtMin {}
 Kinematic bounds for the jet - the units are controlled by m_ptGeV. More...
 
float m_jetPtMax {}
 
float m_jetEtaMax {}
 
bool m_useECFG = false
 Flags controlling whether generalized ECF moments or L-series ratios are needed TODO: Implement the functionality controlled by these. More...
 
bool m_useLSeries = false
 
SG::WriteDecorHandleKey< xAOD::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"}
 
Gaudi::Property< bool > m_suppressInputDependence {this, "SuppressInputDependence", false}
 
Gaudi::Property< bool > m_suppressOutputDependence {this, "SuppressOutputDependence", false}
 
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

double getScore (const xAOD::Jet &jet) const
 Retrieve score for a given ANN type (top/W) More...
 
std::map< std::string, std::map< std::string, double > > getJetProperties (const xAOD::Jet &jet) const
 Update the jet substructure variables for each jet to use in ANN. More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

std::unique_ptr< lwt::LightweightGraph > m_lwnn
 ANN tools. More...
 
std::map< std::string, std::map< std::string, double > > m_ANN_inputValues
 
std::vector< std::string > m_out_names
 Internal stuff to keep track of the output node for the NN. More...
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

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

97 {Unknown, WBoson, ZBoson, TopQuark};

Constructor & Destructor Documentation

◆ JSSWTopTaggerANN()

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

Constructor.

Definition at line 10 of file JSSWTopTaggerANN.cxx.

10  :
11  JSSTaggerBase( name ),
12  m_lwnn(nullptr)
13 {
14 
15 }

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

WTA N-subjettiness ratios

ECF ratios

Definition at line 412 of file JSSTaggerBase.cxx.

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

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

366  {
367 
368  float scale = 1.0;
369  if ( m_ptGeV ) scale = 1.e3;
370 
371  bool passKinRange = true;
372 
375  if ( std::abs(jet.eta()) > m_jetEtaMax ) {
376  ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (|eta| < " << m_jetEtaMax << "). Jet eta = " << jet.eta() );
377  acceptData.setCutResult( "ValidEtaRange", false );
378  passKinRange = false;
379  }
380 
381  if ( jet.pt() < m_jetPtMin * scale ) {
382  ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (pT > " << m_jetPtMin * scale / 1.e3 << "). Jet pT = " << jet.pt() / 1.e3 << " GeV" );
383  acceptData.setCutResult( "ValidPtRangeLow", false );
384  passKinRange = false;
385  }
386 
387  if ( jet.pt() > m_jetPtMax * scale ) {
388  ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (pT < " << m_jetPtMax * scale / 1.e3 << "). Jet pT = " << jet.pt() / 1.e3 << " GeV" );
389  acceptData.setCutResult( "ValidPtRangeHigh", false );
390  passKinRange = false;
391  }
392 
394  SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidPtRangeHigh(m_decValidPtRangeHighKey);
395  SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidPtRangeLow(m_decValidPtRangeLowKey);
396  SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidEtaRange(m_decValidEtaRangeKey);
397  SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidKinRange(m_decValidKinRangeKey);
398 
400  decValidPtRangeHigh(jet) = acceptData.getCutResult( "ValidPtRangeHigh" );
401  decValidPtRangeLow(jet) = acceptData.getCutResult( "ValidPtRangeLow" );
402  decValidEtaRange(jet) = acceptData.getCutResult( "ValidEtaRange" );
403  decValidKinRange(jet) = passKinRange;
404 
405  return StatusCode::SUCCESS;
406 
407 }

◆ declareGaudiProperty() [1/4]

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

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

Definition at line 170 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [2/4]

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

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [3/4]

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

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

Definition at line 184 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [4/4]

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

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

Definition at line 199 of file AthCommonDataStore.h.

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

◆ declareProperty() [1/6]

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

Declare a new Gaudi property.

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

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

Definition at line 245 of file AthCommonDataStore.h.

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

◆ declareProperty() [2/6]

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

Declare a new Gaudi property.

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

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

Definition at line 221 of file AthCommonDataStore.h.

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

◆ declareProperty() [3/6]

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

Definition at line 259 of file AthCommonDataStore.h.

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

◆ declareProperty() [4/6]

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

Declare a new Gaudi property.

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

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

Definition at line 333 of file AthCommonDataStore.h.

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

◆ declareProperty() [5/6]

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

Declare a new Gaudi property.

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

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

Definition at line 352 of file AthCommonDataStore.h.

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

◆ declareProperty() [6/6]

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

Definition at line 145 of file AthCommonDataStore.h.

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

◆ decorate()

StatusCode JSSWTopTaggerANN::decorate ( const xAOD::JetContainer &  jets) const
overridevirtual

Decorate jet collection with tagging info.

Loop over jet collection and decorate each jet.

Reset the AcceptData cut results

Check each kinematic constraint Print warnings using counters

Get the relevant attributes of the jet Mass and pt - note that this will depend on the configuration of the calibration used

Get score

Evaluate the values of the upper and lower mass bounds and the d2 cut

Print cut criteria and jet values

Get SF weight

Inefficiency SF is directly used

Efficiency SF

Calculate inefficiency SF

If inefficiency SF is not available, SF is always 1.0

Decorate cut information if needed

Decorate values

Cut summary

Set the AcceptData depending on whether it is a W/Z or a top tagger

Reimplemented from JSSTaggerBase.

Definition at line 216 of file JSSWTopTaggerANN.cxx.

216  {
217 
218  decorateJSSRatios(jets);
219 
220  ATH_MSG_DEBUG("Using AthAnalysis code. ANN code.");
221  SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidPtRangeHigh(m_decValidPtRangeHighKey);
222  SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidPtRangeLow(m_decValidPtRangeLowKey);
223  SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidEtaRange(m_decValidEtaRangeKey);
224  SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidKinRange(m_decValidKinRangeKey);
225 
226  SG::WriteDecorHandle<xAOD::JetContainer, bool> decPassMass(m_decPassMassKey);
227  SG::WriteDecorHandle<xAOD::JetContainer, bool> decPassScore(m_decPassScoreKey);
228  SG::WriteDecorHandle<xAOD::JetContainer, bool> decTagged(m_decTaggedKey);
229  SG::WriteDecorHandle<xAOD::JetContainer, float> decCutMLow(m_decCutMLowKey);
230  SG::WriteDecorHandle<xAOD::JetContainer, float> decCutMHigh(m_decCutMHighKey);
231  SG::WriteDecorHandle<xAOD::JetContainer, float> decScoreCut(m_decScoreCutKey);
232  SG::WriteDecorHandle<xAOD::JetContainer, float> decScoreValue(m_decScoreValueKey);
233 
234  SG::WriteDecorHandle<xAOD::JetContainer, float> decWeight(m_decWeightKey);
235  SG::WriteDecorHandle<xAOD::JetContainer, float> decEfficiency(m_decEfficiencyKey);
236  SG::WriteDecorHandle<xAOD::JetContainer, float> decEffSF(m_decEffSFKey);
237  SG::WriteDecorHandle<xAOD::JetContainer, float> decSigeffSF(m_decSigeffSFKey);
238 
239 
240  for(const xAOD::Jet* jet : jets){
241 
242  asg::AcceptData acceptData( &m_acceptInfo );
244  ATH_CHECK( resetCuts( acceptData ) );
245 
246  float scale = m_ptGeV ? 1.e3 : 1.0;
247 
248  bool passKinRange = true;
249 
252  if ( std::abs(jet->eta()) > m_jetEtaMax ) {
253  ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (|eta| < " << m_jetEtaMax << "). Jet eta = " << jet->eta() );
254  acceptData.setCutResult( "ValidEtaRange", false );
255  passKinRange = false;
256  }
257 
258  if ( jet->pt() < m_jetPtMin * scale ) {
259  ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (pT > " << m_jetPtMin * scale / 1.e3 << "). Jet pT = " << jet->pt() / 1.e3 << " GeV" );
260  acceptData.setCutResult( "ValidPtRangeLow", false );
261  passKinRange = false;
262  }
263 
264  if ( jet->pt() > m_jetPtMax * scale ) {
265  ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (pT < " << m_jetPtMax * scale / 1.e3 << "). Jet pT = " << jet->pt() / 1.e3 << " GeV" );
266  acceptData.setCutResult( "ValidPtRangeHigh", false );
267  passKinRange = false;
268  }
269  decValidPtRangeHigh(*jet) = acceptData.getCutResult( "ValidPtRangeHigh" );
270  decValidPtRangeLow(*jet) = acceptData.getCutResult( "ValidPtRangeLow" );
271  decValidEtaRange(*jet) = acceptData.getCutResult( "ValidEtaRange" );
272  decValidKinRange(*jet) = passKinRange;
273 
274 
275 
278  float jet_pt = jet->pt()/1000.;
279  float jet_mass = jet->m()/1000.;
280 
282  float jet_score = getScore(*jet);
283 
284 
286  float cut_mass_low = m_funcMassCutLow ->Eval(jet_pt);
287  float cut_mass_high = m_funcMassCutHigh->Eval(jet_pt);
288  float cut_score = m_funcScoreCut ->Eval(jet_pt);
289 
291  ATH_MSG_VERBOSE( "Cut values : Mass window = [" << cut_mass_low << "," << cut_mass_high << "], score cut = " << cut_score );
292  ATH_MSG_VERBOSE( "Jet values : Mass = " << jet_mass << ", score = " << jet_score );
293 
295  float weight = 1.0;
296  float effSF = 1.0;
297  float sigeffSF = 1.0;
298  float efficiency = 1.0;
299 
300  if ( m_isMC && m_calcSF ) {
301 
302  std::string truthLabelStr = getTruthLabelStr( *jet, acceptData );
303  std::tie(effSF, efficiency) = getSF( *jet, truthLabelStr );
304 
305  // calculate signal efficiency SF
306  if ( m_weightHistograms.count("t_qqb") ) {
307  sigeffSF = getSF(*jet, "t_qqb").first;
308  } else if ( m_weightHistograms.count("V_qq") ) {
309  sigeffSF = getSF(*jet, "V_qq").first;
310  } else if ( m_weightHistograms.count("t") ){
311  sigeffSF = getSF(*jet, "t").first;
312  } else {
313  sigeffSF = 1.0;
314  }
315 
317  if ( m_weightFlavors.find("fail") != std::string::npos ) {
318  weight = effSF;
319  }
320 
321  else {
322 
324  if ( jet_score > cut_score ) {
325  weight = effSF;
326  }
327 
329  else {
331  if ( !m_efficiencyHistogramName.empty() && efficiency < 1.0 ) {
332  weight = ( 1. - effSF * efficiency ) / ( 1. - efficiency );
333  }
334 
335  }
336  }
337 
338  }
339 
340 
341  decWeight(*jet) = weight;
342  decEfficiency(*jet) = efficiency;
343  decEffSF(*jet) = effSF;
344  decSigeffSF(*jet) = sigeffSF;
345 
347  ATH_MSG_DEBUG( "Decorating with score" );
348 
350  decCutMLow(*jet) = cut_mass_low;
351  decCutMHigh(*jet) = cut_mass_high;
352  decScoreCut(*jet) = cut_score;
353  decScoreValue(*jet) = jet_score;
354 
356  bool passCuts = true;
357 
359  if ( m_tagClass == TAGCLASS::WBoson || m_tagClass == TAGCLASS::ZBoson ) {
360  ATH_MSG_VERBOSE( "Determining WZ tag return" );
361  if ( jet_mass > cut_mass_low ) acceptData.setCutResult( "PassMassLow", true );
362  if ( jet_mass < cut_mass_high ) acceptData.setCutResult( "PassMassHigh", true );
363  if ( jet_score > cut_score ) acceptData.setCutResult( "PassScore", true );
364  decPassMass(*jet) = acceptData.getCutResult( "PassMassLow" ) && acceptData.getCutResult( "PassMassHigh" );
365  passCuts = passCuts && acceptData.getCutResult( "PassMassLow" ) && acceptData.getCutResult( "PassMassHigh" );
366  }
367  else if ( m_tagClass == TAGCLASS::TopQuark ) {
368  ATH_MSG_VERBOSE( "Determining TopQuark tag return" );
369  if( jet_mass > cut_mass_low ) acceptData.setCutResult( "PassMassLow", true );
370  if( jet_score > cut_score ) acceptData.setCutResult( "PassScore", true );
371  decPassMass(*jet) = acceptData.getCutResult( "PassMassLow" );
372  passCuts = passCuts && acceptData.getCutResult( "PassMassLow" );
373  }
374 
375  decPassScore(*jet) = acceptData.getCutResult( "PassScore" );
376 
377  passCuts = passCuts && acceptData.getCutResult( "PassScore" );
378 
379  decTagged(*jet) = passCuts;
380 
381 
382  }
383 
384  return StatusCode::SUCCESS;
385 
386 }

◆ decorateJSSRatios()

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

Create write decor handles

WTA N-subjettiness ratios

ECF ratios

Definition at line 520 of file JSSTaggerBase.cxx.

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

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

643  {
644 
645  int indexPV = -1;
646 
647  const xAOD::VertexContainer* vxCont = nullptr;
648  if ( evtStore()->retrieve( vxCont, "PrimaryVertices" ).isFailure() ) {
649  ATH_MSG_WARNING( "Unable to retrieve primary vertex container PrimaryVertices" );
650  }
651  else {
652  int vtx_counter = 0;
653  for ( const auto *vx : *vxCont ) {
654  if ( vx->vertexType()==xAOD::VxType::PriVtx ) {
655  indexPV = vtx_counter;
656  break;
657  }
658  vtx_counter++;
659  }
660  }
661 
662  return indexPV;
663 
664 }

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

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

◆ getJetProperties()

std::map< std::string, std::map< std::string, double > > JSSWTopTaggerANN::getJetProperties ( const xAOD::Jet &  jet) const
private

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

Map to store inputs

Create common read decor handles

Mass and pT It is assumed that these are the combined and calibrated mass and pT

Splitting scales

Energy Correlation Functions

Tau21 WTA

Other moments

Create top quark read decor handles

Additional splitting Scales

e3 := normalized ECF3/ECF1**3

N-subjettiness

Qw observable for top tagging

Definition at line 508 of file JSSWTopTaggerANN.cxx.

508  {
509 
511  std::map< std::string, std::map<std::string, double> > ANN_inputs;
512  std::map< std::string, double > ANN_inputValues;
513 
514  ATH_MSG_DEBUG( "Loading variables for common ANN tagger" );
515 
517  SG::ReadDecorHandle<xAOD::JetContainer, float> readSplit12(m_readSplit12Key);
518 
521  ANN_inputValues["CaloTACombinedMassUncorrelated"] = jet.m();
522  ANN_inputValues["JetpTCorrByCombinedMass"] = jet.pt();
523 
525  ANN_inputValues["Split12"] = readSplit12(jet);
526 
528  static const SG::AuxElement::ConstAccessor<float> C2("C2");
529  static const SG::AuxElement::ConstAccessor<float> D2("D2");
530  ANN_inputValues["C2"] = C2(jet);
531  ANN_inputValues["D2"] = D2(jet);
532 
534  static const SG::AuxElement::ConstAccessor<float> Tau21WTA("Tau21_wta");
535  ANN_inputValues["Tau21_wta"] = Tau21WTA(jet);
536 
537  if ( m_tagClass == TAGCLASS::WBoson ) {
538 
539  ATH_MSG_DEBUG( "Loading variables for W boson tagger" );
540 
542  ANN_inputValues["FoxWolfram20"] = jet.getAttribute<float>("FoxWolfram2") / jet.getAttribute<float>("FoxWolfram0");
543  ANN_inputValues["PlanarFlow"] = jet.getAttribute<float>("PlanarFlow");
544  ANN_inputValues["Angularity"] = jet.getAttribute<float>("Angularity");
545  ANN_inputValues["Aplanarity"] = jet.getAttribute<float>("Aplanarity");
546  ANN_inputValues["ZCut12"] = jet.getAttribute<float>("ZCut12");
547  ANN_inputValues["KtDR"] = jet.getAttribute<float>("KtDR");
548 
549  int pv_location = findPV();
550 
551  if(pv_location != -1){
552  int jet_ntrk = GetUnGroomTracks(jet, pv_location);
553  if(jet_ntrk>=0) ANN_inputValues["Ntrk500"] = jet_ntrk;
554  else{
555  ATH_MSG_ERROR("Either the ungroomed parent jet doesn't have 'NumTrkPt500' as an attribute or the parent link is broken");
556  ANN_inputValues["Ntrk500"] = -999;
557  }
558  }
559  else {
560  ATH_MSG_ERROR("Could not find a primary vertex");
561  ANN_inputValues["Ntrk500"] = -999;
562  }
563  }
564 
565  else if ( m_tagClass == TAGCLASS::TopQuark ) {
566 
567  ATH_MSG_DEBUG( "Loading variables for top quark tagger" );
568 
570  SG::ReadDecorHandle<xAOD::JetContainer, float> readTau1WTA(m_readTau1WTAKey);
571  SG::ReadDecorHandle<xAOD::JetContainer, float> readTau2WTA(m_readTau2WTAKey);
572  SG::ReadDecorHandle<xAOD::JetContainer, float> readTau3WTA(m_readTau3WTAKey);
573  SG::ReadDecorHandle<xAOD::JetContainer, float> readSplit23(m_readSplit23Key);
574  SG::ReadDecorHandle<xAOD::JetContainer, float> readQw(m_readQwKey);
575 
577  ANN_inputValues["Split23"] = readSplit23(jet);
578 
580  static const SG::AuxElement::ConstAccessor<float> e3("e3");
581  ANN_inputValues["e3"] = e3(jet);
582 
584  ANN_inputValues["Tau1_wta"] = readTau1WTA(jet);
585  ANN_inputValues["Tau2_wta"] = readTau2WTA(jet);
586  ANN_inputValues["Tau3_wta"] = readTau3WTA(jet);
587 
588  static const SG::AuxElement::ConstAccessor<float> Tau32WTA("Tau32_wta");
589  ANN_inputValues["Tau32_wta"] = Tau32WTA(jet);
590 
592  ANN_inputValues["Qw"] = readQw(jet);
593 
594  }
595 
596  else {
597  ATH_MSG_ERROR( "Loading variables failed because the tagger type is not supported" );
598  }
599 
600  ANN_inputs["node_0"] = ANN_inputValues;
601 
602  return ANN_inputs;
603 
604 }

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

◆ getScore()

double JSSWTopTaggerANN::getScore ( const xAOD::Jet &  jet) const
private

Retrieve score for a given ANN type (top/W)

Create input dictionary map<string,double> for argument to lwtnn

Evaluate the network response

Obtain the output associated with the single output node

Check that input variables are valid

Definition at line 474 of file JSSWTopTaggerANN.cxx.

474  {
475 
477  std::map<std::string, std::map<std::string,double>> ANN_inputs = getJetProperties(jet);
478 
480  auto discriminant = m_lwnn->compute(ANN_inputs);
481 
483  double ANNscore = -666.;
484 
486  bool validVars = true;
487  static const SG::AuxElement::ConstAccessor<float> Tau21WTA("Tau21_wta");
488  if ( Tau21WTA(jet) < 0.0 ) validVars = false;
489  if ( m_tagClass == TAGCLASS::TopQuark ) {
490  static const SG::AuxElement::ConstAccessor<float> Tau32WTA("Tau32_wta");
491  if ( Tau32WTA(jet) < 0.0 ) validVars = false;
492  }
493 
494  if ( !validVars ) {
495 
496  if ( m_nWarnVar++ < m_nWarnMax ) ATH_MSG_INFO( "One (or more) tagger input variable has an out-of-range value, setting score to -666" );
497  else ATH_MSG_DEBUG( "One (or more) tagger input variable has an out-of-range value, setting score to -666" );
498 
499  return ANNscore;
500 
501  }
502 
503  ANNscore = discriminant.at(m_out_names.at(0));
504 
505  return ANNscore;
506 }

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

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

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

837  {
839  std::string truthLabelStr;
840 
842  SG::ReadDecorHandle<xAOD::JetContainer, int> readTruthLabel(m_readTruthLabelKey);
843  LargeRJetTruthLabel::TypeEnum jetContainment = LargeRJetTruthLabel::intToEnum(readTruthLabel(jet));
844 
846  if ( m_weightHistograms.count("t_qqb") ) {
847 
849  if ( jetContainment==LargeRJetTruthLabel::tqqb ) {
850  truthLabelStr = "t_qqb";
851  }
853  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
854  truthLabelStr = "q";
855  }
856 
857  }
859  else if ( m_weightHistograms.count("V_qq_passMpassD2") ) {
860 
862  if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t ) {
863  truthLabelStr = "t_";
864  }
866  else if ( jetContainment==LargeRJetTruthLabel::Wqq || jetContainment==LargeRJetTruthLabel::Zqq || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
867  truthLabelStr = "V_qq_";
868  }
870  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
871  truthLabelStr = "q_";
872  }
873 
875  if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") && acceptData.getCutResult("PassD2") ) {
876  truthLabelStr += "passMpassD2";
877  }
879  else if ( !(acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh")) && acceptData.getCutResult("PassD2") ) {
880  truthLabelStr += "failMpassD2";
881  }
883  else if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") && !acceptData.getCutResult("PassD2") ) {
884  truthLabelStr += "passMfailD2";
885  }
887  else{
888  truthLabelStr += "failMfailD2";
889  }
890 
891  }
892 
894  else if ( m_weightHistograms.count("V_qq") ) {
895 
897  if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t ) {
898  truthLabelStr = "t";
899  }
901  else if ( jetContainment==LargeRJetTruthLabel::Wqq || jetContainment==LargeRJetTruthLabel::Zqq || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
902  truthLabelStr = "V_qq";
903  }
905  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
906  truthLabelStr = "q";
907  }
908  }
909 
910  // inclusive top tagger
911  else {
913  if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
914  truthLabelStr = "t";
915  }
917  else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
918  truthLabelStr = "q";
919  }
920 
921  }
922 
923  return truthLabelStr;
924 }

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

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

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

705  {
706  if ( !m_calcSF ) return StatusCode::SUCCESS;
707 
708  float weight = 1.0;
709  float effSF = 1.0;
710  float sigeffSF = 1.0;
711  float efficiency = 1.0;
712 
713  if ( m_isMC ) {
714 
715  std::string truthLabelStr = getTruthLabelStr( jet, acceptData );
716  std::tie(effSF, efficiency) = getSF( jet, truthLabelStr );
717 
718  // calculate signal efficiency SF
719  if ( m_weightHistograms.count("t_qqb") ) {
720  sigeffSF = getSF(jet, "t_qqb").first;
721  } else if ( m_weightHistograms.count("V_qq") ) {
722  sigeffSF = getSF(jet, "V_qq").first;
723  } else if ( m_weightHistograms.count("t") ){
724  sigeffSF = getSF(jet, "t").first;
725  } else {
726  sigeffSF = 1.0;
727  }
728 
730  if ( m_weightFlavors.find("fail") != std::string::npos ) {
731  weight = effSF;
732  }
733 
734  else {
735 
737  if ( passSel ) {
738  weight = effSF;
739  }
740 
742  else {
744  if ( m_efficiencyHistogramName.empty() ) {
745  weight = 1.0;
746  }
747  else if ( efficiency < 1.0 ) {
748  weight = ( 1. - effSF * efficiency ) / ( 1. - efficiency );
749  }
750  else {
751  weight = 1.0;
752  }
753  }
754  }
755 
756  }
757 
758  else {
759  weight = 1.0;
760  }
761 
763  SG::WriteDecorHandle<xAOD::JetContainer, float> decWeight(m_decWeightKey);
764  SG::WriteDecorHandle<xAOD::JetContainer, float> decEfficiency(m_decEfficiencyKey);
765  SG::WriteDecorHandle<xAOD::JetContainer, float> decEffSF(m_decEffSFKey);
766  SG::WriteDecorHandle<xAOD::JetContainer, float> decSigeffSF(m_decSigeffSFKey);
767 
769  decWeight(jet) = weight;
770  decEfficiency(jet) = efficiency;
771  decEffSF(jet) = effSF;
772  decSigeffSF(jet) = sigeffSF;
773 
774  return StatusCode::SUCCESS;
775 
776 }

◆ initialize()

StatusCode JSSWTopTaggerANN::initialize ( )
overridevirtual

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

Initialize the tagger.

Pt values are defined in GeV

Use mass cut

Use discriminant score cut

Get configReader

Get tagger type

Get the CVMFS calib area where stuff is stored If this is set to "Local" then it will look for the config file in the share space

Get the name/path of the JSON config

Get the name of the Keras output node

Get min and max jet mass. The unit is GeV now. Need to be consistent with ATLAS convention in the future

Get min and max jet pt. The unit is GeV now. Need to be consistent with ATLAS convention in the future

Get cut for ANN score

Get the decoration name

Get the scale factor configuration

Get truth label name information

Print out the configuration parameters for viewing

No config file

Assume the cut functions have been set through properties. Check they are non empty

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

Convert the JSON config file name to the full path

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

Read json file for ANN weights

Build the network

Set internal tagger type

Set the possible states that the tagger can be left in after the JSSTaggerBase::tag() function is called

Loop over and print out the cuts that have been configured

Call base class initialize

Reimplemented from JSSTaggerBase.

Definition at line 18 of file JSSWTopTaggerANN.cxx.

18  {
19 
20  ATH_MSG_INFO( "Initializing JSSWTopTaggerANN tool" );
21 
23  m_ptGeV = true;
24 
26  m_useMassCut = true;
27 
29  m_useScoreCut = true;
30 
31  if( ! m_configFile.empty() ) {
32 
34  ATH_CHECK( getConfigReader() );
35 
37  m_tagType = m_configReader.GetValue("TaggerType" ,"");
38 
41  m_kerasCalibArea = m_configReader.GetValue("CalibAreaKeras" ,"");
42 
44  m_kerasConfigFileName = m_configReader.GetValue("KerasConfigFile" ,"");
45 
47  m_kerasConfigOutputName = m_configReader.GetValue("KerasOutput" ,"");
48 
50  m_strMassCutLow = m_configReader.GetValue("MassCutLow_in_GeV" ,"");
51  m_strMassCutHigh = m_configReader.GetValue("MassCutHigh_in_GeV" ,"");
52 
54  m_jetPtMin = m_configReader.GetValue("pTCutLow_in_GeV", 200.0);
55  m_jetPtMax = m_configReader.GetValue("pTCutHigh_in_GeV", 2000.0);
56 
58  m_strScoreCut = m_configReader.GetValue("ScoreCut" ,"");
59 
61  m_decorationName = m_configReader.GetValue("DecorationName" ,"");
62 
64  m_calcSF = m_configReader.GetValue("CalcSF", false);
65  if ( m_calcSF ) {
66  m_weightDecorationName = m_configReader.GetValue("WeightDecorationName", "");
67  m_weightFileName = m_configReader.GetValue("WeightFile", "");
68  m_weightHistogramName = m_configReader.GetValue("WeightHistogramName", "");
69  m_efficiencyHistogramName = m_configReader.GetValue("EfficiencyHistogramName", "");
70  m_weightFlavors = m_configReader.GetValue("WeightFlavors", "");
71 
73  m_truthLabelName = m_configReader.GetValue("TruthLabelName" , "R10TruthLabel_R21Consolidated");
74  }
75 
77  ATH_MSG_INFO( "Configurations Loaded :");
78  ATH_MSG_INFO( "tagType : " << m_tagType );
79  ATH_MSG_INFO( "calibarea_keras : " << m_kerasCalibArea );
80  ATH_MSG_INFO( "kerasConfigFileName : " << m_kerasConfigFileName );
81  ATH_MSG_INFO( "kerasConfigOutputName : " << m_kerasConfigOutputName );
82  ATH_MSG_INFO( "strMassCutLow : " << m_strMassCutLow );
83  ATH_MSG_INFO( "strMassCutHigh : " << m_strMassCutHigh );
84  ATH_MSG_INFO( "pTCutLow : " << m_jetPtMin );
85  ATH_MSG_INFO( "pTCutHigh : " << m_jetPtMax );
86  ATH_MSG_INFO( "strScoreCut : " << m_strScoreCut );
87  ATH_MSG_INFO( "decorationName : " << m_decorationName );
88  if ( m_calcSF ) {
89  ATH_MSG_INFO( "weightDecorationName : " << m_weightDecorationName );
90  ATH_MSG_INFO( "weightFile : " << m_weightFileName );
91  ATH_MSG_INFO( "weightHistogramName : " << m_weightHistogramName );
92  ATH_MSG_INFO( "efficiencyHistogramName : "<<m_efficiencyHistogramName );
93  ATH_MSG_INFO( "weightFlavors : " << m_weightFlavors );
94  ATH_MSG_INFO( "TruthLabelName : " << m_truthLabelName );
95  }
96  }
97  else {
98  if ( (m_kerasConfigFileName.empty() ||
101  m_kerasConfigOutputName.empty() ||
102  m_strScoreCut.empty() ||
103  m_strMassCutLow.empty() ||
104  m_strMassCutHigh.empty() ||
105  m_decorationName.empty() ||
106  m_weightFileName.empty()) ||
107  ((m_weightDecorationName.empty() ||
108  m_weightHistogramName.empty() ||
109  m_weightFlavors.empty()) && m_calcSF) ) {
110  ATH_MSG_ERROR( "No config file provided OR you haven't manually specified all needed parameters" ) ;
111  ATH_MSG_ERROR( "Please read the TWiki for this tool" );
112  return StatusCode::FAILURE;
113  }
114 
115  }
116 
117  ATH_MSG_INFO( "Mass cut low : " << m_strMassCutLow );
118  ATH_MSG_INFO( "Mass cut High : " << m_strMassCutHigh );
119  ATH_MSG_INFO( "Score cut low : " << m_strScoreCut );
120 
122  if ( m_kerasCalibArea.empty() ) {
123  ATH_MSG_ERROR( "You need to specify where the calibarea is as either being Local or on CVMFS" );
124  return StatusCode::FAILURE;
125  }
126  else if( !m_kerasCalibArea.compare("Local") ) {
127  std::string localCalibArea = "BoostedJetTaggers/JSSWTopTaggerANN/";
128  ATH_MSG_INFO( "Using Local calibarea " << localCalibArea );
130  m_kerasConfigFilePath = PathResolverFindCalibFile(localCalibArea+m_kerasConfigFileName);
131  if ( m_calcSF )
132  m_weightConfigPath = PathResolverFindCalibFile(localCalibArea+m_weightFileName);
133  }
134  else {
135  ATH_MSG_INFO( "Using CVMFS calibarea" );
138  m_kerasConfigFilePath = PathResolverFindCalibFile( (m_kerasCalibArea+m_kerasConfigFileName).c_str() );
139  if ( m_calcSF )
140  m_weightConfigPath = PathResolverFindCalibFile( (m_kerasCalibArea+m_weightFileName).c_str());
141  }
142 
144  ATH_MSG_INFO( "ANN Tagger configured with: " << m_kerasConfigFilePath );
145 
146  std::ifstream input_cfg( m_kerasConfigFilePath.c_str() );
147 
148  if ( !input_cfg.is_open() ) {
149  ATH_MSG_ERROR( "Error openning config file: " << m_kerasConfigFilePath );
150  ATH_MSG_ERROR( "Are you sure that the file exists at this path?" );
151  return StatusCode::FAILURE;
152  }
153 
154  lwt::GraphConfig config = lwt::parse_json_graph( input_cfg );
155 
156  for ( auto& input_node: config.inputs ) {
157  ATH_MSG_INFO( " input node: " << input_node.name );
158  for ( auto& input: input_node.variables ) {
159  ATH_MSG_INFO( " " << input );
160  }
161  }
162 
163  auto output_node_name = config.outputs.begin()->first;
164  m_out_names = config.outputs.at(output_node_name).labels;
165 
166  ATH_MSG_INFO( "Keras Network NLayers: " << config.layers.size() );
167 
168  m_lwnn = std::make_unique< lwt::LightweightGraph >(config, output_node_name);
169 
171  try {
172  m_lwnn = std::make_unique<lwt::LightweightGraph>(config, output_node_name);
173  } catch (lwt::NNConfigurationException& exc) {
174  ATH_MSG_ERROR( "NN configuration problem: " << exc.what() );
175  return StatusCode::FAILURE;
176  }
177 
179  if ( !m_tagType.compare("TopQuark") ) {
180  ATH_MSG_DEBUG( "This is a top quark tagger" );
181  m_tagClass = TAGCLASS::TopQuark;
182  }
183  else if ( !m_tagType.compare("WBoson") ) {
184  ATH_MSG_DEBUG( "This is a W boson tagger" );
185  m_tagClass = TAGCLASS::WBoson;
186  }
187  else if ( !m_tagType.compare("ZBoson") ) {
188  ATH_MSG_DEBUG( "This is a Z boson tagger" );
189  m_tagClass = TAGCLASS::ZBoson;
190  }
191  else {
192  ATH_MSG_ERROR( "I can't tell what kind of tagger your configuration is for." );
193  return StatusCode::FAILURE;
194  }
195 
197  m_acceptInfo.addCut( "PassMassLow" , "mJet > mCutLow" );
198  m_acceptInfo.addCut( "PassScore" , "ScoreJet > ScoreCut" );
199  if ( m_tagClass == TAGCLASS::WBoson || m_tagClass == TAGCLASS::ZBoson ) {
200  m_acceptInfo.addCut( "PassMassHigh", "mJet < mCutHigh" );
201  }
202 
204  printCuts();
205 
207  ATH_CHECK( JSSTaggerBase::initialize() );
208 
209  ATH_MSG_INFO( "ANN Tagger tool initialized" );
210 
211  return StatusCode::SUCCESS;
212 
213 }

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

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

◆ 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, LundVariablesTool, 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 926 of file JSSTaggerBase.cxx.

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

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

334  {
335 
337  acceptData.clear();
338 
340  acceptData.setCutResult( "ValidJetContent", true );
341  acceptData.setCutResult( "ValidEventContent", true );
342 
343  acceptData.setCutResult( "ValidPtRangeHigh", true );
344  acceptData.setCutResult( "ValidPtRangeLow" , true );
345  acceptData.setCutResult( "ValidEtaRange" , true );
346 
347  return StatusCode::SUCCESS;
348 
349 }

◆ 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 JSSWTopTaggerANN::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

Get the relevant attributes of the jet Mass and pt - note that this will depend on the configuration of the calibration used

Calculate NSubjettiness and ECF ratios

Get ANN score for the jet

Evaluate the values of the upper and lower mass bounds and the d2 cut

Print cut criteria and jet values

Get SF weight

Decorate cut information if needed

Create WriteDecorHandles

Decorate values

Cut summary

Set the AcceptData depending on whether it is a W/Z or a top tagger

Implements JSSTaggerBase.

Definition at line 390 of file JSSWTopTaggerANN.cxx.

390  {
391 
392  ATH_MSG_DEBUG("Using AnalysisBase code. ANN code.");
393 
395  asg::AcceptData acceptData( &m_acceptInfo );
396 
398  ATH_CHECK( resetCuts( acceptData ) );
400  ATH_CHECK( checkKinRange( jet, acceptData ) );
401 
404  float jet_pt = jet.pt()/1000.;
405  float jet_mass = jet.m()/1000.;
406 
408  calculateJSSRatios(jet);
409 
411  float jet_score = getScore(jet);
412 
414  float cut_mass_low = m_funcMassCutLow ->Eval(jet_pt);
415  float cut_mass_high = m_funcMassCutHigh->Eval(jet_pt);
416  float cut_score = m_funcScoreCut ->Eval(jet_pt);
417 
419  ATH_MSG_VERBOSE( "Cut values : Mass window = [" << cut_mass_low << "," << cut_mass_high << "], score cut = " << cut_score );
420  ATH_MSG_VERBOSE( "Jet values : Mass = " << jet_mass << ", score = " << jet_score );
421 
423  ATH_CHECK( getWeight( jet, jet_score > cut_score, acceptData ) );
424 
426  ATH_MSG_DEBUG( "Decorating with score" );
427 
429  SG::WriteDecorHandle<xAOD::JetContainer, bool> decPassMass(m_decPassMassKey);
430  SG::WriteDecorHandle<xAOD::JetContainer, bool> decPassScore(m_decPassScoreKey);
431  SG::WriteDecorHandle<xAOD::JetContainer, bool> decTagged(m_decTaggedKey);
432  SG::WriteDecorHandle<xAOD::JetContainer, float> decCutMLow(m_decCutMLowKey);
433  SG::WriteDecorHandle<xAOD::JetContainer, float> decCutMHigh(m_decCutMHighKey);
434  SG::WriteDecorHandle<xAOD::JetContainer, float> decScoreCut(m_decScoreCutKey);
435  SG::WriteDecorHandle<xAOD::JetContainer, float> decScoreValue(m_decScoreValueKey);
436 
437 
439  decCutMLow(jet) = cut_mass_low;
440  decCutMHigh(jet) = cut_mass_high;
441  decScoreCut(jet) = cut_score;
442  decScoreValue(jet) = jet_score;
443 
445  bool passCuts = true;
446 
448  if ( m_tagClass == TAGCLASS::WBoson || m_tagClass == TAGCLASS::ZBoson ) {
449  ATH_MSG_VERBOSE( "Determining WZ tag return" );
450  if ( jet_mass > cut_mass_low ) acceptData.setCutResult( "PassMassLow", true );
451  if ( jet_mass < cut_mass_high ) acceptData.setCutResult( "PassMassHigh", true );
452  if ( jet_score > cut_score ) acceptData.setCutResult( "PassScore", true );
453  decPassMass(jet) = acceptData.getCutResult( "PassMassLow" ) && acceptData.getCutResult( "PassMassHigh" );
454  passCuts = passCuts && acceptData.getCutResult( "PassMassLow" ) && acceptData.getCutResult( "PassMassHigh" );
455  }
456  else if ( m_tagClass == TAGCLASS::TopQuark ) {
457  ATH_MSG_VERBOSE( "Determining TopQuark tag return" );
458  if( jet_mass > cut_mass_low ) acceptData.setCutResult( "PassMassLow", true );
459  if( jet_score > cut_score ) acceptData.setCutResult( "PassScore", true );
460  decPassMass(jet) = acceptData.getCutResult( "PassMassLow" );
461  passCuts = passCuts && acceptData.getCutResult( "PassMassLow" );
462  }
463 
464  decPassScore(jet) = acceptData.getCutResult( "PassScore" );
465 
466  passCuts = passCuts && acceptData.getCutResult( "PassScore" );
467 
468  decTagged(jet) = passCuts;
469 
470  return StatusCode::SUCCESS;
471 
472 }

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

◆ m_ANN_inputValues

std::map<std::string, std::map<std::string,double> > JSSWTopTaggerANN::m_ANN_inputValues
private

Definition at line 35 of file JSSWTopTaggerANN.h.

◆ m_calcSF

bool JSSTaggerBase::m_calcSF {}
protectedinherited

Flag to calculate scale factor.

Definition at line 192 of file JSSTaggerBase.h.

◆ m_calibArea

std::string JSSTaggerBase::m_calibArea
protectedinherited

Location where config files live on cvmfs.

Definition at line 113 of file JSSTaggerBase.h.

◆ m_configFile

std::string JSSTaggerBase::m_configFile
protectedinherited

Configuration file name.

Definition at line 110 of file JSSTaggerBase.h.

◆ m_configReader

TEnv JSSTaggerBase::m_configReader
protectedinherited

TEnv instance to read config files.

Definition at line 63 of file JSSTaggerBase.h.

◆ m_containerName

std::string JSSTaggerBase::m_containerName
protectedinherited

Configurable members.

Jet container name

Definition at line 103 of file JSSTaggerBase.h.

◆ m_decC2Key

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

Definition at line 145 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 219 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 218 of file JSSTaggerBase.h.

◆ m_decD2Key

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

Definition at line 146 of file JSSTaggerBase.h.

◆ m_decE3Key

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

Definition at line 147 of file JSSTaggerBase.h.

◆ m_decEfficiencyKey

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

Definition at line 225 of file JSSTaggerBase.h.

◆ m_decEffSFKey

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

Definition at line 226 of file JSSTaggerBase.h.

◆ m_decL2Key

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

Definition at line 148 of file JSSTaggerBase.h.

◆ m_decL3Key

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

Definition at line 149 of file JSSTaggerBase.h.

◆ m_decorationName

std::string JSSTaggerBase::m_decorationName
protectedinherited

Decoration name.

Definition at line 189 of file JSSTaggerBase.h.

◆ m_decPassMassKey

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

Definition at line 77 of file JSSTaggerBase.h.

◆ m_decPassScoreKey

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

Definition at line 78 of file JSSTaggerBase.h.

◆ m_decScoreCutKey

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

Definition at line 220 of file JSSTaggerBase.h.

◆ m_decScoreValueKey

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

Definition at line 221 of file JSSTaggerBase.h.

◆ m_decSigeffSFKey

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

Definition at line 227 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 69 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 141 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 142 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 143 of file JSSTaggerBase.h.

◆ m_decValidEtaRangeKey

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

Definition at line 72 of file JSSTaggerBase.h.

◆ m_decValidEventContentKey

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

Definition at line 75 of file JSSTaggerBase.h.

◆ m_decValidJetContentKey

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

Definition at line 74 of file JSSTaggerBase.h.

◆ m_decValidKinRangeKey

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

Definition at line 73 of file JSSTaggerBase.h.

◆ m_decValidPtRangeHighKey

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

Definition at line 70 of file JSSTaggerBase.h.

◆ m_decValidPtRangeLowKey

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

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

◆ m_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_efficiencyHistogramName

std::string JSSTaggerBase::m_efficiencyHistogramName
protectedinherited

Definition at line 199 of file JSSTaggerBase.h.

◆ m_efficiencyHistograms

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

Definition at line 205 of file JSSTaggerBase.h.

◆ m_evtStore

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

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_funcMassCutHigh

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

Definition at line 185 of file JSSTaggerBase.h.

◆ m_funcMassCutLow

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

TF1 for cut functions.

Definition at line 184 of file JSSTaggerBase.h.

◆ m_funcScoreCut

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

Definition at line 186 of file JSSTaggerBase.h.

◆ m_isMC

bool JSSTaggerBase::m_isMC {}
protectedinherited

Definition at line 193 of file JSSTaggerBase.h.

◆ m_isSmallRJet

bool JSSTaggerBase::m_isSmallRJet = false
protectedinherited

Definition at line 104 of file JSSTaggerBase.h.

◆ m_jetEtaMax

float JSSTaggerBase::m_jetEtaMax {}
protectedinherited

Definition at line 133 of file JSSTaggerBase.h.

◆ m_jetPtMax

float JSSTaggerBase::m_jetPtMax {}
protectedinherited

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

◆ m_kerasCalibArea

std::string JSSTaggerBase::m_kerasCalibArea
protectedinherited

Definition at line 119 of file JSSTaggerBase.h.

◆ m_kerasConfigFileName

std::string JSSTaggerBase::m_kerasConfigFileName
protectedinherited

Keras configurations for ML taggers.

Definition at line 116 of file JSSTaggerBase.h.

◆ m_kerasConfigFilePath

std::string JSSTaggerBase::m_kerasConfigFilePath
protectedinherited

Definition at line 117 of file JSSTaggerBase.h.

◆ m_kerasConfigOutputName

std::string JSSTaggerBase::m_kerasConfigOutputName
protectedinherited

Definition at line 118 of file JSSTaggerBase.h.

◆ m_lwnn

std::unique_ptr<lwt::LightweightGraph> JSSWTopTaggerANN::m_lwnn
private

ANN tools.

Definition at line 34 of file JSSWTopTaggerANN.h.

◆ m_nWarnMax

const int JSSTaggerBase::m_nWarnMax = 10
protectedinherited

Maximum number of warnings.

Definition at line 81 of file JSSTaggerBase.h.

◆ m_nWarnVar

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

Warning counters.

Definition at line 84 of file JSSTaggerBase.h.

◆ m_out_names

std::vector<std::string> JSSWTopTaggerANN::m_out_names
private

Internal stuff to keep track of the output node for the NN.

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

◆ m_readECF1Key

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

Definition at line 157 of file JSSTaggerBase.h.

◆ m_readECF2Key

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

Definition at line 158 of file JSSTaggerBase.h.

◆ m_readECF3Key

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

Definition at line 159 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 171 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 170 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 169 of file JSSTaggerBase.h.

◆ m_readParentKey

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

Definition at line 173 of file JSSTaggerBase.h.

◆ m_readQwKey

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

Definition at line 164 of file JSSTaggerBase.h.

◆ m_readSphericityKey

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

Definition at line 167 of file JSSTaggerBase.h.

◆ m_readSplit12Key

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

Definition at line 161 of file JSSTaggerBase.h.

◆ m_readSplit23Key

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

Definition at line 162 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 152 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 153 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 154 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 155 of file JSSTaggerBase.h.

◆ m_readThrustMajKey

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

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

◆ m_strMassCutHigh

std::string JSSTaggerBase::m_strMassCutHigh
protectedinherited

Definition at line 180 of file JSSTaggerBase.h.

◆ m_strMassCutLow

std::string JSSTaggerBase::m_strMassCutLow
protectedinherited

Strings for cut functions.

Definition at line 179 of file JSSTaggerBase.h.

◆ m_strScoreCut

std::string JSSTaggerBase::m_strScoreCut
protectedinherited

Definition at line 181 of file JSSTaggerBase.h.

◆ m_suppressInputDependence

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

Definition at line 175 of file JSSTaggerBase.h.

◆ m_suppressOutputDependence

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

Definition at line 176 of file JSSTaggerBase.h.

◆ m_tagClass

TAGCLASS JSSTaggerBase::m_tagClass {Unknown}
protectedinherited

Definition at line 98 of file JSSTaggerBase.h.

◆ m_tagType

std::string JSSTaggerBase::m_tagType
protectedinherited

Definition at line 128 of file JSSTaggerBase.h.

◆ m_tmvaCalibArea

std::string JSSTaggerBase::m_tmvaCalibArea
protectedinherited

Definition at line 124 of file JSSTaggerBase.h.

◆ m_tmvaConfigFileName

std::string JSSTaggerBase::m_tmvaConfigFileName
protectedinherited

TMVA configurations for BDT taggers.

Definition at line 122 of file JSSTaggerBase.h.

◆ m_tmvaConfigFilePath

std::string JSSTaggerBase::m_tmvaConfigFilePath
protectedinherited

Definition at line 123 of file JSSTaggerBase.h.

◆ m_truthBosonContainerName

std::string JSSTaggerBase::m_truthBosonContainerName
protectedinherited

Definition at line 210 of file JSSTaggerBase.h.

◆ m_truthLabelName

std::string JSSTaggerBase::m_truthLabelName
protectedinherited

Definition at line 212 of file JSSTaggerBase.h.

◆ m_truthLabelUseTRUTH3

bool JSSTaggerBase::m_truthLabelUseTRUTH3 {}
protectedinherited

Truth label options.

Definition at line 208 of file JSSTaggerBase.h.

◆ m_truthParticleContainerName

std::string JSSTaggerBase::m_truthParticleContainerName
protectedinherited

Definition at line 209 of file JSSTaggerBase.h.

◆ m_truthTopQuarkContainerName

std::string JSSTaggerBase::m_truthTopQuarkContainerName
protectedinherited

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

◆ m_useLSeries

bool JSSTaggerBase::m_useLSeries = false
protectedinherited

Definition at line 138 of file JSSTaggerBase.h.

◆ m_useMassCut

bool JSSTaggerBase::m_useMassCut = false
protectedinherited

Flag to indicate if mass window cut is used.

Definition at line 91 of file JSSTaggerBase.h.

◆ m_useScoreCut

bool JSSTaggerBase::m_useScoreCut = false
protectedinherited

Flag to indicate if a discriminant score is used.

Definition at line 94 of file JSSTaggerBase.h.

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.

◆ m_weightConfig

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

Histograms for scale factors.

Definition at line 203 of file JSSTaggerBase.h.

◆ m_weightConfigPath

std::string JSSTaggerBase::m_weightConfigPath
protectedinherited

Path to the SF configuration root file.

Definition at line 107 of file JSSTaggerBase.h.

◆ m_weightDecorationName

std::string JSSTaggerBase::m_weightDecorationName
protectedinherited

String for scale factor decoration names.

Definition at line 196 of file JSSTaggerBase.h.

◆ m_weightFileName

std::string JSSTaggerBase::m_weightFileName
protectedinherited

Definition at line 197 of file JSSTaggerBase.h.

◆ m_weightFlavors

std::string JSSTaggerBase::m_weightFlavors
protectedinherited

Definition at line 200 of file JSSTaggerBase.h.

◆ m_weightHistogramName

std::string JSSTaggerBase::m_weightHistogramName
protectedinherited

Definition at line 198 of file JSSTaggerBase.h.

◆ m_weightHistograms

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

Definition at line 204 of file JSSTaggerBase.h.

◆ m_wkpt

std::string JSSTaggerBase::m_wkpt
protectedinherited

Tagger information.

Definition at line 127 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:173
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
LargeRJetTruthLabel::TypeEnum
TypeEnum
Definition: LargeRJetLabelEnum.h:14
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
JSSTaggerBase::m_readTau4WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau4WTAKey
Definition: JSSTaggerBase.h:155
JSSTaggerBase::m_readSplit23Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSplit23Key
Definition: JSSTaggerBase.h:162
Trk::L2
@ L2
Definition: AlignModuleList.h:32
JSSTaggerBase::m_readECF1Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF1Key
Definition: JSSTaggerBase.h:157
JSSTaggerBase::getTruthLabelStr
std::string getTruthLabelStr(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Definition: JSSTaggerBase.cxx:837
JSSTaggerBase::m_decScoreCutKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decScoreCutKey
Definition: JSSTaggerBase.h:220
JSSTaggerBase::m_strMassCutLow
std::string m_strMassCutLow
Strings for cut functions.
Definition: JSSTaggerBase.h:179
get_generator_info.result
result
Definition: get_generator_info.py:21
StateLessPT_NewConfig.proxy
proxy
Definition: StateLessPT_NewConfig.py:407
JSSTaggerBase::m_funcScoreCut
std::unique_ptr< TF1 > m_funcScoreCut
Definition: JSSTaggerBase.h:186
JSSTaggerBase::ZBoson
@ ZBoson
Definition: JSSTaggerBase.h:97
JSSTaggerBase::resetCuts
StatusCode resetCuts(asg::AcceptData &acceptData) const
Reset cuts.
Definition: JSSTaggerBase.cxx:334
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
JSSTaggerBase::m_readECFG311Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG311Key
Definition: JSSTaggerBase.h:170
FlavorTagInference::SaltModelGraphConfig::parse_json_graph
GraphConfig parse_json_graph(const nlohmann::json &metadata)
Definition: SaltModelGraphConfig.cxx:40
JSSTaggerBase::m_jetEtaMax
float m_jetEtaMax
Definition: JSSTaggerBase.h:133
JSSTaggerBase::m_nWarnVar
std::atomic< int > m_nWarnVar
Warning counters.
Definition: JSSTaggerBase.h:84
JSSTaggerBase::m_useMassCut
bool m_useMassCut
Flag to indicate if mass window cut is used.
Definition: JSSTaggerBase.h:91
JSSTaggerBase::m_decValidEtaRangeKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidEtaRangeKey
Definition: JSSTaggerBase.h:72
JSSTaggerBase::m_kerasCalibArea
std::string m_kerasCalibArea
Definition: JSSTaggerBase.h:119
JSSTaggerBase::m_efficiencyHistograms
std::map< std::string, std::unique_ptr< TH2D > > m_efficiencyHistograms
Definition: JSSTaggerBase.h:205
JSSTaggerBase::m_configFile
std::string m_configFile
Configuration file name.
Definition: JSSTaggerBase.h:110
JSSTaggerBase::m_useScoreCut
bool m_useScoreCut
Flag to indicate if a discriminant score is used.
Definition: JSSTaggerBase.h:94
JSSTaggerBase::m_configReader
TEnv m_configReader
TEnv instance to read config files.
Definition: JSSTaggerBase.h:63
RootHelpers::FindBin
Int_t FindBin(const TAxis *axis, const double x)
Definition: RootHelpers.cxx:14
JSSTaggerBase::TopQuark
@ TopQuark
Definition: JSSTaggerBase.h:97
JSSTaggerBase::m_nWarnMax
const int m_nWarnMax
Maximum number of warnings.
Definition: JSSTaggerBase.h:81
JSSTaggerBase::m_weightFileName
std::string m_weightFileName
Definition: JSSTaggerBase.h:197
JSSTaggerBase::m_decPassScoreKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassScoreKey
Definition: JSSTaggerBase.h:78
JSSTaggerBase::getWeight
StatusCode getWeight(const xAOD::Jet &jet, bool passSel, asg::AcceptData &acceptData) const
Get SF weight.
Definition: JSSTaggerBase.cxx:705
python.SystemOfUnits.second
float second
Definition: SystemOfUnits.py:135
TRT_PAI_gasdata::SF
const float SF[NF]
Cross sections for Fluor.
Definition: TRT_PAI_gasdata.h:285
LargeRJetTruthLabel::other_From_t
@ other_From_t
Definition: LargeRJetLabelEnum.h:22
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
defineDB.jets
jets
Definition: JetTagCalibration/share/defineDB.py:24
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
JSSTaggerBase::m_readTau2WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau2WTAKey
Definition: JSSTaggerBase.h:153
JSSTaggerBase::m_readSplit12Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSplit12Key
Definition: JSSTaggerBase.h:161
JSSTaggerBase::m_jetPtMin
float m_jetPtMin
Kinematic bounds for the jet - the units are controlled by m_ptGeV.
Definition: JSSTaggerBase.h:131
JSSTaggerBase::m_decL3Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL3Key
Definition: JSSTaggerBase.h:149
JSSWTopTaggerANN::m_lwnn
std::unique_ptr< lwt::LightweightGraph > m_lwnn
ANN tools.
Definition: JSSWTopTaggerANN.h:34
SG::ConstAccessor
Helper class to provide constant type-safe access to aux data.
Definition: ConstAccessor.h:55
JSSTaggerBase::m_decCutMLowKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutMLowKey
WriteDecorHandle keys for cut values.
Definition: JSSTaggerBase.h:218
JSSTaggerBase::m_decC2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decC2Key
Definition: JSSTaggerBase.h:145
JSSTaggerBase::m_decorationName
std::string m_decorationName
Decoration name.
Definition: JSSTaggerBase.h:189
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
JSSTaggerBase::m_readECFG331Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG331Key
Definition: JSSTaggerBase.h:169
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
ZBoson
@ ZBoson
Definition: TruthClasses.h:67
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:70
JSSTaggerBase::m_funcMassCutHigh
std::unique_ptr< TF1 > m_funcMassCutHigh
Definition: JSSTaggerBase.h:185
JSSTaggerBase::m_weightHistogramName
std::string m_weightHistogramName
Definition: JSSTaggerBase.h:198
config
Definition: PhysicsAnalysis/AnalysisCommon/AssociationUtils/python/config.py:1
JSSWTopTaggerANN::getScore
double getScore(const xAOD::Jet &jet) const
Retrieve score for a given ANN type (top/W)
Definition: JSSWTopTaggerANN.cxx:474
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:132
JSSTaggerBase::m_calibArea
std::string m_calibArea
Location where config files live on cvmfs.
Definition: JSSTaggerBase.h:113
JSSTaggerBase::m_decWeightKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decWeightKey
WriteDecorHandle keys for SF.
Definition: JSSTaggerBase.h:224
JSSTaggerBase::m_readTruthLabelKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTruthLabelKey
Truth label ReadDecorHandle key.
Definition: JSSTaggerBase.h:215
LargeRJetTruthLabel::notruth
@ notruth
Definition: LargeRJetLabelEnum.h:24
WBoson
@ WBoson
Definition: TruthClasses.h:66
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
JSSTaggerBase::m_decSigeffSFKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decSigeffSFKey
Definition: JSSTaggerBase.h:227
JSSWTopTaggerANN::getJetProperties
std::map< std::string, std::map< std::string, double > > getJetProperties(const xAOD::Jet &jet) const
Update the jet substructure variables for each jet to use in ANN.
Definition: JSSWTopTaggerANN.cxx:508
dqt_zlumi_pandas.weight
int weight
Definition: dqt_zlumi_pandas.py:190
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
JSSTaggerBase::m_truthLabelName
std::string m_truthLabelName
Definition: JSSTaggerBase.h:212
JSSTaggerBase::m_decEfficiencyKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decEfficiencyKey
Definition: JSSTaggerBase.h:225
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:85
JSSTaggerBase::m_decE3Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decE3Key
Definition: JSSTaggerBase.h:147
JSSTaggerBase::m_decCutMHighKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutMHighKey
Definition: JSSTaggerBase.h:219
JSSTaggerBase::m_ptGeV
bool m_ptGeV
Flag to indicate units pT is defined in Set to false by default.
Definition: JSSTaggerBase.h:88
efficiency
void efficiency(std::vector< double > &bins, std::vector< double > &values, const std::vector< std::string > &files, const std::string &histname, const std::string &tplotname, const std::string &label="")
Definition: dependence.cxx:128
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:209
JSSTaggerBase::m_decTau32WTAKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau32WTAKey
Definition: JSSTaggerBase.h:142
jet
Definition: JetCalibTools_PlotJESFactors.cxx:23
AthCommonDataStore
Definition: AthCommonDataStore.h:52
LargeRJetTruthLabel::Wqq
@ Wqq
Definition: LargeRJetLabelEnum.h:17
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
CheckAppliedSFs.e3
e3
Definition: CheckAppliedSFs.py:264
JSSTaggerBase::m_readQwKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readQwKey
Definition: JSSTaggerBase.h:164
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:779
JSSTaggerBase::m_tagClass
TAGCLASS m_tagClass
Definition: JSSTaggerBase.h:98
JSSTaggerBase::m_kerasConfigOutputName
std::string m_kerasConfigOutputName
Definition: JSSTaggerBase.h:118
Trk::L3
@ L3
Definition: AlignModuleList.h:32
JSSTaggerBase::m_readECFG212Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG212Key
Definition: JSSTaggerBase.h:171
asg::AcceptInfo::getNCuts
unsigned int getNCuts() const
Get the number of cuts defined.
Definition: AcceptInfo.h:46
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition: AthCommonDataStore.h:145
JSSTaggerBase::checkKinRange
StatusCode checkKinRange(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Check and record if jet passes kinematic constraints.
Definition: JSSTaggerBase.cxx:366
SG::WriteDecorHandle
Handle class for adding a decoration to an object.
Definition: StoreGate/StoreGate/WriteDecorHandle.h:100
PlotPulseshapeFromCool.input
input
Definition: PlotPulseshapeFromCool.py:106
xAOD::VxType::PriVtx
@ PriVtx
Primary vertex.
Definition: TrackingPrimitives.h:572
CalibDbCompareRT.dummy
dummy
Definition: CalibDbCompareRT.py:59
JSSTaggerBase::Unknown
@ Unknown
Definition: JSSTaggerBase.h:97
JSSTaggerBase::m_efficiencyHistogramName
std::string m_efficiencyHistogramName
Definition: JSSTaggerBase.h:199
JSSTaggerBase::m_strScoreCut
std::string m_strScoreCut
Definition: JSSTaggerBase.h:181
JSSTaggerBase::m_decTau42WTAKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau42WTAKey
Definition: JSSTaggerBase.h:143
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:412
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:226
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
JSSTaggerBase::m_calcSF
bool m_calcSF
Flag to calculate scale factor.
Definition: JSSTaggerBase.h:192
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
JSSTaggerBase::findPV
int findPV() const
Find the PV (to be used for Ntrk)
Definition: JSSTaggerBase.cxx:643
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_kerasConfigFilePath
std::string m_kerasConfigFilePath
Definition: JSSTaggerBase.h:117
JSSTaggerBase::m_decD2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decD2Key
Definition: JSSTaggerBase.h:146
LargeRJetTruthLabel::Wqq_From_t
@ Wqq_From_t
Definition: LargeRJetLabelEnum.h:21
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
TauJetParameters::discriminant
@ discriminant
Definition: TauJetParameters.h:166
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
JSSTaggerBase::initialize
virtual StatusCode initialize() override
Initialize the tool.
Definition: JSSTaggerBase.cxx:72
JSSTaggerBase::m_acceptInfo
asg::AcceptInfo m_acceptInfo
Object that stores the results for a jet.
Definition: JSSTaggerBase.h:66
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:141
JSSTaggerBase::decorateJSSRatios
void decorateJSSRatios(const xAOD::JetContainer &jets) const
Definition: JSSTaggerBase.cxx:520
PathResolverFindCalibFile
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
Definition: PathResolver.cxx:283
JSSTaggerBase::m_decScoreValueKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decScoreValueKey
Definition: JSSTaggerBase.h:221
JSSTaggerBase::printCuts
void printCuts() const
Print configured cuts.
Definition: JSSTaggerBase.cxx:926
JSSTaggerBase::m_isMC
bool m_isMC
Definition: JSSTaggerBase.h:193
JSSTaggerBase::m_weightDecorationName
std::string m_weightDecorationName
String for scale factor decoration names.
Definition: JSSTaggerBase.h:196
JSSTaggerBase::WBoson
@ WBoson
Definition: JSSTaggerBase.h:97
JSSWTopTaggerANN::m_out_names
std::vector< std::string > m_out_names
Internal stuff to keep track of the output node for the NN.
Definition: JSSWTopTaggerANN.h:38
xAOD::Jet_v1
Class describing a jet.
Definition: Jet_v1.h:57
JSSTaggerBase::m_decValidKinRangeKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidKinRangeKey
Definition: JSSTaggerBase.h:73
JSSTaggerBase::m_readTau1WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau1WTAKey
ReadDecorHandle keys for JSS moments.
Definition: JSSTaggerBase.h:152
jet::CompScaleVar::Tau21WTA
@ Tau21WTA
Definition: UncertaintyEnum.h:103
JSSTaggerBase::m_weightHistograms
std::map< std::string, std::unique_ptr< TH2D > > m_weightHistograms
Definition: JSSTaggerBase.h:204
a
TList * a
Definition: liststreamerinfos.cxx:10
asg::AcceptData::setCutResult
void setCutResult(const std::string &cutName, bool cutResult)
Set the result of a cut, based on the cut name (safer)
Definition: AcceptData.h:134
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:352
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:71
JSSTaggerBase::m_readTau3WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau3WTAKey
Definition: JSSTaggerBase.h:154
JSSTaggerBase::m_readECF2Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF2Key
Definition: JSSTaggerBase.h:158
JSSTaggerBase::m_readECF3Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF3Key
Definition: JSSTaggerBase.h:159
AthCommonMsg< AlgTool >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
JSSTaggerBase::getConfigReader
StatusCode getConfigReader()
Get configReader StatusCode.
Definition: JSSTaggerBase.cxx:295
LargeRJetTruthLabel::intToEnum
TypeEnum intToEnum(const int type)
Definition: LargeRJetLabelEnum.h:57
JSSTaggerBase::m_decPassMassKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassMassKey
Definition: JSSTaggerBase.h:77
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
JSSTaggerBase::m_isSmallRJet
bool m_isSmallRJet
Definition: JSSTaggerBase.h:104
JSSTaggerBase::GetUnGroomTracks
int GetUnGroomTracks(const xAOD::Jet &jet, int indexPV) const
Retrieve Ntrk variable from the ungroomed parent jet.
Definition: JSSTaggerBase.cxx:667
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:629
jet::CompScaleVar::Tau32WTA
@ Tau32WTA
Definition: UncertaintyEnum.h:104
JSSTaggerBase::m_tagType
std::string m_tagType
Definition: JSSTaggerBase.h:128
JSSTaggerBase::JSSTaggerBase
JSSTaggerBase(const std::string &name)
Default constructor - to be used in all derived classes.
Definition: JSSTaggerBase.cxx:11
JSSTaggerBase::m_funcMassCutLow
std::unique_ptr< TF1 > m_funcMassCutLow
TF1 for cut functions.
Definition: JSSTaggerBase.h:184
dqt_zlumi_alleff_HIST.eff
int eff
Definition: dqt_zlumi_alleff_HIST.py:113
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:801
LargeRJetTruthLabel::tqqb
@ tqqb
Definition: LargeRJetLabelEnum.h:16
JSSTaggerBase::m_strMassCutHigh
std::string m_strMassCutHigh
Definition: JSSTaggerBase.h:180
JSSTaggerBase::m_kerasConfigFileName
std::string m_kerasConfigFileName
Keras configurations for ML taggers.
Definition: JSSTaggerBase.h:116
JSSTaggerBase::m_decL2Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL2Key
Definition: JSSTaggerBase.h:148
JSSTaggerBase::m_weightConfigPath
std::string m_weightConfigPath
Path to the SF configuration root file.
Definition: JSSTaggerBase.h:107
JSSTaggerBase::m_weightFlavors
std::string m_weightFlavors
Definition: JSSTaggerBase.h:200
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
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:69
asg::AcceptInfo::addCut
int addCut(const std::string &cutName, const std::string &cutDescription)
Add a cut; returning the cut position.
Definition: AcceptInfo.h:53
Generated on Thu Sep 4 2025 21:38:20 for ATLAS Offline Software by   doxygen 1.8.18