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

#include <MvaTESVariableDecorator.h>

Inheritance diagram for MvaTESVariableDecorator:
Collaboration diagram for MvaTESVariableDecorator:

Public Member Functions

 MvaTESVariableDecorator (const std::string &name="MvaTESVariableDecorator")
virtual ~MvaTESVariableDecorator ()=default
virtual StatusCode initialize () override
 Tool initializer.
virtual StatusCode execute (xAOD::TauJet &xTau) const override
 Execute - called for each tau candidate.
virtual StatusCode eventInitialize () override
 Event initializer - called at the beginning of each event.
virtual StatusCode executePi0CreateROI (xAOD::TauJet &pTau, CaloConstCellContainer &caloCellContainer, boost::dynamic_bitset<> &map) const override
virtual StatusCode executeVertexFinder (xAOD::TauJet &pTau, const xAOD::VertexContainer *vertexContainer=nullptr) const override
virtual StatusCode executeTrackFinder (xAOD::TauJet &pTau, xAOD::TauTrackContainer &tauTrackContainer) const override
virtual StatusCode executeTrackClassifier (xAOD::TauJet &pTau, xAOD::TauTrackContainer &tauTrackContainer) const override
virtual StatusCode executeShotFinder (xAOD::TauJet &pTau, xAOD::CaloClusterContainer &shotClusterContainer, xAOD::PFOContainer &PFOContainer) const override
virtual StatusCode executePi0ClusterCreator (xAOD::TauJet &pTau, xAOD::PFOContainer &neutralPFOContainer, xAOD::PFOContainer &hadronicPFOContainer, const xAOD::CaloClusterContainer &pCaloClusterContainer) const override
virtual StatusCode executeVertexVariables (xAOD::TauJet &pTau, xAOD::VertexContainer &vertexContainer) const override
virtual StatusCode executePi0ClusterScaler (xAOD::TauJet &pTau, xAOD::PFOContainer &neutralPFOContainer, xAOD::PFOContainer &chargedPFOContainer) const override
virtual StatusCode executePi0nPFO (xAOD::TauJet &pTau, xAOD::PFOContainer &neutralPFOContainer) const override
virtual StatusCode executePanTau (xAOD::TauJet &pTau, xAOD::ParticleContainer &particleContainer, xAOD::PFOContainer &neutralPFOContainer) const override
virtual StatusCode eventFinalize () override
 Event finalizer - called at the end of each event.
virtual StatusCode finalize () override
 Finalizer.
std::string find_file (const std::string &fname) const
virtual StatusCode readConfig () override
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
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 Member Functions

bool inTrigger () const
bool inAOD () const
bool inEleRM () const
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

Gaudi::Property< bool > m_in_trigger {this, "inTrigger", false, "Indicate if the tool is running on trigger"}
Gaudi::Property< bool > m_in_AOD {this, "inAOD", false, "Indicate if the tool is running on AOD"}
Gaudi::Property< bool > m_in_EleRM {this, "inEleRM", false, "Indicate if the tool is running on EleRM routine"}
Gaudi::Property< std::string > m_tauRecToolsTag {this, "calibFolder", "tauRecTools/R22_preprod", "CVMFS path to the tau calibration folder"}

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

Gaudi::Property< bool > m_doVertexCorrection {this, "VertexCorrection", true}
SG::ReadDecorHandleKey< xAOD::EventInfom_aveIntPerXKey
SG::ReadHandleKey< xAOD::VertexContainerm_vertexContainerKey
SG::ReadHandleKey< xAOD::EventShapem_eventShapeKey
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 16 of file MvaTESVariableDecorator.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ MvaTESVariableDecorator()

MvaTESVariableDecorator::MvaTESVariableDecorator ( const std::string & name = "MvaTESVariableDecorator")

Definition at line 13 of file MvaTESVariableDecorator.cxx.

14 : TauRecToolBase(name) {
15}
TauRecToolBase(const std::string &name)

◆ ~MvaTESVariableDecorator()

virtual MvaTESVariableDecorator::~MvaTESVariableDecorator ( )
virtualdefault

