DTDCSByLumiTask.cc

Go to the documentation of this file.

/*
 * \file DTDCSByLumiTask.cc
 *
 * \author C. Battilana - CIEMAT
 * \author P. Bellan - INFN PD
 * \author A. Branca = INFN PD
 
 *
 */
 
#include <DQM/DTMonitorModule/src/DTDCSByLumiTask.h>
 
// Framework
#include <FWCore/Framework/interface/EventSetup.h>
 
// Geometry
#include "Geometry/DTGeometry/interface/DTLayer.h"
#include "Geometry/DTGeometry/interface/DTTopology.h"
 
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include <FWCore/Framework/interface/EventSetupRecord.h>
#include <FWCore/Framework/interface/EventSetupRecordKey.h>
#include "FWCore/Utilities/interface/Transition.h"
 
#include <iostream>
 
using namespace edm;
using namespace std;
 
DTDCSByLumiTask::DTDCSByLumiTask(const edm::ParameterSet& ps)
    : theEvents(0),
      theLumis(0),
      dtGeometryToken_(esConsumes<DTGeometry, MuonGeometryRecord, edm::Transition::BeginRun>()),
      dtHVStatusToken_(esConsumes<DTHVStatus, DTHVStatusRcd, edm::Transition::EndLuminosityBlock>()) {
  LogTrace("DTDQM|DTMonitorModule|DTDCSByLumiTask") << "[DTDCSByLumiTask]: Constructor" << endl;
 
  // If needed put getParameter here
  // dtDCSByLumiLabel = ps.getParameter<InputTag>("dtDCSByLumiLabel");
}
 
DTDCSByLumiTask::~DTDCSByLumiTask() {
  LogTrace("DTDQM|DTMonitorModule|DTDCSByLumiTask")
      << "DTDCSByLumiTask: processed " << theEvents << " events in " << theLumis << " lumi sections" << endl;
}
 
void DTDCSByLumiTask::dqmBeginRun(const edm::Run& run, const edm::EventSetup& context) {
  LogTrace("DTDQM|DTMonitorModule|DTDCSByLumiTask") << "[DTDCSByLumiTask]: begin run" << endl;
 
  theDTGeom = context.getHandle(dtGeometryToken_);
 
  DTHVRecordFound = true;
 
  eventsetup::EventSetupRecordKey recordKey(eventsetup::EventSetupRecordKey::TypeTag::findType("DTHVStatusRcd"));
 
  std::vector<eventsetup::EventSetupRecordKey> recordKeys;
  context.fillAvailableRecordKeys(recordKeys);
  vector<eventsetup::EventSetupRecordKey>::iterator it = find(recordKeys.begin(), recordKeys.end(), recordKey);
 
  if (it == recordKeys.end()) {
    //record not found
    LogTrace("DTDQM|DTMonitorModule|DTDCSByLumiTask") << "Record DTHVStatusRcd does not exist " << std::endl;
 
    DTHVRecordFound = false;
  }
}
 
void DTDCSByLumiTask::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iRun, edm::EventSetup const&) {
  // Book bylumi histo (# of bins as reduced as possible)
  ibooker.setCurrentFolder(topFolder());
 
  for (int wheel = -2; wheel <= 2; wheel++) {
    stringstream wheel_str;
    wheel_str << wheel;
 
    {
      // Set Lumi scope in order to save histo every LS
      auto scope = DQMStore::IBooker::UseLumiScope(ibooker);
      MonitorElement* ME =
          ibooker.book1D("hActiveUnits" + wheel_str.str(), "Active Untis x LS Wh" + wheel_str.str(), 2, 0.5, 2.5);
      hActiveUnits.push_back(ME);
    }
  }
}
 
void DTDCSByLumiTask::dqmBeginLuminosityBlock(LuminosityBlock const& lumiSeg, EventSetup const&) {
  theLumis++;
 
  LogTrace("DTDQM|DTMonitorModule|DTDCSByLumiTask")
      << "[DTDCSByLumiTask]: Begin of processed lumi # " << lumiSeg.id().luminosityBlock() << " " << theLumis
      << " lumi processed by this job" << endl;
 
  for (int wheel = 0; wheel < 5; wheel++) {
    hActiveUnits[wheel]->Reset();  // Cb by lumi histo need to be resetted in between lumi boundaries
  }
}
 
void DTDCSByLumiTask::dqmEndLuminosityBlock(const edm::LuminosityBlock& lumiSeg, const edm::EventSetup& context) {
  const DTHVStatus* dtHVStatus = nullptr;
  if (DTHVRecordFound) {
    dtHVStatus = &context.getData(dtHVStatusToken_);
  }
 
  vector<const DTLayer*>::const_iterator layersIt = theDTGeom->layers().begin();
  vector<const DTLayer*>::const_iterator layersEnd = theDTGeom->layers().end();
 
  for (; layersIt != layersEnd; ++layersIt) {
    int wheel = (*layersIt)->id().wheel();
 
    int nWiresLayer = (*layersIt)->specificTopology().channels();
    hActiveUnits[wheel + 2]->Fill(1, nWiresLayer);  // CB first bin is # of layers
    int nActiveWires = nWiresLayer;
 
    int flagA = -100;
    int flagC = -100;
    int flagS = -100;
    int first = -100;
    int last = -100;
 
    // CB info is not stored if HV is ON -> in this case get returns 1
    // process all other cases and removed wires with "BAD HV" from active
    // wires list
 
    if (DTHVRecordFound) {
      if (!dtHVStatus->get((*layersIt)->id(), 0, first, last, flagA, flagC, flagS) && (flagA || flagC || flagS)) {
        nActiveWires -= (last - first + 1);
      }
 
      if (!dtHVStatus->get((*layersIt)->id(), 1, first, last, flagA, flagC, flagS) && (flagA || flagC || flagS)) {
        nActiveWires -= (last - first + 1);
      }
    } else {
      nActiveWires = -1.;
    }
 
    hActiveUnits[wheel + 2]->Fill(2, nActiveWires);  // CB 2nd bin is the list of wires wit HV ON
  }
}
 
void DTDCSByLumiTask::analyze(const edm::Event&, const edm::EventSetup&) { theEvents++; }
 
string DTDCSByLumiTask::topFolder() const { return string("DT/EventInfo/DCSContents"); }
 
// Local Variables:
// show-trailing-whitespace: t
// truncate-lines: t
// End: