#include <DTDataIntegrityTask.h>

Inheritance diagram for DTDataIntegrityTask:

Public Member Functions
	DTDataIntegrityTask (const edm::ParameterSet &ps, edm::ActivityRegistry &reg)

bool	eventHasErrors () const

void	fedEntry (int dduID)

void	fedFatal (int dduID)

void	fedNonFatal (int dduID)

void	postBeginJob ()

void	postEndJob ()

void	preBeginLumi (const edm::LuminosityBlockID &ls, const edm::Timestamp &iTime)

void	preEndLumi (const edm::LuminosityBlockID &ls, const edm::Timestamp &iTime)

void	preProcessEvent (const edm::EventID &iEvtid, const edm::Timestamp &iTime)

void	processFED (DTDDUData &dduData, const std::vector< DTROS25Data > &rosData, int dduID)

void	processROS25 (DTROS25Data &data, int dduID, int ros)

void	TimeHistos (std::string histoType)

virtual	~DTDataIntegrityTask ()

Public Member Functions inherited from DTDataMonitorInterface
	DTDataMonitorInterface ()

virtual	~DTDataMonitorInterface ()

Private Member Functions
void	bookHistos (const int fedMin, const int fedMax)

void	bookHistos (std::string folder, DTROChainCoding code)

void	bookHistosROS25 (DTROChainCoding code)

void	channelsInCEROS (int cerosId, int chMask, std::vector< int > &channels)

void	channelsInROS (int cerosMask, std::vector< int > &channels)

std::string	topFolder (bool isFEDIntegrity) const

Private Attributes
DTROChainCoding	coding

DQMStore *	dbe

std::map< std::string, std::map< int, MonitorElement * > >	dduHistos

std::map< std::string, std::map< int, DTTimeEvolutionHisto * > >	dduTimeHistos

bool	doTimeHisto

bool	eventErrorFlag

std::set< int >	fedBXIds

std::string	fedIntegrityFolder

bool	getSCInfo

MonitorElement *	hCorruptionSummary

MonitorElement *	hFEDEntry

MonitorElement *	hFEDFatal

MonitorElement *	hFEDNonFatal

MonitorElement *	hTTSSummary

std::multimap< std::string, std::string >::iterator	it

int	mode

int	myPrevEv

int	myPrevFifoVal [7]

int	myPrevRosVal

int	myPrevTtsVal

std::multimap< std::string, std::string >	names

MonitorElement *	nEventMonitor

int	nevents

int	neventsDDU

int	nEventsLS

int	neventsROS25

std::string	outputFile

edm::ParameterSet	parameters

double	rob_max [25]

std::map< std::string, std::map< int, MonitorElement * > >	robHistos

std::map< int, std::set< int > >	rosBxIdsPerFED

std::map< std::string, std::map< int, MonitorElement * > >	rosHistos

std::map< int, std::set< int > >	rosL1AIdsPerFED

std::map< std::string, std::map< int, MonitorElement * > >	rosSHistos

std::map< std::string, std::map< int, DTTimeEvolutionHisto * > >	rosTimeHistos

float	trigger_counter

Detailed Description

Class for DT Data Integrity.

Author: Marco Zanetti (INFN Padova), Gianluca Cerminara (INFN Torino)

Definition at line 37 of file DTDataIntegrityTask.h.

Constructor & Destructor Documentation

DTDataIntegrityTask::DTDataIntegrityTask	(	const edm::ParameterSet &	ps,
		edm::ActivityRegistry &	reg
	)

explicit

Definition at line 35 of file DTDataIntegrityTask.cc.

References edm::hlt::Exception, fedIntegrityFolder, getSCInfo, edm::ParameterSet::getUntrackedParameter(), LogTrace, mode, neventsDDU, neventsROS25, postBeginJob(), preBeginLumi(), preEndLumi(), preProcessEvent(), dtDQMClient_cfg::processingMode, edm::ActivityRegistry::watchPostBeginJob(), edm::ActivityRegistry::watchPreBeginLumi(), edm::ActivityRegistry::watchPreEndLumi(), and edm::ActivityRegistry::watchPreProcessEvent().

                                                                                            : nevents(0) , dbe(0) {
 
   // Register the methods that we want to schedule
   //   reg.watchPostEndJob(this,&DTDataIntegrityTask::postEndJob);
   reg.watchPreBeginLumi(this,&DTDataIntegrityTask::preBeginLumi);
   reg.watchPreEndLumi(this,&DTDataIntegrityTask::preEndLumi);
 
   reg.watchPostBeginJob(this,&DTDataIntegrityTask::postBeginJob);
   reg.watchPreProcessEvent(this,&DTDataIntegrityTask::preProcessEvent);
 
   LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     << "[DTDataIntegrityTask]: Constructor" <<endl;
 
   neventsDDU = 0;
   neventsROS25 = 0;
 
 //   //If you want info VS time histos
 //   doTimeHisto =  ps.getUntrackedParameter<bool>("doTimeHisto", false);
   // Plot quantities about SC
   getSCInfo = ps.getUntrackedParameter<bool>("getSCInfo", false);
 
   fedIntegrityFolder    = ps.getUntrackedParameter<string>("fedIntegrityFolder","DT/FEDIntegrity");
 
   string processingMode = ps.getUntrackedParameter<string>("processingMode","Online");
 
   // processing mode flag to select plots to be produced and basedirs CB vedi se farlo meglio...
   if (processingMode == "Online") {
     mode = 0;
   } else if(processingMode == "SM") {
     mode = 1;
   } else if (processingMode == "Offline") {
     mode = 2;
   } else if (processingMode == "HLT") {
     mode = 3;
   } else {
     throw cms::Exception("MissingParameter")
       << "[DTDataIntegrityTask]: processingMode :" << processingMode
       << " invalid! Must be Online, SM, Offline or HLT !" << endl;
   }
 
 }

