TOPGainEfficiencyCalculatorModule.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 <framework/core/HistoModule.h>
#include <top/modules/TOPGainEfficiencyMonitor/TOPLaserHitSelectorModule.h>
 
#include <string>
 
#include <TTree.h>
#include <TH1.h>
#include <TH2.h>
#include <TF1.h>
 
namespace Belle2 {

  class TOPGainEfficiencyCalculatorModule : public HistoModule {
 
  public:

    enum { c_NParameterGainFit = 6 };

    enum { c_NPlotsPerChannel = 3, c_NChannelPerPage = 4 };

    enum EHistogramType { c_LoadForFitHeight = 1, c_LoadHitRateHeight = 2,
                          c_LoadForFitIntegral = 3, c_LoadHitRateIntegral = 4
                        };

    TOPGainEfficiencyCalculatorModule();

    virtual void initialize() override;

    virtual void terminate() override;

    virtual void defineHisto() override;

    void LoadHistograms(const std::string& histotype);//histotype {Height_gain,Height_efficiency,Integral_gain:}

    void FitHistograms(EHistogramType LoadHisto);

    void DummyFillBranch(EHistogramType LoadHisto);

    void DrawResult(const std::string& histotype, EHistogramType LoadHisto);

    static double TOPGainFunc(double* var, double* par);

    static double FindPeakForSmallerXThan(TH1* histo, double xmax = 0);
 
  private:
 
    TTree* m_tree = 0; 
    std::vector<TBranch*> m_branch[4]; 
 
    TH2F* m_timeChargeHistogram[c_NChannelPerPMT] = {}; 
    TH1D* m_timeHistogram[c_NChannelPerPMT] = {}; 
    TH1D* m_chargeHistogram[c_NChannelPerPMT] = {}; 
    TH1F* m_nCalPulseHistogram =
      0; 
    TF1* m_funcForLaser[c_NChannelPerPMT] = {}; 
    TF1* m_funcForFitRange[c_NChannelPerPMT] = {}; 
    TF1* m_funcForFullRange[c_NChannelPerPMT] = {}; 
 
    std::string m_inputFile = ""; 
    std::string m_outputPDFFile =
      ""; 
 
    std::string m_fitoption =
      "L"; 
 
    short m_targetSlotId = 0; 
    short m_targetPmtId = 0; 
    short m_targetPmtChId = -1; 
    short m_hvDiff = 0; 
 
    float m_fitHalfWidth = 1.; 
    float m_threshold = 100; 
    float m_thresholdForIntegral =
      550; 
    float m_p0HeightIntegral = -50.0; 
    float m_p1HeightIntegral = 6.0; 
    float m_fitMax = 0; 
    float m_fracFit = 0.99; 
    float m_initialP0 = (float)(-1.); 
    float m_initialP1 = (float)(-1.); 
    float m_initialP2 = (float)(-1.); 
    float m_initialX0 = (float)(-1.); 
    float m_pedestalSigma = 10.; 
 
    short m_pixelId = 0; 
    short m_pmtChId = 0; 
    float m_hitTiming = 0; 
    float m_hitTimingSigma = 0; 
    int m_nEntries = 0; 
    int m_nCalPulse = 0; 
    int m_nOverflowEvents = 0; 
    float m_meanPulseHeight = 0; 
    float m_meanPulseHeightError = 0; 
    float m_gain = 0; 
    float m_efficiency = 0; 
    float m_p0 = 0; 
    float m_p1 = 0; 
    float m_p2 = 0; 
    float m_x0 = 0; 
    float m_p0Error = 0; 
    float m_p1Error = 0; 
    float m_p2Error = 0; 
    float m_x0Error = 0; 
    float m_chisquare = 0; 
    int m_ndf = 0; 
    float m_funcFullRangeIntegral = 0; 
    float m_funcFitRangeIntegral = 0; 
    float m_histoFitRangeIntegral = 0; 
    float m_histoMeanAboveThre = 0; 
  };

}  //namespace Belle2