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

#include <TauShotFinder.h>

Inheritance diagram for TauShotFinder:
Collaboration diagram for TauShotFinder:

Classes

struct  ptSort

Public Member Functions

 ASG_TOOL_CLASS2 (TauShotFinder, TauRecToolBase, ITauToolBase)
 TauShotFinder (const std::string &name)
virtual ~TauShotFinder ()=default
virtual StatusCode initialize () override
 Tool initializer.
virtual StatusCode executeShotFinder (xAOD::TauJet &pTau, xAOD::CaloClusterContainer &tauShotCaloClusContainer, xAOD::PFOContainer &tauShotPFOContainer) const override
virtual StatusCode eventInitialize () override
 Event initializer - called at the beginning of each event.
virtual StatusCode execute (xAOD::TauJet &pTau) const override
 Execute - called for each tau candidate.
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 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

StatusCode selectCells (const xAOD::TauJet &tau, std::vector< const CaloCell * > &cells) const
 Apply preselection of the cells Cells within dR < 0.4, in EM1, and pt > 100 MeV are selected.
StatusCode selectSeedCells (const xAOD::TauJet &tau, const CaloCellContainer &cellContainer, std::vector< const CaloCell * > &seedCells) const
 Select the seed cells used to construct the shot Cells must sastisfy:
bool isPhiNeighbour (IdentifierHash cell1Hash, IdentifierHash cell2Hash) const
 Check whether two cells are neighbours in the phi direction.
const CaloCellgetPhiNeighbour (const CaloCell &seedCell, const std::vector< const CaloCell * > &seedCells) const
 Get the hottest neighbour cell in the phi direction.
std::vector< const CaloCell * > getEtaNeighbours (const CaloCell &cell, const CaloCellContainer &cellContainer, int maxDepth) const
 Get neighbour cells in the eta direction.
void addEtaNeighbours (const CaloCell &cell, const CaloCellContainer &cellContainer, std::vector< const CaloCell * > &cells, int depth, int maxDepth, bool next) const
 Get neighbour cells in the eta direction.
xAOD::CaloClustercreateShotCluster (const CaloCell *cell, const CaloCell *phiNeighCell, const CaloCellContainer &cellContainer, xAOD::CaloClusterContainer *clusterContainer) const
 Create the shot cluster Shot cluster contains 5x1 cells from the seed cell and hottestneighbour cell in the phi direction.
int getEtaBin (float eta) const
 Get eta bin.
int getNPhotons (float eta, float energy) const
 Get NPhotons in shot.
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

Gaudi::Property< int > m_nCellsInEta {this, "NCellsInEta"}
Gaudi::Property< std::vector< float > > m_minPtCut {this, "MinPtCut"}
Gaudi::Property< std::vector< float > > m_doubleShotCut {this, "AutoDoubleShotCut"}
Gaudi::Property< bool > m_removeElectronCells {this, "RemoveElectronCells", false}
Gaudi::Property< float > m_dRThreshold {this, "DRThreshold", 0.4, "dR(cell,tau) to select cells"}
Gaudi::Property< float > m_energyThreshold {this, "EnergyThreshold", 100., "energy threshold (in MeV) to select cells"}
SG::ReadHandleKey< CaloCellContainerm_caloCellInputContainer {this,"Key_caloCellInputContainer", "AllCalo", "input calo cell container key"}
ToolHandle< IHadronicCalibrationToolm_caloWeightTool {this, "CaloWeightTool", "H1WeightToolCSC12Generic"}
SG::ReadHandleKey< xAOD::CaloClusterContainerm_removedClusterInputContainer {this,"Key_RemovedClusterInputContainer", "", "input removed cluster key"}
const CaloCell_IDm_calo_id = nullptr
 calo cell navigation
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 30 of file TauShotFinder.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

◆ TauShotFinder()

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

Definition at line 16 of file TauShotFinder.cxx.

16 :
17 TauRecToolBase(name) {}
TauRecToolBase(const std::string &name)

