#include <ESIntegrityTask.h>

Inheritance diagram for ESIntegrityTask:

Public Member Functions
	ESIntegrityTask (const edm::ParameterSet &ps)

	~ESIntegrityTask () override

Public Member Functions inherited from DQMOneEDAnalyzer< edm::LuminosityBlockCache< ESIntLSCache > >
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

	DQMOneEDAnalyzer ()

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

void	calculateDIFraction (const edm::LuminosityBlock &lumi, const edm::EventSetup &c)
	Calculate Data Integrity Fraction. More...

void	dqmEndRun (const edm::Run &r, const edm::EventSetup &c) override
	EndRun. More...

void	endJob (void) override
	EndJob. More...

std::shared_ptr< ESIntLSCache >	globalBeginLuminosityBlock (const edm::LuminosityBlock &lumi, const edm::EventSetup &c) const override
	Begin Lumi. More...

void	globalEndLuminosityBlock (const edm::LuminosityBlock &lumi, const edm::EventSetup &c) override
	End Lumi. More...

Protected Member Functions inherited from DQMOneEDAnalyzer< edm::LuminosityBlockCache< ESIntLSCache > >
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Private Attributes
int	dac_

edm::EDGetTokenT< ESRawDataCollection >	dccCollections_

bool	doLumiAnalysis_

int	eCount_

int	fed_ [2][2][40][40]

int	fiber_ [2][2][40][40]

int	firstDAC_

int	gain_

int	ievt_

int	isPed_

int	kchip_ [2][2][40][40]

edm::EDGetTokenT< ESLocalRawDataCollection >	kchipCollections_

edm::FileInPath	lookup_

MonitorElement *	meDCCCRCErr_

MonitorElement *	meDCCCRCErrByLS_

MonitorElement *	meDCCErr_

MonitorElement *	meDIErrors_ [2][2]

MonitorElement *	meDIErrorsByLS_ [2][2]

MonitorElement *	meDIFraction_

MonitorElement *	meEVDR_

MonitorElement *	meFED_

MonitorElement *	meFiberBadStatus_

MonitorElement *	meFiberErrCode_

MonitorElement *	meFiberErrCodeByLS_

MonitorElement *	meFiberOff_

MonitorElement *	meFiberOffByLS_

MonitorElement *	meGain_

MonitorElement *	meKBC_

MonitorElement *	meKEC_

MonitorElement *	meKF1_

MonitorElement *	meKF2_

MonitorElement *	meOptoBC_

MonitorElement *	meOptoBCByLS_

MonitorElement *	meOptoRX_

MonitorElement *	meSLinkCRCErr_

MonitorElement *	meSLinkCRCErrByLS_

int	nDAC_

int	precision_

std::string	prefixME_

int	runNum_

int	runtype_

int	seqtype_

int	vDAC_ [5]

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer< edm::LuminosityBlockCache< ESIntLSCache > >
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Attributes inherited from DQMOneEDAnalyzer< edm::LuminosityBlockCache< ESIntLSCache > >
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 17 of file ESIntegrityTask.h.

Constructor & Destructor Documentation

◆ ESIntegrityTask()

ESIntegrityTask::ESIntegrityTask ( const edm::ParameterSet & ps )

Definition at line 25 of file ESIntegrityTask.cc.

