ATLAS Offline Software
Public Member Functions | Private Member Functions | Private Attributes | List of all members
PFClusterSelectorTool Class Reference

#include <PFClusterSelectorTool.h>

Inheritance diagram for PFClusterSelectorTool:
Collaboration diagram for PFClusterSelectorTool:

Public Member Functions

 PFClusterSelectorTool (const std::string &type, const std::string &name, const IInterface *parent)
 Default constructor. More...
 
 ~PFClusterSelectorTool ()
 Default destructor. More...
 
StatusCode initialize ()
 Gaudi AthAlgorithm hooks. More...
 
StatusCode execute (eflowRecClusterContainer &theEFlowRecClusterContainer, xAOD::CaloClusterContainer &theCaloClusterContainer) const
 
StatusCode finalize ()
 

Private Member Functions

void retrieveLCCalCellWeight (const double &energy, const unsigned &index, std::map< IdentifierHash, double > &cellsWeight, const xAOD::CaloClusterContainer &caloCalClustersContainer, const CaloDetDescrManager &calo_dd_man) const
 for EM mode, LC weight for cells are retrieved before doing any subtraction; they will be used after subtraction More...
 

Private Attributes

SG::ReadHandleKey< xAOD::CaloClusterContainer > m_caloClustersReadHandleKey
 ReadHandleKey for the CaloClusterContainer to be used as input. More...
 
SG::ReadHandleKey< xAOD::CaloClusterContainer > m_caloCalClustersReadHandleKey
 ReadHandleKey for the CaloClusterContainer, at LC scale, to be used as input. More...
 
SG::ReadDecorHandleKey< xAOD::CaloClusterContainer > m_caloClusterReadDecorHandleKeyNLeadingTruthParticles {this,"CaloClusterReadDecorHandleKey_NLeadingTruthParticles",""}
 Read handle key to decorate CaloCluster with threeN leading truth particle barcode and energy. More...
 
SG::ReadCondHandleKey< CaloDetDescrManager > m_caloMgrKey
 ReadCondHandleKey for CaloDetDescrManager. More...
 

Detailed Description

Definition at line 16 of file PFClusterSelectorTool.h.

Constructor & Destructor Documentation

◆ PFClusterSelectorTool()

PFClusterSelectorTool::PFClusterSelectorTool ( const std::string &  type,
const std::string &  name,
const IInterface *  parent 
)

Default constructor.

Definition at line 10 of file PFClusterSelectorTool.cxx.

10  :
11  base_class(type, name, parent)
12 {
13 }

◆ ~PFClusterSelectorTool()

PFClusterSelectorTool::~PFClusterSelectorTool ( )
inline

Default destructor.

Definition at line 22 of file PFClusterSelectorTool.h.

22 {};

Member Function Documentation

◆ execute()

StatusCode PFClusterSelectorTool::execute ( eflowRecClusterContainer &  theEFlowRecClusterContainer,
xAOD::CaloClusterContainer &  theCaloClusterContainer 
) const

Definition at line 32 of file PFClusterSelectorTool.cxx.

