#include <PresampleClient.h>

Inheritance diagram for ecaldqm::PresampleClient:

Public Member Functions
	PresampleClient ()

void	producePlots (ProcessType) override

	~PresampleClient () override

Public Member Functions inherited from ecaldqm::DQWorkerClient
void	bookMEs (DQMStore::IBooker &) override

	DQWorkerClient ()

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

void	releaseMEs () override

void	releaseSource ()

virtual void	resetMEs ()

void	resetPerLumi ()

bool	retrieveSource (DQMStore::IGetter &, ProcessType)

bool	runsOn (ProcessType _type) const

void	setStatusManager (StatusManager const &_manager)

virtual void	setTokens (edm::ConsumesCollector &)

	~DQWorkerClient () override

Public Member Functions inherited from ecaldqm::DQWorker
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

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

bool	checkElectronicsMap (bool=true)

bool	checkGeometry (bool=true)

bool	checkTopology (bool=true)

bool	checkTrigTowerMap (bool=true)

	DQWorker ()

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

EcalDQMSetupObjects const	getEcalDQMSetupObjects ()

EcalElectronicsMapping const *	GetElectronicsMap ()

CaloGeometry const *	GetGeometry ()

std::string const &	getName () const

CaloTopology const *	GetTopology ()

EcalTrigTowerConstituentsMap const *	GetTrigTowerMap ()

bool	onlineMode () const

void	setEventNumber (edm::EventNumber_t _e)

void	setLumiNumber (edm::LuminosityBlockNumber_t _l)

void	setRunNumber (edm::RunNumber_t _r)

void	setSetupObjects (edm::EventSetup const &)

void	setSetupObjectsEndLumi (edm::EventSetup const &)

void	setTime (time_t _t)

void	setTokens (edm::ConsumesCollector &)

virtual	~DQWorker () noexcept(false)

Private Member Functions
void	setParams (edm::ParameterSet const &) override

Private Attributes
float	expectedMean_

int	minChannelEntries_

float	toleranceHigh_

float	toleranceHighFwd_

float	toleranceLow_

float	toleranceRMS_

float	toleranceRMSFwd_

Additional Inherited Members
Public Types inherited from ecaldqm::DQWorkerClient
enum	ProcessType { kLumi, kJob, nProcessType }

enum	Quality { kBad = 0, kGood = 1, kUnknown = 2, kMBad = 3, kMGood = 4, kMUnknown = 5 }

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

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

Public Attributes inherited from ecaldqm::DQWorker
edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	elecMapHandle

edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	elecMapHandleEndLumi

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geomHandle

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geomHandleEndLumi

edm::ESGetToken< CaloTopology, CaloTopologyRecord >	topoHandle

edm::ESGetToken< CaloTopology, CaloTopologyRecord >	topoHandleEndLumi

edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord >	ttMapHandle

edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord >	ttMapHandleEndLumi

Protected Types inherited from ecaldqm::DQWorker
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

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

void	setSource (edm::ParameterSet const &) override

void	towerAverage_ (MESet &, MESet const &, float)

bool	using_ (std::string const &_name, ProcessType _type=kJob) const

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

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

void	setVerbosity (int _verbosity)

Protected Attributes inherited from ecaldqm::DQWorkerClient
bool	hasLumiPlots_

std::set< std::string >	qualitySummaries_

MESetCollection	sources_

StatusManager const *	statusManager_

Protected Attributes inherited from ecaldqm::DQWorker
bool	booked_

MESetCollection	MEs_

std::string	name_

bool	onlineMode_

Timestamp	timestamp_

int	verbosity_

bool	willConvertToEDM_

Detailed Description

Definition at line 7 of file PresampleClient.h.

Constructor & Destructor Documentation

◆ PresampleClient()

ecaldqm::PresampleClient::PresampleClient ( )

Definition at line 13 of file PresampleClient.cc.

References ecaldqm::DQWorkerClient::qualitySummaries_.

       : DQWorkerClient(),
         minChannelEntries_(0),
         expectedMean_(0.),
         toleranceLow_(0.),
         toleranceHigh_(0.),
         toleranceHighFwd_(0.),
         toleranceRMS_(0.),
         toleranceRMSFwd_(0.) {
     qualitySummaries_.insert("Quality");
     qualitySummaries_.insert("QualitySummary");
   }

◆ ~PresampleClient()

ecaldqm::PresampleClient::~PresampleClient ( )

inlineoverride

Definition at line 10 of file PresampleClient.h.

10 {}

Member Function Documentation

◆ producePlots()

void ecaldqm::PresampleClient::producePlots ( ProcessType )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 36 of file PresampleClient.cc.