◆ ~TauShotFinder()

virtual TauShotFinder::~TauShotFinder ( )
virtualdefault

Member Function Documentation

◆ addEtaNeighbours()

void TauShotFinder::addEtaNeighbours ( const CaloCell & cell,
const CaloCellContainer & cellContainer,
std::vector< const CaloCell * > & cells,
int depth,
int maxDepth,
bool next ) const
private

Get neighbour cells in the eta direction.

Definition at line 329 of file TauShotFinder.cxx.

334 {
335 ++depth;
336
337 if (depth > maxDepth) return;
338
339 const IdentifierHash cellHash = cell.caloDDE()->calo_hash();
340
341 std::vector<IdentifierHash> neigHashes;
342 if (next) {
343 m_calo_id->get_neighbours(cellHash,LArNeighbours::nextInEta,neigHashes);
344 }
345 else {
346 m_calo_id->get_neighbours(cellHash,LArNeighbours::prevInEta,neigHashes);
347 }
348
349 for (const IdentifierHash& hash : neigHashes) {
350 const CaloCell* newCell = cellContainer.findCell(hash);
351
352 if (!newCell) continue;
353
354 cells.push_back(newCell);
355 this->addEtaNeighbours(*newCell, cellContainer, cells, depth, maxDepth, next);
356
357 if (neigHashes.size() > 1) {
358 ATH_MSG_DEBUG(cellHash << " has " << neigHashes.size() << " neighbours in the eta direction !");
359 break;
360 }
361 }
362}
#define ATH_MSG_DEBUG(x)
const CaloCell * findCell(const IdentifierHash theHash) const
fast find method given identifier hash.
void addEtaNeighbours(const CaloCell &cell, const CaloCellContainer &cellContainer, std::vector< const CaloCell * > &cells, int depth, int maxDepth, bool next) const
Get neighbour cells in the eta direction.
const CaloCell_ID * m_calo_id
calo cell navigation
std::string depth
tag string for intendation
Definition fastadd.cxx:46

◆ ASG_TOOL_CLASS2()

TauShotFinder::ASG_TOOL_CLASS2 ( TauShotFinder ,
TauRecToolBase ,
ITauToolBase  )

◆ createShotCluster()

xAOD::CaloCluster * TauShotFinder::createShotCluster ( const CaloCell * cell,
const CaloCell * phiNeighCell,
const CaloCellContainer & cellContainer,
xAOD::CaloClusterContainer * clusterContainer ) const
private

Create the shot cluster Shot cluster contains 5x1 cells from the seed cell and hottestneighbour cell in the phi direction.

Definition at line 366 of file TauShotFinder.cxx.

369 {
370
371 xAOD::CaloCluster* shotCluster = CaloClusterStoreHelper::makeCluster(clusterContainer,&cellContainer);
372
373 int maxDepth = (m_nCellsInEta - 1) / 2;
374
375 std::vector<const CaloCell*> windowNeighbours = this->getEtaNeighbours(*cell, cellContainer, maxDepth);
376 if (phiNeigCell) {
377 std::vector<const CaloCell*> mergeCells = this->getEtaNeighbours(*phiNeigCell, cellContainer, maxDepth);
378 windowNeighbours.push_back(phiNeigCell);
379 windowNeighbours.insert(windowNeighbours.end(), mergeCells.begin(), mergeCells.end());
380 }
381
382 shotCluster->getOwnCellLinks()->reserve(windowNeighbours.size()+1);
383 const IdentifierHash seedHash = cell->caloDDE()->calo_hash();
384 shotCluster->addCell(cellContainer.findIndex(seedHash), 1.);
385
386 for (const CaloCell* cell : windowNeighbours) {
387 shotCluster->addCell(cellContainer.findIndex(cell->caloDDE()->calo_hash()), 1.0);
388 }
389
390 CaloClusterKineHelper::calculateKine(shotCluster,true,true);
391
392 return shotCluster;
393}
int findIndex(const IdentifierHash theHash) const
Return index of the cell with a given hash.
static void calculateKine(xAOD::CaloCluster *clu, const bool useweight=true, const bool updateLayers=true, const bool useGPUCriteria=false)
Helper class to calculate cluster kinematics based on cells.
static std::unique_ptr< xAOD::CaloCluster > makeCluster(const CaloCellContainer *cellCont)
Creates a valid CaloCluster with a private Aux-Store and CellLink container.
std::vector< const CaloCell * > getEtaNeighbours(const CaloCell &cell, const CaloCellContainer &cellContainer, int maxDepth) const
Get neighbour cells in the eta direction.
Gaudi::Property< int > m_nCellsInEta
CaloClusterCellLink * getOwnCellLinks()
Get a pointer to the owned CaloClusterCellLink object (non-const version)
bool addCell(const unsigned index, const double weight)
Method to add a cell to the cluster (Beware: Kinematics not updated!)
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.

◆ 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 TauRecToolBase::execute ( xAOD::TauJet & pTau) const
overridevirtualinherited

◆ 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 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 TauShotFinder::executeShotFinder ( xAOD::TauJet & pTau,
xAOD::CaloClusterContainer & tauShotCaloClusContainer,
xAOD::PFOContainer & tauShotPFOContainer ) const
overridevirtual

Reimplemented from TauRecToolBase.

Definition at line 31 of file TauShotFinder.cxx.

32 {
33
34 // Any tau needs to have shot PFO vectors. Set empty vectors before nTrack cut
35 std::vector<ElementLink<xAOD::PFOContainer>> empty;
36 tau.setShotPFOLinks(empty);
37
38 // Only run on 0-5 prong taus
40 return StatusCode::SUCCESS;
41 }
42
43 SG::ReadHandle<CaloCellContainer> caloCellInHandle( m_caloCellInputContainer );
44 if (!caloCellInHandle.isValid()) {
45 ATH_MSG_ERROR ("Could not retrieve HiveDataObj with key " << caloCellInHandle.key());
46 return StatusCode::FAILURE;
47 }
48 const CaloCellContainer *cellContainer = caloCellInHandle.cptr();
49
50 // Select seed cells:
51 // -- dR < 0.4, EM1, pt > 100
52 // -- largest pt among the neighbours in eta direction
53 // -- no other seed cell as neighbour in eta direction
54 std::vector<const CaloCell*> seedCells;
55 ATH_CHECK(selectSeedCells(tau, *cellContainer, seedCells));
56 ATH_MSG_DEBUG("seedCells.size() = " << seedCells.size());
57
58 // Construt shot by merging neighbour cells in phi direction
59 while (!seedCells.empty()) {
60 // Find the neighbour in phi direction, and choose the one with highest pt
61 const CaloCell* cell = seedCells.front();
62 const CaloCell* phiNeigCell = getPhiNeighbour(*cell, seedCells);
63
64 // Construct shot PFO candidate
65 xAOD::PFO* shot = new xAOD::PFO();
66 shotPFOContainer.push_back(shot);
67
68 // -- Construct the shot cluster
69 xAOD::CaloCluster* shotCluster = createShotCluster(cell, phiNeigCell, *cellContainer, &shotClusterContainer);
70
71 ElementLink<xAOD::CaloClusterContainer> clusElementLink;
72 clusElementLink.toContainedElement( shotClusterContainer, shotCluster );
73 shot->setClusterLink( clusElementLink );
74
75 // -- Calculate the four momentum
76 // TODO: simplify the calculation
77 if (phiNeigCell) {
78 // interpolate position
79 double wtCell_phiNeigCell = m_caloWeightTool->wtCell(phiNeigCell);
80 double wtCell_cell = m_caloWeightTool->wtCell(cell);
81 double dPhi = TVector2::Phi_mpi_pi( phiNeigCell->phi() - cell->phi());
82 double ratio = phiNeigCell->pt()*wtCell_phiNeigCell/(cell->pt()*wtCell_cell + phiNeigCell->pt()*wtCell_phiNeigCell);
83 float phi = cell->phi()+dPhi*ratio;
84 float pt = cell->pt()*wtCell_cell+phiNeigCell->pt()*wtCell_phiNeigCell;
85
86 shot->setP4( static_cast<float>(pt), static_cast<float>(cell->eta()), static_cast<float>(phi), static_cast<float>(cell->m()));
87 }
88 else {
89 shot->setP4( static_cast<float>(cell->pt()*m_caloWeightTool->wtCell(cell)), static_cast<float>(cell->eta()), static_cast<float>(cell->phi()), static_cast<float>(cell->m()));
90 }
91
92 // -- Set the Attribute
93 shot->setBDTPi0Score(-9999.);
94 shot->setCharge(0);
95 shot->setCenterMag(0.0);
96
97 shot->setAttribute<int>(xAOD::PFODetails::PFOAttributes::tauShots_nCellsInEta, m_nCellsInEta);
98
99 const IdentifierHash seedHash = cell->caloDDE()->calo_hash();
100 shot->setAttribute<int>(xAOD::PFODetails::PFOAttributes::tauShots_seedHash, seedHash);
101
102 std::vector<std::vector<const CaloCell*>> cellBlock = TauShotVariableHelpers::getCellBlock(*shot, m_calo_id);
103
104 float pt1 = TauShotVariableHelpers::ptWindow(cellBlock, 1, m_caloWeightTool);
105 shot->setAttribute<float>(xAOD::PFODetails::PFOAttributes::tauShots_pt1, pt1);
106
107 float pt3 = TauShotVariableHelpers::ptWindow(cellBlock, 3, m_caloWeightTool);
108 shot->setAttribute<float>(xAOD::PFODetails::PFOAttributes::tauShots_pt3, pt3);
109
110 float pt5 = TauShotVariableHelpers::ptWindow(cellBlock, 5, m_caloWeightTool);
111 shot->setAttribute<float>(xAOD::PFODetails::PFOAttributes::tauShots_pt5, pt5);
112
113 int nPhotons = getNPhotons(cell->eta(), pt1);
114 shot->setAttribute<int>(xAOD::PFODetails::PFOAttributes::tauShots_nPhotons, nPhotons);
115
116 // Add Element link to the shot PFO container
117 ElementLink<xAOD::PFOContainer> PFOElementLink;
118 PFOElementLink.toContainedElement(shotPFOContainer, shot);
119 tau.addShotPFOLink(PFOElementLink);
120
121 // Remove used cells from list
122 auto cellIndex = std::find(seedCells.begin(), seedCells.end(), cell);
123 seedCells.erase(cellIndex);
124 if (phiNeigCell) {
125 cellIndex = std::find(seedCells.begin(), seedCells.end(), phiNeigCell);
126 seedCells.erase(cellIndex);
127 }
128 } // Loop over seed cells
129
130 return StatusCode::SUCCESS;
131}
Scalar phi() const
phi method
#define ATH_CHECK
Evaluate an expression and check for errors.
static const Attributes_t empty
virtual double phi() const override final
get phi (through CaloDetDescrElement)
Definition CaloCell.h:375
virtual double pt() const
transverse momentum
StatusCode selectSeedCells(const xAOD::TauJet &tau, const CaloCellContainer &cellContainer, std::vector< const CaloCell * > &seedCells) const
Select the seed cells used to construct the shot Cells must sastisfy:
xAOD::CaloCluster * createShotCluster(const CaloCell *cell, const CaloCell *phiNeighCell, const CaloCellContainer &cellContainer, xAOD::CaloClusterContainer *clusterContainer) const
Create the shot cluster Shot cluster contains 5x1 cells from the seed cell and hottestneighbour cell ...
SG::ReadHandleKey< CaloCellContainer > m_caloCellInputContainer
ToolHandle< IHadronicCalibrationTool > m_caloWeightTool
int getNPhotons(float eta, float energy) const
Get NPhotons in shot.
const CaloCell * getPhiNeighbour(const CaloCell &seedCell, const std::vector< const CaloCell * > &seedCells) const
Get the hottest neighbour cell in the phi direction.
void setAttribute(PFODetails::PFOAttributes AttributeType, const T &anAttribute)
Set a PFO Variable via enum - overwrite is allowed.
bool setClusterLink(const ElementLink< xAOD::CaloClusterContainer > &theCluster)
Set a cluster constituent - does NOT append to existing container.
Definition PFO_v1.cxx:549
void setCenterMag(float CenterMag)
set CenterMag moment needed for vertex correction
void setP4(const FourMom_t &vec)
set the 4-vec
Definition PFO_v1.cxx:107
void setBDTPi0Score(float BDTPi0Score)
set BDT Score used to classify clusters as Pi0 like or not
void setCharge(float charge)
set charge of PFO
bool dPhi(const xAOD::TauJet &tau, const xAOD::CaloVertexedTopoCluster &cluster, float &out)
std::vector< std::vector< const CaloCell * > > getCellBlock(const xAOD::PFO &shot, const CaloCell_ID *calo_id)
Get cell block with (currently) 2 x 5 cells in correct order for variable calculations.
float ptWindow(const std::vector< std::vector< const CaloCell * > > &shotCells, int windowSize, const ToolHandle< IHadronicCalibrationTool > &caloWeightTool)
pt in a window of (currently) 2 x windowSize cells
bool doPi0andShots(const xAOD::TauJet &tau)
Determines whether pi0s and shots should be built for a tau candidate.
PFO_v1 PFO
Definition of the current "pfo version".
Definition PFO.h:17

