CDCDedxWireGainAlgorithm.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: jvbennett                                                *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <reconstruction/dbobjects/CDCDedxWireGain.h>
#include <reconstruction/dbobjects/CDCDedxBadWires.h>
#include <calibration/CalibrationAlgorithm.h>
#include <framework/database/DBObjPtr.h>
#include <string>
#include <vector>
#include "TH1D.h"
#include "TH2F.h"
 
namespace Belle2 {
  class CDCDedxWireGainAlgorithm : public CalibrationAlgorithm {
 
  public:
 
    CDCDedxWireGainAlgorithm();
 
    virtual ~CDCDedxWireGainAlgorithm() {}
 
    void setMergePayload(bool value = true) {isMergePayload = value;}
 
    void setMonitoringPlots(bool value = false) {isMakePlots = value;}
 
    void setLocalTrucation(bool value = false) {isLTruc = value;}
 
    void setLayerScaling(bool value = false) {isLayerScale = value;}
 
    void generateNewPayloads(std::vector<double> dedxTruncmean);
 
    void getTrucationBins(TH1D* htemp, int& binlow, int& binhigh);
 
    void plotBadWires(int nDeadwires, int oDeadwires, int Badwires);
 
    TH2F* getHistoPattern(TString badFileName, TString suffix);
 
    double getIndexVal(int iWire, TString what);
 
    double getLayerAverage(std::vector<double> tempWire);
 
  protected:
 
    virtual EResult calibrate() override;
 
 
  private:
 
    std::string m_badWireFPath; 
    std::string m_badWireFName; 
    bool isMakePlots; 
    bool isMergePayload; 
    bool isLTruc; 
    bool isLayerScale; 
    int fdEdxBins; 
    double fdEdxMin; 
    double fdEdxMax; 
    double fTrucMin; 
    double fTrucMax; 
    int fStartRun; 
    std::vector<double> flayerAvg; 
    DBObjPtr<CDCDedxWireGain> m_DBWireGains; 
    DBObjPtr<CDCDedxBadWires> m_DBBadWires; 
  };
} // namespace Belle2