ATLAS Offline Software
Loading...
Searching...
No Matches
SmoothedWZTagger Class Reference

#include <SmoothedWZTagger.h>

Inheritance diagram for SmoothedWZTagger:
Collaboration diagram for SmoothedWZTagger:

Public Member Functions

 SmoothedWZTagger (const std::string &name)
 Constructor.
virtual StatusCode initialize () override
 Run once at the start of the job to setup everything.
virtual StatusCode tag (const xAOD::Jet &jet) const override
 Decorate single jet with tagging info.
virtual StatusCode decorate (const xAOD::JetContainer &jets) const override
 Decorate jet collection with tagging info.
virtual void print () const
 Print the state of the tool.
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysInitialize () override
 Perform system initialization for an algorithm.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const
virtual StatusCode modify (xAOD::JetContainer &jets) const override final
 Concrete implementation of the function inherited from IJetModifier.
Additional helper functions, not directly mimicking Athena
template<class T>
const T * getProperty (const std::string &name) const
 Get one of the tool's properties.
const std::string & msg_level_name () const __attribute__((deprecated))
 A deprecated function for getting the message level's name.
const std::string & getName (const void *ptr) const
 Get the name of an object that is / should be in the event store.
SG::sgkey_t getKey (const void *ptr) const
 Get the (hashed) key of an object that is in the event store.

Protected Types

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

Protected Member Functions

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

Protected Attributes

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

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

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

Private Attributes

bool m_useNtrk {}
 Flag to indicate if Ntrk is used.
bool m_useScore {}
 Flag to indicate if a ML score is used.
std::string m_strD2Cut
 Store functional form of cuts.
std::string m_strNtrkCut
std::string m_strScoreCut
std::unique_ptr< TF1 > m_funcD2Cut
 Functions that are configurable for specific cut values.
std::unique_ptr< TF1 > m_funcNtrkCut
std::unique_ptr< TF1 > m_funcScoreCut
SG::WriteDecorHandleKey< xAOD::JetContainerm_decPassD2Key {this, "PassD2Name", "PassD2", "SG key for PassD2"}
 WriteDecorHandle keys.
SG::WriteDecorHandleKey< xAOD::JetContainerm_decPassNtrkKey {this, "PassNtrkName", "PassNtrk", "SG key for PassNtrk"}
SG::WriteDecorHandleKey< xAOD::JetContainerm_decPassScoreKey {this, "PassScoreName", "PassScore", "SG key for PassScore"}
SG::WriteDecorHandleKey< xAOD::JetContainerm_decCutD2Key {this, "CutD2Name", "Cut_D2", "SG key for Cut_D2"}
SG::WriteDecorHandleKey< xAOD::JetContainerm_decCutNtrkKey {this, "CutNtrkName", "Cut_Ntrk", "SG key for Cut_Ntrk"}
SG::WriteDecorHandleKey< xAOD::JetContainerm_decCutScoreKey {this, "CutScoreName", "Cut_Score", "SG key for Cut_Score"}
SG::WriteDecorHandleKey< xAOD::JetContainerm_decAcceptKey {this, "acceptName", "accept", "SG key for accept"}
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default)
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared

Detailed Description

Definition at line 10 of file SmoothedWZTagger.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.

Constructor & Destructor Documentation

◆ SmoothedWZTagger()

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

Constructor.

Functional forms for cuts

Definition at line 9 of file SmoothedWZTagger.cxx.

9 :
10 JSSTaggerBase( name )
11{
12
14 declareProperty( "MassCutLowFunc", m_strMassCutLow = "", "Lower mass cut");
15 declareProperty( "MassCutHighFunc", m_strMassCutHigh = "", "Higher mass cut");
16 declareProperty( "D2CutFunc", m_strD2Cut = "", "Upper cut on D2");
17 declareProperty( "NtrkCutFunc", m_strNtrkCut = "", "Upper cut on Ntrk");
18 declareProperty( "ScoreCutFunc", m_strScoreCut = "", "Lower cut on score");
19
20}
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
std::string m_strMassCutHigh
std::string m_strMassCutLow
Strings for cut functions.
JSSTaggerBase(const std::string &name)
Default constructor - to be used in all derived classes.
std::string m_strScoreCut
std::string m_strNtrkCut
std::string m_strD2Cut
Store functional form of cuts.

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

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

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