◆ 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

◆ getEtaBin()

int TauShotFinder::getEtaBin ( float eta) const
private

Get eta bin.

Definition at line 135 of file TauShotFinder.cxx.

135 {
136 float absEta=std::abs(eta);
137
138 if (absEta < 0.80) {
139 return 0; // Central Barrel
140 }
141 if (absEta<1.39) {
142 return 1; // Outer Barrel
143 }
144 if (absEta<1.51) {
145 return 2; // Crack region
146 }
147 if (absEta<1.80) {
148 return 3; // Endcap, fine granularity
149 }
150 return 4; // Endcap, coarse granularity
151}
Scalar eta() const
pseudorapidity method
bool absEta(const xAOD::TauJet &tau, float &out)

◆ getEtaNeighbours()

std::vector< const CaloCell * > TauShotFinder::getEtaNeighbours ( const CaloCell & cell,
const CaloCellContainer & cellContainer,
int maxDepth ) const
private

Get neighbour cells in the eta direction.

Definition at line 314 of file TauShotFinder.cxx.

316 {
317 std::vector<const CaloCell*> cells;
318
319 // Add neighbours in next eta direction
320 this->addEtaNeighbours(cell, cellContainer, cells, 0, maxDepth, true);
321 // Add neighbours in previous eta direction
322 this->addEtaNeighbours(cell, cellContainer, cells, 0, maxDepth, false);
323
324 return cells;
325}

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

