DTOfflineSummaryClients.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  \author M. Pelliccioni - INFN Torino
 *
 *  threadsafe version (//-) oct/nov 2014 - WATWanAbdullah -ncpp-um-my
 *
 */
 
#include "DQM/DTMonitorClient/src/DTOfflineSummaryClients.h"
 
// Framework
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
 
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include <string>
#include <cmath>
 
using namespace edm;
using namespace std;
 
DTOfflineSummaryClients::DTOfflineSummaryClients(const ParameterSet& ps) : nevents(0) {
  LogVerbatim("DTDQM|DTMonitorClient|DTOfflineSummaryClients") << "[DTOfflineSummaryClients]: Constructor";
 
  bookingdone = false;
}
 
DTOfflineSummaryClients::~DTOfflineSummaryClients() {
  LogVerbatim("DTDQM|DTMonitorClient|DTOfflineSummaryClients")
      << "DTOfflineSummaryClients: analyzed " << nevents << " events";
}
 
void DTOfflineSummaryClients::beginRun(const edm::Run& r, const edm::EventSetup& c){};
 
void DTOfflineSummaryClients::dqmEndLuminosityBlock(DQMStore::IBooker& ibooker,
                                                    DQMStore::IGetter& igetter,
                                                    edm::LuminosityBlock const& lumiSeg,
                                                    edm::EventSetup const& context) {
  LogVerbatim("DTDQM|DTMonitorClient|DTOfflineSummaryClients")
      << "[DTOfflineSummaryClients]: End of LS transition" << endl;
}
 
void DTOfflineSummaryClients::dqmEndJob(DQMStore::IBooker& ibooker, DQMStore::IGetter& igetter) {
  LogVerbatim("DTDQM|DTMonitorClient|DTOfflineSummaryClients")
      << "[DTOfflineSummaryClients]: end job. Performin client operation";
 
  // book the summary histos
 
  ibooker.setCurrentFolder("DT/EventInfo");
  summaryReport = ibooker.bookFloat("reportSummary");
  // Initialize to 1 so that no alarms are thrown at the beginning of the run
  summaryReport->Fill(1.);
 
  summaryReportMap = ibooker.book2D("reportSummaryMap", "DT Report Summary Map", 12, 1, 13, 5, -2, 3);
  summaryReportMap->setAxisTitle("sector", 1);
  summaryReportMap->setAxisTitle("wheel", 2);
 
  ibooker.setCurrentFolder("DT/EventInfo/reportSummaryContents");
 
  for (int wheel = -2; wheel != 3; ++wheel) {
    stringstream streams;
    streams << "DT_Wheel" << wheel;
    string meName = streams.str();
 
    theSummaryContents.push_back(ibooker.bookFloat(meName));
    // Initialize to 1 so that no alarms are thrown at the beginning of the run
    theSummaryContents[wheel + 2]->Fill(1.);
  }
 
  // reset the monitor elements
  summaryReportMap->Reset();
  summaryReport->Fill(0.);
  for (int ii = 0; ii != 5; ++ii) {
    theSummaryContents[ii]->Fill(0.);
  }
 
  // Fill the map using, at the moment, only the information from DT chamber efficiency
  // problems at a granularity smaller than the chamber are ignored
  for (int wheel = -2; wheel <= 2; wheel++) {  // loop over wheels
    // retrieve the chamber efficiency summary
    stringstream str;
    str << "DT/05-ChamberEff/EfficiencyMap_All_W" << wheel;
    MonitorElement* segmentWheelSummary = igetter.get(str.str());
    if (segmentWheelSummary != nullptr) {
      float nFailingChambers = 0.;
 
      for (int sector = 1; sector <= 12; sector++) {  // loop over sectors
 
        double meaneff = 0.;
        double errorsum = 0.;
 
        for (int station = 1; station != 5; ++station) {  // loop over stations
 
          const double tmpefficiency = segmentWheelSummary->getBinContent(sector, station);
          const double tmpvariance = pow(segmentWheelSummary->getBinError(sector, station), 2);
 
          if (tmpefficiency == 0 || tmpvariance == 0) {
            nFailingChambers++;
            continue;
          }
 
          meaneff += tmpefficiency / tmpvariance;
          errorsum += 1. / tmpvariance;
 
          if (tmpefficiency < 0.2)
            nFailingChambers++;
 
          LogTrace("DTDQM|DTMonitorClient|DTOfflineSummaryClients")
              << "Wheel: " << wheel << " Stat: " << station << " Sect: " << sector << " status: " << meaneff / errorsum
              << endl;
        }
 
        const double eff_result = meaneff / errorsum;
 
        if (eff_result > 0.7)
          summaryReportMap->Fill(sector, wheel, 1.);
        else if (eff_result < 0.7 && eff_result > 0.5)
          summaryReportMap->Fill(sector, wheel, 0.6);
        else if (eff_result < 0.5 && eff_result > 0.3)
          summaryReportMap->Fill(sector, wheel, 0.4);
        else if (eff_result < 0.3 && eff_result > 0.)
          summaryReportMap->Fill(sector, wheel, 0.15);
      }
      theSummaryContents[wheel + 2]->Fill((48. - nFailingChambers) / 48.);
      summaryReport->Fill(summaryReport->getFloatValue() + theSummaryContents[wheel + 2]->getFloatValue() / 5.);
    } else {
      LogWarning("DTDQM|DTMonitorClient|DTOfflineSummaryClients")
          << " [DTOfflineSummaryClients] Segment Summary not found with name: " << str.str() << endl;
    }
  }
}