DQMHistAnalysisTOP.cc

/**************************************************************************
 * 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.                  *
 **************************************************************************/
 
#include <dqm/analysis/modules/DQMHistAnalysisTOP.h>
#include <boost/format.hpp>
#include <daq/slc/base/StringUtil.h>
#include <TClass.h>
#include <TH2.h>
#include "TROOT.h"
 
using namespace std;
using namespace Belle2;
using boost::format;
 
//-----------------------------------------------------------------
//                 Register the Module
//-----------------------------------------------------------------
REG_MODULE(DQMHistAnalysisTOP)
 
//-----------------------------------------------------------------
//                 Implementation
//-----------------------------------------------------------------
 
DQMHistAnalysisTOPModule::DQMHistAnalysisTOPModule()
  : DQMHistAnalysisModule()
{
  //Parameter definition
  B2DEBUG(20, "DQMHistAnalysisTOP: Constructor done.");
}
 
 
DQMHistAnalysisTOPModule::~DQMHistAnalysisTOPModule() { }
 
void DQMHistAnalysisTOPModule::initialize()
{
  gROOT->cd();
  B2DEBUG(20, "DQMHistAnalysisTOP: initialized.");
 
  m_c_goodHitsMean = new TCanvas("TOP/c_good_hits_mean");
  m_c_goodHitsRMS = new TCanvas("TOP/c_good_hits_rms");
  m_c_badHitsMean = new TCanvas("TOP/c_bad_hits_mean");
  m_c_badHitsRMS = new TCanvas("TOP/c_bad_hits_rms");
 
  //using c2_Name to avoid an overlap on default c_Name
  for (int i = 1; i <= 16; i++) {
    m_c_good_hits_xy_[i] = new TCanvas(Form("TOP/c2_good_hits_xy_%d", i));
    m_c_bad_hits_xy_[i] = new TCanvas(Form("TOP/c2_bad_hits_xy_%d", i));
    m_c_good_hits_asics_[i] = new TCanvas(Form("TOP/c2_good_hits_asics_%d", i));
    m_c_bad_hits_asics_[i] = new TCanvas(Form("TOP/c2_bad_hits_asics_%d", i));
  }
 
  m_h_goodHitsMean = new TH1F("TOP/good_hits_mean", "Mean of good hits per event", 16, 0.5, 16.5);
  m_h_goodHitsRMS = new TH1F("TOP/good_hits_rms", "RMS of good hits per event", 16, 0.5, 16.5);
  m_h_badHitsMean = new TH1F("TOP/bad_hits_mean", "Mean of bad hits per event", 16, 0.5, 16.5);
  m_h_badHitsRMS = new TH1F("TOP/bad_hits_rms", "RMS of bad hits per event", 16, 0.5, 16.5);
 
  m_h_goodHitsMean->GetXaxis()->SetTitle("slot number");
  m_h_goodHitsMean->GetYaxis()->SetTitle("hits per event");
  m_h_goodHitsRMS->GetXaxis()->SetTitle("slot number");
  m_h_goodHitsRMS->GetYaxis()->SetTitle("hits per event");
 
  m_h_badHitsMean->GetXaxis()->SetTitle("slot number");
  m_h_badHitsMean->GetYaxis()->SetTitle("hits per event");
  m_h_badHitsRMS->GetXaxis()->SetTitle("slot number");
  m_h_badHitsRMS->GetYaxis()->SetTitle("hits per event");
 
  m_h_goodHitsMean->SetStats(kFALSE);
  m_h_goodHitsRMS->SetStats(kFALSE);
  m_h_badHitsMean->SetStats(kFALSE);
  m_h_badHitsRMS->SetStats(kFALSE);
 
  m_h_goodHitsMean->SetMinimum(0);
  m_h_goodHitsRMS->SetMinimum(0);
  m_h_badHitsMean->SetMinimum(0);
  m_h_badHitsRMS->SetMinimum(0);
 
  m_line1 = new TLine(0.5, 215, 16.5, 215);
  m_line2 = new TLine(0.5, 235, 16.5, 235);
  m_line1->SetLineWidth(2);
  m_line2->SetLineWidth(2);
  m_line1->SetLineColor(kRed);
  m_line2->SetLineColor(kRed);
 
  m_text1 = new TPaveText(1, 435, 10, 500, "NB");
  m_text1->SetFillColorAlpha(kWhite, 0);
  m_text1->SetBorderSize(0);
 
  m_text2 = new TPaveText(0.2, 0.70, 0.50, 0.9, "NB");
  m_text2->SetFillColorAlpha(kWhite, 0);
  m_text2->SetBorderSize(0);
}
 
 
void DQMHistAnalysisTOPModule::beginRun()
{
  //B2DEBUG(20, "DQMHistAnalysisTOP: beginRun called.");
}
 
