CDCDedxBadWireAlgorithm.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 <algorithm>
#include <string>
#include <iostream>
#include <fstream>
 
#include <TCanvas.h>
#include <TH2F.h>
#include <TH1I.h>
#include <TMath.h>
#include <TLine.h>
#include <TStyle.h>
#include <TPaveText.h>
#include <TText.h>
#include <TLegend.h>
 
#include <framework/database/DBObjPtr.h>
#include <framework/database/IntervalOfValidity.h>
 
#include <reconstruction/dbobjects/CDCDedxBadWires.h>
#include <reconstruction/dbobjects/CDCDedxWireGain.h>
 
#include <calibration/CalibrationAlgorithm.h>
 
#include <cdc/dbobjects/CDCGeometry.h>
#include <cdc/geometry/CDCGeometryPar.h>
#include <cdc/geometry/CDCGeometryParConstants.h>
 
 
namespace Belle2 {
  class CDCDedxBadWireAlgorithm : public CalibrationAlgorithm {
 
  public:
 
    CDCDedxBadWireAlgorithm();
 
    virtual ~CDCDedxBadWireAlgorithm() {}
 
    void setMonitoringPlots(bool value = false) {m_isMakePlots = value;}
 
    void setHighFracThres(double value) {m_fracThres = value;}
 
    void setRMSThres(double value) {m_rmsThres = value;}
 
    void setMeanThres(double value) {m_meanThres = value;}
 
    void setADC(bool value = false)
    {
      m_isADC = value;
      if (m_isADC) m_varMax = 1000.0;
      else m_varMax = 7.0;
    }
 
    void setHistPars(int nbin, double min, double max)
    {
      m_varBins = nbin;
      m_varMin = min;
      m_varMax = max;
    }
 
    void getExpRunInfo();
 
    void plotWireDist(const std::vector<double>& inwires, std::map<int, std::vector<double>>& vhitvar);
 
    void printCanvas(TList* list, TList* hflist, Color_t color);
 
    void plotBadWireMap(const std::vector<double>& vbadwires, const std::vector<double>& vdeadwires);
 
    TH2F* getHistoPattern(const std::vector<double>& inwires, const std::string& suffix, int& total);
 
    void plotQaPars(std::map<int, std::vector<double>>& qapars);
 
    void plotEventStats();
 
    void setTextCosmetics(TPaveText* pt, double size)
    {
      pt->SetTextAlign(11);
      pt->SetFillStyle(3001);
      pt->SetLineColor(2);
      pt->SetTextFont(82);
      pt->SetTextSize(size);
    }
 
    void setHistCosmetics(TH2F* hist, Color_t color)
    {
      hist->SetMarkerStyle(20);
      hist->SetMarkerSize(0.3);
      hist->SetMarkerColor(color);
      hist->SetFillColor(color);
      hist->SetStats(0);
    }
 
  protected:
 
    virtual EResult calibrate() override;
 
  private:
 
    unsigned int c_nwireCDC; 
    bool m_isMakePlots; 
    bool m_isADC; 
    int m_varBins; 
    double m_varMin; 
    double m_varMax; 
    double m_meanThres; 
    double m_rmsThres; 
    double m_fracThres; 
    double m_amean; 
    double m_arms; 
    std::string m_varName; 
    std::string m_suffix; 
    DBObjPtr<CDCDedxBadWires> m_DBBadWires; 
    DBObjPtr<CDCDedxWireGain> m_DBWireGains; 
    DBObjPtr<CDCGeometry> m_cdcGeo; 
  };
} // namespace Belle2