Member Function Documentation

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

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

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

◆ eventFinalize()

StatusCode TauRecToolBase::eventFinalize ( )
overridevirtualinherited

Event finalizer - called at the end of each event.

Implements ITauToolBase.

Definition at line 202 of file TauRecToolBase.cxx.

202 {
203 return StatusCode::SUCCESS;
204}

◆ eventInitialize()

StatusCode TauRecToolBase::eventInitialize ( )
overridevirtualinherited

Event initializer - called at the beginning of each event.

Implements ITauToolBase.

Definition at line 133 of file TauRecToolBase.cxx.

133 {
134 return StatusCode::SUCCESS;
135}

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

◆ execute()

StatusCode MvaTESVariableDecorator::execute ( xAOD::TauJet & pTau) const
overridevirtual

Execute - called for each tau candidate.

Reimplemented from TauRecToolBase.

Definition at line 29 of file MvaTESVariableDecorator.cxx.

29 {
30 // Tell clang to optimize assuming that FP operations may trap.
32
33 int mu = 0;
34 SG::ReadDecorHandle<xAOD::EventInfo, float> eventInfoDecorHandle( m_aveIntPerXKey );
35 if (!eventInfoDecorHandle.isPresent()) {
36 ATH_MSG_WARNING ( "EventInfo decoration not available! Will set mu=0." );
37 }
38 else {
39 // convert from float to int to ignore peculiar values used in MC
40 mu = static_cast<int>(eventInfoDecorHandle(0));
41 }
42 static const SG::Accessor<float> acc_mu("mu");
43 acc_mu(xTau) = mu;
44
45 if (!m_vertexContainerKey.empty()) {
46 int nVtxPU = 0;
47 SG::ReadHandle<xAOD::VertexContainer> vertexInHandle( m_vertexContainerKey );
48 if (!vertexInHandle.isValid()) {
49 ATH_MSG_WARNING ("Could not retrieve HiveDataObj with key " << vertexInHandle.key() << ", will set nVtxPU=0.");
50 }
51 else {
52 const xAOD::VertexContainer* vertexContainer = vertexInHandle.cptr();
53 for (const auto *xVertex : *vertexContainer){
54 if (xVertex->vertexType() == xAOD::VxType::PileUp)
55 ++nVtxPU;
56 }
57 }
58 static const SG::Accessor<int> acc_nVtxPU("nVtxPU");
59 acc_nVtxPU(xTau) = nVtxPU;
60 }
61
62 if (!m_eventShapeKey.empty()) {
63 double rho = 0.;
64 SG::ReadHandle<xAOD::EventShape> eventShape(m_eventShapeKey);
65 if (!eventShape.isValid()) {
66 ATH_MSG_WARNING ("Could not retrieve EventShape with key " << m_eventShapeKey );
67 }
68 else if (!eventShape->getDensity(xAOD::EventShape::Density, rho)) {
69 ATH_MSG_WARNING ("Could not retrieve rho.");
70 }
71 static const SG::Accessor<float> acc_rho("rho");
72 acc_rho(xTau) = static_cast<float>(rho);
73 }
74
75 double center_lambda=0. , first_eng_dens=0. , em_probability=0. , second_lambda=0. ;
76 double mean_center_lambda=0. , mean_first_eng_dens=0. , mean_em_probability=0. , mean_second_lambda=0. ;
77 double mean_presampler_frac=0., lead_cluster_frac=0. , second_cluster_frac=0. , third_cluster_frac=0. ;
78 double clE=0., Etot=0.;
79
80 // approximate Upsilon based on clusters, not PFOs (for online TES)
81 TLorentzVector clusters_EM_P4;
82 clusters_EM_P4.SetPtEtaPhiM(0,0,0,0);
83 TLorentzVector clusters_had_P4;
84 clusters_had_P4.SetPtEtaPhiM(0,0,0,0);
85 TLorentzVector tauIntermediateAxisEM;
86 tauIntermediateAxisEM.SetPtEtaPhiM(0,0,0,0);
87
88 TLorentzVector tauAxis = tauRecTools::getTauAxis(xTau, m_doVertexCorrection);
89
90 // in the trigger, we must ignore the tau vertex, otherwise ptFinalCalib computed at calo-only step and precision step will differ
91 const xAOD::Vertex* vertex = nullptr;
92 if (m_doVertexCorrection) vertex = xTau.vertex();
93
94 // loop over tau clusters
95 std::vector<xAOD::CaloVertexedTopoCluster> vertexedClusters = xTau.vertexedClusters();
96
97 for (const xAOD::CaloVertexedTopoCluster& vertexedCluster : vertexedClusters) {
98 const xAOD::CaloCluster& cluster = vertexedCluster.clust();
99
100 TLorentzVector clusterP4 = m_doVertexCorrection? vertexedCluster.p4() : cluster.p4();
101 if (clusterP4.DeltaR(tauAxis) > 0.2) continue;
102
103 // FIXME: should we use calE for EMTopo clusters ?
104 // what's the energy scale when calculating the cluster momentum
105 clE = cluster.calE();
106 Etot += clE;
107
108 if (clE > lead_cluster_frac) {
109 third_cluster_frac = second_cluster_frac;
110 second_cluster_frac = lead_cluster_frac;
111 lead_cluster_frac = clE;
112 }
113 else if (clE > second_cluster_frac) {
114 third_cluster_frac = second_cluster_frac;
115 second_cluster_frac = clE;
116 }
117 else if (clE > third_cluster_frac) {
118 third_cluster_frac = clE;
119 }
120
122 mean_center_lambda += clE*center_lambda;
123 else ATH_MSG_WARNING("Failed to retrieve moment: CENTER_LAMBDA");
124
126 mean_first_eng_dens += clE*first_eng_dens;
127 else ATH_MSG_WARNING("Failed to retrieve moment: FIRST_ENG_DENS");
128
130 mean_em_probability += clE*em_probability;
131
132 // FIXME: should we use calE for EMTopo clusters ?
133 // what's the energy scale when calculating the cluster momentum
134 if (em_probability>0.5) clusters_EM_P4 += cluster.p4(xAOD::CaloCluster::State::CALIBRATED);
135 else clusters_had_P4 += cluster.p4(xAOD::CaloCluster::State::CALIBRATED);
136 }
137 else ATH_MSG_WARNING("Failed to retrieve moment: EM_PROBABILITY");
138
140 mean_second_lambda += clE*second_lambda;
141 else ATH_MSG_WARNING("Failed to retrieve moment: SECOND_LAMBDA");
142
143 mean_presampler_frac += (cluster.eSample(CaloSampling::PreSamplerB) + cluster.eSample(CaloSampling::PreSamplerE));
144
145 // EM-scale equivalent of IntermediateAxis p4
146 if (vertex) {
147 xAOD::CaloVertexedTopoCluster vertexedClusterEM(cluster, xAOD::CaloCluster::State::UNCALIBRATED, vertex->position());
148 tauIntermediateAxisEM += vertexedClusterEM.p4();
149 }
150 else {
151 tauIntermediateAxisEM += cluster.p4(xAOD::CaloCluster::State::UNCALIBRATED);
152 }
153 }
154
155 // calculate mean values
156 if (Etot>0.) {
157 mean_center_lambda /= Etot;
158 mean_first_eng_dens /= Etot;
159 mean_em_probability /= Etot;
160 mean_second_lambda /= Etot;
161 mean_presampler_frac /= Etot;
162 lead_cluster_frac /= Etot;
163 second_cluster_frac /= Etot;
164 third_cluster_frac /= Etot;
165
166 if (mean_first_eng_dens>0.) mean_first_eng_dens = TMath::Log10(mean_first_eng_dens/Etot);
167 }
168
169 // cluster-based upsilon, ranges from -1 to 1
170 double upsilon_cluster = -2.;
171 if (clusters_had_P4.E()+clusters_EM_P4.E()!=0.)
172 upsilon_cluster = (clusters_had_P4.E()-clusters_EM_P4.E())/(clusters_had_P4.E()+clusters_EM_P4.E());
173
174 xTau.setDetail(xAOD::TauJetParameters::ClustersMeanCenterLambda, static_cast<float>(mean_center_lambda));
175 xTau.setDetail(xAOD::TauJetParameters::ClustersMeanFirstEngDens, static_cast<float>(mean_first_eng_dens));
176 xTau.setDetail(xAOD::TauJetParameters::ClustersMeanEMProbability, static_cast<float>(mean_em_probability));
177 xTau.setDetail(xAOD::TauJetParameters::ClustersMeanSecondLambda, static_cast<float>(mean_second_lambda));
178 xTau.setDetail(xAOD::TauJetParameters::ClustersMeanPresamplerFrac, static_cast<float>(mean_presampler_frac));
179
180 static const SG::Accessor<float> acc_ClusterTotalEnergy("ClusterTotalEnergy");
181 acc_ClusterTotalEnergy(xTau) = static_cast<float>(Etot);
182
183 static const SG::Accessor<float> acc_ptIntermediateAxisEM("ptIntermediateAxisEM");
184 acc_ptIntermediateAxisEM(xTau) = static_cast<float>(tauIntermediateAxisEM.Pt());
185
186 // online-specific, not defined in TauDefs enum
187 static const SG::Accessor<float> acc_LeadClusterFrac("LeadClusterFrac");
188 static const SG::Accessor<float> acc_UpsilonCluster("UpsilonCluster");
189 acc_LeadClusterFrac(xTau) = static_cast<float>(lead_cluster_frac);
190 acc_UpsilonCluster(xTau) = static_cast<float>(upsilon_cluster);
191
192 if (inTrigger()) {
193 // for now only used by trigger, but could be used by offline 0p in the 2022 reprocessing
194 static const SG::Accessor<float> acc_SecondClusterFrac("SecondClusterFrac");
195 static const SG::Accessor<float> acc_ThirdClusterFrac("ThirdClusterFrac");
196 acc_SecondClusterFrac(xTau) = static_cast<float>(second_cluster_frac);
197 acc_ThirdClusterFrac(xTau) = static_cast<float>(third_cluster_frac);
198
199 return StatusCode::SUCCESS;
200 }
201
202 // summing corrected Pi0 PFO energies
203 TLorentzVector Pi0_totalP4;
204 Pi0_totalP4.SetPtEtaPhiM(0,0,0,0);
205
206 for (size_t i=0; i<xTau.nPi0PFOs(); i++){
207 Pi0_totalP4 += xTau.pi0PFO(i)->p4();
208 }
209
210 double Pi0_totalE = Pi0_totalP4.E();
211
212 // summing tau track momenta
213 TLorentzVector charged_totalP4;
214 charged_totalP4.SetPtEtaPhiM(0,0,0,0);
215
216 for (const xAOD::TauTrack* track : xTau.tracks()) {
217 charged_totalP4 += track->p4();
218 }
219
220 double charged_totalE = charged_totalP4.E();
221
222 // calculate relative difference and decorate onto tau
223 double relDiff=0.;
224 if (Pi0_totalE+charged_totalE != 0.){
225 relDiff = (charged_totalE - Pi0_totalE) / (charged_totalE + Pi0_totalE) ;
226 }
227 xTau.setDetail(xAOD::TauJetParameters::PFOEngRelDiff, static_cast<float>(relDiff));
228
229 return StatusCode::SUCCESS;
230}
#define ATH_MSG_WARNING(x)
SG::ReadHandleKey< xAOD::VertexContainer > m_vertexContainerKey
SG::ReadDecorHandleKey< xAOD::EventInfo > m_aveIntPerXKey
SG::ReadHandleKey< xAOD::EventShape > m_eventShapeKey
Gaudi::Property< bool > m_doVertexCorrection
bool inTrigger() const
bool retrieveMoment(MomentType type, double &value) const
Retrieve individual moment.
float eSample(const CaloSample sampling) const
@ SECOND_LAMBDA
Second Moment in .
@ EM_PROBABILITY
Classification probability to be em-like.
@ FIRST_ENG_DENS
First Moment in E/V.
@ CENTER_LAMBDA
Shower depth at Cluster Centroid.
flt_t calE() const
Geet Energy in signal state CALIBRATED.
virtual FourMom_t p4() const
The full 4-momentum of the particle.
TLorentzVector getTauAxis(const xAOD::TauJet &tau, bool doVertexCorrection=true)
Return the four momentum of the tau axis The tau axis is widely used to select clusters and cells in ...
std::vector< xAOD::CaloVertexedTopoCluster > vertexedClusters(const xAOD::TauJet &tau, double dRMax)
@ PileUp
Pile-up vertex.
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
VertexContainer_v1 VertexContainer
Definition of the current "Vertex container version".
Vertex_v1 Vertex
Define the latest version of the vertex class.
TauTrack_v1 TauTrack
Definition of the current version.
Definition TauTrack.h:16
#define CXXUTILS_TRAPPING_FP
Definition trapping_fp.h:24