35 {
36 
37  SG::ReadHandle<xAOD::CaloClusterContainer> caloClustersReadHandle(m_caloClustersReadHandleKey);
38 
39  /* Verify the read handle has a valid pointer, and if not return */
40  if (!caloClustersReadHandle.isValid()){
41  ATH_MSG_WARNING(" Invalid ReadHandle for xAOD::CaloCluster with key: " << caloClustersReadHandle.key());
42  return StatusCode::SUCCESS;
43  }
44 
45  const xAOD::CaloClusterContainer* calclusters = nullptr;
46  if (!m_caloCalClustersReadHandleKey.key().empty()) {
47  SG::ReadHandle<xAOD::CaloClusterContainer> caloCalClustersReadHandle(m_caloCalClustersReadHandleKey);
48  calclusters = caloCalClustersReadHandle.get();
49  }
50 
51  SG::ReadCondHandle<CaloDetDescrManager> caloMgrHandle{m_caloMgrKey};
52  ATH_CHECK(caloMgrHandle.isValid());
53 
54  /* Fill the vector of eflowRecClusters */
55  unsigned int nClusters = caloClustersReadHandle->size();
56  for (unsigned int iCluster = 0; iCluster < nClusters; ++iCluster) {
57  /* Create the eflowRecCluster and put it in the container */
58  std::unique_ptr<eflowRecCluster> thisEFRecCluster = std::make_unique<eflowRecCluster>(ElementLink<xAOD::CaloClusterContainer>(*caloClustersReadHandle, iCluster), theCaloClusterContainer);
59 
60  if (calclusters){
61  std::map<IdentifierHash,double> cellsWeightMap;
62  retrieveLCCalCellWeight(caloClustersReadHandle->at(iCluster)->e(), iCluster, cellsWeightMap, *calclusters,**caloMgrHandle);
63 
64  thisEFRecCluster->setCellsWeight(std::move(cellsWeightMap));
65  }
66 
67  if (!m_caloClusterReadDecorHandleKeyNLeadingTruthParticles.empty()){
68 
69  std::string decorHandleName = m_caloClusterReadDecorHandleKeyNLeadingTruthParticles.key();
70  std::string::size_type pos = decorHandleName.find(".");
71  std::string decorName = decorHandleName.substr(pos+1);
72 
73  SG::AuxElement::Accessor< std::vector< std::pair<unsigned int, double> > > accessor(decorName);
74  accessor(*theCaloClusterContainer[iCluster]) = accessor(*caloClustersReadHandle->at(iCluster));
75  }
76 
77  thisEFRecCluster->setClusterId(iCluster);
78  theEFlowRecClusterContainer.push_back(std::move(thisEFRecCluster));
79  }
80 
81  return StatusCode::SUCCESS;
82 }

◆ finalize()

StatusCode PFClusterSelectorTool::finalize ( )

Definition at line 84 of file PFClusterSelectorTool.cxx.

84  {
85  return StatusCode::SUCCESS;
86 }

◆ initialize()

StatusCode PFClusterSelectorTool::initialize ( )

Gaudi AthAlgorithm hooks.

Definition at line 15 of file PFClusterSelectorTool.cxx.

15  {
16 
17  ATH_CHECK(m_caloClustersReadHandleKey.initialize());
18  if(!m_caloCalClustersReadHandleKey.key().empty()) {
19  ATH_CHECK(m_caloCalClustersReadHandleKey.initialize());
20  }
21 
22  if (!m_caloClusterReadDecorHandleKeyNLeadingTruthParticles.empty()){
23  ATH_CHECK(m_caloClusterReadDecorHandleKeyNLeadingTruthParticles.initialize());
24  }
25 
26  ATH_CHECK(m_caloMgrKey.initialize());
27 
28  return StatusCode::SUCCESS;
29 }

◆ retrieveLCCalCellWeight()

void PFClusterSelectorTool::retrieveLCCalCellWeight ( const double &  energy,
const unsigned &  index,
std::map< IdentifierHash, double > &  cellsWeight,
const xAOD::CaloClusterContainer &  caloCalClustersContainer,
const CaloDetDescrManager &  calo_dd_man 
) const
private

for EM mode, LC weight for cells are retrieved before doing any subtraction; they will be used after subtraction

Definition at line 89 of file PFClusterSelectorTool.cxx.

