ESRawDataTask.cc

Go to the documentation of this file.

#include <iostream>
#include <fstream>
#include <vector>

#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"

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

#include "DataFormats/FEDRawData/interface/FEDRawData.h"
#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
#include "DataFormats/FEDRawData/interface/FEDHeader.h"
#include "DataFormats/FEDRawData/interface/FEDTrailer.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"
#include "DataFormats/FEDRawData/src/fed_header.h"

#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalDetId/interface/ESDetId.h"
#include "DataFormats/EcalDigi/interface/ESDataFrame.h"
#include "DataFormats/EcalRawData/interface/ESDCCHeaderBlock.h"
#include "DataFormats/EcalRawData/interface/EcalRawDataCollections.h"

#include "DQM/EcalPreshowerMonitorModule/interface/ESRawDataTask.h"

using namespace cms;
using namespace edm;
using namespace std;

ESRawDataTask::ESRawDataTask(const ParameterSet& ps) {

   init_ = false;

   dqmStore_ = Service<DQMStore>().operator->();

   prefixME_      = ps.getUntrackedParameter<string>("prefixME", "");
   enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
   mergeRuns_     = ps.getUntrackedParameter<bool>("mergeRuns", false);

   FEDRawDataCollection_ = ps.getParameter<InputTag>("FEDRawDataCollection");
   dccCollections_       = ps.getParameter<InputTag>("ESDCCCollections");

}

ESRawDataTask::~ESRawDataTask() {
}

void ESRawDataTask::beginJob(void) {

   ievt_ = 0;

   if ( dqmStore_ ) {
      dqmStore_->setCurrentFolder(prefixME_ + "/ESRawDataTask");
      dqmStore_->rmdir(prefixME_ + "/ESRawDataTask");
   }

}

void ESRawDataTask::beginRun(const Run& r, const EventSetup& c) {

   if ( ! mergeRuns_ ) this->reset();

}

void ESRawDataTask::endRun(const Run& r, const EventSetup& c) {

}

void ESRawDataTask::reset(void) {

}

void ESRawDataTask::setup(void){

   init_ = true;

   char histo[200];

   if (dqmStore_) {
      dqmStore_->setCurrentFolder(prefixME_ + "/ESRawDataTask");

      //sprintf(histo, "ES run number errors");
      //meRunNumberErrors_ = dqmStore_->book1D(histo, histo, 56, 519.5, 575.5); 
      //meRunNumberErrors_->setAxisTitle("ES FED", 1);
      //meRunNumberErrors_->setAxisTitle("Num of Events", 2);

      sprintf(histo, "ES L1A DCC errors");
      meL1ADCCErrors_ = dqmStore_->book1D(histo, histo, 56, 519.5, 575.5); 
      meL1ADCCErrors_->setAxisTitle("ES FED", 1);
      meL1ADCCErrors_->setAxisTitle("Num of Events", 2);

      sprintf(histo, "ES BX DCC errors");
      meBXDCCErrors_ = dqmStore_->book1D(histo, histo, 56, 519.5, 575.5); 
      meBXDCCErrors_->setAxisTitle("ES FED", 1);
      meBXDCCErrors_->setAxisTitle("Num of Events", 2);

      sprintf(histo, "ES Orbit Number DCC errors");
      meOrbitNumberDCCErrors_ = dqmStore_->book1D(histo, histo, 56, 519.5, 575.5); 
      meOrbitNumberDCCErrors_->setAxisTitle("ES FED", 1);
      meOrbitNumberDCCErrors_->setAxisTitle("Num of Events", 2);
      
      sprintf(histo, "Difference between ES and GT L1A");
      meL1ADiff_ = dqmStore_->book1D(histo, histo, 201, -100.5, 100.5);
      meL1ADiff_->setAxisTitle("ES - GT L1A", 1);
      meL1ADiff_->setAxisTitle("Num of Events", 2);

      sprintf(histo, "Difference between ES and GT BX");
      meBXDiff_ = dqmStore_->book1D(histo, histo, 201, -100.5, 100.5);
      meBXDiff_->setAxisTitle("ES - GT BX", 1);
      meBXDiff_->setAxisTitle("Num of Events", 2);

      sprintf(histo, "Difference between ES and GT Orbit Number");
      meOrbitNumberDiff_ = dqmStore_->book1D(histo, histo, 201, -100.5, 100.5);
      meOrbitNumberDiff_->setAxisTitle("ES - GT orbit number", 1);
      meOrbitNumberDiff_->setAxisTitle("Num of Events", 2);
   }

}

void ESRawDataTask::cleanup(void){

   if ( ! init_ ) return;

   if ( dqmStore_ ) {
     //if ( meRunNumberErrors_ ) dqmStore_->removeElement( meRunNumberErrors_->getName() );
     //meRunNumberErrors_ = 0;

     if ( meL1ADCCErrors_ ) dqmStore_->removeElement( meL1ADCCErrors_->getName() );
     meL1ADCCErrors_ = 0;

     if ( meBXDCCErrors_ ) dqmStore_->removeElement( meBXDCCErrors_->getName() );
     meBXDCCErrors_ = 0;

     if ( meOrbitNumberDCCErrors_ ) dqmStore_->removeElement( meOrbitNumberDCCErrors_->getName() );
     meOrbitNumberDCCErrors_ = 0;

     if ( meL1ADiff_ ) dqmStore_->removeElement( meL1ADiff_->getName() );
     meL1ADiff_ = 0;

     if ( meBXDiff_ ) dqmStore_->removeElement( meBXDiff_->getName() );
     meBXDiff_ = 0;

     if ( meOrbitNumberDiff_ ) dqmStore_->removeElement( meOrbitNumberDiff_->getName() );
     meOrbitNumberDiff_ = 0;
   }

   init_ = false;

}

void ESRawDataTask::endJob(void){

   LogInfo("ESRawDataTask") << "analyzed " << ievt_ << " events";

   if ( enableCleanup_ ) this->cleanup();

}

void ESRawDataTask::analyze(const Event& e, const EventSetup& c){

   if ( ! init_ ) this->setup();

   ievt_++;
   runNum_ = e.id().run();

   int gt_L1A = 0, gt_OrbitNumber = 0, gt_BX = 0;
   int esDCC_L1A_MostFreqCounts = 0;
   int esDCC_BX_MostFreqCounts = 0;
   int esDCC_OrbitNumber_MostFreqCounts = 0;

   Handle<ESRawDataCollection> dccs;
   Handle<FEDRawDataCollection> allFedRawData;

   int gtFedDataSize = 0;
   
   if ( e.getByLabel(FEDRawDataCollection_, allFedRawData) ) {
     
     // GT FED data
     const FEDRawData& gtFedData = allFedRawData->FEDData(812);
     
     gtFedDataSize = gtFedData.size()/sizeof(uint64_t);
     
     if ( gtFedDataSize > 0 ) {
       
       FEDHeader header(gtFedData.data());
       
       gt_L1A         = header.lvl1ID();
       gt_OrbitNumber = e.orbitNumber();
       gt_BX          = e.bunchCrossing();
     } else {

       map<int, int> esDCC_L1A_FreqMap;
       map<int, int> esDCC_BX_FreqMap;
       map<int, int> esDCC_OrbitNumber_FreqMap;

       if ( e.getByLabel(dccCollections_, dccs) ) {
     for (ESRawDataCollection::const_iterator dccItr = dccs->begin(); dccItr != dccs->end(); ++dccItr) {
       ESDCCHeaderBlock esdcc = (*dccItr);
       
       esDCC_L1A_FreqMap[esdcc.getLV1()]++;
       esDCC_BX_FreqMap[esdcc.getBX()]++;
       esDCC_OrbitNumber_FreqMap[esdcc.getOrbitNumber()]++;
       
       if (esDCC_L1A_FreqMap[esdcc.getLV1()] > esDCC_L1A_MostFreqCounts) {
         esDCC_L1A_MostFreqCounts = esDCC_L1A_FreqMap[esdcc.getLV1()];
         gt_L1A = esdcc.getLV1();
       } 
       
       if (esDCC_BX_FreqMap[esdcc.getBX()] > esDCC_BX_MostFreqCounts) {
         esDCC_BX_MostFreqCounts = esDCC_BX_FreqMap[esdcc.getBX()];
         gt_BX = esdcc.getBX();
       } 
       
       if (esDCC_OrbitNumber_FreqMap[esdcc.getOrbitNumber()] > esDCC_OrbitNumber_MostFreqCounts) {
         esDCC_OrbitNumber_MostFreqCounts = esDCC_OrbitNumber_FreqMap[esdcc.getOrbitNumber()];
         gt_OrbitNumber = esdcc.getOrbitNumber();
       } 
       
     }
       } else {
     LogWarning("ESRawDataTask") << dccCollections_ << " not available";
       }

     }
   } else {
     LogWarning("ESRawDataTask") << FEDRawDataCollection_ << " not available";
   }

   // DCC 
   vector<int> fiberStatus;
   if ( e.getByLabel(dccCollections_, dccs) ) {
     
     for (ESRawDataCollection::const_iterator dccItr = dccs->begin(); dccItr != dccs->end(); ++dccItr) {
       ESDCCHeaderBlock dcc = (*dccItr);
       
       //if (dcc.getRunNumber() != runNum_) {
       //meRunNumberErrors_->Fill(dcc.fedId());
       //cout<<"Run # err : "<<dcc.getRunNumber()<<" "<<runNum_<<endl;
       //}
       
       if (dcc.getLV1() != gt_L1A) {
     meL1ADCCErrors_->Fill(dcc.fedId());
     //cout<<"L1A err : "<<dcc.getLV1()<<" "<<gt_L1A<<endl;
     Float_t l1a_diff = dcc.getLV1() - gt_L1A;
     if (l1a_diff > 100) l1a_diff = 100;
     else if (l1a_diff < -100) l1a_diff = -100;
     meL1ADiff_->Fill(l1a_diff);
       }
       
       if (dcc.getBX() != gt_BX) {
     meBXDCCErrors_->Fill(dcc.fedId());
     //cout<<"BX err : "<<dcc.getBX()<<" "<<gt_BX<<endl;
     Float_t bx_diff = dcc.getBX() - gt_BX;
     if (bx_diff > 100) bx_diff = 100;
     else if (bx_diff < -100) bx_diff = -100;
     meBXDiff_->Fill(bx_diff);
       }
       if (dcc.getOrbitNumber() != gt_OrbitNumber) {
     meOrbitNumberDCCErrors_->Fill(dcc.fedId());
     //cout<<"Orbit err : "<<dcc.getOrbitNumber()<<" "<<gt_OrbitNumber<<endl;
     Float_t orbitnumber_diff = dcc.getOrbitNumber() - gt_OrbitNumber;
     if (orbitnumber_diff > 100) orbitnumber_diff = 100;
     else if (orbitnumber_diff < -100) orbitnumber_diff = -100;
     meOrbitNumberDiff_->Fill(orbitnumber_diff);
       }
     }
   }

}

DEFINE_FWK_MODULE(ESRawDataTask);