367 {
368
369 float scale = 1.0;
370 if ( m_ptGeV ) scale = 1.e3;
371
372 bool passKinRange = true;
373
376 if ( std::abs(jet.eta()) > m_jetEtaMax ) {
377 ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (|eta| < " << m_jetEtaMax << "). Jet eta = " << jet.eta() );
378 acceptData.setCutResult( "ValidEtaRange", false );
379 passKinRange = false;
380 }
381
382 if ( jet.pt() < m_jetPtMin * scale ) {
383 ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (pT > " << m_jetPtMin * scale / 1.e3 << "). Jet pT = " << jet.pt() / 1.e3 << " GeV" );
384 acceptData.setCutResult( "ValidPtRangeLow", false );
385 passKinRange = false;
386 }
387
388 if ( jet.pt() > m_jetPtMax * scale ) {
389 ATH_MSG_VERBOSE( "Jet does not pass basic kinematic selection (pT < " << m_jetPtMax * scale / 1.e3 << "). Jet pT = " << jet.pt() / 1.e3 << " GeV" );
390 acceptData.setCutResult( "ValidPtRangeHigh", false );
391 passKinRange = false;
392 }
393
395 SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidPtRangeHigh(m_decValidPtRangeHighKey);
396 SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidPtRangeLow(m_decValidPtRangeLowKey);
397 SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidEtaRange(m_decValidEtaRangeKey);
398 SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidKinRange(m_decValidKinRangeKey);
399
401 decValidPtRangeHigh(jet) = acceptData.getCutResult( "ValidPtRangeHigh" );
402 decValidPtRangeLow(jet) = acceptData.getCutResult( "ValidPtRangeLow" );
403 decValidEtaRange(jet) = acceptData.getCutResult( "ValidEtaRange" );
404 decValidKinRange(jet) = passKinRange;
405
406 return StatusCode::SUCCESS;
407
408}
#define ATH_MSG_VERBOSE(x)
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeHighKey
float m_jetPtMin
Kinematic bounds for the jet - the units are controlled by m_ptGeV.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeLowKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidKinRangeKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidEtaRangeKey
bool passKinRange(const xAOD::Jet &jet) const
Check if jet passes kinematic constraints.
bool m_ptGeV
Flag to indicate units pT is defined in Set to false by default.
void setCutResult(const std::string &cutName, bool cutResult)
Set the result of a cut, based on the cut name (safer)
Definition AcceptData.h:134
bool getCutResult(const std::string &cutName) const
Get the result of a cut, based on the cut name (safer)
Definition AcceptData.h:98
virtual double pt() const
The transverse momentum ( ) of the particle.
Definition Jet_v1.cxx:44
virtual double eta() const
The pseudorapidity ( ) of the particle.
Definition Jet_v1.cxx:49

◆ declareGaudiProperty()

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

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

Definition at line 156 of file AthCommonDataStore.h.

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }

◆ declareProperty()

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

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

◆ decorate()

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

Decorate jet collection with tagging info.

Loop over jet collection and decorate each jet.

Implements IJetDecorator.

Reimplemented in JSSWTopTaggerANN, and JSSWTopTaggerDNN.

Definition at line 285 of file JSSTaggerBase.cxx.

285 {
286
287 for ( const auto *jet : jets ) {
288 ATH_CHECK( tag(*jet) );
289 }
290
291 return StatusCode::SUCCESS;
292
293}
#define ATH_CHECK
Evaluate an expression and check for errors.
virtual StatusCode tag(const xAOD::Jet &jet) const =0
Decorate single jet with tagging info.

◆ decorateJSSRatios()

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

Create write decor handles

WTA N-subjettiness ratios

ECF ratios

Definition at line 521 of file JSSTaggerBase.cxx.

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

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

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

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

644 {
645
646 int indexPV = -1;
647
648 const xAOD::VertexContainer* vxCont = nullptr;
649 if ( evtStore()->retrieve( vxCont, "PrimaryVertices" ).isFailure() ) {
650 ATH_MSG_WARNING( "Unable to retrieve primary vertex container PrimaryVertices" );
651 }
652 else {
653 int vtx_counter = 0;
654 for ( const auto *vx : *vxCont ) {
655 if ( vx->vertexType()==xAOD::VxType::PriVtx ) {
656 indexPV = vtx_counter;
657 break;
658 }
659 vtx_counter++;
660 }
661 }
662
663 return indexPV;
664
665}
#define ATH_MSG_WARNING(x)
ServiceHandle< StoreGateSvc > & evtStore()
retrieve(aClass, aKey=None)
Definition PyKernel.py:110
@ PriVtx
Primary vertex.
VertexContainer_v1 VertexContainer
Definition of the current "Vertex container version".

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

296 {
297
298 ATH_MSG_INFO( "Using config file : " << m_configFile );
299
301 std::string configPath;
302
303 if ( m_calibArea.compare("Local") == 0 ) {
305 }
306 else if ( m_calibArea.find("eos") != std::string::npos) {
307 configPath = PathResolverFindCalibFile((m_calibArea+"/"+m_configFile).c_str());
308 configPath = (m_calibArea+"/"+m_configFile).c_str();
309 }
310 else {
311 configPath = PathResolverFindCalibFile(("BoostedJetTaggers/"+m_calibArea+"/"+m_configFile).c_str());
312 }
313
315 FileStat_t fStats;
316 int fSuccess = gSystem->GetPathInfo(configPath.c_str(), fStats);
317 if ( fSuccess ) {
318 ATH_MSG_ERROR( "Recommendations file " << m_configFile << " could not be found" );
319 return StatusCode::FAILURE;
320 }
321 else {
322 ATH_MSG_DEBUG( "Recommendations file was found : " << configPath );
323 }
324
325 if ( m_configReader.ReadFile( configPath.c_str(), EEnvLevel(0) ) ) {
326 ATH_MSG_ERROR( "Error while reading config file : "<< configPath );
327 return StatusCode::FAILURE;
328 }
329
330 return StatusCode::SUCCESS;
331
332}
#define ATH_MSG_ERROR(x)
#define ATH_MSG_INFO(x)
#define ATH_MSG_DEBUG(x)
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
std::string m_calibArea
Location where config files live on cvmfs.
std::string m_configFile
Configuration file name.
TEnv m_configReader
TEnv instance to read config files.