References ecaldqm::MESetCollection::at(), ecaldqm::dccId(), EcalBarrel, EcalEndcap, expectedMean_, ecaldqm::MESet::fill(), ecaldqm::MESet::ConstBin::getBinContent(), ecaldqm::MESet::ConstBin::getBinEntries(), ecaldqm::MESet::ConstBin::getBinError(), ecaldqm::DQWorker::getEcalDQMSetupObjects(), ecaldqm::DQWorker::GetElectronicsMap(), ecaldqm::DQWorker::GetTrigTowerMap(), triggerObjects_cff::id, ecaldqm::DQWorker::Timestamp::iLumi, ecaldqm::isForward(), ecaldqm::DQWorkerClient::kBad, ecaldqm::DQWorkerClient::kGood, ecaldqm::DQWorkerClient::kMBad, ecaldqm::DQWorkerClient::kMGood, ecaldqm::DQWorkerClient::kMUnknown, EcalChannelStatusCode::kOk, ecaldqm::DQWorkerClient::kUnknown, gpuClustering::pixelStatus::mask, SiStripPI::mean, ecaldqm::DQWorker::MEs_, minChannelEntries_, EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR, EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR, SiStripPI::rms, ecaldqm::DQWorkerClient::sources_, mathSSE::sqrt(), ecaldqm::DQWorkerClient::statusManager_, ecaldqm::DQWorker::timestamp_, toleranceHigh_, toleranceHighFwd_, toleranceLow_, toleranceRMS_, toleranceRMSFwd_, ecaldqm::MESet::iterator::toNextChannel(), and ecaldqm::DQWorkerClient::towerAverage_().

                                                 {
     MESet& meQualitySummary(MEs_.at("QualitySummary"));
     MESet& meQuality(MEs_.at("Quality"));
     MESet& meErrorsSummary(MEs_.at("ErrorsSummary"));
     MESet& meMean(MEs_.at("Mean"));
     MESet& meRMS(MEs_.at("RMS"));
     MESet& meRMSMap(MEs_.at("RMSMap"));
     MESet& meRMSMapAll(MEs_.at("RMSMapAll"));
     MESet& meRMSMapAllByLumi(MEs_.at("RMSMapAllByLumi"));
     MESet& meMeanMapAll(MEs_.at("MeanMapAll"));
 
     MESet const& sPedestal(sources_.at("Pedestal"));
     MESet const& sPedestalByLS(sources_.at("PedestalByLS"));
     MESet const& sChStatus(sources_.at("ChStatus"));
 
     uint32_t mask(1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR |
                   1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR);
 
     MESet::iterator qEnd(meQuality.end(GetElectronicsMap()));
 
     MESet::const_iterator pItr(GetElectronicsMap(), sPedestal);
     MESet::const_iterator pLSItr(GetElectronicsMap(), sPedestalByLS);
     double maxEB(0.), minEB(0.), maxEE(0.), minEE(0.);
     double rmsMaxEB(0.), rmsMaxEE(0.);
     for (MESet::iterator qItr(meQuality.beginChannel(GetElectronicsMap())); qItr != qEnd;
          qItr.toNextChannel(GetElectronicsMap())) {
       pItr = qItr;
       pLSItr = qItr;
 
       DetId id(qItr->getId());
 
       bool doMask(meQuality.maskMatches(id, mask, statusManager_, GetTrigTowerMap()));
 
       double rmsThresh(toleranceRMS_);
       double meanThreshHigh(toleranceHigh_);
 
       if (isForward(id)) {
         rmsThresh = toleranceRMSFwd_;
         meanThreshHigh = toleranceHighFwd_;
       }
 
       double entries(pItr->getBinEntries());
       double entriesLS(pLSItr->getBinEntries());
 
       if (entries < minChannelEntries_) {
         qItr->setBinContent(doMask ? kMUnknown : kUnknown);
         meQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, doMask ? kMUnknown : kUnknown);
         meRMSMap.setBinContent(getEcalDQMSetupObjects(), id, -1.);
         continue;
       }
 
       double mean(pItr->getBinContent());
       double meanLS(pLSItr->getBinContent());
       double rms(pItr->getBinError() * std::sqrt(entries));
       double rmsLS(pLSItr->getBinError() * std::sqrt(entriesLS));
 
       int dccid(dccId(id, GetElectronicsMap()));
 
       meMean.fill(getEcalDQMSetupObjects(), dccid, mean);
       meRMS.fill(getEcalDQMSetupObjects(), dccid, rms);
       meRMSMap.setBinContent(getEcalDQMSetupObjects(), id, rms);
       meRMSMapAllByLumi.setBinContent(getEcalDQMSetupObjects(), id, rmsLS);
 
       if (((mean > expectedMean_ + meanThreshHigh) || (mean < expectedMean_ - toleranceLow_)) || rms > rmsThresh) {
         qItr->setBinContent(doMask ? kMBad : kBad);
         meQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, doMask ? kMBad : kBad);
         if (!doMask)
           meErrorsSummary.fill(getEcalDQMSetupObjects(), id);
       } else {
         qItr->setBinContent(doMask ? kMGood : kGood);
         meQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, doMask ? kMGood : kGood);
       }
 
       // Fill Presample Trend plots:
       // Use PedestalByLS which only contains digis from "current" LS
       float chStatus(sChStatus.getBinContent(getEcalDQMSetupObjects(), id));
       if (entriesLS < minChannelEntries_)
         continue;
       if (chStatus != EcalChannelStatusCode::kOk)
         continue;  // exclude problematic channels
 
       // Get max/min
       // Min is effectively just 0
       if (id.subdetId() == EcalBarrel) {
         if (meanLS > maxEB)
           maxEB = meanLS;
         if (meanLS < minEB)
           minEB = meanLS;
         if (rmsLS > rmsMaxEB)
           rmsMaxEB = rmsLS;
       } else {
         if (meanLS > maxEE)
           maxEE = meanLS;
         if (meanLS < minEE)
           minEE = meanLS;
         if (rmsLS > rmsMaxEE)
           rmsMaxEE = rmsLS;
       }
 
     }  // qItr
 
     towerAverage_(meRMSMapAll, meRMSMap, -1.);
     towerAverage_(meMeanMapAll, sPedestal, -1.);
 
     MESet& meTrendMean(MEs_.at("TrendMean"));
     MESet& meTrendRMS(MEs_.at("TrendRMS"));
     meTrendMean.fill(getEcalDQMSetupObjects(), EcalBarrel, double(timestamp_.iLumi), maxEB - minEB);
     meTrendMean.fill(getEcalDQMSetupObjects(), EcalEndcap, double(timestamp_.iLumi), maxEE - minEE);
     meTrendRMS.fill(getEcalDQMSetupObjects(), EcalBarrel, double(timestamp_.iLumi), rmsMaxEB);
     meTrendRMS.fill(getEcalDQMSetupObjects(), EcalEndcap, double(timestamp_.iLumi), rmsMaxEE);
   }

