RPCDeadChannelTest.cc

Go to the documentation of this file.

/* *  \author Anna Cimmino*/
#include "DQM/RPCMonitorDigi/interface/utils.h"
#include <DQM/RPCMonitorClient/interface/RPCDeadChannelTest.h>
// Framework
#include "FWCore/MessageLogger/interface/MessageLogger.h"
// Geometry
#include "Geometry/RPCGeometry/interface/RPCGeomServ.h"
#include "Geometry/RPCGeometry/interface/RPCGeometry.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"

#include <sstream>

RPCDeadChannelTest::RPCDeadChannelTest(const edm::ParameterSet& ps) {
  edm::LogVerbatim("rpcdeadchanneltest") << "[RPCDeadChannelTest]: Constructor";

  useRollInfo_ = ps.getUntrackedParameter<bool>("UseRollInfo", false);

  prescaleFactor_ = ps.getUntrackedParameter<int>("DiagnosticPrescale", 1);
  numberOfDisks_ = ps.getUntrackedParameter<int>("NumberOfEndcapDisks", 4);
  numberOfRings_ = ps.getUntrackedParameter<int>("NumberOfEndcapRings", 2);
}

RPCDeadChannelTest::~RPCDeadChannelTest() {}

void RPCDeadChannelTest::beginJob(std::string& workingFolder) {
  edm::LogVerbatim("rpcdeadchanneltest") << "[RPCDeadChannelTest]: Begin Job";
  globalFolder_ = workingFolder;
}

void RPCDeadChannelTest::getMonitorElements(std::vector<MonitorElement*>& meVector,
                                            std::vector<RPCDetId>& detIdVector,
                                            std::string& clientHistoName) {
  for (unsigned int i = 0; i < meVector.size(); i++) {
    std::string meName = meVector[i]->getName();

    if (meName.find(clientHistoName) != std::string::npos) {
      myOccupancyMe_.push_back(meVector[i]);
      myDetIds_.push_back(detIdVector[i]);
    }
  }
}

void RPCDeadChannelTest::clientOperation() {
  edm::LogVerbatim("rpcdeadchanneltest") << "[RPCDeadChannelTest]:Client Operation";

  MonitorElement* DEAD = nullptr;

  //Loop on chambers
  for (unsigned int i = 0; i < myOccupancyMe_.size(); i++) {
    RPCDetId& detId = myDetIds_[i];
    MonitorElement* myMe = myOccupancyMe_[i];

    if (!myMe)
      continue;

    const QReport* theOccupancyQReport = myMe->getQReport("DeadChannel_0");

    float deadFraction = 0.0;

    if (theOccupancyQReport) {
      float qtresult = theOccupancyQReport->getQTresult();
      // std::vector<dqm::me_util::Channel> badChannels = theOccupancyQReport->getBadChannels();
      deadFraction = 1.0 - qtresult;

    } else {
      int xBins = myMe->getNbinsX();
      float emptyBins = 0.0;
      for (int x = 1; x <= xBins; x++) {
        if (myMe->getBinContent(x) == 0) {
          emptyBins++;
        }
      }
      if (xBins != 0) {
        deadFraction = emptyBins / xBins;
      }
    }

    if (detId.region() == 0)
      DEAD = DEADWheel[detId.ring() + 2];
    else {
      if (-detId.station() + numberOfDisks_ >= 0) {
        if (detId.region() < 0) {
          DEAD = DEADDisk[-detId.station() + numberOfDisks_];
        } else {
          DEAD = DEADDisk[detId.station() + numberOfDisks_ - 1];
        }
      }
    }

    if (DEAD) {
      int xBin, yBin;
      if (detId.region() == 0) {  //Barrel
        xBin = detId.sector();
        rpcdqm::utils rollNumber;
        yBin = rollNumber.detId2RollNr(detId);
      } else {  //Endcap
        //get segment number
        RPCGeomServ RPCServ(detId);
        xBin = RPCServ.segment();
        (numberOfRings_ == 3 ? yBin = detId.ring() * 3 - detId.roll() + 1
                             : yBin = (detId.ring() - 1) * 3 - detId.roll() + 1);
      }
      DEAD->setBinContent(xBin, yBin, deadFraction);
    }

  }  //End loop on rolls in given chambers
}

void RPCDeadChannelTest::myBooker(DQMStore::IBooker& ibooker) {
  ibooker.setCurrentFolder(globalFolder_);

  std::stringstream histoName;

  rpcdqm::utils rpcUtils;

  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
      histoName.str("");
      histoName << "DeadChannelFraction_Roll_vs_Sector_Wheel" << i;
      DEADWheel[i + 2] = ibooker.book2D(histoName.str().c_str(), histoName.str().c_str(), 12, 0.5, 12.5, 21, 0.5, 21.5);

      for (int x = 1; x <= 12; x++) {
        for (int y = 1; y <= 21; y++) {
          DEADWheel[i + 2]->setBinContent(x, y, -1);
        }
      }

      rpcUtils.labelXAxisSector(DEADWheel[i + 2]);
      rpcUtils.labelYAxisRoll(DEADWheel[i + 2], 0, i, useRollInfo_);
    }  //end wheels

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

    int offset = numberOfDisks_;
    if (i > 0) {
      offset--;
    }  //used to skip case equale to zero

    histoName.str("");
    histoName << "DeadChannelFraction_Ring_vs_Segment_Disk" << i;
    DEADDisk[i + offset] = ibooker.book2D(histoName.str().c_str(),
                                          histoName.str().c_str(),
                                          36,
                                          0.5,
                                          36.5,
                                          3 * numberOfRings_,
                                          0.5,
                                          3 * numberOfRings_ + 0.5);

    rpcUtils.labelXAxisSegment(DEADDisk[i + offset]);
    rpcUtils.labelYAxisRing(DEADDisk[i + offset], numberOfRings_, useRollInfo_);

  }  //end loop on wheels and disks
}