TOPPulseHeightAlgorithm.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 <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