TH1* DQMHistAnalysisTOPModule::find_histo_in_canvas(TString histo_name)
{
  StringList s = StringUtil::split(histo_name.Data(), '/');
  std::string dirname = s[0];
  std::string hname = s[1];
  std::string canvas_name = dirname + "/c_" + hname;
 
  auto cobj = findCanvas(canvas_name);
  if (cobj == nullptr) return nullptr;
 
  TIter nextkey(((TCanvas*)cobj)->GetListOfPrimitives());
  TObject* obj{};
  while ((obj = dynamic_cast<TObject*>(nextkey()))) {
    if (obj->IsA()->InheritsFrom("TH1")) {
      if (obj->GetName() == histo_name)
        return  dynamic_cast<TH1*>(obj);
    }
  }
  return nullptr;
}
 
void DQMHistAnalysisTOPModule::event()
{
  for (int i = 1; i <= 16; i++) {
    string hname1 = str(format("TOP/good_hits_per_event%1%") % (i));;
    string hname2 = str(format("TOP/bad_hits_per_event%1%") % (i));;
    TH1* h1 = findHist(hname1);
    TH1* h2 = findHist(hname2);
    //TH1* h1 = find_histo_in_canvas(hname1);
    //TH1* h2 = find_histo_in_canvas(hname2);
    if (h1 != NULL) {
      m_h_goodHitsMean->SetBinContent(i, h1->GetMean());
      m_h_goodHitsRMS->SetBinContent(i, h1->GetRMS());
    } else {
      m_h_goodHitsMean->SetBinContent(i, 0);
      m_h_goodHitsRMS->SetBinContent(i, 0);
    }
    if (h2 != NULL) {
      m_h_badHitsMean->SetBinContent(i, h2->GetMean());
      m_h_badHitsRMS->SetBinContent(i, h2->GetRMS());
    } else {
      m_h_badHitsMean->SetBinContent(i, 0);
      m_h_badHitsRMS->SetBinContent(i, 0);
    }
  }
 
  TH2F* hraw = (TH2F*)findHist("TOP/window_vs_slot");
  double exRatio(0.0);
  double totalraw(0.0);
  double totalbadraw(0.0);
  if (hraw != NULL) {
    totalraw = hraw->GetEntries();
    for (int i = 1; i <= 16; i++) {
      for (int j = 1; j <= 512; j++) {
        double nhraw = hraw->GetBinContent(i, j);
        if (j < 215 || j > 235) totalbadraw += nhraw ;
      }
    }
  }
  if (totalraw > 0) exRatio = totalbadraw * 1.0 / totalraw;
 
  m_text1->Clear();
  m_text1->AddText(Form("Ratio of entries outside of red lines: %.2f %%", exRatio * 100.0));
  if (exRatio > 0.01) {
    m_text1->AddText(">1% bad, report to TOP experts!");
  } else {
    m_text1->AddText("<0.1% good, 0.1-1% recoverable.");
  }
 
  //addHist("", m_h_goodHitsMean->GetName(), m_h_goodHitsMean);
 
  TCanvas* c1 = findCanvas("TOP/c_hitsPerEvent");
  //TH1* h1=find_histo_in_canvas("TOP/hitsPerEvent");
  if (c1 != NULL) {
    c1->SetName("TOP/c_hitsPerEvent_top");
  }
  //if (h1!=NULL) {
  //  h1->SetName("TOP/hitsPerEvent_top");
  //}
 
  TCanvas* c2 = findCanvas("TOP/c_window_vs_slot");
  if (c2 != NULL) {
    c2->cd();
    if (exRatio > 0.01) c2->Pad()->SetFillColor(kRed);
    else c2->Pad()->SetFillColor(kWhite);
    m_line1->Draw();
    m_line2->Draw();
    m_text1->Draw();
  }
 
  TH1F* hBoolEvtMonitor = (TH1F*)findHist("TOP/BoolEvtMonitor");
  double badRatio(0.0);
  int totalBadEvts(0);
  int totalEvts(0);
  if (hBoolEvtMonitor != NULL) {
    totalEvts = hBoolEvtMonitor->GetEntries();
    totalBadEvts = hBoolEvtMonitor->GetBinContent(2);
  }
  if (totalEvts > 0) badRatio = totalBadEvts * 1.0 / totalEvts;
 
  m_text2->Clear();
  m_text2->AddText(Form("fraction of deviating hits: %.4f %%", badRatio * 100.0));
 
  TCanvas* c3 = findCanvas("TOP/c_BoolEvtMonitor");
  if (c3 != NULL) {
    c3->cd();
    if (badRatio > 0.0001) c3->Pad()->SetFillColor(kRed);
    else c3->Pad()->SetFillColor(kWhite);
    m_text2->Draw();
  }
 
  //obtaining the total yield for 16 2D-plots
  double Ntotal_good_hits_xy(0.0);
  double Ntotal_bad_hits_xy(0.0);
  double Ntotal_good_hits_asics(0.0);
  double Ntotal_bad_hits_asics(0.0);
  for (int module = 1; module <= 16; module++) {
    string hnameTmp1 = str(format("TOP/good_hits_xy_%1%") % (module));;
    TH1* h2DTmp1 = findHist(hnameTmp1);
    if (h2DTmp1 != NULL)   Ntotal_good_hits_xy += h2DTmp1->Integral();
 
    string hnameTmp2 = str(format("TOP/bad_hits_xy_%1%") % (module));;
    TH1* h2DTmp2 = findHist(hnameTmp2);
    if (h2DTmp2 != NULL) Ntotal_bad_hits_xy += h2DTmp2->Integral();
 
    string hnameTmp3 = str(format("TOP/good_hits_asics_%1%") % (module));;
    TH1* h2DTmp3 = findHist(hnameTmp3);
    if (h2DTmp3 != NULL) Ntotal_good_hits_asics += h2DTmp3->Integral();
 
    string hnameTmp4 = str(format("TOP/bad_hits_asics_%1%") % (module));;
    TH1* h2DTmp4 = findHist(hnameTmp4);
    if (h2DTmp4 != NULL) Ntotal_bad_hits_asics += h2DTmp4->Integral();
  }
 
  // reset the maximum z-axis of 16 2D plots: 3 times of average for good hits; 30 times of average for bad hits
  for (int i = 1; i <= 16; i++) {
    m_c_good_hits_xy_[i]->Clear();
    m_c_good_hits_xy_[i]->cd();
    TH2F* h2Dscale_xy = (TH2F*)findHist(Form("TOP/good_hits_xy_%d", i));
    if (h2Dscale_xy != NULL && Ntotal_good_hits_xy > 0) {
      h2Dscale_xy->GetZaxis()->SetRangeUser(0, Ntotal_good_hits_xy / 2500.0);
      h2Dscale_xy->Draw();
    }
    m_c_good_hits_xy_[i]->Modified();
  }
 
 
  for (int i = 1; i <= 16; i++) {
    m_c_bad_hits_xy_[i]->Clear();
    m_c_bad_hits_xy_[i]->cd();
    TH2F* h2Dscale_xy = (TH2F*)findHist(Form("TOP/bad_hits_xy_%d", i));
    if (h2Dscale_xy != NULL && Ntotal_bad_hits_xy > 0) {
      h2Dscale_xy->GetZaxis()->SetRangeUser(0, Ntotal_bad_hits_xy / 250.0);
      h2Dscale_xy->Draw();
    }
    m_c_bad_hits_xy_[i]->Modified();
  }
 
  for (int i = 1; i <= 16; i++) {
    m_c_good_hits_asics_[i]->Clear();
    m_c_good_hits_asics_[i]->cd();
    TH2F* h2Dscale_asics = (TH2F*)findHist(Form("TOP/good_hits_asics_%d", i));
    if (h2Dscale_asics != NULL && Ntotal_good_hits_asics > 0) {
      h2Dscale_asics->GetZaxis()->SetRangeUser(0, Ntotal_good_hits_asics / 2500.0);
      h2Dscale_asics->Draw();
    }
    m_c_good_hits_asics_[i]->Modified();
  }
 
  for (int i = 1; i <= 16; i++) {
    m_c_bad_hits_asics_[i]->Clear();
    m_c_bad_hits_asics_[i]->cd();
    TH2F* h2Dscale_asics = (TH2F*)findHist(Form("TOP/bad_hits_asics_%d", i));
    if (h2Dscale_asics != NULL && Ntotal_bad_hits_asics > 0) {
      h2Dscale_asics->GetZaxis()->SetRangeUser(0, Ntotal_bad_hits_asics / 250.0);
      h2Dscale_asics->Draw();
    }
    m_c_bad_hits_asics_[i]->Modified();
  }
 
  m_c_goodHitsMean->Clear();
  m_c_goodHitsMean->cd();
  m_h_goodHitsMean->Draw();
  m_c_goodHitsMean->Modified();
 
  m_c_goodHitsRMS->Clear();
  m_c_goodHitsRMS->cd();
  m_h_goodHitsRMS->Draw();
  m_c_goodHitsRMS->Modified();
 
  m_c_badHitsMean->Clear();
  m_c_badHitsMean->cd();
  m_h_badHitsMean->Draw();
  m_c_badHitsMean->Modified();
 
  m_c_badHitsRMS->Clear();
  m_c_badHitsRMS->cd();
  m_h_badHitsRMS->Draw();
  m_c_badHitsRMS->Modified();
}
 
void DQMHistAnalysisTOPModule::endRun()
{
  B2DEBUG(20, "DQMHistAnalysisTOP : endRun called");
}
 
 
void DQMHistAnalysisTOPModule::terminate()
{
  B2DEBUG(20, "terminate called");
}