#include <RawDataTask.h>

Inheritance diagram for ecaldqm::RawDataTask:

Public Types
enum	Constants { nEventTypes = 25 }

Public Types inherited from ecaldqm::DQWorkerTask
typedef EcalDCCHeaderBlock::EcalDCCEventSettings	EventSettings

Public Member Functions
void	addDependencies (DependencySet &) override

bool	analyze (void const *, Collections) override

void	beginEvent (edm::Event const &, edm::EventSetup const &) override

void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	beginRun (edm::Run const &, edm::EventSetup const &) override

	RawDataTask ()

void	runOnRawData (EcalRawDataCollection const &)

void	runOnSource (FEDRawDataCollection const &)

	~RawDataTask ()

Public Member Functions inherited from ecaldqm::DQWorkerTask
	DQWorkerTask ()

virtual void	endEvent (edm::Event const &, edm::EventSetup const &)

virtual bool	filterRunType (short const *)

virtual bool	filterTrigger (edm::TriggerResultsByName const &)

void	recoverStats ()

virtual void	setTokens (edm::ConsumesCollector &)

void	softReset ()

virtual	~DQWorkerTask ()

Public Member Functions inherited from ecaldqm::DQWorker
virtual void	bookMEs (DQMStore::IBooker &)

	DQWorker ()

virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	endRun (edm::Run const &, edm::EventSetup const &)

std::string const &	getName () const

bool	onlineMode () const

virtual void	releaseMEs ()

void	setEventNumber (edm::EventNumber_t _e)

void	setLumiNumber (edm::LuminosityBlockNumber_t _l)

void	setRunNumber (edm::RunNumber_t _r)

void	setTime (time_t _t)

virtual	~DQWorker ()

Private Attributes
int	bx_

int	feL1Offset_

int	l1A_

int	orbit_

edm::RunNumber_t	runNumber_

short	triggerType_

Additional Inherited Members
Static Public Member Functions inherited from ecaldqm::DQWorkerTask
static void	fillDescriptions (edm::ParameterSetDescription &)

Static Public Member Functions inherited from ecaldqm::DQWorker
static void	fillDescriptions (edm::ParameterSetDescription &_desc)

Protected Member Functions inherited from ecaldqm::DQWorkerTask
void	setME (edm::ParameterSet const &) final

Protected Member Functions inherited from ecaldqm::DQWorker
void	initialize (std::string const &_name, edm::ParameterSet const &)

void	print_ (std::string const &, int=0) const

virtual void	setParams (edm::ParameterSet const &)

virtual void	setSource (edm::ParameterSet const &)

void	setVerbosity (int _verbosity)

Protected Attributes inherited from ecaldqm::DQWorkerTask
std::set< std::string >	resettable_

Protected Attributes inherited from ecaldqm::DQWorker
bool	booked_

MESetCollection	MEs_

std::string	name_

bool	onlineMode_

Timestamp	timestamp_

int	verbosity_

bool	willConvertToEDM_

Detailed Description

Definition at line 13 of file RawDataTask.h.

Member Enumeration Documentation

enum ecaldqm::RawDataTask::Constants

Enumerator
nEventTypes

Definition at line 29 of file RawDataTask.h.

                    {
       nEventTypes = 25
     };

Constructor & Destructor Documentation

ecaldqm::RawDataTask::RawDataTask ( )

Definition at line 14 of file RawDataTask.cc.

                            :
     DQWorkerTask(),
     runNumber_(0),
     l1A_(0),
     orbit_(0),
     bx_(0),
     triggerType_(0),
     feL1Offset_(0)
   {
   }

ecaldqm::RawDataTask::~RawDataTask ( )

inline

Definition at line 16 of file RawDataTask.h.

16 {}

Member Function Documentation

void ecaldqm::RawDataTask::addDependencies ( DependencySet & _dependencies )

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 26 of file RawDataTask.cc.

