TrackerPedestals.h

#ifndef TRACKING_RECO_TRACKERPEDESTALS_H_
#define TRACKING_RECO_TRACKERPEDESTALS_H_
 
#include <array>
#include <cstdint>
#include <numeric>
#include <string>
#include <unordered_map>
 
#include "Framework/ConditionsObject.h"
#include "Tracking/Reco/SiStripChannelMap.h"
 
namespace tracking::reco {
 
class TrackerPedestals : public framework::ConditionsObject {
 public:
  static const std::string CONDITIONS_NAME;
 
  struct Channel {
    std::array<float, channelmap::K_SAMPLES_PER_APV_TRIGGER> mean_{};
    std::array<float, channelmap::K_SAMPLES_PER_APV_TRIGGER> noise_{};
  };
 
  TrackerPedestals() : framework::ConditionsObject(CONDITIONS_NAME) {}
 
  void add(uint8_t feb, uint8_t hybrid, uint8_t apv, uint8_t channel,
           const Channel& ped) {
    channels_[channelmap::channelId(feb, hybrid, apv, channel)] = ped;
  }
 
  const Channel* find(uint8_t feb, uint8_t hybrid, uint8_t apv,
                      uint8_t channel) const {
    auto it = channels_.find(channelmap::channelId(feb, hybrid, apv, channel));
    return it == channels_.end() ? nullptr : &it->second;
  }
 
  float noise(uint8_t feb, uint8_t hybrid, uint8_t apv, uint8_t channel) const {
    const Channel* c = find(feb, hybrid, apv, channel);
    if (!c) return 0.f;
    const float sum =
        std::accumulate(c->noise_.begin(), c->noise_.end(), 0.f);
    return sum / static_cast<float>(c->noise_.size());
  }
 
  std::size_t size() const { return channels_.size(); }
 
 private:
  std::unordered_map<uint32_t, Channel> channels_;
};
 
}  // namespace tracking::reco
 
#endif  // TRACKING_RECO_TRACKERPEDESTALS_H_