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

Creates Pi0 clusters (Pi0 Finder). More...

#include <TauPi0ClusterCreator.h>

Inheritance diagram for TauPi0ClusterCreator:
Collaboration diagram for TauPi0ClusterCreator:

Public Member Functions

 TauPi0ClusterCreator (const std::string &name)
 
virtual ~TauPi0ClusterCreator ()=default
 
virtual StatusCode executePi0ClusterCreator (xAOD::TauJet &pTau, xAOD::PFOContainer &neutralPFOContainer, xAOD::PFOContainer &hadronicClusterPFOContainer, const xAOD::CaloClusterContainer &pi0CaloClusContainer) const override
 
virtual StatusCode initialize () override
 Tool initializer. More...
 
virtual StatusCode eventInitialize () override
 Event initializer - called at the beginning of each event. More...
 
virtual StatusCode execute (xAOD::TauJet &pTau) const override
 Execute - called for each tau candidate. More...
 
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 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. More...
 
virtual StatusCode finalize () override
 Finalizer. More...
 
std::string find_file (const std::string &fname) const
 
virtual StatusCode readConfig () override
 
virtual void print () const
 Print the state of the tool. More...
 
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & evtStore () const
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
 
virtual StatusCode sysInitialize () override
 Perform system initialization for an algorithm. More...
 
virtual StatusCode sysStart () override
 Handle START transition. More...
 
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles. More...
 
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles. More...
 
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T > &t)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc="none")
 Declare a new Gaudi property. More...
 
void updateVHKA (Gaudi::Details::PropertyBase &)
 
MsgStream & msg () const
 
MsgStream & msg (const MSG::Level lvl) const
 
bool msgLvl (const MSG::Level lvl) const
 

Protected Member Functions

bool inTrigger () const
 
bool inAOD () const
 
bool inEleRM () const
 
void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
 
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed. More...
 

Protected Attributes

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< StoreGateSvc > StoreGateSvc_t
 

Private Member Functions

StatusCode configureNeutralPFO (const xAOD::CaloCluster &cluster, const xAOD::CaloClusterContainer &pi0ClusterContainer, const xAOD::TauJet &tau, const std::vector< const xAOD::PFO * > &shotPFOs, const std::map< unsigned, const xAOD::CaloCluster * > &shotsInCluster, xAOD::PFO &neutralPFO) const
 Configure the neutral PFO. More...
 
StatusCode configureHadronicPFO (const xAOD::CaloVertexedTopoCluster &cluster, double clusterEnergyHad, xAOD::PFO &hadronicPFO) const
 Configure the haronic PFO. More...
 
std::map< unsigned, const xAOD::CaloCluster * > getShotToClusterMap (const std::vector< const xAOD::PFO * > &shotVector, const xAOD::CaloClusterContainer &pi0ClusterContainer, const xAOD::TauJet &pTau) const
 
std::vector< unsigned > getShotsMatchedToCluster (const std::vector< const xAOD::PFO * > &shotVector, const std::map< unsigned, const xAOD::CaloCluster * > &clusterToShotMap, const xAOD::CaloCluster &pi0Cluster) const
 
int getNPhotons (const std::vector< const xAOD::PFO * > &shotVector, const std::vector< unsigned > &shotsInCluster) const
 
float getEM1CoreFrac (const xAOD::CaloCluster &cluster) const
 fraction of cluster enegry in central EM1 cells More...
 
std::vector< int > getNPosECells (const xAOD::CaloCluster &cluster) const
 number of cells from cluster with positive energy in PS, EM1 and EM2 More...
 
std::vector< float > get1stEtaMomWRTCluster (const xAOD::CaloCluster &cluster) const
 first eta moment in PS, EM1 and EM2 w.r.t cluster eta More...
 
std::vector< float > get2ndEtaMomWRTCluster (const xAOD::CaloCluster &cluster) const
 second eta moment in PS, EM1 and EM2 w.r.t cluster eta More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

Gaudi::Property< double > m_clusterEtCut {this, "ClusterEtCut", 0.5 * Gaudi::Units::GeV, "Et threshould for pi0 candidate clusters"}
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Creates Pi0 clusters (Pi0 Finder).

Author
Will Davey will..nosp@m.dave.nosp@m.y@cer.nosp@m.n.ch
Benedict Winter bened.nosp@m.ict..nosp@m.tobia.nosp@m.s.wi.nosp@m.nter@.nosp@m.cern.nosp@m..ch
Stephanie Yuen steph.nosp@m.anie.nosp@m..yuen.nosp@m.@cer.nosp@m.n.ch

Definition at line 28 of file TauPi0ClusterCreator.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

◆ TauPi0ClusterCreator()

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

Definition at line 15 of file TauPi0ClusterCreator.cxx.

15  :
16  TauRecToolBase(name) {
17 }

◆ ~TauPi0ClusterCreator()

virtual TauPi0ClusterCreator::~TauPi0ClusterCreator ( )
virtualdefault

Member Function Documentation

◆ configureHadronicPFO()

StatusCode TauPi0ClusterCreator::configureHadronicPFO ( const xAOD::CaloVertexedTopoCluster &  cluster,
double  clusterEnergyHad,
xAOD::PFO &  hadronicPFO 
) const
private

Configure the haronic PFO.

Definition at line 479 of file TauPi0ClusterCreator.cxx.

481  {
482  double clusterPtHad = clusterEnergyHad/std::cosh(cluster.eta());
483  hadronicPFO.setP4(clusterPtHad, cluster.eta(), cluster.phi(), 0.);
484 
485  return StatusCode::SUCCESS;
486 }

◆ configureNeutralPFO()

StatusCode TauPi0ClusterCreator::configureNeutralPFO ( const xAOD::CaloCluster &  cluster,
const xAOD::CaloClusterContainer &  pi0ClusterContainer,
const xAOD::TauJet &  tau,
const std::vector< const xAOD::PFO * > &  shotPFOs,
const std::map< unsigned, const xAOD::CaloCluster * > &  shotsInCluster,
xAOD::PFO &  neutralPFO 
) const
private

Configure the neutral PFO.

Definition at line 378 of file TauPi0ClusterCreator.cxx.

