#include <EcalSelectiveReadoutValidation.h>

Inheritance diagram for EcalSelectiveReadoutValidation:

Classes
struct	energiesEb_t

struct	energiesEe_t

Public Member Functions
void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override

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

	EcalSelectiveReadoutValidation (const edm::ParameterSet &ps)
	Constructor. More...

	~EcalSelectiveReadoutValidation () override
	Destructor. More...

Public Member Functions inherited from DQMOneEDAnalyzer<>
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 final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~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 (edm::Event const &e, edm::EventSetup const &c) override
	Analyzes the event. More...

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

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 Types
typedef EcalRecHit	RecHit

typedef EcalRecHitCollection	RecHitCollection

enum	subdet_t { EB, EE }
	distinguishes barral and endcap of ECAL. More...

Private Member Functions
template<class T , class U >
void	anaDigi (const T &frame, const U &srFlagColl)

void	anaDigiInit ()

void	analyzeDataVolume (const edm::Event &e, const edm::EventSetup &es)

void	analyzeEB (const edm::Event &event, const edm::EventSetup &es)

void	analyzeEE (const edm::Event &event, const edm::EventSetup &es)

void	analyzeTP (const edm::Event &event, const edm::EventSetup &es)

template<class T >
void	checkSrApplication (const edm::Event &event, T &srfs)

int	cIndex2iEta (int i) const

int	cIndex2iPhi (int i) const

int	cIndex2iTtEta (int i) const

int	cIndex2iTtPhi (int i) const

int	cIndex2iXY (int iX0) const

template<class T >
void	compareSrfColl (const edm::Event &event, T &srfFromData, T &computedSrf)

void	configFirWeights (const std::vector< double > &weightsForZsFIR)

std::pair< int, int >	dccCh (const DetId &xtalId) const

double	frame2Energy (const EcalDataFrame &frame) const

template<class T >
double	frame2EnergyForTp (const T &frame, int offset=0) const

double	getBytesPerCrystal () const

int	getCrystalCount (int iDcc, int iDccCh)

double	getDccEventSize (int iDcc0, double nReadXtals) const

double	getDccOverhead (subdet_t subdet) const

double	getDccSrDependentPayload (int iDcc0, double nReadRus, double nReadXtals) const

double	getEbEventSize (double nReadXtals) const

double	getEeEventSize (double nReadXtals) const

double	getL1aRate () const

int	getRuCount (int iDcc0) const

int	iEta2cIndex (int iEta) const

void	initAsciiFile ()

int	iPhi2cIndex (int iPhi) const

int	iTtEta2cIndex (int iEta) const

int	iTtPhi2cIndex (int iPhi) const

int	iXY2cIndex (int iX) const

void	printAvailableHists ()

void	readAllCollections (const edm::Event &e)

bool	registerHist (const std::string &name, const std::string &title)

void	selectFedsForLog ()

void	setTtEtSums (const edm::EventSetup &es, const EBDigiCollection &ebDigis, const EEDigiCollection &eeDigis)

void	updateL1aRate (const edm::Event &event)


EcalTrigTowerDetId	readOutUnitOf (const EBDetId &xtalId) const

EcalScDetId	readOutUnitOf (const EEDetId &xtalId) const


MonitorElement *	bookFloat (DQMStore::IBooker &, const std::string &name)

MonitorElement *	book1D (DQMStore::IBooker &, const std::string &name, const std::string &title, int nbins, double xmin, double xmax)

MonitorElement *	book2D (DQMStore::IBooker &, const std::string &name, const std::string &title, int nxbins, double xmin, double xmax, int nybins, double ymin, double ymax)

MonitorElement *	bookProfile (DQMStore::IBooker &, const std::string &name, const std::string &title, int nbins, double xmin, double xmax)

MonitorElement *	bookProfile2D (DQMStore::IBooker &, const std::string &name, const std::string &title, int nbinx, double xmin, double xmax, int nbiny, double ymin, double ymax, const char *option="")


void	fill (MonitorElement *me, float x)

void	fill (MonitorElement *me, float x, float yw)

void	fill (MonitorElement *me, float x, float y, float zw)

void	fill (MonitorElement *me, float x, float y, float z, float w)


int	ruGraphX (const EcalScDetId &id) const

int	ruGraphY (const EcalScDetId &id) const

int	ruGraphX (const EcalTrigTowerDetId &id) const

int	ruGraphY (const EcalTrigTowerDetId &id) const

int	xtalGraphX (const EEDetId &id) const

int	xtalGraphY (const EEDetId &id) const

int	xtalGraphX (const EBDetId &id) const

int	xtalGraphY (const EBDetId &id) const


int	dccId (const EcalScDetId &detId) const

int	dccId (const EcalTrigTowerDetId &detId) const

Static Private Member Functions
static int	dccZsFIR (const EcalDataFrame &frame, const std::vector< int > &firWeights, int firstFIRSample, bool *saturated=nullptr)

static std::vector< int >	getFIRWeights (const std::vector< double > &normalizedWeights)

Private Attributes
bool	allHists_

std::map< std::string, std::string >	availableHistList_

bool	collNotFoundWarn_
	Switch for collection-not-found warning. More...

int	ebZsThr_

edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	ecalmapping

int	eeZsThr_

const EcalElectronicsMapping *	elecMap_

int	firstFIRSample_

std::vector< int >	firWeights_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geoToken

std::string	histDir_

std::set< std::string >	histList_

edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord >	hTriggerTowerMap

int	ievt_

bool	isRuComplete_ [nDccs_][nDccChs_]

bool	localReco_

std::vector< bool >	logErrForDccs_

bool	logSrApplicationErrors_

bool	logSrpAlgoErrors_

edm::ESGetToken< EcalTPGLutGroup, EcalTPGLutGroupRcd >	lutGrpHandle

edm::ESGetToken< EcalTPGLutIdMap, EcalTPGLutIdMapRcd >	lutMapHandle

MonitorElement *	meL1aRate_

int	nCompleteZS_
	Counter of ZS-flagged RU fully read out. More...

int	nDroppedFRO_
	Counter of FRO-flagged RU dropped from data. More...

int	nEb_

int	nEbFROCnt_
	Counter of EB FRO-flagged RUs. More...

int	nEbHI_

int	nEbLI_

int	nEbZsErrors_
	Counter of EB ZS errors (LI channel below ZS threshold) More...

int	nEbZsErrorsType1_

int	nEe_

int	nEeFROCnt_
	Counter of EE FRO-flagged RUs. More...

int	nEeHI_

int	nEeLI_

int	nEeZsErrors_
	Counter of EE ZS errors (LI channel below ZS threshold) More...

int	nEeZsErrorsType1_

int	nHiPerDcc_ [nDccs_]

int	nHiRuPerDcc_ [nDccs_]

int	nIncompleteFRO_
	Counter of FRO-flagged RU only partial data. More...

int	nLiPerDcc_ [nDccs_]

int	nLiRuPerDcc_ [nDccs_]

int	nPerDcc_ [nDccs_]

int	nPerRu_ [nDccs_][nDccChs_]

int	nRuPerDcc_ [nDccs_]

std::string	outputFile_
	Output file for histograms. More...

edm::ESGetToken< EcalTPGPhysicsConst, EcalTPGPhysicsConstRcd >	physHandle

bool	SkipInnerSC_

std::ofstream	srApplicationErrorLog_
	Output ascii file for unconsistency between Xtals and RU Flags. More...

std::string	srApplicationErrorLogFileName_

std::ofstream	srpAlgoErrorLog_
	Output ascii file for unconsistency on SR flags. More...

std::string	srpAlgoErrorLogFileName_

bool	tpInGeV_

const EcalTrigTowerConstituentsMap *	triggerTowerMap_

double	ttEtSums [nTtEta][nTtPhi]

bool	useEventRate_

bool	verbose_
	Verbosity switch. More...

std::vector< double >	weights_

bool	withEbSimHit_

bool	withEeSimHit_

std::ofstream	zsErrorLog_
	File to log ZS and other errors. More...


CollHandle< EBDigiCollection >	ebDigis_

CollHandle< EEDigiCollection >	eeDigis_

CollHandle< EBDigiCollection >	ebNoZsDigis_

CollHandle< EEDigiCollection >	eeNoZsDigis_

CollHandle< EBSrFlagCollection >	ebSrFlags_

CollHandle< EESrFlagCollection >	eeSrFlags_

CollHandle< EBSrFlagCollection >	ebComputedSrFlags_

CollHandle< EESrFlagCollection >	eeComputedSrFlags_

CollHandle< std::vector< PCaloHit > >	ebSimHits_

CollHandle< std::vector< PCaloHit > >	eeSimHits_

CollHandle< EcalTrigPrimDigiCollection >	tps_

CollHandle< RecHitCollection >	ebRecHits_

CollHandle< RecHitCollection >	eeRecHits_

CollHandle< FEDRawDataCollection >	fedRaw_


int64_t	tmax

int64_t	tmin

int64_t	l1aOfTmin

int64_t	l1aOfTmax

bool	l1aRateErr


MonitorElement *	meDccVol_

MonitorElement *	meDccLiVol_

MonitorElement *	meDccHiVol_

MonitorElement *	meDccVolFromData_

MonitorElement *	meVol_

MonitorElement *	meVolB_

MonitorElement *	meVolE_

MonitorElement *	meVolBLI_

MonitorElement *	meVolELI_

MonitorElement *	meVolLI_

MonitorElement *	meVolBHI_

MonitorElement *	meVolEHI_

MonitorElement *	meVolHI_

MonitorElement *	meChOcc_

MonitorElement *	meTp_

MonitorElement *	meTtf_

MonitorElement *	meTtfVsTp_

MonitorElement *	meTtfVsEtSum_

MonitorElement *	meTpVsEtSum_

MonitorElement *	meEbRecE_

MonitorElement *	meEbEMean_

MonitorElement *	meEbNoise_

MonitorElement *	meEbSimE_

MonitorElement *	meEbRecEHitXtal_

MonitorElement *	meEbRecVsSimE_

MonitorElement *	meEbNoZsRecVsSimE_

MonitorElement *	meEeRecE_

MonitorElement *	meEeEMean_

MonitorElement *	meEeNoise_

MonitorElement *	meEeSimE_

MonitorElement *	meEeRecEHitXtal_

MonitorElement *	meEeRecVsSimE_

MonitorElement *	meEeNoZsRecVsSimE_

MonitorElement *	meFullRoRu_

MonitorElement *	meZs1Ru_

MonitorElement *	meForcedRu_

MonitorElement *	meLiTtf_

MonitorElement *	meMiTtf_

MonitorElement *	meHiTtf_

MonitorElement *	meForcedTtf_

MonitorElement *	meTpMap_

MonitorElement *	meFullRoCnt_

MonitorElement *	meEbFullRoCnt_

MonitorElement *	meEeFullRoCnt_

MonitorElement *	meEbLiZsFir_

MonitorElement *	meEbHiZsFir_

MonitorElement *	meEbIncompleteRUZsFir_

MonitorElement *	meEeLiZsFir_

MonitorElement *	meEeHiZsFir_

MonitorElement *	meSRFlagsFromData_

MonitorElement *	meSRFlagsComputed_

MonitorElement *	meSRFlagsConsistency_

MonitorElement *	meIncompleteFRO_

MonitorElement *	meDroppedFRO_

MonitorElement *	meCompleteZS_

MonitorElement *	meIncompleteFROMap_

MonitorElement *	meDroppedFROMap_

MonitorElement *	meCompleteZSMap_

MonitorElement *	meIncompleteFRORateMap_

MonitorElement *	meDroppedFRORateMap_

MonitorElement *	meCompleteZSRateMap_

MonitorElement *	meIncompleteFROCnt_

MonitorElement *	meDroppedFROCnt_

MonitorElement *	meCompleteZSCnt_

MonitorElement *	meEbZsErrCnt_

MonitorElement *	meEeZsErrCnt_

MonitorElement *	meZsErrCnt_

MonitorElement *	meEbZsErrType1Cnt_

MonitorElement *	meEeZsErrType1Cnt_

MonitorElement *	meZsErrType1Cnt_


MonitorElement *	meEbFixedPayload_

MonitorElement *	meEeFixedPayload_

MonitorElement *	meFixedPayload_


bool	ebRuActive_ [nEbEta/ebTtEdge][nEbPhi/ebTtEdge]

bool	eeRuActive_ [nEndcaps][nEeX/scEdge][nEeY/scEdge]

Static Private Attributes
static const int	ebTtEdge = 5
	Number of crystals along an EB TT. More...

static const int	kByte_ = 1024
	number of bytes in 1 kByte: More...

static const int	maxDccId_ = minDccId_ + nDccs_ - 1

static const int	minDccId_ = 1

static const unsigned	nDccChs_ = 68
	Number of input channels of a DCC. More...

static const int	nDccRus_ [nDccs_]
	number of RUs for each DCC More...

static const unsigned	nDccs_ = 54
	Total number of DCCs. More...

static const int	nEbDccs = 36
	number of DCCs for EB More...

static const int	nEbEta = 170
	number of crystals along Eta in EB More...

static const int	nEbPhi = 360
	number of crystals along Phi in EB More...

static const int	nEbRus = 36 * 68
	number of RUs for EB More...

static const int	nEbTtEta = 34
	Number of Trigger Towers in barrel along Eta. More...

static const int	nEeDccs = 18
	number of DCCs for EE More...

static const int	nEeRus = 2 * (34 + 32 + 33 + 33 + 32 + 34 + 33 + 34 + 33)
	number of RUs for EE More...

static const int	nEeX = 100
	EE crystal grid size along X. More...

static const int	nEeY = 100
	EE crystal grid size along Y. More...

static const int	nEndcaps = 2
	number of endcaps More...

static const int	nMaxXtalPerRu = 25
	Number of crystals per Readout Unit excepted partial SCs. More...

static const int	nOneEeTtEta = 11
	Number of Trigger Towers in an endcap along Eta. More...

static const int	nTtEta = 2 * nOneEeTtEta + nEbTtEta
	Number of Trigger Towers along Eta. More...

static const int	nTtPhi = 72
	Number of Trigger Towers along Phi. More...

static const double	rad2deg = 45. / atan(1.)
	Conversion factor from radian to degree. More...

static const int	scEdge = 5
	Number of crystals along a supercrystal edge. More...

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer<>
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::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 46 of file EcalSelectiveReadoutValidation.h.

Member Typedef Documentation

◆ RecHit

typedef EcalRecHit EcalSelectiveReadoutValidation::RecHit

private

Definition at line 48 of file EcalSelectiveReadoutValidation.h.

◆ RecHitCollection

typedef EcalRecHitCollection EcalSelectiveReadoutValidation::RecHitCollection

private

Definition at line 47 of file EcalSelectiveReadoutValidation.h.

Member Enumeration Documentation

◆ subdet_t

enum EcalSelectiveReadoutValidation::subdet_t

private

distinguishes barral and endcap of ECAL.

Enumerator
EB
EE

Definition at line 67 of file EcalSelectiveReadoutValidation.h.

67 { EB, EE };

EcalSelectiveReadoutValidation::EE

Definition: EcalSelectiveReadoutValidation.h:67

EcalSelectiveReadoutValidation::EB

Definition: EcalSelectiveReadoutValidation.h:67

Constructor & Destructor Documentation

◆ EcalSelectiveReadoutValidation()

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

Constructor.

Definition at line 98 of file EcalSelectiveReadoutValidation.cc.

References allHists_, edm::BeginRun, configFirWeights(), edm::EDConsumerBase::consumesCollector(), ebComputedSrFlags_, ebDigis_, ebNoZsDigis_, ebRecHits_, ebSimHits_, ebSrFlags_, ebZsThr_, eeComputedSrFlags_, eeDigis_, eeNoZsDigis_, eeRecHits_, eeSimHits_, eeSrFlags_, eeZsThr_, fedRaw_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), histList_, compare::hists, logSrApplicationErrors_, logSrpAlgoErrors_, outputFile_, CollHandle< T >::setToken(), srApplicationErrorLogFileName_, srpAlgoErrorLogFileName_, tps_, and verbose_.

     : geoToken(esConsumes()),
       ecalmapping(esConsumes<edm::Transition::BeginRun>()),
       hTriggerTowerMap(esConsumes<edm::Transition::BeginRun>()),
       physHandle(esConsumes()),
       lutGrpHandle(esConsumes()),
       lutMapHandle(esConsumes()),
       collNotFoundWarn_(ps.getUntrackedParameter<bool>("warnIfCollectionNotFound", true)),
       ebDigis_(ps.getParameter<edm::InputTag>("EbDigiCollection"), false, collNotFoundWarn_),
       eeDigis_(ps.getParameter<edm::InputTag>("EeDigiCollection"), false, collNotFoundWarn_),
       ebNoZsDigis_(ps.getParameter<edm::InputTag>("EbUnsuppressedDigiCollection"), false, false /*collNotFoundWarn_*/),
       eeNoZsDigis_(ps.getParameter<edm::InputTag>("EeUnsuppressedDigiCollection"), false, false /*collNotFoundWarn_*/),
       ebSrFlags_(ps.getParameter<edm::InputTag>("EbSrFlagCollection"), false, collNotFoundWarn_),
       eeSrFlags_(ps.getParameter<edm::InputTag>("EeSrFlagCollection"), false, collNotFoundWarn_),
       ebComputedSrFlags_(
           ps.getParameter<edm::InputTag>("EbSrFlagFromTTCollection"), false, false /*collNotFoundWarn_*/),
       eeComputedSrFlags_(
           ps.getParameter<edm::InputTag>("EeSrFlagFromTTCollection"), false, false /*collNotFoundWarn_*/),
       ebSimHits_(ps.getParameter<edm::InputTag>("EbSimHitCollection"), false, false /*collNotFoundWarn_*/),
       eeSimHits_(ps.getParameter<edm::InputTag>("EeSimHitCollection"), false, false /*collNotFoundWarn_*/),
       tps_(ps.getParameter<edm::InputTag>("TrigPrimCollection"), false, collNotFoundWarn_),
       ebRecHits_(ps.getParameter<edm::InputTag>("EbRecHitCollection"), false, false /*collNotFoundWarn_*/),
       eeRecHits_(ps.getParameter<edm::InputTag>("EeRecHitCollection"), false, false /*collNotFoundWarn_*/),
       fedRaw_(ps.getParameter<edm::InputTag>("FEDRawCollection"), false, false /*collNotFoundWarn_*/),
       tmax(0),
       tmin(numeric_limits<int64_t>::max()),
       l1aOfTmin(0),
       l1aOfTmax(0),
       triggerTowerMap_(nullptr),
       localReco_(ps.getParameter<bool>("LocalReco")),
       weights_(ps.getParameter<vector<double>>("weights")),
       tpInGeV_(ps.getParameter<bool>("tpInGeV")),
       firstFIRSample_(ps.getParameter<int>("ecalDccZs1stSample")),
       useEventRate_(ps.getParameter<bool>("useEventRate")),
       logErrForDccs_(nDccs_, false),
       ievt_(0),
       allHists_(false),
       histDir_(ps.getParameter<string>("histDir")),
       withEeSimHit_(false),
       withEbSimHit_(false) {
   edm::ConsumesCollector collector(consumesCollector());
   ebDigis_.setToken(collector);
   eeDigis_.setToken(collector);
   ebNoZsDigis_.setToken(collector);
   eeNoZsDigis_.setToken(collector);
   ebSrFlags_.setToken(collector);
   eeSrFlags_.setToken(collector);
   ebComputedSrFlags_.setToken(collector);
   eeComputedSrFlags_.setToken(collector);
   ebSimHits_.setToken(collector);
   eeSimHits_.setToken(collector);
   tps_.setToken(collector);
   ebRecHits_.setToken(collector);
   eeRecHits_.setToken(collector);
   fedRaw_.setToken(collector);
 
   double ebZsThr = ps.getParameter<double>("ebZsThrADCCount");
   double eeZsThr = ps.getParameter<double>("eeZsThrADCCount");
 
   ebZsThr_ = lround(ebZsThr * 4);
   eeZsThr_ = lround(eeZsThr * 4);
 
   //File to log SRP algorithem inconsistency
   srpAlgoErrorLogFileName_ = ps.getUntrackedParameter<string>("srpAlgoErrorLogFile", "");
   logSrpAlgoErrors_ = (!srpAlgoErrorLogFileName_.empty());
 
   //File to log SRP decision application inconsistency
   srApplicationErrorLogFileName_ = ps.getUntrackedParameter<string>("srApplicationErrorLogFile", "");
   logSrApplicationErrors_ = (!srApplicationErrorLogFileName_.empty());
 
   //FIR ZS weights
   configFirWeights(ps.getParameter<vector<double>>("dccWeights"));
 
   // DQM ROOT output
   outputFile_ = ps.getUntrackedParameter<string>("outputFile", "");
 
   if (!outputFile_.empty()) {
     LogInfo("OutputInfo") << " Ecal Digi Task histograms will be saved to '" << outputFile_.c_str() << "'";
   } else {
     LogInfo("OutputInfo") << " Ecal Digi Task histograms will NOT be saved";
   }
 
   // verbosity switch
   verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
 
   // get hold of back-end interface
 
   vector<string> hists(ps.getUntrackedParameter<vector<string>>("histograms", vector<string>(1, "all")));
 
   for (vector<string>::iterator it = hists.begin(); it != hists.end(); ++it)
     histList_.insert(*it);
   if (histList_.find("all") != histList_.end())
     allHists_ = true;
 }

◆ ~EcalSelectiveReadoutValidation()

EcalSelectiveReadoutValidation::~EcalSelectiveReadoutValidation ( )

override

Destructor.

Definition at line 840 of file EcalSelectiveReadoutValidation.cc.

840 {}

Member Function Documentation

◆ anaDigi()

template<class T , class U >

void EcalSelectiveReadoutValidation::anaDigi	(	const T &	frame,
		const U &	srFlagColl
	)

private

Accumulates statitics for data volume analysis. To be called for each ECAL digi. See anaDigiInit().

Definition at line 1658 of file EcalSelectiveReadoutValidation.cc.

References Reference_intrackfit_cff::barrel, dccCh(), ebRuActive_, ebTtEdge, EcalBarrel, eeRuActive_, Exception, RemoveAddSevLevel::flag, hcalRecHitTable_cff::ieta, iEta2cIndex(), hcalRecHitTable_cff::iphi, iPhi2cIndex(), iXY2cIndex(), minDccId_, nEb_, nEbHI_, nEbLI_, nEe_, nEeHI_, nEeLI_, nHiPerDcc_, nHiRuPerDcc_, nLiPerDcc_, nLiRuPerDcc_, nPerDcc_, nPerRu_, nRuPerDcc_, DetId::rawId(), readOutUnitOf(), scEdge, EcalSrFlag::SRF_FULL, DetId::subdetId(), and ecaldqm::zside().

