TOPModuleT0DeltaTAlgorithm.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 TOPModuleT0DeltaTAlgorithm : public CalibrationAlgorithm {
    public:
 
      TOPModuleT0DeltaTAlgorithm();
 
      virtual ~TOPModuleT0DeltaTAlgorithm() {}
 
      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;
 
      int fitSingleGaus(std::shared_ptr<TH1F> h);
 
      int fitDoubleGaus(std::shared_ptr<TH1F> h);
 
      int fitHistogram(std::shared_ptr<TH1F> h);
 
      // algorithm parameters
      int m_minEntries = 10; 
      double m_sigmaCoreInit = 0.120; 
      double m_sigmaTailInit = 0.240; 
      double m_tailFractInit = 0.20; 
      int m_cutoffEntries = 100; 
      double m_minError = 0.100; 
      // temporary store for results of a single fit
      double m_delT0 = 0; 
      double m_error = 0; 
      double m_chi2 = 0;  
      double m_ndf = 0;  
    };
 
  } // end namespace TOP
} // end namespace Belle2