References gather_cfg::cout, FCDTask_cfi::fiber, geometryDiff::file, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), dqmiolumiharvest::j, dqmdumpme::k, and visualization-live-secondInstance_cfg::m.

                                                        {
   prefixME_ = ps.getUntrackedParameter<string>("prefixME", "");
   lookup_ = ps.getUntrackedParameter<FileInPath>("LookupTable");
 
   dccCollections_ = consumes<ESRawDataCollection>(ps.getParameter<InputTag>("ESDCCCollections"));
   kchipCollections_ = consumes<ESLocalRawDataCollection>(ps.getParameter<InputTag>("ESKChipCollections"));
 
   // read in look-up table
   for (int i = 0; i < 2; ++i)
     for (int j = 0; j < 2; ++j)
       for (int k = 0; k < 40; ++k)
         for (int m = 0; m < 40; ++m) {
           fed_[i][j][k][m] = -1;
           kchip_[i][j][k][m] = -1;
           fiber_[i][j][k][m] = -1;
         }
 
   int nLines_, z, iz, ip, ix, iy, fed, kchip, pace, bundle, fiber, optorx;
   ifstream file;
   file.open(lookup_.fullPath().c_str());
   if (file.is_open()) {
     file >> nLines_;
 
     for (int i = 0; i < nLines_; ++i) {
       file >> iz >> ip >> ix >> iy >> fed >> kchip >> pace >> bundle >> fiber >> optorx;
 
       z = (iz == -1) ? 2 : iz;
       fed_[z - 1][ip - 1][ix - 1][iy - 1] = fed;
       kchip_[z - 1][ip - 1][ix - 1][iy - 1] = kchip;
       fiber_[z - 1][ip - 1][ix - 1][iy - 1] = (fiber - 1) + (optorx - 1) * 12;
     }
   } else {
     cout << "ESIntegrityTask : Look up table file can not be found in " << lookup_.fullPath().c_str() << endl;
   }
 
   ievt_ = 0;
 }

◆ ~ESIntegrityTask()

ESIntegrityTask::~ESIntegrityTask ( )

inlineoverride

Definition at line 20 of file ESIntegrityTask.h.

20 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Analyze.

Reimplemented from DQMOneEDAnalyzer< edm::LuminosityBlockCache< ESIntLSCache > >.

Definition at line 231 of file ESIntegrityTask.cc.