◆ executePanTau()

StatusCode TauRecToolBase::executePanTau ( xAOD::TauJet & pTau,
xAOD::ParticleContainer & particleContainer,
xAOD::PFOContainer & neutralPFOContainer ) const
overridevirtualinherited

Implements ITauToolBase.

Reimplemented in PanTau::PanTauProcessor.

Definition at line 197 of file TauRecToolBase.cxx.

197 {
198 ATH_MSG_ERROR("function not implemented");
199 return StatusCode::FAILURE;
200}
#define ATH_MSG_ERROR(x)

◆ executePi0ClusterCreator()

StatusCode TauRecToolBase::executePi0ClusterCreator ( xAOD::TauJet & pTau,
xAOD::PFOContainer & neutralPFOContainer,
xAOD::PFOContainer & hadronicPFOContainer,
const xAOD::CaloClusterContainer & pCaloClusterContainer ) const
overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauPi0ClusterCreator.

Definition at line 175 of file TauRecToolBase.cxx.

177 {
178 ATH_MSG_ERROR("function not implemented");
179 return StatusCode::FAILURE;
180}

◆ executePi0ClusterScaler()

StatusCode TauRecToolBase::executePi0ClusterScaler ( xAOD::TauJet & pTau,
xAOD::PFOContainer & neutralPFOContainer,
xAOD::PFOContainer & chargedPFOContainer ) const
overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauPi0ClusterScaler.

