WireEfficiencyAlgorithm.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 <TEfficiency.h>
#include <calibration/CalibrationAlgorithm.h>
#include <framework/database/DBObjPtr.h>
 
#include <cdc/dbobjects/CDCGeometry.h>
#include <cdc/dbobjects/CDCBadWires.h>
 
#include <string>
 
namespace Belle2 {
  namespace CDC {
    class WireEfficiencyAlgorithm: public CalibrationAlgorithm {
    public:
      WireEfficiencyAlgorithm();
      ~WireEfficiencyAlgorithm() override { delete m_efficiencyList; delete m_badWireList; }
      WireEfficiencyAlgorithm(const WireEfficiencyAlgorithm&) = delete;

      WireEfficiencyAlgorithm& operator=(const WireEfficiencyAlgorithm&) = delete;
      void setHistFileName(const std::string& name) {m_outputFileName = "histWireEff_" + name + ".root";}
      void setAverageOccupancyThreshold(const float threshold) { m_averageOccupancyThreshold = threshold; };
    protected:
      EResult calibrate() override;
      bool hasEnoughData();
      void buildEfficiencies();
      void detectBadWires();
      static double chiTest(TGraphAsymmErrors* graph1, TGraphAsymmErrors* graph2, double minVale, double maxValue);
    private:
      float m_averageOccupancyThreshold = 4000.0; 
      TList* m_efficiencyList = new TList(); 
      std::string m_outputFileName = "wire_efficiencies.root"; 
      DBObjPtr<CDCGeometry> m_cdcGeo; 
      CDCBadWires* m_badWireList = nullptr; 
    };
 
 
  }
}