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

#include <tauMonitorAlgorithm.h>

Inheritance diagram for tauMonitorAlgorithm:
Collaboration diagram for tauMonitorAlgorithm:

Public Types

enum class  Environment_t {
  user = 0 , online , tier0 , tier0Raw ,
  tier0ESD , AOD , altprod
}
 Specifies the processing environment. More...
enum class  DataType_t {
  userDefined = 0 , monteCarlo , collisions , cosmics ,
  heavyIonCollisions
}
 Specifies what type of input data is being monitored. More...

Public Member Functions

 tauMonitorAlgorithm (const std::string &name, ISvcLocator *pSvcLocator)
virtual ~tauMonitorAlgorithm ()
virtual StatusCode initialize () override
 initialize
virtual StatusCode fillHistograms (const EventContext &ctx) const override
 adds event to the monitoring histograms
virtual StatusCode execute (const EventContext &ctx) const override
 Applies filters and trigger requirements.
void fill (const ToolHandle< GenericMonitoringTool > &groupHandle, std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &&variables) const
 Fills a vector of variables to a group by reference.
void fill (const ToolHandle< GenericMonitoringTool > &groupHandle, const std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &variables) const
 Fills a vector of variables to a group by reference.
template<typename... T>
void fill (const ToolHandle< GenericMonitoringTool > &groupHandle, T &&... variables) const
 Fills a variadic list of variables to a group by reference.
void fill (const std::string &groupName, std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &&variables) const
 Fills a vector of variables to a group by name.
void fill (const std::string &groupName, const std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &variables) const
 Fills a vector of variables to a group by name.
template<typename... T>
void fill (const std::string &groupName, T &&... variables) const
 Fills a variadic list of variables to a group by name.
Environment_t environment () const
 Accessor functions for the environment.
Environment_t envStringToEnum (const std::string &str) const
 Convert the environment string from the python configuration to an enum object.
DataType_t dataType () const
 Accessor functions for the data type.
DataType_t dataTypeStringToEnum (const std::string &str) const
 Convert the data type string from the python configuration to an enum object.
const ToolHandle< GenericMonitoringTool > & getGroup (const std::string &name) const
 Get a specific monitoring tool from the tool handle array.
const ToolHandle< Trig::TrigDecisionTool > & getTrigDecisionTool () const
 Get the trigger decision tool member.
bool trigChainsArePassed (const std::vector< std::string > &vTrigNames) const
 Check whether triggers are passed.
SG::ReadHandle< xAOD::EventInfoGetEventInfo (const EventContext &) const
 Return a ReadHandle for an EventInfo object (get run/event numbers, etc.)