Definition at line 187 of file TauRecToolBase.cxx.

187 {
188 ATH_MSG_ERROR("function not implemented");
189 return StatusCode::FAILURE;
190}

◆ executePi0CreateROI()

StatusCode TauRecToolBase::executePi0CreateROI ( xAOD::TauJet & pTau,
CaloConstCellContainer & caloCellContainer,
boost::dynamic_bitset<> & map ) const
overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauPi0CreateROI.

Definition at line 149 of file TauRecToolBase.cxx.

149 {
150 ATH_MSG_ERROR("function not implemented");
151 return StatusCode::FAILURE;
152}

◆ executePi0nPFO()

StatusCode TauRecToolBase::executePi0nPFO ( xAOD::TauJet & pTau,
xAOD::PFOContainer & neutralPFOContainer ) const
overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauPi0ScoreCalculator, and TauPi0Selector.

Definition at line 192 of file TauRecToolBase.cxx.

192 {
193 ATH_MSG_ERROR("function not implemented");
194 return StatusCode::FAILURE;
195}

◆ executeShotFinder()

StatusCode TauRecToolBase::executeShotFinder ( xAOD::TauJet & pTau,
xAOD::CaloClusterContainer & shotClusterContainer,
xAOD::PFOContainer & PFOContainer ) const
overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauShotFinder.

