TOPGainEfficiencyCalculatorModule.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2017 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Maeda Yosuke, Okuto Rikuya                               *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#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 ~TOPGainEfficiencyCalculatorModule();
 
    virtual void initialize() override;
 
    virtual void beginRun() override;
 
    virtual void event() override;
 
    virtual void endRun() 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