◆ getNPhotons()

int TauShotFinder::getNPhotons ( float eta,
float energy ) const
private

Get NPhotons in shot.

Definition at line 155 of file TauShotFinder.cxx.

155 {
156 int etaBin = getEtaBin(eta);
157
158 // No photons in crack region
159 if(etaBin==2) return 0;
160
161 const std::vector<float>& minPtCut = m_minPtCut.value();
162 const std::vector<float>& doubleShotCut = m_doubleShotCut.value();
163 ATH_MSG_DEBUG("etaBin = " << etaBin << ", energy = " << energy);
164 ATH_MSG_DEBUG("MinPtCut: " << minPtCut.at(etaBin) << "DoubleShotCut: " << doubleShotCut.at(etaBin));
165
166 if (energy < minPtCut.at(etaBin)) return 0;
167 if (energy > doubleShotCut.at(etaBin)) return 2;
168 return 1;
169}
Gaudi::Property< std::vector< float > > m_minPtCut
Gaudi::Property< std::vector< float > > m_doubleShotCut
int getEtaBin(float eta) const
Get eta bin.
setSAddress setEtaMS setDirPhiMS setDirZMS setBarrelRadius setEndcapAlpha setEndcapRadius setInterceptInner setEtaMap etaBin

◆ getPhiNeighbour()

