#include <CastorDataIntegrityMonitor.h>

Inheritance diagram for CastorDataIntegrityMonitor:

Public Member Functions
void	beginRun (const edm::EventSetup &iSetup)

	CastorDataIntegrityMonitor ()

void	cleanup ()

void	processEvent (const FEDRawDataCollection &RawData, const HcalUnpackerReport &report, const CastorElectronicsMap &emap)

void	reset ()

void	setup (const edm::ParameterSet &ps, DQMStore *dbe)

void	unpack (const FEDRawData &raw, const CastorElectronicsMap &emap)

	~CastorDataIntegrityMonitor ()

Public Member Functions inherited from CastorBaseMonitor
	CastorBaseMonitor ()

virtual void	clearME ()

virtual void	done ()

bool	getDiagnostics () const

int	getVerbosity () const

void	setDiagnostics (bool myval)

void	setVerbosity (int verb)

bool	vetoCell (HcalCastorDetId id)

virtual	~CastorBaseMonitor ()

Public Attributes
std::map< uint32_t, std::vector< HcalCastorDetId > >	DCCtoCell

std::map< uint32_t, std::vector< HcalCastorDetId > >::iterator	thisDCC

Private Member Functions
void	fillzoos (int bin, int dccid)

void	mapDCCproblem (int dcc)

void	mapHTRproblem (int dcc, int spigot)

void	UpdateMap ()

Private Attributes
std::map< int, short >	CDFEventType_list

std::map< int, short >::iterator	CDFEvT_it

bool	CDFProbThisDCC

std::map< int, short >::iterator	CDFReservedBits_it

std::map< int, short >	CDFReservedBits_list

std::map< int, short >::iterator	CDFvers_it

std::map< int, short >	CDFversionNumber_list

std::map< int, short >::iterator	DCCEvtFormat_it

std::map< int, short >	DCCEvtFormat_list

int	dccid

std::map< int, short >::iterator	DCCRsvdBits_it

std::map< int, short >	DCCRsvdBits_list

MonitorElement *	fedEntries

MonitorElement *	fedFatal

MonitorElement *	fedNonFatal

std::vector< int >	fedUnpackList_

int	ievt_

MonitorElement *	meDCCVersion

MonitorElement *	meEVT_

std::vector< std::vector< bool > >	problemCASTOR

int	problemsSpigot [15][3]

int	spigotStatus

MonitorElement *	spigotStatusMap

double	statusSpigotDCC

Additional Inherited Members
Protected Attributes inherited from CastorBaseMonitor
std::string	baseFolder_

edm::CPUTimer	cpu_timer

int	fVerbosity

DQMStore *	m_dbe

bool	makeDiagnostics

std::string	rootFolder_

bool	showTiming

Detailed Description

Definition at line 10 of file CastorDataIntegrityMonitor.h.

Constructor & Destructor Documentation

CastorDataIntegrityMonitor::CastorDataIntegrityMonitor ( )

Definition at line 25 of file CastorDataIntegrityMonitor.cc.

25 {

26 }

CastorDataIntegrityMonitor::~CastorDataIntegrityMonitor ( )

Definition at line 32 of file CastorDataIntegrityMonitor.cc.

32 {

33 }

Member Function Documentation

void CastorDataIntegrityMonitor::beginRun ( const edm::EventSetup & iSetup )

Definition at line 88 of file CastorDataIntegrityMonitor.cc.

References CastorBaseMonitor::baseFolder_, DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookInt(), gather_cfg::cout, fedEntries, fedFatal, CastorBaseMonitor::fVerbosity, CastorBaseMonitor::m_dbe, meEVT_, NULL, DQMStore::setCurrentFolder(), and spigotStatusMap.