Definition at line 170 of file TauRecToolBase.cxx.

170 {
171 ATH_MSG_ERROR("function not implemented");
172 return StatusCode::FAILURE;
173}

◆ executeTrackClassifier()

StatusCode TauRecToolBase::executeTrackClassifier ( xAOD::TauJet & pTau,
xAOD::TauTrackContainer & tauTrackContainer ) const
overridevirtualinherited

Implements ITauToolBase.

Reimplemented in tauRecTools::TauTrackRNNClassifier.

Definition at line 165 of file TauRecToolBase.cxx.

165 {
166 ATH_MSG_ERROR("function not implemented");
167 return StatusCode::FAILURE;
168}

◆ executeTrackFinder()

StatusCode TauRecToolBase::executeTrackFinder ( xAOD::TauJet & pTau,
xAOD::TauTrackContainer & tauTrackContainer ) const
overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauTrackFinder.

Definition at line 160 of file TauRecToolBase.cxx.

160 {
161 ATH_MSG_ERROR("function not implemented");
162 return StatusCode::FAILURE;
163}

◆ executeVertexFinder()

StatusCode TauRecToolBase::executeVertexFinder ( xAOD::TauJet & pTau,
const xAOD::VertexContainer * vertexContainer = nullptr ) const
overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauVertexFinder.