virtual float lbAverageInteractionsPerCrossing (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate the average mu, i.e.
virtual float lbInteractionsPerCrossing (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate instantaneous number of interactions, i.e.
virtual float lbAverageLuminosity (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate average luminosity (in ub-1 s-1 => 10^30 cm-2 s-1).
virtual float lbLuminosityPerBCID (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate the instantaneous luminosity per bunch crossing.
virtual double lbDuration (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate the duration of the luminosity block (in seconds)
virtual float lbAverageLivefraction (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate the average luminosity livefraction.
virtual float livefractionPerBCID (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate the live fraction per bunch crossing ID.
virtual double lbLumiWeight (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate the average integrated luminosity multiplied by the live fraction.
virtual StatusCode parseList (const std::string &line, std::vector< std::string > &result) const
 Parse a string into a vector.
virtual StatusCode sysInitialize () override
 Override sysInitialize.
virtual bool isClonable () const override
 Specify if the algorithm is clonable.
virtual unsigned int cardinality () const override
 Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
virtual bool filterPassed (const EventContext &ctx) const
virtual void setFilterPassed (bool state, const EventContext &ctx) const
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 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

Public Attributes

 flags = initConfigFlags()
str nightly = '/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/CommonInputs/'
str file = 'data16_13TeV.00311321.physics_Main.recon.AOD.r9264/AOD.11038520._000001.pool.root.1'
 Files
 isMC
 HISTFileName
 cfg = MainServicesCfg(flags)
 exampleMonitorAcc = tauMonitoringConfig(flags)
 OutputLevel
 withDetails

Protected Member Functions

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

ToolHandleArray< GenericMonitoringToolm_tools {this,"GMTools",{}}
 Array of Generic Monitoring Tools.
PublicToolHandle< Trig::TrigDecisionToolm_trigDecTool
 Tool to tell whether a specific trigger is passed.
ToolHandleArray< IDQFilterToolm_DQFilterTools {this,"FilterTools",{}}
 Array of Data Quality filter tools.
SG::ReadCondHandleKey< LuminosityCondDatam_lumiDataKey {this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"}
SG::ReadCondHandleKey< LBDurationCondDatam_lbDurationDataKey {this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"}
SG::ReadCondHandleKey< TrigLiveFractionCondDatam_trigLiveFractionDataKey {this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"}
AthMonitorAlgorithm::Environment_t m_environment
 Instance of the Environment_t enum.
AthMonitorAlgorithm::DataType_t m_dataType
 Instance of the DataType_t enum.
Gaudi::Property< std::string > m_environmentStr {this,"Environment","user"}
 Environment string pulled from the job option and converted to enum.
Gaudi::Property< std::string > m_dataTypeStr {this,"DataType","userDefined"}
 DataType string pulled from the job option and converted to enum.
Gaudi::Property< std::string > m_triggerChainString {this,"TriggerChain",""}
 Trigger chain string pulled from the job option and parsed into a vector.
std::vector< std::string > m_vTrigChainNames
 Vector of trigger chain names parsed from trigger chain string.
Gaudi::Property< std::string > m_fileKey {this,"FileKey",""}
 Internal Athena name for file.
Gaudi::Property< bool > m_useLumi {this,"EnableLumi",false}
 Allows use of various luminosity functions.
Gaudi::Property< float > m_defaultLBDuration {this,"DefaultLBDuration",60.}
 Default duration of one lumi block.
Gaudi::Property< int > m_detailLevel {this,"DetailLevel",0}
 Sets the level of detail used in the monitoring.
SG::ReadHandleKey< xAOD::EventInfom_EventInfoKey {this,"EventInfoKey","EventInfo"}
 Key for retrieving EventInfo from StoreGate.

Private Types

typedef std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > MonVarVec_t
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

std::vector< int > m_abGroups1
std::vector< std::vector< int > > m_abGroups2
std::map< std::string, int > m_cGroups1
std::map< std::string, std::map< std::string, int > > m_cGroups2
SG::ReadHandleKey< xAOD::TauJetContainerm_TauContainerKey {this, "TauRecContainer", "TauJets"}
SG::ReadDecorHandleKey< xAOD::TauJetContainerm_offlineGNTauScoreDecorKey {this, "OfflineGNTauScoreDecorKey", "GNTauScore_v0prune", "Offline GNTau Score decoration key"}
SG::ReadDecorHandleKey< xAOD::TauJetContainerm_offlineGNTauTransScoreDecorKey {this, "OfflineGNTauTransScoreDecorKey", "GNTauScoreSigTrans_v0prune", "Offline GNTau Trans Score decoration key"}
SG::ReadDecorHandleKey< xAOD::TauJetContainerm_offlineGNTauLooseWPDecorKey {this, "OfflineGNTauLooseWPDecorKey", "GNTauL_v0prune", "Offline GNTau Loose WP decoration key"}
SG::ReadDecorHandleKey< xAOD::TauJetContainerm_offlineGNTauMediumWPDecorKey {this, "OfflineGNTauMediumWPDecorKey", "GNTauM_v0prune", "Offline GNTau Medium WP decoration key"}
SG::ReadDecorHandleKey< xAOD::TauJetContainerm_offlineGNTauTightWPDecorKey {this, "OfflineGNTauTightWPDecorKey", "GNTauT_v0prune", "Offline GNTau Tight WP decoration key"}
Gaudi::Property< float > m_etaMin {this, "etaMin", -1.}
Gaudi::Property< float > m_etaMax {this, "etaMax", 3.0}
Gaudi::Property< std::string > m_kinGroupName {this, "kinGroupName", "tauMonKinGroupBA"}
std::string m_name
std::unordered_map< std::string, size_t > m_toolLookupMap
const ToolHandle< GenericMonitoringToolm_dummy
Gaudi::Property< bool > m_enforceExpressTriggers
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
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 18 of file tauMonitorAlgorithm.h.

Member Typedef Documentation

◆ MonVarVec_t

typedef std::vector<std::reference_wrapper<Monitored::IMonitoredVariable> > AthMonitorAlgorithm::MonVarVec_t
privateinherited

Definition at line 370 of file AthMonitorAlgorithm.h.

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ DataType_t

enum class AthMonitorAlgorithm::DataType_t
stronginherited

Specifies what type of input data is being monitored.

An enumeration of the different types of data the monitoring application may be running over. This can be used to select which histograms to produce, e.g. to prevent the production of colliding-beam histograms when running on cosmic-ray data. Strings of the same names may be given as jobOptions.

Enumerator
userDefined 
monteCarlo 
collisions 
cosmics 
heavyIonCollisions 

Definition at line 194 of file AthMonitorAlgorithm.h.

194 {
195 userDefined = 0,
196 monteCarlo,
197 collisions,
198 cosmics,
199 heavyIonCollisions,
200 };
AthConfigFlags cosmics(AthConfigFlags flags, str instanceName, str recoMode)

◆ Environment_t

enum class AthMonitorAlgorithm::Environment_t
stronginherited

Specifies the processing environment.

The running environment may be used to select which histograms are produced, and where they are located in the output. For example, the output paths of the histograms are different for the "user", "online" and the various offline flags. Strings of the same names may be given as jobOptions.

Enumerator
user 
online 
tier0 
tier0Raw 
tier0ESD 
AOD 
altprod 

Definition at line 175 of file AthMonitorAlgorithm.h.

175 {
176 user = 0,
177 online,
178 tier0,
179 tier0Raw,
180 tier0ESD,
181 AOD,
182 altprod,
183 };

Constructor & Destructor Documentation

◆ tauMonitorAlgorithm()

tauMonitorAlgorithm::tauMonitorAlgorithm ( const std::string & name,
ISvcLocator * pSvcLocator )

Definition at line 15 of file tauMonitorAlgorithm.cxx.

17 : AthMonitorAlgorithm(name, pSvcLocator)
18{}
Base class for Athena Monitoring Algorithms.

◆ ~tauMonitorAlgorithm()

tauMonitorAlgorithm::~tauMonitorAlgorithm ( )
virtual

Definition at line 20 of file tauMonitorAlgorithm.cxx.

20{}

Member Function Documentation

◆ cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality ( ) const
overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

64{
65 return 0;
66}

◆ dataType()

DataType_t AthMonitorAlgorithm::dataType ( ) const
inlineinherited

Accessor functions for the data type.

Returns
the current value of the class's DataType_t instance.

Definition at line 224 of file AthMonitorAlgorithm.h.

224{ return m_dataType; }
AthMonitorAlgorithm::DataType_t m_dataType
Instance of the DataType_t enum.

◆ dataTypeStringToEnum()

AthMonitorAlgorithm::DataType_t AthMonitorAlgorithm::dataTypeStringToEnum ( const std::string & str) const
inherited

Convert the data type string from the python configuration to an enum object.

Returns
a value in the DataType_t enumeration which matches the input string.

Definition at line 144 of file AthMonitorAlgorithm.cxx.

144 {
145 // convert the string to all lowercase
146 std::string lowerCaseStr = str;
147 std::transform(lowerCaseStr.begin(), lowerCaseStr.end(), lowerCaseStr.begin(), ::tolower);
148
149 // check if it matches one of the enum choices
150 if( lowerCaseStr == "userdefined" ) {
152 } else if( lowerCaseStr == "montecarlo" ) {
154 } else if( lowerCaseStr == "collisions" ) {
156 } else if( lowerCaseStr == "cosmics" ) {
157 return DataType_t::cosmics;
158 } else if( lowerCaseStr == "heavyioncollisions" ) {
160 } else { // otherwise, warn the user and return "userDefined"
161 ATH_MSG_WARNING("AthMonitorAlgorithm::dataTypeStringToEnum(): Unknown data type "
162 <<str<<", returning userDefined.");
164 }
165}
#define ATH_MSG_WARNING(x)
void tolower(std::string &s)

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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>

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

◆ environment()

Environment_t AthMonitorAlgorithm::environment ( ) const
inlineinherited

Accessor functions for the environment.

Returns
the current value of the class's Environment_t instance.

Definition at line 208 of file AthMonitorAlgorithm.h.

208{ return m_environment; }
AthMonitorAlgorithm::Environment_t m_environment
Instance of the Environment_t enum.

◆ envStringToEnum()

AthMonitorAlgorithm::Environment_t AthMonitorAlgorithm::envStringToEnum ( const std::string & str) const
inherited

Convert the environment string from the python configuration to an enum object.

Returns
a value in the Environment_t enumeration which matches the input string.

Definition at line 116 of file AthMonitorAlgorithm.cxx.

116 {
117 // convert the string to all lowercase
118 std::string lowerCaseStr = str;
119 std::transform(lowerCaseStr.begin(), lowerCaseStr.end(), lowerCaseStr.begin(), ::tolower);
120
121 // check if it matches one of the enum choices
122 if( lowerCaseStr == "user" ) {
123 return Environment_t::user;
124 } else if( lowerCaseStr == "online" ) {
126 } else if( lowerCaseStr == "tier0" ) {
128 } else if( lowerCaseStr == "tier0raw" ) {
130 } else if( lowerCaseStr == "tier0esd" ) {
132 } else if( lowerCaseStr == "aod" ) {
133 return Environment_t::AOD;
134 } else if( lowerCaseStr == "altprod" ) {
136 } else { // otherwise, warn the user and return "user"
137 ATH_MSG_WARNING("AthMonitorAlgorithm::envStringToEnum(): Unknown environment "
138 <<str<<", returning user.");
139 return Environment_t::user;
140 }
141}

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

StatusCode AthMonitorAlgorithm::execute ( const EventContext & ctx) const
overridevirtualinherited

Applies filters and trigger requirements.

Then, calls fillHistograms().

Parameters
ctxevent context for reentrant Athena call
Returns
StatusCode

Definition at line 77 of file AthMonitorAlgorithm.cxx.

77 {
78
79 // Checks that all of the DQ filters are passed. If any one of the filters
80 // fails, return SUCCESS code and do not fill the histograms with the event.
81 for ( const auto& filterItr : m_DQFilterTools ) {
82 if (!filterItr->accept()) {
83 ATH_MSG_DEBUG("Event rejected due to filter tool.");
84 return StatusCode::SUCCESS;
85 }
86 }
87
88 // Trigger: If there is a decision tool and the chains fail, skip the event.
90 ATH_MSG_DEBUG("Event rejected due to trigger filter.");
91 return StatusCode::SUCCESS;
92 }
93
94 ATH_MSG_DEBUG("Event accepted!");
95 return fillHistograms(ctx);
96}
#define ATH_MSG_DEBUG(x)
virtual StatusCode fillHistograms(const EventContext &ctx) const =0
adds event to the monitoring histograms
bool trigChainsArePassed(const std::vector< std::string > &vTrigNames) const
Check whether triggers are passed.
std::vector< std::string > m_vTrigChainNames
Vector of trigger chain names parsed from trigger chain string.
PublicToolHandle< Trig::TrigDecisionTool > m_trigDecTool
Tool to tell whether a specific trigger is passed.
ToolHandleArray< IDQFilterTool > m_DQFilterTools
Array of Data Quality filter tools.

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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

◆ extraOutputDeps()

const DataObjIDColl & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

90{
91 // If we didn't find any symlinks to add, just return the collection
92 // from the base class. Otherwise, return the extended collection.
93 if (!m_extendedExtraObjects.empty()) {
95 }
97}
An algorithm that can be simultaneously executed in multiple threads.

◆ fillHistograms()

StatusCode tauMonitorAlgorithm::fillHistograms ( const EventContext & ctx) const
overridevirtual

adds event to the monitoring histograms

User will overwrite this function. Histogram booking is no longer done in C++. This function is called in execute once the filters are all passed.

Parameters
ctxforwarded from execute
Returns
StatusCode

Implements AthMonitorAlgorithm.

Definition at line 44 of file tauMonitorAlgorithm.cxx.

44 {
45
46 SG::ReadHandle<xAOD::TauJetContainer> taus(m_TauContainerKey, ctx);
47
48 if (!taus.isValid()) {
49 ATH_MSG_ERROR("evtStore() does not contain tau Collection with name "
50 << m_TauContainerKey);
51 return StatusCode::FAILURE;
52 }
53
54 auto shallowCopy = xAOD::shallowCopyContainer(*taus);
55 std::unique_ptr<xAOD::TauJetContainer> shallowTaus(shallowCopy.first);
56 std::unique_ptr<xAOD::ShallowAuxContainer> shallowTausAux(shallowCopy.second);
57
58 const int lowerEtThreshold = 15;
59 const int higherEtThreshold = 75;
60 auto tool = getGroup(m_kinGroupName);
61
62 auto tauEta = Monitored::Scalar<float>("tauEta", 0.0);
63
64 auto tauPhi = Monitored::Scalar<float>("tauPhi", 0.0);
65 auto tauEt = Monitored::Scalar<float>("tauEt", 0.0);
66
67 std::string tauid = "";
68 if(m_TauContainerKey.key().find("TMTauJets") != std::string::npos)
69 {
70 tauid = "GNTau";
71 } else {
72 tauid = "RNN";
73 }
74
75 auto tauEtEt15Loose = Monitored::Scalar<float>("tauEtEt15"+tauid+"Loose", 0.0);
76 auto panModeEt15Loose =
77 Monitored::Scalar<float>("panModeEt15"+tauid+"Loose", 0.0);
78 auto panModeSubstructure =
79 Monitored::Scalar<float>("panModeSubstructure", 0.0);
80
81 auto tauPhiEt15 = Monitored::Scalar<float>("tauPhiEt15", 0.0);
82 auto tauEtaEt15 = Monitored::Scalar<float>("tauEtaEt15", 0.0);
83
84 auto tauPhiEt15Loose = Monitored::Scalar<float>("tauPhiEt15"+tauid+"Loose", 0.0);
85 auto tauEtaEt15Loose = Monitored::Scalar<float>("tauEtaEt15"+tauid+"Loose", 0.0);
86
87 auto tauCharge = Monitored::Scalar<int>("tauCharge", 0.0);
88 auto JetScore = Monitored::Scalar<float>(tauid+"JetScore", 0.0);
89 auto JetScoreSigTrans =
90 Monitored::Scalar<float>(tauid+"JetScoreSigTrans", 0.0);
91
92 auto RNNEleScore = Monitored::Scalar<float>("RNNEleScore", 0.0);
94 Monitored::Scalar<float>("RNNEleScoreSigTrans", 0.0);
95
96 auto NumTracks = Monitored::Scalar<int>("NumTracks", 0.0);
97 auto NumTracksEt15Loose =
98 Monitored::Scalar<int>("NumTracksEt15"+tauid+"Loose", 0.0);
99
100 auto nTauCandidates = Monitored::Scalar<int>("nTauCandidates", 0.0);
101 auto nHighPtTauCandidates =
102 Monitored::Scalar<int>("nHighPtTauCandidates", 0.0);
103 auto nClusters = Monitored::Scalar<int>("nClusters", 0.0);
104 auto nClustersEt15Loose =
105 Monitored::Scalar<int>("nClustersEt15"+tauid+"Loose", 0.0);
106
107 auto tauEtLoose = Monitored::Scalar<float>("tauEt"+tauid+"Loose", 0.0);
108 auto tauEtaLoose = Monitored::Scalar<float>("tauEta"+tauid+"Loose", 0.0);
109 auto tauPhiLoose = Monitored::Scalar<float>("tauPhi"+tauid+"Loose", 0.0);
110 auto NumTracksLoose = Monitored::Scalar<float>("NumTracks"+tauid+"Loose", 0.0);
111
112 auto tauEtMedium = Monitored::Scalar<float>("tauEt"+tauid+"Medium", 0.0);
113 auto tauEtaMedium = Monitored::Scalar<float>("tauEta"+tauid+"Medium", 0.0);
114 auto tauPhiMedium = Monitored::Scalar<float>("tauPhi"+tauid+"Medium", 0.0);
115 auto NumTracksMedium = Monitored::Scalar<float>("NumTracks"+tauid+"Medium", 0.0);
116
117 auto tauEtTight = Monitored::Scalar<float>("tauEt"+tauid+"Tight", 0.0);
118 auto tauEtaTight = Monitored::Scalar<float>("tauEta"+tauid+"Tight", 0.0);
119 auto tauPhiTight = Monitored::Scalar<float>("tauPhi"+tauid+"Tight", 0.0);
120 auto NumTracksTight = Monitored::Scalar<float>("NumTracks"+tauid+"Tight", 0.0);
121
122 auto LB = Monitored::Scalar<int>("LB", 0.0);
123
124 auto EMRadius = Monitored::Scalar<float>("EMRadius", 0.0);
125 auto hadRadius = Monitored::Scalar<float>("hadRadius", 0.0);
126 auto isolFrac = Monitored::Scalar<float>("isolFrac", 0.0);
127 auto etEMAtEMScale = Monitored::Scalar<float>("etEMAtEMScale", 0.0);
128
129 auto etHadAtEMScale = Monitored::Scalar<float>("etHadAtEMScale", 0.0);
130 auto centFrac = Monitored::Scalar<float>("centFrac", 0.0);
131 auto jetSeedEta = Monitored::Scalar<float>("jetSeedEta", 0.0);
132 auto jetSeedPhi = Monitored::Scalar<float>("jetSeedPhi", 0.0);
133 auto jetSeedPt = Monitored::Scalar<float>("jetSeedPt", 0.0);
134
135 auto muonVeto = Monitored::Scalar<float>("muonVeto", 0.0);
136
137 auto tauLoose = Monitored::Scalar<float>("tau"+tauid+"Loose", 0.0);
138 auto tauMedium = Monitored::Scalar<float>("tau"+tauid+"Medium", 0.0);
139 auto tauTight = Monitored::Scalar<float>("tau"+tauid+"Tight", 0.0);
140
141 auto PSSFrac = Monitored::Scalar<float>("PSSFrac", 0.0);
142 auto EMFrac = Monitored::Scalar<float>("EMFrac", 0.0);
143
144 auto EMFracTrk = Monitored::Scalar<float>("EMFracTrk", 0.0);
145 auto nNeutPFO = Monitored::Scalar<float>("nNeutPFO", 0.0);
146 auto nShot = Monitored::Scalar<float>("nShot", 0.0);
147
148 auto BDTScoreAsP0 = Monitored::Scalar<float>("BDTScoreAsP0", 0.0);
149 auto dRmax = Monitored::Scalar<float>("dRmax", 0.0);
150
151 auto ipSigLeadTrk = Monitored::Scalar<float>("ipSigLeadTrk", 0.0);
152 auto massTrkSys = Monitored::Scalar<float>("massTrkSys", 0.0);
153 auto etOverPtLeadTrack = Monitored::Scalar<float>("etOverPtLeadTrack", 0.0);
154 auto ptRatioEflowApprox = Monitored::Scalar<float>("ptRatioEflowApprox", 0.0);
155 auto trFlightPathSig = Monitored::Scalar<float>("trFlightPathSig", 0.0);
156 auto trkAvgDist = Monitored::Scalar<float>("trkAvgDist", 0.0);
157
158 auto panEta = Monitored::Scalar<float>("panEta", 0.0);
159 auto panPhi = Monitored::Scalar<float>("panPhi", 0.0);
160 auto panPt = Monitored::Scalar<float>("panPt", 0.0);
161 auto d0 = Monitored::Scalar<float>("d0", 0.0);
162 auto dRJetSeedAxis = Monitored::Scalar<float>("dRJetSeedAxis", 0.0);
163 auto z0 = Monitored::Scalar<float>("z0", 0.0);
164
165 auto etaTrack = Monitored::Scalar<float>("etaTrack", 0.0);
166 auto ptTrack = Monitored::Scalar<float>("ptTrack", 0.0);
167 auto phiTrack = Monitored::Scalar<float>("phiTrack", 0.0);
168 auto leadTrkPt = Monitored::Scalar<float>("leadTrkPt", 0.0);
169 auto nHighPtTaus = Monitored::Scalar<float>("nHighPtTaus", 0.0);
171 Monitored::Scalar<float>("numberOfTRTHighThresholdHits", 0.0);
173 Monitored::Scalar<float>("numberOfTRTHighThresholdOutliers", 0.0);
174 auto numberOfTRTHits = Monitored::Scalar<float>("numberOfTRTHits", 0.0);
176 Monitored::Scalar<float>("numberOfTRTOutliers", 0.0);
178 Monitored::Scalar<float>("ipZ0SinThetaSigLeadTrk", 0.0);
179 auto numberOfPixelHits = Monitored::Scalar<float>("numberOfPixelHits", 0.0);
181 Monitored::Scalar<float>("numberOfPixelSharedHits", 0.0);
182 auto numberOfSCTHits = Monitored::Scalar<float>("numberOfSCTHits", 0.0);
184 Monitored::Scalar<float>("numberOfSCTSharedHits", 0.0);
185 auto rConv = Monitored::Scalar<float>("rConv", 0.0);
186 auto rConvII = Monitored::Scalar<float>("rConvII", 0.0);
187
188 // trackLogSeedJetPt
189 auto trackLogSeedJetPt = Monitored::Scalar<float>("trackLogSeedJetPt", 0.0);
190 auto trackLogPt = Monitored::Scalar<float>("trackLogPt", 0.0);
191 auto trackEta = Monitored::Scalar<float>("trackEta", 0.0);
192 auto trackd0TJVA = Monitored::Scalar<float>("trackd0TJVA", 0.0);
193 auto trackZ0SinthetaTJVA =
194 Monitored::Scalar<float>("trackZ0SinthetaTJVA", 0.0);
195 auto trackD0SigTJVA = Monitored::Scalar<float>("trackD0SigTJVA", 0.0);
196 auto trackZ0sinthetaSigTJVA =
197 Monitored::Scalar<float>("trackZ0sinthetaSigTJVA", 0.0);
198 auto trackCharge = Monitored::Scalar<float>("trackCharge", 0.0);
199 auto trackqOverP = Monitored::Scalar<float>("trackqOverP", 0.0);
200 auto trackLogRConv = Monitored::Scalar<float>("trackLogRConv", 0.0);
201 auto trackTanhRConvII = Monitored::Scalar<float>("trackTanhRConvII", 0.0);
202 auto trackPtRatioSeedJet =
203 Monitored::Scalar<float>("trackPtRatioSeedJet", 0.0);
204 auto trackdRJetSeedAxis = Monitored::Scalar<float>("trackdRJetSeedAxis", 0.0);
205 auto trackNInnermostPixHits =
206 Monitored::Scalar<float>("trackNInnermostPixHits", 0.0);
207 auto trackNPixHits = Monitored::Scalar<float>("trackNPixHits", 0.0);
208 auto trackNSiHits = Monitored::Scalar<float>("trackNSiHits", 0.0);
209 auto trackeProbabilityHT =
210 Monitored::Scalar<float>("trackeProbabilityHT", 0.0);
211 auto trackeProbabilityNN =
212 Monitored::Scalar<float>("trackeProbabilityNN", 0.0);
213 auto trackeProbabilityHTorNN =
214 Monitored::Scalar<float>("trackeProbabilityHTorNN", 0.0);
215 auto trackIdScoreCharged = Monitored::Scalar<float>("track", 0.0);
216 auto trackIdScoreIso = Monitored::Scalar<float>("track", 0.0);
217 auto trackIdScoreConv = Monitored::Scalar<float>("track", 0.0);
218 auto trackIdScoreFake = Monitored::Scalar<float>("track", 0.0);
219
220 auto clusterLogEt = Monitored::Scalar<float>("clusterLogEt", 0.0);
221 auto clusterEta = Monitored::Scalar<float>("clusterEta", 0.0);
222 auto clusterPhi = Monitored::Scalar<float>("clusterPhi", 0.0);
223 auto clusterSecondR = Monitored::Scalar<float>("clusterSecondR", 0.0);
224 auto clusterSecondLambda =
225 Monitored::Scalar<float>("clusterSecondLambda", 0.0);
226 auto clusterCenterLambda =
227 Monitored::Scalar<float>("clusterCenterLambda", 0.0);
228 auto clusterFirstEngDens =
229 Monitored::Scalar<float>("clusterFirstEngDens", 0.0);
230 auto clusterEMproba = Monitored::Scalar<float>("clusterEMproba", 0.0);
231 auto clustersMeanCenterLambda =
232 Monitored::Scalar<float>("clustersMeanCenterLambda", 0.0);
233 auto clustersMeanFirstEngDens =
234 Monitored::Scalar<float>("clustersMeanFirstEngDens", 0.0);
235 auto clustersMeanEMProbability =
236 Monitored::Scalar<float>("clustersMeanEMProbability", 0.0);
237 auto clustersMeanSecondLambda =
238 Monitored::Scalar<float>("clustersMeanSecondLambda", 0.0);
239 auto clustersMeanPresamplerFrac =
240 Monitored::Scalar<float>("clustersMeanPresamplerFrac", 0.0);
241 auto clustersPFOEngRelDiff =
242 Monitored::Scalar<float>("clustersPFOEngRelDiff", 0.0);
243
244 nTauCandidates = 0;
245
246 static const SG::ConstAccessor<char> passThinningAcc("passThinning");
247
248 for (const auto tau : *shallowTaus) {
249
250 if (!passThinningAcc.withDefault(*tau, true)) continue;
251
252 tauEta = tau->eta();
253 tauPhi = tau->phi();
254 tauEt = tau->pt() / GeV;
255 tauCharge = tau->charge();
256 NumTracks = tau->nTracks();
257 nClusters = tau->detail<int>(xAOD::TauJetParameters::numTopoClusters);
258 LB = GetEventInfo(ctx)->lumiBlock();
259
260 // calo
261 EMRadius = tau->detail<float>(xAOD::TauJetParameters::EMRadius);
262 hadRadius = tau->detail<float>(xAOD::TauJetParameters::hadRadius);
263 isolFrac = tau->detail<float>(xAOD::TauJetParameters::isolFrac);
266
267 centFrac = tau->detail<float>(xAOD::TauJetParameters::centFrac);
268 jetSeedEta = tau->etaJetSeed();
269 jetSeedPhi = tau->phiJetSeed();
270 jetSeedPt = tau->ptJetSeed() / GeV;
271
275
276 muonVeto = tau->isTau(xAOD::TauJetParameters::MuonVeto);
277
278 // check to understand which TauID we can use
279 SG::ReadDecorHandle<xAOD::TauJetContainer, float> trans_score{m_offlineGNTauTransScoreDecorKey, ctx};
280 if( trans_score.isAvailable()) {
281
282 // we are using a special container decorated with GNTau
283 SG::ReadDecorHandle<xAOD::TauJetContainer, float> score{m_offlineGNTauScoreDecorKey, ctx};
284 JetScore = score.isAvailable() ? score(*tau) : -1234;
285
286 JetScoreSigTrans = trans_score(*tau);
287
288 SG::ReadDecorHandle<xAOD::TauJetContainer, char> tauid_loose{m_offlineGNTauLooseWPDecorKey, ctx};
289 tauLoose = tauid_loose.isAvailable() ? tauid_loose(*tau) : -1234;
290
291 SG::ReadDecorHandle<xAOD::TauJetContainer, char> tauid_medium{m_offlineGNTauMediumWPDecorKey, ctx};
292 tauMedium = tauid_medium.isAvailable() ? tauid_medium(*tau) : -1234;
293
294 SG::ReadDecorHandle<xAOD::TauJetContainer, char> tauid_tight{m_offlineGNTauTightWPDecorKey, ctx};
295 tauTight = tauid_tight.isAvailable() ? tauid_tight(*tau) : -1234;
296
297 } else{
298 // GNTau is not present -> use RNN based TauID
299 JetScore = tau->discriminant(xAOD::TauJetParameters::TauID::RNNJetScore);
300 JetScoreSigTrans =
302
303 tauLoose = tau->isTau(xAOD::TauJetParameters::JetRNNSigLoose);
304 tauMedium = tau->isTau(xAOD::TauJetParameters::JetRNNSigMedium);
305 tauTight = tau->isTau(xAOD::TauJetParameters::JetRNNSigTight);
306 }
307
308 dRmax = tau->detail<float>(xAOD::TauJetParameters::dRmax);
312
313 panEta = tau->etaPanTauCellBased();
314 panPhi = tau->phiPanTauCellBased();
315 panPt = tau->ptPanTauCellBased() / GeV; // GeV ;
316
317 // TauB/Identification/EleVetoBDTinputs
318 PSSFrac = tau->detail<float>(xAOD::TauJetParameters::PSSFraction);
319 static const SG::ConstAccessor<float> EMFracFixedAcc("EMFracFixed");
320 EMFrac = EMFracFixedAcc(*tau);
321
322 // TauB/SubStructure
323 EMFracTrk = tau->detail<float>(xAOD::TauJetParameters::ChPiEMEOverCaloEME);
324 nNeutPFO = tau->nProtoNeutralPFOs();
325 nShot = tau->nShotPFOs();
326
327 int panModeDummy = -1;
328 int panModeSubstructureDummy = -1;
329
330 if (m_etaMin < std::abs(tauEta) && std::abs(tauEta) < m_etaMax) {
331
332 if( tauEt > lowerEtThreshold) {
333 nTauCandidates += 1;
334 }
335
336 if (tauEt > higherEtThreshold) {
337 nHighPtTauCandidates += 1;
338 nHighPtTaus += 1;
339 }
340
341 if (m_kinGroupName != "tauMonKinGroupGlobal" &&
342 tauEt > lowerEtThreshold && tauLoose) {
343
344 tauPhiEt15Loose = tau->phi();
345 tauEtaEt15Loose = tau->eta();
346 tauEtEt15Loose = tau->pt() / GeV;
347 nClustersEt15Loose =
349 NumTracksEt15Loose = tau->nTracks();
350
352 panModeDummy);
353 panModeEt15Loose = panModeDummy;
354 fill(tool, tauPhiEt15Loose, tauEtaEt15Loose,
355 nClustersEt15Loose, NumTracksEt15Loose, tauEtEt15Loose,
356 panModeEt15Loose);
357 }
358
359 if (m_kinGroupName != "tauMonKinGroupGlobal" && tauLoose) {
360 tauPhiLoose = tau->phi();
361 tauEtaLoose = tau->eta();
362 tauEtLoose = tau->pt() / GeV;
363 NumTracksLoose = tau->nTracks();
364
365 fill(tool, tauPhiLoose, tauEtaLoose, NumTracksLoose,
366 tauEtLoose);
367 }
368
369 if (m_kinGroupName != "tauMonKinGroupGlobal" && tauMedium) {
370 tauPhiMedium = tau->phi();
371 tauEtaMedium = tau->eta();
372 tauEtMedium = tau->pt() / GeV;
373 NumTracksMedium = tau->nTracks();
374
375 fill(tool, tauPhiMedium, tauEtaMedium, NumTracksMedium,
376 tauEtMedium);
377 }
378
379 if (m_kinGroupName != "tauMonKinGroupGlobal" && tauTight) {
380 tauPhiTight = tau->phi();
381 tauEtaTight = tau->eta();
382 tauEtTight = tau->pt() / GeV;
383 NumTracksTight = tau->nTracks();
384
385 fill(tool, tauPhiTight, tauEtaTight, NumTracksTight,
386 tauEtTight);
387 }
388
389
390 // tracks - old
391 if (tau->nTracks() != 0) {
392
393 massTrkSys =
394 tau->detail<float>(xAOD::TauJetParameters::massTrkSys) / GeV; // GeV
396 tau->detail<float>(xAOD::TauJetParameters::trFlightPathSig);
397
398 static const SG::ConstAccessor<float> d0SigTJVAAcc("d0SigTJVA");
399 if (d0SigTJVAAcc.isAvailable(*tau->track(0))) {
400 ipSigLeadTrk = tau->track(0)->d0SigTJVA();
401 }
402 else {
404 }
405
406 static const SG::ConstAccessor<float> z0sinthetaSigTJVAAcc("z0sinthetaSigTJVA");
407 if (z0sinthetaSigTJVAAcc.isAvailable(*tau->track(0))) {
408 ipZ0SinThetaSigLeadTrk = tau->track(0)->z0sinthetaSigTJVA();
409 }
410 else {
412 }
413
414 etOverPtLeadTrack =
415 tau->detail<float>(xAOD::TauJetParameters::etOverPtLeadTrk);
416 leadTrkPt = tau->detail<float>(xAOD::TauJetParameters::leadTrkPt) / GeV;
417
418 fill(tool, massTrkSys, etOverPtLeadTrack, trFlightPathSig,
419 ipSigLeadTrk, ipZ0SinThetaSigLeadTrk, leadTrkPt);
420
421 if (environment() != Environment_t::AOD) {
422 const xAOD::TrackParticle *track = tau->track(0)->track();
423 const Trk::Perigee perigee = track->perigeeParameters();
424
425 uint8_t dummy(0);
426
427 if (track->summaryValue(dummy, xAOD::numberOfSCTSharedHits)) {
429 fill(tool, numberOfSCTSharedHits);
430 }
431
432 if (track->summaryValue(dummy, xAOD::numberOfSCTHits)) {
434 fill(tool, numberOfSCTHits);
435 }
436
437 if (track->summaryValue(dummy, xAOD::numberOfPixelSharedHits)) {
439 fill(tool, numberOfPixelSharedHits);
440 }
441
442 if (track->summaryValue(dummy, xAOD::numberOfPixelHits)) {
444 fill(tool, numberOfPixelHits);
445 }
446 if (track->summaryValue(dummy, xAOD::numberOfTRTHighThresholdHits)) {
448 fill(tool, numberOfTRTHighThresholdHits);
449 }
450 if (track->summaryValue(dummy,
453 fill(tool, numberOfTRTHighThresholdOutliers);
454 }
455 if (track->summaryValue(dummy, xAOD::numberOfTRTHits)) {
457 fill(tool, numberOfTRTHits);
458 }
459 if (track->summaryValue(dummy, xAOD::numberOfTRTOutliers)) {
461 fill(tool, numberOfTRTOutliers);
462 }
463
464 d0 = perigee.parameters()[Trk::d0];
465 z0 = perigee.parameters()[Trk::z0];
466
467 phiTrack = perigee.parameters()[Trk::phi];
468 etaTrack = perigee.eta();
469 ptTrack = perigee.pT() / GeV;
470
471 fill(tool, d0, z0, phiTrack, etaTrack, ptTrack);
472 }
473 }
474 // this else can be removed, but it sets any track variable to 0 if
475 // there are no tracks this solution makes entry numbers match calo which
476 // is desired but there are too many zeros.
477 else {
478 leadTrkPt = 0;
479 fill(tool, leadTrkPt);
480 }
481
482 // Code for All Tracks and Mean track variables.
483
484 // pre loop variables:
485 double tauSeedPt = tau->ptJetSeed();
486 double logTauSeedPt = std::log(tauSeedPt);
487 trackLogSeedJetPt = logTauSeedPt;
488
489 fill(tool, trackLogSeedJetPt);
490
491 static const SG::ConstAccessor<float> acc_trackScoreCharged("rnn_chargedScore");
492 static const SG::ConstAccessor<float> acc_trackScoreIso("rnn_isolationScore");
493 static const SG::ConstAccessor<float> acc_trackScoreConv("rnn_conversionScore");
494 // rnn_fakeScore may not be available (it is not provided by the TauJets smart slimming list), it can be obtained from unitarity
495
496 for (const xAOD::TauTrack *track : tau->allTracks()) {
497
498 static const SG::Accessor<
500 trackAcc("trackLinks");
501 if (!trackAcc(*track)[0]) {
502 continue;
503 }
504
505 const xAOD::TrackParticle *trackParticle = track->track();
506
507 float d0TJVA = track->track()->d0();
508 trackZ0SinthetaTJVA = track->z0sinThetaTJVA(*tau);
509 trackD0SigTJVA = 999.;
510 trackZ0sinthetaSigTJVA = 999.;
511 float rConv = 999.;
512 float rConvII = 999.;
513 static const SG::ConstAccessor<float> z0sinthetaTJVAAcc("z0sinthetaTJVA");
514 if (z0sinthetaTJVAAcc.isAvailable(*track)) {
515 d0TJVA = track->d0TJVA();
516 trackZ0SinthetaTJVA = track->z0sinthetaTJVA();
517 trackD0SigTJVA = track->d0SigTJVA();
518 trackZ0sinthetaSigTJVA = track->z0sinthetaSigTJVA();
519 rConv = track->rConv();
520 rConvII = track->rConvII();
521 }
522 trackdRJetSeedAxis =
523 track->p4().DeltaR(tau->p4(xAOD::TauJetParameters::JetSeed));
524 double qOverP = trackParticle->qOverP();
525 double trackPt = trackParticle->pt();
526 uint8_t nInnermostPixelLayerHits = 0;
527 trackParticle->summaryValue(nInnermostPixelLayerHits,
530 trackParticle->summaryValue(nPixelHits, xAOD::numberOfPixelHits);
532 trackParticle->summaryValue(nPixelDeadSensors,
534 uint8_t nSCTHits = 0;
535 trackParticle->summaryValue(nSCTHits, xAOD::numberOfSCTHits);
537 trackParticle->summaryValue(nSCTDeadSensors,
540 trackParticle->summaryValue(nTRTHighThresholdHits,
543 trackParticle->summaryValue(numberOfPixelHoles,
546 trackParticle->summaryValue(numberOfSCTHoles, xAOD::numberOfSCTHoles);
547 float eProbabilityHT = 0.;
548 trackParticle->summaryValue(eProbabilityHT, xAOD::eProbabilityHT);
549 static const SG::ConstAccessor<float> eProbabilityNNAcc("eProbabilityNN");
550 float eProbabilityNN = eProbabilityNNAcc.withDefault(*trackParticle, -1);
551 // hybrid variable (eProbabilityNN is not computed for tracks with pt
552 // < 2 GeV)
553 trackeProbabilityHTorNN =
554 (trackPt > 2000.) ? eProbabilityNN : eProbabilityHT;
555
556 trackeProbabilityNN = eProbabilityNN;
557 trackeProbabilityHT = eProbabilityHT;
558 trackLogPt = std::log(trackPt);
559 trackEta = track->eta();
560 trackd0TJVA = std::tanh(d0TJVA / 10.);
561 trackCharge = trackParticle->charge();
562 trackqOverP = qOverP * 1000.;
563 trackLogRConv = std::log(rConv);
564 trackTanhRConvII = std::tanh(rConvII / 500.0);
565 trackPtRatioSeedJet = trackPt / tauSeedPt;
566 trackNInnermostPixHits = nInnermostPixelLayerHits;
567 trackNPixHits = nPixelHits + nPixelDeadSensors;
568 trackNSiHits =
570
571 fill(tool, trackLogPt, trackEta, trackd0TJVA, trackZ0SinthetaTJVA,
572 trackD0SigTJVA, trackZ0sinthetaSigTJVA, trackCharge, trackqOverP,
573 trackLogRConv, trackTanhRConvII, trackPtRatioSeedJet,
574 trackdRJetSeedAxis, trackNInnermostPixHits,
575 trackNPixHits,
576 trackNSiHits,
577 trackeProbabilityHT, trackeProbabilityNN, trackeProbabilityHTorNN);
578
579 if (acc_trackScoreCharged.isAvailable(*track)) {
580 float chargedScore = acc_trackScoreCharged(*track);
581 float isolationScore = acc_trackScoreIso(*track);
582 float conversionScore = acc_trackScoreConv(*track);
583 float fakeScore = 1. - chargedScore - isolationScore - conversionScore;
584 // ensure the probability is within [0.,1.]
585 fakeScore = std::max(0.f, fakeScore);
586 fakeScore = std::min(1.f, fakeScore);
587
588 trackIdScoreCharged = chargedScore;
589 trackIdScoreIso = isolationScore;
590 trackIdScoreConv = conversionScore;
591 trackIdScoreFake = fakeScore;
592
593 fill(tool, trackIdScoreCharged, trackIdScoreIso, trackIdScoreConv,
594 trackIdScoreFake);
595 }
596 }
597
598 // clusters
599 std::vector<const xAOD::IParticle *> particleList = tau->clusters();
600 std::vector<xAOD::CaloVertexedTopoCluster> clusters;
601 const xAOD::Vertex *vertex = nullptr;
602 if (tau->vertexLink().isValid())
603 vertex = tau->vertex();
604
605 for (const xAOD::IParticle *particle : particleList) {
606 const xAOD::CaloCluster *cluster =
607 static_cast<const xAOD::CaloCluster *>(particle);
608 if (vertex) {
609 clusters.emplace_back(*cluster, xAOD::CaloCluster::State::CALIBRATED,
610 vertex->position());
611 } else {
612 clusters.emplace_back(*cluster, xAOD::CaloCluster::State::CALIBRATED);
613 }
614 }
615
616 // sort by decreasing Et
617 auto et_cmp = [](const xAOD::CaloVertexedTopoCluster &lhs,
618 const xAOD::CaloVertexedTopoCluster &rhs) {
619 return lhs.p4().Et() > rhs.p4().Et();
620 };
621 std::sort(clusters.begin(), clusters.end(), et_cmp);
622
623 // keep first 6 leading clusters as in RNN ID
624 if (clusters.size() > 6) {
625 clusters.resize(6, clusters[0]);
626 }
627
628 double moment;
629 for (const auto& vertexedCluster : clusters) {
630 const xAOD::CaloCluster &cluster = vertexedCluster.clust();
631
632 clusterLogEt = std::log10(vertexedCluster.p4().Et());
633 clusterEta = vertexedCluster.eta();
634 clusterPhi = vertexedCluster.phi();
635
637 clusterSecondR = std::log10(moment + 0.1);
638
640 moment);
641 clusterSecondLambda = std::log10(moment + 0.1);
642
644 moment);
645 clusterCenterLambda = std::log10(moment + 1e-6);
646
648 moment);
649 if (moment != 0.)
650 moment = std::log10(std::abs(moment));
651 clusterFirstEngDens = moment;
652
654 moment);
655 clusterEMproba = moment;
656
657 fill(tool, clusterLogEt, clusterEta, clusterPhi, clusterSecondR,
658 clusterSecondLambda, clusterCenterLambda, clusterFirstEngDens,
659 clusterEMproba);
660 }
661
662 // Et-weighted average of cluster moments (MVA TES)
663 float avariable = 0.;
665 avariable);
666 if (test)
667 clustersMeanCenterLambda = avariable;
668
670 avariable);
671 if (test)
672 clustersMeanFirstEngDens = avariable;
673
675 avariable);
676 if (test)
677 clustersMeanEMProbability = avariable;
678
680 avariable);
681 if (test)
682 clustersMeanSecondLambda = avariable;
683
685 avariable);
686 if (test)
687 clustersMeanPresamplerFrac = avariable;
688
689 test = tau->detail(xAOD::TauJetParameters::PFOEngRelDiff, avariable);
690 if (test)
691 clustersPFOEngRelDiff = avariable;
692
693 fill(tool, clustersMeanCenterLambda, clustersMeanFirstEngDens,
694 clustersMeanEMProbability, clustersMeanSecondLambda,
695 clustersMeanPresamplerFrac, clustersPFOEngRelDiff);
696
697 for (unsigned int np = 0; np < nNeutPFO; np++) {
698 const xAOD::PFO *npfo = tau->protoNeutralPFO(np);
699 BDTScoreAsP0 = npfo->bdtPi0Score();
700 fill(tool, BDTScoreAsP0);
701 }
702
703 if (tauEt > lowerEtThreshold) {
704 tauPhiEt15 = tau->phi();
705 tauEtaEt15 = tau->eta();
706 fill(tool, tauPhiEt15, tauEtaEt15);
707 }
708
709 fill(tool, tauPhi, tauEta, LB, tauEt, centFrac, isolFrac, EMRadius,
710 hadRadius, etEMAtEMScale, etHadAtEMScale,
711 tauCharge, JetScore, JetScoreSigTrans, RNNEleScore, RNNEleScoreSigTrans,
712 muonVeto, tauLoose, tauMedium, tauTight, PSSFrac, EMFrac, EMFracTrk, nNeutPFO,
713 nShot, NumTracks, nClusters, jetSeedEta, jetSeedPhi, jetSeedPt,
714 dRmax, ptRatioEflowApprox, trkAvgDist);
715
717 panModeSubstructureDummy);
718 panModeSubstructure = panModeSubstructureDummy;
719
720 fill(tool, panModeSubstructure);
721 if (panPhi > -100) {
722 fill(tool, panEta, panPhi, panPt);
723 }
724 }
725 }
726
727 fill(tool, nHighPtTauCandidates, nTauCandidates, nHighPtTaus);
728
729 return StatusCode::SUCCESS;
730}
#define ATH_MSG_ERROR(x)
static Double_t taus
Environment_t environment() const
Accessor functions for the environment.
SG::ReadHandle< xAOD::EventInfo > GetEventInfo(const EventContext &) const
Return a ReadHandle for an EventInfo object (get run/event numbers, etc.)
Declare a monitored scalar variable.
bool isAvailable()
Test to see if this variable exists in the store, for the referenced object.
double eta() const
Access method for pseudorapidity - from momentum.
double pT() const
Access method for transverse momentum.
SG::ReadDecorHandleKey< xAOD::TauJetContainer > m_offlineGNTauScoreDecorKey
SG::ReadDecorHandleKey< xAOD::TauJetContainer > m_offlineGNTauMediumWPDecorKey
SG::ReadDecorHandleKey< xAOD::TauJetContainer > m_offlineGNTauLooseWPDecorKey
SG::ReadDecorHandleKey< xAOD::TauJetContainer > m_offlineGNTauTransScoreDecorKey
SG::ReadDecorHandleKey< xAOD::TauJetContainer > m_offlineGNTauTightWPDecorKey
bool retrieveMoment(MomentType type, double &value) const
Retrieve individual moment.
@ SECOND_LAMBDA
Second Moment in .
@ EM_PROBABILITY
Classification probability to be em-like.
@ SECOND_R
Second Moment in .
@ FIRST_ENG_DENS
First Moment in E/V.
@ CENTER_LAMBDA
Shower depth at Cluster Centroid.
virtual FourMom_t p4() const final
The full 4-momentum of the particle.
float bdtPi0Score() const
get BDT Score used to classify clusters as Pi0 like or not
std::vector< ElementLink< xAOD::TrackParticleContainer > > TrackParticleLinks_t
Definition TauTrack_v1.h:91
bool summaryValue(uint8_t &value, const SummaryType &information) const
Accessor for TrackSummary values.
float qOverP() const
Returns the parameter.
virtual double pt() const override final
The transverse momentum ( ) of the particle.
float charge() const
Returns the charge.
float nTRTHighThresholdHits(const U &p)
float nPixelDeadSensors(const U &p)
float nSCTHits(const U &p)
float nSCTDeadSensors(const U &p)
float nPixelHits(const U &p)
@ qOverP
perigee
void fill(const ToolHandle< GenericMonitoringTool > &tool, T &&... variables)
double tauEta(const xAOD::TauJet &xTau)
return tau eta
bool trackEta(const xAOD::TauJet &, const xAOD::TauTrack &track, float &out)
bool d0TJVA(const xAOD::TauJet &, const xAOD::TauTrack &track, float &out)
bool trackPt(const xAOD::TauJet &, const xAOD::TauTrack &track, float &out)
constexpr ParticleHypothesis particle[PARTICLEHYPOTHESES]
the array of masses
ParametersT< TrackParametersDim, Charged, PerigeeSurface > Perigee
@ phi
Definition ParamDefs.h:75
@ d0
Definition ParamDefs.h:63
@ z0
Definition ParamDefs.h:64
void sort(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end)
Specialization of sort for DataVector/List.
@ hadRadius
Get hadron calorimeter radius.
Definition TauDefs.h:192
@ EMRadius
Get E_T radius.
Definition TauDefs.h:190
@ etHadAtEMScale
Get Hadronic energy at EM scale.
Definition TauDefs.h:196
@ isolFrac
Get isolation fraction.
Definition TauDefs.h:198
@ trkAvgDist
Get calibrated EM transverse energy (DEPRECATED since r19)
Definition TauDefs.h:214
@ etEMAtEMScale
Get EM energy at EM scale.
Definition TauDefs.h:194
@ centFrac
Get centrality fraction.
Definition TauDefs.h:200
@ numTopoClusters
get number of topocluster constituents of jet associated to tau candidate
Definition TauDefs.h:173
@ dRmax
Get maximal dR of tracks associated to calo-seeded tau.
Definition TauDefs.h:226
@ RNNEleScore
RNN score for Ele rejection (not transformed)
Definition TauDefs.h:94
@ RNNEleScoreSigTrans
RNN score which is signal transformed/flattened.
Definition TauDefs.h:96
@ RNNJetScore
RNN score for Jet rejection (not transformed)
Definition TauDefs.h:90
@ RNNJetScoreSigTrans
RNN score which is signal transformed/flattened.
Definition TauDefs.h:92
PFO_v1 PFO
Definition of the current "pfo version".
Definition PFO.h:17
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
std::pair< std::unique_ptr< T >, std::unique_ptr< ShallowAuxContainer > > shallowCopyContainer(const T &cont, const EventContext &ctx)
Function making a shallow copy of a constant container.
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Vertex_v1 Vertex
Define the latest version of the vertex class.
TauTrack_v1 TauTrack
Definition of the current version.
Definition TauTrack.h:16
@ numberOfPixelHoles
number of pixel layers on track with absence of hits [unit8_t].
@ numberOfTRTHighThresholdOutliers
number of TRT high threshold outliers (only xenon counted) [unit8_t].
@ numberOfTRTHits
number of TRT hits [unit8_t].
@ numberOfSCTDeadSensors
number of dead SCT sensors crossed [unit8_t].
@ eProbabilityHT
Electron probability from High Threshold (HT) information [float].
@ numberOfSCTHits
number of hits in SCT [unit8_t].
@ numberOfInnermostPixelLayerHits
these are the hits in the 0th pixel barrel layer
@ numberOfPixelHits
these are the pixel hits, including the b-layer [unit8_t].
@ numberOfPixelSharedHits
number of Pixel all-layer hits shared by several tracks [unit8_t].
@ numberOfSCTSharedHits
number of SCT hits shared by several tracks [unit8_t].
@ numberOfTRTHighThresholdHits
number of TRT hits which pass the high threshold (only xenon counted) [unit8_t].
@ numberOfTRTOutliers
number of TRT outliers [unit8_t].
@ numberOfPixelDeadSensors
number of dead pixel sensors crossed [unit8_t].
@ numberOfSCTHoles
number of SCT holes [unit8_t].

◆ filterPassed()

virtual bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx) const
inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

96 {
97 return execState( ctx ).filterPassed();
98 }
virtual bool filterPassed(const EventContext &ctx) const

◆ GetEventInfo()

SG::ReadHandle< xAOD::EventInfo > AthMonitorAlgorithm::GetEventInfo ( const EventContext & ctx) const
inherited

Return a ReadHandle for an EventInfo object (get run/event numbers, etc.)

Parameters
ctxEventContext for the event
Returns
a SG::ReadHandle<xAOD::EventInfo>

Definition at line 111 of file AthMonitorAlgorithm.cxx.

111 {
112 return SG::ReadHandle<xAOD::EventInfo>(m_EventInfoKey, ctx);
113}
SG::ReadHandleKey< xAOD::EventInfo > m_EventInfoKey
Key for retrieving EventInfo from StoreGate.

◆ getGroup()

const ToolHandle< GenericMonitoringTool > & AthMonitorAlgorithm::getGroup ( const std::string & name) const
inherited

Get a specific monitoring tool from the tool handle array.

Finds a specific GenericMonitoringTool instance from the list of monitoring tools (a ToolHandleArray). Throws a FATAL warning if the object found is empty.

Parameters
namestring name of the desired tool
Returns
reference to the desired monitoring tool

Definition at line 168 of file AthMonitorAlgorithm.cxx.

168 {
169 // get the pointer to the tool, and check that it exists
170 auto idx = m_toolLookupMap.find(name);
171 if (ATH_LIKELY(idx != m_toolLookupMap.end())) {
172 return m_tools[idx->second];
173 }
174 else {
175 // treat empty tool handle case as in Monitored::Group
176 if (m_toolLookupMap.empty()) {
177 return m_dummy;
178 }
179
180 if (!isInitialized()) {
182 "It seems that the AthMonitorAlgorithm::initialize was not called "
183 "in derived class initialize method");
184 } else {
185 std::string available = std::accumulate(
186 m_toolLookupMap.begin(), m_toolLookupMap.end(), std::string(""),
187 [](const std::string& s, auto h) { return s + "," + h.first; });
188 ATH_MSG_FATAL("The tool " << name << " could not be found in the tool array of the "
189 << "monitoring algorithm " << m_name << ". This probably reflects a discrepancy between "
190 << "your python configuration and c++ filling code. Note: your available groups are {"
191 << available << "}.");
192 }
193 }
194 return m_dummy;
195}
#define ATH_MSG_FATAL(x)
#define ATH_LIKELY(x)
std::unordered_map< std::string, size_t > m_toolLookupMap
const ToolHandle< GenericMonitoringTool > m_dummy
ToolHandleArray< GenericMonitoringTool > m_tools
Array of Generic Monitoring Tools.

◆ getTrigDecisionTool()

const ToolHandle< Trig::TrigDecisionTool > & AthMonitorAlgorithm::getTrigDecisionTool ( ) const
inherited

Get the trigger decision tool member.

The trigger decision tool is used to check whether a specific trigger is passed by an event.

Returns
m_trigDecTool

Definition at line 198 of file AthMonitorAlgorithm.cxx.

198 {
199 return m_trigDecTool;
200}

◆ initialize()

StatusCode tauMonitorAlgorithm::initialize ( )
overridevirtual

initialize

Returns
StatusCode

Reimplemented from AthMonitorAlgorithm.

Definition at line 22 of file tauMonitorAlgorithm.cxx.

22 {
23
24 ATH_CHECK(m_TauContainerKey.initialize());
25
26 m_offlineGNTauScoreDecorKey = m_TauContainerKey.key() + "." + m_offlineGNTauScoreDecorKey.key();
27 ATH_CHECK(m_offlineGNTauScoreDecorKey.initialize());
28
29 m_offlineGNTauTransScoreDecorKey = m_TauContainerKey.key() + "." + m_offlineGNTauTransScoreDecorKey.key();
30 ATH_CHECK(m_offlineGNTauTransScoreDecorKey.initialize());
31
32 m_offlineGNTauLooseWPDecorKey = m_TauContainerKey.key() + "." + m_offlineGNTauLooseWPDecorKey.key();
33 ATH_CHECK(m_offlineGNTauLooseWPDecorKey.initialize());
34
35 m_offlineGNTauMediumWPDecorKey = m_TauContainerKey.key() + "." + m_offlineGNTauMediumWPDecorKey.key();
36 ATH_CHECK(m_offlineGNTauMediumWPDecorKey.initialize());
37
38 m_offlineGNTauTightWPDecorKey = m_TauContainerKey.key() + "." + m_offlineGNTauTightWPDecorKey.key();
39 ATH_CHECK(m_offlineGNTauTightWPDecorKey.initialize());
40
42}
#define ATH_CHECK
Evaluate an expression and check for errors.
virtual StatusCode initialize() override
initialize

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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.

◆ isClonable()

◆ msg()

MsgStream & AthCommonMsg< Gaudi::Algorithm >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

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

◆ msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::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< Gaudi::Algorithm > >::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.

◆ parseList()

StatusCode AthMonitorAlgorithm::parseList ( const std::string & line,
std::vector< std::string > & result ) const
virtualinherited

Parse a string into a vector.

The input string is a single long string of all of the trigger names. It parses this string and turns it into a vector, where each element is one trigger or trigger category.

Parameters
lineThe input string.
resultThe parsed output vector of strings.
Returns
StatusCode

Definition at line 345 of file AthMonitorAlgorithm.cxx.

345 {
346 std::string item;
347 std::stringstream ss(line);
348
349 ATH_MSG_DEBUG( "AthMonitorAlgorithm::parseList()" );
350
351 while ( std::getline(ss, item, ',') ) {
352 std::stringstream iss(item); // remove whitespace
353 iss >> item;
354 result.push_back(item);
355 }
356
357 return StatusCode::SUCCESS;
358}
static Double_t ss

◆ 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< Gaudi::Algorithm > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
384 }
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)

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ setFilterPassed()

virtual void AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state,
const EventContext & ctx ) const
inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

