ESDaqInfoTask.cc

Go to the documentation of this file.

#include <iostream>

#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"

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

#include <DataFormats/EcalDetId/interface/ESDetId.h>

#include "CondFormats/DataRecord/interface/RunSummaryRcd.h"
#include "CondFormats/RunInfo/interface/RunSummary.h"
#include "CondFormats/RunInfo/interface/RunInfo.h"

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

#include "DataFormats/FEDRawData/interface/FEDNumbering.h"

#include "Geometry/EcalMapping/interface/ESElectronicsMapper.h"

#include "DQM/EcalPreshowerMonitorModule/interface/ESDaqInfoTask.h"

using namespace cms;
using namespace edm;
using namespace std;

ESDaqInfoTask::ESDaqInfoTask(const ParameterSet& ps) {
  usesResource("DQMStore");
  dqmStore_ = Service<DQMStore>().operator->();
  runInfoToken_ = esConsumes<edm::Transition::BeginLuminosityBlock>();
  prefixME_ = ps.getUntrackedParameter<string>("prefixME", "");

  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);

  ESFedRangeMin_ = ps.getUntrackedParameter<int>("ESFedRangeMin", 520);
  ESFedRangeMax_ = ps.getUntrackedParameter<int>("ESFedRangeMax", 575);

  meESDaqFraction_ = nullptr;
  meESDaqActiveMap_ = nullptr;
  meESDaqError_ = nullptr;

  for (int i = 0; i < 56; i++) {
    meESDaqActive_[i] = nullptr;
  }

  if (ps.exists("esMapping")) {
    edm::ParameterSet esMap = ps.getParameter<edm::ParameterSet>("esMapping");
    es_mapping_ = new ESElectronicsMapper(esMap);
  } else {
    edm::LogError("ESDaqInfoTask") << "preshower mapping pointer not initialized. Temporary.";
    es_mapping_ = nullptr;
  }
}

ESDaqInfoTask::~ESDaqInfoTask() { delete es_mapping_; }

void ESDaqInfoTask::beginJob(void) {
  char histo[200];

  if (dqmStore_) {
    dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo");

    sprintf(histo, "DAQSummary");
    meESDaqFraction_ = dqmStore_->bookFloat(histo);
    meESDaqFraction_->Fill(0.0);

    sprintf(histo, "DAQSummaryMap");
    meESDaqActiveMap_ = dqmStore_->book2D(histo, histo, 80, 0.5, 80.5, 80, 0.5, 80.5);
    meESDaqActiveMap_->setAxisTitle("Si X", 1);
    meESDaqActiveMap_->setAxisTitle("Si Y", 2);

    dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo/DAQContents");

    for (int i = 0; i < 56; i++) {
      sprintf(histo, "EcalPreshower_%d", ESFedRangeMin_ + i);
      meESDaqActive_[i] = dqmStore_->bookFloat(histo);
      meESDaqActive_[i]->Fill(0.0);

      ESOnFed_[i] = false;
      for (int x = 0; x < 80; x++) {
        for (int y = 0; y < 80; y++) {
          if (getFEDNumber(x, y) == ESFedRangeMin_ + i) {
            ESOnFed_[i] = true;
            break;
          }
        }
        if (ESOnFed_[i] == true)
          break;
      }
    }

    dqmStore_->setCurrentFolder(prefixME_ + "/ESIntegrityTask");
    sprintf(histo, "DAQError");
    meESDaqError_ = dqmStore_->book1D(histo, histo, 56, ESFedRangeMin_ - 0.5, ESFedRangeMax_ + 0.5);
    meESDaqError_->setAxisTitle("FedID", 1);
  }
}

void ESDaqInfoTask::endJob(void) {}

void ESDaqInfoTask::beginLuminosityBlock(const edm::LuminosityBlock& lumiBlock, const edm::EventSetup& iSetup) {
  this->reset();

  for (int x = 0; x < 80; ++x) {
    for (int y = 0; y < 80; ++y) {
      if (getFEDNumber(x, y) > 0)
        meESDaqActiveMap_->setBinContent(x + 1, y + 1, 0.0);
      else
        meESDaqActiveMap_->setBinContent(x + 1, y + 1, -1.0);
    }
  }

  for (int i = 0; i < 56; i++) {
    if (meESDaqError_)
      meESDaqError_->setBinContent(i, 0.0);
  }

  if (auto runInfoRec = iSetup.tryToGet<RunInfoRcd>()) {
    const auto& sumFED = runInfoRec->getHandle(runInfoToken_);

    std::vector<int> FedsInIds = sumFED->m_fed_in;

    float ESFedCount = 0.;

    for (unsigned int fedItr = 0; fedItr < FedsInIds.size(); ++fedItr) {
      int fedID = FedsInIds[fedItr];

      if (fedID >= ESFedRangeMin_ && fedID <= ESFedRangeMax_) {
        if (ESOnFed_[fedID - ESFedRangeMin_])
          ESFedCount++;

        if (meESDaqActive_[fedID - ESFedRangeMin_])
          meESDaqActive_[fedID - ESFedRangeMin_]->Fill(1.0);

        if (meESDaqActiveMap_) {
          for (int x = 0; x < 80; x++) {
            for (int y = 0; y < 80; y++) {
              if (fedID == getFEDNumber(x, y))
                meESDaqActiveMap_->setBinContent(x + 1, y + 1, 1.0);
            }
          }
        }

        if (meESDaqFraction_)
          meESDaqFraction_->Fill(ESFedCount / 40.);

        if (meESDaqError_) {
          for (int i = 0; i < 56; i++) {
            if (ESOnFed_[fedID - ESFedRangeMin_])
              meESDaqError_->setBinContent(i + 1, 1.0);
            else
              meESDaqError_->setBinContent(i + 1, 2.0);
          }
        }
      }
    }

  } else {
    LogWarning("ESDaqInfoTask") << "Cannot find any RunInfoRcd" << endl;
  }
}

void ESDaqInfoTask::endLuminosityBlock(const edm::LuminosityBlock&, const edm::EventSetup&) {}

void ESDaqInfoTask::reset(void) {
  if (meESDaqFraction_)
    meESDaqFraction_->Reset();

  for (int i = 0; i < 56; i++) {
    if (meESDaqActive_[i])
      meESDaqActive_[i]->Reset();
  }

  if (meESDaqActiveMap_)
    meESDaqActiveMap_->Reset();

  if (meESDaqError_)
    meESDaqError_->Reset();
}

void ESDaqInfoTask::analyze(const Event& e, const EventSetup& c) {}

DEFINE_FWK_MODULE(ESDaqInfoTask);