Referenced by analyzeDataVolume().

                                                                                 {
   const DetId& xtalId = frame.id();
   typedef typename U::key_type RuDetId;
   const RuDetId& ruId = readOutUnitOf(frame.id());
   typename U::const_iterator srf = srFlagColl.find(ruId);
 
   bool highInterest = false;
   int flag = 0;
 
   if (srf != srFlagColl.end()) {
     flag = srf->value() & ~EcalSrFlag::SRF_FORCED_MASK;
 
     highInterest = (flag == EcalSrFlag::SRF_FULL);
   }
 
   bool barrel = (xtalId.subdetId() == EcalBarrel);
 
   pair<int, int> ch = dccCh(xtalId);
 
   if (barrel) {
     ++nEb_;
     if (highInterest) {
       ++nEbHI_;
     } else {  //low interest
       ++nEbLI_;
     }
     int iEta0 = iEta2cIndex(static_cast<const EBDetId&>(xtalId).ieta());
     int iPhi0 = iPhi2cIndex(static_cast<const EBDetId&>(xtalId).iphi());
     if (!ebRuActive_[iEta0 / ebTtEdge][iPhi0 / ebTtEdge]) {
       ++nRuPerDcc_[ch.first - minDccId_];
       if (highInterest) {
         ++nHiRuPerDcc_[ch.first - minDccId_];
       } else {
         ++nLiRuPerDcc_[ch.first - minDccId_];
       }
 
       ebRuActive_[iEta0 / ebTtEdge][iPhi0 / ebTtEdge] = true;
     }
   } else {  //endcap
     ++nEe_;
     if (highInterest) {
       ++nEeHI_;
     } else {  //low interest
       ++nEeLI_;
     }
     int iX0 = iXY2cIndex(static_cast<const EEDetId&>(frame.id()).ix());
     int iY0 = iXY2cIndex(static_cast<const EEDetId&>(frame.id()).iy());
     int iZ0 = static_cast<const EEDetId&>(frame.id()).zside() > 0 ? 1 : 0;
 
     if (!eeRuActive_[iZ0][iX0 / scEdge][iY0 / scEdge]) {
       ++nRuPerDcc_[ch.first - minDccId_];
       if (highInterest) {
         ++nHiRuPerDcc_[ch.first - minDccId_];
       } else {
         ++nLiRuPerDcc_[ch.first - minDccId_];
       }
 
       eeRuActive_[iZ0][iX0 / scEdge][iY0 / scEdge] = true;
     }
   }
 
   if (ch.second < 1 || ch.second > 68) {
     throw cms::Exception("EcalSelectiveReadoutValidation")
         << "Error in DCC channel retrieval for crystal with detId " << xtalId.rawId()
         << "DCC channel out of allowed range [1..68]\n";
   }
   ++nPerDcc_[ch.first - minDccId_];
   ++nPerRu_[ch.first - minDccId_][ch.second - 1];
   if (highInterest) {
     ++nHiPerDcc_[ch.first - minDccId_];
   } else {  //low interest channel
     ++nLiPerDcc_[ch.first - minDccId_];
   }
 }

◆ anaDigiInit()

void EcalSelectiveReadoutValidation::anaDigiInit ( )

private

Initializes statistics accumalator for data volume analysis. To be call at start of each event analysis.

Definition at line 1733 of file EcalSelectiveReadoutValidation.cc.

References ebRuActive_, eeRuActive_, nEb_, nEbHI_, nEbLI_, nEe_, nEeHI_, nEeLI_, nHiPerDcc_, nHiRuPerDcc_, nLiPerDcc_, nLiRuPerDcc_, nPerDcc_, nPerRu_, and nRuPerDcc_.

Referenced by analyzeDataVolume().

                                                  {
   nEb_ = 0;
   nEe_ = 0;
   nEeLI_ = 0;
   nEeHI_ = 0;
   nEbLI_ = 0;
   nEbHI_ = 0;
   bzero(nPerDcc_, sizeof(nPerDcc_));
   bzero(nLiPerDcc_, sizeof(nLiPerDcc_));
   bzero(nHiPerDcc_, sizeof(nHiPerDcc_));
   bzero(nRuPerDcc_, sizeof(nRuPerDcc_));
   bzero(ebRuActive_, sizeof(ebRuActive_));
   bzero(eeRuActive_, sizeof(eeRuActive_));
   bzero(nPerRu_, sizeof(nPerRu_));
   bzero(nLiRuPerDcc_, sizeof(nLiRuPerDcc_));
   bzero(nHiRuPerDcc_, sizeof(nHiRuPerDcc_));
 }

◆ analyze()

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

overrideprotectedvirtual

Analyzes the event.

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 218 of file EcalSelectiveReadoutValidation.cc.

References analyzeDataVolume(), analyzeEB(), analyzeEE(), analyzeTP(), checkSrApplication(), compareSrfColl(), ebComputedSrFlags_, ebNoZsDigis_, ebSimHits_, ebSrFlags_, eeComputedSrFlags_, eeNoZsDigis_, eeSimHits_, eeSrFlags_, fill(), ievt_, meCompleteZSCnt_, meDroppedFROCnt_, meEbFullRoCnt_, meEbZsErrCnt_, meEbZsErrType1Cnt_, meEeFullRoCnt_, meEeZsErrCnt_, meEeZsErrType1Cnt_, meFullRoCnt_, meIncompleteFROCnt_, meZsErrCnt_, meZsErrType1Cnt_, nCompleteZS_, nDroppedFRO_, nEbFROCnt_, nEbZsErrors_, nEbZsErrorsType1_, nEeFROCnt_, nEeZsErrors_, nEeZsErrorsType1_, nIncompleteFRO_, readAllCollections(), selectFedsForLog(), setTtEtSums(), tps_, updateL1aRate(), withEbSimHit_, and withEeSimHit_.

                                                                                      {
   updateL1aRate(event);
 
   //retrieves event products:
   readAllCollections(event);
 
   withEeSimHit_ = (!eeSimHits_->empty());
   withEbSimHit_ = (!ebSimHits_->empty());
 
   if (ievt_ < 10) {
     edm::LogInfo("EcalSrValid") << "Size of TP collection: " << tps_->size() << std::endl
                                 << "Size of EB SRF collection read from data: " << ebSrFlags_->size() << std::endl
                                 << "Size of EB SRF collection computed from data TTFs: " << ebComputedSrFlags_->size()
                                 << std::endl
                                 << "Size of EE SRF collection read from data: " << eeSrFlags_->size() << std::endl
                                 << "Size of EE SRF collection computed from data TTFs: " << eeComputedSrFlags_->size()
                                 << std::endl;
   }
 
   if (ievt_ == 0) {
     selectFedsForLog();  //note: must be called after readAllCollection
   }
 
   //computes Et sum trigger tower crystals:
   setTtEtSums(es, *ebNoZsDigis_, *eeNoZsDigis_);
 
   //Data Volume
   analyzeDataVolume(event, es);
 
   //EB digis
   //must be called after analyzeDataVolume because it uses
   //isRuComplete_ array that this method fills
   analyzeEB(event, es);
 
   //EE digis
   //must be called after analyzeDataVolume because it uses
   //isRuComplete_ array that this method fills
   analyzeEE(event, es);
 
   fill(meFullRoCnt_, nEeFROCnt_ + nEbFROCnt_);
   fill(meEbFullRoCnt_, nEbFROCnt_);
   fill(meEeFullRoCnt_, nEeFROCnt_);
 
   fill(meEbZsErrCnt_, nEbZsErrors_);
   fill(meEeZsErrCnt_, nEeZsErrors_);
   fill(meZsErrCnt_, nEbZsErrors_ + nEeZsErrors_);
 
   fill(meEbZsErrType1Cnt_, nEbZsErrorsType1_);
   fill(meEeZsErrType1Cnt_, nEeZsErrorsType1_);
   fill(meZsErrType1Cnt_, nEbZsErrorsType1_ + nEeZsErrorsType1_);
 
   //TP
   analyzeTP(event, es);
 
   if (!ebComputedSrFlags_->empty()) {
     compareSrfColl(event, *ebSrFlags_, *ebComputedSrFlags_);
   }
   if (!eeComputedSrFlags_->empty()) {
     compareSrfColl(event, *eeSrFlags_, *eeComputedSrFlags_);
   }
   nDroppedFRO_ = 0;
   nIncompleteFRO_ = 0;
   nCompleteZS_ = 0;
   checkSrApplication(event, *ebSrFlags_);
   checkSrApplication(event, *eeSrFlags_);
   fill(meDroppedFROCnt_, nDroppedFRO_);
   fill(meIncompleteFROCnt_, nIncompleteFRO_);
   fill(meCompleteZSCnt_, nCompleteZS_);
   ++ievt_;
 }

◆ analyzeDataVolume()

void EcalSelectiveReadoutValidation::analyzeDataVolume	(	const edm::Event &	e,
		const edm::EventSetup &	es
	)

private

Data volume analysis. To be called for each event.

Parameters

event	EDM event
es	event setup

Definition at line 1604 of file EcalSelectiveReadoutValidation.cc.

References a, anaDigi(), anaDigiInit(), b, ebDigis_, ebSrFlags_, eeDigis_, eeSrFlags_, FEDRawDataCollection::FEDData(), fedRaw_, fill(), getBytesPerCrystal(), getCrystalCount(), getDccEventSize(), getDccSrDependentPayload(), getEbEventSize(), getEeEventSize(), isRuComplete_, kByte_, maxDccId_, meDccHiVol_, meDccLiVol_, meDccVol_, meDccVolFromData_, meVol_, meVolB_, meVolBHI_, meVolBLI_, meVolE_, meVolEHI_, meVolELI_, meVolHI_, meVolLI_, minDccId_, nDccRus_, nDccs_, nEb_, nEbHI_, nEbLI_, nEe_, nEeHI_, nEeLI_, nHiPerDcc_, nHiRuPerDcc_, nLiPerDcc_, nLiRuPerDcc_, nPerDcc_, nPerRu_, FEDRawData::size(), and edm::DataFrameContainer::size().

Referenced by analyze().

                                                                                            {
   anaDigiInit();
 
   //Complete RU, i.e. RU actually fully readout
   for (int iDcc = minDccId_; iDcc <= maxDccId_; ++iDcc) {
     for (int iCh = 1; iCh < nDccRus_[iDcc - minDccId_]; ++iCh) {
       isRuComplete_[iDcc - minDccId_][iCh - 1] = (nPerRu_[iDcc - minDccId_][iCh - 1] == getCrystalCount(iDcc, iCh));
     }
   }
 
   //Barrel
   for (unsigned int digis = 0; digis < ebDigis_->size(); ++digis) {
     EBDataFrame ebdf = (*ebDigis_)[digis];
     anaDigi(ebdf, *ebSrFlags_);
   }
 
   // Endcap
   for (unsigned int digis = 0; digis < eeDigis_->size(); ++digis) {
     EEDataFrame eedf = (*eeDigis_)[digis];
     anaDigi(eedf, *eeSrFlags_);
   }
 
   //histos
   for (unsigned iDcc0 = 0; iDcc0 < nDccs_; ++iDcc0) {
     fill(meDccVol_, iDcc0 + 1, getDccEventSize(iDcc0, nPerDcc_[iDcc0]) / kByte_);
     fill(meDccLiVol_, iDcc0 + 1, getDccSrDependentPayload(iDcc0, nLiRuPerDcc_[iDcc0], nLiPerDcc_[iDcc0]) / kByte_);
     fill(meDccHiVol_, iDcc0 + 1, getDccSrDependentPayload(iDcc0, nHiRuPerDcc_[iDcc0], nHiPerDcc_[iDcc0]) / kByte_);
     const FEDRawDataCollection& raw = *fedRaw_;
     fill(meDccVolFromData_, iDcc0 + 1, ((double)raw.FEDData(601 + iDcc0).size()) / kByte_);
   }
 
   //low interesest channels:
   double a = nEbLI_ * getBytesPerCrystal() / kByte_;  //getEbEventSize(nEbLI_)/kByte_;
   fill(meVolBLI_, a);
   double b = nEeLI_ * getBytesPerCrystal() / kByte_;  //getEeEventSize(nEeLI_)/kByte_;
   fill(meVolELI_, b);
   fill(meVolLI_, a + b);
 
   //high interest chanels:
   a = nEbHI_ * getBytesPerCrystal() / kByte_;  //getEbEventSize(nEbHI_)/kByte_;
   fill(meVolBHI_, a);
   b = nEeHI_ * getBytesPerCrystal() / kByte_;  //getEeEventSize(nEeHI_)/kByte_;
   fill(meVolEHI_, b);
   fill(meVolHI_, a + b);
 
   //any-interest channels:
   a = getEbEventSize(nEb_) / kByte_;
   fill(meVolB_, a);
   b = getEeEventSize(nEe_) / kByte_;
   fill(meVolE_, b);
   fill(meVol_, a + b);
 }

◆ analyzeEB()

void EcalSelectiveReadoutValidation::analyzeEB	(	const edm::Event &	event,
		const edm::EventSetup &	es
	)

private

ECAL barrel data analysis. To be called for each event.

Parameters

event	EDM event
es	event setup

Energy deposited in ECAL barrel crystals. Eta index starts from 0 at eta minimum and phi index starts at phi=0+ in CMS std coordinate system.

Definition at line 563 of file EcalSelectiveReadoutValidation.cc.

References cms::cuda::assert(), edm::DataFrameContainer::begin(), gather_cfg::cout, dccCh(), dccZsFIR(), hcalRecHitTable_cff::detId, ebDigis_, ebNoZsDigis_, ebRecHits_, ebSimHits_, ebSrFlags_, ebZsThr_, DetId::Ecal, EcalBarrel, edm::DataFrameContainer::end(), PV3DBase< T, PVType, FrameType >::eta(), Exception, fill(), dqmdumpme::first, firstFIRSample_, firWeights_, RemoveAddSevLevel::flag, frame2Energy(), ecalLiteDTU::gainId(), geoToken, CaloSubdetectorGeometry::getGeometry(), edm::EventSetup::getHandle(), mps_fire::i, EBSrFlag::id(), hit::id, hcalRecHitTable_cff::ieta, EcalTrigTowerDetId::ieta(), l1tTowerCalibrationProducer_cfi::iEta, iEta2cIndex(), EcalTrigTowerDetId::ietaAbs(), ievt_, hcalRecHitTable_cff::iphi, EcalTrigTowerDetId::iphi(), iPhi2cIndex(), isRuComplete_, localReco_, SiStripPI::max, meChOcc_, meEbEMean_, meEbHiZsFir_, meEbLiZsFir_, meEbNoise_, meEbNoZsRecVsSimE_, meEbRecE_, meEbRecEHitXtal_, meEbRecVsSimE_, meEbSimE_, meForcedRu_, meFullRoRu_, meZs1Ru_, nEbEta, nEbFROCnt_, nEbPhi, nEbZsErrors_, nEbZsErrorsType1_, EcalSelectiveReadoutValidation::energiesEb_t::noZsRecE, PV3DBase< T, PVType, FrameType >::phi(), rad2deg, readOutUnitOf(), EcalSelectiveReadoutValidation::energiesEb_t::recE, alignCSCRings::s, edm::second(), EcalSelectiveReadoutValidation::energiesEb_t::simE, rpcPointValidation_cfi::simHit, EcalSelectiveReadoutValidation::energiesEb_t::simHit, srApplicationErrorLog_, EcalSrFlag::SRF_FORCED_MASK, EcalSrFlag::SRF_FULL, EcalSrFlag::SRF_ZS1, findQualityFiles::v, EcalSrFlag::value(), withEbSimHit_, x, xtalGraphX(), xtalGraphY(), y, and EcalTrigTowerDetId::zside().

Referenced by analyze().

                                                                                              {
   bool eventError = false;
   nEbZsErrors_ = 0;
   nEbZsErrorsType1_ = 0;
   vector<pair<int, int>> xtalEtaPhi;
 
   xtalEtaPhi.reserve(nEbPhi * nEbEta);
 
   std::unique_ptr<std::array<std::array<energiesEb_t, nEbPhi>, nEbEta>> ebEnergies =
       std::make_unique<std::array<std::array<energiesEb_t, nEbPhi>, nEbEta>>();
 
   for (int iEta0 = 0; iEta0 < nEbEta; ++iEta0) {
     for (int iPhi0 = 0; iPhi0 < nEbPhi; ++iPhi0) {
       (*ebEnergies)[iEta0][iPhi0].noZsRecE = -numeric_limits<double>::max();
       (*ebEnergies)[iEta0][iPhi0].recE = -numeric_limits<double>::max();
       (*ebEnergies)[iEta0][iPhi0].simE = 0;  //must be zero.
       (*ebEnergies)[iEta0][iPhi0].simHit = 0;
       (*ebEnergies)[iEta0][iPhi0].gain12 = false;
       xtalEtaPhi.push_back(pair<int, int>(iEta0, iPhi0));
     }
   }
 
   // get the barrel geometry:
   auto geoHandle = es.getHandle(geoToken);
   const CaloSubdetectorGeometry* geometry_p = (*geoHandle).getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
   //CaloSubdetectorGeometry const& geometry = *geometry_p;
 
   //EB unsuppressed digis:
   for (EBDigiCollection::const_iterator it = ebNoZsDigis_->begin(); it != ebNoZsDigis_->end(); ++it) {
     const EBDataFrame& frame = *it;
     int iEta0 = iEta2cIndex(static_cast<const EBDetId&>(frame.id()).ieta());
     int iPhi0 = iPhi2cIndex(static_cast<const EBDetId&>(frame.id()).iphi());
     if (iEta0 < 0 || iEta0 >= nEbEta) {
       stringstream s;
       s << "EcalSelectiveReadoutValidation: "
         << "iEta0 (= " << iEta0 << ") is out of range ("
         << "[0," << nEbEta - 1 << "]\n";
       throw cms::Exception(s.str());
     }
     if (iPhi0 < 0 || iPhi0 >= nEbPhi) {
       stringstream s;
       s << "EcalSelectiveReadoutValidation: "
         << "iPhi0 (= " << iPhi0 << ") is out of range ("
         << "[0," << nEbPhi - 1 << "]\n";
       throw cms::Exception(s.str());
     }
 
     (*ebEnergies)[iEta0][iPhi0].noZsRecE = frame2Energy(frame);
     (*ebEnergies)[iEta0][iPhi0].gain12 = true;
     for (int i = 0; i < frame.size(); ++i) {
       const int gain12Code = 0x1;
       if (frame[i].gainId() != gain12Code)
         (*ebEnergies)[iEta0][iPhi0].gain12 = false;
     }
 
     const GlobalPoint xtalPos = geometry_p->getGeometry(frame.id())->getPosition();
 
     (*ebEnergies)[iEta0][iPhi0].phi = rad2deg * ((double)xtalPos.phi());
     (*ebEnergies)[iEta0][iPhi0].eta = xtalPos.eta();
   }  //next non-zs digi
 
   //EB sim hits
   for (vector<PCaloHit>::const_iterator it = ebSimHits_->begin(); it != ebSimHits_->end(); ++it) {
     const PCaloHit& simHit = *it;
     EBDetId detId(simHit.id());
     int iEta = detId.ieta();
     int iEta0 = iEta2cIndex(iEta);
     int iPhi = detId.iphi();
     int iPhi0 = iPhi2cIndex(iPhi);
     (*ebEnergies)[iEta0][iPhi0].simE += simHit.energy();
     ++(*ebEnergies)[iEta0][iPhi0].simHit;
   }
 
   bool crystalShot[nEbEta][nEbPhi];
   pair<int, int> EBxtalCoor[nEbEta][nEbPhi];
 
   for (int iEta0 = 0; iEta0 < nEbEta; ++iEta0) {
     for (int iPhi0 = 0; iPhi0 < nEbPhi; ++iPhi0) {
       crystalShot[iEta0][iPhi0] = false;
       EBxtalCoor[iEta0][iPhi0] = make_pair(0, 0);
     }
   }
 
   for (EBDigiCollection::const_iterator it = ebDigis_->begin(); it != ebDigis_->end(); ++it) {
     const EBDataFrame& frame = *it;
     int iEta = static_cast<const EBDetId&>(frame.id()).ieta();
     int iPhi = static_cast<const EBDetId&>(frame.id()).iphi();
     int iEta0 = iEta2cIndex(iEta);
     int iPhi0 = iPhi2cIndex(iPhi);
     if (iEta0 < 0 || iEta0 >= nEbEta) {
       throw(cms::Exception("EcalSelectiveReadoutValidation") << "iEta0 (= " << iEta0 << ") is out of range ("
                                                              << "[0," << nEbEta - 1 << "]");
     }
     if (iPhi0 < 0 || iPhi0 >= nEbPhi) {
       throw(cms::Exception("EcalSelectiveReadoutValidation") << "iPhi0 (= " << iPhi0 << ") is out of range ("
                                                              << "[0," << nEbPhi - 1 << "]");
     }
     assert(iEta0 >= 0 && iEta0 < nEbEta);
     assert(iPhi0 >= 0 && iPhi0 < nEbPhi);
     if (!crystalShot[iEta0][iPhi0]) {
       crystalShot[iEta0][iPhi0] = true;
       EBxtalCoor[iEta0][iPhi0] = make_pair(xtalGraphX(frame.id()), xtalGraphY(frame.id()));
     } else {
       cout << "Error: several digi for same crystal!";
       abort();
     }
     if (localReco_) {
       (*ebEnergies)[iEta0][iPhi0].recE = frame2Energy(frame);
     }
 
     (*ebEnergies)[iEta0][iPhi0].gain12 = true;
     for (int i = 0; i < frame.size(); ++i) {
       const int gain12Code = 0x1;
       if (frame[i].gainId() != gain12Code) {
         (*ebEnergies)[iEta0][iPhi0].gain12 = false;
       }
     }
 
     EBSrFlagCollection::const_iterator srf = ebSrFlags_->find(readOutUnitOf(frame.id()));
 
     bool highInterest = false;
 
     // if(srf == ebSrFlags_->end()){
     //       throw cms::Exception("EcalSelectiveReadoutValidation")
     //  << __FILE__ << ":" << __LINE__ << ": SR flag not found";
     //}
 
     if (srf != ebSrFlags_->end()) {
       highInterest = ((srf->value() & ~EcalSrFlag::SRF_FORCED_MASK) == EcalSrFlag::SRF_FULL);
     }
 
     if (highInterest) {
       fill(meEbHiZsFir_, dccZsFIR(frame, firWeights_, firstFIRSample_, nullptr));
     } else {
       int v = dccZsFIR(frame, firWeights_, firstFIRSample_, nullptr);
       fill(meEbLiZsFir_, v);
       if (v < ebZsThr_) {
         eventError = true;
         ++nEbZsErrors_;
         pair<int, int> ru = dccCh(frame.id());
         if (isRuComplete_[ru.first - 1][ru.second - 1])
           ++nEbZsErrorsType1_;
         if (nEbZsErrors_ < 3) {
           srApplicationErrorLog_ << event.id() << ", "
                                  << "RU " << frame.id() << ", "
                                  << "DCC " << ru.first << " Ch : " << ru.second << ": "
                                  << "LI channel under ZS threshold.\n";
         }
         if (nEbZsErrors_ == 3) {
           srApplicationErrorLog_ << event.id() << ": "
                                  << "more ZS errors for this event...\n";
         }
       }
     }
   }  //next EB digi
 
   for (int iEta0 = 0; iEta0 < nEbEta; ++iEta0) {
     for (int iPhi0 = 0; iPhi0 < nEbPhi; ++iPhi0)
       fill(meChOcc_,
            EBxtalCoor[iEta0][iPhi0].first,
            EBxtalCoor[iEta0][iPhi0].second,
            crystalShot[iEta0][iPhi0] ? 1. : 0.);
   }
 
   if (!localReco_) {
     for (RecHitCollection::const_iterator it = ebRecHits_->begin(); it != ebRecHits_->end(); ++it) {
       const RecHit& hit = *it;
       int iEta = static_cast<const EBDetId&>(hit.id()).ieta();
       int iPhi = static_cast<const EBDetId&>(hit.id()).iphi();
       int iEta0 = iEta2cIndex(iEta);
       int iPhi0 = iPhi2cIndex(iPhi);
       if (iEta0 < 0 || iEta0 >= nEbEta) {
         LogError("EcalSrValid") << "iEta0 (= " << iEta0 << ") is out of range ("
                                 << "[0," << nEbEta - 1 << "]\n";
       }
       if (iPhi0 < 0 || iPhi0 >= nEbPhi) {
         LogError("EcalSrValid") << "iPhi0 (= " << iPhi0 << ") is out of range ("
                                 << "[0," << nEbPhi - 1 << "]\n";
       }
       (*ebEnergies)[iEta0][iPhi0].recE = hit.energy();
     }
   }
 
   for (unsigned int i = 0; i < xtalEtaPhi.size(); ++i) {
     int iEta0 = xtalEtaPhi[i].first;
     int iPhi0 = xtalEtaPhi[i].second;
     energiesEb_t& energies = (*ebEnergies)[iEta0][iPhi0];
 
     double recE = energies.recE;
     if (recE != -numeric_limits<double>::max()) {  //not zero suppressed
       fill(meEbRecE_, (*ebEnergies)[iEta0][iPhi0].recE);
       fill(meEbEMean_, ievt_ + 1, recE);
     }  //not zero suppressed
 
     if (withEbSimHit_) {
       if (!energies.simHit) {  //noise only crystal channel
         fill(meEbNoise_, energies.noZsRecE);
       } else {
         fill(meEbSimE_, energies.simE);
         fill(meEbRecEHitXtal_, energies.recE);
       }
       fill(meEbRecVsSimE_, energies.simE, energies.recE);
       fill(meEbNoZsRecVsSimE_, energies.simE, energies.noZsRecE);
     }
   }
 
   int EBZs1RuCount[2][17][72];
   int EBFullRuCount[2][17][72];
   int EBForcedRuCount[2][17][72];
   std::pair<int, int> EBtowerCoor[2][17][72];
   for (int iZ(0); iZ < 2; iZ++) {
     for (int iEta(0); iEta < 17; iEta++) {
       for (int iPhi(0); iPhi < 72; iPhi++) {
         EBZs1RuCount[iZ][iEta][iPhi] = 0;
         EBFullRuCount[iZ][iEta][iPhi] = 0;
         EBForcedRuCount[iZ][iEta][iPhi] = 0;
       }
     }
   }
 
   //SRF
   nEbFROCnt_ = 0;
   char ebSrfMark[2][17][72];
   bzero(ebSrfMark, sizeof(ebSrfMark));
   //      int idbg = 0;
   for (EBSrFlagCollection::const_iterator it = ebSrFlags_->begin(); it != ebSrFlags_->end(); ++it) {
     const EBSrFlag& srf = *it;
     int iEtaAbs = srf.id().ietaAbs();
     int iPhi = srf.id().iphi();
     int iZ = srf.id().zside();
 
     //  cout << "--> " << ++idbg << iEtaAbs << " " << iPhi << " "  << iZ
     //       << " " << srf.id() << "\n";
 
     if (iEtaAbs < 1 || iEtaAbs > 17 || iPhi < 1 || iPhi > 72)
       throw cms::Exception("EcalSelectiveReadoutValidation")
           << "Found a barrel SRF with an invalid det ID: " << srf.id() << ".\n";
     ++ebSrfMark[iZ > 0 ? 1 : 0][iEtaAbs - 1][iPhi - 1];
     if (ebSrfMark[iZ > 0 ? 1 : 0][iEtaAbs - 1][iPhi - 1] > 1)
       throw cms::Exception("EcalSelectiveReadoutValidation") << "Duplicate SRF for RU " << srf.id() << ".\n";
 
     EBtowerCoor[iZ > 0 ? 1 : 0][iEtaAbs - 1][iPhi - 1] = std::pair<int, int>(srf.id().ieta(), srf.id().iphi());
 
     int flag = srf.value() & ~EcalSrFlag::SRF_FORCED_MASK;
     if (flag == EcalSrFlag::SRF_ZS1) {
       EBZs1RuCount[iZ > 0 ? 1 : 0][iEtaAbs - 1][iPhi - 1] += 1;
     }
     if (flag == EcalSrFlag::SRF_FULL) {
       EBFullRuCount[iZ > 0 ? 1 : 0][iEtaAbs - 1][iPhi - 1] += 1;
       ++nEbFROCnt_;
     }
     if (srf.value() & EcalSrFlag::SRF_FORCED_MASK) {
       EBForcedRuCount[iZ > 0 ? 1 : 0][iEtaAbs - 1][iPhi - 1] += 1;
     }
   }
   for (int iZ(0); iZ < 2; iZ++) {
     for (int iEta(0); iEta < 17; iEta++) {
       for (int iPhi(0); iPhi < 72; iPhi++) {
         float x(EBtowerCoor[iZ][iEta][iPhi].first);
         float y(EBtowerCoor[iZ][iEta][iPhi].second);
         fill(meZs1Ru_, x, y, EBZs1RuCount[iZ][iEta][iPhi]);
         fill(meFullRoRu_, x, y, EBFullRuCount[iZ][iEta][iPhi]);
         fill(meForcedRu_, x, y, EBForcedRuCount[iZ][iEta][iPhi]);
       }
     }
   }
 
   if (eventError)
     srApplicationErrorLog_ << event.id() << ": " << nEbZsErrors_
                            << " ZS-flagged EB channels under "
                               "the ZS threshold, whose "
                            << nEbZsErrorsType1_ << " in a complete RU.\n";
 }