100 {
102 }
virtual void setFilterPassed(bool state, const EventContext &ctx) const

◆ sysExecute()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx)
overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

77{
78 return BaseAlg::sysExecute (ctx);
79}

◆ sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysInitialize ( )
overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

107 {
109
110 if (sc.isFailure()) {
111 return sc;
112 }
113
114 ServiceHandle<ICondSvc> cs("CondSvc",name());
115 for (auto h : outputHandles()) {
116 if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
117 // do this inside the loop so we don't create the CondSvc until needed
118 if ( cs.retrieve().isFailure() ) {
119 ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
120 return StatusCode::SUCCESS;
121 }
122 if (cs->regHandle(this,*h).isFailure()) {
124 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
125 << " with CondSvc");
126 }
127 }
128 }
129 return sc;
130}
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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.

◆ trigChainsArePassed()

bool AthMonitorAlgorithm::trigChainsArePassed ( const std::vector< std::string > & vTrigNames) const
inherited

Check whether triggers are passed.

For the event, use the trigger decision tool to check that at least one of the triggers listed in the supplied vector is passed.

Parameters
vTrigNamesList of trigger names.
Returns
If empty input, default to true. If at least one trigger is specified, returns whether at least one trigger was passed.

Definition at line 203 of file AthMonitorAlgorithm.cxx.