◆ getKey()

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

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

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

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

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

Definition at line 119 of file AsgTool.cxx.

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

◆ getName()

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

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

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

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

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

Definition at line 106 of file AsgTool.cxx.

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

◆ getProperty()

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

Get one of the tool's properties.

◆ getSF()

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

Get scale factor and efficiency.

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

Definition at line 780 of file JSSTaggerBase.cxx.

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

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

838 {
840 std::string truthLabelStr;
841
843 SG::ReadDecorHandle<xAOD::JetContainer, int> readTruthLabel(m_readTruthLabelKey);
844 LargeRJetTruthLabel::TypeEnum jetContainment = LargeRJetTruthLabel::intToEnum(readTruthLabel(jet));
845
847 if ( m_weightHistograms.count("t_qqb") ) {
848
850 if ( jetContainment==LargeRJetTruthLabel::tqqb ) {
851 truthLabelStr = "t_qqb";
852 }
854 else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
855 truthLabelStr = "q";
856 }
857
858 }
860 else if ( m_weightHistograms.count("V_qq_passMpassD2") ) {
861
863 if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t ) {
864 truthLabelStr = "t_";
865 }
867 else if ( jetContainment==LargeRJetTruthLabel::Wqq || jetContainment==LargeRJetTruthLabel::Zqq || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
868 truthLabelStr = "V_qq_";
869 }
871 else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
872 truthLabelStr = "q_";
873 }
874
876 if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") && acceptData.getCutResult("PassD2") ) {
877 truthLabelStr += "passMpassD2";
878 }
880 else if ( !(acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh")) && acceptData.getCutResult("PassD2") ) {
881 truthLabelStr += "failMpassD2";
882 }
884 else if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") && !acceptData.getCutResult("PassD2") ) {
885 truthLabelStr += "passMfailD2";
886 }
888 else{
889 truthLabelStr += "failMfailD2";
890 }
891
892 }
893
895 else if ( m_weightHistograms.count("V_qq") ) {
896
898 if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t ) {
899 truthLabelStr = "t";
900 }
902 else if ( jetContainment==LargeRJetTruthLabel::Wqq || jetContainment==LargeRJetTruthLabel::Zqq || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
903 truthLabelStr = "V_qq";
904 }
906 else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
907 truthLabelStr = "q";
908 }
909 }
910
911 // inclusive top tagger
912 else {
914 if ( jetContainment==LargeRJetTruthLabel::tqqb || jetContainment==LargeRJetTruthLabel::other_From_t || jetContainment==LargeRJetTruthLabel::Wqq_From_t ) {
915 truthLabelStr = "t";
916 }
918 else if ( jetContainment==LargeRJetTruthLabel::notruth || jetContainment==LargeRJetTruthLabel::qcd ) {
919 truthLabelStr = "q";
920 }
921
922 }
923
924 return truthLabelStr;
925}
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTruthLabelKey
Truth label ReadDecorHandle key.
TypeEnum intToEnum(const int type)

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

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

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

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

◆ initialize()

StatusCode SmoothedWZTagger::initialize ( void )
overridevirtual

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

Pt values are defined in GeV

Use mass cut

Get configReader

Get min and max jet pt

Get the decoration name

Set flag to indicate if a ML score cut is used

Get the scale factor configuration

Get truth label name information

No config file

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

Set flag to indicate if Ntrk cut is used

Transform these strings into functions

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

Initialize additional decorators

Reimplemented from JSSTaggerBase.

Definition at line 22 of file SmoothedWZTagger.cxx.