References edm::SortedCollection< T, SORT >::begin(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::dcc(), MillePedeFileConverter_cfg::e, edm::SortedCollection< T, SORT >::end(), sistrip::extrainfo::gain_, ESKCHIPBlock::getBC(), ESKCHIPBlock::getEC(), ESKCHIPBlock::getFlag1(), ESKCHIPBlock::getFlag2(), ESKCHIPBlock::getOptoBC(), ESKCHIPBlock::getOptoEC(), mps_fire::i, ESKCHIPBlock::id(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), and dqmiolumiharvest::j.

                                                                  {
   ievt_++;
   auto lumiCache = luminosityBlockCache(e.getLuminosityBlock().index());
   ++(lumiCache->ievtLS_);
 
   Handle<ESRawDataCollection> dccs;
   Handle<ESLocalRawDataCollection> kchips;
 
   // Fill # of events
   meDCCErr_->Fill(575, 2, 1);
   meDCCCRCErr_->Fill(575, 2, 1);
   meDCCCRCErrByLS_->Fill(575, 2, 1);
   meOptoRX_->Fill(575, 2, 1);
   meOptoBC_->Fill(575, 2, 1);
   meOptoBCByLS_->Fill(575, 2, 1);
   meFiberBadStatus_->Fill(575, 36, 1);
   meFiberOff_->Fill(575, 36, 1);
   meFiberOffByLS_->Fill(575, 36, 1);
   meEVDR_->Fill(575, 36, 1);
 
   // # of DI errors
   Double_t nDIErr[56][36];
   for (int i = 0; i < 56; ++i)
     for (int j = 0; j < 36; ++j)
       nDIErr[i][j] = 0;
 
   // DCC
   vector<int> fiberStatus;
   if (e.getByToken(dccCollections_, dccs)) {
     for (ESRawDataCollection::const_iterator dccItr = dccs->begin(); dccItr != dccs->end(); ++dccItr) {
       const ESDCCHeaderBlock& dcc = (*dccItr);
 
       meFED_->Fill(dcc.fedId());
 
       meDCCErr_->Fill(dcc.fedId(), dcc.getDCCErrors());
 
       // SLink CRC error
       if (dcc.getDCCErrors() == 101) {
         meSLinkCRCErr_->Fill(dcc.fedId());
         meSLinkCRCErrByLS_->Fill(dcc.fedId());
         for (int j = 0; j < 36; ++j)
           nDIErr[dcc.fedId() - 520][j]++;
       }
 
       if (dcc.getOptoRX0() == 129) {
         meOptoRX_->Fill(dcc.fedId(), 0);
         if (((dcc.getOptoBC0() - 15) & 0x0fff) != dcc.getBX()) {
           meOptoBC_->Fill(dcc.fedId(), 0);
           meOptoBCByLS_->Fill(dcc.fedId(), 0);
         }
       }
       if (dcc.getOptoRX1() == 129) {
         meOptoRX_->Fill(dcc.fedId(), 1);
         if (((dcc.getOptoBC1() - 15) & 0x0fff) != dcc.getBX()) {
           meOptoBC_->Fill(dcc.fedId(), 1);
           meOptoBCByLS_->Fill(dcc.fedId(), 1);
         }
       }
       if (dcc.getOptoRX2() == 129) {
         meOptoRX_->Fill(dcc.fedId(), 2);
         if (((dcc.getOptoBC2() - 15) & 0x0fff) != dcc.getBX()) {
           meOptoBC_->Fill(dcc.fedId(), 2);
           meOptoBCByLS_->Fill(dcc.fedId(), 2);
         }
       }
 
       if (dcc.getOptoRX0() == 128) {
         meDCCCRCErr_->Fill(dcc.fedId(), 0);
         meDCCCRCErrByLS_->Fill(dcc.fedId(), 0);
         for (int j = 0; j < 12; ++j)
           nDIErr[dcc.fedId() - 520][j]++;
       }
       if (dcc.getOptoRX1() == 128) {
         meDCCCRCErr_->Fill(dcc.fedId(), 1);
         meDCCCRCErrByLS_->Fill(dcc.fedId(), 1);
         for (int j = 12; j < 24; ++j)
           nDIErr[dcc.fedId() - 520][j]++;
       }
       if (dcc.getOptoRX2() == 128) {
         meDCCCRCErr_->Fill(dcc.fedId(), 2);
         meDCCCRCErrByLS_->Fill(dcc.fedId(), 2);
         for (int j = 24; j < 36; ++j)
           nDIErr[dcc.fedId() - 520][j]++;
       }
 
       fiberStatus = dcc.getFEChannelStatus();
 
       for (unsigned int i = 0; i < fiberStatus.size(); ++i) {
         if (fiberStatus[i] == 4 || fiberStatus[i] == 8 || fiberStatus[i] == 10 || fiberStatus[i] == 11 ||
             fiberStatus[i] == 12) {
           meFiberBadStatus_->Fill(dcc.fedId(), i + 1, 1);
           meFiberErrCode_->Fill(fiberStatus[i]);
           meFiberErrCodeByLS_->Fill(fiberStatus[i]);
           nDIErr[dcc.fedId() - 520][i]++;
         }
         if (fiberStatus[i] == 7) {
           meFiberOff_->Fill(dcc.fedId(), i + 1, 1);
           meFiberOffByLS_->Fill(dcc.fedId(), i + 1, 1);
         }
         if (fiberStatus[i] == 6) {
           meFiberErrCode_->Fill(fiberStatus[i]);
           meFiberErrCodeByLS_->Fill(fiberStatus[i]);
           meEVDR_->Fill(dcc.fedId(), i + 1, 1);
         }
       }
 
       runtype_ = dcc.getRunType();
       seqtype_ = dcc.getSeqType();
       dac_ = dcc.getDAC();
       gain_ = dcc.getGain();
       precision_ = dcc.getPrecision();
 
       meGain_->Fill(gain_);
     }
   } else {
     LogWarning("ESIntegrityTask") << "dccCollections not available";
   }
 
   // KCHIP's
   if (e.getByToken(kchipCollections_, kchips)) {
     for (ESLocalRawDataCollection::const_iterator kItr = kchips->begin(); kItr != kchips->end(); ++kItr) {
       ESKCHIPBlock kchip = (*kItr);
 
       meKF1_->Fill(kchip.id(), kchip.getFlag1());
       meKF2_->Fill(kchip.id(), kchip.getFlag2());
       if (kchip.getBC() != kchip.getOptoBC())
         meKBC_->Fill(kchip.id());
       if (kchip.getEC() != kchip.getOptoEC())
         meKEC_->Fill(kchip.id());
     }
   } else {
     LogWarning("ESIntegrityTask") << "kchipCollections not available";
   }
 
   // Fill # of DI errors
   for (int iz = 0; iz < 2; ++iz)
     for (int ip = 0; ip < 2; ++ip)
       for (int ix = 0; ix < 40; ++ix)
         for (int iy = 0; iy < 40; ++iy) {
           if (fed_[iz][ip][ix][iy] == -1)
             continue;
 
           if (nDIErr[fed_[iz][ip][ix][iy] - 520][fiber_[iz][ip][ix][iy]] > 0) {
             meDIErrors_[iz][ip]->Fill(ix + 1, iy + 1, 1);
             (lumiCache->DIErrorsByLS_)[iz][ip][ix][iy] += 1;
             meDIErrorsByLS_[iz][ip]->Fill(ix + 1, iy + 1, 1);
           }
         }
 }

◆ bookHistograms()

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

overrideprotectedvirtual

Implements DQMOneEDAnalyzer< edm::LuminosityBlockCache< ESIntLSCache > >.

Definition at line 99 of file ESIntegrityTask.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), timingPdfMaker::histo, mps_fire::i, dqmiolumiharvest::j, dqm::impl::MonitorElement::setAxisTitle(), and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                                   {
   char histo[200];
 
   iBooker.setCurrentFolder(prefixME_ + "/ESIntegrityTask");
 
   sprintf(histo, "ES FEDs used for data taking");
   meFED_ = iBooker.book1D(histo, histo, 56, 519.5, 575.5);
   meFED_->setAxisTitle("ES FED", 1);
   meFED_->setAxisTitle("Num of Events", 2);
 
   sprintf(histo, "ES Gain used for data taking");
   meGain_ = iBooker.book1D(histo, histo, 2, -0.5, 1.5);
   meGain_->setAxisTitle("ES Gain", 1);
   meGain_->setAxisTitle("Num of Events", 2);
 
   sprintf(histo, "ES DCC Error codes");
   meDCCErr_ = iBooker.book2D(histo, histo, 56, 519.5, 575.5, 6, -0.5, 5.5);
   meDCCErr_->setAxisTitle("ES FED", 1);
   meDCCErr_->setAxisTitle("ES DCC Error code", 2);
 
   sprintf(histo, "ES SLink CRC Errors");
   meSLinkCRCErr_ = iBooker.book1D(histo, histo, 56, 519.5, 575.5);
   meSLinkCRCErr_->setAxisTitle("ES FED", 1);
   meSLinkCRCErr_->setAxisTitle("Num of Events", 2);
 
   sprintf(histo, "ES DCC CRC Errors");
   meDCCCRCErr_ = iBooker.book2D(histo, histo, 56, 519.5, 575.5, 3, -0.5, 2.5);
   meDCCCRCErr_->setAxisTitle("ES FED", 1);
   meDCCCRCErr_->setAxisTitle("ES OptoRX", 2);
 
   sprintf(histo, "ES OptoRX used for data taking");
   meOptoRX_ = iBooker.book2D(histo, histo, 56, 519.5, 575.5, 3, -0.5, 2.5);
   meOptoRX_->setAxisTitle("ES FED", 1);
   meOptoRX_->setAxisTitle("ES OptoRX", 2);
 
   sprintf(histo, "ES OptoRX BC mismatch");
   meOptoBC_ = iBooker.book2D(histo, histo, 56, 519.5, 575.5, 3, -0.5, 2.5);
   meOptoBC_->setAxisTitle("ES FED", 1);
   meOptoBC_->setAxisTitle("ES OptoRX", 2);
 
   sprintf(histo, "ES Fiber Bad Status");
   meFiberBadStatus_ = iBooker.book2D(histo, histo, 56, 519.5, 575.5, 36, 0.5, 36.5);
   meFiberBadStatus_->setAxisTitle("ES FED", 1);
   meFiberBadStatus_->setAxisTitle("Fiber Number", 2);
 
   sprintf(histo, "ES Fiber Error Code");
   meFiberErrCode_ = iBooker.book1D(histo, histo, 17, -0.5, 16.5);
   meFiberErrCode_->setAxisTitle("Fiber Error Code", 1);
 
   sprintf(histo, "ES Fiber Off");
 
   meFiberOff_ = iBooker.book2D(histo, histo, 56, 519.5, 575.5, 36, 0.5, 36.5);
   meFiberOff_->setAxisTitle("ES FED", 1);
   meFiberOff_->setAxisTitle("Fiber Number", 2);
 
   sprintf(histo, "ES Warning Event Dropped");
   meEVDR_ = iBooker.book2D(histo, histo, 56, 519.5, 575.5, 36, 0.5, 36.5);
   meEVDR_->setAxisTitle("ES FED", 1);
   meEVDR_->setAxisTitle("Fiber Number", 2);
 
   sprintf(histo, "ES KChip Flag 1 Error codes");
   meKF1_ = iBooker.book2D(histo, histo, 1550, -0.5, 1549.5, 16, -0.5, 15.5);
   meKF1_->setAxisTitle("ES KChip", 1);
   meKF1_->setAxisTitle("ES KChip F1 Error Code ", 2);
 
   sprintf(histo, "ES KChip Flag 2 Error codes");
   meKF2_ = iBooker.book2D(histo, histo, 1550, -0.5, 1549.5, 16, -0.5, 15.5);
   meKF2_->setAxisTitle("ES KChip", 1);
   meKF2_->setAxisTitle("ES KChip F1 Error Code ", 2);
 
   sprintf(histo, "ES KChip BC mismatch with OptoRX");
   meKBC_ = iBooker.book1D(histo, histo, 1550, -0.5, 1549.5);
   meKBC_->setAxisTitle("ES KChip", 1);
   meKBC_->setAxisTitle("Num of BC mismatch", 2);
 
   sprintf(histo, "ES KChip EC mismatch with OptoRX");
   meKEC_ = iBooker.book1D(histo, histo, 1550, -0.5, 1549.5);
   meKEC_->setAxisTitle("ES KChip", 1);
   meKEC_->setAxisTitle("Num of EC mismatch", 2);
 
   for (int i = 0; i < 2; ++i)
     for (int j = 0; j < 2; ++j) {
       int iz = (i == 0) ? 1 : -1;
       sprintf(histo, "ES Integrity Errors Z %d P %d", iz, j + 1);
       meDIErrors_[i][j] = iBooker.book2D(histo, histo, 40, 0.5, 40.5, 40, 0.5, 40.5);
       meDIErrors_[i][j]->setAxisTitle("Si X", 1);
       meDIErrors_[i][j]->setAxisTitle("Si Y", 2);
     }
 
   // LS-based histos
   sprintf(histo, "ES Good Channel Fraction");
   meDIFraction_ = iBooker.book2D(histo, histo, 3, 1.0, 3.0, 3, 1.0, 3.0);
 
   iBooker.setCurrentFolder(prefixME_ + "/ByLumiSection");
 
   sprintf(histo, "ES SLink CRC Errors");
   meSLinkCRCErrByLS_ = iBooker.book1D(histo, histo, 56, 519.5, 575.5);
   meSLinkCRCErrByLS_->setAxisTitle("ES FED", 1);
   meSLinkCRCErrByLS_->setAxisTitle("Num of Events", 2);
 
   sprintf(histo, "ES DCC CRC Errors");
   meDCCCRCErrByLS_ = iBooker.book2D(histo, histo, 56, 519.5, 575.5, 3, -0.5, 2.5);
   meDCCCRCErrByLS_->setAxisTitle("ES FED", 1);
   meDCCCRCErrByLS_->setAxisTitle("ES OptoRX", 2);
 
   sprintf(histo, "ES OptoRX BC mismatch");
   meOptoBCByLS_ = iBooker.book2D(histo, histo, 56, 519.5, 575.5, 3, -0.5, 2.5);
   meOptoBCByLS_->setAxisTitle("ES FED", 1);
   meOptoBCByLS_->setAxisTitle("ES OptoRX", 2);
 
   sprintf(histo, "ES Fiber Error Code");
   meFiberErrCodeByLS_ = iBooker.book1D(histo, histo, 17, -0.5, 16.5);
   meFiberErrCodeByLS_->setAxisTitle("Fiber Error Code", 1);
 
   sprintf(histo, "ES Fiber Off");
   meFiberOffByLS_ = iBooker.book2D(histo, histo, 56, 519.5, 575.5, 36, 0.5, 36.5);
   meFiberOffByLS_->setAxisTitle("ES FED", 1);
   meFiberOffByLS_->setAxisTitle("Fiber Number", 2);
 
   for (int i = 0; i < 2; ++i) {
     for (int j = 0; j < 2; ++j) {
       int iz = (i == 0) ? 1 : -1;
       sprintf(histo, "ES Integrity Errors Z %d P %d", iz, j + 1);
       meDIErrorsByLS_[i][j] = iBooker.book2D(histo, histo, 40, 0.5, 40.5, 40, 0.5, 40.5);
       meDIErrorsByLS_[i][j]->setAxisTitle("Si X", 1);
       meDIErrorsByLS_[i][j]->setAxisTitle("Si Y", 2);
     }
   }
 }

