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.

    : base_class(type, name, parent)
  {}

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.

  {
    // Unpack RDOs (would need a proper function here)
    CellCollection cells;
    cells.reserve(RDOs.size());
    for (const HGTD_RDO* rdo : RDOs) {
      Identifier id = rdo->identify();
      cells.emplace_back(-1,
             m_hgtd_id->phi_index(id),
             m_hgtd_id->eta_index(id),
             rdo->getTOA(),
             rdo->getTOT(),
             id);
    }
 
    ATH_MSG_DEBUG("Clustering on " << RDOs.size() << " RDOs using time information");
    ClusterCollection clusters =
      Acts::Ccl::createClusters<CellCollection, ClusterCollection, 2>
      (cells, Acts::Ccl::TimedConnect<Cell, 2ul>(m_timeTollerance.value(), m_addCorners.value()));
    ATH_MSG_DEBUG("   \\_ " << clusters.size() << " clusters reconstructed");
 
    // Fast insertion trick
    std::size_t previousSizeContainer = container.size();
    std::vector<xAOD::HGTDCluster*> toAdd;
    toAdd.reserve(clusters.size());
    for (std::size_t i(0), n(clusters.size()); i < n; ++i)
      toAdd.push_back( new xAOD::HGTDCluster() );
    container.insert(container.end(), toAdd.begin(), toAdd.end());
 
    for (std::size_t i(0); i<clusters.size(); ++i) {
      const typename HgtdTimedClusteringTool::Cluster& cluster = clusters[i];
      ATH_CHECK(makeCluster(ctx, cluster, *container[previousSizeContainer+i]));
    }
 
    return StatusCode::SUCCESS;
  }

initialize()

StatusCode ActsTrk::HgtdTimedClusteringTool::initialize ( )

overridevirtual

Definition at line 37 of file HgtdTimedClusteringTool.cxx.

  {
    ATH_MSG_DEBUG("Initializing " << name() << "...");
 
    ATH_CHECK(detStore()->retrieve(m_hgtd_det_mgr, "HGTD"));
    ATH_CHECK(detStore()->retrieve(m_hgtd_id, "HGTD_ID"));
 
    ATH_MSG_DEBUG(m_timeTollerance);
 
    return StatusCode::SUCCESS;
  }

makeCluster()

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

private

Definition at line 88 of file HgtdTimedClusteringTool.cxx.

  {
    if (cluster.ids.empty()) return StatusCode::SUCCESS;
 
    InDetDD::SiLocalPosition pos_acc(0,0);
    double tot_time = 0;
 
    for (size_t i = 0; i < cluster.ids.size(); i++) {
      Identifier rdo_id = cluster.ids[i];
 
      const InDetDD::HGTD_DetectorElement* element = m_hgtd_det_mgr->getDetectorElement(rdo_id);
      InDetDD::SiCellId si_cell_id = element->cellIdFromIdentifier(rdo_id);
      InDetDD::SiLocalPosition si_pos = element->design().localPositionOfCell(si_cell_id);
 
      pos_acc += si_pos;
      tot_time += cluster.times[i];
    }
 
    pos_acc /= cluster.ids.size();
    tot_time /= cluster.ids.size();
 
    // Create the cluster
    Eigen::Matrix<float, 3, 1> loc_pos(pos_acc.xPhi(), pos_acc.xEta(), tot_time);
    Eigen::Matrix<float, 3, 3> cov_matrix= Eigen::Matrix<float, 3, 3>::Zero();
 
    constexpr float xWidth = 1.3;
    constexpr float yWidth = 1.3;
    cov_matrix(0,0) = xWidth * xWidth / 12 / cluster.ids.size(); // i.e. Cov XX
    cov_matrix(1,1) = yWidth * yWidth / 12 / cluster.ids.size(); // i.e. Cov YY
    constexpr float time_of_arrival_err = 0.035;
    cov_matrix(2,2) = time_of_arrival_err * time_of_arrival_err / cluster.ids.size(); // i.e. Cov TT
 
    IdentifierHash id_hash = m_hgtd_det_mgr->getDetectorElement(cluster.ids.front())->identifyHash();
 
    // Fill
    xaodcluster.setMeasurement<3>(id_hash,loc_pos,cov_matrix);
    xaodcluster.setIdentifier(cluster.ids.front().get_compact());
    xaodcluster.setRDOlist(std::move(cluster.ids));
    xaodcluster.setToTlist(std::move(cluster.tots));
 
    return StatusCode::SUCCESS;
  }

Member Data Documentation

m_addCorners

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

private

Definition at line 68 of file HgtdTimedClusteringTool.h.

m_hgtd_det_mgr

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

private

Definition at line 64 of file HgtdTimedClusteringTool.h.

m_hgtd_id

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

private

Definition at line 65 of file HgtdTimedClusteringTool.h.

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.

---

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

• HgtdTimedClusteringTool.h
• HgtdTimedClusteringTool.cxx