TOPAsicShiftsBS13dAlgorithm.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2019 - Belle II Collaboration                             *
 *                                                                        *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Marko Staric                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <calibration/CalibrationAlgorithm.h>
#include <TH1F.h>
#include <vector>
 
namespace Belle2 {
  namespace TOP {
 
    class TOPAsicShiftsBS13dAlgorithm : public CalibrationAlgorithm {
 
    public:
 
      TOPAsicShiftsBS13dAlgorithm();
 
      virtual ~TOPAsicShiftsBS13dAlgorithm() {}
 
      void setMinValue(double minVal) {m_minVal = minVal;}
 
      void setWindowSize(int size) {m_winSize = size;}
 
      void setShiftRange(int begin, int end) {m_shiftBegin = begin; m_shiftEnd = end;}
 
      void setMinimalSignificance(double s) {m_minSignificance = s;}
 
    private:
 
      virtual EResult calibrate() final;
 
      void setWindow();
 
      double fun(int x);
 
      std::vector<double> getPDF(int shift);
 
      double logL(std::shared_ptr<TH1F> h, int shift);
 
      // parameters
      double m_minVal = 1.0; 
      int m_shiftBegin = -20; 
      int m_shiftEnd = 20;  
      double m_minSignificance = 5; 
      // other
      std::vector<double> m_timeReference; 
      int m_winSize = 21; 
      int m_i0 = 0; 
      int m_i1 = 0; 
      std::vector<int> m_lastEntries; 
    };
 
  } // end namespace TOP
} // end namespace Belle2