203 {
204
205
206 // If no triggers were given, return true.
207 if (vTrigNames.empty()) return true;
208
209
210 // Trigger: Check if this Algorithm is being run as an Express Stream job.
211 // Events are entering the express stream are chosen randomly, and by chain,
212 // Hence an additional check should be aplied to see if the chain(s)
213 // monitored here are responsible for the event being selected for
214 // the express stream.
215
216 const auto group = m_trigDecTool->getChainGroup(vTrigNames);
218 const auto passedBits = m_trigDecTool->isPassedBits(group);
219 bool expressPass = passedBits & TrigDefs::Express_passed; //bitwise AND
220 if(!expressPass) {
221 return false;
222 }
223 }
224
225 // monitor the event if any of the chains in the chain group passes the event.
226 return group->isPassed();
227
228}
Gaudi::Property< bool > m_enforceExpressTriggers

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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 }

Member Data Documentation

◆ cfg

tauMonitorAlgorithm.cfg = MainServicesCfg(flags)

Definition at line 640 of file tauMonitorAlgorithm.py.

◆ exampleMonitorAcc

tauMonitorAlgorithm.exampleMonitorAcc = tauMonitoringConfig(flags)

Definition at line 643 of file tauMonitorAlgorithm.py.

◆ file

str tauMonitorAlgorithm.file = 'data16_13TeV.00311321.physics_Main.recon.AOD.r9264/AOD.11038520._000001.pool.root.1'