383  {
384  // Set the properties of the PFO
385  // -- Default value
386  neutralPFO.setBDTPi0Score(-9999.);
387  neutralPFO.setCharge(0);
388  neutralPFO.setAttribute<int>(xAOD::PFODetails::PFOAttributes::nPi0Proto, -1);
389 
390  // -- CENTER_MAG
391  double CENTER_MAG = 0.0;
392  if (!cluster.retrieveMoment(xAOD::CaloCluster::MomentType::CENTER_MAG, CENTER_MAG)) {
393  ATH_MSG_WARNING("Couldn't retrieve CENTER_MAG moment. Set it to 0.");
394  }
395  neutralPFO.setCenterMag( (float) CENTER_MAG);
396 
397  // -- Number of photons
398  std::vector<unsigned> shotsInCluster = getShotsMatchedToCluster(shotPFOs, shotToClusterMap, cluster);
399  int NHitsInEM1 = getNPhotons(shotPFOs, shotsInCluster);
400  neutralPFO.setAttribute<int>(xAOD::PFODetails::PFOAttributes::cellBased_NHitsInEM1, NHitsInEM1);
401 
402  // -- Energy at each layer
403  float eEM1 = cluster.eSample(CaloSampling::EMB1) + cluster.eSample(CaloSampling::EME1);
404  neutralPFO.setAttribute<float>(xAOD::PFODetails::PFOAttributes::cellBased_energy_EM1, eEM1);
405 
406  float eEM2 = cluster.eSample(CaloSampling::EMB2) + cluster.eSample(CaloSampling::EME2);
407  neutralPFO.setAttribute<float>(xAOD::PFODetails::PFOAttributes::cellBased_energy_EM2, eEM2);
408 
409  // -- Number of positive cells in each layer
410  std::vector<int> nPosECells = getNPosECells(cluster);
411  neutralPFO.setAttribute<int>(xAOD::PFODetails::PFOAttributes::cellBased_NPosECells_PS, nPosECells.at(0));
412  neutralPFO.setAttribute<int>(xAOD::PFODetails::PFOAttributes::cellBased_NPosECells_EM1, nPosECells.at(1));
413  neutralPFO.setAttribute<int>(xAOD::PFODetails::PFOAttributes::cellBased_NPosECells_EM2, nPosECells.at(2));
414 
415  // -- Core Fraction of the energy in EM1
416  float EM1CoreFrac = getEM1CoreFrac(cluster);
417  neutralPFO.setAttribute<float>(xAOD::PFODetails::PFOAttributes::cellBased_EM1CoreFrac, EM1CoreFrac);
418 
419  // -- First moment of deltaEta(cluster, cell) in EM1 and EM2
420  std::vector<float> deltaEtaFirstMom = get1stEtaMomWRTCluster(cluster);
421  neutralPFO.setAttribute<float>(xAOD::PFODetails::PFOAttributes::cellBased_firstEtaWRTClusterPosition_EM1, deltaEtaFirstMom.at(1));
422  neutralPFO.setAttribute<float>(xAOD::PFODetails::PFOAttributes::cellBased_firstEtaWRTClusterPosition_EM2, deltaEtaFirstMom.at(2));
423 
424  // -- Second moment of deltaEta(cluster, cell) in EM1 and EM2
425  std::vector<float> secondEtaWRTClusterPositionInLayer = get2ndEtaMomWRTCluster(cluster);
426  neutralPFO.setAttribute<float>(xAOD::PFODetails::PFOAttributes::cellBased_secondEtaWRTClusterPosition_EM1, secondEtaWRTClusterPositionInLayer.at(1));
427  neutralPFO.setAttribute<float>(xAOD::PFODetails::PFOAttributes::cellBased_secondEtaWRTClusterPosition_EM2, secondEtaWRTClusterPositionInLayer.at(2));
428 
429  // -- Retrieve cluster moments
430  using Moment = xAOD::CaloCluster::MomentType;
431  using Attribute = xAOD::PFODetails::PFOAttributes;
432  const std::array< std::pair<Moment, Attribute>, 12> momentAttributePairs {{
433  {Moment::FIRST_ETA, Attribute::cellBased_FIRST_ETA},
434  {Moment::SECOND_R, Attribute::cellBased_SECOND_R},
435  {Moment::SECOND_LAMBDA, Attribute::cellBased_SECOND_LAMBDA},
436  {Moment::DELTA_PHI, Attribute::cellBased_DELTA_PHI},
437  {Moment::DELTA_THETA, Attribute::cellBased_DELTA_THETA},
438  {Moment::CENTER_LAMBDA, Attribute::cellBased_CENTER_LAMBDA},
439  {Moment::LATERAL, Attribute::cellBased_LATERAL},
440  {Moment::LONGITUDINAL, Attribute::cellBased_LONGITUDINAL},
441  {Moment::ENG_FRAC_EM, Attribute::cellBased_ENG_FRAC_EM},
442  {Moment::ENG_FRAC_MAX, Attribute::cellBased_ENG_FRAC_MAX},
443  {Moment::ENG_FRAC_CORE, Attribute::cellBased_ENG_FRAC_CORE},
444  {Moment::SECOND_ENG_DENS, Attribute::cellBased_SECOND_ENG_DENS}
445  }};
446 
447  for (const auto& [moment, attribute] : momentAttributePairs) {
448  double value = 0.0;
449  if (! cluster.retrieveMoment(moment, value)) {
450  ATH_MSG_WARNING("Cound not retrieve " << moment);
451  }
452  neutralPFO.setAttribute(attribute, static_cast<float>(value));
453  }
454 
455  // -- Element link to the cluster
456  ElementLink<xAOD::CaloClusterContainer> clusElementLink;
457  clusElementLink.toContainedElement(pi0ClusterContainer, &cluster);
458  neutralPFO.setClusterLink( clusElementLink );
459 
460  // -- Element link to the shots
461  std::vector<ElementLink<xAOD::IParticleContainer>> shotlinks;
462  for (unsigned index = 0; index < shotsInCluster.size(); ++index) {
463  auto& shotPFOElementLink = tau.shotPFOLinks().at(shotsInCluster.at(index));
464  ElementLink<xAOD::IParticleContainer> shotElementLink( shotPFOElementLink.dataID(), shotPFOElementLink.index() );
465  if (!shotElementLink.isValid()) {
466  ATH_MSG_WARNING("Created an invalid element link to xAOD::PFO");
467  }
468  shotlinks.push_back(shotElementLink);
469  }
470  if(!neutralPFO.setAssociatedParticleLinks( xAOD::PFODetails::TauShot,shotlinks)) {
471  ATH_MSG_WARNING("Couldn't add shot links to neutral PFO!");
472  }
473 
474  return StatusCode::SUCCESS;
475 }