References ecaldqm::kEcalRawData, ecaldqm::kSource, and ecaldqm::DependencySet::push_back().

   {
     _dependencies.push_back(Dependency(kEcalRawData, kSource));
   }

bool ecaldqm::RawDataTask::analyze	(	void const *	_p,
		Collections	_collection
	)

inlineoverridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 43 of file RawDataTask.h.

References ecaldqm::kEcalRawData, ecaldqm::kSource, runOnRawData(), and runOnSource().

                                                                          {
     switch(_collection){
     case kSource:
       if(_p) runOnSource(*static_cast<FEDRawDataCollection const*>(_p));
       return true;
       break;
     case kEcalRawData:
       if(_p) runOnRawData(*static_cast<EcalRawDataCollection const*>(_p));
       return true;
       break;
     default:
       break;
     }
     return false;
   }

void ecaldqm::RawDataTask::beginEvent	(	edm::Event const &	_evt,
		edm::EventSetup const &
	)

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 45 of file RawDataTask.cc.

References edm::EventBase::bunchCrossing(), bx_, edm::EventBase::experimentType(), feL1Offset_, edm::EventBase::isRealData(), l1A_, orbit_, edm::EventBase::orbitNumber(), and triggerType_.

   {
     orbit_ = _evt.orbitNumber() & 0xffffffff;
     bx_ = _evt.bunchCrossing() & 0xfff;
     triggerType_ = _evt.experimentType() & 0xf;
     l1A_ = 0;
     feL1Offset_ = _evt.isRealData() ? 1 : 0;
   }

void ecaldqm::RawDataTask::beginLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

overridevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 38 of file RawDataTask.cc.

References ecaldqm::DQWorker::MEs_.

   {
     MEs_.at("DesyncByLumi").reset();
     MEs_.at("FEByLumi").reset();
   }

void ecaldqm::RawDataTask::beginRun	(	edm::Run const &	_run,
		edm::EventSetup const &
	)

overridevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 32 of file RawDataTask.cc.

References edm::RunBase::run(), and runNumber_.

   {
     runNumber_ = _run.run();
   }

void ecaldqm::RawDataTask::runOnRawData ( EcalRawDataCollection const & _dcchs )

Definition at line 77 of file RawDataTask.cc.

References edm::SortedCollection< T, SORT >::begin(), ecaldqm::BlockSize, bx_, ecaldqm::BXDesync, ecaldqm::ccuExists(), ecaldqm::ChannelId, ecaldqm::dccId(), edm::SortedCollection< T, SORT >::end(), feL1Offset_, ecaldqm::getElectronicsMap(), ecaldqm::HeaderError, ecaldqm::HParity, ecaldqm::DQWorker::Timestamp::iLumi, ecaldqm::kEEpLow, l1A_, ecaldqm::L1ABXDesync, ecaldqm::L1ADesync, ecaldqm::LinkError, ecaldqm::DQWorker::MEs_, orbit_, runNumber_, fileinputsource_cfi::runType, ntuplemaker::status, ecaldqm::Timeout, ecaldqm::DQWorker::timestamp_, triggerType_, and ecaldqm::VParity.

Referenced by analyze().

   {
     using namespace std;
 
     MESet& meRunNumber(MEs_.at("RunNumber"));
     MESet& meOrbit(MEs_.at("Orbit"));
     MESet& meOrbitDiff(MEs_.at("OrbitDiff"));
     MESet& meTriggerType(MEs_.at("TriggerType"));
     MESet& meL1ADCC(MEs_.at("L1ADCC"));
     MESet& meBXDCC(MEs_.at("BXDCC"));
     MESet& meBXDCCDiff(MEs_.at("BXDCCDiff"));
     MESet& meBXFE(MEs_.at("BXFE"));
     MESet& meBXFEDiff(MEs_.at("BXFEDiff"));
     MESet& meBXFEInvalid(MEs_.at("BXFEInvalid"));
     MESet& meL1AFE(MEs_.at("L1AFE"));
     MESet& meFEStatus(MEs_.at("FEStatus"));
     MESet& meDesyncByLumi(MEs_.at("DesyncByLumi"));
     MESet& meDesyncTotal(MEs_.at("DesyncTotal"));
     MESet& meFEByLumi(MEs_.at("FEByLumi"));
     MESet& meBXTCC(MEs_.at("BXTCC"));
     MESet& meL1ATCC(MEs_.at("L1ATCC"));
     MESet& meBXSRP(MEs_.at("BXSRP"));
     MESet& meL1ASRP(MEs_.at("L1ASRP"));
     MESet& meTrendNSyncErrors(MEs_.at("L1ATCC"));
     MESet& meEventTypePreCalib(MEs_.at("EventTypePreCalib"));
     MESet& meEventTypeCalib(MEs_.at("EventTypeCalib"));
     MESet& meEventTypePostCalib(MEs_.at("EventTypePostCalib"));
 
     if(!l1A_){
       // majority vote on L1A.. is there no better implementation?
       map<int, int> l1aCounts;
       for(EcalRawDataCollection::const_iterator dcchItr(_dcchs.begin()); dcchItr != _dcchs.end(); ++dcchItr){
         l1aCounts[dcchItr->getLV1()]++;
       }
       int maxVote(0);
       for(map<int, int>::iterator l1aItr(l1aCounts.begin()); l1aItr != l1aCounts.end(); ++l1aItr){
         if(l1aItr->second > maxVote){
           maxVote = l1aItr->second;
           l1A_ = l1aItr->first;
         }
       }
     }
 
     for(EcalRawDataCollection::const_iterator dcchItr(_dcchs.begin()); dcchItr != _dcchs.end(); ++dcchItr){
       int dccId(dcchItr->id());
 
       int dccL1A(dcchItr->getLV1());
       short dccL1AShort(dccL1A & 0xfff);
       int dccBX(dcchItr->getBX());
 
       meOrbitDiff.fill(dccId, dcchItr->getOrbit() - orbit_);
       meBXDCCDiff.fill(dccId, dccBX - bx_);
       if(dccBX == -1) meBXFEInvalid.fill(dccId, 68.5);
 
       if(dcchItr->getRunNumber() != int(runNumber_)) meRunNumber.fill(dccId);
       if(dcchItr->getOrbit() != orbit_) meOrbit.fill(dccId);
       if(dcchItr->getBasicTriggerType() != triggerType_) meTriggerType.fill(dccId);
       if(dccL1A != l1A_) meL1ADCC.fill(dccId);
       if(dccBX != bx_) meBXDCC.fill(dccId);
 
       const vector<short> &feStatus(dcchItr->getFEStatus());
       const vector<short> &feBxs(dcchItr->getFEBxs());
       const vector<short> &feL1s(dcchItr->getFELv1());
 
       double feDesync(0.);
       double statusError(0.);
 
       for(unsigned iFE(0); iFE < feStatus.size(); iFE++){
         if(!ccuExists(dccId, iFE + 1)) continue;
 
         short status(feStatus[iFE]);
 
         if(feBxs[iFE] != -1 && dccBX != -1){
           meBXFEDiff.fill(dccId, feBxs[iFE] - dccBX);
         }
         if(feBxs[iFE] == -1) meBXFEInvalid.fill(dccId, iFE + 0.5);
 
         if(status != BXDesync && status != L1ABXDesync){ // BX desync not detected in the DCC
           if(feBxs[iFE] != dccBX && feBxs[iFE] != -1 && dccBX != -1){
             meBXFE.fill(dccId, iFE + 0.5);
             feDesync += 1.;
           }
         }
 
         if(status != L1ADesync && status != L1ABXDesync){
           if(feL1s[iFE] + feL1Offset_ != dccL1AShort && feL1s[iFE] != -1 && dccL1AShort != 0){
             meL1AFE.fill(dccId, iFE + 0.5);
             feDesync += 1.;
           }
         }
 
         if(iFE >= 68) continue;
 
         DetId id(getElectronicsMap()->dccTowerConstituents(dccId, iFE + 1).at(0));
         meFEStatus.fill(id, status);
 
         switch(status){
         case Timeout:
         case HeaderError:
         case ChannelId:
         case LinkError:
         case BlockSize:
         case L1ADesync:
         case BXDesync:
         case L1ABXDesync:
         case HParity:
         case VParity:
           statusError += 1.;
           break;
         default:
           continue;
         }
       }
 
       if(feDesync > 0.){
         meDesyncByLumi.fill(dccId, feDesync);
         meDesyncTotal.fill(dccId, feDesync);
         meTrendNSyncErrors.fill(double(timestamp_.iLumi), feDesync);
       }
       if(statusError > 0.)
         meFEByLumi.fill(dccId, statusError);
 
       const vector<short> &tccBx(dcchItr->getTCCBx());
       const vector<short> &tccL1(dcchItr->getTCCLv1());
 
       if(tccBx.size() == 4){ // EB uses tccBx[0]; EE uses all
         if(dccId <= kEEmHigh + 1 || dccId >= kEEpLow + 1){
           for(int iTCC(0); iTCC < 4; iTCC++){
 
             if(tccBx[iTCC] != dccBX && tccBx[iTCC] != -1 && dccBX != -1)
               meBXTCC.fill(dccId);
 
             if(tccL1[iTCC] != dccL1AShort && tccL1[iTCC] != -1 && dccL1AShort != 0)
               meL1ATCC.fill(dccId);
 
           }
         }else{
 
           if(tccBx[0] != dccBX && tccBx[0] != -1 && dccBX != -1)
             meBXTCC.fill(dccId);
 
           if(tccL1[0] != dccL1AShort && tccL1[0] != -1 && dccL1AShort != 0)
             meL1ATCC.fill(dccId);
 
         }
       }
 
       short srpBx(dcchItr->getSRPBx());
       short srpL1(dcchItr->getSRPLv1());
 
       if(srpBx != dccBX && srpBx != -1 && dccBX != -1)
         meBXSRP.fill(dccId);
 
       if(srpL1 != dccL1AShort && srpL1 != -1 && dccL1AShort != 0)
         meL1ASRP.fill(dccId);
 
       const int calibBX(3490);
 
       short runType(dcchItr->getRunType() + 1);
       if(runType < 0 || runType > 22) runType = 0;
       if(dccBX < calibBX) meEventTypePreCalib.fill(dccId, runType, 1. / 54.);
       else if(dccBX == calibBX) meEventTypeCalib.fill(dccId, runType, 1. / 54.);
       else meEventTypePostCalib.fill(dccId, runType, 1. / 54.);
 
     }
   }

void ecaldqm::RawDataTask::runOnSource ( FEDRawDataCollection const & _fedRaw )

Definition at line 55 of file RawDataTask.cc.

References FEDRawDataCollection::FEDData(), l1A_, and ecaldqm::DQWorker::MEs_.

Referenced by analyze().

   {
     MESet& meCRC(MEs_.at("CRC"));
 
     // Get GT L1 info
     const FEDRawData &gtFED(_fedRaw.FEDData(812));
     if(gtFED.size() > sizeof(uint64_t)){ // FED header is one 64 bit word
       uint32_t *halfHeader((uint32_t *)gtFED.data());
       l1A_ = *(halfHeader + 1) & 0xffffff;
     }
 
     for(int iFED(601); iFED <= 654; iFED++){
       const FEDRawData& fedData(_fedRaw.FEDData(iFED));
       unsigned length(fedData.size() / sizeof(uint64_t));
       if(length > 1){ // FED header is one 64 bit word
         const uint64_t* pData(reinterpret_cast<uint64_t const*>(fedData.data()));
         if((pData[length - 1] & 0x4) != 0) meCRC.fill(iFED - 600);
       }
     }
   }