DTCertificationSummary.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  \author G. Cerminara - INFN Torino
 */

#include "DQM/DTMonitorClient/src/DTCertificationSummary.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"

#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

using namespace std;
using namespace edm;

DTCertificationSummary::DTCertificationSummary(const ParameterSet& pset) {}

DTCertificationSummary::~DTCertificationSummary() {}

void DTCertificationSummary::dqmEndLuminosityBlock(DQMStore::IBooker& ibooker,
                                                   DQMStore::IGetter& igetter,
                                                   edm::LuminosityBlock const& lumiSeg,
                                                   edm::EventSetup const& context) {
  //FR moved here from beginJob

  // book the ME
  ibooker.setCurrentFolder("DT/EventInfo");
  // global fraction
  if (!igetter.get("CertificationSummary")) {
    totalCertFraction = ibooker.bookFloat("CertificationSummary");
    totalCertFraction->Fill(-1);

    // certification map
    certMap = ibooker.book2D("CertificationSummaryMap", "DT Certification Summary Map", 12, 1, 13, 5, -2, 3);
    certMap->setAxisTitle("sector", 1);
    certMap->setAxisTitle("wheel", 2);

    ibooker.setCurrentFolder("DT/EventInfo/CertificationContents");
    // Wheel "fractions" -> will be 0 or 1
    for (int wheel = -2; wheel != 3; ++wheel) {
      stringstream streams;
      streams << "DT_Wheel" << wheel;
      certFractions[wheel] = ibooker.bookFloat(streams.str());
      certFractions[wheel]->Fill(-1);
    }
  }
}

void DTCertificationSummary::beginRun(const Run& run, const EventSetup& setup) {}

void DTCertificationSummary::endRun(const Run& run, const EventSetup& setup) {}

void DTCertificationSummary::dqmEndJob(DQMStore::IBooker& ibooker, DQMStore::IGetter& igetter) {
  // get the relevant summary histos
  MonitorElement* effSummary = igetter.get("DT/05-ChamberEff/EfficiencyGlbSummary");
  MonitorElement* resSummary = igetter.get("DT/02-Segments/ResidualsGlbSummary");
  MonitorElement* segQualSummary = igetter.get("DT/02-Segments/segmentSummary");

  // check that all needed histos are there
  if (effSummary == nullptr || resSummary == nullptr || segQualSummary == nullptr) {
    LogWarning("DQM|DTMonitorClient|DTCertificationSummary")
        << "*** Warning: not all needed summaries are present!" << endl;
    return;
  }

  // reset the MEs
  totalCertFraction->Fill(0.);
  certFractions[-2]->Fill(0.);
  certFractions[-1]->Fill(0.);
  certFractions[-0]->Fill(0.);
  certFractions[1]->Fill(0.);
  certFractions[2]->Fill(0.);
  certMap->Reset();

  // loop over all sectors and wheels
  for (int wheel = -2; wheel != 3; ++wheel) {
    for (int sector = 1; sector != 13; ++sector) {
      double eff = effSummary->getBinContent(sector, wheel + 3);
      double res = resSummary->getBinContent(sector, wheel + 3);
      double segQual = segQualSummary->getBinContent(sector, wheel + 3);

      double total = 0;
      if (segQual != 0) {
        total = min(res, eff);
      } else {
        total = eff;
      }

      certMap->Fill(sector, wheel, total);
      // can use variable weight depending on the sector
      double weight = 1. / 12.;
      certFractions[wheel]->Fill(certFractions[wheel]->getFloatValue() + weight * total);
      double totalWeight = 1. / 60.;
      totalCertFraction->Fill(totalCertFraction->getFloatValue() + totalWeight * total);
    }
  }
}