◆ calculateDIFraction()

void ESIntegrityTask::calculateDIFraction	(	const edm::LuminosityBlock &	lumi,
		const edm::EventSetup &	c
	)

protected

Calculate Data Integrity Fraction.

Definition at line 383 of file ESIntegrityTask.cc.

References mps_fire::i, dqmiolumiharvest::j, heppy_batch::val, and x.

                                                                                                 {
   float nValidChannels = 0;
   float nGlobalErrors = 0;
   auto lumiCache = luminosityBlockCache(lumi.index());
 
   for (int i = 0; i < 2; ++i) {
     for (int j = 0; j < 2; ++j) {
       float nValidChannelsES = 0;
       float nGlobalErrorsES = 0;
       float reportSummaryES = -1;
       for (int x = 0; x < 40; ++x) {
         for (int y = 0; y < 40; ++y) {
           float val = 1.0 * ((lumiCache->DIErrorsByLS_)[i][j][x][y]);
           if (fed_[i][j][x][y] == -1)
             continue;
           if ((lumiCache->ievtLS_) != 0)
             nGlobalErrors += val / (lumiCache->ievtLS_);
           nValidChannels++;
           if ((lumiCache->ievtLS_) != 0)
             nGlobalErrorsES += val / (lumiCache->ievtLS_);
           nValidChannelsES++;
         }
       }
       if (nValidChannelsES != 0)
         reportSummaryES = 1 - nGlobalErrorsES / nValidChannelsES;
       meDIFraction_->setBinContent(i + 1, j + 1, reportSummaryES);
     }
   }
   float reportSummary = -1.0;
   if (nValidChannels != 0)
     reportSummary = 1.0 - nGlobalErrors / nValidChannels;
   meDIFraction_->setBinContent(3, 3, reportSummary);
 }

◆ dqmEndRun()

void ESIntegrityTask::dqmEndRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

EndRun.

Reimplemented from DQMOneEDAnalyzer< edm::LuminosityBlockCache< ESIntLSCache > >.

Definition at line 63 of file ESIntegrityTask.cc.

                                                                  {
   // In case of Lumi based analysis Disable SoftReset from Integrity histogram to get full statistics
   // TODO: no longer possible, clone histo beforehand if full statisticcs at end of run are required.
 }

◆ endJob()

void ESIntegrityTask::endJob ( void )

overrideprotectedvirtual

EndJob.

Reimplemented from edm::one::EDProducerBase.

Definition at line 229 of file ESIntegrityTask.cc.

229 { LogInfo("ESIntegrityTask") << "analyzed " << ievt_ << " events"; }

ESIntegrityTask::ievt_

int ievt_

Definition: ESIntegrityTask.h:45

edm::LogInfo

Log< level::Info, false > LogInfo

Definition: MessageLogger.h:131

◆ globalBeginLuminosityBlock()

std::shared_ptr< ESIntLSCache > ESIntegrityTask::globalBeginLuminosityBlock	(	const edm::LuminosityBlock &	lumi,
		const edm::EventSetup &	c
	)		const

overrideprotected

Begin Lumi.