◆ analyzeEE()

void EcalSelectiveReadoutValidation::analyzeEE	(	const edm::Event &	event,
		const edm::EventSetup &	es
	)

private

ECAL endcap data analysis. To be called for each event.

Parameters

event	EDM event
es	event setup

Energy deposited in ECAL endcap crystals. Endcap index is 0 for EE- and 1 for EE+. X and Y index starts at x and y minimum in std CMS coordinate system.

Definition at line 289 of file EcalSelectiveReadoutValidation.cc.

References edm::DataFrameContainer::begin(), gather_cfg::cout, dccCh(), dccZsFIR(), hcalRecHitTable_cff::detId, DetId::Ecal, EcalEndcap, eeDigis_, eeNoZsDigis_, eeRecHits_, eeSimHits_, eeSrFlags_, eeZsThr_, edm::DataFrameContainer::end(), PV3DBase< T, PVType, FrameType >::eta(), Exception, fill(), dqmdumpme::first, firstFIRSample_, firWeights_, RemoveAddSevLevel::flag, frame2Energy(), ecalLiteDTU::gainId(), geoToken, CaloSubdetectorGeometry::getGeometry(), edm::EventSetup::getHandle(), mps_fire::i, EESrFlag::id(), hit::id, ievt_, isRuComplete_, EcalScDetId::ix(), iXY2cIndex(), EcalScDetId::iy(), localReco_, SiStripPI::max, meChOcc_, meEeEMean_, meEeHiZsFir_, meEeLiZsFir_, meEeNoise_, meEeNoZsRecVsSimE_, meEeRecE_, meEeRecEHitXtal_, meEeRecVsSimE_, meEeSimE_, meForcedRu_, meFullRoRu_, meZs1Ru_, nEeFROCnt_, nEeX, nEeY, nEeZsErrors_, nEeZsErrorsType1_, nEndcaps, PV3DBase< T, PVType, FrameType >::phi(), rad2deg, readOutUnitOf(), edm::second(), rpcPointValidation_cfi::simHit, edm::DataFrameContainer::size(), srApplicationErrorLog_, EcalSrFlag::SRF_FORCED_MASK, EcalSrFlag::SRF_FULL, EcalSrFlag::SRF_ZS1, findQualityFiles::v, EcalSrFlag::value(), withEeSimHit_, xtalGraphX(), xtalGraphY(), EcalScDetId::zside(), and ecaldqm::zside().

Referenced by analyze().

                                                                                              {
   bool eventError = false;
   nEeZsErrors_ = 0;
   nEeZsErrorsType1_ = 0;
 
   std::unique_ptr<std::array<std::array<std::array<energiesEe_t, nEeY>, nEeX>, nEndcaps>> eeEnergies =
       std::make_unique<std::array<std::array<std::array<energiesEe_t, nEeY>, nEeX>, nEndcaps>>();
 
   for (int iZ0 = 0; iZ0 < nEndcaps; ++iZ0) {
     for (int iX0 = 0; iX0 < nEeX; ++iX0) {
       for (int iY0 = 0; iY0 < nEeY; ++iY0) {
         (*eeEnergies)[iZ0][iX0][iY0].noZsRecE = -numeric_limits<double>::max();
         (*eeEnergies)[iZ0][iX0][iY0].recE = -numeric_limits<double>::max();
         (*eeEnergies)[iZ0][iX0][iY0].simE = 0;  //must be set to zero.
         (*eeEnergies)[iZ0][iX0][iY0].simHit = 0;
         (*eeEnergies)[iZ0][iX0][iY0].gain12 = false;
       }
     }
   }
 
   // gets the endcap geometry:
   auto geoHandle = es.getHandle(geoToken);
   const CaloSubdetectorGeometry* geometry_p = (*geoHandle).getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
   //CaloSubdetectorGeometry const& geometry = *geometry_p;
 
   //EE unsupressed digis:
   for (unsigned int digis = 0; digis < eeNoZsDigis_->size(); ++digis) {
     EEDataFrame frame = (*eeNoZsDigis_)[digis];
     int iX0 = iXY2cIndex(frame.id().ix());
     int iY0 = iXY2cIndex(frame.id().iy());
     int iZ0 = frame.id().zside() > 0 ? 1 : 0;
 
     if (iX0 < 0 || iX0 >= nEeX) {
       edm::LogError("EcalSrValid") << "iX0 (= " << iX0 << ") is out of range ("
                                    << "[0," << nEeX - 1 << "]\n";
     }
     if (iY0 < 0 || iY0 >= nEeY) {
       edm::LogError("EcalSrValid") << "iY0 (= " << iY0 << ") is out of range ("
                                    << "[0," << nEeY - 1 << "]\n";
     }
     //    cout << "EE no ZS energy computation..." ;
     (*eeEnergies)[iZ0][iX0][iY0].noZsRecE = frame2Energy(frame);
 
     (*eeEnergies)[iZ0][iX0][iY0].gain12 = true;
     for (int i = 0; i < frame.size(); ++i) {
       const int gain12Code = 0x1;
       if (frame[i].gainId() != gain12Code)
         (*eeEnergies)[iZ0][iX0][iY0].gain12 = false;
     }
 
     const GlobalPoint xtalPos = geometry_p->getGeometry(frame.id())->getPosition();
 
     (*eeEnergies)[iZ0][iX0][iY0].phi = rad2deg * ((double)xtalPos.phi());
     (*eeEnergies)[iZ0][iX0][iY0].eta = xtalPos.eta();
   }
 
   //EE rec hits:
   if (!localReco_) {
     for (RecHitCollection::const_iterator it = eeRecHits_->begin(); it != eeRecHits_->end(); ++it) {
       const RecHit& hit = *it;
       int iX0 = iXY2cIndex(static_cast<const EEDetId&>(hit.id()).ix());
       int iY0 = iXY2cIndex(static_cast<const EEDetId&>(hit.id()).iy());
       int iZ0 = static_cast<const EEDetId&>(hit.id()).zside() > 0 ? 1 : 0;
 
       if (iX0 < 0 || iX0 >= nEeX) {
         LogError("EcalSrValid") << "iX0 (= " << iX0 << ") is out of range ("
                                 << "[0," << nEeX - 1 << "]\n";
       }
       if (iY0 < 0 || iY0 >= nEeY) {
         LogError("EcalSrValid") << "iY0 (= " << iY0 << ") is out of range ("
                                 << "[0," << nEeY - 1 << "]\n";
       }
       //    cout << "EE no ZS energy computation..." ;
       (*eeEnergies)[iZ0][iX0][iY0].recE = hit.energy();
     }
   }
 
   //EE sim hits:
   for (vector<PCaloHit>::const_iterator it = eeSimHits_->begin(); it != eeSimHits_->end(); ++it) {
     const PCaloHit& simHit = *it;
     EEDetId detId(simHit.id());
     int iX = detId.ix();
     int iX0 = iXY2cIndex(iX);
     int iY = detId.iy();
     int iY0 = iXY2cIndex(iY);
     int iZ0 = detId.zside() > 0 ? 1 : 0;
     (*eeEnergies)[iZ0][iX0][iY0].simE += simHit.energy();
     ++(*eeEnergies)[iZ0][iX0][iY0].simHit;
   }
 
   bool EEcrystalShot[nEeX][nEeY][2];
   pair<int, int> EExtalCoor[nEeX][nEeY][2];
 
   for (int iEeZ = 0; iEeZ < 2; ++iEeZ) {
     for (int iEeX = 0; iEeX < nEeX; ++iEeX) {
       for (int iEeY = 0; iEeY < nEeY; ++iEeY) {
         EEcrystalShot[iEeX][iEeY][iEeZ] = false;
         EExtalCoor[iEeX][iEeY][iEeZ] = make_pair(0, 0);
       }
     }
   }
 
   //EE suppressed digis
   for (EEDigiCollection::const_iterator it = eeDigis_->begin(); it != eeDigis_->end(); ++it) {
     const EEDataFrame& frame = *it;
     int iX0 = iXY2cIndex(static_cast<const EEDetId&>(frame.id()).ix());
     int iY0 = iXY2cIndex(static_cast<const EEDetId&>(frame.id()).iy());
     int iZ0 = static_cast<const EEDetId&>(frame.id()).zside() > 0 ? 1 : 0;
     if (iX0 < 0 || iX0 >= nEeX) {
       LogError("EcalSrValid") << "iX0 (= " << iX0 << ") is out of range ("
                               << "[0," << nEeX - 1 << "]\n";
     }
     if (iY0 < 0 || iY0 >= nEeY) {
       LogError("EcalSrValid") << "iY0 (= " << iY0 << ") is out of range ("
                               << "[0," << nEeY - 1 << "]\n";
     }
 
     if (!EEcrystalShot[iX0][iY0][iZ0]) {
       EEcrystalShot[iX0][iY0][iZ0] = true;
       EExtalCoor[iX0][iY0][iZ0] = make_pair(xtalGraphX(frame.id()), xtalGraphY(frame.id()));
     } else {
       cout << "Error: several digi for same crystal!";
       abort();
     }
 
     if (localReco_) {
       (*eeEnergies)[iZ0][iX0][iY0].recE = frame2Energy(frame);
     }
 
     (*eeEnergies)[iZ0][iX0][iY0].gain12 = true;
     for (int i = 0; i < frame.size(); ++i) {
       const int gain12Code = 0x1;
       if (frame[i].gainId() != gain12Code) {
         (*eeEnergies)[iZ0][iX0][iY0].gain12 = false;
       }
     }
 
     EESrFlagCollection::const_iterator srf = eeSrFlags_->find(readOutUnitOf(frame.id()));
 
     bool highInterest = false;
 
     if (srf == eeSrFlags_->end())
       continue;
 
     if (srf != eeSrFlags_->end()) {
       highInterest = ((srf->value() & ~EcalSrFlag::SRF_FORCED_MASK) == EcalSrFlag::SRF_FULL);
     }
 
     if (highInterest) {
       fill(meEeHiZsFir_, dccZsFIR(frame, firWeights_, firstFIRSample_, nullptr));
     } else {
       int v = dccZsFIR(frame, firWeights_, firstFIRSample_, nullptr);
       fill(meEeLiZsFir_, v);
       if (v < eeZsThr_) {
         eventError = true;
         ++nEeZsErrors_;
         pair<int, int> ru = dccCh(frame.id());
         if (isRuComplete_[ru.first - 1][ru.second - 1])
           ++nEeZsErrorsType1_;
         if (nEeZsErrors_ < 3) {
           srApplicationErrorLog_ << event.id() << ", "
                                  << "RU " << frame.id() << ", "
                                  << "DCC " << ru.first << " Ch : " << ru.second << ": "
                                  << "LI channel under ZS threshold.\n";
         }
         if (nEeZsErrors_ == 3) {
           srApplicationErrorLog_ << event.id() << ": "
                                  << "more ZS errors for this event...\n";
         }
       }
     }
   }  //next ZS digi.
 
   for (int iEeZ = 0; iEeZ < 2; ++iEeZ) {
     for (int iEeX = 0; iEeX < nEeX; ++iEeX) {
       for (int iEeY = 0; iEeY < nEeY; ++iEeY) {
         fill(meChOcc_,
              EExtalCoor[iEeX][iEeY][iEeZ].first,
              EExtalCoor[iEeX][iEeY][iEeZ].second,
              EEcrystalShot[iEeX][iEeY][iEeZ] ? 1 : 0);
       }
     }
   }
 
   for (int iZ0 = 0; iZ0 < nEndcaps; ++iZ0) {
     for (int iX0 = 0; iX0 < nEeX; ++iX0) {
       for (int iY0 = 0; iY0 < nEeY; ++iY0) {
         double recE = (*eeEnergies)[iZ0][iX0][iY0].recE;
         if (recE == -numeric_limits<double>::max())
           continue;  //not a crystal or ZS
         fill(meEeRecE_, (*eeEnergies)[iZ0][iX0][iY0].recE);
 
         fill(meEeEMean_, ievt_ + 1, (*eeEnergies)[iZ0][iX0][iY0].recE);
 
         if (withEeSimHit_) {
           if (!(*eeEnergies)[iZ0][iX0][iY0].simHit) {  //noise only crystal channel
             fill(meEeNoise_, (*eeEnergies)[iZ0][iX0][iY0].noZsRecE);
           } else {
             fill(meEeSimE_, (*eeEnergies)[iZ0][iX0][iY0].simE);
             fill(meEeRecEHitXtal_, (*eeEnergies)[iZ0][iX0][iY0].recE);
           }
           fill(meEeRecVsSimE_, (*eeEnergies)[iZ0][iX0][iY0].simE, (*eeEnergies)[iZ0][iX0][iY0].recE);
           fill(meEeNoZsRecVsSimE_, (*eeEnergies)[iZ0][iX0][iY0].simE, (*eeEnergies)[iZ0][iX0][iY0].noZsRecE);
         }
       }
     }
   }
 
   int EEZs1RuCount[2][20][20];
   int EEFullRuCount[2][20][20];
   int EEForcedRuCount[2][20][20];
   for (int iZ(0); iZ < 2; iZ++) {
     for (int iX(0); iX < 20; iX++) {
       for (int iY(0); iY < 20; iY++) {
         EEZs1RuCount[iZ][iX][iY] = 0;
         EEFullRuCount[iZ][iX][iY] = 0;
         EEForcedRuCount[iZ][iX][iY] = 0;
       }
     }
   }
 
   nEeFROCnt_ = 0;
   char eeSrfMark[2][20][20];
   bzero(eeSrfMark, sizeof(eeSrfMark));
   //Filling RU histo
   for (EESrFlagCollection::const_iterator it = eeSrFlags_->begin(); it != eeSrFlags_->end(); ++it) {
     const EESrFlag& srf = *it;
     // srf.id() is EcalScDetId; 1 <= ix <= 20 1 <= iy <= 20
     int iX = srf.id().ix();
     int iY = srf.id().iy();
     int zside = srf.id().zside();  //-1 for EE-, +1 for EE+
     if (iX < 1 || iY > 100)
       throw cms::Exception("EcalSelectiveReadoutValidation")
           << "Found an endcap SRF with an invalid det ID: " << srf.id() << ".\n";
     ++eeSrfMark[zside > 0 ? 1 : 0][iX - 1][iY - 1];
     if (eeSrfMark[zside > 0 ? 1 : 0][iX - 1][iY - 1] > 1)
       throw cms::Exception("EcalSelectiveReadoutValidation") << "Duplicate SRF for supercrystal " << srf.id() << ".\n";
     int flag = srf.value() & ~EcalSrFlag::SRF_FORCED_MASK;
     if (flag == EcalSrFlag::SRF_ZS1) {
       EEZs1RuCount[zside > 0 ? 1 : 0][iX - 1][iY - 1] += 1;
     }
 
     if (flag == EcalSrFlag::SRF_FULL) {
       EEFullRuCount[zside > 0 ? 1 : 0][iX - 1][iY - 1] += 1;
       ++nEeFROCnt_;
     }
 
     if (srf.value() & EcalSrFlag::SRF_FORCED_MASK) {
       EEForcedRuCount[zside > 0 ? 1 : 0][iX - 1][iY - 1] += 1;
     }
   }
   for (int iZ(0); iZ < 2; iZ++) {
     int xOffset(iZ == 0 ? -40 : 20);
     for (int iX(0); iX < 20; iX++) {
       for (int iY(0); iY < 20; iY++) {
         int GraphX = (iX + 1) + xOffset;
         int GraphY = (iY + 1);
         fill(meZs1Ru_, GraphX, GraphY, EEZs1RuCount[iZ][iX][iY]);
         fill(meFullRoRu_, GraphX, GraphY, EEFullRuCount[iZ][iX][iY]);
         fill(meForcedRu_, GraphX, GraphY, EEForcedRuCount[iZ][iX][iY]);
       }
     }
   }
 
   if (eventError)
     srApplicationErrorLog_ << event.id() << ": " << nEeZsErrors_
                            << " ZS-flagged EE channels under "
                               "the ZS threshold, whose "
                            << nEeZsErrorsType1_ << " in a complete RU.\n";
 }  //end of analyzeEE

◆ analyzeTP()

void EcalSelectiveReadoutValidation::analyzeTP	(	const edm::Event &	event,
		const edm::EventSetup &	es
	)

private

Trigger primitive analysis. To be called for each event.

Parameters

event	EDM event
es	event setup

Definition at line 1487 of file EcalSelectiveReadoutValidation.cc.

References cIndex2iTtEta(), cIndex2iTtPhi(), DetId::Ecal, EcalBarrel, EcalEndcap, fill(), edm::EventSetup::getData(), dqm::impl::MonitorElement::getTProfile(), mps_fire::i, l1tTowerCalibrationProducer_cfi::iEta, iTtEta2cIndex(), iTtPhi2cIndex(), lutGrpHandle, lutMapHandle, meForcedTtf_, meHiTtf_, meLiTtf_, meMiTtf_, meTp_, meTpMap_, meTpVsEtSum_, meTtf_, meTtfVsEtSum_, meTtfVsTp_, nTtEta, nTtPhi, physHandle, nano_mu_digi_cff::rawId, ecaldqm::towerId(), tpInGeV_, tps_, and ttEtSums.

Referenced by analyze().

                                                                                              {
   int TTFlagCount[8];
   int LiTTFlagCount[nTtEta][nTtPhi];
   int MiTTFlagCount[nTtEta][nTtPhi];
   int HiTTFlagCount[nTtEta][nTtPhi];
   for (int iTTFlag(0); iTTFlag < 8; iTTFlag++) {
     TTFlagCount[iTTFlag] = 0;
   }
   for (int iTtEta(0); iTtEta < nTtEta; iTtEta++) {
     for (int iTtPhi(0); iTtPhi < nTtPhi; iTtPhi++) {
       LiTTFlagCount[iTtEta][iTtPhi] = 0;
       MiTTFlagCount[iTtEta][iTtPhi] = 0;
       HiTTFlagCount[iTtEta][iTtPhi] = 0;
     }
   }
   int tpEtCount[100];
   for (int iEt(0); iEt < 100; iEt++) {
     tpEtCount[iEt] = 0;
   }
 
   const EcalTPGPhysicsConstMap& physMap = es.getData(physHandle).getMap();
 
   const EcalTPGGroups::EcalTPGGroupsMap& lutGrpMap = es.getData(lutGrpHandle).getMap();
 
   const EcalTPGLutIdMap::EcalTPGLutMap& lutMap = es.getData(lutMapHandle).getMap();
 
   EcalTPGPhysicsConstMapIterator ebItr(physMap.find(DetId(DetId::Ecal, EcalBarrel).rawId()));
   double lsb10bitsEB(ebItr == physMap.end() ? 0. : ebItr->second.EtSat / 1024.);
   EcalTPGPhysicsConstMapIterator eeItr(physMap.find(DetId(DetId::Ecal, EcalEndcap).rawId()));
   double lsb10bitsEE(eeItr == physMap.end() ? 0. : eeItr->second.EtSat / 1024.);
 
   for (EcalTrigPrimDigiCollection::const_iterator it = tps_->begin(); it != tps_->end(); ++it) {
     double tpEt;
     if (tpInGeV_) {
       EcalTrigTowerDetId const& towerId(it->id());
       unsigned int ADC = it->compressedEt();
 
       double lsb10bits(0.);
       if (towerId.subDet() == EcalBarrel)
         lsb10bits = lsb10bitsEB;
       else if (towerId.subDet() == EcalEndcap)
         lsb10bits = lsb10bitsEE;
 
       int tpg10bits = 0;
       EcalTPGGroups::EcalTPGGroupsMapItr itgrp = lutGrpMap.find(towerId.rawId());
       uint32_t lutGrp = 999;
       if (itgrp != lutGrpMap.end())
         lutGrp = itgrp->second;
 
       EcalTPGLutIdMap::EcalTPGLutMapItr itLut = lutMap.find(lutGrp);
       if (itLut != lutMap.end()) {
         const unsigned int* lut = (itLut->second).getLut();
         for (unsigned int i = 0; i < 1024; i++)
           if (ADC == (0xff & lut[i])) {
             tpg10bits = i;
             break;
           }
       }
 
       tpEt = lsb10bits * tpg10bits;
     } else {
       tpEt = it->compressedEt();
     }
     int iEta = it->id().ieta();
     int iEta0 = iTtEta2cIndex(iEta);
     int iPhi = it->id().iphi();
     int iPhi0 = iTtPhi2cIndex(iPhi);
     double etSum = ttEtSums[iEta0][iPhi0];
 
     int iE = meTp_->getTProfile()->FindFixBin(tpEt);
     if ((iE >= 0) && (iE < 100)) {
       ++tpEtCount[iE];
     } else {
       // FindFixBin might return an overflow bin (outside tpEtCount range).
       // To prevent a memory overflow / segfault, these values are ignored.
       //std::cout << "EcalSelectiveReadoutValidation: Invalid iE value: " << iE << std::endl;
     }
 
     fill(meTpVsEtSum_, etSum, tpEt);
     ++TTFlagCount[it->ttFlag()];
     if ((it->ttFlag() & 0x3) == 0) {
       LiTTFlagCount[iEta0][iPhi0] += 1;
     } else if ((it->ttFlag() & 0x3) == 1) {
       MiTTFlagCount[iEta0][iPhi0] += 1;
     } else if ((it->ttFlag() & 0x3) == 3) {
       HiTTFlagCount[iEta0][iPhi0] += 1;
     }
     if ((it->ttFlag() & 0x4)) {
       fill(meForcedTtf_, iEta, iPhi);
     }
 
     fill(meTtfVsTp_, tpEt, it->ttFlag());
     fill(meTtfVsEtSum_, etSum, it->ttFlag());
     fill(meTpMap_, iEta, iPhi, tpEt, 1.);
   }
 
   for (int ittflag(0); ittflag < 8; ittflag++) {
     fill(meTtf_, ittflag, TTFlagCount[ittflag]);
   }
   for (int iTtEta(0); iTtEta < nTtEta; iTtEta++) {
     for (int iTtPhi(0); iTtPhi < nTtPhi; iTtPhi++) {
       fill(meLiTtf_, cIndex2iTtEta(iTtEta), cIndex2iTtPhi(iTtPhi), LiTTFlagCount[iTtEta][iTtPhi]);
       fill(meMiTtf_, cIndex2iTtEta(iTtEta), cIndex2iTtPhi(iTtPhi), MiTTFlagCount[iTtEta][iTtPhi]);
       fill(meHiTtf_, cIndex2iTtEta(iTtEta), cIndex2iTtPhi(iTtPhi), HiTTFlagCount[iTtEta][iTtPhi]);
     }
   }
   if (tpInGeV_) {
     for (int iE(0); iE < 100; iE++) {
       fill(meTp_, iE, tpEtCount[iE]);
     }
   } else {
     for (int iE(0); iE < 40; iE++) {
       fill(meTp_, iE, tpEtCount[iE]);
     }
   }
 }