const CaloCell * TauShotFinder::getPhiNeighbour ( const CaloCell & seedCell,
const std::vector< const CaloCell * > & seedCells ) const
private

Get the hottest neighbour cell in the phi direction.

Definition at line 286 of file TauShotFinder.cxx.

287 {
288
289 const IdentifierHash seedHash = seedCell.caloDDE()->calo_hash();
290
291 // Obtain the neighbour cells in the phi direction
292 std::vector<const CaloCell*> neighCells;
293 for (const CaloCell* neighCell : seedCells) {
294 if (neighCell == &seedCell) continue;
295
296 IdentifierHash neighHash = neighCell->caloDDE()->calo_hash();
297 if (this->isPhiNeighbour(seedHash, neighHash)) {
298 neighCells.push_back(neighCell);
299 }
300 }
301 std::sort(neighCells.begin(),neighCells.end(),ptSort(*this));
302
303 // Select the one with largest pt
304 const CaloCell* phiNeigCell = nullptr;
305 if (!neighCells.empty()) {
306 phiNeigCell = neighCells[0];
307 }
308
309 return phiNeigCell;
310}
const CaloDetDescrElement * caloDDE() const
get pointer to CaloDetDescrElement (data member)
Definition CaloCell.h:321
bool isPhiNeighbour(IdentifierHash cell1Hash, IdentifierHash cell2Hash) const
Check whether two cells are neighbours in the phi direction.
void sort(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end)
Specialization of sort for DataVector/List.

◆ 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 TauShotFinder::initialize ( void )
overridevirtual

Tool initializer.

Reimplemented from TauRecToolBase.

Definition at line 19 of file TauShotFinder.cxx.

19 {
20
21 ATH_CHECK(m_caloWeightTool.retrieve());
24 ATH_CHECK(detStore()->retrieve (m_calo_id, "CaloCell_ID"));
25 ATH_MSG_INFO("Find TauShot in context: " << (inEleRM() ? "`EleRM`" : "`Standard`") << ", with Electron cell removal Flag: " << m_removeElectronCells);
26 return StatusCode::SUCCESS;
27}
#define ATH_MSG_INFO(x)
const ServiceHandle< StoreGateSvc > & detStore() const
bool inEleRM() const
SG::ReadHandleKey< xAOD::CaloClusterContainer > m_removedClusterInputContainer
Gaudi::Property< bool > m_removeElectronCells

◆ 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

◆ isPhiNeighbour()

bool TauShotFinder::isPhiNeighbour ( IdentifierHash cell1Hash,
IdentifierHash cell2Hash ) const
private

Check whether two cells are neighbours in the phi direction.

Definition at line 264 of file TauShotFinder.cxx.

264 {
265 std::vector<IdentifierHash> neigHashes;
266
267 // Next cell in phi direction
268 m_calo_id->get_neighbours(cell1Hash,LArNeighbours::nextInPhi,neigHashes);
269 if (neigHashes.size() > 1) {
270 ATH_MSG_DEBUG(cell1Hash << " has " << neigHashes.size() << " neighbours in the next phi direction !");
271 }
272 if (std::find(neigHashes.begin(), neigHashes.end(), cell2Hash) != neigHashes.end()) {
273 return true;
274 }
275
276 // Previous cell in phi direction
277 m_calo_id->get_neighbours(cell1Hash,LArNeighbours::prevInPhi,neigHashes);
278 if (neigHashes.size() > 1) {
279 ATH_MSG_DEBUG(cell1Hash << " has " << neigHashes.size() << " neighbours in the previous phi direction !");
280 }
281 return std::find(neigHashes.begin(), neigHashes.end(), cell2Hash) != neigHashes.end();
282}

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

