ATLAS Offline Software
Loading...
Searching...
No Matches
Classes | Public Types | Public Member Functions | Private Member Functions | Private Attributes | List of all members
ActsTrk::HgtdTimedClusteringTool Class Reference

#include <HgtdTimedClusteringTool.h>

Inheritance diagram for ActsTrk::HgtdTimedClusteringTool:
Collaboration diagram for ActsTrk::HgtdTimedClusteringTool:

Classes

struct  Cluster

Public Types

using Cell = Hgtd::UnpackedHgtdRDO
using CellCollection = std::vector<Cell>
using ClusterCollection = std::vector<Cluster>

Public Member Functions

 HgtdTimedClusteringTool (const std::string &type, const std::string &name, const IInterface *parent)
virtual StatusCode initialize () override
virtual StatusCode clusterize (const EventContext &ctx, const RawDataCollection &RDOs, ClusterContainer &container) const override

Private Member Functions

StatusCode makeCluster (const EventContext &ctx, const typename HgtdTimedClusteringTool::Cluster &cluster, xAOD::HGTDCluster &xaodcluster) const

Private Attributes

const HGTD_DetectorManager * m_hgtd_det_mgr {nullptr}
const HGTD_ID * m_hgtd_id {nullptr}
Gaudi::Property< double > m_timeTollerance {this, "TimeTollerance", 0.035 * Acts::UnitConstants::ns}
Gaudi::Property< bool > m_addCorners {this, "AddCorners", true}

Detailed Description

Definition at line 33 of file HgtdTimedClusteringTool.h.

Member Typedef Documentation

◆ Cell

using ActsTrk::HgtdTimedClusteringTool::Cell = Hgtd::UnpackedHgtdRDO

Definition at line 42 of file HgtdTimedClusteringTool.h.

◆ CellCollection

using ActsTrk::HgtdTimedClusteringTool::CellCollection = std::vector<Cell>

Definition at line 43 of file HgtdTimedClusteringTool.h.

◆ ClusterCollection

using ActsTrk::HgtdTimedClusteringTool::ClusterCollection = std::vector<Cluster>

Definition at line 44 of file HgtdTimedClusteringTool.h.

Constructor & Destructor Documentation

◆ HgtdTimedClusteringTool()

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

Definition at line 31 of file HgtdTimedClusteringTool.cxx.

34 : base_class(type, name, parent)
35 {}

Member Function Documentation

◆ clusterize()

StatusCode ActsTrk::HgtdTimedClusteringTool::clusterize ( const EventContext & ctx,
const RawDataCollection & RDOs,
ClusterContainer & container ) const
overridevirtual

Definition at line 49 of file HgtdTimedClusteringTool.cxx.

52 {
53 // Unpack RDOs (would need a proper function here)
54 CellCollection cells;
55 cells.reserve(RDOs.size());
56 for (const HGTD_RDO* rdo : RDOs) {
57 Identifier id = rdo->identify();
58 cells.emplace_back(-1,
59 m_hgtd_id->phi_index(id),
60 m_hgtd_id->eta_index(id),
61 rdo->getTOA(),
62 rdo->getTOT(),
63 id);
64 }
65
66 ATH_MSG_DEBUG("Clustering on " << RDOs.size() << " RDOs using time information");
67 Acts::Ccl::ClusteringData data;
68 ClusterCollection clusters;
69 Acts::Ccl::createClusters<CellCollection, ClusterCollection, 2>
70 (data, cells, clusters, Acts::Ccl::TimedConnect<Cell, 2ul>(m_timeTollerance.value(), m_addCorners.value()));
71 ATH_MSG_DEBUG(" \\_ " << clusters.size() << " clusters reconstructed");
72
73 // Fast insertion trick
74 std::size_t previousSizeContainer = container.size();
75 std::vector<xAOD::HGTDCluster*> toAdd;
76 toAdd.reserve(clusters.size());
77 for (std::size_t i(0), n(clusters.size()); i < n; ++i)
78 toAdd.push_back( new xAOD::HGTDCluster() );
79 container.insert(container.end(), toAdd.begin(), toAdd.end());
80
81 for (std::size_t i(0); i<clusters.size(); ++i) {
82 const typename HgtdTimedClusteringTool::Cluster& cluster = clusters[i];
83 ATH_CHECK(makeCluster(ctx, cluster, *container[previousSizeContainer+i]));
84 }
85
86 return StatusCode::SUCCESS;
87 }
ATH_CHECK
#define ATH_CHECK
Evaluate an expression and check for errors.
Definition AthCheckMacros.h:40
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition AthMsgStreamMacros.h:29
data
char data[hepevt_bytes_allocation_ATLAS]
Definition HepEvt.cxx:11
ActsTrk::HgtdTimedClusteringTool::m_timeTollerance
Gaudi::Property< double > m_timeTollerance
Definition HgtdTimedClusteringTool.h:67
ActsTrk::HgtdTimedClusteringTool::makeCluster
StatusCode makeCluster(const EventContext &ctx, const typename HgtdTimedClusteringTool::Cluster &cluster, xAOD::HGTDCluster &xaodcluster) const
Definition HgtdTimedClusteringTool.cxx:89
xAOD::HGTDCluster
HGTDCluster_v1 HGTDCluster
Define the version of the pixel cluster class.
Definition HGTDCluster.h:13

◆ initialize()

StatusCode ActsTrk::HgtdTimedClusteringTool::initialize ( )
overridevirtual

Definition at line 37 of file HgtdTimedClusteringTool.cxx.