◆ book1D()

EcalSelectiveReadoutValidation::MonitorElement * EcalSelectiveReadoutValidation::book1D	(	DQMStore::IBooker &	ibook,
		const std::string &	name,
		const std::string &	title,
		int	nbins,
		double	xmin,
		double	xmax
	)

private

Definition at line 1981 of file EcalSelectiveReadoutValidation.cc.

References dqm::implementation::IBooker::book1D(), Exception, Skims_PA_cff::name, LaserClient_cfi::nbins, registerHist(), mps_fire::result, runGCPTkAlMap::title, TrackerOfflineValidation_Dqm_cff::xmax, and TrackerOfflineValidation_Dqm_cff::xmin.

Referenced by bookHistograms().

                                                                                                                 {
   if (!registerHist(name, title))
     return nullptr;  //this histo is disabled
   MonitorElement* result = ibook.book1D(name, title, nbins, xmin, xmax);
   if (result == nullptr) {
     throw cms::Exception("Histo") << "Failed to book histogram " << name;
   }
   return result;
 }

◆ book2D()

EcalSelectiveReadoutValidation::MonitorElement * EcalSelectiveReadoutValidation::book2D	(	DQMStore::IBooker &	ibook,
		const std::string &	name,
		const std::string &	title,
		int	nxbins,
		double	xmin,
		double	xmax,
		int	nybins,
		double	ymin,
		double	ymax
	)

private

Definition at line 1992 of file EcalSelectiveReadoutValidation.cc.

References dqm::implementation::IBooker::book2D(), Exception, Skims_PA_cff::name, registerHist(), mps_fire::result, runGCPTkAlMap::title, TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, L1TOccupancyClient_cfi::ymax, and L1TOccupancyClient_cfi::ymin.

Referenced by bookHistograms().

                                                                                                     {
   if (!registerHist(name, title))
     return nullptr;  //this histo is disabled
   MonitorElement* result = ibook.book2D(name, title, nxbins, xmin, xmax, nybins, ymin, ymax);
   if (result == nullptr) {
     throw cms::Exception("Histo") << "Failed to book histogram " << name;
   }
   return result;
 }

◆ bookFloat()

EcalSelectiveReadoutValidation::MonitorElement * EcalSelectiveReadoutValidation::bookFloat	(	DQMStore::IBooker &	ibook,
		const std::string &	name
	)

private

Wrappers to the book methods of the DQMStore DQM histogramming interface.

Definition at line 1970 of file EcalSelectiveReadoutValidation.cc.

References dqm::implementation::IBooker::bookFloat(), Exception, Skims_PA_cff::name, registerHist(), and mps_fire::result.

Referenced by bookHistograms().

                                                                                                                  {
   if (!registerHist(name, ""))
     return nullptr;  //this histo is disabled
   MonitorElement* result = ibook.bookFloat(name);
   if (result == nullptr) {
     throw cms::Exception("DQM") << "Failed to book integer DQM monitor element" << name;
   }
   return result;
 }

◆ bookHistograms()

void EcalSelectiveReadoutValidation::bookHistograms	(	DQMStore::IBooker &	i,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Implements DQMOneEDAnalyzer<>.

Definition at line 856 of file EcalSelectiveReadoutValidation.cc.

References availableHistList_, book1D(), book2D(), bookFloat(), bookProfile(), bookProfile2D(), visDQMUpload::buf, ebZsThr_, eeZsThr_, histList_, meChOcc_, meCompleteZSCnt_, meCompleteZSMap_, meCompleteZSRateMap_, meDccHiVol_, meDccLiVol_, meDccVol_, meDccVolFromData_, meDroppedFROCnt_, meDroppedFROMap_, meDroppedFRORateMap_, meEbEMean_, meEbFullRoCnt_, meEbHiZsFir_, meEbLiZsFir_, meEbNoise_, meEbNoZsRecVsSimE_, meEbRecE_, meEbRecEHitXtal_, meEbRecVsSimE_, meEbSimE_, meEbZsErrCnt_, meEbZsErrType1Cnt_, meEeEMean_, meEeFullRoCnt_, meEeHiZsFir_, meEeLiZsFir_, meEeNoise_, meEeNoZsRecVsSimE_, meEeRecE_, meEeRecEHitXtal_, meEeRecVsSimE_, meEeSimE_, meEeZsErrCnt_, meEeZsErrType1Cnt_, meForcedRu_, meForcedTtf_, meFullRoCnt_, meFullRoRu_, meHiTtf_, meIncompleteFROCnt_, meIncompleteFROMap_, meIncompleteFRORateMap_, meL1aRate_, meLiTtf_, meMiTtf_, meSRFlagsConsistency_, meTp_, meTpMap_, meTpVsEtSum_, meTtf_, meTtfVsEtSum_, meTtfVsTp_, meVol_, meVolB_, meVolBHI_, meVolBLI_, meVolE_, meVolEHI_, meVolELI_, meVolHI_, meVolLI_, meZs1Ru_, meZsErrCnt_, meZsErrType1Cnt_, nDccs_, printAvailableHists(), alignCSCRings::s, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, runGCPTkAlMap::title, and tpInGeV_.

                                                                           {
   ibooker.setCurrentFolder("EcalDigisV/SelectiveReadout");
 
   {
     auto scope = DQMStore::IBooker::UseRunScope(ibooker);
     meL1aRate_ = bookFloat(ibooker, "l1aRate_");
   }
 
   meDccVol_ = bookProfile(ibooker,
                           "hDccVol",  //"EcalDccEventSizeComputed",
                           "ECAL DCC event fragment size;Dcc id; "
                           "<Event size> (kB)",
                           nDccs_,
                           .5,
                           .5 + nDccs_);
 
   meDccLiVol_ = bookProfile(ibooker,
                             "hDccLiVol",
                             "LI channel payload per DCC;Dcc id; "
                             "<Event size> (kB)",
                             nDccs_,
                             .5,
                             .5 + nDccs_);
 
   meDccHiVol_ = bookProfile(ibooker,
                             "hDccHiVol",
                             "HI channel payload per DCC;Dcc id; "
                             "<Event size> (kB)",
                             nDccs_,
                             .5,
                             .5 + nDccs_);
 
   meDccVolFromData_ = bookProfile(ibooker,
                                   "hDccVolFromData",  //"EcalDccEventSize",
                                   "ECAL DCC event fragment size;Dcc id; "
                                   "<Event size> (kB)",
                                   nDccs_,
                                   .5,
                                   .5 + nDccs_);
 
   meVolBLI_ = book1D(ibooker,
                      "hVolBLI",  // "EBLowInterestPayload",
                      "ECAL Barrel low interest crystal data payload;"
                      "Event size (kB);Nevts",
                      100,
                      0.,
                      200.);
 
   meVolELI_ = book1D(ibooker,
                      "hVolELI",  //"EELowInterestPayload",
                      "Endcap low interest crystal data payload;"
                      "Event size (kB);Nevts",
                      100,
                      0.,
                      200.);
 
   meVolLI_ = book1D(ibooker,
                     "hVolLI",  //"EcalLowInterestPayload",
                     "ECAL low interest crystal data payload;"
                     "Event size (kB);Nevts",
                     100,
                     0.,
                     200.);
 
   meVolBHI_ = book1D(ibooker,
                      "hVolBHI",  //"EBHighInterestPayload",
                      "Barrel high interest crystal data payload;"
                      "Event size (kB);Nevts",
                      100,
                      0.,
                      200.);
 
   meVolEHI_ = book1D(ibooker,
                      "hVolEHI",  //"EEHighInterestPayload",
                      "Endcap high interest crystal data payload;"
                      "Event size (kB);Nevts",
                      100,
                      0.,
                      200.);
 
   meVolHI_ = book1D(ibooker,
                     "hVolHI",  //"EcalHighInterestPayload",
                     "ECAL high interest crystal data payload;"
                     "Event size (kB);Nevts",
                     100,
                     0.,
                     200.);
 
   meVolB_ = book1D(ibooker,
                    "hVolB",  //"EBEventSize",
                    "Barrel data volume;Event size (kB);Nevts",
                    100,
                    0.,
                    200.);
 
   meVolE_ = book1D(ibooker,
                    "hVolE",  //"EEEventSize",
                    "Endcap data volume;Event size (kB);Nevts",
                    100,
                    0.,
                    200.);
 
   meVol_ = book1D(ibooker,
                   "hVol",  //"EcalEventSize",
                   "ECAL data volume;Event size (kB);Nevts",
                   100,
                   0.,
                   200.);
 
   meChOcc_ = bookProfile2D(ibooker,
                            "h2ChOcc",  //"EcalChannelOccupancy",
                            "ECAL crystal channel occupancy after zero suppression;"
                            "iX -200 / iEta / iX + 100;"
                            "iY / iPhi (starting from -10^{o}!);"
                            "Event count rate",
                            401,
                            -200.5,
                            200.5,
                            360,
                            .5,
                            360.5);
 
   //TP
   string tpUnit;
   if (tpInGeV_)
     tpUnit = string("GeV");
   else
     tpUnit = string("TP hw unit");
   string title;
   title = string("Trigger primitive TT E_{T};E_{T} ") + tpUnit + string(";Event Count");
   meTp_ = bookProfile(ibooker,
                       "hTp",  //"EcalTriggerPrimitiveEt",
                       title,
                       (tpInGeV_ ? 100 : 40),
                       0.,
                       (tpInGeV_ ? 10. : 40.));
 
   meTtf_ = bookProfile(ibooker,
                        "hTtf",  //"EcalTriggerTowerFlag",
                        "Trigger primitive TT flag;Flag number;Event count",
                        8,
                        -.5,
                        7.5);
 
   title = string("Trigger tower flag vs TP;E_{T}(TT) (") + tpUnit + string(");Flag number");
   meTtfVsTp_ = book2D(ibooker, "h2TtfVsTp", title, 100, 0., (tpInGeV_ ? 10. : 40.), 8, -.5, 7.5);
 
   meTtfVsEtSum_ = book2D(ibooker,
                          "h2TtfVsEtSum",
                          "Trigger tower flag vs #sumE_{T};"
                          "E_{T}(TT) (GeV);"
                          "TTF",
                          100,
                          0.,
                          10.,
                          8,
                          -.5,
                          7.5);
   title = string(
               "Trigger primitive Et (TP) vs #sumE_{T};"
               "E_{T} (sum) (GeV);"
               "E_{T} (TP) (") +
           tpUnit + string(")");
 
   meTpVsEtSum_ = book2D(ibooker, "h2TpVsEtSum", title, 100, 0., 10., 100, 0., (tpInGeV_ ? 10. : 40.));
 
   title = string(
               "Trigger primitive E_{T};"
               "iEta;"
               "iPhi;"
               "E_{T} (TP) (") +
           tpUnit + string(")");
   meTpMap_ = bookProfile2D(ibooker, "h2Tp", title, 57, -28.5, 28.5, 72, .5, 72.5);
 
   //SRF
   meFullRoRu_ = book2D(ibooker,
                        "h2FRORu",  //"EcalFullReadoutSRFlagMap",
                        "Full Read-out readout unit;"
                        "iX - 40 / iEta / iX + 20;"
                        "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                        "Event count",
                        80,
                        -39.5,
                        40.5,
                        72,
                        .5,
                        72.5);
 
   meFullRoCnt_ = book1D(ibooker,
                         "hFROCnt",
                         "Number of Full-readout-flagged readout units;"
                         "FRO RU count;Event count",
                         300,
                         -.5,
                         299.5);
 
   meEbFullRoCnt_ = book1D(ibooker,
                           "hEbFROCnt",
                           "Number of EB Full-readout-flagged readout units;"
                           "FRO RU count;Event count",
                           200,
                           -.5,
                           199.5);
 
   meEeFullRoCnt_ = book1D(ibooker,
                           "hEeFROCnt",
                           "Number of EE Full-readout-flagged readout units;"
                           "FRO RU count;Event count",
                           200,
                           -.5,
                           199.5);
 
   meZs1Ru_ = book2D(ibooker,
                     "h2Zs1Ru",  //"EbZeroSupp1SRFlagMap",
                     "Readout unit with ZS-thr-1 flag;"
                     "iX - 40 / iEta / iX + 20;"
                     "iY0 / iPhi0 (iPhi = 1 at phi = 0 rad);"
                     "Event count",
                     80,
                     -39.5,
                     40.5,
                     72,
                     .5,
                     72.5);
 
   meForcedRu_ = book2D(ibooker,
                        "h2ForcedRu",  //"EcalReadoutUnitForcedBitMap",
                        "ECAL readout unit with forced bit of SR flag on;"
                        "iX - 40 / iEta / iX + 20;"
                        "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                        "Event count",
                        80,
                        -39.5,
                        40.5,
                        72,
                        .5,
                        72.5);
 
   meLiTtf_ = bookProfile2D(ibooker,
                            "h2LiTtf",  //"EcalLowInterestTriggerTowerFlagMap",
                            "Low interest trigger tower flags;"
                            "iEta;"
                            "iPhi;"
                            "Event count",
                            57,
                            -28.5,
                            28.5,
                            72,
                            .5,
                            72.5);
 
   meMiTtf_ = bookProfile2D(ibooker,
                            "h2MiTtf",  //"EcalMidInterestTriggerTowerFlagMap",
                            "Mid interest trigger tower flags;"
                            "iEta;"
                            "iPhi;"
                            "Event count",
                            57,
                            -28.5,
                            28.5,
                            72,
                            .5,
                            72.5);
 
   meHiTtf_ = bookProfile2D(ibooker,
                            "h2HiTtf",  //"EcalHighInterestTriggerTowerFlagMap",
                            "High interest trigger tower flags;"
                            "iEta;"
                            "iPhi;"
                            "Event count",
                            57,
                            -28.5,
                            28.5,
                            72,
                            .5,
                            72.5);
 
   meForcedTtf_ = book2D(ibooker,
                         "h2ForcedTtf",  //"EcalTtfForcedBitMap",
                         "Trigger tower flags with forced bit set;"
                         "iEta;"
                         "iPhi;"
                         "Event count",
                         57,
                         -28.5,
                         28.5,
                         72,
                         .5,
                         72.5);
 
   const float ebMinNoise = -1.;
   const float ebMaxNoise = 1.;
 
   const float eeMinNoise = -1.;
   const float eeMaxNoise = 1.;
 
   const float ebMinE = ebMinNoise;
   const float ebMaxE = ebMaxNoise;
 
   const float eeMinE = eeMinNoise;
   const float eeMaxE = eeMaxNoise;
 
   const int evtMax = 500;
 
   meEbRecE_ = book1D(ibooker, "hEbRecE", "Crystal reconstructed energy;E (GeV);Event count", 100, ebMinE, ebMaxE);
 
   meEbEMean_ = bookProfile(ibooker, "hEbEMean", "EE <E_hit>;event #;<E_hit> (GeV)", evtMax, .5, evtMax + .5);
 
   meEbNoise_ = book1D(ibooker,
                       "hEbNoise",
                       "Crystal noise "
                       "(rec E of crystal without deposited energy)"
                       ";Rec E (GeV);Event count",
                       100,
                       ebMinNoise,
                       ebMaxNoise);
 
   meEbLiZsFir_ = book1D(ibooker,
                         "zsEbLiFIRemu",
                         "Emulated ouput of ZS FIR filter for EB "
                         "low interest crystals;"
                         "ADC count*4;"
                         "Event count",
                         60,
                         -30,
                         30);
 
   meEbHiZsFir_ = book1D(ibooker,
                         "zsEbHiFIRemu",
                         "Emulated ouput of ZS FIR filter for EB "
                         "high interest crystals;"
                         "ADC count*4;"
                         "Event count",
                         60,
                         -30,
                         30);
 
   //TODO: Fill this histogram...
   //   meEbIncompleteRUZsFir_ = book1D(ibooker, "zsEbIncompleteRUFIRemu",
   //                                   "Emulated ouput of ZS FIR filter for EB "
   //                                   "incomplete FRO-flagged RU;"
   //                                   "ADC count*4;"
   //                                   "Event count",
   //                                   60, -30, 30);
 
   meEbSimE_ = book1D(ibooker, "hEbSimE", "EB hit crystal simulated energy", 100, ebMinE, ebMaxE);
 
   meEbRecEHitXtal_ = book1D(ibooker, "hEbRecEHitXtal", "EB rec energy of hit crystals", 100, ebMinE, ebMaxE);
 
   meEbRecVsSimE_ = book2D(ibooker,
                           "hEbRecVsSimE",
                           "Crystal simulated vs reconstructed energy;"
                           "Esim (GeV);Erec GeV);Event count",
                           100,
                           ebMinE,
                           ebMaxE,
                           100,
                           ebMinE,
                           ebMaxE);
 
   meEbNoZsRecVsSimE_ = book2D(ibooker,
                               "hEbNoZsRecVsSimE",
                               "Crystal no-zs simulated vs reconstructed "
                               "energy;"
                               "Esim (GeV);Erec GeV);Event count",
                               100,
                               ebMinE,
                               ebMaxE,
                               100,
                               ebMinE,
                               ebMaxE);
 
   meEeRecE_ = book1D(ibooker,
                      "hEeRecE",
                      "EE crystal reconstructed energy;E (GeV);"
                      "Event count",
                      100,
                      eeMinE,
                      eeMaxE);
 
   meEeEMean_ = bookProfile(ibooker, "hEeEMean", "<E_{EE hit}>;event;<E_{hit}> (GeV)", evtMax, .5, evtMax + .5);
 
   meEeNoise_ = book1D(ibooker,
                       "hEeNoise",
                       "EE crystal noise "
                       "(rec E of crystal without deposited energy);"
                       "E (GeV);Event count",
                       200,
                       eeMinNoise,
                       eeMaxNoise);
 
   meEeLiZsFir_ = book1D(ibooker,
                         "zsEeLiFIRemu",
                         "Emulated ouput of ZS FIR filter for EE "
                         "low interest crystals;"
                         "ADC count*4;"
                         "Event count",
                         60,
                         -30,
                         30);
 
   meEeHiZsFir_ = book1D(ibooker,
                         "zsEeHiFIRemu",
                         "Emulated ouput of ZS FIR filter for EE "
                         "high interest crystals;"
                         "ADC count*4;"
                         "Event count",
                         60,
                         -30,
                         30);
 
   //TODO: Fill this histogram...
   //   meEeIncompleteRUZsFir_ = book1D(ibooker, "zsEeIncompleteRUFIRemu",
   //                                     "Emulated ouput of ZS FIR filter for EE "
   //                                     "incomplete FRO-flagged RU;"
   //                                     "ADC count*4;"
   //                                     "Event count",
   //                                     60, -30, 30);
 
   meEeSimE_ = book1D(ibooker, "hEeSimE", "EE hit crystal simulated energy", 100, eeMinE, eeMaxE);
 
   meEeRecEHitXtal_ = book1D(ibooker, "hEeRecEHitXtal", "EE rec energy of hit crystals", 100, eeMinE, eeMaxE);
 
   meEeRecVsSimE_ = book2D(ibooker,
                           "hEeRecVsSimE",
                           "EE crystal simulated vs reconstructed energy;"
                           "Esim (GeV);Erec GeV);Event count",
                           100,
                           eeMinE,
                           eeMaxE,
                           100,
                           eeMinE,
                           eeMaxE);
 
   meEeNoZsRecVsSimE_ = book2D(ibooker,
                               "hEeNoZsRecVsSimE",
                               "EE crystal no-zs simulated vs "
                               "reconstructed "
                               "energy;Esim (GeV);Erec GeV);Event count",
                               100,
                               eeMinE,
                               eeMaxE,
                               100,
                               eeMinE,
                               eeMaxE);
 
   meSRFlagsConsistency_ = book2D(ibooker,
                                  "hSRAlgoErrorMap",
                                  "TTFlags and SR Flags mismatch;"
                                  "iX - 40 / iEta / iX + 20;"
                                  "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                                  "Event count",
                                  80,
                                  -39.5,
                                  40.5,
                                  72,
                                  .5,
                                  72.5);
 
   //Readout Units histos (interest/Ncrystals)
   meIncompleteFROMap_ = book2D(ibooker,
                                "hIncompleteFROMap",
                                "Incomplete full-readout-flagged readout units;"
                                "iX - 40 / iEta / iX + 20;"
                                "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                                "Event count",
                                80,
                                -39.5,
                                40.5,
                                72,
                                .5,
                                72.5);
 
   meIncompleteFROCnt_ = book1D(ibooker,
                                "hIncompleteFROCnt",
                                "Number of incomplete full-readout-flagged "
                                "readout units;"
                                "Number of RUs;Event count;",
                                200,
                                -.5,
                                199.5);
 
   meIncompleteFRORateMap_ = bookProfile2D(ibooker,
                                           "hIncompleteFRORateMap",
                                           "Incomplete full-readout-flagged readout units;"
                                           "iX - 40 / iEta / iX + 20;"
                                           "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                                           "Incomplete error rate",
                                           80,
                                           -39.5,
                                           40.5,
                                           72,
                                           .5,
                                           72.5);
 
   meDroppedFROMap_ = book2D(ibooker,
                             "hDroppedFROMap",
                             "Dropped full-readout-flagged readout units;"
                             "iX - 40 / iEta / iX + 20;"
                             "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                             "Event count",
                             80,
                             -39.5,
                             40.5,
                             72,
                             .5,
                             72.5);
 
   meDroppedFROCnt_ = book1D(ibooker,
                             "hDroppedFROCnt",
                             "Number of dropped full-readout-flagged "
                             "RU count;RU count;Event count",
                             200,
                             -.5,
                             199.5);
 
   meCompleteZSCnt_ = book1D(ibooker,
                             "hCompleteZsCnt",
                             "Number of zero-suppressed-flagged RU fully "
                             "readout;"
                             "RU count;Event count",
                             200,
                             -.5,
                             199.5);
 
   stringstream buf;
   buf << "Number of LI EB channels below the " << ebZsThr_ / 4.
       << " ADC count ZS threshold;"
          "Channel count;Event count",
       meEbZsErrCnt_ = book1D(ibooker, "hEbZsErrCnt", buf.str(), 200, -.5, 199.5);
 
   buf.str("");
   buf << "Number of LI EE channels below the " << eeZsThr_ / 4.
       << " ADC count ZS theshold;"
          "Channel count;Event count",
       meEeZsErrCnt_ = book1D(ibooker, "hEeZsErrCnt", buf.str(), 200, -.5, 199.5);
 
   meZsErrCnt_ = book1D(ibooker,
                        "hZsErrCnt",
                        "Number of LI channels below the ZS threshold;"
                        "Channel count;Event count",
                        200,
                        -.5,
                        199.5);
 
   meEbZsErrType1Cnt_ = book1D(ibooker,
                               "hEbZsErrType1Cnt",
                               "Number of EB channels below the ZS "
                               "threshold in a LI but fully readout RU;"
                               "Channel count;Event count;",
                               200,
                               -.5,
                               199.5);
 
   meEeZsErrType1Cnt_ = book1D(ibooker,
                               "hEeZsErrType1Cnt",
                               "Number EE channels below the ZS threshold"
                               " in a LI but fully readout RU;"
                               "Channel count;Event count",
                               200,
                               -.5,
                               199.5);
 
   meZsErrType1Cnt_ = book1D(ibooker,
                             "hZsErrType1Cnt",
                             "Number of LI channels below the ZS threshold "
                             "in a LI but fully readout RU;"
                             "Channel count;Event count",
                             200,
                             -.5,
                             199.5);
 
   meDroppedFRORateMap_ = bookProfile2D(ibooker,
                                        "hDroppedFRORateMap",
                                        "Dropped full-readout-flagged readout units"
                                        "iX - 40 / iEta / iX + 20;"
                                        "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                                        "Dropping rate",
                                        80,
                                        -39.5,
                                        40.5,
                                        72,
                                        .5,
                                        72.5);
 
   meCompleteZSMap_ = book2D(ibooker,
                             "hCompleteZSMap",
                             "Complete zero-suppressed-flagged readout units;"
                             "iX - 40 / iEta / iX + 20;"
                             "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                             "Event count",
                             80,
                             -39.5,
                             40.5,
                             72,
                             .5,
                             72.5);
 
   meCompleteZSRateMap_ = bookProfile2D(ibooker,
                                        "hCompleteZSRate",
                                        "Complete zero-suppressed-flagged readout units;"
                                        "iX - 40 / iEta / iX + 20;"
                                        "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                                        "Completeness rate",
                                        80,
                                        -39.5,
                                        40.5,
                                        72,
                                        .5,
                                        72.5);
 
   //print list of available histograms (must be called after
   //the bookXX methods):
   printAvailableHists();
 
   //check the histList parameter:
   stringstream s;
   for (set<string>::iterator it = histList_.begin(); it != histList_.end(); ++it) {
     if (*it != string("all") && availableHistList_.find(*it) == availableHistList_.end()) {
       s << (s.str().empty() ? "" : ", ") << *it;
     }
   }
   if (!s.str().empty()) {
     LogWarning("Configuration") << "Parameter 'histList' contains some unknown histogram(s). "
                                    "Check spelling. Following name were not found: "
                                 << s.str();
   }
 }

◆ bookProfile()

EcalSelectiveReadoutValidation::MonitorElement * EcalSelectiveReadoutValidation::bookProfile	(	DQMStore::IBooker &	ibook,
		const std::string &	name,
		const std::string &	title,
		int	nbins,
		double	xmin,
		double	xmax
	)

private

Definition at line 2010 of file EcalSelectiveReadoutValidation.cc.

References dqm::implementation::IBooker::bookProfile(), Exception, Skims_PA_cff::name, LaserClient_cfi::nbins, registerHist(), mps_fire::result, runGCPTkAlMap::title, TrackerOfflineValidation_Dqm_cff::xmax, and TrackerOfflineValidation_Dqm_cff::xmin.

Referenced by bookHistograms().

                                                                                                                 {
   if (!registerHist(name, title))
     return nullptr;  //this histo is disabled
   MonitorElement* result = ibook.bookProfile(name, title, nbins, xmin, xmax, 0, 0, 0);
   if (result == nullptr) {
     throw cms::Exception("Histo") << "Failed to book histogram " << name;
   }
   return result;
 }

◆ bookProfile2D()

EcalSelectiveReadoutValidation::MonitorElement * EcalSelectiveReadoutValidation::bookProfile2D	(	DQMStore::IBooker &	ibook,
		const std::string &	name,
		const std::string &	title,
		int	nbinx,
		double	xmin,
		double	xmax,
		int	nbiny,
		double	ymin,
		double	ymax,
		const char *	option = `""`
	)

private

Definition at line 2021 of file EcalSelectiveReadoutValidation.cc.

References dqm::implementation::IBooker::bookProfile2D(), Exception, Skims_PA_cff::name, fileinputsource_cfi::option, registerHist(), mps_fire::result, runGCPTkAlMap::title, TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, L1TOccupancyClient_cfi::ymax, and L1TOccupancyClient_cfi::ymin.

Referenced by bookHistograms().

                                                                                                                   {
   if (!registerHist(name, title))
     return nullptr;  //this histo is disabled
   MonitorElement* result = ibook.bookProfile2D(name, title, nbinx, xmin, xmax, nbiny, ymin, ymax, 0, 0, 0, option);
   if (result == nullptr) {
     throw cms::Exception("Histo") << "Failed to book histogram " << name;
   }
   return result;
 }

◆ checkSrApplication()

template<class T >

void EcalSelectiveReadoutValidation::checkSrApplication	(	const edm::Event &	event,
		T &	srfs
	)

private

Checks application of SR decision by the DCC.

Parameters

event	event currently analyzed.
srfs	Selective readou flags

Definition at line 2318 of file EcalSelectiveReadoutValidation.cc.

References dccCh(), fill(), RemoveAddSevLevel::flag, getCrystalCount(), meCompleteZSMap_, meCompleteZSRateMap_, meDroppedFROMap_, meDroppedFRORateMap_, meIncompleteFROMap_, meIncompleteFRORateMap_, minDccId_, nCompleteZS_, nDroppedFRO_, nIncompleteFRO_, nPerRu_, ruGraphX(), ruGraphY(), srApplicationErrorLog_, EcalSrFlag::SRF_FULL, EcalSrFlag::SRF_ZS1, and EcalSrFlag::SRF_ZS2.

Referenced by analyze().

                                                                                       {
   typedef typename T::const_iterator SrFlagCollectionConstIt;
   typedef typename T::key_type MyRuDetIdType;
 
   for (SrFlagCollectionConstIt itSrf = srfs.begin(); itSrf != srfs.end(); ++itSrf) {
     int flag = itSrf->value() & ~EcalSrFlag::SRF_FORCED_MASK;
     pair<int, int> ru = dccCh(itSrf->id());
 
     if (flag == EcalSrFlag::SRF_FULL) {
       if (nPerRu_[ru.first - minDccId_][ru.second - 1] == getCrystalCount(ru.first, ru.second)) {  //no error
         fill(meIncompleteFRORateMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 0);
         fill(meDroppedFRORateMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 0);
       } else if (nPerRu_[ru.first - minDccId_][ru.second - 1] == 0) {  //tower dropped!
         fill(meIncompleteFRORateMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 0);
         fill(meDroppedFRORateMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 1);
         fill(meDroppedFROMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 1);
         ++nDroppedFRO_;
         srApplicationErrorLog_ << event.id() << ": Flag of RU " << itSrf->id() << " (DCC " << ru.first << " ch "
                                << ru.second << ") is 'Full readout' "
                                << "while none of its channel was read out\n";
       } else {  //tower partially read out
         fill(meIncompleteFRORateMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 1);
         fill(meDroppedFRORateMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 0);
         fill(meIncompleteFROMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 1);
         ++nIncompleteFRO_;
         srApplicationErrorLog_ << event.id() << ": Flag of RU" << itSrf->id() << " (DCC " << ru.first << " ch "
                                << ru.second << ") is 'Full readout' "
                                << "while only " << nPerRu_[ru.first - minDccId_][ru.second - 1] << " / "
                                << getCrystalCount(ru.first, ru.second) << " channels were read out.\n";
       }
     }
 
     if (flag == EcalSrFlag::SRF_ZS1 || flag == EcalSrFlag::SRF_ZS2) {
       if (nPerRu_[ru.first - minDccId_][ru.second - 1] == getCrystalCount(ru.first, ru.second)) {
         //ZS readout unit whose every channel was read
 
         fill(meCompleteZSMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()));
         fill(meCompleteZSRateMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 1);
 
         ++nCompleteZS_;
       } else {
         fill(meCompleteZSRateMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 0);
       }
     }
   }
 }