Referenced by CastorMonitorModule::beginRun().

                                                                      {
 
   if(fVerbosity>0) std::cout << "CastorDataIntegrityMonitor::beginRun (start)" << std::endl;
   
   if ( m_dbe !=NULL ) {
     m_dbe->setCurrentFolder(baseFolder_);
     meEVT_      = m_dbe->bookInt("Digi Task Event Number");
     fedEntries  = m_dbe->book1D("FEDEntries_CASTOR" ,"Number of Entries per CASTOR FED",3,690,693); //-- FED 693 take out for the time being
     fedFatal    = m_dbe->book1D("CASTOR FEDFatal errors"   ,"Number of Fatal Errors per CASTOR FED",3,690,693); //-- FED 693 take out for the time being
     //fedNonFatal = m_dbe->book1D("FEDNonFatal_CASTOR","Number of non-fatal errors CASTOR FED",3,690,693); //-- take it out for the time being
 
     //meDCCVersion =  m_dbe->bookProfile("DCC Firmware Version","DCC Firmware Version", 3, 690, 693, 256, -0.5, 255.5);
     //meDCCVersion->setAxisTitle("FED ID", 1);
     spigotStatusMap = m_dbe->book2D("CASTOR spigot status","CASTOR spigot status",15,0,15,4,690,694); //-- get some space for the legend
 
   }
   
    else{ 
    if(fVerbosity>0) std::cout << "CastorDataIntegrityMonitor::setup - NO DQMStore service" << std::endl; 
   }
  
  
  if(fVerbosity>0) std::cout << "CastorDataIntegrityMonitor::beginRun (end)" << std::endl;
 
 
  return;
 }

void CastorDataIntegrityMonitor::cleanup ( )

Definition at line 47 of file CastorDataIntegrityMonitor.cc.

                                         {
  
 } 

void CastorDataIntegrityMonitor::fillzoos	(	int	bin,
		int	dccid
	)

private

void CastorDataIntegrityMonitor::mapDCCproblem ( int dcc )

private

void CastorDataIntegrityMonitor::mapHTRproblem	(	int	dcc,
		int	spigot
	)

private

void CastorDataIntegrityMonitor::processEvent	(	const FEDRawDataCollection &	RawData,
		const HcalUnpackerReport &	report,
		const CastorElectronicsMap &	emap
	)

— check the spigot arrays

Definition at line 123 of file CastorDataIntegrityMonitor.cc.

References gather_cfg::cout, FEDRawDataCollection::FEDData(), fedUnpackList_, MonitorElement::Fill(), CastorBaseMonitor::fVerbosity, MonitorElement::getTH2F(), i, ievt_, meEVT_, problemsSpigot, FEDRawData::size(), spigotStatusMap, statusSpigotDCC, and unpack().

Referenced by CastorMonitorModule::analyze().

   {
   if(fVerbosity>0) std::cout << "CastorDataIntegrityMonitor::processEvent (begin)" << std::endl;
 
   meEVT_->Fill(ievt_); 
 
   ievt_++;
 
   for (std::vector<int>::const_iterator i=fedUnpackList_.begin();i!=fedUnpackList_.end(); i++)  {
       const FEDRawData& fed = RawData.FEDData(*i);
       if (fed.size()<12) continue;
       unpack(fed,emap);
     }
 
   for (int spigot=0; spigot<15; spigot++){
     for (int dcc=0; dcc<3; dcc++){
 
       if( problemsSpigot[spigot][dcc] == 0)  statusSpigotDCC=1.0;
       else if( double(problemsSpigot[spigot][dcc])/double(ievt_)  < 0.05) statusSpigotDCC=0.;       
       else statusSpigotDCC=-1.0;
       spigotStatusMap->getTH2F()->SetBinContent(spigot+1,dcc+1,statusSpigotDCC);
       if(fVerbosity>1) 
       std::cout<< "==> SpigotNr:"<< spigot+1 <<" DCC_ID:"<< dcc+690 << " # problems=" << problemsSpigot[spigot][dcc]
            << "==> ievt_:"<< ievt_ << " ratio=" << double(problemsSpigot[spigot][dcc])/double(ievt_) << " STATUS=" << statusSpigotDCC << std::endl;
     }
   }
 
   if(fVerbosity>0) std::cout << "CastorDataIntegrityMonitor::processEvent (end)" << std::endl;
 
   return;
 } 