Definition at line 155 of file TauRecToolBase.cxx.

155 {
156 ATH_MSG_ERROR("function not implemented");
157 return StatusCode::FAILURE;
158}

◆ executeVertexVariables()

StatusCode TauRecToolBase::executeVertexVariables ( xAOD::TauJet & pTau,
xAOD::VertexContainer & vertexContainer ) const
overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauVertexVariables.

Definition at line 182 of file TauRecToolBase.cxx.

182 {
183 ATH_MSG_ERROR("function not implemented");
184 return StatusCode::FAILURE;
185}

◆ 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

◆ finalize()

StatusCode TauRecToolBase::finalize ( )
overridevirtualinherited

Finalizer.

Implements ITauToolBase.

Reimplemented in PanTau::PanTauProcessor.

Definition at line 206 of file TauRecToolBase.cxx.

206 {
207 return StatusCode::SUCCESS;
208}

◆ find_file()

std::string TauRecToolBase::find_file ( const std::string & fname) const
inherited

Definition at line 19 of file TauRecToolBase.cxx.

19 {
20 std::string full_path;
21 //offline calib files are in GroupData
22 //online calib files are in release
23 full_path = PathResolverFindCalibFile(m_tauRecToolsTag+"/"+fname);
24 if(full_path.empty()) full_path = PathResolverFindCalibFile(fname);
25 return full_path;
26}
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
Gaudi::Property< std::string > m_tauRecToolsTag

◆ 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 }
ServiceHandle< StoreGateSvc > & evtStore()

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

◆ inAOD()

bool TauRecToolBase::inAOD ( ) const
inlineprotectedinherited

Definition at line 88 of file TauRecToolBase.h.

88{ return m_in_AOD; }
Gaudi::Property< bool > m_in_AOD

◆ inEleRM()

bool TauRecToolBase::inEleRM ( ) const
inlineprotectedinherited

Definition at line 89 of file TauRecToolBase.h.

89{ return m_in_EleRM; }
Gaudi::Property< bool > m_in_EleRM

◆ initialize()

StatusCode MvaTESVariableDecorator::initialize ( void )
overridevirtual

Tool initializer.

Reimplemented from TauRecToolBase.

Definition at line 18 of file MvaTESVariableDecorator.cxx.

18 {
19
20 ATH_CHECK(m_aveIntPerXKey.initialize());
23
24 return StatusCode::SUCCESS;
25}
#define ATH_CHECK
Evaluate an expression and check for errors.

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

◆ inTrigger()

bool TauRecToolBase::inTrigger ( ) const
inlineprotectedinherited

Definition at line 87 of file TauRecToolBase.h.

87{ return m_in_trigger; }
Gaudi::Property< bool > m_in_trigger

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

◆ print()

◆ readConfig()

StatusCode TauRecToolBase::readConfig ( )
overridevirtualinherited

Implements ITauToolBase.

Definition at line 27 of file TauRecToolBase.cxx.