◆ cIndex2iEta()

int EcalSelectiveReadoutValidation::cIndex2iEta ( int i ) const

inlineprivate

converse of iEta2cIndex() method.

Definition at line 248 of file EcalSelectiveReadoutValidation.h.

References mps_fire::i.

Referenced by dccCh().

248 { return (i < 85) ? i - 85 : i - 84; }

mps_fire.i

i

Definition: mps_fire.py:429

◆ cIndex2iPhi()

int EcalSelectiveReadoutValidation::cIndex2iPhi ( int i ) const

inlineprivate

converse of iPhi2cIndex() method.

Definition at line 252 of file EcalSelectiveReadoutValidation.h.

References mps_fire::i.

Referenced by dccCh().

252 { return (i + 11) % 360; }

mps_fire.i

i

Definition: mps_fire.py:429

◆ cIndex2iTtEta()

int EcalSelectiveReadoutValidation::cIndex2iTtEta ( int i ) const

inlineprivate

converse of iTtEta2cIndex() method.

Definition at line 275 of file EcalSelectiveReadoutValidation.h.

References mps_fire::i.

Referenced by analyzeTP().

275 { return (i < 27) ? i - 28 : i - 27; }

mps_fire.i

i

Definition: mps_fire.py:429

◆ cIndex2iTtPhi()

int EcalSelectiveReadoutValidation::cIndex2iTtPhi ( int i ) const

inlineprivate

converse of iTtPhi2cIndex() method.

Definition at line 279 of file EcalSelectiveReadoutValidation.h.

References mps_fire::i.

Referenced by analyzeTP().

279 { return i + 1; }

mps_fire.i

i

Definition: mps_fire.py:429

◆ cIndex2iXY()

int EcalSelectiveReadoutValidation::cIndex2iXY ( int iX0 ) const

inlineprivate

converse of iXY2cIndex() method.

Definition at line 244 of file EcalSelectiveReadoutValidation.h.

244 { return iX0 + 1; }

◆ compareSrfColl()

template<class T >

void EcalSelectiveReadoutValidation::compareSrfColl	(	const edm::Event &	event,
		T &	srfFromData,
		T &	computedSrf
	)

private

Compares two SR flag collection, flags read from data and computed flags. Descripencies are recorded in relevant histogram and log file.

Template Parameters

T	collection type. Must be either an EESrFlagCollection or an EBSrFlagCollection.

Parameters

event	event currently analyzed. Used in logs.
srfFromData	SR flag collection read from data
compareSrf	SR flag collection computed from TTF by SRP emulation

Definition at line 2224 of file EcalSelectiveReadoutValidation.cc.

References dccCh(), fill(), logErrForDccs_, meSRFlagsConsistency_, minDccId_, ruGraphX(), ruGraphY(), and srpAlgoErrorLog_.

Referenced by analyze().

                                                                                                          {
   typedef typename T::const_iterator SrFlagCollectionConstIt;
   typedef typename T::key_type MyRuDetIdType;
   SrFlagCollectionConstIt itSrfFromData = srfFromData.begin();
   SrFlagCollectionConstIt itComputedSr = computedSrf.begin();
 
   while (itSrfFromData != srfFromData.end() || itComputedSr != computedSrf.end()) {
     MyRuDetIdType inconsistentRu = 0;
     bool inconsistent = false;
     if (itComputedSr == computedSrf.end() ||
         (itSrfFromData != srfFromData.end() && itSrfFromData->id() < itComputedSr->id())) {
       //computedSrf is missig a detid found in srfFromData
       pair<int, int> ch = dccCh(itSrfFromData->id());
       srpAlgoErrorLog_ << event.id() << ": " << itSrfFromData->id() << ", DCC " << ch.first << " ch " << ch.second
                        << " found in data (SRF:" << itSrfFromData->flagName()
                        << ") but not in the set of SRFs computed from the data TTF.\n";
       inconsistentRu = itSrfFromData->id();
       inconsistent = true;
       ++itSrfFromData;
     } else if (itSrfFromData == srfFromData.end() ||
                (itComputedSr != computedSrf.end() && itComputedSr->id() < itSrfFromData->id())) {
       //ebSrFlags is missing a detid found in computedSrf
       pair<int, int> ch = dccCh(itComputedSr->id());
       if (logErrForDccs_[ch.first - minDccId_]) {
         srpAlgoErrorLog_ << event.id() << ": " << itComputedSr->id() << ", DCC " << ch.first << " ch " << ch.second
                          << " not found in data. Computed SRF: " << itComputedSr->flagName() << ".\n";
         inconsistentRu = itComputedSr->id();
         inconsistent = true;
       }
       ++itComputedSr;
     } else {
       //*itSrfFromData and *itComputedSr has same detid
       if (itComputedSr->value() != itSrfFromData->value()) {
         pair<int, int> ch = dccCh(itSrfFromData->id());
         srpAlgoErrorLog_ << event.id() << ", " << itSrfFromData->id() << ", DCC " << ch.first << " ch " << ch.second
                          << ", SRF inconsistency: "
                          << "from data: " << itSrfFromData->flagName()
                          << ", computed from TTF: " << itComputedSr->flagName() << "\n";
         inconsistentRu = itComputedSr->id();
         inconsistent = true;
       }
       if (itComputedSr != computedSrf.end())
         ++itComputedSr;
       if (itSrfFromData != srfFromData.end())
         ++itSrfFromData;
     }
 
     if (inconsistent)
       fill(meSRFlagsConsistency_, ruGraphX(inconsistentRu), ruGraphY(inconsistentRu));
   }
 }

◆ configFirWeights()

void EcalSelectiveReadoutValidation::configFirWeights ( const std::vector< double > & weightsForZsFIR )

private

Configure DCC ZS FIR weights. Heuristic is used to determine if input weights are normalized weights or integer weights in the hardware representation.

Parameters

weightsForZsFIR weights from configuration file

Definition at line 2156 of file EcalSelectiveReadoutValidation.cc.

References Exception, firstFIRSample_, firWeights_, getFIRWeights(), mps_fire::i, createfilelist::int, dqm-mbProfile::log, and mathSSE::sqrt().

Referenced by EcalSelectiveReadoutValidation().

                                                                                            {
   bool notNormalized = false;
   bool notInt = false;
   for (unsigned i = 0; i < weightsForZsFIR.size(); ++i) {
     if (weightsForZsFIR[i] > 1.)
       notNormalized = true;
     if ((int)weightsForZsFIR[i] != weightsForZsFIR[i])
       notInt = true;
   }
   if (notInt && notNormalized) {
     throw cms::Exception("InvalidParameter") << "weigtsForZsFIR paramater values are not valid: they "
                                              << "must either be integer and uses the hardware representation "
                                              << "of the weights or less or equal than 1 and used the normalized "
                                              << "representation.";
   }
   LogInfo log("DccFir");
   if (notNormalized) {
     firWeights_ = vector<int>(weightsForZsFIR.size());
     for (unsigned i = 0; i < weightsForZsFIR.size(); ++i) {
       firWeights_[i] = (int)weightsForZsFIR[i];
     }
   } else {
     firWeights_ = getFIRWeights(weightsForZsFIR);
   }
 
   log << "Input weights for FIR: ";
   for (unsigned i = 0; i < weightsForZsFIR.size(); ++i) {
     log << weightsForZsFIR[i] << "\t";
   }
 
   double s2 = 0.;
   log << "\nActual FIR weights: ";
   for (unsigned i = 0; i < firWeights_.size(); ++i) {
     log << firWeights_[i] << "\t";
     s2 += firWeights_[i] * firWeights_[i];
   }
 
   s2 = sqrt(s2);
   log << "\nNormalized FIR weights after hw representation rounding: ";
   for (unsigned i = 0; i < firWeights_.size(); ++i) {
     log << firWeights_[i] / (double)(1 << 10) << "\t";
   }
 
   log << "\nFirst FIR sample: " << firstFIRSample_;
 }

◆ dccCh()

pair< int, int > EcalSelectiveReadoutValidation::dccCh ( const DetId & xtalId ) const

private

Retrieves the logical number of the DCC reading a given crystal channel.

Parameters

xtarId crystal channel identifier

Returns: the DCC logical number starting from 1. Retrieves the DCC channel reading out a crystal, the crystals of a barrel trigger tower or the crystals, of an endcap supercrystal.

Parameters

xtarId crystal channel, barrel trigger tower or endcap supercrystal identifier

Returns: pair of (DCC ID, DCC channel)

Definition at line 1787 of file EcalSelectiveReadoutValidation.cc.

References cIndex2iEta(), cIndex2iPhi(), EcalElectronicsId::dccId(), TauDecayModes::dec, hcalRecHitTable_cff::detId, DetId::Ecal, EcalBarrel, EcalEndcap, EcalTriggerTower, elecMap_, Exception, EcalElectronicsMapping::getDCCandSC(), EcalElectronicsMapping::getElectronicsId(), iTtEta2cIndex(), iTtPhi2cIndex(), maxDccId_, minDccId_, nOneEeTtEta, DetId::rawId(), mps_fire::result, and EcalElectronicsId::towerId().

Referenced by anaDigi(), analyzeEB(), analyzeEE(), checkSrApplication(), compareSrfColl(), and dccId().

                                                                              {
   if (detId.det() != DetId::Ecal) {
     throw cms::Exception("InvalidParameter") << "Wrong type of DetId passed to the "
                                                 "EcalSelectiveReadoutValidation::dccCh(const DetId&). "
                                                 "An ECAL DetId was expected.\n";
   }
 
   DetId xtalId;
   switch (detId.subdetId()) {
     case EcalTriggerTower:  //Trigger tower
     {
       const EcalTrigTowerDetId tt = detId;
       //pick up one crystal of the trigger tower: they are however all readout by
       //the same DCC channel in the barrel.
       //Arithmetic is easier on the "c" indices:
       const int iTtPhi0 = iTtPhi2cIndex(tt.iphi());
       const int iTtEta0 = iTtEta2cIndex(tt.ieta());
       const int oneXtalPhi0 = iTtPhi0 * 5;
       const int oneXtalEta0 = (iTtEta0 - nOneEeTtEta) * 5;
 
       xtalId = EBDetId(cIndex2iEta(oneXtalEta0), cIndex2iPhi(oneXtalPhi0));
     } break;
     case EcalEndcap:
       if (detId.rawId() & 0x8000) {  //Supercrystal
         return elecMap_->getDCCandSC(EcalScDetId(detId));
       } else {  //EE crystal
         xtalId = detId;
       }
       break;
     case EcalBarrel:  //EB crystal
       xtalId = detId;
       break;
     default:
       throw cms::Exception("InvalidParameter")
           << "Wrong type of DetId passed to the method "
              "EcalSelectiveReadoutValidation::dccCh(const DetId&). "
              "A valid EcalTriggerTower, EcalBarrel or EcalEndcap DetId was expected. "
              "detid = "
           << xtalId.rawId() << ".\n";
   }
 
   const EcalElectronicsId& EcalElecId = elecMap_->getElectronicsId(xtalId);
 
   pair<int, int> result;
   result.first = EcalElecId.dccId();
 
   if (result.first < minDccId_ || result.second > maxDccId_) {
     throw cms::Exception("OutOfRange") << "Got an invalid DCC ID, DCCID = " << result.first << " for DetId 0x" << hex
                                        << detId.rawId() << " and 0x" << xtalId.rawId() << dec << "\n";
   }
 
   result.second = EcalElecId.towerId();
 
   if (result.second < 1 || result.second > 68) {
     throw cms::Exception("OutOfRange") << "Got an invalid DCC channel ID, DCC_CH = " << result.second << " for DetId 0x"
                                        << hex << detId.rawId() << " and 0x" << xtalId.rawId() << dec << "\n";
   }
 
   return result;
 }

◆ dccId() [1/2]

int EcalSelectiveReadoutValidation::dccId ( const EcalScDetId & detId ) const

private

Retrieves the ID of the DCC reading a readout unit

Parameters

detId detid of the readout unit

Definition at line 2276 of file EcalSelectiveReadoutValidation.cc.

References hcalRecHitTable_cff::detId, elecMap_, and EcalElectronicsMapping::getDCCandSC().

Referenced by selectFedsForLog().

2276 { return elecMap_->getDCCandSC(detId).first; }

EcalElectronicsMapping::getDCCandSC

std::pair< int, int > getDCCandSC(EcalScDetId id) const

Definition: EcalElectronicsMapping.cc:762

hcalRecHitTable_cff.detId

detId

Definition: hcalRecHitTable_cff.py:12

EcalSelectiveReadoutValidation::elecMap_

const EcalElectronicsMapping * elecMap_

Definition: EcalSelectiveReadoutValidation.h:675

◆ dccId() [2/2]

int EcalSelectiveReadoutValidation::dccId ( const EcalTrigTowerDetId & detId ) const

private

Definition at line 2278 of file EcalSelectiveReadoutValidation.cc.

References dccCh(), hcalRecHitTable_cff::detId, and Exception.

                                                                                {
   if (detId.ietaAbs() > 17) {
     throw cms::Exception("InvalidArgument")
         << "Argument of EcalSelectiveReadoutValidation::dccId(const EcalTrigTowerDetId&) "
         << "must be a barrel trigger tower Id\n";
   }
   return dccCh(detId).first;
 }

◆ dccZsFIR()

int EcalSelectiveReadoutValidation::dccZsFIR	(	const EcalDataFrame &	frame,
		const std::vector< int > &	firWeights,
		int	firstFIRSample,
		bool *	saturated = `nullptr`
	)

staticprivate

Emulates the DCC zero suppression FIR filter. If one of the time sample is not in gain 12, numeric_limits<int>::max() is returned.

Parameters

frame	data frame
firWeights	TAP weights
firstFIRSample	index (starting from 1) of the first time sample to be used in the filter
saturated	if not null, *saturated is set to true if all the time sample are not in gain 12 and set to false otherwise.

Returns: FIR output or numeric_limits<int>::max().

Definition at line 2096 of file EcalSelectiveReadoutValidation.cc.

References LogTrace, SiStripPI::max, ecalGpuTask_cfi::sample, and w().

Referenced by analyzeEB(), and analyzeEE().

                                                               {
   const int nFIRTaps = 6;
   //FIR filter weights:
   const vector<int>& w = firWeights;
 
   //accumulator used to compute weighted sum of samples
   int acc = 0;
   bool gain12saturated = false;
   const int gain12 = 0x01;
   const int lastFIRSample = firstFIRSample + nFIRTaps - 1;
   //LogDebug("DccFir") << "DCC FIR operation: ";
   int iWeight = 0;
   for (int iSample = firstFIRSample - 1; iSample < lastFIRSample; ++iSample, ++iWeight) {
     if (iSample >= 0 && iSample < frame.size()) {
       EcalMGPASample sample(frame[iSample]);
       if (sample.gainId() != gain12)
         gain12saturated = true;
       LogTrace("DccFir") << (iSample >= firstFIRSample ? "+" : "") << sample.adc() << "*(" << w[iWeight] << ")";
       acc += sample.adc() * w[iWeight];
     } else {
       edm::LogWarning("DccFir") << __FILE__ << ":" << __LINE__
                                 << ": Not enough samples in data frame or 'ecalDccZs1stSample' module "
                                    "parameter is not valid...";
     }
   }
   LogTrace("DccFir") << "\n";
   //discards the 8 LSBs
   //(shift operator cannot be used on negative numbers because
   // the result depends on compilator implementation)
   acc = (acc >= 0) ? (acc >> 8) : -(-acc >> 8);
   //ZS passed if weighted sum acc above ZS threshold or if
   //one sample has a lower gain than gain 12 (that is gain 12 output
   //is saturated)
 
   LogTrace("DccFir") << "acc: " << acc << "\n"
                      << "saturated: " << (gain12saturated ? "yes" : "no") << "\n";
 
   if (saturated) {
     *saturated = gain12saturated;
   }
 
   return gain12saturated ? numeric_limits<int>::max() : acc;
 }

