TOPCommonT0BFAlgorithm.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>
 
namespace Belle2 {
  namespace TOP {
 
    class TOPCommonT0BFAlgorithm : public CalibrationAlgorithm {
    public:
 
      TOPCommonT0BFAlgorithm();
 
      virtual ~TOPCommonT0BFAlgorithm() {}
 
      void setMinEntries(int minEntries) {m_minEntries = minEntries;}
 
      void setFitInitializers(double sigmaCore, double sigmaTail, double tailFract)
      {
        m_sigmaCoreInit = sigmaCore;
        m_sigmaTailInit = sigmaTail;
        m_tailFractInit = tailFract;
      }
 
      void setCutoffEntries(int cutoffEntries) {m_cutoffEntries = cutoffEntries;}
 
      void setMinError(double minError) {m_minError = minError;}
 
    private:
 
      virtual EResult calibrate() final;
 
      std::shared_ptr<TH1F> getHistogram();
 
      int fitSingleGaus(std::shared_ptr<TH1F> h);
 
      int fitDoubleGaus(std::shared_ptr<TH1F> h);
 
      // algorithm parameters
      int m_minEntries = 10; 
      double m_sigmaCoreInit = 0.060; 
      double m_sigmaTailInit = 0.120; 
      double m_tailFractInit = 0.20; 
      int m_cutoffEntries = 100; 
      double m_minError = 0.020; 
      // ntuple variables
      int m_expNo = 0; 
      int m_runNo = 0; 
      int m_runLast = 0; 
      float m_fittedOffset = 0; 
      float m_offset = 0; 
      float m_offsetError = 0; 
      float m_sigmaCore = 0; 
      float m_sigmaTail = 0; 
      float m_tailFract = 0; 
      float m_chi2 = 0;  
      float m_integral = 0; 
      int m_numEvents = 0; 
      int m_fitStatus = -1; 
      // other
      double m_bunchTimeSep = 0; 
    };
 
  } // end namespace TOP
} // end namespace Belle2