◆ declareGaudiProperty() [1/4]

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

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

Definition at line 170 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [2/4]

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

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [3/4]

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

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

Definition at line 184 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [4/4]

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

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

Definition at line 199 of file AthCommonDataStore.h.

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

◆ declareProperty() [1/6]

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

Declare a new Gaudi property.

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

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

Definition at line 245 of file AthCommonDataStore.h.

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

◆ declareProperty() [2/6]

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

Declare a new Gaudi property.

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

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

Definition at line 221 of file AthCommonDataStore.h.

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

◆ declareProperty() [3/6]

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

Definition at line 259 of file AthCommonDataStore.h.

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

◆ declareProperty() [4/6]

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

Declare a new Gaudi property.

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

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

Definition at line 333 of file AthCommonDataStore.h.

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

◆ declareProperty() [5/6]

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

Declare a new Gaudi property.

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

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

Definition at line 352 of file AthCommonDataStore.h.

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

◆ declareProperty() [6/6]

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

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ 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() [1/2]

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

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

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ execute()

StatusCode TauRecToolBase::execute ( xAOD::TauJet &  pTau) const
overridevirtualinherited

Execute - called for each tau candidate.

Implements ITauToolBase.

Reimplemented in TauDecayModeNNClassifier, MvaTESVariableDecorator, TauCombinedTES, MvaTESEvaluator, TauWPDecorator, TauAxisSetter, TauEleOverlapChecker, TauJetRNNEvaluator, TauClusterFinder, TauCalibrateLC, TauGNNEvaluator, TauSubstructureVariables, TauVertexedClusterDecorator, TauAODSelector, TauElectronVetoVariables, JetSeedBuilder, TauCellVariables, TauCommonCalcVars, TauIDVarCalculator, and TauAODMuonRemovalTool.

Definition at line 138 of file TauRecToolBase.cxx.

138  {
139  ATH_MSG_ERROR("function not implemented");
140  return StatusCode::FAILURE;
141 }

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

◆ executePi0ClusterCreator()

StatusCode TauPi0ClusterCreator::executePi0ClusterCreator ( xAOD::TauJet &  pTau,
xAOD::PFOContainer &  neutralPFOContainer,
xAOD::PFOContainer &  hadronicClusterPFOContainer,
const xAOD::CaloClusterContainer &  pi0CaloClusContainer 
) const
overridevirtual

Reimplemented from TauRecToolBase.

Definition at line 21 of file TauPi0ClusterCreator.cxx.

24  {
25  // Any tau needs to have PFO vectors. Set empty vectors before nTrack cut
26  std::vector<ElementLink<xAOD::PFOContainer>> empty;
27  tau.setProtoNeutralPFOLinks(empty);
28  tau.setHadronicPFOLinks(empty);
29 
30  // only need pi0s and shots for 0-5 prong taus
31  if (!tauRecTools::doPi0andShots(tau)) {
32  return StatusCode::SUCCESS;
33  }
34 
35  // Retrieve Ecal1 shots and match them to clusters
36  std::vector<const xAOD::PFO*> shotPFOs;
37 
38  unsigned nShots = tau.nShotPFOs();
39  for (unsigned index=0; index<nShots; ++index) {
40  const xAOD::PFO* shotPFO = tau.shotPFO(index);
41  shotPFOs.push_back(shotPFO);
42  }
43 
44  // Map shot to the pi0 cluster
45  std::map<unsigned, const xAOD::CaloCluster*> shotToClusterMap = getShotToClusterMap(shotPFOs, pi0ClusterContainer, tau);
46 
47  // We will always perform the vertex correction
48  const xAOD::Vertex* vertex = tau.vertex();
49 
50  // Tau custom PFO reconstruction is only used in offline reconstrution
51  TLorentzVector tauAxis = tauRecTools::getTauAxis(tau);
52 
53  // Loop over custom pi0 clusters, and create neutral PFOs
54  for (const xAOD::CaloCluster* cluster: pi0ClusterContainer) {
55  // correct pi0 clusters w.r.t. the tau vertex
56  TLorentzVector clusterP4;
57  if (vertex) {
58  xAOD::CaloVertexedTopoCluster vertexedCluster(*cluster, vertex->position());
59  clusterP4 = vertexedCluster.p4();
60  }
61  else {
62  clusterP4 = cluster->p4();
63  }
64 
65  // Clusters must have enough energy, and within 0.4 cone of the tau candidate
66  if (clusterP4.Pt() < m_clusterEtCut) continue;
67  if (clusterP4.DeltaR(tauAxis) > 0.4) continue;
68 
69  // Create the neutral PFOs
70  xAOD::PFO* neutralPFO = new xAOD::PFO();
71  neutralPFOContainer.push_back(neutralPFO);
72 
73  // -- Set the PFO four momentum corrected w.r.t. the tau vertex
74  neutralPFO->setP4(clusterP4.Pt(), clusterP4.Eta(), clusterP4.Phi(), 0.);
75 
76  // Add the link to the tau candidate
77  ElementLink<xAOD::PFOContainer> PFOElementLink;
78  PFOElementLink.toContainedElement(neutralPFOContainer, neutralPFO);
79  tau.addProtoNeutralPFOLink(PFOElementLink);
80 
81  ATH_CHECK(configureNeutralPFO(*cluster, pi0ClusterContainer, tau, shotPFOs, shotToClusterMap, *neutralPFO));
82  }
83 
84  // Loop over clusters from jet seed, and create hadronic PFOs
85  std::vector<xAOD::CaloVertexedTopoCluster> vertexedClusterList = tau.vertexedClusters();
86  for (const xAOD::CaloVertexedTopoCluster& vertexedCluster : vertexedClusterList){
87  TLorentzVector clusterP4 = vertexedCluster.p4();
88 
89  // Clusters must have positive energy, and within 0.2 cone of the tau candidate
90  if(clusterP4.E()<=0.) continue;
91  if(clusterP4.DeltaR(tauAxis) > 0.2) continue;
92 
93  double clusterEnergyHad = 0.;
94 
95  const xAOD::CaloCluster& cluster = vertexedCluster.clust();
96  const CaloClusterCellLink* cellLinks = cluster.getCellLinks();
97  if (!cellLinks) {
98  ATH_MSG_WARNING("The cell links of the tau cluster is unavailable.");
99  continue;
100  }
101  CaloClusterCellLink::const_iterator cellLink = cellLinks->begin();
102  for (; cellLink != cellLinks->end(); ++cellLink) {
103  const CaloCell* cell = static_cast<const CaloCell*>(*cellLink);
104 
105  int sampling = cell->caloDDE()->getSampling();
106  if (sampling < 8) continue;
107 
108  // TODO: what is the weight for EMTopo
109  double cellEnergy = cell->e() * cellLink.weight();
110  clusterEnergyHad += cellEnergy;
111  }
112 
113  // Energy in Had Calorimeter must be positive
114  if(clusterEnergyHad <= 0.) continue;
115 
116  // Create the hadrnic PFO
117  xAOD::PFO* hadronicPFO = new xAOD::PFO();
118  hadronicPFOContainer.push_back(hadronicPFO);
119 
120  // Add element link from tau to hadronic PFO
121  ElementLink<xAOD::PFOContainer> PFOElementLink;
122  PFOElementLink.toContainedElement( hadronicPFOContainer, hadronicPFO );
123  tau.addHadronicPFOLink( PFOElementLink );
124 
125  ATH_CHECK(configureHadronicPFO(vertexedCluster, clusterEnergyHad, *hadronicPFO));
126  }
127 
128  return StatusCode::SUCCESS;
129 }

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