Definition at line 68 of file ESIntegrityTask.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), and ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy().

                                                                                                         {
   LogInfo("ESIntegrityTask") << "analyzed " << ievt_ << " events";
   // In case of Lumi based analysis SoftReset the Integrity histogram
   auto lumiCache = std::make_shared<ESIntLSCache>();
   lumiCache->ievtLS_ = 0;
 
   for (int iz = 0; iz < 2; ++iz) {
     for (int ip = 0; ip < 2; ++ip) {
       meDIErrorsByLS_[iz][ip]->Reset();
       for (int ix = 0; ix < 40; ++ix) {
         for (int iy = 0; iy < 40; ++iy) {
           (lumiCache->DIErrorsByLS_)[iz][ip][ix][iy] = 0;
         }
       }
     }
   }
 
   meDCCCRCErrByLS_->Reset();
   meFiberErrCodeByLS_->Reset();
   meFiberOffByLS_->Reset();
   meOptoBCByLS_->Reset();
   meSLinkCRCErrByLS_->Reset();
 
   return lumiCache;
 }

◆ globalEndLuminosityBlock()

void ESIntegrityTask::globalEndLuminosityBlock	(	const edm::LuminosityBlock &	lumi,
		const edm::EventSetup &	c
	)

overrideprotected

End Lumi.