void CastorDataIntegrityMonitor::reset ( void )

Definition at line 39 of file CastorDataIntegrityMonitor.cc.

Referenced by CastorMonitorModule::reset().

                                       {
  
 }

void CastorDataIntegrityMonitor::setup	(	const edm::ParameterSet &	ps,
		DQMStore *	dbe
	)

virtual

Reimplemented from CastorBaseMonitor.

Definition at line 55 of file CastorDataIntegrityMonitor.cc.

References CastorBaseMonitor::baseFolder_, CDFProbThisDCC, cuy::col, gather_cfg::cout, fedUnpackList_, CastorBaseMonitor::fVerbosity, i, ievt_, FEDNumbering::MAXCASTORFEDID, FEDNumbering::MINCASTORFEDID, problemsSpigot, CastorBaseMonitor::rootFolder_, CastorBaseMonitor::setup(), spigotStatus, and statusSpigotDCC.

Referenced by CastorMonitorModule::CastorMonitorModule().

 {
  
  if(fVerbosity>0) std::cout << "CastorDataIntegrityMonitor::setup (start)" << std::endl;
  CastorBaseMonitor::setup(ps,dbe);
  baseFolder_ = rootFolder_+"CastorDataIntegrityMonitor";
  ievt_=0;
  spigotStatus=-99;
  statusSpigotDCC=-99;
  CDFProbThisDCC = false; 
 
 
   for (int row=0; row<15; row++) 
     for (int col=0; col<3; col++)
         problemsSpigot[row][col] = 0;
  
 
 
  for (int i=FEDNumbering::MINCASTORFEDID; i<=FEDNumbering:: MAXCASTORFEDID; ++i)       
       fedUnpackList_.push_back(i);
 
   if(fVerbosity>0) std::cout << "CastorDataIntegrityMonitor::setup (end)" << std::endl;
 
  return;
 
 }

void CastorDataIntegrityMonitor::unpack	(	const FEDRawData &	raw,
		const CastorElectronicsMap &	emap
	)

Definition at line 164 of file CastorDataIntegrityMonitor.cc.