◆ get1stEtaMomWRTCluster()

std::vector< float > TauPi0ClusterCreator::get1stEtaMomWRTCluster ( const xAOD::CaloCluster &  cluster) const
private

first eta moment in PS, EM1 and EM2 w.r.t cluster eta

Definition at line 304 of file TauPi0ClusterCreator.cxx.

304  {
305  std::vector<float> deltaEtaFirstMom (3, 0.);
306  std::vector<float> totalEnergy (3, 0.);
307 
308  const CaloClusterCellLink* cellLinks = cluster.getCellLinks();
309  CaloClusterCellLink::const_iterator cellLink = cellLinks->begin();
310  for (; cellLink != cellLinks->end(); ++cellLink) {
311  const CaloCell* cell = static_cast<const CaloCell*>(*cellLink);
312 
313  // Only consider PS, EM1, and EM2
314  int sampling = cell->caloDDE()->getSampling();
315  int layer = sampling%4;
316  if (layer >= 3) continue;
317 
318  // Only consider positive cells
319  float cellEnergy = cell->e();
320  if (cellEnergy <= 0) continue;
321 
322  float deltaEta = cell->eta() - cluster.eta();
323  deltaEtaFirstMom[layer] += deltaEta * cellEnergy;
324  totalEnergy[layer] += cellEnergy;
325  }
326 
327  for (int layer=0; layer < 3; ++layer) {
328  if (totalEnergy[layer] != 0.) {
329  deltaEtaFirstMom[layer]/=std::abs(totalEnergy[layer]);
330  }
331  else {
332  deltaEtaFirstMom[layer]=0.;
333  }
334  }
335 
336  return deltaEtaFirstMom;
337 }

◆ get2ndEtaMomWRTCluster()

std::vector< float > TauPi0ClusterCreator::get2ndEtaMomWRTCluster ( const xAOD::CaloCluster &  cluster) const
private

second eta moment in PS, EM1 and EM2 w.r.t cluster eta

Definition at line 341 of file TauPi0ClusterCreator.cxx.

341  {
342  std::vector<float> deltaEtaSecondMom (3, 0.);
343  std::vector<float> totalEnergy (3, 0.);
344 
345  const CaloClusterCellLink* cellLinks = cluster.getCellLinks();
346  CaloClusterCellLink::const_iterator cellLink = cellLinks->begin();
347  for (; cellLink != cellLinks->end(); ++cellLink) {
348  const CaloCell* cell = static_cast<const CaloCell*>(*cellLink);
349 
350  // Only consider PS, EM1, and EM2
351  int sampling = cell->caloDDE()->getSampling();
352  int layer = sampling%4;
353  if (layer >= 3) continue;
354 
355  // Only consider positive cells
356  float cellEnergy=cell->e();
357  if (cellEnergy <= 0) continue;
358 
359  float deltaEta = cell->eta() - cluster.eta();
360  deltaEtaSecondMom[layer] += deltaEta * deltaEta * cellEnergy;
361  totalEnergy[layer] += cellEnergy;
362  }
363 
364  for (int layer=0; layer < 3; ++layer) {
365  if (totalEnergy[layer] != 0.) {
366  deltaEtaSecondMom[layer]/=std::abs(totalEnergy[layer]);
367  }
368  else {
369  deltaEtaSecondMom[layer]=0.;
370  }
371  }
372 
373  return deltaEtaSecondMom;
374 }

◆ getEM1CoreFrac()

float TauPi0ClusterCreator::getEM1CoreFrac ( const xAOD::CaloCluster &  cluster) const
private

fraction of cluster enegry in central EM1 cells

Definition at line 269 of file TauPi0ClusterCreator.cxx.

269  {
270  float coreEnergyEM1 = 0.;
271  float totalEnergyEM1 = 0.;
272 
273  const CaloClusterCellLink* cellLinks = cluster.getCellLinks();
274  CaloClusterCellLink::const_iterator cellLink = cellLinks->begin();
275  for (; cellLink != cellLinks->end(); ++cellLink) {
276  const CaloCell* cell = static_cast<const CaloCell*>(*cellLink);
277 
278  // Only consider EM1
279  int sampling = cell->caloDDE()->getSampling();
280  if (sampling != 1 && sampling != 5) continue;
281 
282  // Only consider positive cells
283  // FIXME: is the weight needed ?
284  float cellEnergy = cell->e() * cellLink.weight();
285  if (cellEnergy <= 0) continue;
286 
287  totalEnergyEM1 += cellEnergy;
288 
289  float deltaEta = cell->eta() - cluster.eta();
290  float deltaPhi = P4Helpers::deltaPhi(cell->phi(), cluster.phi());
291 
292  // Core region: [0.05, 0.05/8]
293  if(std::abs(deltaPhi) > 0.05 || std::abs(deltaEta) > 2 * 0.025/8.) continue;
294 
295  coreEnergyEM1 += cellEnergy;
296  }
297 
298  if (totalEnergyEM1 <= 0.) return 0.;
299  return coreEnergyEM1/totalEnergyEM1;
300 }