22 {
23
24 ATH_MSG_INFO( "Initializing SmoothedWZTagger tool" );
25
27 m_ptGeV = true;
28
30 m_useMassCut = true;
31
32 if ( ! m_configFile.empty() ) {
33
36
37 if ( m_wkpt.empty() ) {
38 m_strMassCutLow = m_configReader.GetValue("MassCutLow", "");
39 m_strMassCutHigh = m_configReader.GetValue("MassCutHigh", "");
40 m_strD2Cut = m_configReader.GetValue("D2Cut", "");
41 m_strNtrkCut = m_configReader.GetValue("NtrkCut", "");
42 m_strScoreCut = m_configReader.GetValue("ScoreCut", "");
43 }
44 else {
45 m_strMassCutLow = m_configReader.GetValue((m_wkpt+".MassCutLow").c_str(), "");
46 m_strMassCutHigh = m_configReader.GetValue((m_wkpt+".MassCutHigh").c_str(), "");
47 m_strD2Cut = m_configReader.GetValue((m_wkpt+".D2Cut").c_str(), "");
48 m_strNtrkCut = m_configReader.GetValue((m_wkpt+".NtrkCut").c_str(), "");
49 m_strScoreCut = m_configReader.GetValue((m_wkpt+".ScoreCut").c_str(), "");
50 }
51
53 m_jetPtMin = m_configReader.GetValue("pTCutLow", 200.0);
54 m_jetPtMax = m_configReader.GetValue("pTCutHigh", 4000.0);
55
57 m_decorationName = m_configReader.GetValue("DecorationName", "");
58
60 if (!m_strScoreCut.empty())
61 m_useScore = true;
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 if ( m_calibArea.compare("Local") == 0 ) {
76 m_weightConfigPath = PathResolverFindCalibFile(("$WorkDir_DIR/data/BoostedJetTaggers/SmoothedWZTaggers/"+m_weightFileName).c_str());
77 }
78 else if ( m_calibArea.find("eos") != std::string::npos ) {
80 }
81 else {
83 }
84 }
85
86 }
87 else {
90 if( m_strD2Cut.empty() || m_strMassCutLow.empty() || m_strMassCutHigh.empty() ||
91 ((m_weightDecorationName.empty() ||
92 m_weightHistogramName.empty() ||
93 m_weightFlavors.empty()) && m_calcSF) ) {
94 ATH_MSG_ERROR( "No config file provided AND no parameters specified." ) ;
95 return StatusCode::FAILURE;
96 }
97 }
98
100 m_useNtrk = !m_strNtrkCut.empty();
101
103 m_funcD2Cut = std::make_unique<TF1>("strD2Cut", m_strD2Cut.c_str(), 0, 14000);
104 if ( m_useNtrk ) m_funcNtrkCut = std::make_unique<TF1>("strNtrkCut", m_strNtrkCut.c_str(), 0, 14000);
105 if(m_useScore)
106 m_funcScoreCut = std::make_unique<TF1>("strScoreCut", m_strScoreCut.c_str(), 0, 14000);
107
108 ATH_MSG_INFO( "Smoothed WZ Tagger tool initialized" );
109 ATH_MSG_INFO( " Mass cut low : " << m_strMassCutLow );
110 ATH_MSG_INFO( " Mass cut High : " << m_strMassCutHigh );
111 ATH_MSG_INFO( " D2 cut high : " << m_strD2Cut );
112 if ( m_useNtrk )
113 ATH_MSG_INFO( " Ntrk cut high : " << m_strNtrkCut );
114 if( m_useScore)
115 ATH_MSG_INFO( " Score cut low : " << m_strScoreCut );
116 ATH_MSG_INFO( " DecorationName : " << m_decorationName );
117 if ( m_calcSF ) {
118 ATH_MSG_INFO( "weightDecorationName : " << m_weightDecorationName );
119 ATH_MSG_INFO( "weightFile : " << m_weightFileName );
120 ATH_MSG_INFO( "weightHistogramName : " << m_weightHistogramName );
121 ATH_MSG_INFO( "efficiencyHistogramName : " << m_efficiencyHistogramName );
122 ATH_MSG_INFO( "weightFlavors : " << m_weightFlavors );
123 ATH_MSG_INFO( "TruthLabelName : " << m_truthLabelName );
124 }
125 ATH_MSG_INFO( " Pt cut low : " << m_jetPtMin );
126 ATH_MSG_INFO( " Pt cut high : " << m_jetPtMax );
127
129 m_acceptInfo.addCut( "PassMassLow", "mJet > mCutLow" );
130 m_acceptInfo.addCut( "PassMassHigh", "mJet < mCutHigh" );
131 m_acceptInfo.addCut( "PassD2", "D2Jet < D2Cut" );
132 if ( m_useNtrk ) {
133 m_acceptInfo.addCut( "PassNtrk", "NtrkJet < NtrkCut" );
134 }
135 if(m_useScore){
136 m_acceptInfo.addCut( "PassScore", "ScoreJet > ScoreCut" );
137 }
138
140 printCuts();
141
144
146 ATH_MSG_INFO( "Additional decorators that will be attached to jet :" );
147
150
151 ATH_CHECK( m_decPassD2Key.initialize() );
152 ATH_CHECK( m_decCutD2Key.initialize() );
153
154 ATH_MSG_INFO( " " << m_decPassD2Key.key() << " : pass D2 cut" );
155 ATH_MSG_INFO( " " << m_decCutD2Key.key() << " : D2 cut" );
156
159
160 ATH_CHECK( m_decPassNtrkKey.initialize() );
161 ATH_CHECK( m_decCutNtrkKey.initialize() );
162
163 ATH_MSG_INFO( " " << m_decPassNtrkKey.key() << " : pass Ntrk cut" );
164 ATH_MSG_INFO( " " << m_decCutNtrkKey.key() << " : Ntrk cut" );
165
168
169 ATH_CHECK( m_decPassScoreKey.initialize() );
170 ATH_CHECK( m_decCutScoreKey.initialize() );
171
172 ATH_MSG_INFO( " " << m_decPassScoreKey.key() << " : pass Score cut" );
173 ATH_MSG_INFO( " " << m_decCutScoreKey.key() << " : Score cut" );
174
176 ATH_CHECK( m_decAcceptKey.initialize() );
177
178#ifndef XAOD_STANDALONE
190 }
191#endif
192
194 ATH_CHECK( m_decAcceptKey.initialize() );
195
196 return StatusCode::SUCCESS;
197
198}
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)
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidJetContentKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutMHighKey
Gaudi::Property< bool > m_suppressOutputDependence
std::string m_truthLabelName
std::string m_weightDecorationName
String for scale factor decoration names.
bool m_useMassCut
Flag to indicate if mass window cut is used.
std::string m_wkpt
Tagger information.
std::string m_weightConfigPath
Path to the SF configuration root file.
StatusCode getConfigReader()
Get configReader StatusCode.
std::string m_weightFileName
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutMLowKey
WriteDecorHandle keys for cut values.
std::string m_containerName
Configurable members.
std::string m_weightHistogramName
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassMassKey
void printCuts() const
Print configured cuts.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidEventContentKey
asg::AcceptInfo m_acceptInfo
Object that stores the results for a jet.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTaggedKey
WriteDecorHandle keys for tagging bools.
virtual StatusCode initialize() override
Initialize the tool.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutNtrkKey
std::unique_ptr< TF1 > m_funcScoreCut
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decAcceptKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassD2Key
WriteDecorHandle keys.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutD2Key
bool m_useScore
Flag to indicate if a ML score is used.
std::unique_ptr< TF1 > m_funcD2Cut
Functions that are configurable for specific cut values.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassNtrkKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decPassScoreKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decCutScoreKey
bool m_useNtrk
Flag to indicate if Ntrk is used.
std::unique_ptr< TF1 > m_funcNtrkCut

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

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg_level_name()

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

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

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

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

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

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