◆ setParams()

void ecaldqm::PresampleClient::setParams ( edm::ParameterSet const & _params )

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 26 of file PresampleClient.cc.

References l1extraParticles_cfi::_params, expectedMean_, minChannelEntries_, toleranceHigh_, toleranceHighFwd_, toleranceLow_, toleranceRMS_, and toleranceRMSFwd_.

                                                                 {
     minChannelEntries_ = _params.getUntrackedParameter<int>("minChannelEntries");
     expectedMean_ = _params.getUntrackedParameter<double>("expectedMean");
     toleranceLow_ = _params.getUntrackedParameter<double>("toleranceLow");
     toleranceHigh_ = _params.getUntrackedParameter<double>("toleranceHigh");
     toleranceHighFwd_ = _params.getUntrackedParameter<double>("toleranceHighFwd");
     toleranceRMS_ = _params.getUntrackedParameter<double>("toleranceRMS");
     toleranceRMSFwd_ = _params.getUntrackedParameter<double>("toleranceRMSFwd");
   }

Member Data Documentation

◆ expectedMean_

float ecaldqm::PresampleClient::expectedMean_

private

Definition at line 18 of file PresampleClient.h.

Referenced by producePlots(), and setParams().

◆ minChannelEntries_

int ecaldqm::PresampleClient::minChannelEntries_

private

Definition at line 17 of file PresampleClient.h.

Referenced by producePlots(), and setParams().

◆ toleranceHigh_

float ecaldqm::PresampleClient::toleranceHigh_

private

Definition at line 20 of file PresampleClient.h.

Referenced by producePlots(), and setParams().

◆ toleranceHighFwd_

float ecaldqm::PresampleClient::toleranceHighFwd_

private

Definition at line 21 of file PresampleClient.h.

Referenced by producePlots(), and setParams().

◆ toleranceLow_

float ecaldqm::PresampleClient::toleranceLow_

private

Definition at line 19 of file PresampleClient.h.

Referenced by producePlots(), and setParams().

◆ toleranceRMS_

float ecaldqm::PresampleClient::toleranceRMS_

private

Definition at line 22 of file PresampleClient.h.

Referenced by producePlots(), and setParams().

◆ toleranceRMSFwd_

float ecaldqm::PresampleClient::toleranceRMSFwd_

private

Definition at line 23 of file PresampleClient.h.

Referenced by producePlots(), and setParams().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PresampleClient()

◆ ~PresampleClient()

Member Function Documentation

◆ producePlots()

◆ setParams()

Member Data Documentation

◆ expectedMean_

◆ minChannelEntries_

◆ toleranceHigh_

◆ toleranceHighFwd_

◆ toleranceLow_

◆ toleranceRMS_

◆ toleranceRMSFwd_