◆ getKey()

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

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

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

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

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

Definition at line 119 of file AsgTool.cxx.

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

◆ getName()

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

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

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

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

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

Definition at line 106 of file AsgTool.cxx.

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

◆ getNPhotons()

int TauPi0ClusterCreator::getNPhotons ( const std::vector< const xAOD::PFO * > &  shotVector,
const std::vector< unsigned > &  shotsInCluster 
) const
private

Definition at line 230 of file TauPi0ClusterCreator.cxx.

231  {
232  int totalPhotons = 0;
233 
234  for (unsigned index = 0; index < shotsInCluster.size(); ++index) {
235  int nPhotons = 0;
236  const xAOD::PFO* shotPFO = shotPFOs.at(shotsInCluster.at(index));
237  if (! shotPFO->attribute(xAOD::PFODetails::PFOAttributes::tauShots_nPhotons, nPhotons)) {
238  ATH_MSG_WARNING("Can't find NHitsInEM1. Set it to 0.");
239  }
240  totalPhotons += nPhotons;
241  }
242 
243  return totalPhotons;
244 }

◆ getNPosECells()

std::vector< int > TauPi0ClusterCreator::getNPosECells ( const xAOD::CaloCluster &  cluster) const
private

number of cells from cluster with positive energy in PS, EM1 and EM2

Definition at line 248 of file TauPi0ClusterCreator.cxx.

248  {
249  std::vector<int> nPosECells(3, 0);
250 
251  const CaloClusterCellLink* cellLinks = cluster.getCellLinks();
252  CaloClusterCellLink::const_iterator cellLink = cellLinks->begin();
253  for (; cellLink != cellLinks->end(); ++cellLink) {
254  const CaloCell* cell = static_cast<const CaloCell*>(*cellLink);
255  int sampling = cell->caloDDE()->getSampling();
256 
257  // layer0: PS, layer1: EM1, layer2: EM2
258  int layer = sampling%4;
259  if (layer < 3 && cell->e() > 0) {
260  ++nPosECells[layer];
261  }
262  }
263 
264  return nPosECells;
265 }

◆ getProperty()

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

Get one of the tool's properties.

◆ getShotsMatchedToCluster()

std::vector< unsigned > TauPi0ClusterCreator::getShotsMatchedToCluster ( const std::vector< const xAOD::PFO * > &  shotVector,
const std::map< unsigned, const xAOD::CaloCluster * > &  clusterToShotMap,
const xAOD::CaloCluster &  pi0Cluster 
) const
private

Definition at line 211 of file TauPi0ClusterCreator.cxx.

213  {
214  std::vector<unsigned> shotsMatchedToCluster;
215 
216  // Loop over the shots, and select those matched to the cluster
217  for (unsigned index = 0; index < shotPFOs.size(); ++index) {
218  auto iterator = shotToClusterMap.find(index);
219  if (iterator == shotToClusterMap.end()) continue;
220  if (iterator->second != &pi0Cluster) continue;
221 
222  shotsMatchedToCluster.push_back(index);
223  }
224 
225  return shotsMatchedToCluster;
226 }

◆ getShotToClusterMap()

std::map< unsigned, const xAOD::CaloCluster * > TauPi0ClusterCreator::getShotToClusterMap ( const std::vector< const xAOD::PFO * > &  shotVector,
const xAOD::CaloClusterContainer &  pi0ClusterContainer,
const xAOD::TauJet &  pTau 
) const
private

Definition at line 133 of file TauPi0ClusterCreator.cxx.

135  {
136  std::map<unsigned, const xAOD::CaloCluster*> shotToClusterMap;
137  for (unsigned index = 0; index < shotPFOs.size(); ++index) {
138  const xAOD::PFO* shotPFO = shotPFOs.at(index);
139 
140  int seedHashInt = -1;
141  if (!shotPFO->attribute(xAOD::PFODetails::PFOAttributes::tauShots_seedHash, seedHashInt)) {
142  ATH_MSG_WARNING("Couldn't find seed hash. Set it to -1, no cluster will be associated to shot.");
143  }
144  const IdentifierHash seedHash = static_cast<const IdentifierHash>(seedHashInt);
145 
146  // We will always perform the vertex correction
147  const xAOD::Vertex* vertex = tau.vertex();
148 
149  // Tau custom PFO reconstruction is only used in offline reconstrution
150  TLorentzVector tauAxis = tauRecTools::getTauAxis(tau);
151 
152  float weightInCluster = -1.;
153  float weightInPreviousCluster = -1;
154 
155  // Loop over custom pi0 clusters, and map shot to the cluster
156  for (const xAOD::CaloCluster* cluster: pi0ClusterContainer) {
157  // custom clusters could be correctd directly using the tau vertex
158  TLorentzVector clusterP4;
159  if (vertex) {
160  xAOD::CaloVertexedTopoCluster vertexedCluster(*cluster, vertex->position());
161  clusterP4 = vertexedCluster.p4();
162  }
163  else {
164  clusterP4 = cluster->p4();
165  }
166 
167  weightInCluster = -1.;
168  if (clusterP4.Et() < m_clusterEtCut) continue;
169  if (clusterP4.DeltaR(tauAxis) > 0.4) continue;
170 
171  const CaloClusterCellLink* cellLinks = cluster->getCellLinks();
172  CaloClusterCellLink::const_iterator cellLink = cellLinks->begin();
173  for (; cellLink != cellLinks->end(); ++cellLink) {
174  const CaloCell* cell = static_cast<const CaloCell*>(*cellLink);
175 
176  // Check if seed cell is in cluster.
177  if (cell->caloDDE()->calo_hash() != seedHash) continue;
178 
179  weightInCluster = cellLink.weight();
180  // found cell, no need to loop over other cells
181  break;
182  }
183 
184  if (weightInCluster < 0) continue;
185 
186  // Check if cell was already found in a previous cluster
187  if (weightInPreviousCluster < 0) {
188  // Cell not found in a previous cluster.
189  // Have to check whether cell is shared with other cluster
190  shotToClusterMap[index] = cluster;
191  weightInPreviousCluster = weightInCluster;
192  }
193  else {
194  // Cell has been found in a previous cluster
195  // assign shot to this cluster if it has larger weight for the cell
196  // otherwise the shots keeps assigned to the previous cluster
197  if (weightInCluster > weightInPreviousCluster) {
198  shotToClusterMap[index] = cluster;
199  }
200  // FIXME: why break here ? Should loop all the cluster, and find the largest weight
201  break;
202  }
203  }
204  }
205 
206  return shotToClusterMap;
207 }