◆ dqmBeginRun()

void EcalSelectiveReadoutValidation::dqmBeginRun	(	edm::Run const &	r,
		edm::EventSetup const &	es
	)

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 842 of file EcalSelectiveReadoutValidation.cc.

References ecalmapping, elecMap_, edm::EventSetup::getData(), hTriggerTowerMap, initAsciiFile(), and triggerTowerMap_.

                                                                                          {
   // endcap mapping
   triggerTowerMap_ = &es.getData(hTriggerTowerMap);
 
   //electronics map
   elecMap_ = &es.getData(ecalmapping);
 
   initAsciiFile();
 }

◆ dqmEndRun()

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

overrideprotectedvirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 852 of file EcalSelectiveReadoutValidation.cc.

References dqm::impl::MonitorElement::Fill(), getL1aRate(), and meL1aRate_.

                                                                                        {
   meL1aRate_->Fill(getL1aRate());
 }

◆ fill() [1/4]

void EcalSelectiveReadoutValidation::fill	(	MonitorElement *	me,
		float	x
	)

inlineprivate

Wrapper to fill methods of DQM monitor elements.

Definition at line 350 of file EcalSelectiveReadoutValidation.h.

References hlt_dqm_clientPB-live_cfg::me, and x.

Referenced by analyze(), analyzeDataVolume(), analyzeEB(), analyzeEE(), analyzeTP(), checkSrApplication(), and compareSrfColl().

                                          {
     if (me)
       me->Fill(x);
   }

◆ fill() [2/4]

void EcalSelectiveReadoutValidation::fill	(	MonitorElement *	me,
		float	x,
		float	yw
	)

inlineprivate

Definition at line 354 of file EcalSelectiveReadoutValidation.h.

References hlt_dqm_clientPB-live_cfg::me, and x.

                                                    {
     if (me)
       me->Fill(x, yw);
   }

◆ fill() [3/4]

void EcalSelectiveReadoutValidation::fill	(	MonitorElement *	me,
		float	x,
		float	y,
		float	zw
	)

inlineprivate

Definition at line 358 of file EcalSelectiveReadoutValidation.h.

References hlt_dqm_clientPB-live_cfg::me, x, y, and zw().

                                                             {
     if (me)
       me->Fill(x, y, zw);
   }

◆ fill() [4/4]

void EcalSelectiveReadoutValidation::fill	(	MonitorElement *	me,
		float	x,
		float	y,
		float	z,
		float	w
	)

inlineprivate

Definition at line 362 of file EcalSelectiveReadoutValidation.h.

References hlt_dqm_clientPB-live_cfg::me, w(), x, y, and z.

                                                                     {
     if (me)
       me->Fill(x, y, z, w);
   }

◆ frame2Energy()

double EcalSelectiveReadoutValidation::frame2Energy ( const EcalDataFrame & frame ) const

private

Selective Readout decisions Validation

Parameters

event	EDM event
es	event setup Energy reconstruction from ADC samples.
frame	the ADC sample of an ECA channel

Definition at line 1751 of file EcalSelectiveReadoutValidation.cc.

References cms::cuda::assert(), visDQMUpload::buf, mps_fire::i, createfilelist::int, SiStripPI::min, dqmiodumpmetadata::n, and weights_.

Referenced by analyzeEB(), and analyzeEE().

                                                                                     {
   static std::atomic<bool> firstCall{true};
   bool expected = true;
   if (firstCall.compare_exchange_strong(expected, false)) {
     stringstream buf;
     buf << "Weights:";
     for (unsigned i = 0; i < weights_.size(); ++i) {
       buf << "\t" << weights_[i];
     }
     edm::LogInfo("EcalSrValid") << buf.str() << "\n";
     firstCall = false;
   }
   double adc2GeV = 0.;
 
   if (typeid(EBDataFrame) == typeid(frame)) {  //barrel APD
     adc2GeV = .035;
   } else if (typeid(EEDataFrame) == typeid(frame)) {  //endcap VPT
     adc2GeV = 0.06;
   } else {
     assert(false);
   }
 
   double acc = 0;
 
   const int n = min(frame.size(), (int)weights_.size());
 
   double gainInv[] = {12., 1., 6., 12.};
 
   for (int i = 0; i < n; ++i) {
     acc += weights_[i] * frame[i].adc() * gainInv[frame[i].gainId()] * adc2GeV;
   }
   return acc;
 }

◆ frame2EnergyForTp()

template<class T >

double EcalSelectiveReadoutValidation::frame2EnergyForTp	(	const T &	frame,
		int	offset = `0`
	)		const

private

Energy reconstruction from ADC samples to be used for trigger primitive estimate.

Parameters

frame	the ADC sample of an ECA channel
offset	time offset. To be used to evaluate energy of the event previous (offset=-1) and next (offset=+1) to the triggered one.

Definition at line 1938 of file EcalSelectiveReadoutValidation.cc.

References Exception, mps_fire::i, dqmiodumpmetadata::n, hltrates_dqm_sourceclient-live_cfg::offset, and hltDeepSecondaryVertexTagInfosPFPuppi_cfi::weights.

Referenced by setTtEtSums().

                                                                                          {
   //we have to start by 0 in order to handle offset=-1
   //(however Fenix FIR has AFAK only 5 taps)
   double weights[] = {0., -1 / 3., -1 / 3., -1 / 3., 0., 1.};
 
   double adc2GeV = 0.;
   if (typeid(frame) == typeid(EBDataFrame)) {
     adc2GeV = 0.035;
   } else if (typeid(frame) == typeid(EEDataFrame)) {
     adc2GeV = 0.060;
   } else {  //T is an invalid type!
     //TODO: replace message by a cms exception
     throw cms::Exception("Severe Error") << __FILE__ << ":" << __LINE__ << ": "
                                          << "this is a bug. Please report it.\n";
   }
 
   double acc = 0;
 
   const int n = min<int>(frame.size(), sizeof(weights) / sizeof(weights[0]));
 
   double gainInv[] = {12., 1., 6., 12};
 
   for (int i = offset; i < n; ++i) {
     int iframe = i + offset;
     if (iframe >= 0 && iframe < frame.size()) {
       acc += weights[i] * frame[iframe].adc() * gainInv[frame[iframe].gainId()] * adc2GeV;
     }
   }
   //cout << "\n";
   return acc;
 }

◆ getBytesPerCrystal()

double EcalSelectiveReadoutValidation::getBytesPerCrystal ( ) const

inlineprivate

Gets the number of bytes per crystal channel of the event part depending on the number of read crystal channels.

Returns: the number of bytes.

Definition at line 154 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), getDccSrDependentPayload(), getEbEventSize(), and getEeEventSize().

154 { return 3 * 8; }

◆ getCrystalCount()

int EcalSelectiveReadoutValidation::getCrystalCount	(	int	iDcc,
		int	iDccCh
	)

private

Retrieves number of crystal channel read out by a DCC channel

Parameters

iDcc	DCC ID starting from 1
iDccCh	DCC channel starting from 1

Returns: crystal count

Definition at line 2365 of file EcalSelectiveReadoutValidation.cc.

References maxDccId_.

Referenced by analyzeDataVolume(), and checkSrApplication().

                                                                         {
   if (iDcc < minDccId_ || iDcc > maxDccId_) {  //invalid DCC
     return 0;
   } else if (10 <= iDcc && iDcc <= 45) {  //EB
     return 25;
   } else {  //EE
     int iDccPhi;
     if (iDcc < 10)
       iDccPhi = iDcc;
     else
       iDccPhi = iDcc - 45;
     switch (iDccPhi * 100 + iDccCh) {
       case 110:
       case 232:
       case 312:
       case 412:
       case 532:
       case 610:
       case 830:
       case 806:
         //inner partials at 12, 3, and 9 o'clock
         return 20;
       case 134:
       case 634:
       case 827:
       case 803:
         return 10;
       case 330:
       case 430:
         return 20;
       case 203:
       case 503:
       case 721:
       case 921:
         return 21;
       default:
         return 25;
     }
   }
 }

◆ getDccEventSize()

double EcalSelectiveReadoutValidation::getDccEventSize	(	int	iDcc0,
		double	nReadXtals
	)		const

inlineprivate

Gets the size of an DCC event fragment.

Parameters

iDcc0	the DCC logical number starting from 0.
nReadXtals	number of read crystal channels.

Returns: the DCC event fragment size in bytes.

Definition at line 161 of file EcalSelectiveReadoutValidation.h.

References EB, EE, getDccOverhead(), getDccSrDependentPayload(), and getRuCount().

Referenced by analyzeDataVolume().

                                                              {
     subdet_t subdet;
     if (iDcc0 < 9 || iDcc0 >= 45) {
       subdet = EE;
     } else {
       subdet = EB;
     }
     //     return getDccOverhead(subdet)+nReadXtals*getBytesPerCrystal()
     //       + getRuCount(iDcc0)*8;
     return getDccOverhead(subdet) + getDccSrDependentPayload(iDcc0, getRuCount(iDcc0), nReadXtals);
   }

◆ getDccOverhead()

double EcalSelectiveReadoutValidation::getDccOverhead ( subdet_t subdet ) const

inlineprivate

Gets the size in bytes fixed-size part of a DCC event fragment.

Returns: the fixed size in bytes.

Definition at line 145 of file EcalSelectiveReadoutValidation.h.

References EB.

Referenced by getDccEventSize(), getEbEventSize(), and getEeEventSize().

                                                {
     //  return (subdet==EB?34:25)*8;
     return (subdet == EB ? 34 : 52) * 8;
   }

◆ getDccSrDependentPayload()

double EcalSelectiveReadoutValidation::getDccSrDependentPayload	(	int	iDcc0,
		double	nReadRus,
		double	nReadXtals
	)		const

inlineprivate

Gets DCC event fragment payload depending on the channel selection made by the selective readout.

Parameters

iDcc0	the DCC logical number starting from 0.
nReadRus	number of read-out RUs
nReadXtals	number of read-out crystal channels.

Returns: the DCC event fragment payload in bytes.

Definition at line 180 of file EcalSelectiveReadoutValidation.h.

References getBytesPerCrystal().

Referenced by analyzeDataVolume(), and getDccEventSize().

                                                                                        {
     return nReadXtals * getBytesPerCrystal() + nReadRus * 8;
   }

◆ getEbEventSize()

double EcalSelectiveReadoutValidation::getEbEventSize ( double nReadXtals ) const

private

Computes the size of an ECAL barrel event fragment.

Parameters

nReadXtals number of read crystal channels

Returns: the event fragment size in bytes

Definition at line 2072 of file EcalSelectiveReadoutValidation.cc.

References EB, getBytesPerCrystal(), getDccOverhead(), getRuCount(), nEbDccs, and nEeDccs.

Referenced by analyzeDataVolume().

                                                                              {
   double ruHeaderPayload = 0.;
   const int firstEbDcc0 = nEeDccs / 2;
   for (int iDcc0 = firstEbDcc0; iDcc0 < firstEbDcc0 + nEbDccs; ++iDcc0) {
     ruHeaderPayload += getRuCount(iDcc0) * 8.;
   }
 
   return getDccOverhead(EB) * nEbDccs + nReadXtals * getBytesPerCrystal() + ruHeaderPayload;
 }

◆ getEeEventSize()

double EcalSelectiveReadoutValidation::getEeEventSize ( double nReadXtals ) const

private

Computes the size of an ECAL endcap event fragment.

Parameters

nReadXtals number of read crystal channels

Returns: the event fragment size in bytes

Definition at line 2082 of file EcalSelectiveReadoutValidation.cc.

References EE, getBytesPerCrystal(), getDccOverhead(), getRuCount(), nDccs_, nEbDccs, and nEeDccs.

Referenced by analyzeDataVolume().

                                                                              {
   double ruHeaderPayload = 0.;
   const unsigned firstEbDcc0 = nEeDccs / 2;
   for (unsigned iDcc0 = 0; iDcc0 < nDccs_; ++iDcc0) {
     //skip barrel:
     if (iDcc0 == firstEbDcc0)
       iDcc0 += nEbDccs;
     ruHeaderPayload += getRuCount(iDcc0) * 8.;
   }
   return getDccOverhead(EE) * nEeDccs + nReadXtals * getBytesPerCrystal() + ruHeaderPayload;
 }

◆ getFIRWeights()

std::vector< int > EcalSelectiveReadoutValidation::getFIRWeights ( const std::vector< double > & normalizedWeights )

staticprivate

Computes the ZS FIR filter weights from the normalized weights.

Parameters

normalizedWeights the normalized weights

Returns: the computed ZS filter weights.

Definition at line 2143 of file EcalSelectiveReadoutValidation.cc.

References funct::abs(), mps_fire::i, and SiStripPI::min.

Referenced by configFirWeights().

                                                                                                        {
   const int nFIRTaps = 6;
   vector<int> firWeights(nFIRTaps, 0);  //default weight: 0;
   const static int maxWeight = 0xEFF;   //weights coded on 11+1 signed bits
   for (unsigned i = 0; i < min((size_t)nFIRTaps, normalizedWeights.size()); ++i) {
     firWeights[i] = lround(normalizedWeights[i] * (1 << 10));
     if (abs(firWeights[i]) > maxWeight) {  //overflow
       firWeights[i] = firWeights[i] < 0 ? -maxWeight : maxWeight;
     }
   }
   return firWeights;
 }

◆ getL1aRate()

double EcalSelectiveReadoutValidation::getL1aRate ( ) const

private

Gets L1A rate estimate.

See also: updateL1aRate(const edm::Event&)

Returns: L1A rate estimate

Definition at line 211 of file EcalSelectiveReadoutValidation.cc.

References MillePedeFileConverter_cfg::e, l1aOfTmax, l1aOfTmin, LogDebug, tmax, and tmin.

Referenced by dqmEndRun().

                                                         {
   LogDebug("EcalSrValid") << __FILE__ << ":" << __LINE__ << ": "
                           << "Tmax = " << tmax << " x 25ns; Tmin = " << tmin << " x 25ns; L1A(Tmax) = " << l1aOfTmax
                           << "; L1A(Tmin) = " << l1aOfTmin << "\n";
   return (double)(l1aOfTmax - l1aOfTmin) / ((tmax - tmin) * 25e-9);
 }

◆ getRuCount()

int EcalSelectiveReadoutValidation::getRuCount ( int iDcc0 ) const

private

Gets the number of readout unit read by a DCC. A readout unit correspond to an active DCC input channel.

Parameters

iDcc0 DCC logical number starting from 0.

Definition at line 1785 of file EcalSelectiveReadoutValidation.cc.

References nRuPerDcc_.

Referenced by getDccEventSize(), getEbEventSize(), and getEeEventSize().

1785 { return nRuPerDcc_[iDcc0]; }

EcalSelectiveReadoutValidation::nRuPerDcc_

int nRuPerDcc_[nDccs_]

Definition: EcalSelectiveReadoutValidation.h:756

◆ iEta2cIndex()

int EcalSelectiveReadoutValidation::iEta2cIndex ( int iEta ) const

inlineprivate

Converts a std CMSSW crystal eta index to a c-array index (starting from zero and without hole).

Definition at line 224 of file EcalSelectiveReadoutValidation.h.

References l1tTowerCalibrationProducer_cfi::iEta.

Referenced by anaDigi(), and analyzeEB().

224 { return (iEta < 0) ? iEta + 85 : iEta + 84; }

l1tTowerCalibrationProducer_cfi.iEta

iEta

Definition: l1tTowerCalibrationProducer_cfi.py:60

◆ initAsciiFile()

void EcalSelectiveReadoutValidation::initAsciiFile ( )

private

Definition at line 2202 of file EcalSelectiveReadoutValidation.cc.

References Exception, logSrApplicationErrors_, logSrpAlgoErrors_, MillePedeFileConverter_cfg::out, srApplicationErrorLog_, srApplicationErrorLogFileName_, srpAlgoErrorLog_, srpAlgoErrorLogFileName_, and pileupReCalc_HLTpaths::trunc.

Referenced by dqmBeginRun().

                                                    {
   if (logSrpAlgoErrors_) {
     srpAlgoErrorLog_.open(srpAlgoErrorLogFileName_.c_str(), ios::out | ios::trunc);
     if (!srpAlgoErrorLog_.good()) {
       throw cms::Exception("Output") << "Failed to open the log file '" << srpAlgoErrorLogFileName_
                                      << "' for SRP algorithm result check.\n";
     }
   }
 
   if (logSrApplicationErrors_) {
     srApplicationErrorLog_.open(srApplicationErrorLogFileName_.c_str(), ios::out | ios::trunc);
     if (!srApplicationErrorLog_.good()) {
       throw cms::Exception("Output") << "Failed to open the log file '" << srApplicationErrorLogFileName_
                                      << "' for Selective Readout decision application check.\n";
     }
   }
 }

◆ iPhi2cIndex()

int EcalSelectiveReadoutValidation::iPhi2cIndex ( int iPhi ) const

inlineprivate

Converts a std CMSSW crystal phi index to a c-array index (starting from zero and without hole).

Definition at line 229 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), and analyzeEB().

                                   {
     //    return iPhi-1;
     int iPhi0 = iPhi - 11;
     if (iPhi0 < 0)
       iPhi0 += 360;
     return iPhi0;
   }

◆ iTtEta2cIndex()

int EcalSelectiveReadoutValidation::iTtEta2cIndex ( int iEta ) const

inlineprivate

Transforms CMSSW eta ECAL TT indices to indices starting at 0 to use for c-array or vector.

Parameters

iEta	CMSSW eta index (numbering -28...-1,28...56)

Returns: index in numbering from 0 to 55

Definition at line 259 of file EcalSelectiveReadoutValidation.h.

References l1tTowerCalibrationProducer_cfi::iEta.

Referenced by analyzeTP(), dccCh(), and setTtEtSums().

259 { return (iEta < 0) ? iEta + 28 : iEta + 27; }

l1tTowerCalibrationProducer_cfi.iEta

iEta

Definition: l1tTowerCalibrationProducer_cfi.py:60

◆ iTtPhi2cIndex()

int EcalSelectiveReadoutValidation::iTtPhi2cIndex ( int iPhi ) const

inlineprivate

Transforms CMSSW phi ECAL crystal indices to indices starting at 0 to use for c-array or vector.

Parameters

iPhi	CMSSW phi index (numbering 1...72)

Returns: index in numbering 0...71

Definition at line 266 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), dccCh(), and setTtEtSums().

                                     {
     return iPhi - 1;
     //int iPhi0 = iPhi - 3;
     //if(iPhi0<0) iPhi0 += 72;
     //return iPhi0;
   }

◆ iXY2cIndex()

int EcalSelectiveReadoutValidation::iXY2cIndex ( int iX ) const

inlineprivate

Converts a std CMSSW crystal x or y index to a c-array index (starting from zero and without hole).

Definition at line 240 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), and analyzeEE().

240 { return iX - 1; }

◆ printAvailableHists()

void EcalSelectiveReadoutValidation::printAvailableHists ( )

private

Prints the list of available histograms (registered by the registerHist method), including disabled one.

Definition at line 2062 of file EcalSelectiveReadoutValidation.cc.

References availableHistList_, and dqm-mbProfile::log.

Referenced by bookHistograms().

                                                          {
   LogInfo log("HistoList");
   log << "Avalailable histograms (DQM monitor elements): \n";
   for (map<string, string>::iterator it = availableHistList_.begin(); it != availableHistList_.end(); ++it) {
     log << it->first << ": " << it->second << "\n";
   }
   log << "\nTo include an histogram add its name in the vstring parameter "
          "'histograms' of the EcalSelectiveReadoutValidation module\n";
 }

◆ readAllCollections()

void EcalSelectiveReadoutValidation::readAllCollections ( const edm::Event & e )

private

Reads the data collections from the event. Called at start of each event analysis.

Parameters

event the EDM event.

Definition at line 2045 of file EcalSelectiveReadoutValidation.cc.

References ebComputedSrFlags_, ebDigis_, ebNoZsDigis_, ebRecHits_, ebSimHits_, ebSrFlags_, eeComputedSrFlags_, eeDigis_, eeNoZsDigis_, eeRecHits_, eeSimHits_, eeSrFlags_, fedRaw_, CollHandle< T >::read(), and tps_.

Referenced by analyze().

                                                                              {
   ebRecHits_.read(event);
   eeRecHits_.read(event);
   ebDigis_.read(event);
   eeDigis_.read(event);
   ebNoZsDigis_.read(event);
   eeNoZsDigis_.read(event);
   ebSrFlags_.read(event);
   eeSrFlags_.read(event);
   ebComputedSrFlags_.read(event);
   eeComputedSrFlags_.read(event);
   ebSimHits_.read(event);
   eeSimHits_.read(event);
   tps_.read(event);
   fedRaw_.read(event);
 }

◆ readOutUnitOf() [1/2]

EcalTrigTowerDetId EcalSelectiveReadoutValidation::readOutUnitOf ( const EBDetId & xtalId ) const

private

Retrives the readout unit, a trigger tower in the barrel case, and a supercrystal in the endcap case, a given crystal belongs to.

Parameters

xtalId identifier of the crystal

Returns: identifer of the supercrystal or of the trigger tower.

Definition at line 1848 of file EcalSelectiveReadoutValidation.cc.

References EcalTrigTowerConstituentsMap::towerOf(), and triggerTowerMap_.

Referenced by anaDigi(), analyzeEB(), and analyzeEE().

                                                                                             {
   return triggerTowerMap_->towerOf(xtalId);
 }

◆ readOutUnitOf() [2/2]

EcalScDetId EcalSelectiveReadoutValidation::readOutUnitOf ( const EEDetId & xtalId ) const

private

Definition at line 1852 of file EcalSelectiveReadoutValidation.cc.

References EcalElectronicsId::dccId(), elecMap_, EcalElectronicsMapping::getEcalScDetId(), EcalElectronicsMapping::getElectronicsId(), and EcalElectronicsId::towerId().

                                                                                      {
   const EcalElectronicsId& EcalElecId = elecMap_->getElectronicsId(xtalId);
   int iDCC = EcalElecId.dccId();
   int iDccChan = EcalElecId.towerId();
   const bool ignoreSingle = true;
   const vector<EcalScDetId> id = elecMap_->getEcalScDetId(iDCC, iDccChan, ignoreSingle);
   return !id.empty() ? id[0] : EcalScDetId();
 }

◆ registerHist()

bool EcalSelectiveReadoutValidation::registerHist	(	const std::string &	name,
		const std::string &	title
	)

private

Register a histogram in the available histogram list and check if the histogram is enabled. Called by the histogram booking methods.

Returns: true if the histogram is enable, false otherwise

Definition at line 2040 of file EcalSelectiveReadoutValidation.cc.

References allHists_, availableHistList_, histList_, Skims_PA_cff::name, and runGCPTkAlMap::title.

Referenced by book1D(), book2D(), bookFloat(), bookProfile(), and bookProfile2D().

                                                                                                {
   availableHistList_.insert(pair<string, string>(name, title));
   return allHists_ || histList_.find(name) != histList_.end();
 }

◆ ruGraphX() [1/2]

int EcalSelectiveReadoutValidation::ruGraphX ( const EcalScDetId & id ) const

inlineprivate

Functions to compute x and y coordinates of RU maps grouping endcap and barrel.

Definition at line 870 of file EcalSelectiveReadoutValidation.h.

Referenced by checkSrApplication(), and compareSrfColl().

870 { return id.ix() + (id.zside() > 0 ? 20 : -40); }

◆ ruGraphX() [2/2]