DTDataIntegrityTask::~DTDataIntegrityTask ( )

virtual

Definition at line 79 of file DTDataIntegrityTask.cc.

References LogTrace, and neventsDDU.

                                           {
   LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     <<"[DTDataIntegrityTask]: Destructor. Analyzed "<< neventsDDU <<" events"<<endl;
 }

Member Function Documentation

void DTDataIntegrityTask::bookHistos	(	const int	fedMin,
		const int	fedMax
	)

private

Definition at line 104 of file DTDataIntegrityTask.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookFloat(), dbe, hCorruptionSummary, hFEDEntry, hFEDFatal, hFEDNonFatal, hTTSSummary, nEventMonitor, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), and topFolder().

Referenced by bookHistosROS25(), and postBeginJob().

                                                                        {
 
   dbe->setCurrentFolder("DT/EventInfo/Counters");
   nEventMonitor = dbe->bookFloat("nProcessedEventsDataIntegrity");
 
   // Standard FED integrity histos
   dbe->setCurrentFolder(topFolder(true));
 
   int nFED = (fedMax - fedMin)+1;
 
   hFEDEntry = dbe->book1D("FEDEntries","# entries per DT FED",nFED,fedMin,fedMax+1);
   hFEDFatal = dbe->book1D("FEDFatal","# fatal errors DT FED",nFED,fedMin,fedMax+1);
   hFEDNonFatal = dbe->book1D("FEDNonFatal","# NON fatal errors DT FED",nFED,fedMin,fedMax+1);
 
 
   dbe->setCurrentFolder(topFolder(false));
   hTTSSummary = dbe->book2D("TTSSummary","Summary Status TTS",nFED,fedMin,fedMax+1,9,1,10);
   hTTSSummary->setAxisTitle("FED",1);
   hTTSSummary->setBinLabel(1,"ROS PAF",2);
   hTTSSummary->setBinLabel(2,"DDU PAF",2);
   hTTSSummary->setBinLabel(3,"ROS PAF",2);
   hTTSSummary->setBinLabel(4,"DDU PAF",2);
   hTTSSummary->setBinLabel(5,"DDU Full",2);
   hTTSSummary->setBinLabel(6,"L1A Mism.",2);
   hTTSSummary->setBinLabel(7,"ROS Error",2);
   hTTSSummary->setBinLabel(8,"BX Mism.",2);
   hTTSSummary->setBinLabel(9,"DDU Logic Err.",2);
 
   // bookkeeping of the
 
   hCorruptionSummary =  dbe->book2D("DataCorruptionSummary", "Data Corruption Sources",
                    nFED,fedMin,fedMax+1, 8, 1, 9);
   hCorruptionSummary->setAxisTitle("FED",1);
   hCorruptionSummary->setBinLabel(1,"Miss Ch.",2);
   hCorruptionSummary->setBinLabel(2,"ROS BX mism",2);
   hCorruptionSummary->setBinLabel(3,"DDU BX mism",2);
   hCorruptionSummary->setBinLabel(4,"ROS L1A mism",2);
   hCorruptionSummary->setBinLabel(5,"Miss Payload",2);
   hCorruptionSummary->setBinLabel(6,"FCRC bit",2);
   hCorruptionSummary->setBinLabel(7,"Header check",2);
   hCorruptionSummary->setBinLabel(8,"Triler Check",2);
 
 }

void DTDataIntegrityTask::bookHistos	(	std::string	folder,
		DTROChainCoding	code
	)

private

void DTDataIntegrityTask::bookHistosROS25 ( DTROChainCoding code )

private

Definition at line 526 of file DTDataIntegrityTask.cc.

References bookHistos(), getSCInfo, and mode.

Referenced by postBeginJob().

                                                               {
     bookHistos( string("ROS"), code);
 //     for(int robId = 0; robId != 25; ++robId) {
 //       code.setROB(robId);
 //       bookHistos( string("TDCError"), code);
 //     }
     if(mode <= 1)
       if(getSCInfo)
     bookHistos( string("SC"), code);
 }

void DTDataIntegrityTask::channelsInCEROS	(	int	cerosId,
		int	chMask,
		std::vector< int > &	channels
	)

private

Definition at line 1184 of file DTDataIntegrityTask.cc.

Referenced by processROS25().

                                                                                         {
   for (int iCh=0; iCh<6;++iCh) {
     if ((chMask >> iCh) & 0x1){
       channels.push_back(cerosId*6+iCh);
     }
   }
   return;
 }

void DTDataIntegrityTask::channelsInROS	(	int	cerosMask,
		std::vector< int > &	channels
	)

private

Definition at line 1193 of file DTDataIntegrityTask.cc.

Referenced by processROS25().

                                                                            {
   for (int iCeros=0; iCeros<5;++iCeros) {
     if ((cerosMask >> iCeros) & 0x1){
       for (int iCh=0; iCh<6;++iCh) {
     channels.push_back(iCeros*6+iCh);
       }
     }
   }
   return;
 }

bool DTDataIntegrityTask::eventHasErrors ( ) const

Definition at line 1148 of file DTDataIntegrityTask.cc.

References eventErrorFlag.

Referenced by DTDataErrorFilter::hltFilter().

                                                {
   return eventErrorFlag;
 }

void DTDataIntegrityTask::fedEntry ( int dduID )

virtual

Implements DTDataMonitorInterface.

Definition at line 1155 of file DTDataIntegrityTask.cc.

References MonitorElement::Fill(), and hFEDEntry.

                                             {
   hFEDEntry->Fill(dduID);
 }

void DTDataIntegrityTask::fedFatal ( int dduID )

virtual

Implements DTDataMonitorInterface.

Definition at line 1162 of file DTDataIntegrityTask.cc.

