LaserCalibratorFit.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2017 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Wenlong Yuan                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <TH1F.h>
#include <TF1.h>
#include <string>
 
namespace Belle2 {
  namespace TOP {
 
    double fcnCB(double* x, double* par);
 
    double fcnCB2(double* x, double* par);
 
    class LaserCalibratorFit {
    public:
      explicit LaserCalibratorFit(unsigned moduleID);
 
      ~LaserCalibratorFit();
 
      void setHist(const std::vector<TH1F*>& hist);
 
      void setFitMethod(const std::string& method) {m_fitMethod = method;}
 
      void setFitRange(double xmin = -200, double xmax = 200) {m_xmin = xmin; m_xmax = xmax;}
 
      int  fitChannel(unsigned channel);
 
      void writeFile(const std::string& outfile);
 
      int getModuleID() const {return m_moduleID;}
 
      double getFitChisq(unsigned channel);
 
      std::string getFitMethod() const {return m_fitMethod;}
 
      std::vector<double> getMaxPos()  {return m_maxpos;}
 
      double getFitT()  {return m_fitT;}
 
      double getFitTErr() {return m_fitTErr;}
 
 
    private:
 
      TF1* makeGFit(unsigned channel);
 
      TF1* makeCBFit(unsigned channel);
 
      TF1* makeCB2Fit(unsigned channel, bool minOut);
 
      unsigned m_moduleID = 0; 
      std::vector<double> m_maxpos; 
      std::vector<double> m_maxpos_error; 
      double m_xmin = 0; 
      double m_xmax = 0; 
      double m_fitT = 0; 
      double m_fitTErr = 0.; 
      std::vector<TH1F*> m_hist; 
      std::vector<TF1*> m_func; 
      std::string m_fitMethod; 
    };
  }// TOP namespace
}//Belle2 namespace