Definition at line 629 of file tauMonitorAlgorithm.py.

◆ Files

tauMonitorAlgorithm.Files

Definition at line 630 of file tauMonitorAlgorithm.py.

◆ flags

tauMonitorAlgorithm.flags = initConfigFlags()

Definition at line 626 of file tauMonitorAlgorithm.py.

◆ HISTFileName

tauMonitorAlgorithm.HISTFileName

Definition at line 633 of file tauMonitorAlgorithm.py.

◆ isMC

tauMonitorAlgorithm.isMC

Definition at line 632 of file tauMonitorAlgorithm.py.

◆ m_abGroups1

std::vector<int> tauMonitorAlgorithm::m_abGroups1
private

Definition at line 25 of file tauMonitorAlgorithm.h.

◆ m_abGroups2

std::vector<std::vector<int> > tauMonitorAlgorithm::m_abGroups2
private

Definition at line 26 of file tauMonitorAlgorithm.h.

◆ m_cGroups1

std::map<std::string,int> tauMonitorAlgorithm::m_cGroups1
private

Definition at line 27 of file tauMonitorAlgorithm.h.

◆ m_cGroups2

std::map<std::string,std::map<std::string,int> > tauMonitorAlgorithm::m_cGroups2
private

Definition at line 28 of file tauMonitorAlgorithm.h.