Definition at line 101 of file AsgTool.cxx.

101 {
102
103 return MSG::name( msg().level() );
104 }
MsgStream & msg() const
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition MsgLevel.cxx:19

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

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

◆ outputHandles()

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

Return this algorithm's output handles.

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

◆ passKinRange()

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

Check if jet passes kinematic constraints.

Definition at line 353 of file JSSTaggerBase.cxx.

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

◆ print()

◆ printCuts()

void JSSTaggerBase::printCuts ( ) const
protectedinherited

Print configured cuts.

Definition at line 927 of file JSSTaggerBase.cxx.

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

◆ 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();
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 {
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 335 of file JSSTaggerBase.cxx.

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

◆ sysInitialize()

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

Perform system initialization for an algorithm.

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

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

◆ sysStart()

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

Handle START transition.

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

◆ tag()

StatusCode SmoothedWZTagger::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 D2 value

Get Score value

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

Evaluate the values of the lower score cut

Decorate the cut values

Create WriteDecorHandles

Decorate values

Evaluate the cut criteria on mass and d2

Check if it's a smooth three-variable tagger (ntrk)

Decorate Ntrk cut value

Create WriteDecorHandles

Decorate values

check if it is a 1D score tagger

Decorate jet with tagging summary

Get enum to decorate acceptData state if only using 2-var tagger

Pass mass cut

Fail mass cut

Get SF weight

Create WriteDecorHandles

Decorate values

Implements JSSTaggerBase.

Definition at line 200 of file SmoothedWZTagger.cxx.

200 {
201
202 ATH_MSG_DEBUG( "Obtaining Smooth WZ result" );
203
205 asg::AcceptData acceptData( &m_acceptInfo );
206
208 ATH_CHECK( resetCuts( acceptData ) );
209
211 ATH_CHECK( checkKinRange( jet, acceptData ) );
212
215 float jet_pt = jet.pt()/1000.0;
216 float jet_mass = jet.m()/1000.0;
217
220
222 static const SG::AuxElement::ConstAccessor<float> D2("D2");
223 float jet_d2 = D2(jet);
224
226 static const SG::AuxElement::ConstAccessor<float> Score(m_decorationName + "_ConstScore");
227 float jet_score = m_useScore ? Score(jet) : -99;
228
230 float cut_mass_low (-99), cut_mass_high (-99), cut_d2 (-99);
231 if(!m_useScore){
232 cut_mass_low = m_funcMassCutLow ->Eval(jet_pt);
233 cut_mass_high = m_funcMassCutHigh->Eval(jet_pt);
234 cut_d2 = m_funcD2Cut ->Eval(jet_pt);
235 }
236
238 float cut_score = m_useScore ? m_funcScoreCut -> Eval(jet_pt) : -99;
239
241
243 SG::WriteDecorHandle<xAOD::JetContainer, bool> decPassMass(m_decPassMassKey);
244 SG::WriteDecorHandle<xAOD::JetContainer, bool> decPassD2(m_decPassD2Key);
245 SG::WriteDecorHandle<xAOD::JetContainer, bool> decPassScore(m_decPassScoreKey);
246 SG::WriteDecorHandle<xAOD::JetContainer, bool> decTagged(m_decTaggedKey);
247 SG::WriteDecorHandle<xAOD::JetContainer, float> decCutMLow(m_decCutMLowKey);
248 SG::WriteDecorHandle<xAOD::JetContainer, float> decCutMHigh(m_decCutMHighKey);
249 SG::WriteDecorHandle<xAOD::JetContainer, float> decCutD2(m_decCutD2Key);
250 SG::WriteDecorHandle<xAOD::JetContainer, float> decCutScore(m_decCutScoreKey);
251
253 decCutMLow(jet) = cut_mass_low;
254 decCutMHigh(jet) = cut_mass_high;
255 decCutD2(jet) = cut_d2;
256 decCutScore(jet) = cut_score;
257
259 ATH_MSG_VERBOSE( "Cut Values : MassWindow = [" << cut_mass_low << "," << cut_mass_high << "], D2Cut = " << cut_d2 );
260 ATH_MSG_VERBOSE( "Cut Values : JetMass = " << jet_mass << ", D2 = " << jet_d2 );
261 ATH_MSG_VERBOSE( "Cut Values : score = " << jet_score);
262
263 if ( jet_mass >= cut_mass_low ) acceptData.setCutResult( "PassMassLow", true );
264
265 if ( jet_mass <= cut_mass_high ) acceptData.setCutResult( "PassMassHigh", true );
266
267 if ( jet_d2 < cut_d2 ) acceptData.setCutResult( "PassD2", true );
268
269 if ( jet_score > cut_score ) acceptData.setCutResult( "PassScore", true );
270
271 decPassMass(jet) = acceptData.getCutResult( "PassMassLow" ) && acceptData.getCutResult( "PassMassHigh" );
272 decPassD2(jet) = acceptData.getCutResult( "PassD2" );
273 decPassScore(jet) = acceptData.getCutResult( "PassScore" );
274
275 bool passCuts = acceptData.getCutResult( "PassMassLow" ) && acceptData.getCutResult( "PassMassHigh" );
276 passCuts = passCuts && acceptData.getCutResult( "PassD2" );
277
279 if ( m_useNtrk ) {
280
281 float cut_ntrk = m_funcNtrkCut->Eval(jet_pt);
282
284
286 SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidJetContent(m_decValidJetContentKey);
287 SG::WriteDecorHandle<xAOD::JetContainer, bool> decValidEventContent(m_decValidEventContentKey);
288 SG::WriteDecorHandle<xAOD::JetContainer, bool> decPassNtrk(m_decPassNtrkKey);
289 SG::WriteDecorHandle<xAOD::JetContainer, float> decCutNtrk(m_decCutNtrkKey);
290
292 decCutNtrk(jet) = cut_ntrk;
293
294 int pv_location = findPV();
295
296 if(pv_location != -1){
297 int jet_ntrk = GetUnGroomTracks(jet, pv_location);
298 if(jet_ntrk>=0){
299 if ( jet_ntrk < cut_ntrk ) acceptData.setCutResult( "PassNtrk", true );
300 decPassNtrk(jet) = acceptData.getCutResult( "PassNtrk" );
301 passCuts = passCuts && acceptData.getCutResult( "PassNtrk" );
302 }
303 else{
304 acceptData.setCutResult( "ValidJetContent", false );
305 decValidJetContent(jet) = false;
306 return StatusCode::FAILURE;
307 }
308 }
309 else {
310 acceptData.setCutResult( "ValidEventContent", false );
311 }
312
313 decValidJetContent(jet) = acceptData.getCutResult( "ValidJetContent" );
314 decValidEventContent(jet) = acceptData.getCutResult( "ValidEventContent" );
315
316 }
317
319 if(m_useScore){
320 passCuts = acceptData.getCutResult( "PassScore" );
321 }
322
324 decTagged(jet) = passCuts;
325
327 TagResult::TypeEnum myCutResultForSF = TagResult::UNKNOWN;
328 if ( !m_useNtrk && !m_useScore ) {
330 if ( acceptData.getCutResult("PassMassLow") && acceptData.getCutResult("PassMassHigh") ) {
331 if ( acceptData.getCutResult("PassD2") ) {
332 myCutResultForSF = TagResult::passMpassD2_2Var;
333 }
334 else {
335 myCutResultForSF = TagResult::passMfailD2_2Var;
336 }
337 }
339 else {
340 if ( acceptData.getCutResult("PassD2") ) {
341 myCutResultForSF = TagResult::failMpassD2_2Var;
342 }
343 else {
344 myCutResultForSF = TagResult::failMfailD2_2Var;
345 }
346 }
347 }
348
350 ATH_CHECK( getWeight( jet, (bool)acceptData, acceptData ) );
351
352 if ( m_calcSF ) {
353
355 SG::WriteDecorHandle<xAOD::JetContainer, float> decAccept(m_decAcceptKey);
356
358 decAccept(jet) = myCutResultForSF;
359
360 }
361
362 return StatusCode::SUCCESS;
363
364}
std::unique_ptr< TF1 > m_funcMassCutLow
TF1 for cut functions.
int calculateJSSRatios(const xAOD::Jet &jet) const
Calculate JSS moment ratios in case they are not already saved TODO: Remove this once JSSMomentTools ...
StatusCode checkKinRange(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
Check and record if jet passes kinematic constraints.
int GetUnGroomTracks(const xAOD::Jet &jet, int indexPV) const
Retrieve Ntrk variable from the ungroomed parent jet.
StatusCode resetCuts(asg::AcceptData &acceptData) const
Reset cuts.
int findPV() const
Find the PV (to be used for Ntrk)
StatusCode getWeight(const xAOD::Jet &jet, bool passSel, asg::AcceptData &acceptData) const
Get SF weight.
std::unique_ptr< TF1 > m_funcMassCutHigh

◆ 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) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }
std::vector< SG::VarHandleKeyArray * > m_vhka

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_calcSF