27 {
28 // Sanity check to see if property ConfigPath is declared for a tool. Might be
29 // removed once all tools are updated to have a config path declared.
30 // in athena getProperties returns std::vector<Gaudi::Details::PropertyBase*>
31 // in rc getProperties returns std::map<std::string,Property*>
32#ifndef XAOD_STANDALONE
33 bool configPathDeclared = false;
34 for (Gaudi::Details::PropertyBase* property : getProperties())
35 {
36 if (property->name() == "ConfigPath")
37 {
38 configPathDeclared = true;
39 break;
40 }
41 }
42 if (!configPathDeclared)
43#elif defined(XAOD_STANDALONE)
44 PropertyMgr::PropMap_t property_map = getPropertyMgr()->getProperties();
45 if (property_map.find("ConfigPath") == property_map.end())
46#else
47# error "What environment are we in?!?"
48#endif // XAOD_STANDALONE
49 {
50 ATH_MSG_INFO("No config file path property declared yet, this is not recommended");
51 return StatusCode::SUCCESS;
52 }
53
54 // get configured config path and load file via TEnv
55 const std::string* config_file_path_property;
56 // if (getProperty("ConfigPath", config_file_path).isFailure())
57 // return StatusCode::FAILURE;
58 config_file_path_property = getProperty<std::string>("ConfigPath");
59 std::string config_file_path = find_file(*config_file_path_property);
60 TEnv env;
61 env.ReadFile(PathResolverFindCalibFile(config_file_path).c_str(),kEnvAll);
62
63 THashList* lList = env.GetTable();
64 for( Int_t i = 0; lList && i < lList->GetEntries(); ++i )
65 {
67 // types of properties are handled differently as well
68#ifndef XAOD_STANDALONE
69 // get type of variable with the entry name
70 const std::type_info* type = getProperty(lList->At( i )->GetName()).type_info();
71
72 // search for type is needed by env.GetValue function (needs a variable of the correct type as 2nd argument)
73 if (*type == typeid(bool))
74 sc = this->setProperty(lList->At( i )->GetName(),
75 bool(env.GetValue(lList->At( i )->GetName(),bool(true))));
76 else if (*type == typeid(int))
77 sc = this->setProperty(lList->At( i )->GetName(),
78 env.GetValue(lList->At( i )->GetName(),int(0)));
79 else if (*type == typeid(float))
80 sc = this->setProperty(lList->At( i )->GetName(),
81 env.GetValue(lList->At( i )->GetName(),float(0)));
82 else if (*type == typeid(double))
83 sc = this->setProperty(lList->At( i )->GetName(),
84 env.GetValue(lList->At( i )->GetName(),double(0)));
85 else if (*type == typeid(std::string))
86 sc = this->setProperty(lList->At( i )->GetName(),
87 env.GetValue(lList->At( i )->GetName(),""));
88#else
89 // get type of variable with the entry name
90 Property::Type type = getPropertyMgr()->getProperty(lList->At( i )->GetName())->type();
91
92 if (type == Property::BOOL)
93 sc = this->setProperty(lList->At( i )->GetName(),
94 bool(env.GetValue(lList->At( i )->GetName(),bool(true))));
95 else if (type == Property::INT)
96 sc = this->setProperty(lList->At( i )->GetName(),
97 env.GetValue(lList->At( i )->GetName(),int(0)));
98 else if (type == Property::FLOAT)
99 sc = this->setProperty(lList->At( i )->GetName(),
100 env.GetValue(lList->At( i )->GetName(),float(0)));
101 else if (type == Property::DOUBLE)
102 sc = this->setProperty(lList->At( i )->GetName(),
103 env.GetValue(lList->At( i )->GetName(),double(0)));
104 else if (type == Property::STRING)
105 sc = this->setProperty(lList->At( i )->GetName(),
106 env.GetValue(lList->At( i )->GetName(),""));
107#endif // XAOD_STANDALONE
108 else
109 {
110#ifndef XAOD_STANDALONE
111 ATH_MSG_FATAL("there was a problem to find the correct type enum: "<<type->name());
112#else
113 ATH_MSG_FATAL("there was a problem to find the correct type enum: "<<type);
114#endif // XAOD_STANDALONE
115 return StatusCode::FAILURE;
116 }
117 if (!sc.isSuccess()) {
118 ATH_MSG_FATAL("failed to set property: " << lList->At( i )->GetName());
119 return StatusCode::FAILURE;
120 }
121 }
122 return StatusCode::SUCCESS;
123}
#define ATH_MSG_FATAL(x)
static Double_t sc
void setProperty(columnar::PythonToolHandle &self, const std::string &key, nb::object value)
Type
Property type enumeration.
Definition Property.h:27
std::string find_file(const std::string &fname) const
const T * getProperty(const std::string &name) const
Get one of the tool's properties.
::StatusCode StatusCode
StatusCode definition for legacy code.

