TOPSampleTimes.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Marko Staric                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <TObject.h>
#include <vector>
#include <top/dbobjects/TOPASICPedestals.h>
 
namespace Belle2 {
  class TOPSampleTimes: public TObject {
  public:
 
    enum EStatus {
      c_Default = 0,    
      c_Calibrated = 1, 
      c_Unusable = 2    
    };
 
    enum {c_WindowSize = TOPASICPedestals::c_WindowSize,
          c_TimeAxisSize = c_WindowSize * 4
         };
 
    TOPSampleTimes()
    {}
 
    TOPSampleTimes(unsigned scrodID, unsigned channel, double syncTimeBase):
      m_scrodID(scrodID), m_channel(channel)
    {
      setTimeAxis(syncTimeBase);
    }
 
    ~TOPSampleTimes()
    {
    }
 
    void setTimeAxis(double syncTimeBase);
 
    void setTimeAxis(const std::vector<double>& sampleTimes, double syncTimeBase);
 
    void setTimeError(double error) {m_timeError = error;}
 
    void setUnusable() {m_calibrated = c_Unusable;}
 
    unsigned getScrodID() const {return m_scrodID;}
 
    unsigned getChannel() const {return m_channel;}
 
    double getTimeRange() const {return m_timeAxis[c_TimeAxisSize];}
 
    std::vector<double> getTimeAxis() const;
 
    double getTimeError() const {return m_timeError;}
 
    double getFullTime(int window, double sample) const;
 
    double getTime(int window, double sample) const
    {
      return getFullTime(window, sample) - window * getTimeRange() / 4.0;
    }
 
    double getDeltaTime(int window, double sample2, double sample1) const
    {
      return getFullTime(window, sample2) - getFullTime(window, sample1);
    }
 
    double getTimeBin(int window, int sampleNumber) const;
 
    double getSample(int window, double time) const;
 
    bool isCalibrated() const {return m_calibrated == c_Calibrated;}
 
    bool isDefault() const {return m_calibrated == c_Default;}
 
    bool isUnusable() const {return m_calibrated == c_Unusable;}
 
  private:
 
    unsigned short m_scrodID = 0;          
    unsigned short m_channel = 0;          
    float m_timeAxis[c_TimeAxisSize + 1] = {0};  
    float m_timeError = 0;                 
    EStatus m_calibrated = c_Default;      
    ClassDef(TOPSampleTimes, 2); 
  };
 
} // end namespace Belle2