DTDataIntegrityTest.cc

Go to the documentation of this file.

/*
 * \file DTDataIntegrityTest.cc
 * 
 * $Date: 2011/07/29 16:23:07 $
 * $Revision: 1.39 $
 * \author S. Bolognesi - CERN
 *
 */

#include <DQM/DTMonitorClient/src/DTDataIntegrityTest.h>

//Framework
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"
#include "CondFormats/DTObjects/interface/DTReadOutMapping.h"
#include "CondFormats/DataRecord/interface/DTReadOutMappingRcd.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <iostream>
#include <string>


using namespace std;
using namespace edm;


DTDataIntegrityTest::DTDataIntegrityTest(const ParameterSet& ps) : nevents(0) {
  
  LogTrace ("DTDQM|DTRawToDigi|DTMonitorClient|DTDataIntegrityTest") << "[DTDataIntegrityTest]: Constructor";

  // prescale on the # of LS to update the test
  prescaleFactor = ps.getUntrackedParameter<int>("diagnosticPrescale", 1);


}


DTDataIntegrityTest::~DTDataIntegrityTest(){

  LogTrace ("DTDQM|DTRawToDigi|DTMonitorClient|DTDataIntegrityTest") << "DataIntegrityTest: analyzed " << nupdates << " updates";

}


void DTDataIntegrityTest::beginJob(){

  LogTrace ("DTDQM|DTRawToDigi|DTMonitorClient|DTDataIntegrityTest") << "[DTDataIntegrityTest]: BeginJob";

  //nSTAEvents = 0;
  nupdates = 0;
  run=0;

  dbe = Service<DQMStore>().operator->();
  
  // book the summary histogram
  dbe->setCurrentFolder("DT/00-DataIntegrity");
  summaryHisto = dbe->book2D("DataIntegritySummary","Summary Data Integrity",12,1,13,5,-2,3);
  summaryHisto->setAxisTitle("Sector",1);
  summaryHisto->setAxisTitle("Wheel",2);

  dbe->setCurrentFolder("DT/00-DataIntegrity");
  summaryTDCHisto = dbe->book2D("DataIntegrityTDCSummary","TDC Summary Data Integrity",12,1,13,5,-2,3);
  summaryTDCHisto->setAxisTitle("Sector",1);
  summaryTDCHisto->setAxisTitle("Wheel",2);

  dbe->setCurrentFolder("DT/00-DataIntegrity");
  glbSummaryHisto = dbe->book2D("DataIntegrityGlbSummary","Summary Data Integrity",12,1,13,5,-2,3);
  glbSummaryHisto->setAxisTitle("Sector",1);
  glbSummaryHisto->setAxisTitle("Wheel",2);

}

void DTDataIntegrityTest::beginRun(const Run& run, const EventSetup& context){

  context.get<DTReadOutMappingRcd>().get(mapping);

}



void DTDataIntegrityTest::beginLuminosityBlock(LuminosityBlock const& lumiSeg, EventSetup const& context) {

  LogTrace ("DTDQM|DTRawToDigi|DTMonitorClient|DTDataIntegrityTest") <<"[DTDataIntegrityTest]: Begin of LS transition";

  // Get the run number
  run = lumiSeg.run();

}



void DTDataIntegrityTest::analyze(const Event& e, const EventSetup& context){
  // count the analyzed events
  nevents++;
  if(nevents%1000 == 0)
    LogTrace ("DTDQM|DTRawToDigi|DTMonitorClient|DTDataIntegrityTest")
      << "[DTDataIntegrityTest]: "<<nevents<<" events";
}



void DTDataIntegrityTest::endLuminosityBlock(LuminosityBlock const& lumiSeg, EventSetup const& context) {

  // counts number of lumiSegs 
  nLumiSegs = lumiSeg.id().luminosityBlock();
  stringstream nLumiSegs_s; nLumiSegs_s << nLumiSegs;
  
  // prescale factor
  if (nLumiSegs%prescaleFactor != 0) return;
  
  LogTrace ("DTDQM|DTRawToDigi|DTMonitorClient|DTDataIntegrityTest")
    <<"[DTDataIntegrityTest]: End of LS " << nLumiSegs << ", performing client operations";


  // counts number of updats 
  nupdates++;

  
  //Counter for x bin in the timing histos
  counter++;

  //Loop on FED id
  for (int dduId=FEDNumbering::MINDTFEDID; dduId<=FEDNumbering::MAXDTFEDID; ++dduId){
    LogTrace ("DTDQM|DTRawToDigi|DTMonitorClient|DTDataIntegrityTest")
      <<"[DTDataIntegrityTest]:FED Id: "<<dduId;
 
    //Each nTimeBin onUpdate remove timing histos and book a new bunch of them
    stringstream dduId_s; dduId_s << dduId;
    
    string histoType;
    
    //Check if the list of ROS is compatible with the channels enabled
    string rosStatusName = "DT/00-DataIntegrity/FED" + dduId_s.str() + "/FED" + dduId_s.str() + "_ROSStatus";
    MonitorElement * FED_ROSStatus = dbe->get(rosStatusName);
     
    // Get the error summary histo
    string fedSummaryName = "DT/00-DataIntegrity/FED" + dduId_s.str() + "_ROSSummary";
    MonitorElement * FED_ROSSummary = dbe->get(fedSummaryName);

    // Get the event lenght plot (used to counr # of processed evts)
    string fedEvLenName = "DT/00-DataIntegrity/FED" + dduId_s.str() + "/FED" + dduId_s.str() + "_EventLenght";
    MonitorElement * FED_EvLenght = dbe->get(fedEvLenName);

    // Get the histos for FED integrity
    string fedIntegrityFolder = "DT/FEDIntegrity/";
    MonitorElement * hFEDEntry = dbe->get(fedIntegrityFolder+"FEDEntries");
    MonitorElement * hFEDFatal = dbe->get(fedIntegrityFolder+"FEDFatal");
    MonitorElement * hFEDNonFatal = dbe->get(fedIntegrityFolder+"FEDNonFatal");

    if(hFEDEntry && hFEDFatal && hFEDNonFatal) {

      if(FED_ROSSummary && FED_ROSStatus && FED_EvLenght) {
    TH2F * histoFEDSummary = FED_ROSSummary->getTH2F();
    TH2F * histoROSStatus  = FED_ROSStatus->getTH2F();
    TH1F * histoEvLenght   = FED_EvLenght->getTH1F();
    // Check that the FED is in the ReadOut using the FEDIntegrity histos
    bool fedNotReadout = (hFEDEntry->getBinContent(dduId-769) == 0 &&
                  hFEDFatal->getBinContent(dduId-769) == 0 &&
                  hFEDNonFatal->getBinContent(dduId-769) == 0);
    int nFEDEvts = histoEvLenght->Integral();
    for(int rosNumber = 1; rosNumber <= 12; ++rosNumber) { // loop on the ROS
      int wheelNumber, sectorNumber;
      if (!readOutToGeometry(dduId,rosNumber,wheelNumber,sectorNumber)) {
        int result = -2;
        float nErrors  = histoFEDSummary->Integral(1,14,rosNumber,rosNumber);
        nErrors += histoROSStatus->Integral(2,8,rosNumber,rosNumber);
        if(nErrors == 0) { // no errors
          result = 0;
        } else { // there are errors
          result = 2;
        }
        summaryHisto->setBinContent(sectorNumber,wheelNumber+3,result);
        int tdcResult = -2;
        float nTDCErrors = histoFEDSummary->Integral(15,15,rosNumber,rosNumber); 
        if(nTDCErrors == 0) { // no errors
          tdcResult = 0;
        } else { // there are errors
          tdcResult = 2;
        }
        summaryTDCHisto->setBinContent(sectorNumber,wheelNumber+3,tdcResult);
        // FIXME: different errors should have different weights
        float sectPerc = max((float)0., ((float)nFEDEvts-nErrors)/(float)nFEDEvts);
        glbSummaryHisto->setBinContent(sectorNumber,wheelNumber+3,sectPerc);
       
        if(fedNotReadout) {
          // no data in this FED: it is off
          summaryHisto->setBinContent(sectorNumber,wheelNumber+3,1);
          summaryTDCHisto->setBinContent(sectorNumber,wheelNumber+3,1);
          glbSummaryHisto->setBinContent(sectorNumber,wheelNumber+3,0);
        }
      }
    }
      
      } else { // no data in this FED: it is off
    for(int rosNumber = 1; rosNumber <= 12; ++rosNumber) {
      int wheelNumber, sectorNumber;
      if (!readOutToGeometry(dduId,rosNumber,wheelNumber,sectorNumber)) {
        summaryHisto->setBinContent(sectorNumber,wheelNumber+3,1);
        summaryTDCHisto->setBinContent(sectorNumber,wheelNumber+3,1);
        glbSummaryHisto->setBinContent(sectorNumber,wheelNumber+3,0);
      } 
    }
      }

    }
    
  }
  
}



void DTDataIntegrityTest::endJob(){

  LogTrace ("DTDQM|DTRawToDigi|DTMonitorClient|DTDataIntegrityTest") <<"[DTDataIntegrityTest] endjob called!";

//   dbe->rmdir("DT/DTDataIntegrity");
}



string DTDataIntegrityTest::getMEName(string histoType, int FEDId){
  //Use the DDU name to find the ME
  stringstream dduID_s; dduID_s << FEDId;

  string folderName = "DT/00-DataIntegrity/FED" + dduID_s.str(); 

  string histoName = folderName + "/FED" + dduID_s.str() + "_" + histoType;
  return histoName;
}



void DTDataIntegrityTest::bookHistos(string histoType, int dduId){
  stringstream dduId_s; dduId_s << dduId;
  dbe->setCurrentFolder("DT/00-DataIntegrity/FED" + dduId_s.str());
  string histoName;

}


int DTDataIntegrityTest::readOutToGeometry(int dduId, int ros, int& wheel, int& sector){

  int dummy;
  return mapping->readOutToGeometry(dduId,ros,2,2,2,wheel,dummy,sector,dummy,dummy,dummy);

}