ChannelMerger.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file ChannelMerger.hpp

// This file is part of the Acts project.
//
// Copyright (C) 2020 CERN for the benefit of the Acts project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at http://mozilla.org/MPL/2.0/.

#pragma once

#include "Acts/Surfaces/Surface.hpp"
#include "Acts/Surfaces/SurfaceError.hpp"
#include "ActsFatras/Digitization/DigitizationData.hpp"

#include <array>
#include <utility>

namespace ActsFatras {

template <typename signal_t, size_t kSize>
const std::vector<Channel<signal_t, kSize>> mergeChannels(
    const std::vector<Channel<signal_t, kSize>>& channels) {
  using Channel = Channel<signal_t, kSize>;
  using ChannelKey = std::array<unsigned int, kSize>;

  // Fill a channel map - use the channel identification
  auto extractChannelKey = [&](const Channel& ch) -> ChannelKey {
    ChannelKey cKey;
    for (unsigned int ik = 0; ik < kSize; ++ik) {
      cKey[ik] = ch.cellId[ik].first;
    }
    return cKey;
  };

  std::map<ChannelKey, Channel> channelMap;
  for (const auto& ch : channels) {
    ChannelKey key = extractChannelKey(ch);
    auto chItr = channelMap.find(key);
    if (chItr != channelMap.end()) {
      chItr->second.value += ch.value;
      chItr->second.links.insert(ch.links.begin(), ch.links.end());
    } else {
      channelMap.insert(std::pair<ChannelKey, Channel>(key, ch));
    }
  }
  // Unroll the channels after merging
  std::vector<Channel> mergedChannels;
  mergedChannels.reserve(channelMap.size());
  for (auto& [key, value] : channelMap) {
    mergedChannels.push_back(value);
  }
  return mergedChannels;
}

}  // namespace ActsFatras