Referenced by processEvent().

                                                                                                {
 
 
   const HcalDCCHeader* dccHeader=(const HcalDCCHeader*)(raw.data());
   if(!dccHeader) return;
 
   unsigned char* trailer_ptr = (unsigned char*) (raw.data()+raw.size()-sizeof(uint64_t));
   FEDTrailer trailer = FEDTrailer(trailer_ptr);
 
   dccid=dccHeader->getSourceId();
 
   HcalHTRData htr;
 
   CDFProbThisDCC = false; 
 
 
   //================ BASIC CHECKS OF DATA INTEGRITY ====================//
 
  
 
   /* 1 */ //There should always be a second CDF header word indicated.
   if (!dccHeader->thereIsASecondCDFHeaderWord()) 
     {
       CDFProbThisDCC = true;  
     }
 
   /* 2 */ //Make sure a reference CDF Version value has been recorded for this dccid
   CDFvers_it = CDFversionNumber_list.find(dccid);
   if (CDFvers_it  == CDFversionNumber_list.end()) 
     {
       CDFversionNumber_list.insert(std::pair<int,short>
                    (dccid,dccHeader->getCDFversionNumber() ) );
       CDFvers_it = CDFversionNumber_list.find(dccid);
     } // then check against it.
 
   if (dccHeader->getCDFversionNumber()!= CDFvers_it->second) 
     {
       CDFProbThisDCC = true;  
     }
   
   /* 3 */ //Make sure a reference CDF EventType value has been recorded for this dccid
   CDFEvT_it = CDFEventType_list.find(dccid);
   if (CDFEvT_it  == CDFEventType_list.end()) 
     {
       CDFEventType_list.insert(std::pair<int,short>
                    (dccid,dccHeader->getCDFEventType() ) );
       CDFEvT_it = CDFEventType_list.find(dccid);
     } // then check against it.
   
   if (dccHeader->getCDFEventType()!= CDFEvT_it->second) 
     {
       // On probation until safe against Orbit Gap Calibration Triggers...
        CDFProbThisDCC = true;  
     }
 
   /* 4 */ //There should always be a '5' in CDF Header word 0, bits [63:60]
   if (dccHeader->BOEshouldBe5Always()!=5) 
     {
       CDFProbThisDCC = true;  
     }
 
   /* 5 */ //There should never be a third CDF Header word indicated.
   if (dccHeader->thereIsAThirdCDFHeaderWord()) 
     {
       CDFProbThisDCC = true;  
     }
 
   /* 6 */ //Make sure a reference value of Reserved Bits has been recorded for this dccid
 
   CDFReservedBits_it = CDFReservedBits_list.find(dccid);
   if (CDFReservedBits_it  == CDFReservedBits_list.end()) {
     CDFReservedBits_list.insert(std::pair<int,short>
                 (dccid,dccHeader->getSlink64ReservedBits() ) );
     CDFReservedBits_it = CDFReservedBits_list.find(dccid);
   } // then check against it.
   
   if ((int) dccHeader->getSlink64ReservedBits()!= CDFReservedBits_it->second) 
     {
     // On probation until safe against Orbit Gap Calibration Triggers...
     //       CDFProbThisDCC = true;
     }
 
   /* 7 */ //There should always be 0x0 in CDF Header word 1, bits [63:60]
   if (dccHeader->BOEshouldBeZeroAlways() !=0) 
     {
       CDFProbThisDCC = true;
     }
   
   /* 8 */ //There should only be one trailer
   if (trailer.moreTrailers()) 
     {
       CDFProbThisDCC = true; 
     }
   //  if trailer.
 
   /* 9 */ //CDF Trailer [55:30] should be the # 64-bit words in the EvFragment
   if ((uint64_t) raw.size() != ( (uint64_t) trailer.lenght()*sizeof(uint64_t)) )  //The function name is a typo! Awesome.
     {
       CDFProbThisDCC = true; 
     }
   /*10 */ //There is a rudimentary sanity check built into the FEDTrailer class
   if (!trailer.check()) 
     {
       CDFProbThisDCC = true; 
     }
 
   if (CDFProbThisDCC) fedFatal->Fill(dccid);
 
   fedEntries->Fill(dccid);
     
 
 
   //================== do similar what HCALRawDataMonitor does   
   //uint64_t* dccfw= (uint64_t*) (raw.data()+(sizeof(uint64_t)*2)); //64-bit DAQ word number 2 (from 0)
   //int dcc_fw =  ( ((*dccfw)>>(6*8))&0x00000000000000FF );         //Shift right 6 bytes, get that low byte.
   //meDCCVersion->Fill(dccid,dcc_fw);
   //char TTS_state = (char)trailer.ttsBits();
 
   //errors per-Spigot bits from the DCC Header
   int WholeErrorList=0; 
 
   for(int spigot=0; spigot<HcalDCCHeader::SPIGOT_COUNT; spigot++) {
 
   if (!( dccHeader->getSpigotEnabled((unsigned int) spigot)) ) continue; //-- skip when not enabled
 
   //-- set OK here 
   spigotStatus=1;
   //-- get DCC id
   dccid=dccHeader->getSourceId();
 
   if (dccid==693) continue; //-- skip this one
 
    WholeErrorList=dccHeader->getLRBErrorBits((unsigned int) spigot);
     if (WholeErrorList!=0) {
 
       if ((WholeErrorList>>0)&0x01) {
     if(fVerbosity>0)  std::cout << "CastorDataIntegrityMonitor:: error==> " << std::endl;
         spigotStatus=-1; 
         ++problemsSpigot[spigot][dccid-690];
       }
 
       if (((WholeErrorList>>1)&0x01)!=0)  {
         if(fVerbosity>0)  std::cout << "CastorDataIntegrityMonitor:: error ==> Uncorrected Error" << std::endl; 
            if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];  
            spigotStatus=-1;
        }
 
       if (((WholeErrorList>>2)&0x01)!=0)  {
     if(fVerbosity>0)  std::cout << "CastorDataIntegrityMonitor:: error ==> Truncated data coming into LRB" << std::endl; 
           if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];  
           spigotStatus=-1;
       }
 
       if (((WholeErrorList>>3)&0x01)!=0)  {
     if(fVerbosity>0)  std::cout << "CastorDataIntegrityMonitor: error ==>: FIFO Overflow" << std::endl; 
         if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];  
         spigotStatus=-1;
       }
 
       if (((WholeErrorList>>4)&0x01)!=0)  {
     if(fVerbosity>0)  std::cout << "CastorDataIntegrityMonitor:: error ==> (EvN Mismatch), htr payload metadeta" << std::endl; 
          if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];  
          spigotStatus=-1;
       }
 
       if (((WholeErrorList>>5)&0x01)!=0)  {
     if(fVerbosity>0)  std::cout << "CastorDataIntegrityMonitor:: error ==> STatus: hdr/data/trlr error" << std::endl; 
         if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];  
         spigotStatus=-1;
       }
 
       if (((WholeErrorList>>6)&0x01)!=0)  {
     if(fVerbosity>0)  std::cout << "CastorDataIntegrityMonitor:: error ==> ODD 16-bit word count from HT error" << std::endl;  
         if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];  
         spigotStatus=-1;
       }
 
     }
 
 
     if (!dccHeader->getSpigotPresent((unsigned int) spigot)){
           if(fVerbosity>0)  std::cout <<"CastorDataIntegrityMonitor:: HTR Problem: Spigot Not Present"<<std::endl;
           if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];  
           spigotStatus=-1;
     }
    else {
       if ( dccHeader->getSpigotDataTruncated((unsigned int) spigot)) {
         if(fVerbosity>0)  std::cout <<"CastorDataIntegrityMonitor:: HTR Problem: Spigot Data Truncated"<<std::endl;
         if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];  
         spigotStatus=-1;
       }
       if ( dccHeader->getSpigotCRCError((unsigned int) spigot)) {
     if(fVerbosity>0)  std::cout <<"CastorDataIntegrityMonitor:: HTR Problem: Spigot CRC Error"<<std::endl; 
         if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];  
         spigotStatus=-1;
       }
      if (dccHeader->getSpigotDataLength(spigot) <(unsigned long)4) {
       if(fVerbosity>0)  std::cout <<"CastorDataIntegrityMonitor:: HTR Problem: Spigot Data Length too small"<<std::endl; 
       if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];  
       spigotStatus=-1;
      }  
 
      if (dccHeader->getSpigotData(spigot,htr,raw.size())==-1) {
     if(fVerbosity>0)  std::cout<< "CastorDataIntegrityMonitor:: Invalid HTR data (data beyond payload size) observed on spigot " << spigot 
          << " of DCC with source id " << dccHeader->getSourceId()<< std::endl;
       if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];
        spigotStatus=-1;
      }
     
     if (!htr.check()) {
       if(fVerbosity>0)  std::cout << "CastorDataIntegrityMonitor:: Invalid HTR data observed on spigot " << spigot << " of DCC with source id " << dccHeader->getSourceId() << std::endl;
      if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];
      spigotStatus=-1;
     }
 
     if (htr.isHistogramEvent()) {
      if(fVerbosity>0)  std::cout << "CastorDataIntegrityMonitor:: Histogram data passed to non-histogram unpacker on spigot " << spigot << " of DCC with source id " 
      << dccHeader->getSourceId() << std::endl; 
      if(spigotStatus>0) ++problemsSpigot[spigot][dccid-690];
      spigotStatus=-1;
     }
   
    }
     
   } //-- end of loop over spigots
 
   return;
 } 