◆ inAOD()

bool TauRecToolBase::inAOD ( ) const
inlineprotectedinherited

Definition at line 88 of file TauRecToolBase.h.

88 { return m_in_AOD; }

◆ inEleRM()

bool TauRecToolBase::inEleRM ( ) const
inlineprotectedinherited

Definition at line 89 of file TauRecToolBase.h.

89 { return m_in_EleRM; }

◆ initialize()

StatusCode TauRecToolBase::initialize ( )
overridevirtualinherited

Tool initializer.

Reimplemented from asg::AsgTool.

Reimplemented in tauRecTools::TrackRNN, TauTrackFinder, tauRecTools::TauTrackRNNClassifier, TauWPDecorator, TauVertexFinder, TauDecayModeNNClassifier, TauPi0CreateROI, TauShotFinder, TauElectronVetoVariables, TauEleOverlapChecker, TauJetRNNEvaluator, TauCalibrateLC, TauGNNEvaluator, TauPi0ScoreCalculator, TauVertexVariables, TauVertexedClusterDecorator, MvaTESVariableDecorator, TauAODMuonRemovalTool, TauCombinedTES, MvaTESEvaluator, and PanTau::PanTauProcessor.

Definition at line 129 of file TauRecToolBase.cxx.

129  {
130  return StatusCode::SUCCESS;
131 }

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

◆ msg() [1/2]

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

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

Definition at line 27 of file AthCommonMsg.h.

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

◆ msg_level_name()

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

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

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

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

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

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

Definition at line 101 of file AsgTool.cxx.

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

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

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

◆ outputHandles()

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

Return this algorithm's output handles.

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

◆ print()

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

Print the state of the tool.

Implements asg::IAsgTool.

Reimplemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.

Definition at line 131 of file AsgTool.cxx.

131  {
132 
133  ATH_MSG_INFO( "AsgTool " << name() << " @ " << this );
134  return;
135  }

◆ 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  {
66  StatusCode sc;
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 }

◆ renounce()

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

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ sysInitialize()

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

Perform system initialization for an algorithm.

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

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

◆ sysStart()

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

Handle START transition.

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

◆ updateVHKA()

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

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_clusterEtCut

Gaudi::Property<double> TauPi0ClusterCreator::m_clusterEtCut {this, "ClusterEtCut", 0.5 * Gaudi::Units::GeV, "Et threshould for pi0 candidate clusters"}
private

Definition at line 82 of file TauPi0ClusterCreator.h.

◆ m_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

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

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

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

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

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.