◆ selectCells()

StatusCode TauShotFinder::selectCells ( const xAOD::TauJet & tau,
std::vector< const CaloCell * > & cells ) const
private

Apply preselection of the cells Cells within dR < 0.4, in EM1, and pt > 100 MeV are selected.

Definition at line 173 of file TauShotFinder.cxx.

174 {
175 // if in EleRM tau reco, do electron cell removal
176 std::vector<const CaloCell*> removed_cells;
178 SG::ReadHandle<xAOD::CaloClusterContainer> removedClustersHandle( m_removedClusterInputContainer );
179 if (!removedClustersHandle.isValid()){
180 ATH_MSG_ERROR ("Could not retrieve HiveDataObj with key " << removedClustersHandle.key());
181 return StatusCode::FAILURE;
182 }
183 const xAOD::CaloClusterContainer *removed_clusters_cont = removedClustersHandle.cptr();
184
185 for (auto cluster : *removed_clusters_cont){
186 for(auto cell_it = cluster->cell_cbegin(); cell_it != cluster->cell_cend(); cell_it++){
187 removed_cells.push_back(*cell_it);
188 }
189 }
190 }
191
192 // retrieve EM1 cells within dR=0.4, pre-selected by TauPi0CreateROI
193 static const SG::ConstAccessor<std::vector<const CaloCell*>> acc_shotCells("shotCells");
194 std::vector<const CaloCell*> shotCells = acc_shotCells(tau);
195
196 for (const CaloCell* cell : shotCells) {
197 // Require cells above threshold (100 MeV by default)
198 // FIXME: cells are not corrected to point at tau vertex
199 if (cell->pt() * m_caloWeightTool->wtCell(cell) < m_energyThreshold) continue;
200 // if in EleRM, check the clusters do not include electron activities
201 if (m_removeElectronCells && inEleRM() && std::find(removed_cells.cbegin(), removed_cells.cend(), cell) != removed_cells.cend()) continue;
202
203 cells.push_back(cell);
204 }
205 return StatusCode::SUCCESS;
206}
Gaudi::Property< float > m_energyThreshold
CaloClusterContainer_v1 CaloClusterContainer
Define the latest version of the calorimeter cluster container.

◆ selectSeedCells()

StatusCode TauShotFinder::selectSeedCells ( const xAOD::TauJet & tau,
const CaloCellContainer & cellContainer,
std::vector< const CaloCell * > & seedCells ) const
private

Select the seed cells used to construct the shot Cells must sastisfy:

  1. pre-selction: dR < 0.4, in EM1, and pt > 100 MeV
  2. have largest pt among the neighbours in the eta direction
  3. no other seed cells as neighbors in the eta direction

Definition at line 210 of file TauShotFinder.cxx.

212 {
213
214 // Apply pre-selection of the cells
215 assert(seedCells.empty());
216 std::vector<const CaloCell*> cells;
217 ATH_CHECK(selectCells(tau, cells));
218 std::sort(cells.begin(),cells.end(),ptSort(*this));
219
220 std::set<IdentifierHash> seedCellHashes;
221
222 // Loop the pt sorted cells, and select the seed cells
223 for (const CaloCell* cell: cells) {
224 const IdentifierHash cellHash = cell->caloDDE()->calo_hash();
225
226 std::vector<IdentifierHash> nextEtaHashes;
227 m_calo_id->get_neighbours(cellHash, LArNeighbours::nextInEta, nextEtaHashes);
228 std::vector<IdentifierHash> prevEtaHashes;
229 m_calo_id->get_neighbours(cellHash, LArNeighbours::prevInEta, prevEtaHashes);
230
231 std::vector<IdentifierHash> neighHashes = nextEtaHashes;
232 neighHashes.insert(neighHashes.end(),prevEtaHashes.begin(),prevEtaHashes.end());
233
234 // Check whether it is a seed cell
235 bool status = true;
236 for (const IdentifierHash& neighHash : neighHashes) {
237 // Seed cells must not have seed cells as neighbours
238 // TODO: maybe this requirement can be removed
239 if (seedCellHashes.find(neighHash) != seedCellHashes.end()) {
240 status = false;
241 break;
242 }
243
244 // Pt of seed cells must be larger than neighbours'
245 const CaloCell* neighCell = cellContainer.findCell(neighHash);
246 if (!neighCell) continue;
247 if (neighCell->pt() * m_caloWeightTool->wtCell(neighCell) >= cell->pt() * m_caloWeightTool->wtCell(cell)) {
248 status = false;
249 break;
250 }
251 } // End of the loop of neighbour cells
252
253 if (!status) continue;
254
255 seedCells.push_back(cell);
256 seedCellHashes.insert(cellHash);
257 } // End of the loop of cells
258
259 return StatusCode::SUCCESS;
260}
StatusCode selectCells(const xAOD::TauJet &tau, std::vector< const CaloCell * > &cells) const
Apply preselection of the cells Cells within dR < 0.4, in EM1, and pt > 100 MeV are selected.
status
Definition merge.py:16

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

