HLTTauDQMPathSummaryPlotter.cc

Go to the documentation of this file.

#include "DQMOffline/Trigger/interface/HLTTauDQMPathSummaryPlotter.h"
#include "DQMOffline/Trigger/interface/HLTTauDQMPath.h"

HLTTauDQMPathSummaryPlotter::HLTTauDQMPathSummaryPlotter(const edm::ParameterSet& pset,
                                                         bool doRefAnalysis,
                                                         const std::string& dqmBaseFolder,
                                                         double hltMatchDr)
    : HLTTauDQMPlotter(pset, dqmBaseFolder), hltMatchDr_(hltMatchDr), doRefAnalysis_(doRefAnalysis) {}

HLTTauDQMPathSummaryPlotter::~HLTTauDQMPathSummaryPlotter() = default;

void HLTTauDQMPathSummaryPlotter::bookHistograms(HistoWrapper& iWrapper, DQMStore::IBooker& iBooker) {
  if (!isValid() || pathObjects_.empty())
    return;

  //Create the histograms
  iBooker.setCurrentFolder(triggerTag() + "/helpers");

  all_events = iWrapper.book1D(iBooker, "RefEvents", "All events", pathObjects_.size(), 0, pathObjects_.size(), kVital);
  accepted_events = iWrapper.book1D(iBooker,
                                    "PathTriggerBits",
                                    "Accepted Events per Path;;entries",
                                    pathObjects_.size(),
                                    0,
                                    pathObjects_.size(),
                                    kVital);
  for (size_t i = 0; i < pathObjects_.size(); ++i) {
    if (all_events)
      all_events->setBinLabel(i + 1, pathObjects_[i]->getPathName());
    if (accepted_events)
      accepted_events->setBinLabel(i + 1, pathObjects_[i]->getPathName());
  }

  iBooker.setCurrentFolder(triggerTag());
}

void HLTTauDQMPathSummaryPlotter::analyze(const edm::TriggerResults& triggerResults,
                                          const trigger::TriggerEvent& triggerEvent,
                                          const HLTTauDQMOfflineObjects& refCollection) {
  if (doRefAnalysis_) {
    std::vector<HLTTauDQMPath::Object> triggerObjs;
    std::vector<HLTTauDQMPath::Object> matchedTriggerObjs;
    HLTTauDQMOfflineObjects matchedOfflineObjs;

    for (size_t i = 0; i < pathObjects_.size(); ++i) {
      const HLTTauDQMPath* path = pathObjects_[i];
      const int lastFilter = path->filtersSize() - 1;

      if (path->goodOfflineEvent(lastFilter, refCollection)) {
        if (all_events)
          all_events->Fill(i + 0.5);
      }
      if (path->fired(triggerResults)) {
        triggerObjs.clear();
        matchedTriggerObjs.clear();
        matchedOfflineObjs.clear();
        path->getFilterObjects(triggerEvent, lastFilter, triggerObjs);
        if (path->offlineMatching(
                lastFilter, triggerObjs, refCollection, hltMatchDr_, matchedTriggerObjs, matchedOfflineObjs)) {
          if (accepted_events)
            accepted_events->Fill(i + 0.5);
        }
      }
    }
  } else {
    for (size_t i = 0; i < pathObjects_.size(); ++i) {
      const HLTTauDQMPath* path = pathObjects_[i];
      if (all_events)
        all_events->Fill(i + 0.5);
      if (path->fired(triggerResults)) {
        if (accepted_events)
          accepted_events->Fill(i + 0.5);
      }
    }
  }
}