PedestalProcessorData.h

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
#ifndef EMCAL_PEDESTAL_PROCESSOR_DATA_H_
#define EMCAL_PEDESTAL_PROCESSOR_DATA_H_
#include <array>
#include <exception>
#include <iosfwd>
#include <string>
#include <tuple>
#include "MathUtils/Utils.h"
 
#include "Rtypes.h"
 
namespace o2::emcal
{
 
class PedestalProcessorData
{
 private:
 public:
  class ChannelIndexException : public std::exception
  {
   private:
    unsigned short mChannelIndex; 
    unsigned short mMaxChannels;  
    std::string mErrorMessage;    
 
   public:
    ChannelIndexException(unsigned short channelIndex, unsigned short maxChannels) : std::exception(), mChannelIndex(channelIndex), mMaxChannels(maxChannels)
    {
      mErrorMessage = "Channel index " + std::to_string(mChannelIndex) + " not found (max " + std::to_string(mMaxChannels) + ")";
    }
 
    ~ChannelIndexException() noexcept final = default;
 
    const char* what() const noexcept final { return mErrorMessage.data(); }
 
    unsigned short getChannelIndex() const noexcept { return mChannelIndex; }
 
    unsigned short getMaxChannels() const noexcept { return mMaxChannels; }
  };
 
  using ProfileHistFEC = std::array<o2::math_utils::StatAccumulator, 17664>;
  using ProfileHistLEDMON = std::array<o2::math_utils::StatAccumulator, 480>;
  using PedestalValue = std::tuple<double, double>;
 
  PedestalProcessorData() = default;
 
  ~PedestalProcessorData() = default;
 
  PedestalProcessorData& operator+=(const PedestalProcessorData& other);
 
  void fillADC(unsigned short adc, unsigned short tower, bool lowGain, bool LEDMON);
 
  PedestalValue getValue(unsigned short tower, bool lowGain, bool LEDMON) const;
 
  int getEntriesForChannel(unsigned short tower, bool lowGain, bool LEDMON) const;
 
  void reset();
 
  const ProfileHistFEC& accessFECData(bool lowGain) const { return lowGain ? mDataFECLG : mDataFECHG; }
 
  const ProfileHistLEDMON& accessLEDMONData(bool lowGain) const { return lowGain ? mDataLEDMONLG : mDataLEDMONHG; }
 
 private:
  ProfileHistFEC mDataFECHG;       
  ProfileHistFEC mDataFECLG;       
  ProfileHistLEDMON mDataLEDMONHG; 
  ProfileHistLEDMON mDataLEDMONLG; 
 
  ClassDefNV(PedestalProcessorData, 1);
};
 
PedestalProcessorData operator+(const PedestalProcessorData& lhs, const PedestalProcessorData& rhs);
 
std::ostream& operator<<(std::ostream& stream, const PedestalProcessorData::ChannelIndexException& ex);
 
} // namespace o2::emcal
 
#endif