◆ m_dataType

AthMonitorAlgorithm::DataType_t AthMonitorAlgorithm::m_dataType
protectedinherited

Instance of the DataType_t enum.

Definition at line 356 of file AthMonitorAlgorithm.h.

◆ m_dataTypeStr

Gaudi::Property<std::string> AthMonitorAlgorithm::m_dataTypeStr {this,"DataType","userDefined"}
protectedinherited

DataType string pulled from the job option and converted to enum.

Definition at line 358 of file AthMonitorAlgorithm.h.

358{this,"DataType","userDefined"};

◆ m_defaultLBDuration

Gaudi::Property<float> AthMonitorAlgorithm::m_defaultLBDuration {this,"DefaultLBDuration",60.}
protectedinherited

Default duration of one lumi block.

Definition at line 365 of file AthMonitorAlgorithm.h.

365{this,"DefaultLBDuration",60.};

◆ m_detailLevel

Gaudi::Property<int> AthMonitorAlgorithm::m_detailLevel {this,"DetailLevel",0}
protectedinherited

Sets the level of detail used in the monitoring.

Definition at line 366 of file AthMonitorAlgorithm.h.

366{this,"DetailLevel",0};

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_DQFilterTools

ToolHandleArray<IDQFilterTool> AthMonitorAlgorithm::m_DQFilterTools {this,"FilterTools",{}}
protectedinherited