int EcalSelectiveReadoutValidation::ruGraphX ( const EcalTrigTowerDetId & id ) const

inlineprivate

Definition at line 874 of file EcalSelectiveReadoutValidation.h.

874 { return id.ieta(); }

◆ ruGraphY() [1/2]

int EcalSelectiveReadoutValidation::ruGraphY ( const EcalScDetId & id ) const

inlineprivate

Definition at line 872 of file EcalSelectiveReadoutValidation.h.

Referenced by checkSrApplication(), and compareSrfColl().

872 { return id.iy(); }

◆ ruGraphY() [2/2]

int EcalSelectiveReadoutValidation::ruGraphY ( const EcalTrigTowerDetId & id ) const

inlineprivate

Definition at line 876 of file EcalSelectiveReadoutValidation.h.

876 { return id.iphi(); }

◆ selectFedsForLog()

void EcalSelectiveReadoutValidation::selectFedsForLog ( )

private

Look in events whose DCC has SR flags and enable error logging for them. To be called with the processed first event. List of monitored DCCs is reported in the log file.

Definition at line 2287 of file EcalSelectiveReadoutValidation.cc.

References visDQMUpload::buf, dccId(), ebSrFlags_, eeSrFlags_, dqmdumpme::first, logErrForDccs_, minDccId_, nDccs_, srApplicationErrorLog_, and srpAlgoErrorLog_.

Referenced by analyze().

                                                       {
   logErrForDccs_ = vector<bool>(nDccs_, false);
 
   for (EBSrFlagCollection::const_iterator it = ebSrFlags_->begin(); it != ebSrFlags_->end(); ++it) {
     int iDcc = dccId(it->id()) - minDccId_;
 
     logErrForDccs_.at(iDcc) = true;
   }
 
   for (EESrFlagCollection::const_iterator it = eeSrFlags_->begin(); it != eeSrFlags_->end(); ++it) {
     int iDcc = dccId(it->id()) - minDccId_;
 
     logErrForDccs_.at(iDcc) = true;
   }
 
   stringstream buf;
   buf << "List of DCCs found in the first processed event: ";
   bool first = true;
   for (unsigned iDcc = 0; iDcc < nDccs_; ++iDcc) {
     if (logErrForDccs_[iDcc]) {
       buf << (first ? "" : ", ") << (iDcc + minDccId_);
       first = false;
     }
   }
   buf << "\nOnly DCCs from this list will be considered for error logging\n";
   srpAlgoErrorLog_ << buf.str();
   srApplicationErrorLog_ << buf.str();
   LogInfo("EcalSrValid") << buf.str();
 }

◆ setTtEtSums()

void EcalSelectiveReadoutValidation::setTtEtSums	(	const edm::EventSetup &	es,
		const EBDigiCollection &	ebDigis,
		const EEDigiCollection &	eeDigis
	)

private

Computes trigger primitive estimates. A sum of crystal deposited transverse energy is performed.

Parameters

es	event setup
ebDigis	the ECAL barrel unsuppressed digi to use for the computation
ebDigis	the ECAL endcap unsuppressed digi to use for the computation

Definition at line 1861 of file EcalSelectiveReadoutValidation.cc.

References edm::DataFrameContainer::begin(), MillePedeFileConverter_cfg::e, ebDigis_, DetId::Ecal, EcalBarrel, EcalEndcap, SimL1EmulatorDM_cff::eeDigis, edm::DataFrameContainer::end(), EgHLTOffHistBins_cfi::et, frame2EnergyForTp(), geoToken, edm::EventSetup::getData(), CaloSubdetectorGeometry::getGeometry(), edm::EventSetup::getHandle(), iTtEta2cIndex(), iTtPhi2cIndex(), nTtEta, nTtPhi, funct::sin(), theta(), EcalTrigTowerConstituentsMap::towerOf(), triggerTowerMap_, ttEtSums, and ecaldqm::ttId().

Referenced by analyze().

                                                                                   {
   //ecal geometry:
   const CaloSubdetectorGeometry* eeGeometry = nullptr;
   const CaloSubdetectorGeometry* ebGeometry = nullptr;
   if (eeGeometry == nullptr || ebGeometry == nullptr) {
     es.getData(geoToken);
     auto geoHandle = es.getHandle(geoToken);
     eeGeometry = (*geoHandle).getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
     ebGeometry = (*geoHandle).getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
   }
 
   //init etSum array:
   for (int iEta0 = 0; iEta0 < nTtEta; ++iEta0) {
     for (int iPhi0 = 0; iPhi0 < nTtPhi; ++iPhi0) {
       ttEtSums[iEta0][iPhi0] = 0.;
     }
   }
 
   for (EBDigiCollection::const_iterator it = ebDigis_->begin(); it != ebDigis_->end(); ++it) {
     const EBDataFrame& frame = *it;
     const EcalTrigTowerDetId& ttId = triggerTowerMap_->towerOf(frame.id());
 
     const int iTtEta0 = iTtEta2cIndex(ttId.ieta());
     const int iTtPhi0 = iTtPhi2cIndex(ttId.iphi());
     double theta = ebGeometry->getGeometry(frame.id())->getPosition().theta();
     double e = frame2EnergyForTp(frame);
     if ((frame2EnergyForTp(frame, -1) < e) && (frame2EnergyForTp(frame, 1) < e)) {
       ttEtSums[iTtEta0][iTtPhi0] += e * sin(theta);
     }
   }
 
   for (EEDigiCollection::const_iterator it = eeDigis.begin(); it != eeDigis.end(); ++it) {
     const EEDataFrame& frame = *it;
     const EcalTrigTowerDetId& ttId = triggerTowerMap_->towerOf(frame.id());
     const int iTtEta0 = iTtEta2cIndex(ttId.ieta());
     const int iTtPhi0 = iTtPhi2cIndex(ttId.iphi());
 
     double theta = eeGeometry->getGeometry(frame.id())->getPosition().theta();
     double e = frame2EnergyForTp(frame);
     if ((frame2EnergyForTp(frame, -1) < e) && (frame2EnergyForTp(frame, 1) < e)) {
       ttEtSums[iTtEta0][iTtPhi0] += e * sin(theta);
     }
   }
 
   //dealing with pseudo-TT in two inner EE eta-ring:
   int innerTTEtas[] = {0, 1, 54, 55};
   for (unsigned iRing = 0; iRing < sizeof(innerTTEtas) / sizeof(innerTTEtas[0]); ++iRing) {
     int iTtEta0 = innerTTEtas[iRing];
     //this detector eta-section is divided in only 36 phi bins
     //For this eta regions,
     //current tower eta numbering scheme is inconsistent. For geometry
     //version 133:
     //- TT are numbered from 0 to 72 for 36 bins
     //- some TT have an even index, some an odd index
     //For geometry version 125, there are 72 phi bins.
     //The code below should handle both geometry definition.
     //If there are 72 input trigger primitives for each inner eta-ring,
     //then the average of the trigger primitive of the two pseudo-TT of
     //a pair (nEta, nEta+1) is taken as Et of both pseudo TTs.
     //If there are only 36 input TTs for each inner eta ring, then half
     //of the present primitive of a pseudo TT pair is used as Et of both
     //pseudo TTs.
 
     for (unsigned iTtPhi0 = 0; iTtPhi0 < nTtPhi - 1; iTtPhi0 += 2) {
       double et = .5 * (ttEtSums[iTtEta0][iTtPhi0] + ttEtSums[iTtEta0][iTtPhi0 + 1]);
       //divides the TT into 2 phi bins in order to match with 72 phi-bins SRP
       //scheme or average the Et on the two pseudo TTs if the TT is already
       //divided into two trigger primitives.
       ttEtSums[iTtEta0][iTtPhi0] = et;
       ttEtSums[iTtEta0][iTtPhi0 + 1] = et;
     }
   }
 }

◆ updateL1aRate()

void EcalSelectiveReadoutValidation::updateL1aRate ( const edm::Event & event )

private

Updates estimate of L1A rate

Parameters

event EDM event

Definition at line 193 of file EcalSelectiveReadoutValidation.cc.

References nano_mu_digi_cff::bx, l1aOfTmax, l1aOfTmin, submitPVValidationJobs::t, tmax, and tmin.

Referenced by analyze().

                                                                         {
   const int32_t bx = event.bunchCrossing();
   if (bx < 1 || bx > 3564)
     return;
 
   int64_t t = event.bunchCrossing() + (event.orbitNumber() - 1) * 3564;
 
   if (t < tmin) {
     tmin = t;
     l1aOfTmin = event.id().event();
   }
 
   if (t > tmax) {
     tmax = t;
     l1aOfTmax = event.id().event();
   }
 }

◆ xtalGraphX() [1/2]

int EcalSelectiveReadoutValidation::xtalGraphX ( const EEDetId & id ) const

inlineprivate

Definition at line 878 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and analyzeEE().

878 { return id.ix() + (id.zside() > 0 ? 100 : -200); }

◆ xtalGraphX() [2/2]

int EcalSelectiveReadoutValidation::xtalGraphX ( const EBDetId & id ) const

inlineprivate

Definition at line 882 of file EcalSelectiveReadoutValidation.h.

882 { return id.ieta(); }

◆ xtalGraphY() [1/2]

int EcalSelectiveReadoutValidation::xtalGraphY ( const EEDetId & id ) const

inlineprivate

Definition at line 880 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and analyzeEE().

880 { return id.iy(); }

◆ xtalGraphY() [2/2]

int EcalSelectiveReadoutValidation::xtalGraphY ( const EBDetId & id ) const

inlineprivate

Definition at line 884 of file EcalSelectiveReadoutValidation.h.

884 { return id.iphi(); }

Member Data Documentation

◆ allHists_

bool EcalSelectiveReadoutValidation::allHists_

private

When true, every histogram is enabled.

Definition at line 800 of file EcalSelectiveReadoutValidation.h.

Referenced by EcalSelectiveReadoutValidation(), and registerHist().

◆ availableHistList_

std::map<std::string, std::string> EcalSelectiveReadoutValidation::availableHistList_

private

List of available histograms. Filled by the booking methods. key: name, value: title.

Definition at line 809 of file EcalSelectiveReadoutValidation.h.

Referenced by bookHistograms(), printAvailableHists(), and registerHist().

◆ collNotFoundWarn_

bool EcalSelectiveReadoutValidation::collNotFoundWarn_

private

Switch for collection-not-found warning.

Definition at line 491 of file EcalSelectiveReadoutValidation.h.

◆ ebComputedSrFlags_

CollHandle<EBSrFlagCollection> EcalSelectiveReadoutValidation::ebComputedSrFlags_

private

Definition at line 519 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), EcalSelectiveReadoutValidation(), and readAllCollections().

◆ ebDigis_

CollHandle<EBDigiCollection> EcalSelectiveReadoutValidation::ebDigis_

private

The event product collections.

Definition at line 513 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), analyzeEB(), EcalSelectiveReadoutValidation(), readAllCollections(), and setTtEtSums().

◆ ebNoZsDigis_

CollHandle<EBDigiCollection> EcalSelectiveReadoutValidation::ebNoZsDigis_

private

Definition at line 515 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), analyzeEB(), EcalSelectiveReadoutValidation(), and readAllCollections().

◆ ebRecHits_

CollHandle<RecHitCollection> EcalSelectiveReadoutValidation::ebRecHits_

private

Definition at line 524 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), EcalSelectiveReadoutValidation(), and readAllCollections().

◆ ebRuActive_

bool EcalSelectiveReadoutValidation::ebRuActive_[nEbEta/ebTtEdge][nEbPhi/ebTtEdge]

private

For book keeping of RU actually read out (not fully zero suppressed)

Definition at line 769 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), and anaDigiInit().

◆ ebSimHits_

CollHandle<std::vector<PCaloHit> > EcalSelectiveReadoutValidation::ebSimHits_

private

Definition at line 521 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), analyzeEB(), EcalSelectiveReadoutValidation(), and readAllCollections().

◆ ebSrFlags_

CollHandle<EBSrFlagCollection> EcalSelectiveReadoutValidation::ebSrFlags_

private

Definition at line 517 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), analyzeDataVolume(), analyzeEB(), EcalSelectiveReadoutValidation(), readAllCollections(), and selectFedsForLog().

◆ ebTtEdge

const int EcalSelectiveReadoutValidation::ebTtEdge = 5

staticprivate

Number of crystals along an EB TT.

Definition at line 462 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi().

◆ ebZsThr_

int EcalSelectiveReadoutValidation::ebZsThr_

private

ZS threshold in 1/4th ADC count for EB

Definition at line 692 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), bookHistograms(), and EcalSelectiveReadoutValidation().

◆ ecalmapping

edm::ESGetToken<EcalElectronicsMapping, EcalMappingRcd> EcalSelectiveReadoutValidation::ecalmapping

private

Definition at line 409 of file EcalSelectiveReadoutValidation.h.

Referenced by dqmBeginRun().

◆ eeComputedSrFlags_

CollHandle<EESrFlagCollection> EcalSelectiveReadoutValidation::eeComputedSrFlags_

private

Definition at line 520 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), EcalSelectiveReadoutValidation(), and readAllCollections().

◆ eeDigis_

CollHandle<EEDigiCollection> EcalSelectiveReadoutValidation::eeDigis_

private

Definition at line 514 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), analyzeEE(), EcalSelectiveReadoutValidation(), and readAllCollections().

◆ eeNoZsDigis_

CollHandle<EEDigiCollection> EcalSelectiveReadoutValidation::eeNoZsDigis_

private

Definition at line 516 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), analyzeEE(), EcalSelectiveReadoutValidation(), and readAllCollections().

◆ eeRecHits_

CollHandle<RecHitCollection> EcalSelectiveReadoutValidation::eeRecHits_

private

Definition at line 525 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE(), EcalSelectiveReadoutValidation(), and readAllCollections().

◆ eeRuActive_

bool EcalSelectiveReadoutValidation::eeRuActive_[nEndcaps][nEeX/scEdge][nEeY/scEdge]

private

Definition at line 770 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), and anaDigiInit().

◆ eeSimHits_

CollHandle<std::vector<PCaloHit> > EcalSelectiveReadoutValidation::eeSimHits_

private

Definition at line 522 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), analyzeEE(), EcalSelectiveReadoutValidation(), and readAllCollections().

◆ eeSrFlags_

CollHandle<EESrFlagCollection> EcalSelectiveReadoutValidation::eeSrFlags_

private

Definition at line 518 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), analyzeDataVolume(), analyzeEE(), EcalSelectiveReadoutValidation(), readAllCollections(), and selectFedsForLog().

◆ eeZsThr_

int EcalSelectiveReadoutValidation::eeZsThr_

private

ZS threshold in 1/4th ADC count for EE

Definition at line 696 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE(), bookHistograms(), and EcalSelectiveReadoutValidation().

◆ elecMap_

const EcalElectronicsMapping* EcalSelectiveReadoutValidation::elecMap_

private

Ecal electronics/geometrical mapping.

Definition at line 675 of file EcalSelectiveReadoutValidation.h.

Referenced by dccCh(), dccId(), dqmBeginRun(), and readOutUnitOf().

◆ fedRaw_

CollHandle<FEDRawDataCollection> EcalSelectiveReadoutValidation::fedRaw_

private

Definition at line 526 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), EcalSelectiveReadoutValidation(), and readAllCollections().

◆ firstFIRSample_

int EcalSelectiveReadoutValidation::firstFIRSample_

private

Time position of the first sample to use in zero suppession FIR filter. Numbering starts at 0.

Definition at line 705 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), analyzeEE(), and configFirWeights().

◆ firWeights_

std::vector<int> EcalSelectiveReadoutValidation::firWeights_

private

Weights to be used for the ZS FIR filter

Definition at line 688 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), analyzeEE(), and configFirWeights().

◆ geoToken

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> EcalSelectiveReadoutValidation::geoToken

private

Definition at line 408 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), analyzeEE(), and setTtEtSums().

◆ histDir_

std::string EcalSelectiveReadoutValidation::histDir_

private

Histogram directory PATH in DQM or within the output ROOT file

Definition at line 804 of file EcalSelectiveReadoutValidation.h.

◆ histList_

std::set<std::string> EcalSelectiveReadoutValidation::histList_

private

List of enabled histograms. Special name "all" is used to indicate all available histograms.

Definition at line 796 of file EcalSelectiveReadoutValidation.h.

Referenced by bookHistograms(), EcalSelectiveReadoutValidation(), and registerHist().

◆ hTriggerTowerMap

edm::ESGetToken<EcalTrigTowerConstituentsMap, IdealGeometryRecord> EcalSelectiveReadoutValidation::hTriggerTowerMap

private

Definition at line 410 of file EcalSelectiveReadoutValidation.h.

Referenced by dqmBeginRun().

◆ ievt_

int EcalSelectiveReadoutValidation::ievt_

private

Event sequence number

Definition at line 781 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), analyzeEB(), and analyzeEE().

◆ isRuComplete_

bool EcalSelectiveReadoutValidation::isRuComplete_[nDccs_][nDccChs_]

private

Definition at line 773 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), analyzeEB(), and analyzeEE().

◆ kByte_

const int EcalSelectiveReadoutValidation::kByte_ = 1024

staticprivate

number of bytes in 1 kByte:

Definition at line 416 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume().

◆ l1aOfTmax

int64_t EcalSelectiveReadoutValidation::l1aOfTmax

private

Definition at line 535 of file EcalSelectiveReadoutValidation.h.

Referenced by getL1aRate(), and updateL1aRate().

◆ l1aOfTmin

int64_t EcalSelectiveReadoutValidation::l1aOfTmin

private

Definition at line 534 of file EcalSelectiveReadoutValidation.h.

Referenced by getL1aRate(), and updateL1aRate().

◆ l1aRateErr

bool EcalSelectiveReadoutValidation::l1aRateErr

private

Definition at line 536 of file EcalSelectiveReadoutValidation.h.

◆ localReco_

bool EcalSelectiveReadoutValidation::localReco_

private

Local reconstruction switch: true to reconstruct locally the amplitude insted of using the Rec Hits.

Definition at line 680 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and analyzeEE().

◆ logErrForDccs_

std::vector<bool> EcalSelectiveReadoutValidation::logErrForDccs_

private

List of TCC masks for validation If tccMasks[iTcc-1] is false then TCC is considered to have been out of the run and related validations are skipped.

Definition at line 716 of file EcalSelectiveReadoutValidation.h.

Referenced by compareSrfColl(), and selectFedsForLog().

◆ logSrApplicationErrors_

bool EcalSelectiveReadoutValidation::logSrApplicationErrors_

private

Switch to log SR decision that fails to be applied on data: inconstitencies between SRF and number of read out crystals.

Definition at line 846 of file EcalSelectiveReadoutValidation.h.

Referenced by EcalSelectiveReadoutValidation(), and initAsciiFile().

◆ logSrpAlgoErrors_

bool EcalSelectiveReadoutValidation::logSrpAlgoErrors_

private

Switch to log in an ascii file inconsistencies found between SRFs read from data and SRFs obtained by rerunning SRP algorithm on TTFs.

Definition at line 841 of file EcalSelectiveReadoutValidation.h.

Referenced by EcalSelectiveReadoutValidation(), and initAsciiFile().

◆ lutGrpHandle

edm::ESGetToken<EcalTPGLutGroup, EcalTPGLutGroupRcd> EcalSelectiveReadoutValidation::lutGrpHandle

private

Definition at line 412 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP().

◆ lutMapHandle

edm::ESGetToken<EcalTPGLutIdMap, EcalTPGLutIdMapRcd> EcalSelectiveReadoutValidation::lutMapHandle

private

Definition at line 413 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP().

◆ maxDccId_

const int EcalSelectiveReadoutValidation::maxDccId_ = minDccId_ + nDccs_ - 1

staticprivate

Definition at line 429 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), dccCh(), and getCrystalCount().

◆ meChOcc_

MonitorElement* EcalSelectiveReadoutValidation::meChOcc_

private

Definition at line 555 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), analyzeEE(), and bookHistograms().

◆ meCompleteZS_

MonitorElement* EcalSelectiveReadoutValidation::meCompleteZS_

private

Definition at line 606 of file EcalSelectiveReadoutValidation.h.

◆ meCompleteZSCnt_

MonitorElement* EcalSelectiveReadoutValidation::meCompleteZSCnt_

private

Definition at line 618 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and bookHistograms().

◆ meCompleteZSMap_

MonitorElement* EcalSelectiveReadoutValidation::meCompleteZSMap_

private

Definition at line 610 of file EcalSelectiveReadoutValidation.h.

Referenced by bookHistograms(), and checkSrApplication().

◆ meCompleteZSRateMap_

MonitorElement* EcalSelectiveReadoutValidation::meCompleteZSRateMap_

private

Definition at line 614 of file EcalSelectiveReadoutValidation.h.

Referenced by bookHistograms(), and checkSrApplication().

◆ meDccHiVol_

MonitorElement* EcalSelectiveReadoutValidation::meDccHiVol_

private

Definition at line 544 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meDccLiVol_

MonitorElement* EcalSelectiveReadoutValidation::meDccLiVol_

private

Definition at line 543 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meDccVol_

MonitorElement* EcalSelectiveReadoutValidation::meDccVol_

private

The histograms

Definition at line 542 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meDccVolFromData_

MonitorElement* EcalSelectiveReadoutValidation::meDccVolFromData_

private

Definition at line 545 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meDroppedFRO_

MonitorElement* EcalSelectiveReadoutValidation::meDroppedFRO_

private

Definition at line 605 of file EcalSelectiveReadoutValidation.h.

◆ meDroppedFROCnt_

MonitorElement* EcalSelectiveReadoutValidation::meDroppedFROCnt_

private

Definition at line 617 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and bookHistograms().

◆ meDroppedFROMap_

MonitorElement* EcalSelectiveReadoutValidation::meDroppedFROMap_

private

Definition at line 609 of file EcalSelectiveReadoutValidation.h.

Referenced by bookHistograms(), and checkSrApplication().

◆ meDroppedFRORateMap_

MonitorElement* EcalSelectiveReadoutValidation::meDroppedFRORateMap_

private

Definition at line 613 of file EcalSelectiveReadoutValidation.h.

Referenced by bookHistograms(), and checkSrApplication().

◆ meEbEMean_

MonitorElement* EcalSelectiveReadoutValidation::meEbEMean_

private

Definition at line 564 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and bookHistograms().

◆ meEbFixedPayload_

MonitorElement* EcalSelectiveReadoutValidation::meEbFixedPayload_

private

Event payload that do not depend on the number of crystals passing the SR

Definition at line 631 of file EcalSelectiveReadoutValidation.h.

◆ meEbFullRoCnt_

MonitorElement* EcalSelectiveReadoutValidation::meEbFullRoCnt_

private

Definition at line 591 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and bookHistograms().

◆ meEbHiZsFir_

MonitorElement* EcalSelectiveReadoutValidation::meEbHiZsFir_

private

Definition at line 595 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and bookHistograms().

◆ meEbIncompleteRUZsFir_

MonitorElement* EcalSelectiveReadoutValidation::meEbIncompleteRUZsFir_

private

Definition at line 596 of file EcalSelectiveReadoutValidation.h.

◆ meEbLiZsFir_

MonitorElement* EcalSelectiveReadoutValidation::meEbLiZsFir_

private

Definition at line 594 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and bookHistograms().

◆ meEbNoise_

MonitorElement* EcalSelectiveReadoutValidation::meEbNoise_

private

Definition at line 565 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and bookHistograms().

◆ meEbNoZsRecVsSimE_

MonitorElement* EcalSelectiveReadoutValidation::meEbNoZsRecVsSimE_

private

Definition at line 569 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and bookHistograms().

◆ meEbRecE_

MonitorElement* EcalSelectiveReadoutValidation::meEbRecE_

private

Definition at line 563 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and bookHistograms().

◆ meEbRecEHitXtal_

MonitorElement* EcalSelectiveReadoutValidation::meEbRecEHitXtal_

private

Definition at line 567 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and bookHistograms().

◆ meEbRecVsSimE_

MonitorElement* EcalSelectiveReadoutValidation::meEbRecVsSimE_

private

Definition at line 568 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and bookHistograms().

◆ meEbSimE_