95 {
96  /* match CaloCluster with CaloCalCluster to obtain cell weight */
97  /* first try the position at 'index'. If we are lucky, the loop can be
98  * avoided. */
99  /* Note the read handle has been tested to be valid prior to the call of this
100  * function */
101  const xAOD::CaloCluster* matchedCalCluster = caloCalClustersContainer.at(index);
102  if (matchedCalCluster){
103  if (!(fabs(energy - matchedCalCluster->rawE()) < 0.001)) {
104  matchedCalCluster = nullptr;
105  for (const auto *iCalCalCluster : caloCalClustersContainer) {
106  matchedCalCluster = iCalCalCluster;
107  if (fabs(energy - matchedCalCluster->rawE()) < 0.001) {
108  break;
109  }
110  }
111  if (!matchedCalCluster) ATH_MSG_WARNING("Invalid pointer to matched cluster - failed to find cluster match");
112  }
113  if (not matchedCalCluster){
114  throw std::runtime_error("matchedCluster is a null pointer in PFClusterSelectorTool::retrieveLCCalCellWeight");
115  }
116  /* obtain cell index and cell weight */
117  const CaloCell_ID* calo_id = calo_dd_man.getCaloCell_ID();
118  xAOD::CaloCluster::const_cell_iterator itCell = matchedCalCluster->cell_begin();
119  xAOD::CaloCluster::const_cell_iterator endCell = matchedCalCluster->cell_end();
120  for (; itCell != endCell; ++itCell) {
121  const CaloCell* pCell = *itCell;
122  Identifier myId = pCell->ID();
123  IdentifierHash myHashId = calo_id->calo_cell_hash(myId);
124  cellsWeight[myHashId] = itCell.weight();
125  }
126  } else ATH_MSG_WARNING("Invalid pointer to matched cluster - could not look up local hadron cell weights");
127 }

Member Data Documentation

◆ m_caloCalClustersReadHandleKey

SG::ReadHandleKey<xAOD::CaloClusterContainer> PFClusterSelectorTool::m_caloCalClustersReadHandleKey
private
Initial value:
{
this,
"calClustersName",
"CaloCalTopoClusters",
"ReadHandleKey for the CaloClusterContainer, at LC scale, to be used as "
"input"
}

ReadHandleKey for the CaloClusterContainer, at LC scale, to be used as input.

Definition at line 50 of file PFClusterSelectorTool.h.

◆ m_caloClusterReadDecorHandleKeyNLeadingTruthParticles

SG::ReadDecorHandleKey<xAOD::CaloClusterContainer> PFClusterSelectorTool::m_caloClusterReadDecorHandleKeyNLeadingTruthParticles {this,"CaloClusterReadDecorHandleKey_NLeadingTruthParticles",""}
private

Read handle key to decorate CaloCluster with threeN leading truth particle barcode and energy.

Definition at line 59 of file PFClusterSelectorTool.h.

◆ m_caloClustersReadHandleKey

SG::ReadHandleKey<xAOD::CaloClusterContainer> PFClusterSelectorTool::m_caloClustersReadHandleKey
private
Initial value:
{
this,
"clustersName",
"CaloTopoClusters",
"ReadHandleKey for the CaloClusterContainer to be used as input"
}

ReadHandleKey for the CaloClusterContainer to be used as input.

Definition at line 41 of file PFClusterSelectorTool.h.

◆ m_caloMgrKey

SG::ReadCondHandleKey<CaloDetDescrManager> PFClusterSelectorTool::m_caloMgrKey
private
Initial value:
{
this,
"CaloDetDescrManager",
"CaloDetDescrManager",
"SG Key for CaloDetDescrManager in the Condition Store"
}

ReadCondHandleKey for CaloDetDescrManager.

Definition at line 62 of file PFClusterSelectorTool.h.