References MonitorElement::Fill(), and hFEDFatal.

                                             {
   hFEDFatal->Fill(dduID);
 }

void DTDataIntegrityTask::fedNonFatal ( int dduID )

virtual

Implements DTDataMonitorInterface.

Definition at line 1169 of file DTDataIntegrityTask.cc.

References MonitorElement::Fill(), and hFEDNonFatal.

                                                {
   hFEDNonFatal->Fill(dduID);
 }

void DTDataIntegrityTask::postBeginJob ( )

Definition at line 1262 of file DTDataIntegrityTask.cc.

References bookHistos(), bookHistosROS25(), dbe, LogTrace, FEDNumbering::MAXDTFEDID, FEDNumbering::MINDTFEDID, cppFunctionSkipper::operator, DTROChainCoding::setDDU(), and DTROChainCoding::setROS().

Referenced by DTDataIntegrityTask().

                                        {
   LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask") << "[DTDataIntegrityTask]: postBeginJob" <<endl;
   // get the DQMStore service if needed
   dbe = edm::Service<DQMStore>().operator->();
   LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask") << "[DTDataIntegrityTask] Get DQMStore service" << endl;
 
 
 
   // Loop over the DT FEDs
   int FEDIDmin = FEDNumbering::MINDTFEDID;
   int FEDIDMax = FEDNumbering::MAXDTFEDID;
 
   LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     << " FEDS: " << FEDIDmin  << " to " <<  FEDIDMax << " in the RO" << endl;
 
   // book FED integrity histos
   bookHistos(FEDIDmin, FEDIDMax);
 
   // static booking of the histograms
   for(int fed = FEDIDmin; fed <= FEDIDMax; ++fed) { // loop over the FEDs in the readout
     DTROChainCoding code;
     code.setDDU(fed);
 
     bookHistos( string("ROS_S"), code);
 
     bookHistos( string("DDU"), code);
 
     for(int ros = 1; ros <= 12; ++ros) {// loop over all ROS
       code.setROS(ros);
       bookHistosROS25(code);
     }
   }
 
 }

void DTDataIntegrityTask::postEndJob ( )

Definition at line 95 of file DTDataIntegrityTask.cc.

References LogTrace.

                                     {
   LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     << "[DTDataIntegrityTask]: postEndJob called!" <<endl;
 
 //   if(doTimeHisto) TimeHistos("Event_word_vs_time");
 
 }

void DTDataIntegrityTask::preBeginLumi	(	const edm::LuminosityBlockID &	ls,
		const edm::Timestamp &	iTime
	)

Definition at line 1230 of file DTDataIntegrityTask.cc.

References nEventsLS.

Referenced by DTDataIntegrityTask().

                                                                                                 {
 
   nEventsLS = 0;
 
 }

void DTDataIntegrityTask::preEndLumi	(	const edm::LuminosityBlockID &	ls,
		const edm::Timestamp &	iTime
	)

Definition at line 1236 of file DTDataIntegrityTask.cc.

References dduTimeHistos, edm::LuminosityBlockID::luminosityBlock(), nEventsLS, and rosTimeHistos.

Referenced by DTDataIntegrityTask().

                                                                                               {
 
   int lumiBlock = ls.luminosityBlock();
 
   map<std::string, map<int, DTTimeEvolutionHisto*> >::iterator dduIt  = dduTimeHistos.begin();
   map<std::string, map<int, DTTimeEvolutionHisto*> >::iterator dduEnd = dduTimeHistos.end();
   for(; dduIt!=dduEnd; ++dduIt) {
     map<int, DTTimeEvolutionHisto*>::iterator histoIt  = dduIt->second.begin();
     map<int, DTTimeEvolutionHisto*>::iterator histoEnd = dduIt->second.end();
     for(; histoIt!=histoEnd; ++histoIt) {
       histoIt->second->updateTimeSlot(lumiBlock,nEventsLS);
     }
   }
 
   map<std::string, map<int, DTTimeEvolutionHisto*> >::iterator rosIt  = rosTimeHistos.begin();
   map<std::string, map<int, DTTimeEvolutionHisto*> >::iterator rosEnd = rosTimeHistos.end();
   for(; rosIt!=rosEnd; ++rosIt) {
     map<int, DTTimeEvolutionHisto*>::iterator histoIt  = rosIt->second.begin();
     map<int, DTTimeEvolutionHisto*>::iterator histoEnd = rosIt->second.end();
     for(; histoIt!=histoEnd; ++histoIt) {
       histoIt->second->updateTimeSlot(lumiBlock,nEventsLS);
     }
   }
 
 }

void DTDataIntegrityTask::preProcessEvent	(	const edm::EventID &	iEvtid,
		const edm::Timestamp &	iTime
	)

Definition at line 1204 of file DTDataIntegrityTask.cc.

References eventErrorFlag, fedBXIds, MonitorElement::Fill(), LogTrace, python.multivaluedict::map(), nEventMonitor, nevents, nEventsLS, rosBxIdsPerFED, and rosL1AIdsPerFED.

Referenced by DTDataIntegrityTask().

                                                                                              {
 
   nevents++;
   nEventMonitor->Fill(nevents);
 
   nEventsLS++;
 
   LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask") << "[DTDataIntegrityTask]: preProcessEvent" <<endl;
   // clear the set of BXids from the ROSs
   for(map<int, set<int> >::iterator rosBxIds = rosBxIdsPerFED.begin();
       rosBxIds != rosBxIdsPerFED.end(); ++rosBxIds) {
     (*rosBxIds).second.clear();
   }
 
   fedBXIds.clear();
 
   for(map<int, set<int> >::iterator rosL1AIds = rosL1AIdsPerFED.begin();
       rosL1AIds != rosL1AIdsPerFED.end(); ++rosL1AIds) {
     (*rosL1AIds).second.clear();
   }
 
   // reset the error flag
   eventErrorFlag = false;
 
 }