Array of Data Quality filter tools.

Definition at line 346 of file AthMonitorAlgorithm.h.

346{this,"FilterTools",{}};

◆ m_dummy

const ToolHandle<GenericMonitoringTool> AthMonitorAlgorithm::m_dummy
privateinherited

Definition at line 374 of file AthMonitorAlgorithm.h.

◆ m_enforceExpressTriggers

Gaudi::Property<bool> AthMonitorAlgorithm::m_enforceExpressTriggers
privateinherited
Initial value:
{this,
"EnforceExpressTriggers", false,
"Requires that matched triggers made the event enter the express stream"}

Definition at line 377 of file AthMonitorAlgorithm.h.

377 {this,
378 "EnforceExpressTriggers", false,
379 "Requires that matched triggers made the event enter the express stream"};

◆ m_environment

AthMonitorAlgorithm::Environment_t AthMonitorAlgorithm::m_environment
protectedinherited

Instance of the Environment_t enum.

Definition at line 355 of file AthMonitorAlgorithm.h.

◆ m_environmentStr

Gaudi::Property<std::string> AthMonitorAlgorithm::m_environmentStr {this,"Environment","user"}
protectedinherited

Environment string pulled from the job option and converted to enum.

Definition at line 357 of file AthMonitorAlgorithm.h.