Definition at line 95 of file ESIntegrityTask.cc.

References DummyCfis::c.

                                                                                                      {
   calculateDIFraction(lumi, c);
 }

Member Data Documentation

◆ dac_

int ESIntegrityTask::dac_

private

Definition at line 78 of file ESIntegrityTask.h.

◆ dccCollections_

edm::EDGetTokenT<ESRawDataCollection> ESIntegrityTask::dccCollections_

private

Definition at line 49 of file ESIntegrityTask.h.

◆ doLumiAnalysis_

bool ESIntegrityTask::doLumiAnalysis_

private

Definition at line 81 of file ESIntegrityTask.h.

◆ eCount_

int ESIntegrityTask::eCount_

private

Definition at line 78 of file ESIntegrityTask.h.

◆ fed_

int ESIntegrityTask::fed_[2][2][40][40]

private

Definition at line 80 of file ESIntegrityTask.h.

◆ fiber_

int ESIntegrityTask::fiber_[2][2][40][40]

private

Definition at line 80 of file ESIntegrityTask.h.

◆ firstDAC_

int ESIntegrityTask::firstDAC_

private

Definition at line 79 of file ESIntegrityTask.h.

◆ gain_

int ESIntegrityTask::gain_

private

Definition at line 78 of file ESIntegrityTask.h.

◆ ievt_

int ESIntegrityTask::ievt_

private

Definition at line 45 of file ESIntegrityTask.h.

◆ isPed_

int ESIntegrityTask::isPed_

private

Definition at line 79 of file ESIntegrityTask.h.

◆ kchip_

int ESIntegrityTask::kchip_[2][2][40][40]

private

Definition at line 80 of file ESIntegrityTask.h.

◆ kchipCollections_

edm::EDGetTokenT<ESLocalRawDataCollection> ESIntegrityTask::kchipCollections_

private

Definition at line 50 of file ESIntegrityTask.h.

◆ lookup_

edm::FileInPath ESIntegrityTask::lookup_

private

Definition at line 76 of file ESIntegrityTask.h.

◆ meDCCCRCErr_

MonitorElement* ESIntegrityTask::meDCCCRCErr_

private

Definition at line 57 of file ESIntegrityTask.h.

◆ meDCCCRCErrByLS_

MonitorElement* ESIntegrityTask::meDCCCRCErrByLS_

private

Definition at line 58 of file ESIntegrityTask.h.

◆ meDCCErr_

MonitorElement* ESIntegrityTask::meDCCErr_

private

Definition at line 56 of file ESIntegrityTask.h.

◆ meDIErrors_

MonitorElement* ESIntegrityTask::meDIErrors_[2][2]

private

Definition at line 72 of file ESIntegrityTask.h.

◆ meDIErrorsByLS_

MonitorElement* ESIntegrityTask::meDIErrorsByLS_[2][2]

private

Definition at line 73 of file ESIntegrityTask.h.

◆ meDIFraction_

MonitorElement* ESIntegrityTask::meDIFraction_

private

Definition at line 74 of file ESIntegrityTask.h.

◆ meEVDR_

MonitorElement* ESIntegrityTask::meEVDR_

private

Definition at line 67 of file ESIntegrityTask.h.

◆ meFED_

MonitorElement* ESIntegrityTask::meFED_

private

Definition at line 53 of file ESIntegrityTask.h.

◆ meFiberBadStatus_

MonitorElement* ESIntegrityTask::meFiberBadStatus_

private

Definition at line 62 of file ESIntegrityTask.h.

◆ meFiberErrCode_

MonitorElement* ESIntegrityTask::meFiberErrCode_

private

Definition at line 63 of file ESIntegrityTask.h.

◆ meFiberErrCodeByLS_

MonitorElement* ESIntegrityTask::meFiberErrCodeByLS_

