DQMHistAnalysisSVDOnMiraBelle.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 <dqm/core/DQMHistAnalysis.h>
#include <vxd/dataobjects/VxdID.h>
#include <vxd/geometry/GeoTools.h>
#include <svd/dbobjects/SVDDQMPlotsConfiguration.h>
#include <framework/database/DBObjPtr.h>
 
#include <vector>
 
#include <TROOT.h>
#include <TCanvas.h>
#include <TH1F.h>
#include <TH2F.h>
 
namespace Belle2 {

  class DQMHistAnalysisSVDOnMiraBelleModule final : public DQMHistAnalysisModule {
 
  public:

    DQMHistAnalysisSVDOnMiraBelleModule();
    void initialize() override final;
    void beginRun() override final;
    void event() override final;
    void endRun() override final;
    void terminate() override final;
 
  private:

    TCanvas* m_c_avgEfficiency = nullptr; 
    TCanvas* m_c_avgOffOccupancy = nullptr; 
    TCanvas* m_c_MPVChargeClusterOnTrack = nullptr; 
    TCanvas* m_c_MPVSNRClusterOnTrack = nullptr; 
    TCanvas* m_c_MPVTimeClusterOnTrack = nullptr; 
    TCanvas* m_c_avgMaxBinClusterOnTrack =
      nullptr; 
    TCanvas* m_c_MeanSVDEventT0 = nullptr; 
 
    DBObjPtr<SVDDQMPlotsConfiguration> m_svdPlotsConfig; 

    MonitoringObject* m_monObj = nullptr;

    std::vector<VxdID> m_SVDModules;

    const VXD::GeoTools* m_gTools = nullptr;

    std::vector<std::string> m_listOfSensorsToMonitor;

    std::pair<float, float> avgOccupancyUV(TH1F* hU, TH1F* hV, int nEvents, int layer = -1, int ladder = -1, int sensor = -1) const;

    std::pair<float, float>  avgOccupancyGrpId0UV(int iLayer, int nEvents) const;

    std::pair<float, float> avgEfficiencyUV(TH2F* hMCU, TH2F* hMCV, TH2F* hFTU, TH2F* hFTV, int layer = -1, int ladder = -1,
                                            int sensor = -1) const;

    float xForMaxY(TH1F* h) const;

    float histFWHM(TH1F* h) const;

    void addVariable(std::string name, std::pair<float, float>& varUV);

    Int_t findBinY(Int_t layer, Int_t sensor) const
    {
      // should take the method from SVDSummaryPlot (->put as static)
      if (layer == 3)
        return sensor; //2
      if (layer == 4)
        return 2 + 1 + sensor; //6
      if (layer == 5)
        return 6 + 1 + sensor; // 11
      if (layer == 6)
        return 11 + 1 + sensor; // 17
      else
        return -1;
    }

    Int_t getNumberOfLadders(Int_t layer) const
    {
      if (layer == 3)
        return 7;
      if (layer == 4)
        return 10;
      if (layer == 5)
        return 12;
      if (layer == 6)
        return 16;
      else
        return -1;
    }

    Int_t getNumberOfSensors(Int_t layer) const
    {
      if (layer == 3)
        return 2;
      if (layer == 4)
        return 3;
      if (layer == 5)
        return 4;
      if (layer == 6)
        return 5;
      else
        return -1;
    }
 
  };

#define SetVariable(x) addVariable(#x , x)

} // end namespace Belle2