38 {
39 ATH_MSG_DEBUG("Initializing " << name() << "...");
40
41 ATH_CHECK(detStore()->retrieve(m_hgtd_det_mgr, "HGTD"));
42 ATH_CHECK(detStore()->retrieve(m_hgtd_id, "HGTD_ID"));
43
44 ATH_MSG_DEBUG(m_timeTollerance);
45
46 return StatusCode::SUCCESS;
47 }
ActsTrk::HgtdTimedClusteringTool::m_hgtd_det_mgr
const HGTD_DetectorManager * m_hgtd_det_mgr
Definition HgtdTimedClusteringTool.h:64
python.PyKernel.retrieve
retrieve(aClass, aKey=None)
Definition PyKernel.py:110

◆ makeCluster()

StatusCode ActsTrk::HgtdTimedClusteringTool::makeCluster ( const EventContext & ctx,
const typename HgtdTimedClusteringTool::Cluster & cluster,
xAOD::HGTDCluster & xaodcluster ) const
private

Definition at line 89 of file HgtdTimedClusteringTool.cxx.

92 {
93 if (cluster.ids.empty()) return StatusCode::SUCCESS;
94
95 InDetDD::SiLocalPosition pos_acc(0,0);
96 double tot_time = 0;
97
98 for (size_t i = 0; i < cluster.ids.size(); i++) {
99 Identifier rdo_id = cluster.ids[i];
100
101 const InDetDD::HGTD_DetectorElement* element = m_hgtd_det_mgr->getDetectorElement(rdo_id);
102 InDetDD::SiCellId si_cell_id = element->cellIdFromIdentifier(rdo_id);
103 InDetDD::SiLocalPosition si_pos = element->design().localPositionOfCell(si_cell_id);
104
105 pos_acc += si_pos;
106 tot_time += cluster.times[i];
107 }
108
109 pos_acc /= cluster.ids.size();
110 tot_time /= cluster.ids.size();
111
112 // Create the cluster
113 Eigen::Matrix<float, 3, 1> loc_pos(pos_acc.xPhi(), pos_acc.xEta(), tot_time);
114 Eigen::Matrix<float, 3, 3> cov_matrix= Eigen::Matrix<float, 3, 3>::Zero();
115
116 constexpr float xWidth = 1.3;
117 constexpr float yWidth = 1.3;
118 cov_matrix(0,0) = xWidth * xWidth / 12 / cluster.ids.size(); // i.e. Cov XX
119 cov_matrix(1,1) = yWidth * yWidth / 12 / cluster.ids.size(); // i.e. Cov YY
120 constexpr float time_of_arrival_err = 0.035;
121 cov_matrix(2,2) = time_of_arrival_err * time_of_arrival_err / cluster.ids.size(); // i.e. Cov TT
122
123 IdentifierHash id_hash = m_hgtd_det_mgr->getDetectorElement(cluster.ids.front())->identifyHash();
124
125 // Fill
126 xaodcluster.setMeasurement<3>(id_hash,loc_pos,cov_matrix);
127 xaodcluster.setIdentifier(cluster.ids.front().get_compact());
128 xaodcluster.setRDOlist(std::move(cluster.ids));
129 xaodcluster.setToTlist(std::move(cluster.tots));
130
131 return StatusCode::SUCCESS;
132 }
InDetDD::HGTD_DetectorElement::cellIdFromIdentifier
SiCellId cellIdFromIdentifier(const Identifier &identifier) const override final
SiCellId from Identifier.
Definition HGTD_DetectorElement.cxx:51
InDetDD::HGTD_DetectorElement::design
const HGTD_ModuleDesign & design() const override final
access to the local description:
Definition HGTD_DetectorElement.h:121
InDetDD::HGTD_ModuleDesign::localPositionOfCell
virtual SiLocalPosition localPositionOfCell(const SiCellId &cellId) const
readout or diode id -> position.
Definition HGTD_ModuleDesign.cxx:145
xAOD::HGTDCluster_v1::setToTlist
void setToTlist(const std::vector< int > &tots)
Sets the list of ToT of the channels building the cluster.
xAOD::HGTDCluster_v1::setRDOlist
void setRDOlist(const std::vector< Identifier > &rdolist)
Sets the list of identifiers of the channels building the cluster.
Definition HGTDCluster_v1.cxx:15
xAOD::UncalibratedMeasurement_v1::setMeasurement
void setMeasurement(const DetectorIDHashType idHash, MeasVector< N > locPos, MeasMatrix< N > locCov)
Sets IdentifierHash, local position and local covariance of the measurement.
xAOD::UncalibratedMeasurement_v1::setIdentifier
void setIdentifier(const DetectorIdentType measId)
Sets the full Identifier of the measurement.

Member Data Documentation

◆ m_addCorners

Gaudi::Property<bool> ActsTrk::HgtdTimedClusteringTool::m_addCorners {this, "AddCorners", true}
private

Definition at line 68 of file HgtdTimedClusteringTool.h.

68{this, "AddCorners", true};

◆ m_hgtd_det_mgr

const HGTD_DetectorManager* ActsTrk::HgtdTimedClusteringTool::m_hgtd_det_mgr {nullptr}
private

Definition at line 64 of file HgtdTimedClusteringTool.h.

64{nullptr};

◆ m_hgtd_id

const HGTD_ID* ActsTrk::HgtdTimedClusteringTool::m_hgtd_id {nullptr}
private

Definition at line 65 of file HgtdTimedClusteringTool.h.

65{nullptr};

◆ m_timeTollerance

Gaudi::Property<double> ActsTrk::HgtdTimedClusteringTool::m_timeTollerance {this, "TimeTollerance", 0.035 * Acts::UnitConstants::ns}
private

Definition at line 67 of file HgtdTimedClusteringTool.h.

67{this, "TimeTollerance", 0.035 * Acts::UnitConstants::ns};
The documentation for this class was generated from the following files:
Generated on for ATLAS Offline Software by doxygen 1.14.0