RPCDCSSummary.cc

Go to the documentation of this file.

#include "DQM/RPCMonitorClient/interface/RPCDCSSummary.h"

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

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

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

RPCDCSSummary::RPCDCSSummary(const edm::ParameterSet& ps) {
  numberOfDisks_ = ps.getUntrackedParameter<int>("NumberOfEndcapDisks", 4);

  FEDRange_.first = ps.getUntrackedParameter<unsigned int>("MinimumRPCFEDId", 790);
  FEDRange_.second = ps.getUntrackedParameter<unsigned int>("MaximumRPCFEDId", 792);

  offlineDQM_ = ps.getUntrackedParameter<bool>("OfflineDQM", true);

  NumberOfFeds_ = FEDRange_.second - FEDRange_.first + 1;
  init_ = false;
  defaultValue_ = 1.;
}

RPCDCSSummary::~RPCDCSSummary() {}

void RPCDCSSummary::beginJob() {}

void RPCDCSSummary::dqmEndLuminosityBlock(DQMStore::IBooker& ibooker,
                                          DQMStore::IGetter& igetter,
                                          edm::LuminosityBlock const& lumiB,
                                          edm::EventSetup const& setup) {
  if (!init_) {
    this->checkDCSbit(setup);
    if (!offlineDQM_) {
      this->myBooker(ibooker);
    }
  }
}

void RPCDCSSummary::dqmEndJob(DQMStore::IBooker& ibooker, DQMStore::IGetter& igetter) {
  if (offlineDQM_) {
    this->myBooker(ibooker);
  }
}

void RPCDCSSummary::checkDCSbit(edm::EventSetup const& setup) {
  defaultValue_ = 1.;

  if (auto runInfoRec = setup.tryToGet<RunInfoRcd>()) {
    defaultValue_ = -1.;
    //get fed summary information
    edm::ESHandle<RunInfo> sumFED;
    runInfoRec->get(sumFED);
    std::vector<int> FedsInIds = sumFED->m_fed_in;
    unsigned int f = 0;
    bool flag = false;
    while (!flag && f < FedsInIds.size()) {
      int fedID = FedsInIds[f];
      //make sure fed id is in allowed range
      if (fedID >= FEDRange_.first && fedID <= FEDRange_.second) {
        defaultValue_ = 1.;
        flag = true;
      }
      f++;
    }
  }
  init_ = true;
}

void RPCDCSSummary::myBooker(DQMStore::IBooker& ibooker) {
  ibooker.setCurrentFolder("RPC/EventInfo");
  // global fraction
  totalDCSFraction = ibooker.bookFloat("DCSSummary");
  totalDCSFraction->Fill(defaultValue_);

  DCSMap_ = ibooker.book2D("DCSSummaryMap", "RPC DCS Summary Map", 15, -7.5, 7.5, 12, 0.5, 12.5);

  //customize the 2d histo
  std::stringstream BinLabel;
  for (int i = 1; i < 13; i++) {
    BinLabel.str("");
    BinLabel << "Sec" << i;
    DCSMap_->setBinLabel(i, BinLabel.str(), 2);
  }

  for (int i = -2; i <= 2; i++) {
    BinLabel.str("");
    BinLabel << "Wheel" << i;
    DCSMap_->setBinLabel((i + 8), BinLabel.str(), 1);
  }

  for (int i = 1; i <= numberOfDisks_; i++) {
    BinLabel.str("");
    BinLabel << "Disk" << i;
    DCSMap_->setBinLabel((i + 11), BinLabel.str(), 1);
    BinLabel.str("");
    BinLabel << "Disk" << -i;
    DCSMap_->setBinLabel((-i + 5), BinLabel.str(), 1);
  }

  //fill the histo with "1" --- just for the moment
  for (int i = 1; i <= 15; i++) {
    for (int j = 1; j <= 12; j++) {
      if (i == 5 || i == 11 || (j > 6 && (i < 6 || i > 10)))
        DCSMap_->setBinContent(i, j, -1);  //bins that not correspond to subdetector parts
      else
        DCSMap_->setBinContent(i, j, defaultValue_);
    }
  }

  if (numberOfDisks_ < 4) {
    for (int j = 1; j <= 12; j++) {
      DCSMap_->setBinContent(1, j, -1);  //bins that not correspond to subdetector parts
      DCSMap_->setBinContent(15, j, -1);
    }
  }

  // book the ME
  ibooker.setCurrentFolder("RPC/EventInfo/DCSContents");

  int limit = numberOfDisks_;
  if (numberOfDisks_ < 2)
    limit = 2;

  for (int i = -1 * limit; i <= limit; i++) {  //loop on wheels and disks
    if (i > -3 && i < 3) {                     //wheels
      std::stringstream streams;
      streams << "RPC_Wheel" << i;
      dcsWheelFractions[i + 2] = ibooker.bookFloat(streams.str());
      dcsWheelFractions[i + 2]->Fill(defaultValue_);
    }

    if (i == 0 || i > numberOfDisks_ || i < (-1 * numberOfDisks_))
      continue;

    int offset = numberOfDisks_;
    if (i > 0)
      offset--;  //used to skip case equale to zero
    std::stringstream streams;
    streams << "RPC_Disk" << i;
    dcsDiskFractions[i + 2] = ibooker.bookFloat(streams.str());
    dcsDiskFractions[i + 2]->Fill(defaultValue_);
  }
}