bool JSSTaggerBase::m_calcSF {}
protectedinherited

Flag to calculate scale factor.

Definition at line 192 of file JSSTaggerBase.h.

192{};

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

SG::WriteDecorHandleKey<xAOD::JetContainer> SmoothedWZTagger::m_decAcceptKey {this, "acceptName", "accept", "SG key for accept"}
private

Definition at line 51 of file SmoothedWZTagger.h.

51{this, "acceptName", "accept", "SG key for accept"};

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

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

◆ m_decCutD2Key

SG::WriteDecorHandleKey<xAOD::JetContainer> SmoothedWZTagger::m_decCutD2Key {this, "CutD2Name", "Cut_D2", "SG key for Cut_D2"}
private

Definition at line 48 of file SmoothedWZTagger.h.

48{this, "CutD2Name", "Cut_D2", "SG key for Cut_D2"};

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

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

◆ m_decCutMLowKey

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

WriteDecorHandle keys for cut values.

Definition at line 218 of file JSSTaggerBase.h.

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

◆ m_decCutNtrkKey

SG::WriteDecorHandleKey<xAOD::JetContainer> SmoothedWZTagger::m_decCutNtrkKey {this, "CutNtrkName", "Cut_Ntrk", "SG key for Cut_Ntrk"}
private

Definition at line 49 of file SmoothedWZTagger.h.

49{this, "CutNtrkName", "Cut_Ntrk", "SG key for Cut_Ntrk"};

◆ m_decCutScoreKey

SG::WriteDecorHandleKey<xAOD::JetContainer> SmoothedWZTagger::m_decCutScoreKey {this, "CutScoreName", "Cut_Score", "SG key for Cut_Score"}
private

Definition at line 50 of file SmoothedWZTagger.h.

50{this, "CutScoreName", "Cut_Score", "SG key for Cut_Score"};

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

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

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

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

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

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

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

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

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

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

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

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

◆ m_decorationName

std::string JSSTaggerBase::m_decorationName
protectedinherited

Decoration name.

Definition at line 189 of file JSSTaggerBase.h.

◆ m_decPassD2Key

SG::WriteDecorHandleKey<xAOD::JetContainer> SmoothedWZTagger::m_decPassD2Key {this, "PassD2Name", "PassD2", "SG key for PassD2"}
private

WriteDecorHandle keys.

Definition at line 44 of file SmoothedWZTagger.h.

44{this, "PassD2Name", "PassD2", "SG key for PassD2"};

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

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

◆ m_decPassNtrkKey

SG::WriteDecorHandleKey<xAOD::JetContainer> SmoothedWZTagger::m_decPassNtrkKey {this, "PassNtrkName", "PassNtrk", "SG key for PassNtrk"}
private

Definition at line 45 of file SmoothedWZTagger.h.

45{this, "PassNtrkName", "PassNtrk", "SG key for PassNtrk"};

◆ m_decPassScoreKey

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

Definition at line 46 of file SmoothedWZTagger.h.

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

◆ m_decScoreCutKey

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

Definition at line 220 of file JSSTaggerBase.h.

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

◆ m_decScoreValueKey

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

Definition at line 221 of file JSSTaggerBase.h.

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

◆ m_decSigeffSFKey

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

Definition at line 227 of file JSSTaggerBase.h.

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

◆ m_decTaggedKey

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

WriteDecorHandle keys for tagging bools.

Definition at line 69 of file JSSTaggerBase.h.

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

◆ m_decTau21WTAKey

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

WriteDecorHandle keys for JSS moments.

Definition at line 141 of file JSSTaggerBase.h.

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