const CaloCell_ID* TauShotFinder::m_calo_id = nullptr
private

calo cell navigation

Definition at line 111 of file TauShotFinder.h.

◆ m_caloCellInputContainer

SG::ReadHandleKey<CaloCellContainer> TauShotFinder::m_caloCellInputContainer {this,"Key_caloCellInputContainer", "AllCalo", "input calo cell container key"}
private

Definition at line 106 of file TauShotFinder.h.

106{this,"Key_caloCellInputContainer", "AllCalo", "input calo cell container key"};

◆ m_caloWeightTool

ToolHandle<IHadronicCalibrationTool> TauShotFinder::m_caloWeightTool {this, "CaloWeightTool", "H1WeightToolCSC12Generic"}
private

Definition at line 107 of file TauShotFinder.h.

107{this, "CaloWeightTool", "H1WeightToolCSC12Generic"};

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

Gaudi::Property<std::vector<float> > TauShotFinder::m_doubleShotCut {this, "AutoDoubleShotCut"}
private

Definition at line 101 of file TauShotFinder.h.

101{this, "AutoDoubleShotCut"};

◆ m_dRThreshold

Gaudi::Property<float> TauShotFinder::m_dRThreshold {this, "DRThreshold", 0.4, "dR(cell,tau) to select cells"}
private

Definition at line 103 of file TauShotFinder.h.

103{this, "DRThreshold", 0.4, "dR(cell,tau) to select cells"};

◆ m_energyThreshold

Gaudi::Property<float> TauShotFinder::m_energyThreshold {this, "EnergyThreshold", 100., "energy threshold (in MeV) to select cells"}
private

Definition at line 104 of file TauShotFinder.h.

104{this, "EnergyThreshold", 100., "energy threshold (in MeV) to select cells"};

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

Gaudi::Property<std::vector<float> > TauShotFinder::m_minPtCut {this, "MinPtCut"}
private

Definition at line 100 of file TauShotFinder.h.

100{this, "MinPtCut"};

◆ m_nCellsInEta

Gaudi::Property<int> TauShotFinder::m_nCellsInEta {this, "NCellsInEta"}
private

Definition at line 99 of file TauShotFinder.h.

99{this, "NCellsInEta"};

◆ m_removedClusterInputContainer

SG::ReadHandleKey<xAOD::CaloClusterContainer> TauShotFinder::m_removedClusterInputContainer {this,"Key_RemovedClusterInputContainer", "", "input removed cluster key"}
private

Definition at line 108 of file TauShotFinder.h.

108{this,"Key_RemovedClusterInputContainer", "", "input removed cluster key"};

◆ m_removeElectronCells

Gaudi::Property<bool> TauShotFinder::m_removeElectronCells {this, "RemoveElectronCells", false}
private

Definition at line 102 of file TauShotFinder.h.

102{this, "RemoveElectronCells", false};

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