The documentation for this class was generated from the following files:
xAOD::CaloCluster_v1::rawE
flt_t rawE() const
PFClusterSelectorTool::m_caloClustersReadHandleKey
SG::ReadHandleKey< xAOD::CaloClusterContainer > m_caloClustersReadHandleKey
ReadHandleKey for the CaloClusterContainer to be used as input.
Definition: PFClusterSelectorTool.h:41
CaloCell_Base_ID::calo_cell_hash
IdentifierHash calo_cell_hash(const Identifier cellId) const
create hash id from 'global' cell id
SG::ReadCondHandle
Definition: ReadCondHandle.h:44
SG::Accessor
Helper class to provide type-safe access to aux data.
Definition: Control/AthContainers/AthContainers/Accessor.h:68
SG::ReadHandle
Definition: StoreGate/StoreGate/ReadHandle.h:70
index
Definition: index.py:1
DataVector::get
const T * get(size_type n) const
Access an element, as an rvalue.
SG::VarHandleKey::key
const std::string & key() const
Return the StoreGate ID for the referenced object.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:141
CaloClusterCellLink::const_iterator::weight
weight_t weight() const
Accessor for weight associated to this cell.
Definition: CaloClusterCellLink.h:82
xAOD::CaloCluster_v1
Description of a calorimeter cluster.
Definition: CaloCluster_v1.h:59
ParticleGun_FastCalo_ChargeFlip_Config.energy
energy
Definition: ParticleGun_FastCalo_ChargeFlip_Config.py:78
test_pyathena.parent
parent
Definition: test_pyathena.py:15
eflowRecCluster::setCellsWeight
void setCellsWeight(std::map< IdentifierHash, double > &&cellsWeight)
Definition: eflowRecCluster.h:81
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
CaloCell_ID
Helper class for offline cell identifiers.
Definition: CaloCell_ID.h:34
SG::VarHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:103
DataVector
Derived DataVector<T>.
Definition: DataVector.h:581
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:221
PFClusterSelectorTool::retrieveLCCalCellWeight
void retrieveLCCalCellWeight(const double &energy, const unsigned &index, std::map< IdentifierHash, double > &cellsWeight, const xAOD::CaloClusterContainer &caloCalClustersContainer, const CaloDetDescrManager &calo_dd_man) const
for EM mode, LC weight for cells are retrieved before doing any subtraction; they will be used after ...
Definition: PFClusterSelectorTool.cxx:89
DataVector::push_back
value_type push_back(value_type pElem)
Add an element to the end of the collection.
SG::CondHandleKey::initialize
StatusCode initialize(bool used=true)
CaloCell::ID
Identifier ID() const
get ID (from cached data member) non-virtual and inline for fast access
Definition: CaloCell.h:279
xAOD::JetAttributeAccessor::accessor
const AccessorWrapper< T > * accessor(xAOD::JetAttribute::AttributeID id)
Returns an attribute accessor corresponding to an AttributeID.
Definition: JetAccessorMap.h:26
python.LumiBlobConversion.pos
pos
Definition: LumiBlobConversion.py:18
CaloCell
Data object for each calorimeter readout cell.
Definition: CaloCell.h:57
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
CaloDetDescrManager::getCaloCell_ID
const CaloCell_ID * getCaloCell_ID() const
get calo cell ID helper
Definition: CaloDetDescrManager.cxx:1590
PFClusterSelectorTool::m_caloMgrKey
SG::ReadCondHandleKey< CaloDetDescrManager > m_caloMgrKey
ReadCondHandleKey for CaloDetDescrManager.
Definition: PFClusterSelectorTool.h:62
DataVector::at
const T * at(size_type n) const
Access an element, as an rvalue.
PFClusterSelectorTool::m_caloClusterReadDecorHandleKeyNLeadingTruthParticles
SG::ReadDecorHandleKey< xAOD::CaloClusterContainer > m_caloClusterReadDecorHandleKeyNLeadingTruthParticles
Read handle key to decorate CaloCluster with threeN leading truth particle barcode and energy.
Definition: PFClusterSelectorTool.h:59
PFClusterSelectorTool::m_caloCalClustersReadHandleKey
SG::ReadHandleKey< xAOD::CaloClusterContainer > m_caloCalClustersReadHandleKey
ReadHandleKey for the CaloClusterContainer, at LC scale, to be used as input.
Definition: PFClusterSelectorTool.h:50
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
CaloClusterCellLink::const_iterator
const_iterator to loop over cells belonging to a cluster
Definition: CaloClusterCellLink.h:57
eflowRecCluster::setClusterId
void setClusterId(int clusterId)
Definition: eflowRecCluster.h:80
Identifier
Definition: IdentifierFieldParser.cxx:14
Generated on Thu Nov 7 2024 21:53:54 for ATLAS Offline Software by   doxygen 1.8.18