RPCMonitorRaw.cc

Go to the documentation of this file.

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


#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.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/RPCRawDataCounts.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_(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::beginJob()
{

// Get DQM interface
  DQMStore* theDMBE = edm::Service<DQMStore>().operator->();
  
  theDMBE->setCurrentFolder("RPC/LinkMonitor");
  
  me_t[0]=theDMBE->book1D("recordType_790",RPCRawDataCountsHistoMaker::emptyRecordTypeHisto(790));
  me_t[1]=theDMBE->book1D("recordType_791",RPCRawDataCountsHistoMaker::emptyRecordTypeHisto(791));
  me_t[2]=theDMBE->book1D("recordType_792",RPCRawDataCountsHistoMaker::emptyRecordTypeHisto(792));
  for (int i=0;i<3;++i)me_t[i]->getTH1F()->SetStats(0);
  
  me_e[0]=theDMBE->book1D("readoutErrors_790",RPCRawDataCountsHistoMaker::emptyReadoutErrorHisto(790));
  me_e[1]=theDMBE->book1D("readoutErrors_791",RPCRawDataCountsHistoMaker::emptyReadoutErrorHisto(791));
  me_e[2]=theDMBE->book1D("readoutErrors_792",RPCRawDataCountsHistoMaker::emptyReadoutErrorHisto(792));
  for (int i=0;i<3;++i)me_e[i]->getTH1F()->SetStats(0);

  me_mapGoodEvents=theDMBE->book2D("mapGoodRecords","mapGoodRecords",36,-0.5,35.5, 3, 789.5,792.5);
  me_mapGoodEvents->getTH2F()->SetNdivisions(3,"y");
  me_mapGoodEvents->getTH2F()->SetXTitle("rmb");
  me_mapGoodEvents->getTH2F()->SetYTitle("fed");
  me_mapGoodEvents->getTH2F()->SetStats(0);
  me_mapBadEvents =theDMBE->book2D("mapErrorRecords", "mapErrorRecords", 36,-0.5,35.5, 3, 789.5,792.5);
  me_mapBadEvents->getTH2F()->SetXTitle("fed");
  me_mapBadEvents->getTH2F()->SetYTitle("rmb");
  me_mapBadEvents->getTH2F()->SetNdivisions(3,"y");
  me_mapBadEvents->getTH2F()->SetStats(0);

  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 = theDMBE->book2D(histo->GetName(),histo);
        theWatchedErrorHistos[fed-790].push_back(watched);
        theWatchedErrorHistoPos[i] = theWatchedErrorHistos[fed-790].size();
      }
    } 
  }

}
void RPCMonitorRaw::endJob()
{
bool writeHistos = theConfig.getUntrackedParameter<bool>("writeHistograms", false);
  if (writeHistos) {
    std::string histoFile = theConfig.getUntrackedParameter<std::string>("histoFileName"); 
    TFile f(histoFile.c_str(),"RECREATE");
    for (int i=0; i<3; ++i) {
      me_t[i]->getTH1F()->Write();
      me_e[i]->getTH1F()->Write();
      std::vector<MonitorElement* > & wh = theWatchedErrorHistos[i];
      for (std::vector<MonitorElement* >::const_iterator it=wh.begin(); it != wh.end(); ++it)(*it)->getTH2F()->Write(); 
    }
    me_mapGoodEvents->getTH2F()->Write();
    me_mapBadEvents->getTH2F()->Write();
    edm::LogInfo(" END JOB, histos saved!");
    f.Close();
  }
}


void RPCMonitorRaw::analyze(const  edm::Event& ev, const edm::EventSetup& es) 
{

  edm::Handle<RPCRawDataCounts> rawCounts;
  ev.getByLabel(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();

 

}