TOPASICPedestals.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <TObject.h>
#include <TProfile.h>
#include <vector>
 
namespace Belle2 {
  class TOPASICPedestals : public TObject {
  public:
 
    enum {c_WindowSize = 64, 
          c_Bits       = 9   
         };
 
    TOPASICPedestals(): m_asicWindow(0), m_offset(0)
    {
      for (unsigned i = 0; i < c_WindowSize; i++) m_pedestals[i] = 0;
    }
 
    explicit TOPASICPedestals(unsigned short asicWindow): m_asicWindow(asicWindow),
      m_offset(0)
    {
      for (unsigned i = 0; i < c_WindowSize; i++) m_pedestals[i] = 0;
    }
 
    int setPedestals(const TProfile* profile, double average = 0);
 
    unsigned getASICWindow() const {return m_asicWindow;}
 
    unsigned getSize() const {return c_WindowSize;}
 
    float getValue(unsigned i) const
    {
      if (i < c_WindowSize) {
        unsigned short mask = (1 << c_Bits);
        mask--;
        return (m_pedestals[i] & mask) + m_offset;
      }
      return 0;
    }
 
    float getError(unsigned i) const
    {
      if (i < c_WindowSize) {
        return (m_pedestals[i] >> c_Bits);
      }
      return 0;
    }
 
    bool isValid(unsigned i) const
    {
      if (i < c_WindowSize) {
        return (m_pedestals[i] != 0);
      }
      return false;
    }
 
 
    unsigned getNumofUnvalid() const
    {
      unsigned bad = 0;
      for (int i = 0; i < c_WindowSize; i++) {
        if (!isValid(i)) bad++;
      }
      return bad;
    }
 
  private:
 
    unsigned getOptimizedOffset(const std::vector<unsigned>& values,
                                const std::vector<unsigned>& errors,
                                unsigned maxDif,
                                unsigned maxErr);
 
    unsigned short m_asicWindow;  
    unsigned short m_offset;      
    unsigned short m_pedestals[c_WindowSize]; 
    ClassDef(TOPASICPedestals, 1); 
  };
 
} // end namespace Belle2