MonitorElement* EcalSelectiveReadoutValidation::meEbSimE_

private

Definition at line 566 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and bookHistograms().

◆ meEbZsErrCnt_

MonitorElement* EcalSelectiveReadoutValidation::meEbZsErrCnt_

private

Definition at line 619 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and bookHistograms().

◆ meEbZsErrType1Cnt_

MonitorElement* EcalSelectiveReadoutValidation::meEbZsErrType1Cnt_

private

Definition at line 622 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and bookHistograms().

◆ meEeEMean_

MonitorElement* EcalSelectiveReadoutValidation::meEeEMean_

private

Definition at line 572 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE(), and bookHistograms().

◆ meEeFixedPayload_

MonitorElement* EcalSelectiveReadoutValidation::meEeFixedPayload_

private

Definition at line 632 of file EcalSelectiveReadoutValidation.h.

◆ meEeFullRoCnt_

MonitorElement* EcalSelectiveReadoutValidation::meEeFullRoCnt_

private

Definition at line 592 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and bookHistograms().

◆ meEeHiZsFir_

MonitorElement* EcalSelectiveReadoutValidation::meEeHiZsFir_

private

Definition at line 599 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE(), and bookHistograms().

◆ meEeLiZsFir_

MonitorElement* EcalSelectiveReadoutValidation::meEeLiZsFir_

private

Definition at line 598 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE(), and bookHistograms().

◆ meEeNoise_

MonitorElement* EcalSelectiveReadoutValidation::meEeNoise_

private

Definition at line 573 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE(), and bookHistograms().

◆ meEeNoZsRecVsSimE_

MonitorElement* EcalSelectiveReadoutValidation::meEeNoZsRecVsSimE_

private

Definition at line 577 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE(), and bookHistograms().

◆ meEeRecE_

MonitorElement* EcalSelectiveReadoutValidation::meEeRecE_

private

Definition at line 571 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE(), and bookHistograms().

◆ meEeRecEHitXtal_

MonitorElement* EcalSelectiveReadoutValidation::meEeRecEHitXtal_

private

Definition at line 575 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE(), and bookHistograms().

◆ meEeRecVsSimE_

MonitorElement* EcalSelectiveReadoutValidation::meEeRecVsSimE_

private

Definition at line 576 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE(), and bookHistograms().

◆ meEeSimE_

MonitorElement* EcalSelectiveReadoutValidation::meEeSimE_

private

Definition at line 574 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE(), and bookHistograms().

◆ meEeZsErrCnt_

MonitorElement* EcalSelectiveReadoutValidation::meEeZsErrCnt_

private

Definition at line 620 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and bookHistograms().

◆ meEeZsErrType1Cnt_

MonitorElement* EcalSelectiveReadoutValidation::meEeZsErrType1Cnt_

private

Definition at line 623 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and bookHistograms().

◆ meFixedPayload_

MonitorElement* EcalSelectiveReadoutValidation::meFixedPayload_

private

Definition at line 633 of file EcalSelectiveReadoutValidation.h.

◆ meForcedRu_

MonitorElement* EcalSelectiveReadoutValidation::meForcedRu_

private

Definition at line 581 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), analyzeEE(), and bookHistograms().

◆ meForcedTtf_

MonitorElement* EcalSelectiveReadoutValidation::meForcedTtf_

private

Definition at line 586 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and bookHistograms().

◆ meFullRoCnt_

MonitorElement* EcalSelectiveReadoutValidation::meFullRoCnt_

private

Definition at line 590 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and bookHistograms().

◆ meFullRoRu_

MonitorElement* EcalSelectiveReadoutValidation::meFullRoRu_

private

Definition at line 579 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), analyzeEE(), and bookHistograms().

◆ meHiTtf_

MonitorElement* EcalSelectiveReadoutValidation::meHiTtf_

private

Definition at line 585 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and bookHistograms().

◆ meIncompleteFRO_

MonitorElement* EcalSelectiveReadoutValidation::meIncompleteFRO_

private

Definition at line 604 of file EcalSelectiveReadoutValidation.h.

◆ meIncompleteFROCnt_

MonitorElement* EcalSelectiveReadoutValidation::meIncompleteFROCnt_

private

Definition at line 616 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and bookHistograms().

◆ meIncompleteFROMap_

MonitorElement* EcalSelectiveReadoutValidation::meIncompleteFROMap_

private

Definition at line 608 of file EcalSelectiveReadoutValidation.h.

Referenced by bookHistograms(), and checkSrApplication().

◆ meIncompleteFRORateMap_

MonitorElement* EcalSelectiveReadoutValidation::meIncompleteFRORateMap_

private

Definition at line 612 of file EcalSelectiveReadoutValidation.h.

Referenced by bookHistograms(), and checkSrApplication().

◆ meL1aRate_

MonitorElement* EcalSelectiveReadoutValidation::meL1aRate_

private

Estimate of L1A rate

Definition at line 638 of file EcalSelectiveReadoutValidation.h.

Referenced by bookHistograms(), and dqmEndRun().

◆ meLiTtf_

MonitorElement* EcalSelectiveReadoutValidation::meLiTtf_

private

Definition at line 583 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and bookHistograms().

◆ meMiTtf_

MonitorElement* EcalSelectiveReadoutValidation::meMiTtf_

private

Definition at line 584 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and bookHistograms().

◆ meSRFlagsComputed_

MonitorElement* EcalSelectiveReadoutValidation::meSRFlagsComputed_

private

Definition at line 601 of file EcalSelectiveReadoutValidation.h.

◆ meSRFlagsConsistency_

MonitorElement* EcalSelectiveReadoutValidation::meSRFlagsConsistency_

private

Definition at line 602 of file EcalSelectiveReadoutValidation.h.

Referenced by bookHistograms(), and compareSrfColl().

◆ meSRFlagsFromData_

MonitorElement* EcalSelectiveReadoutValidation::meSRFlagsFromData_

private

Definition at line 600 of file EcalSelectiveReadoutValidation.h.

◆ meTp_

MonitorElement* EcalSelectiveReadoutValidation::meTp_

private

Definition at line 557 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and bookHistograms().

◆ meTpMap_

MonitorElement* EcalSelectiveReadoutValidation::meTpMap_

private

Definition at line 588 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and bookHistograms().

◆ meTpVsEtSum_

MonitorElement* EcalSelectiveReadoutValidation::meTpVsEtSum_

private

Definition at line 561 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and bookHistograms().

◆ meTtf_

MonitorElement* EcalSelectiveReadoutValidation::meTtf_

private

Definition at line 558 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and bookHistograms().

◆ meTtfVsEtSum_

MonitorElement* EcalSelectiveReadoutValidation::meTtfVsEtSum_

private

Definition at line 560 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and bookHistograms().

◆ meTtfVsTp_

MonitorElement* EcalSelectiveReadoutValidation::meTtfVsTp_

private

Definition at line 559 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and bookHistograms().

◆ meVol_

MonitorElement* EcalSelectiveReadoutValidation::meVol_

private

Definition at line 546 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meVolB_

MonitorElement* EcalSelectiveReadoutValidation::meVolB_

private

Definition at line 547 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meVolBHI_

MonitorElement* EcalSelectiveReadoutValidation::meVolBHI_

private

Definition at line 552 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meVolBLI_

MonitorElement* EcalSelectiveReadoutValidation::meVolBLI_

private

Definition at line 549 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meVolE_

MonitorElement* EcalSelectiveReadoutValidation::meVolE_

private

Definition at line 548 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meVolEHI_

MonitorElement* EcalSelectiveReadoutValidation::meVolEHI_

private

Definition at line 553 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meVolELI_

MonitorElement* EcalSelectiveReadoutValidation::meVolELI_

private

Definition at line 550 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meVolHI_

MonitorElement* EcalSelectiveReadoutValidation::meVolHI_

private

Definition at line 554 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meVolLI_

MonitorElement* EcalSelectiveReadoutValidation::meVolLI_

private

Definition at line 551 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), and bookHistograms().

◆ meZs1Ru_

MonitorElement* EcalSelectiveReadoutValidation::meZs1Ru_

private

Definition at line 580 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), analyzeEE(), and bookHistograms().

◆ meZsErrCnt_

MonitorElement* EcalSelectiveReadoutValidation::meZsErrCnt_

private

Definition at line 621 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and bookHistograms().

◆ meZsErrType1Cnt_

MonitorElement* EcalSelectiveReadoutValidation::meZsErrType1Cnt_

private

Definition at line 624 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and bookHistograms().

◆ minDccId_

const int EcalSelectiveReadoutValidation::minDccId_ = 1

staticprivate

Definition at line 426 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), analyzeDataVolume(), checkSrApplication(), compareSrfColl(), dccCh(), and selectFedsForLog().

◆ nCompleteZS_

int EcalSelectiveReadoutValidation::nCompleteZS_

private

Counter of ZS-flagged RU fully read out.

Definition at line 647 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and checkSrApplication().

◆ nDccChs_

const unsigned EcalSelectiveReadoutValidation::nDccChs_ = 68

staticprivate

Number of input channels of a DCC.

Definition at line 423 of file EcalSelectiveReadoutValidation.h.

◆ nDccRus_

const int EcalSelectiveReadoutValidation::nDccRus_

staticprivate

number of RUs for each DCC

Definition at line 444 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume().

◆ nDccs_

const unsigned EcalSelectiveReadoutValidation::nDccs_ = 54

staticprivate

Total number of DCCs.

Definition at line 419 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), bookHistograms(), getEeEventSize(), and selectFedsForLog().

◆ nDroppedFRO_

int EcalSelectiveReadoutValidation::nDroppedFRO_

private

Counter of FRO-flagged RU dropped from data.

Definition at line 641 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and checkSrApplication().

◆ nEb_

int EcalSelectiveReadoutValidation::nEb_

private

ECAL barrel read channel count

Definition at line 720 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), and analyzeDataVolume().

◆ nEbDccs

const int EcalSelectiveReadoutValidation::nEbDccs = 36

staticprivate

number of DCCs for EB

Definition at line 432 of file EcalSelectiveReadoutValidation.h.

Referenced by getEbEventSize(), and getEeEventSize().

◆ nEbEta

const int EcalSelectiveReadoutValidation::nEbEta = 170

staticprivate

number of crystals along Eta in EB

Definition at line 450 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB().

◆ nEbFROCnt_

int EcalSelectiveReadoutValidation::nEbFROCnt_

private

Counter of EB FRO-flagged RUs.

Definition at line 650 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and analyzeEB().

◆ nEbHI_

int EcalSelectiveReadoutValidation::nEbHI_

private

ECAL barrel high interest read channel count

Definition at line 740 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), and analyzeDataVolume().

◆ nEbLI_

int EcalSelectiveReadoutValidation::nEbLI_

private

ECAL barrel low interest read channel count

Definition at line 736 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), and analyzeDataVolume().

◆ nEbPhi

const int EcalSelectiveReadoutValidation::nEbPhi = 360

staticprivate

number of crystals along Phi in EB

Definition at line 453 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB().

◆ nEbRus

const int EcalSelectiveReadoutValidation::nEbRus = 36 * 68

staticprivate

number of RUs for EB

Definition at line 438 of file EcalSelectiveReadoutValidation.h.

◆ nEbTtEta

const int EcalSelectiveReadoutValidation::nEbTtEta = 34

staticprivate

Number of Trigger Towers in barrel along Eta.

Definition at line 471 of file EcalSelectiveReadoutValidation.h.

◆ nEbZsErrors_

int EcalSelectiveReadoutValidation::nEbZsErrors_

private

Counter of EB ZS errors (LI channel below ZS threshold)

Definition at line 656 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and analyzeEB().

◆ nEbZsErrorsType1_

int EcalSelectiveReadoutValidation::nEbZsErrorsType1_

private

Counter of EB ZS errors of type 1: LI channel below ZS threshold and in a RU which was fully readout

Definition at line 663 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and analyzeEB().

◆ nEe_

int EcalSelectiveReadoutValidation::nEe_

private

ECAL endcap read channel count

Definition at line 724 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), and analyzeDataVolume().

◆ nEeDccs

const int EcalSelectiveReadoutValidation::nEeDccs = 18

staticprivate

number of DCCs for EE

Definition at line 435 of file EcalSelectiveReadoutValidation.h.

Referenced by getEbEventSize(), and getEeEventSize().

◆ nEeFROCnt_

int EcalSelectiveReadoutValidation::nEeFROCnt_

private

Counter of EE FRO-flagged RUs.

Definition at line 653 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and analyzeEE().

◆ nEeHI_

int EcalSelectiveReadoutValidation::nEeHI_

private

ECAL endcap high interest read channel count

Definition at line 732 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), and analyzeDataVolume().

◆ nEeLI_

int EcalSelectiveReadoutValidation::nEeLI_

private

ECAL endcap low interest read channel count

Definition at line 728 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), and analyzeDataVolume().

◆ nEeRus

const int EcalSelectiveReadoutValidation::nEeRus = 2 * (34 + 32 + 33 + 33 + 32 + 34 + 33 + 34 + 33)

staticprivate

number of RUs for EE

Definition at line 441 of file EcalSelectiveReadoutValidation.h.

◆ nEeX

const int EcalSelectiveReadoutValidation::nEeX = 100

staticprivate

EE crystal grid size along X.

Definition at line 456 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE().

◆ nEeY

const int EcalSelectiveReadoutValidation::nEeY = 100

staticprivate

EE crystal grid size along Y.

Definition at line 459 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE().

◆ nEeZsErrors_

int EcalSelectiveReadoutValidation::nEeZsErrors_

private

Counter of EE ZS errors (LI channel below ZS threshold)

Definition at line 659 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and analyzeEE().

◆ nEeZsErrorsType1_

int EcalSelectiveReadoutValidation::nEeZsErrorsType1_

private

Counter of EE ZS errors of tyoe 1: LI channel below ZS threshold and in a RU which was fully readout

Definition at line 667 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and analyzeEE().

◆ nEndcaps

const int EcalSelectiveReadoutValidation::nEndcaps = 2

staticprivate

number of endcaps

Definition at line 447 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEE().

◆ nHiPerDcc_

int EcalSelectiveReadoutValidation::nHiPerDcc_[nDccs_]

private

read-out ECAL Hiugh interest channel count for each DCC:

Definition at line 752 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), and analyzeDataVolume().

◆ nHiRuPerDcc_

int EcalSelectiveReadoutValidation::nHiRuPerDcc_[nDccs_]

private

Count for each DCC of HI RUs with at leat one channel read out:

Definition at line 764 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), and analyzeDataVolume().

◆ nIncompleteFRO_

int EcalSelectiveReadoutValidation::nIncompleteFRO_

private

Counter of FRO-flagged RU only partial data.

Definition at line 644 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and checkSrApplication().

◆ nLiPerDcc_

int EcalSelectiveReadoutValidation::nLiPerDcc_[nDccs_]

private

read-out ECAL Low interest channel count for each DCC:

Definition at line 748 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), and analyzeDataVolume().

◆ nLiRuPerDcc_

int EcalSelectiveReadoutValidation::nLiRuPerDcc_[nDccs_]

private

Count for each DCC of LI RUs with at leat one channel read out:

Definition at line 760 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), and analyzeDataVolume().

◆ nMaxXtalPerRu

const int EcalSelectiveReadoutValidation::nMaxXtalPerRu = 25

staticprivate

Number of crystals per Readout Unit excepted partial SCs.

Definition at line 480 of file EcalSelectiveReadoutValidation.h.

◆ nOneEeTtEta

const int EcalSelectiveReadoutValidation::nOneEeTtEta = 11

staticprivate

Number of Trigger Towers in an endcap along Eta.

Definition at line 468 of file EcalSelectiveReadoutValidation.h.

Referenced by dccCh().

◆ nPerDcc_

int EcalSelectiveReadoutValidation::nPerDcc_[nDccs_]

private

read-out ECAL channel count for each DCC:

Definition at line 744 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), and analyzeDataVolume().

◆ nPerRu_

int EcalSelectiveReadoutValidation::nPerRu_[nDccs_][nDccChs_]

private

Number of crystal read for each DCC channel (aka readout unit).

Definition at line 777 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), analyzeDataVolume(), and checkSrApplication().

◆ nRuPerDcc_

int EcalSelectiveReadoutValidation::nRuPerDcc_[nDccs_]

private

Count for each DCC of RUs with at leat one channel read out:

Definition at line 756 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi(), anaDigiInit(), and getRuCount().

◆ nTtEta

const int EcalSelectiveReadoutValidation::nTtEta = 2 * nOneEeTtEta + nEbTtEta

staticprivate

Number of Trigger Towers along Eta.

Definition at line 474 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and setTtEtSums().

◆ nTtPhi

const int EcalSelectiveReadoutValidation::nTtPhi = 72

staticprivate

Number of Trigger Towers along Phi.

Definition at line 477 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and setTtEtSums().

◆ outputFile_

std::string EcalSelectiveReadoutValidation::outputFile_

private

Output file for histograms.

Definition at line 488 of file EcalSelectiveReadoutValidation.h.

Referenced by EcalSelectiveReadoutValidation().

◆ physHandle

edm::ESGetToken<EcalTPGPhysicsConst, EcalTPGPhysicsConstRcd> EcalSelectiveReadoutValidation::physHandle

private

Definition at line 411 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP().

◆ rad2deg

const double EcalSelectiveReadoutValidation::rad2deg = 45. / atan(1.)

staticprivate

Conversion factor from radian to degree.

Definition at line 482 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), and analyzeEE().

◆ scEdge

const int EcalSelectiveReadoutValidation::scEdge = 5

staticprivate

Number of crystals along a supercrystal edge.

Definition at line 465 of file EcalSelectiveReadoutValidation.h.

Referenced by anaDigi().

◆ SkipInnerSC_

bool EcalSelectiveReadoutValidation::SkipInnerSC_

private

Permits to skip inner SC

Definition at line 791 of file EcalSelectiveReadoutValidation.h.

◆ srApplicationErrorLog_

std::ofstream EcalSelectiveReadoutValidation::srApplicationErrorLog_

private

Output ascii file for unconsistency between Xtals and RU Flags.

Definition at line 505 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeEB(), analyzeEE(), checkSrApplication(), initAsciiFile(), and selectFedsForLog().

◆ srApplicationErrorLogFileName_

std::string EcalSelectiveReadoutValidation::srApplicationErrorLogFileName_

private

Output ascii file name for unconsistency between SRFs and actual number of read-out crystals.

Definition at line 499 of file EcalSelectiveReadoutValidation.h.

Referenced by EcalSelectiveReadoutValidation(), and initAsciiFile().

◆ srpAlgoErrorLog_

std::ofstream EcalSelectiveReadoutValidation::srpAlgoErrorLog_

private

Output ascii file for unconsistency on SR flags.

Definition at line 502 of file EcalSelectiveReadoutValidation.h.

Referenced by compareSrfColl(), initAsciiFile(), and selectFedsForLog().

◆ srpAlgoErrorLogFileName_

std::string EcalSelectiveReadoutValidation::srpAlgoErrorLogFileName_

private

Output ascii file name for unconsistency between SRFs read from data and SRF obtained by rerunning SRP algorithm on TTFs.

Definition at line 495 of file EcalSelectiveReadoutValidation.h.

Referenced by EcalSelectiveReadoutValidation(), and initAsciiFile().

◆ tmax

int64_t EcalSelectiveReadoutValidation::tmax

private

For L1A rate estimate

Definition at line 532 of file EcalSelectiveReadoutValidation.h.

Referenced by getL1aRate(), and updateL1aRate().

◆ tmin

int64_t EcalSelectiveReadoutValidation::tmin

private

Definition at line 533 of file EcalSelectiveReadoutValidation.h.

Referenced by getL1aRate(), and updateL1aRate().

◆ tpInGeV_

bool EcalSelectiveReadoutValidation::tpInGeV_

private

Switch for uncompressing TP value

Definition at line 700 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and bookHistograms().

◆ tps_

CollHandle<EcalTrigPrimDigiCollection> EcalSelectiveReadoutValidation::tps_

private

Definition at line 523 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), analyzeTP(), EcalSelectiveReadoutValidation(), and readAllCollections().

◆ triggerTowerMap_

const EcalTrigTowerConstituentsMap* EcalSelectiveReadoutValidation::triggerTowerMap_

private

ECAL trigger tower mapping

Definition at line 671 of file EcalSelectiveReadoutValidation.h.

Referenced by dqmBeginRun(), readOutUnitOf(), and setTtEtSums().

◆ ttEtSums

double EcalSelectiveReadoutValidation::ttEtSums[nTtEta][nTtPhi]

private

Trigger tower Et computed as sum the crystal Et. Indices stands for the eta and phi TT index starting from 0 at eta minimum and at phi=0+ in std CMS coordinate system.

Definition at line 787 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP(), and setTtEtSums().

◆ useEventRate_

bool EcalSelectiveReadoutValidation::useEventRate_

private

Switch to fill histograms with event rate instead of event count. Applies only to some histograms.

Definition at line 710 of file EcalSelectiveReadoutValidation.h.

◆ verbose_

bool EcalSelectiveReadoutValidation::verbose_

private

Verbosity switch.

Definition at line 485 of file EcalSelectiveReadoutValidation.h.

Referenced by EcalSelectiveReadoutValidation().

◆ weights_

std::vector<double> EcalSelectiveReadoutValidation::weights_

private

Weights for amplitude local reconstruction

Definition at line 684 of file EcalSelectiveReadoutValidation.h.

Referenced by frame2Energy().

◆ withEbSimHit_

bool EcalSelectiveReadoutValidation::withEbSimHit_

private

Indicates if EB sim hits are available

Definition at line 817 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and analyzeEB().

◆ withEeSimHit_

bool EcalSelectiveReadoutValidation::withEeSimHit_

private

Indicates if EE sim hits are available

Definition at line 813 of file EcalSelectiveReadoutValidation.h.

Referenced by analyze(), and analyzeEE().

◆ zsErrorLog_

std::ofstream EcalSelectiveReadoutValidation::zsErrorLog_

private

File to log ZS and other errors.

Definition at line 508 of file EcalSelectiveReadoutValidation.h.

Classes

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Static Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ RecHit

◆ RecHitCollection

Member Enumeration Documentation

◆ subdet_t

Constructor & Destructor Documentation

◆ EcalSelectiveReadoutValidation()

◆ ~EcalSelectiveReadoutValidation()

Member Function Documentation

◆ anaDigi()

◆ anaDigiInit()

◆ analyze()

◆ analyzeDataVolume()

◆ analyzeEB()

◆ analyzeEE()

◆ analyzeTP()

◆ book1D()

◆ book2D()

◆ bookFloat()

◆ bookHistograms()

◆ bookProfile()

◆ bookProfile2D()

◆ checkSrApplication()

◆ cIndex2iEta()

◆ cIndex2iPhi()

◆ cIndex2iTtEta()

◆ cIndex2iTtPhi()

◆ cIndex2iXY()

◆ compareSrfColl()

◆ configFirWeights()

◆ dccCh()

◆ dccId() [1/2]

◆ dccId() [2/2]

◆ dccZsFIR()

◆ dqmBeginRun()

◆ dqmEndRun()

◆ fill() [1/4]

◆ fill() [2/4]

◆ fill() [3/4]

◆ fill() [4/4]

◆ frame2Energy()

◆ frame2EnergyForTp()

◆ getBytesPerCrystal()

◆ getCrystalCount()

◆ getDccEventSize()

◆ getDccOverhead()

◆ getDccSrDependentPayload()

◆ getEbEventSize()

◆ getEeEventSize()

◆ getFIRWeights()

◆ getL1aRate()

◆ getRuCount()

◆ iEta2cIndex()

◆ initAsciiFile()

◆ iPhi2cIndex()

◆ iTtEta2cIndex()

◆ iTtPhi2cIndex()

◆ iXY2cIndex()

◆ printAvailableHists()

◆ readAllCollections()

◆ readOutUnitOf() [1/2]

◆ readOutUnitOf() [2/2]

◆ registerHist()

◆ ruGraphX() [1/2]

◆ ruGraphX() [2/2]

◆ ruGraphY() [1/2]

◆ ruGraphY() [2/2]

◆ selectFedsForLog()

◆ setTtEtSums()

◆ updateL1aRate()