private

Definition at line 64 of file ESIntegrityTask.h.

◆ meFiberOff_

MonitorElement* ESIntegrityTask::meFiberOff_

private

Definition at line 65 of file ESIntegrityTask.h.

◆ meFiberOffByLS_

MonitorElement* ESIntegrityTask::meFiberOffByLS_

private

Definition at line 66 of file ESIntegrityTask.h.

◆ meGain_

MonitorElement* ESIntegrityTask::meGain_

private

Definition at line 52 of file ESIntegrityTask.h.

◆ meKBC_

MonitorElement* ESIntegrityTask::meKBC_

private

Definition at line 70 of file ESIntegrityTask.h.

◆ meKEC_

MonitorElement* ESIntegrityTask::meKEC_

private

Definition at line 71 of file ESIntegrityTask.h.

◆ meKF1_

MonitorElement* ESIntegrityTask::meKF1_

private

Definition at line 68 of file ESIntegrityTask.h.

◆ meKF2_

MonitorElement* ESIntegrityTask::meKF2_

private

Definition at line 69 of file ESIntegrityTask.h.

◆ meOptoBC_

MonitorElement* ESIntegrityTask::meOptoBC_

private

Definition at line 60 of file ESIntegrityTask.h.

◆ meOptoBCByLS_

MonitorElement* ESIntegrityTask::meOptoBCByLS_

private

Definition at line 61 of file ESIntegrityTask.h.

◆ meOptoRX_

MonitorElement* ESIntegrityTask::meOptoRX_

private

Definition at line 59 of file ESIntegrityTask.h.

◆ meSLinkCRCErr_

MonitorElement* ESIntegrityTask::meSLinkCRCErr_

private

Definition at line 54 of file ESIntegrityTask.h.

◆ meSLinkCRCErrByLS_

MonitorElement* ESIntegrityTask::meSLinkCRCErrByLS_

private

Definition at line 55 of file ESIntegrityTask.h.

◆ nDAC_

int ESIntegrityTask::nDAC_

private

Definition at line 79 of file ESIntegrityTask.h.

◆ precision_

int ESIntegrityTask::precision_

private

Definition at line 78 of file ESIntegrityTask.h.

◆ prefixME_

std::string ESIntegrityTask::prefixME_

private

Definition at line 47 of file ESIntegrityTask.h.

◆ runNum_

int ESIntegrityTask::runNum_

private

Definition at line 78 of file ESIntegrityTask.h.

◆ runtype_

int ESIntegrityTask::runtype_

private

Definition at line 78 of file ESIntegrityTask.h.

◆ seqtype_

int ESIntegrityTask::seqtype_

private

Definition at line 78 of file ESIntegrityTask.h.

◆ vDAC_

int ESIntegrityTask::vDAC_[5]

private

Definition at line 79 of file ESIntegrityTask.h.

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ESIntegrityTask()

◆ ~ESIntegrityTask()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ calculateDIFraction()

◆ dqmEndRun()

◆ endJob()

◆ globalBeginLuminosityBlock()

◆ globalEndLuminosityBlock()

Member Data Documentation

◆ dac_

◆ dccCollections_

◆ doLumiAnalysis_

◆ eCount_

◆ fed_

◆ fiber_

◆ firstDAC_

◆ gain_

◆ ievt_

◆ isPed_

◆ kchip_

◆ kchipCollections_

◆ lookup_

◆ meDCCCRCErr_

◆ meDCCCRCErrByLS_

◆ meDCCErr_

◆ meDIErrors_

◆ meDIErrorsByLS_

◆ meDIFraction_

◆ meEVDR_

◆ meFED_

◆ meFiberBadStatus_

◆ meFiberErrCode_

◆ meFiberErrCodeByLS_

◆ meFiberOff_

◆ meFiberOffByLS_

◆ meGain_

◆ meKBC_

◆ meKEC_

◆ meKF1_

◆ meKF2_

◆ meOptoBC_

◆ meOptoBCByLS_

◆ meOptoRX_

◆ meSLinkCRCErr_

◆ meSLinkCRCErrByLS_

◆ nDAC_

◆ precision_

◆ prefixME_

◆ runNum_

◆ runtype_

◆ seqtype_

◆ vDAC_