357{this,"Environment","user"};

◆ m_etaMax

Gaudi::Property<float> tauMonitorAlgorithm::m_etaMax {this, "etaMax", 3.0}
private

Definition at line 38 of file tauMonitorAlgorithm.h.

38{this, "etaMax", 3.0};

◆ m_etaMin

Gaudi::Property<float> tauMonitorAlgorithm::m_etaMin {this, "etaMin", -1.}
private

Definition at line 37 of file tauMonitorAlgorithm.h.

37{this, "etaMin", -1.};

◆ m_EventInfoKey

SG::ReadHandleKey<xAOD::EventInfo> AthMonitorAlgorithm::m_EventInfoKey {this,"EventInfoKey","EventInfo"}
protectedinherited

Key for retrieving EventInfo from StoreGate.

Definition at line 367 of file AthMonitorAlgorithm.h.

367{this,"EventInfoKey","EventInfo"};

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

◆ m_fileKey

Gaudi::Property<std::string> AthMonitorAlgorithm::m_fileKey {this,"FileKey",""}
protectedinherited

Internal Athena name for file.

Definition at line 363 of file AthMonitorAlgorithm.h.

363{this,"FileKey",""};

◆ m_kinGroupName

Gaudi::Property<std::string> tauMonitorAlgorithm::m_kinGroupName {this, "kinGroupName", "tauMonKinGroupBA"}
private

Definition at line 40 of file tauMonitorAlgorithm.h.

40{this, "kinGroupName", "tauMonKinGroupBA"};

◆ m_lbDurationDataKey

SG::ReadCondHandleKey<LBDurationCondData> AthMonitorAlgorithm::m_lbDurationDataKey {this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"}
protectedinherited

Definition at line 350 of file AthMonitorAlgorithm.h.

351{this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"};

◆ m_lumiDataKey

SG::ReadCondHandleKey<LuminosityCondData> AthMonitorAlgorithm::m_lumiDataKey {this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"}
protectedinherited

Definition at line 348 of file AthMonitorAlgorithm.h.

349{this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"};

◆ m_name

std::string AthMonitorAlgorithm::m_name
privateinherited

Definition at line 371 of file AthMonitorAlgorithm.h.

◆ m_offlineGNTauLooseWPDecorKey

SG::ReadDecorHandleKey<xAOD::TauJetContainer> tauMonitorAlgorithm::m_offlineGNTauLooseWPDecorKey {this, "OfflineGNTauLooseWPDecorKey", "GNTauL_v0prune", "Offline GNTau Loose WP decoration key"}
private

Definition at line 33 of file tauMonitorAlgorithm.h.

33{this, "OfflineGNTauLooseWPDecorKey", "GNTauL_v0prune", "Offline GNTau Loose WP decoration key"};

◆ m_offlineGNTauMediumWPDecorKey

SG::ReadDecorHandleKey<xAOD::TauJetContainer> tauMonitorAlgorithm::m_offlineGNTauMediumWPDecorKey {this, "OfflineGNTauMediumWPDecorKey", "GNTauM_v0prune", "Offline GNTau Medium WP decoration key"}
private

Definition at line 34 of file tauMonitorAlgorithm.h.

34{this, "OfflineGNTauMediumWPDecorKey", "GNTauM_v0prune", "Offline GNTau Medium WP decoration key"};

◆ m_offlineGNTauScoreDecorKey

SG::ReadDecorHandleKey<xAOD::TauJetContainer> tauMonitorAlgorithm::m_offlineGNTauScoreDecorKey {this, "OfflineGNTauScoreDecorKey", "GNTauScore_v0prune", "Offline GNTau Score decoration key"}
private

Definition at line 31 of file tauMonitorAlgorithm.h.

31{this, "OfflineGNTauScoreDecorKey", "GNTauScore_v0prune", "Offline GNTau Score decoration key"};

◆ m_offlineGNTauTightWPDecorKey

SG::ReadDecorHandleKey<xAOD::TauJetContainer> tauMonitorAlgorithm::m_offlineGNTauTightWPDecorKey {this, "OfflineGNTauTightWPDecorKey", "GNTauT_v0prune", "Offline GNTau Tight WP decoration key"}
private

Definition at line 35 of file tauMonitorAlgorithm.h.

35{this, "OfflineGNTauTightWPDecorKey", "GNTauT_v0prune", "Offline GNTau Tight WP decoration key"};

◆ m_offlineGNTauTransScoreDecorKey

SG::ReadDecorHandleKey<xAOD::TauJetContainer> tauMonitorAlgorithm::m_offlineGNTauTransScoreDecorKey {this, "OfflineGNTauTransScoreDecorKey", "GNTauScoreSigTrans_v0prune", "Offline GNTau Trans Score decoration key"}
private

Definition at line 32 of file tauMonitorAlgorithm.h.

32{this, "OfflineGNTauTransScoreDecorKey", "GNTauScoreSigTrans_v0prune", "Offline GNTau Trans Score decoration key"};

◆ m_TauContainerKey

SG::ReadHandleKey<xAOD::TauJetContainer> tauMonitorAlgorithm::m_TauContainerKey {this, "TauRecContainer", "TauJets"}
private

Definition at line 30 of file tauMonitorAlgorithm.h.

30{this, "TauRecContainer", "TauJets"};

◆ m_toolLookupMap

std::unordered_map<std::string, size_t> AthMonitorAlgorithm::m_toolLookupMap
privateinherited

Definition at line 372 of file AthMonitorAlgorithm.h.

◆ m_tools

ToolHandleArray<GenericMonitoringTool> AthMonitorAlgorithm::m_tools {this,"GMTools",{}}
protectedinherited

Array of Generic Monitoring Tools.

Definition at line 341 of file AthMonitorAlgorithm.h.

341{this,"GMTools",{}};

◆ m_trigDecTool

PublicToolHandle<Trig::TrigDecisionTool> AthMonitorAlgorithm::m_trigDecTool
protectedinherited

Tool to tell whether a specific trigger is passed.

Definition at line 345 of file AthMonitorAlgorithm.h.

◆ m_triggerChainString

Gaudi::Property<std::string> AthMonitorAlgorithm::m_triggerChainString {this,"TriggerChain",""}
protectedinherited

Trigger chain string pulled from the job option and parsed into a vector.

Definition at line 360 of file AthMonitorAlgorithm.h.

360{this,"TriggerChain",""};

◆ m_trigLiveFractionDataKey

SG::ReadCondHandleKey<TrigLiveFractionCondData> AthMonitorAlgorithm::m_trigLiveFractionDataKey {this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"}
protectedinherited

Definition at line 352 of file AthMonitorAlgorithm.h.

353{this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"};

◆ m_useLumi

Gaudi::Property<bool> AthMonitorAlgorithm::m_useLumi {this,"EnableLumi",false}
protectedinherited

Allows use of various luminosity functions.

Definition at line 364 of file AthMonitorAlgorithm.h.

364{this,"EnableLumi",false};

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.

◆ m_vTrigChainNames

std::vector<std::string> AthMonitorAlgorithm::m_vTrigChainNames
protectedinherited

Vector of trigger chain names parsed from trigger chain string.

Definition at line 361 of file AthMonitorAlgorithm.h.

◆ nightly

str tauMonitorAlgorithm.nightly = '/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/CommonInputs/'

Definition at line 628 of file tauMonitorAlgorithm.py.

◆ OutputLevel

tauMonitorAlgorithm.OutputLevel

Definition at line 647 of file tauMonitorAlgorithm.py.

◆ withDetails

tauMonitorAlgorithm.withDetails

Definition at line 662 of file tauMonitorAlgorithm.py.


The documentation for this class was generated from the following files: