ESRawDataTask Class Reference

#include <ESRawDataTask.h>

Inheritance diagram for ESRawDataTask:

Public Member Functions
	ESRawDataTask (const edm::ParameterSet &ps)
	~ESRawDataTask () override
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
	Analyze. More...
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	Setup. More...
void	endJob (void) override
	EndJob. More...

Private Attributes
edm::EDGetTokenT< ESRawDataCollection >	dccCollections_
edm::EDGetTokenT< FEDRawDataCollection >	FEDRawDataCollection_
int	ievt_
MonitorElement *	meBXDCCErrors_
MonitorElement *	meBXDiff_
MonitorElement *	meL1ADCCErrors_
MonitorElement *	meL1ADiff_
MonitorElement *	meOrbitNumberDCCErrors_
MonitorElement *	meOrbitNumberDiff_
std::string	prefixME_
int	runNum_

Detailed Description

Definition at line 15 of file ESRawDataTask.h.

Constructor & Destructor Documentation

ESRawDataTask::ESRawDataTask

(

const edm::ParameterSet &

ps

)

Definition at line 29 of file ESRawDataTask.cc.

References edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

                                                   {

   prefixME_      = ps.getUntrackedParameter<string>("prefixME", "");

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

   ievt_ = 0;
}

ESRawDataTask::~ESRawDataTask ( )

inlineoverride

Definition at line 20 of file ESRawDataTask.h.

References analyze(), bookHistograms(), EnergyCorrector::c, MillePedeFileConverter_cfg::e, and endJob().

{}

Member Function Documentation

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

overrideprotected

Analyze.

Definition at line 87 of file ESRawDataTask.cc.

References edm::SortedCollection< T, SORT >::begin(), edm::EventBase::bunchCrossing(), FEDRawData::data(), DEFINE_FWK_MODULE, edm::SortedCollection< T, SORT >::end(), FEDRawDataCollection::FEDData(), ESDCCHeaderBlock::fedId(), ESDCCHeaderBlock::getBX(), edm::Event::getByToken(), ESDCCHeaderBlock::getLV1(), ESDCCHeaderBlock::getOrbitNumber(), RecoTauValidation_cfi::header, edm::EventBase::id(), edm::EventBase::orbitNumber(), edm::EventID::run(), and FEDRawData::size().

Referenced by ~ESRawDataTask().

                                                              {

   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.getByToken(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.getByToken(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.getByToken(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);
       }
     }
   }

}

void ESRawDataTask::bookHistograms ( DQMStore::IBooker &  iBooker, edm::Run const &  , edm::EventSetup const &    )

overrideprotected

Setup.

Definition at line 39 of file ESRawDataTask.cc.

References DQMStore::IBooker::book1D(), trackerHits::histo, MonitorElement::setAxisTitle(), and DQMStore::IBooker::setCurrentFolder().

Referenced by ~ESRawDataTask().

{
   char histo[200];

   iBooker.setCurrentFolder(prefixME_ + "/ESRawDataTask");

   //sprintf(histo, "ES run number errors");
   //meRunNumberErrors_ = iBooker.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_ = iBooker.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_ = iBooker.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_ = iBooker.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_ = iBooker.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_ = iBooker.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_ = iBooker.book1D(histo, histo, 201, -100.5, 100.5);
   meOrbitNumberDiff_->setAxisTitle("ES - GT orbit number", 1);
   meOrbitNumberDiff_->setAxisTitle("Num of Events", 2);
}

void ESRawDataTask::endJob ( void  )

overrideprotected

EndJob.

Definition at line 81 of file ESRawDataTask.cc.

Referenced by ~ESRawDataTask().

                              {

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

}

Member Data Documentation

edm::EDGetTokenT<ESRawDataCollection> ESRawDataTask::dccCollections_

private

Definition at line 39 of file ESRawDataTask.h.

edm::EDGetTokenT<FEDRawDataCollection> ESRawDataTask::FEDRawDataCollection_

private

Definition at line 40 of file ESRawDataTask.h.

int ESRawDataTask::ievt_

private

Definition at line 35 of file ESRawDataTask.h.

MonitorElement* ESRawDataTask::meBXDCCErrors_

private

Definition at line 44 of file ESRawDataTask.h.

MonitorElement* ESRawDataTask::meBXDiff_

private

Definition at line 47 of file ESRawDataTask.h.

MonitorElement* ESRawDataTask::meL1ADCCErrors_

private

Definition at line 43 of file ESRawDataTask.h.

MonitorElement* ESRawDataTask::meL1ADiff_

private

Definition at line 46 of file ESRawDataTask.h.

MonitorElement* ESRawDataTask::meOrbitNumberDCCErrors_

private

Definition at line 45 of file ESRawDataTask.h.

MonitorElement* ESRawDataTask::meOrbitNumberDiff_

private

Definition at line 48 of file ESRawDataTask.h.

std::string ESRawDataTask::prefixME_

private

Definition at line 37 of file ESRawDataTask.h.

int ESRawDataTask::runNum_

private

Definition at line 50 of file ESRawDataTask.h.