TOPPulseHeightAlgorithm.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
#include <calibration/CalibrationAlgorithm.h>
#include <TH2F.h>
#include <TH1D.h>
#include <TFile.h>
#include <TTree.h>
#include <top/dbobjects/TOPCalChannelPulseHeight.h>
#include <top/dbobjects/TOPCalChannelThresholdEff.h>
 
namespace Belle2 {
  namespace TOP {

    class TOPPulseHeightAlgorithm : public CalibrationAlgorithm {
    public:

      TOPPulseHeightAlgorithm();

      virtual ~TOPPulseHeightAlgorithm() {}

      void setMinEntries(int minEntries) {m_minEntries = minEntries;}

      void setXmin(double xmin) {m_xmin = xmin;}

      void setMinFracCalibrated(double minFract) {m_minCalibrated = minFract;}
 
    private:

      virtual EResult calibrate() final;

      int fitChannels(std::shared_ptr<TH2F> h);

      int fitPulseHeight(TH1D* h);

      double getEfficiency(TH1D* h, TF1* func);
 
      // algorithm parameters
      int m_minEntries = 100; 
      double m_xmin = 100.0; 
      double m_minCalibrated = 0.90; 
 
      // ntuple variables
      int m_slot = 0; 
      unsigned m_channel = 0; 
      int m_numPhot = 0; 
      float m_x0 = 0; 
      float m_p1 = 0; 
      float m_p2 = 0; 
      float m_x0err = 0; 
      float m_p1err = 0; 
      float m_p2err = 0; 
      float m_effi = 0; 
      float m_meanHist = 0; 
      float m_meanFunc = 0; 
      float m_chi2 = 0;  
      int m_ndf = 0; 
      int m_fitStatus = -1; 
      bool m_good = false; 
 
      // output file and tree
      TFile* m_file = 0; 
      TTree* m_tree = 0; 
 
      // output DB objects
      TOPCalChannelPulseHeight* m_pulseHeight = 0; 
      TOPCalChannelThresholdEff* m_thresholdEffi = 0; 
 
    };
 
  } // end namespace TOP
} // end namespace Belle2