DQMHistAnalysisTrackingER.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/DQMHistAnalysisTrackingER.h>
 
#include <TROOT.h>
 
#include <string>
 
using namespace std;
using namespace Belle2;
 
//-----------------------------------------------------------------
//                 Register the Module
//-----------------------------------------------------------------
REG_MODULE(DQMHistAnalysisTrackingER);
 
//-----------------------------------------------------------------
//                 Implementation
//-----------------------------------------------------------------
 
DQMHistAnalysisTrackingERModule::DQMHistAnalysisTrackingERModule()
  : DQMHistAnalysisModule()
{
 
  setDescription("DQM Analysis Module of the Tracking ER Plots.");
 
  addParam("onTimeHalfWidth", m_onTimeHalfWidth, "a cluster is on time if within ± onTimeHalfWidth [ns]", float(m_onTimeHalfWidth));
 
}
 
void DQMHistAnalysisTrackingERModule::initialize()
{
 
  gROOT->cd();
 
  // add MonitoringObject
  m_monObj = getMonitoringObject("trackingER");
 
}
 
 
void DQMHistAnalysisTrackingERModule::event()
{
  // Repeat this for all tracks (no suffix) and tracks from IP (_FromIP suffix)
  for (const string suffix : {"_FromIP", "_NotFromIP"}) {
    //compute fraction of tracks with no PXD hits
    TH1* hNoPXDHits = findHist("TrackingERDQM" + suffix + "/NoOfHitsInTrack_PXD");
    if (hNoPXDHits != nullptr) {
      int nTracks = hNoPXDHits->GetEntries();
      int nTracksNoPXD = hNoPXDHits->GetBinContent(1);
      if (nTracks > 0)
        m_monObj->setVariable("tracksNoPXDHit" + suffix, (double)nTracksNoPXD / nTracks);
    }
 
    // add average number of tracks per event to Mirabelle
    TH1* hnTracks = findHist("TrackingERDQM" + suffix + "/NoOfTracks");
    if (hnTracks != nullptr) {
      double averageNTracks = hnTracks->GetMean();
      m_monObj->setVariable("nTracksPerEvent" + suffix, averageNTracks);
    }
    TH1* hnVXDTracks = findHist("TrackingERDQM" + suffix + "/NoOfTracksInVXDOnly");
    if (hnVXDTracks != nullptr) {
      double averageNVXDTracks = hnVXDTracks->GetMean();
      m_monObj->setVariable("nVXDTracksPerEvent" + suffix, averageNVXDTracks);
    }
    TH1* hnCDCTracks = findHist("TrackingERDQM" + suffix + "/NoOfTracksInCDCOnly");
    if (hnCDCTracks != nullptr) {
      double averageNCDCTracks = hnCDCTracks->GetMean();
      m_monObj->setVariable("nCDCTracksPerEvent" + suffix, averageNCDCTracks);
    }
    TH1* hnVXDCDCTracks = findHist("TrackingERDQM" + suffix + "/NoOfTracksInVXDCDC");
    if (hnVXDCDCTracks != nullptr) {
      double averageNVXDCDCTracks = hnVXDCDCTracks->GetMean();
      m_monObj->setVariable("nVXDCDCTracksPerEvent" + suffix, averageNVXDCDCTracks);
    }
  }
 
  //fraction of offtime SVD hits
  //considering L3V and L456V clusters (separately)
 
  TH1* hSVDL3VTime = findHist("SVDClsTrk/SVDTRK_ClusterTimeV3");
  if (hSVDL3VTime != nullptr) {
    int all = hSVDL3VTime->GetEntries();
    if (all > 0) {
      int bin_min = hSVDL3VTime->GetXaxis()->FindBin(-m_onTimeHalfWidth);
      int bin_max = hSVDL3VTime->GetXaxis()->FindBin(+m_onTimeHalfWidth);
      int offtime = all - hSVDL3VTime->Integral(bin_min, bin_max);
      double offtimeL3Hits = (double)offtime / all;
      m_monObj->setVariable("offtimeL3Hits", offtimeL3Hits);
    }
  }
 
  TH1* hSVDL456VTime = findHist("SVDClsTrk/SVDTRK_ClusterTimeV456");
  if (hSVDL456VTime != nullptr) {
    int all = hSVDL456VTime->GetEntries();
    if (all > 0) {
      int bin_min = hSVDL456VTime->GetXaxis()->FindBin(-m_onTimeHalfWidth);
      int bin_max = hSVDL456VTime->GetXaxis()->FindBin(+m_onTimeHalfWidth);
      int offtime = all - hSVDL456VTime->Integral(bin_min, bin_max);
      double offtimeL456Hits = (double)offtime / all;
      m_monObj->setVariable("offtimeL456Hits", offtimeL456Hits);
    }
  }
 
  // Compute tracking abort rate for Mirabelle
  // Using the _FromIP suffix: these histograms are event-based, so it does not make any difference
  TH1* hAbort = findHist("TrackingERDQM_FromIP/NumberTrackingErrorFlags");
  TH1* hAbortReasons = findHist("TrackingERDQM_FromIP/TrackingErrorFlagsReasons");
  if (hAbort != nullptr) {
    double nEvents = hAbort->GetEntries();
    if (nEvents > 0) {
      m_monObj->setVariable("abortRate", hAbort->GetMean());
      if (hAbortReasons != nullptr) {
        m_monObj->setVariable("abortRateUnspecifiedPR", hAbortReasons->GetBinContent(2) / nEvents);
        m_monObj->setVariable("abortRateVXDTF2", hAbortReasons->GetBinContent(3) / nEvents);
        m_monObj->setVariable("abortRateSVDCKF", hAbortReasons->GetBinContent(4) / nEvents);
        m_monObj->setVariable("abortRatePXDCKF", hAbortReasons->GetBinContent(5) / nEvents);
        m_monObj->setVariable("abortRateSpacePoint", hAbortReasons->GetBinContent(6) / nEvents);
 
        // Scale the bin contents to contain relative instead of absolute values.
        for (int i = 1; i <= 6; i++) {
          hAbortReasons->SetBinContent(i, hAbortReasons->GetBinContent(i) / nEvents);
        }
      }
    }
  }
}