void DTDataIntegrityTask::processFED	(	DTDDUData &	dduData,
		const std::vector< DTROS25Data > &	rosData,
		int	dduID
	)

virtual

Implements DTDataMonitorInterface.

Definition at line 861 of file DTDataIntegrityTask.cc.

                                                                                                       {
 
   neventsDDU++;
   if (neventsDDU%1000 == 0)
     LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
       << "[DTDataIntegrityTask]: " << neventsDDU << " events analyzed by processFED" << endl;
 
 
   DTROChainCoding code;
   code.setDDU(ddu);
 
   hFEDEntry->Fill(code.getDDUID());
 
   FEDTrailer trailer = data.getDDUTrailer();
   FEDHeader header = data.getDDUHeader();
 
   // check consistency of header and trailer
   if(!header.check()) {
     // error code 7
     hFEDFatal->Fill(code.getDDUID());
     hCorruptionSummary->Fill(code.getDDUID(), 7);
   }
 
   if(!trailer.check()) {
     // error code 8
     hFEDFatal->Fill(code.getDDUID());
     hCorruptionSummary->Fill(code.getDDUID(), 8);
   }
 
   // check CRC error bit set by DAQ before sending data on SLink
   if(data.crcErrorBit()) {
     // error code 6
     hFEDFatal->Fill(code.getDDUID());
     hCorruptionSummary->Fill(code.getDDUID(), 6);
   }
 
   DTDDUSecondStatusWord secondWord = data.getSecondStatusWord();
 
   // Fill the status summary of the TTS
 
 
   //1D HISTO WITH TTS VALUES form trailer (7 bins = 7 values)
   int ttsCodeValue = -1;
   int ttsSummaryBin = -1;
 
   switch(trailer.ttsBits()) {
   case 0:{ //disconnected
     ttsCodeValue = 0;
     break;
   }
   case 1:{ //warning overflow
     ttsCodeValue = 1;
     if(secondWord.warningROSPAF()) { // ROS PAF
       ttsSummaryBin = 1;
     } else { // DDU PAF
       ttsSummaryBin = 2;
     }
 
     break;
   }
   case 2:{ //out of sinch
     ttsCodeValue = 2;
     bool knownOrigin = false;
     if(secondWord.outOfSynchROSError()) {// ROS Error
       ttsSummaryBin = 7;
       knownOrigin = true;
     }
     if(secondWord.l1AIDError()) {// L1A Mism.
       ttsSummaryBin = 6;
       knownOrigin = true;
     }
     if(secondWord.bxIDError()) {// BX Mism.
       ttsSummaryBin = 8;
       knownOrigin = true;
     }
     if(secondWord.outputFifoFull() || secondWord.inputFifoFull() || secondWord.fifoFull()) { // DDU Full
       ttsSummaryBin = 5;
       knownOrigin = true;
     }
     if(!knownOrigin) ttsSummaryBin = 9; // Error in DDU logic
 
     break;
   }
   case 4:{ //busy
     ttsCodeValue = 3;
     bool knownOrigin = false;
     if(secondWord.busyROSPAF()) { // ROS PAF
       ttsSummaryBin = 3;
       knownOrigin = true;
     }
     if(secondWord.outputFifoAlmostFull() || secondWord.inputFifoAlmostFull() || secondWord.fifoAlmostFull() ){ // DDU PAF
       ttsSummaryBin = 4;
       knownOrigin = true;
     }
     if(!knownOrigin) ttsSummaryBin = 9; // Error in DDU logic
     break;
   }
   case 8:{ //ready
     ttsCodeValue = 4;
     break;
   }
   case 12:{ //error
     ttsCodeValue = 5;
     break;
   }
   case 16:{ //disconnected
     ttsCodeValue = 6;
     break;
   }
   default:{
     LogError("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
       <<"[DTDataIntegrityTask] DDU control: wrong TTS value "<<trailer.ttsBits()<<endl;
     ttsCodeValue = 7;
   }
   }
   if(mode <= 1) dduHistos["TTSValues"][code.getDDUID()]->Fill(ttsCodeValue);
   if(ttsSummaryBin != -1) {
     hTTSSummary->Fill(ddu, ttsSummaryBin);
   }
 
 
 
 
 
 
   //2D HISTO: ROS VS STATUS (8 BIT = 8 BIN) from 1st-2nd status words (9th BIN FROM LIST OF ROS in 2nd status word)
   MonitorElement* hROSStatus = dduHistos["ROSStatus"][code.getDDUID()];
   //1D HISTO: NUMBER OF ROS IN THE EVENTS from 2nd status word
 
   int rosList = secondWord.rosList();
   set<int> rosPositions;
   for(int i=0;i<12;i++) {
     if(rosList & 0x1) {
       rosPositions.insert(i);
       //9th BIN FROM LIST OF ROS in 2nd status word
       if(mode <= 2) hROSStatus->Fill(8,i,1);
     }
     rosList >>= 1;
   }
 
   int channel=0;
   for (vector<DTDDUFirstStatusWord>::const_iterator fsw_it = data.getFirstStatusWord().begin();
        fsw_it != data.getFirstStatusWord().end(); fsw_it++) {
     // assuming association one-to-one between DDU channel and ROS
     if(mode <= 2) {
       hROSStatus->Fill(0,channel,(*fsw_it).channelEnabled());
       hROSStatus->Fill(1,channel,(*fsw_it).timeout());
       hROSStatus->Fill(2,channel,(*fsw_it).eventTrailerLost());
       hROSStatus->Fill(3,channel,(*fsw_it).opticalFiberSignalLost());
       hROSStatus->Fill(4,channel,(*fsw_it).tlkPropagationError());
       hROSStatus->Fill(5,channel,(*fsw_it).tlkPatternError());
       hROSStatus->Fill(6,channel,(*fsw_it).tlkSignalLost());
       hROSStatus->Fill(7,channel,(*fsw_it).errorFromROS());
     }
     // check that the enabled channel was also in the read-out
     if((*fsw_it).channelEnabled() == 1 &&
        rosPositions.find(channel) == rosPositions.end()) {
       if(mode <= 2) hROSStatus->Fill(9,channel,1);
       // error code 1
       hFEDFatal->Fill(code.getDDUID());
       hCorruptionSummary->Fill(code.getDDUID(), 1);
     }
     channel++;
   }
 
 
   // ---------------------------------------------------------------------
   // cross checks between FED and ROS data
   // check the BX ID against the ROSs
   set<int> rosBXIds = rosBxIdsPerFED[ddu];
   if((rosBXIds.size() > 1 || rosBXIds.find(header.bxID()) == rosBXIds.end()) && rosBXIds.size() != 0) { // in this case look for faulty ROSs
     for(vector<DTROS25Data>::const_iterator rosControlData = rosData.begin();
     rosControlData != rosData.end(); ++rosControlData) { // loop over the ROS data
       for (vector<DTROSDebugWord>::const_iterator debug_it = (*rosControlData).getROSDebugs().begin();
        debug_it != (*rosControlData).getROSDebugs().end(); debug_it++) { // Loop over ROS debug words
     if ((*debug_it).debugType() == 0 && (*debug_it).debugMessage() != header.bxID()) { // check the BX
       int ros = (*rosControlData).getROSID();
       // fill the error bin
       if(mode <= 2) hROSStatus->Fill(11,ros-1);
       // error code 2
       hFEDFatal->Fill(code.getDDUID());
       hCorruptionSummary->Fill(code.getDDUID(), 2);
     }
       }
     }
   }
 
   // check the BX ID against other FEDs
   fedBXIds.insert(header.bxID());
   if(fedBXIds.size() != 1) {
     LogWarning("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
       << "ERROR: FED " << ddu << " BX ID different from other feds: " << header.bxID() << endl;
     // error code 3
     hFEDFatal->Fill(code.getDDUID());
     hCorruptionSummary->Fill(code.getDDUID(), 3);
   }
 
 
   // check the L1A ID against the ROSs
   set<int> rosL1AIds = rosL1AIdsPerFED[ddu];
   if((rosL1AIds.size() > 1 || rosL1AIds.find(header.lvl1ID()-1) == rosL1AIds.end()) && rosL1AIds.size() != 0) { // in this case look for faulty ROSs
     //If L1A_ID error identify which ROS has wrong L1A
     for (vector<DTROS25Data>::const_iterator rosControlData = rosData.begin();
      rosControlData != rosData.end(); rosControlData++) { // loop over the ROS data
       int ROSHeader_TTCCount = ((*rosControlData).getROSHeader().TTCEventCounter() + 1) % 0x1000000; // fix comparison in case of last counting bin in ROS /first one in DDU
       if( ROSHeader_TTCCount != header.lvl1ID() ) {
     int ros = (*rosControlData).getROSID();
     if(mode <= 2) hROSStatus->Fill(10,ros-1);
     // error code 4
     hFEDFatal->Fill(code.getDDUID());
     hCorruptionSummary->Fill(code.getDDUID(), 4);
       }
     }
   }
 
   //1D HISTOS: EVENT LENGHT from trailer
   int fedEvtLenght = trailer.lenght()*8;
   //   if(fedEvtLenght > 16000) fedEvtLenght = 16000; // overflow bin
   dduHistos["EventLenght"][code.getDDUID()]->Fill(fedEvtLenght);
 
   if(mode > 1) return;
 
   dduTimeHistos["FEDAvgEvLenghtvsLumi"][code.getDDUID()]->accumulateValueTimeSlot(fedEvtLenght);
 
   // size of the list of ROS in the Read-Out
   dduHistos["ROSList"][code.getDDUID()]->Fill(rosPositions.size());
 
 
   //2D HISTO: FIFO STATUS from 2nd status word
   MonitorElement *hFIFOStatus = dduHistos["FIFOStatus"][code.getDDUID()];
   int inputFifoFull = secondWord.inputFifoFull();
   int inputFifoAlmostFull = secondWord.inputFifoAlmostFull();
   int fifoFull = secondWord.fifoFull();
   int fifoAlmostFull = secondWord.fifoAlmostFull();
   int outputFifoFull = secondWord.outputFifoFull();
   int outputFifoAlmostFull = secondWord.outputFifoAlmostFull();
   for(int i=0;i<3;i++){
     if(inputFifoFull & 0x1){
       hFIFOStatus->Fill(i,0);
     }
     if(inputFifoAlmostFull & 0x1){
       hFIFOStatus->Fill(i,1);
     }
     if(fifoFull & 0x1){
       hFIFOStatus->Fill(3+i,0);
     }
     if(fifoAlmostFull & 0x1){
       hFIFOStatus->Fill(3+i,1);
     }
     if(!(inputFifoFull & 0x1) && !(inputFifoAlmostFull & 0x1)){
       hFIFOStatus->Fill(i,2);
     }
     if(!(fifoFull & 0x1) && !(fifoAlmostFull & 0x1)){
       hFIFOStatus->Fill(3+i,2);
     }
     inputFifoFull >>= 1;
     inputFifoAlmostFull >>= 1;
     fifoFull >>= 1;
     fifoAlmostFull >>= 1;
   }
 
   if(outputFifoFull){
     hFIFOStatus->Fill(6,0);
   }
   if(outputFifoAlmostFull){
     hFIFOStatus->Fill(6,1);
   }
   if(!outputFifoFull && !outputFifoAlmostFull){
     hFIFOStatus->Fill(6,2);
   }
 
 
 
 
 
 
 
   //1D HISTO: EVENT TYPE from header
   dduHistos["EventType"][code.getDDUID()]->Fill(header.triggerType());
 
   // fill the distribution of the BX ids
   dduHistos["BXID"][code.getDDUID()]->Fill(header.bxID());
 
 
 }

void DTDataIntegrityTask::processROS25	(	DTROS25Data &	data,
		int	dduID,
		int	ros
	)

virtual

Implements DTDataMonitorInterface.

Definition at line 538 of file DTDataIntegrityTask.cc.

References DTROBHeaderWord::bunchID(), channelsInCEROS(), channelsInROS(), eventErrorFlag, DTROSTrailerWord::EventWordCount(), MonitorElement::Fill(), FirstRos, DTROChainCoding::getDDUID(), DTROChainCoding::getROB(), DTROS25Data::getROBHeaders(), DTROS25Data::getROBTrailers(), DTROChainCoding::getROS(), DTROS25Data::getROSDebugs(), DTROS25Data::getROSErrors(), DTROS25Data::getROSHeader(), DTROChainCoding::getROSID(), DTROS25Data::getROSTrailer(), DTROChainCoding::getSCID(), getSCInfo, DTROS25Data::getSCPrivHeader(), DTROS25Data::getSCTrailer(), DTROS25Data::getTDCData(), DTROS25Data::getTDCError(), DTROSTrailerWord::l1AFifoOccupancy(), LogTrace, mode, nevents, neventsROS25, DTLocalTriggerSectorCollectorHeaderWord::NumberOf16bitWords(), DTTDCMeasurementWord::PC(), ResetCount_unfolded, DTROBHeaderWord::robID(), rosBxIdsPerFED, rosHistos, rosL1AIdsPerFED, rosSHistos, rosTimeHistos, ROSWords_t(), DTROChainCoding::setDDU(), DTROChainCoding::setROB(), DTROChainCoding::setROS(), DTROSTrailerWord::TPX(), DTROSHeaderWord::TTCEventCounter(), and DTLocalTriggerTrailerWord::wordCount().

                                                                            {
   neventsROS25++; // FIXME: implement a counter which makes sense
 
 //   if (neventsROS25%1000 == 0)
       LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     << "[DTDataIntegrityTask]: " << neventsROS25 << " events analyzed by processROS25" << endl;
 
   // The ID of the RO board (used to map the histos)
   DTROChainCoding code;
   code.setDDU(ddu);
   code.setROS(ros);
 
   MonitorElement* ROSSummary = rosSHistos["ROSSummary"][code.getDDUID()];
 
   // Summary of all ROB errors
   MonitorElement* ROSError = 0;
   if(mode <= 2) ROSError = rosHistos["ROSError"][code.getROSID()];
 
   if ( (mode<=2) && (!ROSError) ) {
     LogError("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask") <<
     "Trying to access non existing ME at ROSID " << code.getROSID() <<
     std::endl;
     return;
   }
 
   // L1A ids to be checked against FED one
   rosL1AIdsPerFED[ddu].insert(data.getROSHeader().TTCEventCounter());
 
   // ROS errors
 
 
   // check for TPX errors
   if (data.getROSTrailer().TPX() != 0) {
     LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask") << " TXP error en ROS "
                                       << code.getROS() << endl;
     ROSSummary->Fill(9,code.getROS());
   }
 
   // L1 Buffer almost full (non-critical error!)
   if(data.getROSTrailer().l1AFifoOccupancy() > 31) {
      ROSSummary->Fill(10,code.getROS());
    }
 
   // FIXME: what is this about???
   if (neventsROS25 == 1) FirstRos = code.getROSID();
   if (code.getROSID() == FirstRos) nevents++ ;
 
 
   for (vector<DTROSErrorWord>::const_iterator error_it = data.getROSErrors().begin();
        error_it != data.getROSErrors().end(); error_it++) { // Loop over ROS error words
 
     LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask") << " Error in ROS " << code.getROS()
                                       << " ROB Id " << (*error_it).robID()
                                       << " Error type " << (*error_it).errorType() << endl;
 
     // Fill the ROSSummary (1 per FED) histo
     ROSSummary->Fill((*error_it).errorType(), code.getROS());
     if((*error_it).errorType() <= 11) { // set error flag
        eventErrorFlag = true;
     }
 
     if(mode <= 2) {
       // Fill the ROB Summary (1 per ROS) histo
       if ((*error_it).robID() != 31) {
     ROSError->Fill((*error_it).errorType(),(*error_it).robID());
       }
       else if ((*error_it).errorType() == 4) {
     vector<int> channelBins;
     channelsInROS((*error_it).cerosID(),channelBins);
     vector<int>::const_iterator channelIt  = channelBins.begin();
     vector<int>::const_iterator channelEnd = channelBins.end();
     for(;channelIt!=channelEnd;++channelIt) {
       ROSError->Fill(4,(*channelIt));
     }
       }
     }
   }
 
 
   int ROSDebug_BunchNumber = -1;
   int ROSDebug_BcntResCntLow = 0;
   int ROSDebug_BcntResCntHigh = 0;
   int ROSDebug_BcntResCnt = 0;
 
   for (vector<DTROSDebugWord>::const_iterator debug_it = data.getROSDebugs().begin();
        debug_it != data.getROSDebugs().end(); debug_it++) { // Loop over ROS debug words
 
     int debugROSSummary = 0;
     int debugROSError   = 0;
     vector<int> debugBins;
     bool hasEvIdMis = false;
     vector<int> evIdMisBins;
 
     if ((*debug_it).debugType() == 0 ) {
       ROSDebug_BunchNumber = (*debug_it).debugMessage();
     } else if ((*debug_it).debugType() == 1 ) {
       ROSDebug_BcntResCntLow = (*debug_it).debugMessage();
     } else if ((*debug_it).debugType() == 2 ) {
       ROSDebug_BcntResCntHigh = (*debug_it).debugMessage();
     } else if ((*debug_it).debugType() == 3) {
       if ((*debug_it).dontRead()){
     debugROSSummary = 11;
     debugROSError   = 8;
     if (mode <= 2) channelsInCEROS((*debug_it).cerosIdCerosStatus(),(*debug_it).dontRead(),debugBins);
       } if ((*debug_it).evIdMis()){
     hasEvIdMis = true;
     if (mode <= 2) channelsInCEROS((*debug_it).cerosIdCerosStatus(),(*debug_it).evIdMis(),evIdMisBins);
       }
     } else if ((*debug_it).debugType() == 4 &&
            (*debug_it).cerosIdRosStatus()){
       debugROSSummary = 13;
       debugROSError   = 10;
       if (mode <= 2) channelsInROS((*debug_it).cerosIdRosStatus(),debugBins);
     }
 
     if (debugROSSummary) {
       ROSSummary->Fill(debugROSSummary,code.getROS());
       if (mode <= 2) {
     vector<int>::const_iterator channelIt  = debugBins.begin();
     vector<int>::const_iterator channelEnd = debugBins.end();
     for (;channelIt!=channelEnd;++channelIt) {
       ROSError->Fill(debugROSError,(*channelIt));
     }
       }
     }
 
     if (hasEvIdMis) {
       ROSSummary->Fill(12,code.getROS());
       if (mode <= 2) {
     vector<int>::const_iterator channelIt  = evIdMisBins.begin();
     vector<int>::const_iterator channelEnd = evIdMisBins.end();
     for (;channelIt!=channelEnd;++channelIt) {
       ROSError->Fill(9,(*channelIt));
     }
       }
     }
 
   }
 
   ROSDebug_BcntResCnt = (ROSDebug_BcntResCntHigh << 15) + ROSDebug_BcntResCntLow;
   //   LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
   //     << " ROS: " << code.getROS() << " ROSDebug_BunchNumber " << ROSDebug_BunchNumber
   //     << " ROSDebug_BcntResCnt " << ROSDebug_BcntResCnt << endl;
 
 
   //     Event words vs time
   // FIXME: what is this doing???
   ROSWords_t(ResetCount_unfolded,code.getROS(),ROSDebug_BcntResCnt,nevents);
 
   // fill hists it here
   //   histoType = "Event_word_vs_time";
   //   if (rosHistos[histoType].find(code.getROSID()) != rosHistos[histoType].end()){
   //   (rosHistos.find(histoType)->second).find(code.getROSID())->second->
   //        Fill((ResetCount_unfolded),data.getROSTrailer().EventWordCount());
   //   (rosHistos.find(histoType)->second).find(code.getROSID())->second->setAxisTitle("Time(s)",1);
   //    }
   //   else {
   //      (rosHistos.find(histoType)->second).find(code.getROSID())->second->
   //            Fill((ResetCount_unfolded),data.getROSTrailer().EventWordCount());}
 
 
 
 
   // ROB Group Header
   // Check the BX of the ROB headers against the BX of the ROS
   for (vector<DTROBHeader>::const_iterator rob_it = data.getROBHeaders().begin();
        rob_it != data.getROBHeaders().end(); rob_it++) { // loop over ROB headers
 
     code.setROB((*rob_it).first);
     DTROBHeaderWord robheader = (*rob_it).second;
 
     rosBxIdsPerFED[ddu].insert(ROSDebug_BunchNumber);
 
     if (robheader.bunchID() != ROSDebug_BunchNumber) {
       // fill ROS Summary plot
       ROSSummary->Fill(8,code.getROS());
       eventErrorFlag = true;
 
       // fill ROB Summary plot for that particular ROS
       if(mode <= 2) ROSError->Fill(7,robheader.robID());
     }
   }
 
 
   if(mode <= 1) { // produce only when not in HLT
     // ROB Trailer
     for (vector<DTROBTrailerWord>::const_iterator robt_it = data.getROBTrailers().begin();
      robt_it != data.getROBTrailers().end(); robt_it++) { // loop over ROB trailers
       float  wCount = (*robt_it).wordCount()<100. ? (*robt_it).wordCount() : 99.9;
       rosHistos["ROB_mean"][code.getROSID()]->Fill((*robt_it).robID(),wCount);
     }
 
 //     // Trigger frequency
 //     double frequency = 0;
 //     // FIXME: how is the frequency computed
 //     ROS_L1A_Frequency(code.getROS(),ROSDebug_BcntResCnt,neventsROS25,frequency,trigger_counter);
 //     rosHistos["Trigger_frequency"][code.getROSID()]->Fill(frequency);
 
     // Plot the event lenght //NOHLT
     int rosEventLenght = data.getROSTrailer().EventWordCount()*4;
     rosTimeHistos["ROSAvgEventLenghtvsLumi"][code.getROSID()]->accumulateValueTimeSlot(rosEventLenght);
 
     if(rosEventLenght > 1600) rosEventLenght = 1600;
     rosHistos["ROSEventLenght"][code.getROSID()]->Fill(rosEventLenght);
   }
 
 
   // TDC Data
   for (vector<DTTDCData>::const_iterator tdc_it = data.getTDCData().begin();
        tdc_it != data.getTDCData().end(); tdc_it++) { // loop over TDC data
 
     DTTDCMeasurementWord tdcDatum = (*tdc_it).second;
 
     if ( tdcDatum.PC() !=0)  {
       LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     << " PC error in ROS " << code.getROS() << " TDC " << (*tdc_it).first << endl;
       //     fill ROS Summary plot
       ROSSummary->Fill(7,code.getROS());
 
       eventErrorFlag = true;
 
       // fill ROB Summary plot for that particular ROS
       if(mode <= 2) ROSError->Fill(6,(*tdc_it).first);
     }
   }
 
   // TDC Error
   for (vector<DTTDCError>::const_iterator tdc_it = data.getTDCError().begin();
        tdc_it != data.getTDCError().end(); tdc_it++) { // loop over TDC errors
 
     code.setROB((*tdc_it).first);
 
     int tdcError_ROSSummary = 0;
     int tdcError_ROSError = 0;
     int tdcError_TDCHisto = 0;
 
     if(((*tdc_it).second).tdcError() & 0x4000 ) {
       LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     << " ROS " << code.getROS() << " ROB " << code.getROB()
     << " Internal fatal Error 4000 in TDC " << (*tdc_it).first << endl;
 
       tdcError_ROSSummary = 14;
       tdcError_ROSError   = 11;
       tdcError_TDCHisto   = 0;
 
     } else if ( ((*tdc_it).second).tdcError() & 0x0249 ) {
       LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     << " ROS " << code.getROS() << " ROB " << code.getROB()
     << " TDC FIFO overflow in TDC " << (*tdc_it).first << endl;
 
       tdcError_ROSSummary = 15;
       tdcError_ROSError   = 12;
       tdcError_TDCHisto   = 1;
 
     } else if ( ((*tdc_it).second).tdcError() & 0x0492 ) {
       LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     << " ROS " << code.getROS() << " ROB " << code.getROB()
     << " TDC L1 buffer overflow in TDC " << (*tdc_it).first << endl;
 
       tdcError_ROSSummary = 16;
       tdcError_ROSError   = 13;
       tdcError_TDCHisto   = 2;
 
     } else if ( ((*tdc_it).second).tdcError() & 0x2000 ) {
       LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     << " ROS " << code.getROS() << " ROB " << code.getROB()
     << " TDC L1A FIFO overflow in TDC " << (*tdc_it).first << endl;
 
       tdcError_ROSSummary = 17;
       tdcError_ROSError   = 14;
       tdcError_TDCHisto   = 3;
 
     } else if ( ((*tdc_it).second).tdcError() & 0x0924 ) {
       LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     << " ROS " << code.getROS() << " ROB " << code.getROB()
     << " TDC hit error in TDC " << (*tdc_it).first << endl;
 
       tdcError_ROSSummary = 18;
       tdcError_ROSError   = 15;
       tdcError_TDCHisto   = 4;
 
     } else if ( ((*tdc_it).second).tdcError() & 0x1000 ) {
       LogTrace("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     << " ROS " << code.getROS() << " ROB " << code.getROB()
     << " TDC hit rejected in TDC " << (*tdc_it).first << endl;
 
       tdcError_ROSSummary = 19;
       tdcError_ROSError   = 16;
       tdcError_TDCHisto   = 5;
 
     } else {
       LogWarning("DTRawToDigi|DTDQM|DTMonitorModule|DTDataIntegrityTask")
     << " TDC error code not known " << ((*tdc_it).second).tdcError() << endl;
     }
 
     ROSSummary->Fill(tdcError_ROSSummary,code.getROS());
 
     if(tdcError_ROSSummary <= 15) {
       eventErrorFlag = true;
     }
 
     if(mode <= 2) {
       ROSError->Fill(tdcError_ROSError,(*tdc_it).first);
       if(mode <= 1)
     rosHistos["TDCError"][code.getROSID()]->Fill(tdcError_TDCHisto+6*((*tdc_it).second).tdcID(),(*tdc_it).first);
     }
   }
 
   // Read SC data
   if (mode <= 1 && getSCInfo) {
     // SC Data
 
     // NumberOf16bitWords counts the # of words + 1 subheader
     // the SC includes the SC "private header" and the ROS header and trailer (= NumberOf16bitWords +3)
     rosHistos["SCSizeVsROSSize"][code.getSCID()]->Fill(ros,data.getSCPrivHeader().NumberOf16bitWords()+3-data.getSCTrailer().wordCount());
 
   }
 
 
 
 
 }

void DTDataIntegrityTask::TimeHistos ( std::string histoType )

Definition at line 490 of file DTDataIntegrityTask.cc.

References newFWLiteAna::bin, dbe, DQMStore::get(), MonitorElement::getBinContent(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), it, j, names, and MonitorElement::setAxisRange().

                                                     {
 
  if(histoType == "Event_word_vs_time"){
 
   for (it = names.begin(); it != names.end(); it++) {
 
     if ((*it).first==histoType){
 
      MonitorElement * h1 =dbe->get((*it).second);
 
  int first_bin = -1, last_bin=-1;
    for( int bin=1; bin < h1->getNbinsX()+1; bin++ ){
     for( int j=1; j < h1->getNbinsY(); j++ ){
      if( h1->getBinContent(bin,j) > 0 ) {
       if( first_bin == -1 ) { first_bin = bin; }
       last_bin = bin;
    }
   }
  }
 
   if( first_bin > 1 ) { first_bin -= 1; }
   if( last_bin < h1-> getNbinsX() ){ last_bin += 1; }
     h1->setAxisRange(0,last_bin,1);
    }
   }
  }
 }

std::string DTDataIntegrityTask::topFolder ( bool isFEDIntegrity ) const

private

Definition at line 1173 of file DTDataIntegrityTask.cc.

References fedIntegrityFolder, and mode.

Referenced by bookHistos().

                                                                   {
 
   string folder = isFEDIntegrity ? fedIntegrityFolder : "DT/00-DataIntegrity";
 
   if (!isFEDIntegrity)
     folder += (mode==1) ? "_SM/" : (mode==3) ? "_EvF/" : "/";
 
   return folder;
 
 }