RPCMonitorRaw.cc

Go to the documentation of this file.

#include "DQM/RPCMonitorClient/interface/RPCMonitorRaw.h"
 
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Framework/interface/Event.h"
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include "CondFormats/RPCObjects/interface/LinkBoardElectronicIndex.h"
 
#include "DataFormats/RPCDigi/interface/ReadoutError.h"
#include "DQM/RPCMonitorClient/interface/RPCRawDataCountsHistoMaker.h"
 
#include "TH1F.h"
#include "TH2F.h"
#include "TFile.h"
 
#include <fstream>
#include <iostream>
#include <iomanip>
#include <bitset>
 
typedef std::map<std::pair<int, int>, int>::const_iterator IT;
 
RPCMonitorRaw::RPCMonitorRaw(const edm::ParameterSet& cfg) : theConfig(cfg) {
  rpcRawDataCountsTag_ = consumes<RPCRawDataCounts>(cfg.getParameter<edm::InputTag>("rpcRawDataCountsTag"));
 
  for (unsigned int i = 0; i < 10; i++)
    theWatchedErrorHistoPos[i] = 0;
  std::vector<int> algos = cfg.getUntrackedParameter<std::vector<int> >("watchedErrors");
  for (std::vector<int>::const_iterator it = algos.begin(); it != algos.end(); ++it) {
    unsigned int ialgo = *it;
    if (ialgo < 10)
      theWatchedErrorHistoPos[ialgo] = 1;  // real position initialisain is in begin job. here mark just switched on.
  }
}
 
RPCMonitorRaw::~RPCMonitorRaw() { LogTrace("") << "RPCMonitorRaw destructor"; }
 
void RPCMonitorRaw::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iRun, edm::EventSetup const& iSetup) {
  ibooker.cd();
  ibooker.setCurrentFolder("RPC/LinkMonitor");
 
  me_t[0] = ibooker.book1D("recordType_790", RPCRawDataCountsHistoMaker::emptyRecordTypeHisto(790));
  me_t[1] = ibooker.book1D("recordType_791", RPCRawDataCountsHistoMaker::emptyRecordTypeHisto(791));
  me_t[2] = ibooker.book1D("recordType_792", RPCRawDataCountsHistoMaker::emptyRecordTypeHisto(792));
  for (int i = 0; i < 3; ++i)
    me_t[i]->getTH1F()->SetStats(false);
 
  me_e[0] = ibooker.book1D("readoutErrors_790", RPCRawDataCountsHistoMaker::emptyReadoutErrorHisto(790));
  me_e[1] = ibooker.book1D("readoutErrors_791", RPCRawDataCountsHistoMaker::emptyReadoutErrorHisto(791));
  me_e[2] = ibooker.book1D("readoutErrors_792", RPCRawDataCountsHistoMaker::emptyReadoutErrorHisto(792));
  for (int i = 0; i < 3; ++i)
    me_e[i]->getTH1F()->SetStats(false);
 
  me_mapGoodEvents = ibooker.book2D("mapGoodRecords", "mapGoodRecords", 36, -0.5, 35.5, 3, 789.5, 792.5);
  me_mapGoodEvents->getTH2F()->SetNdivisions(3, "y");
  me_mapGoodEvents->setAxisTitle("rmb");
  me_mapGoodEvents->getTH2F()->SetYTitle("fed");
  me_mapGoodEvents->getTH2F()->SetStats(false);
  me_mapBadEvents = ibooker.book2D("mapErrorRecords", "mapErrorRecords", 36, -0.5, 35.5, 3, 789.5, 792.5);
  me_mapBadEvents->setAxisTitle("fed");
  me_mapBadEvents->getTH2F()->SetYTitle("rmb");
  me_mapBadEvents->getTH2F()->SetNdivisions(3, "y");
  me_mapBadEvents->getTH2F()->SetStats(false);
 
  for (unsigned int i = 0; i <= 9; ++i) {
    if (theWatchedErrorHistoPos[i]) {
      for (unsigned int fed = 790; fed <= 792; ++fed) {
        TH2F* histo = RPCRawDataCountsHistoMaker::emptyReadoutErrorMapHisto(fed, i);
        MonitorElement* watched = ibooker.book2D(histo->GetName(), histo);
        theWatchedErrorHistos[fed - 790].push_back(watched);
        theWatchedErrorHistoPos[i] = theWatchedErrorHistos[fed - 790].size();
      }
    }
  }
}
 
void RPCMonitorRaw::analyze(const edm::Event& ev, const edm::EventSetup& es) {
  edm::Handle<RPCRawDataCounts> rawCounts;
  ev.getByToken(rpcRawDataCountsTag_, rawCounts);
  const RPCRawDataCounts& counts = *rawCounts.product();
 
  //
  // record type
  //
  for (IT it = counts.theRecordTypes.begin(); it != counts.theRecordTypes.end(); ++it)
    me_t[it->first.first - 790]->Fill(it->first.second, it->second);
 
  //
  // good events topology
  //
  for (IT it = counts.theGoodEvents.begin(); it != counts.theGoodEvents.end(); ++it)
    me_mapGoodEvents->Fill(it->first.second, it->first.first, it->second);
 
  //
  // bad events topology
  //
  for (IT it = counts.theBadEvents.begin(); it != counts.theBadEvents.end(); ++it)
    me_mapBadEvents->Fill(it->first.second, it->first.first, it->second);
 
  //
  // readout errors
  //
  for (IT it = counts.theReadoutErrors.begin(); it != counts.theReadoutErrors.end(); ++it) {
    rpcrawtodigi::ReadoutError error(it->first.second);
    LinkBoardElectronicIndex ele = error.where();
    rpcrawtodigi::ReadoutError::ReadoutErrorType type = error.type();
 
    int fed = it->first.first;
    me_e[fed - 790]->Fill(type, it->second);
 
    //
    // in addition fill location map for selected errors
    //
    int idx = theWatchedErrorHistoPos[type] - 1;
    if (idx >= 0) {
      std::vector<MonitorElement*>& wh = theWatchedErrorHistos[fed - 790];
      MonitorElement* me = wh[idx];
      me->Fill(ele.dccInputChannelNum, ele.tbLinkInputNum, it->second);
    }
  }
 
  //  for (int i=0; i<3; ++i) {
  //    me_t[i]->update();
  //    me_e[i]->update();
  //    std::vector<MonitorElement* > & wh = theWatchedErrorHistos[i];
  //    for (std::vector<MonitorElement* >::iterator it=wh.begin(); it != wh.end(); ++it) (*it)->update();
  //  }
  //  me_mapGoodEvents->update();
  //  me_mapBadEvents->update();
}