L1TEfficiencyHarvesting.cc

Go to the documentation of this file.

// L1TMonitor includes
#include "DQMOffline/L1Trigger/interface/L1TEfficiencyHarvesting.h"
#include "FWCore/Framework/interface/MakerMacros.h"

using namespace edm;
using namespace std;

namespace dqmoffline {
  namespace l1t {
    L1TEfficiencyPlotHandler::L1TEfficiencyPlotHandler(const ParameterSet &ps, std::string plotName)
        : numeratorDir_(ps.getUntrackedParameter<std::string>("numeratorDir")),
          denominatorDir_(ps.getUntrackedParameter<std::string>("denominatorDir", numeratorDir_)),
          outputDir_(ps.getUntrackedParameter<std::string>("outputDir", numeratorDir_)),
          plotName_(plotName),
          numeratorSuffix_(ps.getUntrackedParameter<std::string>("numeratorSuffix", "Num")),
          denominatorSuffix_(ps.getUntrackedParameter<std::string>("denominatorSuffix", "Den")),
          h_efficiency_() {}

    L1TEfficiencyPlotHandler::L1TEfficiencyPlotHandler(const L1TEfficiencyPlotHandler &handler)
        : numeratorDir_(handler.numeratorDir_),
          denominatorDir_(handler.denominatorDir_),
          outputDir_(handler.outputDir_),
          plotName_(handler.plotName_),
          numeratorSuffix_(handler.numeratorSuffix_),
          denominatorSuffix_(handler.denominatorSuffix_),
          h_efficiency_(handler.h_efficiency_) {}

    void L1TEfficiencyPlotHandler::book(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) {
      edm::LogInfo("L1TEfficiencyPlotHandler")
          << "Booking efficiency histogram for " << outputDir_ << " and " << plotName_ << endl;

      std::string numeratorName = numeratorDir_ + "/" + plotName_ + numeratorSuffix_;
      std::string denominatorName = denominatorDir_ + "/" + plotName_ + denominatorSuffix_;
      MonitorElement *num = igetter.get(numeratorName);
      MonitorElement *den = igetter.get(denominatorName);

      if (!num || !den) {
        edm::LogWarning("L1TEfficiencyPlotHandler") << (!num && !den ? numeratorName + " && " + denominatorName
                                                        : !num       ? numeratorName
                                                                     : denominatorName)
                                                    << " not gettable. Quitting booking" << endl;
        return;
      }

      TH1 *numH = num->getTH1();
      TH1 *denH = den->getTH1();

      if (!numH || !denH) {
        edm::LogWarning("L1TEfficiencyPlotHandler") << (!numH && !denH ? numeratorName + " && " + denominatorName
                                                        : !num         ? numeratorName
                                                                       : denominatorName)
                                                    << " is not TH1F. Quitting booking" << endl;

        return;
      }

      if (numH->GetNbinsX() != denH->GetNbinsX()) {
        edm::LogWarning("L1TEfficiencyPlotHandler") << " # X bins in " << numeratorName << " and " << denominatorName
                                                    << " are different. Quitting booking" << endl;
        return;
      }

      MonitorElement::Kind kind = num->kind();
      bool is1D = kind == MonitorElement::Kind::TH1F || kind == MonitorElement::Kind::TH1D;
      bool is2D = kind == MonitorElement::Kind::TH2F || kind == MonitorElement::Kind::TH2D;

      if (is2D) {
        if (numH->GetNbinsY() != denH->GetNbinsY()) {
          edm::LogWarning("L1TEfficiencyPlotHandler") << " # Y bins in " << numeratorName << " and " << denominatorName
                                                      << " are different. Quitting booking" << endl;
          return;
        }
      }

      ibooker.setCurrentFolder(outputDir_);
      if (is1D) {
        h_efficiency_ = ibooker.book1D(plotName_, den->getTH1F());
      } else if (is2D) {
        h_efficiency_ = ibooker.book2D(plotName_, den->getTH2F());
      }
      h_efficiency_->setEfficiencyFlag();
    }

    void L1TEfficiencyPlotHandler::computeEfficiency(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) {
      if (!h_efficiency_)
        return;

      edm::LogInfo("L1TEfficiencyPlotHandler") << " Computing efficiency for " << plotName_ << endl;

      MonitorElement *num = igetter.get(numeratorDir_ + "/" + plotName_ + numeratorSuffix_);
      MonitorElement *den = igetter.get(denominatorDir_ + "/" + plotName_ + denominatorSuffix_);

      TH1 *numH = num->getTH1();
      TH1 *denH = den->getTH1();
      TH1 *effH = h_efficiency_->getTH1();

      effH->Divide(numH, denH, 1.0, 1.0, "B");
    }

    //___________DQM_analyzer_class________________________________________
    L1TEfficiencyHarvesting::L1TEfficiencyHarvesting(const ParameterSet &ps)
        : verbose_(ps.getUntrackedParameter<bool>("verbose")), plotHandlers_() {
      if (verbose_) {
        edm::LogInfo("L1TEfficiencyHarvesting") << "____________ Storage initialization ____________ " << endl;
      }

      for (const auto &plotConfig : ps.getUntrackedParameter<std::vector<edm::ParameterSet>>("plotCfgs")) {
        vector<string> plots = plotConfig.getUntrackedParameter<vector<string>>("plots");
        for (const auto &plot : plots) {
          plotHandlers_.push_back(L1TEfficiencyPlotHandler(plotConfig, plot));
        }
      }
    }

    //_____________________________________________________________________
    L1TEfficiencyHarvesting::~L1TEfficiencyHarvesting() {}

    //_____________________________________________________________________
    void L1TEfficiencyHarvesting::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) {
      if (verbose_) {
        edm::LogInfo("L1TEfficiencyHarvesting") << "Called endRun." << endl;
      }

      for (auto plotHandler : plotHandlers_) {
        plotHandler.book(ibooker, igetter);
        plotHandler.computeEfficiency(ibooker, igetter);
      }
    }

    //define this as a plug-in
    DEFINE_FWK_MODULE(L1TEfficiencyHarvesting);
  }  // namespace l1t
}  // namespace dqmoffline