◆ m_decTau32WTAKey

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

Definition at line 142 of file JSSTaggerBase.h.

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

◆ m_decTau42WTAKey

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

Definition at line 143 of file JSSTaggerBase.h.

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

◆ m_decValidEtaRangeKey

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

Definition at line 72 of file JSSTaggerBase.h.

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

◆ m_decValidEventContentKey

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

Definition at line 75 of file JSSTaggerBase.h.

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

◆ m_decValidJetContentKey

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

Definition at line 74 of file JSSTaggerBase.h.

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

◆ m_decValidKinRangeKey

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

Definition at line 73 of file JSSTaggerBase.h.

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

◆ m_decValidPtRangeHighKey

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

Definition at line 70 of file JSSTaggerBase.h.

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

◆ m_decValidPtRangeLowKey

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

Definition at line 71 of file JSSTaggerBase.h.

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

◆ m_decWeightKey

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

WriteDecorHandle keys for SF.

Definition at line 224 of file JSSTaggerBase.h.

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

◆ m_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_efficiencyHistogramName

std::string JSSTaggerBase::m_efficiencyHistogramName
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_funcD2Cut

std::unique_ptr<TF1> SmoothedWZTagger::m_funcD2Cut
private

Functions that are configurable for specific cut values.

Definition at line 39 of file SmoothedWZTagger.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_funcNtrkCut

std::unique_ptr<TF1> SmoothedWZTagger::m_funcNtrkCut
private

Definition at line 40 of file SmoothedWZTagger.h.

◆ m_funcScoreCut

std::unique_ptr<TF1> SmoothedWZTagger::m_funcScoreCut
private

Definition at line 41 of file SmoothedWZTagger.h.

◆ m_isMC

bool JSSTaggerBase::m_isMC {}
protectedinherited

Definition at line 193 of file JSSTaggerBase.h.

193{};

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

133{};

◆ m_jetPtMax

float JSSTaggerBase::m_jetPtMax {}
protectedinherited

Definition at line 132 of file JSSTaggerBase.h.

132{};

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

131{};

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

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

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

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

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

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

◆ m_readECFG212Key

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

Definition at line 171 of file JSSTaggerBase.h.

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

◆ m_readECFG311Key

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

Definition at line 170 of file JSSTaggerBase.h.

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

◆ m_readECFG331Key

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

Definition at line 169 of file JSSTaggerBase.h.

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

◆ m_readParentKey

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

Definition at line 173 of file JSSTaggerBase.h.

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

◆ m_readQwKey

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

Definition at line 164 of file JSSTaggerBase.h.

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

◆ m_readSphericityKey

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

Definition at line 167 of file JSSTaggerBase.h.

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

◆ m_readSplit12Key

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

Definition at line 161 of file JSSTaggerBase.h.

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

◆ m_readSplit23Key

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

Definition at line 162 of file JSSTaggerBase.h.

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

◆ m_readTau1WTAKey

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

ReadDecorHandle keys for JSS moments.

Definition at line 152 of file JSSTaggerBase.h.

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

◆ m_readTau2WTAKey

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

Definition at line 153 of file JSSTaggerBase.h.

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

◆ m_readTau3WTAKey

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

Definition at line 154 of file JSSTaggerBase.h.

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

◆ m_readTau4WTAKey

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

Definition at line 155 of file JSSTaggerBase.h.

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

◆ m_readThrustMajKey

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

Definition at line 166 of file JSSTaggerBase.h.

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

◆ m_readTruthLabelKey

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

Truth label ReadDecorHandle key.

Definition at line 215 of file JSSTaggerBase.h.

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

◆ m_strD2Cut

std::string SmoothedWZTagger::m_strD2Cut
private

Store functional form of cuts.

Definition at line 34 of file SmoothedWZTagger.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_strNtrkCut

std::string SmoothedWZTagger::m_strNtrkCut
private

Definition at line 35 of file SmoothedWZTagger.h.

◆ m_strScoreCut

std::string SmoothedWZTagger::m_strScoreCut
private

Definition at line 36 of file SmoothedWZTagger.h.

◆ m_suppressInputDependence

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

Definition at line 175 of file JSSTaggerBase.h.

175{this, "SuppressInputDependence", false};

◆ m_suppressOutputDependence

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

Definition at line 176 of file JSSTaggerBase.h.

176{this, "SuppressOutputDependence", false};

◆ m_tagClass

TAGCLASS JSSTaggerBase::m_tagClass {Unknown}
protectedinherited

Definition at line 98 of file JSSTaggerBase.h.

98{Unknown};

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

208{};

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

bool SmoothedWZTagger::m_useNtrk {}
private

Flag to indicate if Ntrk is used.

Definition at line 28 of file SmoothedWZTagger.h.

28{};

◆ m_useScore

bool SmoothedWZTagger::m_useScore {}
private

Flag to indicate if a ML score is used.

Definition at line 31 of file SmoothedWZTagger.h.

31{};

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