28 using HitIds = std::vector<size_t>;
30 const auto& clusters = ctx.
eventStore.
get<Clusters>(m_cfg.inputClusters);
31 const auto& hitParticlesMap =
32 ctx.
eventStore.
get<HitParticlesMap>(m_cfg.inputHitParticlesMap);
33 const auto& hitIds = ctx.
eventStore.
get<HitIds>(m_cfg.inputHitIds);
37 size_t hitIdEnd = m_cfg.hitIdStart + m_cfg.hitIdLength;
38 for (
size_t ihit = m_cfg.hitIdStart;
ihit < hitIdEnd; ++
ihit) {
39 auto hitId = hitIds[
ihit];
40 auto ic = clusters.nth(
ihit);
41 if (ic == clusters.end()) {
46 ACTS_INFO(
" Cluster " <<
ihit <<
" hitId " << hitId <<
" geoId " << geoId
49 for (
const auto&
p :
makeRange(hitParticlesMap.equal_range(
ihit))) {
55 auto numVolume =
selectVolume(clusters, m_cfg.volumeId).size();
56 auto numLayer =
selectLayer(clusters, m_cfg.volumeId, m_cfg.layerId).size();
58 selectModule(clusters, m_cfg.volumeId, m_cfg.layerId, m_cfg.moduleId);
60 ACTS_INFO(
"Hits total: " << clusters.size());
61 ACTS_INFO(
"Hits in volume " << m_cfg.volumeId <<
": " << numVolume);
62 ACTS_INFO(
"Hits in volume " << m_cfg.volumeId <<
" layer " << m_cfg.layerId
64 ACTS_INFO(
"Hits in volume " << m_cfg.volumeId <<
" layer " << m_cfg.layerId
65 <<
" module " << m_cfg.moduleId <<
": "
66 << rangeModule.size());
70 for (
auto ic = rangeModule.begin(); ic != rangeModule.end(); ++ic) {
71 auto ihit = clusters.index_of(ic);
72 auto hitId = hitIds[
ihit];
76 ACTS_INFO(
" Cluster " <<
ihit <<
" hitId " << hitId <<
" geoId " << geoId