The documentation for this class was generated from the following files:
xAOD::iterator
JetConstituentVector::iterator iterator
Definition: JetConstituentVector.cxx:68
xAOD::CaloCluster_v1::phi
virtual double phi() const
The azimuthal angle ( ) of the particle.
Definition: CaloCluster_v1.cxx:256
Property::BOOL
@ BOOL
Definition: Property.h:29
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
Property::Type
Type
Property type enumeration.
Definition: Property.h:27
TauGNNUtils::Variables::Cluster::CENTER_LAMBDA
bool CENTER_LAMBDA(const xAOD::TauJet &, const xAOD::CaloVertexedTopoCluster &cluster, double &out)
Definition: TauGNNUtils.cxx:840
xAOD::PFODetails::cellBased_firstEtaWRTClusterPosition_EM2
@ cellBased_firstEtaWRTClusterPosition_EM2
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:137
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
xAOD::PFODetails::cellBased_NPosECells_EM1
@ cellBased_NPosECells_EM1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:134
TauPi0ClusterCreator::get2ndEtaMomWRTCluster
std::vector< float > get2ndEtaMomWRTCluster(const xAOD::CaloCluster &cluster) const
second eta moment in PS, EM1 and EM2 w.r.t cluster eta
Definition: TauPi0ClusterCreator.cxx:341
constants.EMB1
int EMB1
Definition: Calorimeter/CaloClusterCorrection/python/constants.py:53
StateLessPT_NewConfig.proxy
proxy
Definition: StateLessPT_NewConfig.py:392
ReadCellNoiseFromCool.cell
cell
Definition: ReadCellNoiseFromCool.py:53
xAOD::PFODetails::cellBased_SECOND_R
@ cellBased_SECOND_R
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:119
TauPi0ClusterCreator::configureHadronicPFO
StatusCode configureHadronicPFO(const xAOD::CaloVertexedTopoCluster &cluster, double clusterEnergyHad, xAOD::PFO &hadronicPFO) const
Configure the haronic PFO.
Definition: TauPi0ClusterCreator.cxx:479
Property::STRING
@ STRING
Definition: Property.h:33
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
xAOD::PFODetails::TauShot
@ TauShot
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:171
xAOD::PFODetails::PFOAttributes
PFOAttributes
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:28
TauPi0ClusterCreator::m_clusterEtCut
Gaudi::Property< double > m_clusterEtCut
Definition: TauPi0ClusterCreator.h:82
index
Definition: index.py:1
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
tauRecTools::getTauAxis
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 ...
Definition: Reconstruction/tauRecTools/Root/HelperFunctions.cxx:33
xAOD::deltaPhi
setSAddress setEtaMS setDirPhiMS setDirZMS setBarrelRadius setEndcapAlpha setEndcapRadius setInterceptInner setEtaMap setEtaBin setIsTgcFailure setDeltaPt deltaPhi
Definition: L2StandAloneMuon_v1.cxx:160
xAOD::PFO_v1::attribute
bool attribute(PFODetails::PFOAttributes AttributeType, T &anAttribute) const
get a PFO Variable via enum
xAOD::PFODetails::cellBased_NPosECells_PS
@ cellBased_NPosECells_PS
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:133
xAOD::PFODetails::cellBased_NHitsInEM1
@ cellBased_NHitsInEM1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:132
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
TauRecToolBase::m_tauRecToolsTag
Gaudi::Property< std::string > m_tauRecToolsTag
Definition: TauRecToolBase.h:79
athena.value
value
Definition: athena.py:124
xAOD::CaloVertexedClusterBase::phi
virtual double phi() const final
The azimuthal angle ( ) of the particle.
Definition: Event/xAOD/xAODCaloEvent/xAODCaloEvent/CaloVertexedClusterBase.h:79
asg::AsgTool::getProperty
const T * getProperty(const std::string &name) const
Get one of the tool's properties.
xAOD::CaloVertexedClusterBase::eta
virtual double eta() const final
The pseudorapidity ( ) of the particle.
Definition: Event/xAOD/xAODCaloEvent/xAODCaloEvent/CaloVertexedClusterBase.h:77
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
dbg::ptr
void * ptr(T *p)
Definition: SGImplSvc.cxx:74
xAOD::PFO_v1::setAssociatedParticleLinks
bool setAssociatedParticleLinks(PFODetails::PFOParticleType ParticleType, const std::vector< ElementLink< IParticleContainer > > &theParticles)
Set a vector of PFO constituent particle types via enum - overwrite is allowed.
Definition: PFO_v1.cxx:600
xAOD::PFODetails::cellBased_secondEtaWRTClusterPosition_EM1
@ cellBased_secondEtaWRTClusterPosition_EM1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:138
empty
bool empty(TH1 *h)
Definition: computils.cxx:295
xAOD::PFODetails::cellBased_DELTA_THETA
@ cellBased_DELTA_THETA
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:122
xAOD::PFODetails::cellBased_ENG_FRAC_MAX
@ cellBased_ENG_FRAC_MAX
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:127
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
CaloClusterCellLink::const_iterator::weight
weight_t weight() const
Accessor for weight associated to this cell.
Definition: CaloClusterCellLink.h:82
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
xAOD::PFODetails::cellBased_firstEtaWRTClusterPosition_EM1
@ cellBased_firstEtaWRTClusterPosition_EM1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:136
xAOD::PFODetails::tauShots_nPhotons
@ tauShots_nPhotons
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:162
xAOD::CaloCluster_v1::MomentType
MomentType
Enums to identify different moments.
Definition: CaloCluster_v1.h:120
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
xAOD::PFO_v1::setClusterLink
bool setClusterLink(const ElementLink< xAOD::CaloClusterContainer > &theCluster)
Set a cluster constituent - does NOT append to existing container
Definition: PFO_v1.cxx:549
Property::DOUBLE
@ DOUBLE
Definition: Property.h:32
xAOD::PFODetails::cellBased_SECOND_LAMBDA
@ cellBased_SECOND_LAMBDA
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:120
AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:85
xAOD::CaloCluster_v1
Description of a calorimeter cluster.
Definition: CaloCluster_v1.h:59
P4Helpers::deltaPhi
double deltaPhi(double phiA, double phiB)
delta Phi in range [-pi,pi[
Definition: P4Helpers.h:34
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
AthCommonDataStore
Definition: AthCommonDataStore.h:52
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
constants.EMB2
int EMB2
Definition: Calorimeter/CaloClusterCorrection/python/constants.py:54
P4Helpers::deltaEta
double deltaEta(const I4Momentum &p1, const I4Momentum &p2)
Computes efficiently .
Definition: P4Helpers.h:66
xAOD::CaloCluster_v1::eta
virtual double eta() const
The pseudorapidity ( ) of the particle.
Definition: CaloCluster_v1.cxx:251
lumiFormat.i
int i
Definition: lumiFormat.py:85
xAOD::PFO
PFO_v1 PFO
Definition of the current "pfo version".
Definition: PFO.h:17
xAOD::TauJet_v3::shotPFOLinks
const PFOLinks_t & shotPFOLinks() const
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
TRT::Hit::layer
@ layer
Definition: HitInfo.h:79
xAOD::PFODetails::cellBased_secondEtaWRTClusterPosition_EM2
@ cellBased_secondEtaWRTClusterPosition_EM2
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:139
tauRecTools::doPi0andShots
bool doPi0andShots(const xAOD::TauJet &tau)
Determines whether pi0s and shots should be built for a tau candidate.
Definition: Reconstruction/tauRecTools/Root/HelperFunctions.cxx:93
xAOD::PFODetails::tauShots_seedHash
@ tauShots_seedHash
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:163
TauRecToolBase::m_in_AOD
Gaudi::Property< bool > m_in_AOD
Definition: TauRecToolBase.h:77
constants.EME1
int EME1
Definition: Calorimeter/CaloClusterCorrection/python/constants.py:55
python.xAODType.dummy
dummy
Definition: xAODType.py:4
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
MSG::name
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition: MsgLevel.cxx:19
xAOD::CaloCluster_v1::retrieveMoment
bool retrieveMoment(MomentType type, double &value) const
Retrieve individual moment.
Definition: CaloCluster_v1.cxx:738
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
xAOD::PFODetails::cellBased_FIRST_ETA
@ cellBased_FIRST_ETA
These variables belong to the cell-based particle flow algorithm.
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:118
TauRecToolBase::TauRecToolBase
TauRecToolBase(const std::string &name)
Definition: TauRecToolBase.cxx:125
TauGNNUtils::Variables::Cluster::SECOND_LAMBDA
bool SECOND_LAMBDA(const xAOD::TauJet &, const xAOD::CaloVertexedTopoCluster &cluster, double &out)
Definition: TauGNNUtils.cxx:835
xAOD::PFO_v1::setBDTPi0Score
void setBDTPi0Score(float BDTPi0Score)
set BDT Score used to classify clusters as Pi0 like or not
xAOD::CaloCluster_v1::getCellLinks
const CaloClusterCellLink * getCellLinks() const
Get a pointer to the CaloClusterCellLink object (const version)
Definition: CaloCluster_v1.cxx:905
TauRecToolBase::m_in_trigger
Gaudi::Property< bool > m_in_trigger
Definition: TauRecToolBase.h:76
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
Property::INT
@ INT
Definition: Property.h:30
xAOD::PFODetails::cellBased_energy_EM2
@ cellBased_energy_EM2
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:141
xAOD::PFODetails::cellBased_energy_EM1
@ cellBased_energy_EM1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:140
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
xAOD::PFO_v1
Class describing a particle flow object.
Definition: PFO_v1.h:35
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
Property::FLOAT
@ FLOAT
Definition: Property.h:31
TauRecToolBase::find_file
std::string find_file(const std::string &fname) const
Definition: TauRecToolBase.cxx:19
TauGNNUtils::Variables::Cluster::SECOND_R
bool SECOND_R(const xAOD::TauJet &, const xAOD::CaloVertexedTopoCluster &cluster, double &out)
Definition: TauGNNUtils.cxx:830
DataVector::push_back
value_type push_back(value_type pElem)
Add an element to the end of the collection.
xAOD::PFODetails::nPi0Proto
@ nPi0Proto
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:31
xAOD::PFODetails::cellBased_DELTA_PHI
@ cellBased_DELTA_PHI
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:121
xAOD::PFODetails::cellBased_ENG_FRAC_CORE
@ cellBased_ENG_FRAC_CORE
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:128
xAOD::PFODetails::cellBased_EM1CoreFrac
@ cellBased_EM1CoreFrac
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:130
TauGNNUtils::Variables::Cluster::CENTER_MAG
bool CENTER_MAG(const xAOD::TauJet &, const xAOD::CaloVertexedTopoCluster &cluster, double &out)
Definition: TauGNNUtils.cxx:916
PathResolverFindCalibFile
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
Definition: PathResolver.cxx:431
python.AthDsoLogger.fname
string fname
Definition: AthDsoLogger.py:67
xAOD::PFO_v1::setP4
void setP4(const FourMom_t &vec)
set the 4-vec
Definition: PFO_v1.cxx:107
xAOD::PFO_v1::setCharge
void setCharge(float charge)
set charge of PFO
Trk::vertex
@ vertex
Definition: MeasurementType.h:21
DeMoScan.index
string index
Definition: DeMoScan.py:364
xAOD::PFODetails::cellBased_SECOND_ENG_DENS
@ cellBased_SECOND_ENG_DENS
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:129
xAOD::PFODetails::cellBased_NPosECells_EM2
@ cellBased_NPosECells_EM2
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:135
xAOD::CaloCluster_v1::eSample
float eSample(const CaloSample sampling) const
Definition: CaloCluster_v1.cxx:521
a
TList * a
Definition: liststreamerinfos.cxx:10
h
xAOD::Vertex_v1
Class describing a Vertex.
Definition: Vertex_v1.h:42
CaloCell
Data object for each calorimeter readout cell.
Definition: CaloCell.h:57
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
xAOD::PFODetails::cellBased_CENTER_LAMBDA
@ cellBased_CENTER_LAMBDA
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:123
TauPi0ClusterCreator::configureNeutralPFO
StatusCode configureNeutralPFO(const xAOD::CaloCluster &cluster, const xAOD::CaloClusterContainer &pi0ClusterContainer, const xAOD::TauJet &tau, const std::vector< const xAOD::PFO * > &shotPFOs, const std::map< unsigned, const xAOD::CaloCluster * > &shotsInCluster, xAOD::PFO &neutralPFO) const
Configure the neutral PFO.
Definition: TauPi0ClusterCreator.cxx:378
xAOD::PFODetails::cellBased_LONGITUDINAL
@ cellBased_LONGITUDINAL
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:125
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
TauPi0ClusterCreator::getShotToClusterMap
std::map< unsigned, const xAOD::CaloCluster * > getShotToClusterMap(const std::vector< const xAOD::PFO * > &shotVector, const xAOD::CaloClusterContainer &pi0ClusterContainer, const xAOD::TauJet &pTau) const
Definition: TauPi0ClusterCreator.cxx:133
xAOD::PFODetails::cellBased_LATERAL
@ cellBased_LATERAL
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:124
AthCommonMsg< AlgTool >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
TauPi0ClusterCreator::getNPhotons
int getNPhotons(const std::vector< const xAOD::PFO * > &shotVector, const std::vector< unsigned > &shotsInCluster) const
Definition: TauPi0ClusterCreator.cxx:230
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
xAOD::PFODetails::cellBased_ENG_FRAC_EM
@ cellBased_ENG_FRAC_EM
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:126
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
TauPi0ClusterCreator::get1stEtaMomWRTCluster
std::vector< float > get1stEtaMomWRTCluster(const xAOD::CaloCluster &cluster) const
first eta moment in PS, EM1 and EM2 w.r.t cluster eta
Definition: TauPi0ClusterCreator.cxx:304
xAOD::CaloVertexedTopoCluster
Evaluate cluster kinematics with a different vertex / signal state.
Definition: Event/xAOD/xAODCaloEvent/xAODCaloEvent/CaloVertexedTopoCluster.h:38
IdentifierHash
This is a "hash" representation of an Identifier. This encodes a 32 bit index which can be used to lo...
Definition: IdentifierHash.h:25
TauPi0ClusterCreator::getEM1CoreFrac
float getEM1CoreFrac(const xAOD::CaloCluster &cluster) const
fraction of cluster enegry in central EM1 cells
Definition: TauPi0ClusterCreator.cxx:269
python.DataFormatRates.env
env
Definition: DataFormatRates.py:32
xAOD::PFO_v1::setAttribute
void setAttribute(PFODetails::PFOAttributes AttributeType, const T &anAttribute)
Set a PFO Variable via enum - overwrite is allowed.
SG::DataProxy
Definition: DataProxy.h:45
TauPi0ClusterCreator::getNPosECells
std::vector< int > getNPosECells(const xAOD::CaloCluster &cluster) const
number of cells from cluster with positive energy in PS, EM1 and EM2
Definition: TauPi0ClusterCreator.cxx:248
CaloClusterCellLink::const_iterator
const_iterator to loop over cells belonging to a cluster
Definition: CaloClusterCellLink.h:57
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
TauRecToolBase::m_in_EleRM
Gaudi::Property< bool > m_in_EleRM
Definition: TauRecToolBase.h:78
constants.EME2
int EME2
Definition: Calorimeter/CaloClusterCorrection/python/constants.py:56
fitman.k
k
Definition: fitman.py:528
TauPi0ClusterCreator::getShotsMatchedToCluster
std::vector< unsigned > getShotsMatchedToCluster(const std::vector< const xAOD::PFO * > &shotVector, const std::map< unsigned, const xAOD::CaloCluster * > &clusterToShotMap, const xAOD::CaloCluster &pi0Cluster) const
Definition: TauPi0ClusterCreator.cxx:211
xAOD::PFO_v1::setCenterMag
void setCenterMag(float CenterMag)
set CenterMag moment needed for vertex correction
Generated on Sun Dec 22 2024 21:45:28 for ATLAS Offline Software by   doxygen 1.8.18