CDCDedxCosEdgeAlgorithm.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 <TH1D.h>
#include <TPaveText.h>
 
#include<map>
#include<vector>
 
#include <reconstruction/dbobjects/CDCDedxCosineEdge.h>
#include <calibration/CalibrationAlgorithm.h>
#include <framework/database/DBObjPtr.h>
 
enum fitstatus {FitOK, FitFailed, LowStats};
 
namespace Belle2 {
  class CDCDedxCosEdgeAlgorithm : public CalibrationAlgorithm {
 
  public:
 
    CDCDedxCosEdgeAlgorithm();
 
    virtual ~CDCDedxCosEdgeAlgorithm() {}
 
    void setMergePayload(bool value = true) {m_isMerge = value;}
 
    void setMonitoringPlots(bool value = false) {m_isMakePlots = value;}
 
    void setFitWidth(double value = 2.5) {m_sigLim = value;}
 
    void setCosineBins(unsigned int value) {m_npBins = value;}
 
    void setNegCosRange(double min, double max) {m_negMin = min; m_negMax = max;}
 
    void setPosCosRange(double min, double max) {m_posMin = min; m_posMax = max;}
 
    void setDedxHistPars(int value, double min, double max)
    {
      m_dedxBins = value;
      m_dedxMin = min;
      m_dedxMax = max;
    }
 
    void getExpRunInfo();
 
    void fitGaussianWRange(TH1D*& temphist, fitstatus& status);
 
    void createPayload(std::vector<std::vector<double>>& vfinalconst);
 
    void plotHist(std::vector<TH1D*>& hdedx, std::map<int, std::vector<double>>& fPars, std::string type);
 
    void plotFitPars(std::map<int, std::vector<double>>& fPars_Neg, std::map<int, std::vector<double>>& fPars_Pos);
 
    void plotConstants(std::vector<std::vector<double>>& vfinalconst);
 
    void plotStats();
 
    void setTextCosmetics(TPaveText* pt, double size)
    {
      pt->SetTextAlign(11);
      pt->SetFillStyle(3001);
      pt->SetLineColor(2);
      pt->SetTextFont(82);
      pt->SetTextSize(size);
    }
 
    void setHistCosmetics(TH1D& hist, Color_t color, double min, double max, double size)
    {
      hist.SetMarkerStyle(20);
      hist.SetMarkerSize(0.60);
      hist.SetMarkerColor(color);
      hist.SetLineColor(color);
      hist.SetStats(0);
      hist.GetYaxis()->SetRangeUser(min, max);
      hist.GetXaxis()->SetLabelSize(size);
      hist.GetYaxis()->SetLabelSize(size);
    }
 
  protected:
 
    virtual EResult calibrate() override;
 
  private:
 
    bool m_isMakePlots; 
    bool m_isMerge; 
    double m_sigLim; 
    unsigned int m_npBins; 
    double m_negMin; 
    double m_negMax; 
    double m_posMin; 
    double m_posMax; 
    int m_dedxBins; 
    double m_dedxMin; 
    double m_dedxMax; 
    std::string m_suffix; 
    DBObjPtr<CDCDedxCosineEdge> m_DBCosineCor; 
  };
} // namespace Belle2