◆ 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 }
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< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

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

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

SG::ReadDecorHandleKey<xAOD::EventInfo> MvaTESVariableDecorator::m_aveIntPerXKey
private
Initial value:
{this,
"averageInteractionsPerCrossingKey",
"EventInfo.averageInteractionsPerCrossing",
"Decoration for Average Interaction Per Crossing"}

Definition at line 34 of file MvaTESVariableDecorator.h.

34 {this,
35 "averageInteractionsPerCrossingKey",
36 "EventInfo.averageInteractionsPerCrossing",
37 "Decoration for Average Interaction Per Crossing"};

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

Gaudi::Property<bool> MvaTESVariableDecorator::m_doVertexCorrection {this, "VertexCorrection", true}
private

Definition at line 32 of file MvaTESVariableDecorator.h.

32{this, "VertexCorrection", true};

◆ m_eventShapeKey

SG::ReadHandleKey<xAOD::EventShape> MvaTESVariableDecorator::m_eventShapeKey
private
Initial value:
{this,
"EventShapeKey",
"Kt4LCTopoOriginEventShape",
"EventShape key"}

Definition at line 44 of file MvaTESVariableDecorator.h.

44 {this,
45 "EventShapeKey",
46 "Kt4LCTopoOriginEventShape",
47 "EventShape key"};

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

Gaudi::Property<bool> TauRecToolBase::m_in_AOD {this, "inAOD", false, "Indicate if the tool is running on AOD"}
protectedinherited

Definition at line 77 of file TauRecToolBase.h.

77{this, "inAOD", false, "Indicate if the tool is running on AOD"};

◆ m_in_EleRM

Gaudi::Property<bool> TauRecToolBase::m_in_EleRM {this, "inEleRM", false, "Indicate if the tool is running on EleRM routine"}
protectedinherited

Definition at line 78 of file TauRecToolBase.h.

78{this, "inEleRM", false, "Indicate if the tool is running on EleRM routine"};

◆ m_in_trigger

Gaudi::Property<bool> TauRecToolBase::m_in_trigger {this, "inTrigger", false, "Indicate if the tool is running on trigger"}
protectedinherited

Definition at line 76 of file TauRecToolBase.h.

76{this, "inTrigger", false, "Indicate if the tool is running on trigger"};

◆ m_tauRecToolsTag

Gaudi::Property<std::string> TauRecToolBase::m_tauRecToolsTag {this, "calibFolder", "tauRecTools/R22_preprod", "CVMFS path to the tau calibration folder"}
protectedinherited

Definition at line 79 of file TauRecToolBase.h.

79{this, "calibFolder", "tauRecTools/R22_preprod", "CVMFS path to the tau calibration folder"};

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vertexContainerKey

SG::ReadHandleKey<xAOD::VertexContainer> MvaTESVariableDecorator::m_vertexContainerKey
private
Initial value:
{this,
"Key_vertexInputContainer",
"PrimaryVertices",
"input vertex container key"}

Definition at line 39 of file MvaTESVariableDecorator.h.

39 {this,
40 "Key_vertexInputContainer",
41 "PrimaryVertices",
42 "input vertex container key"};

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.


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