#include <TimingClient.h>

Inheritance diagram for ecaldqm::TimingClient:

Public Member Functions
void	producePlots (ProcessType) override

	TimingClient ()

	~TimingClient () 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	ebtoleranceMean_

float	eetoleranceMean_

int	minChannelEntries_

int	minChannelEntriesFwd_

int	minTowerEntries_

int	minTowerEntriesFwd_

float	tailPopulThreshold_

float	toleranceMeanFwd_

float	toleranceRMS_

float	toleranceRMSFwd_

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

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

Static Public Attributes inherited from ecaldqm::DQWorkerClient
static constexpr int	kBad = 0

static constexpr int	kGood = 1

static constexpr int	kMBad = 3

static constexpr int	kMGood = 4

static constexpr int	kMUnknown = 5

static constexpr int	kUnknown = 2

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 8 of file TimingClient.h.

Constructor & Destructor Documentation

◆ TimingClient()

ecaldqm::TimingClient::TimingClient ( )

Definition at line 13 of file TimingClient.cc.

References ecaldqm::DQWorkerClient::qualitySummaries_.

       : DQWorkerClient(),
         ebtoleranceMean_(0.),
         eetoleranceMean_(0.),
         toleranceMeanFwd_(0.),
         toleranceRMS_(0.),
         toleranceRMSFwd_(0.),
         minChannelEntries_(0),
         minChannelEntriesFwd_(0),
         minTowerEntries_(0),
         minTowerEntriesFwd_(0),
         tailPopulThreshold_(0.) {
     qualitySummaries_.insert("Quality");
     qualitySummaries_.insert("QualitySummary");
   }

◆ ~TimingClient()

ecaldqm::TimingClient::~TimingClient ( )

inlineoverride

Definition at line 11 of file TimingClient.h.

11 {}

Member Function Documentation

◆ producePlots()

void ecaldqm::TimingClient::producePlots ( ProcessType )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 42 of file TimingClient.cc.

References funct::abs(), ecaldqm::MESetCollection::at(), EcalTrigTowerConstituentsMap::constituentsOf(), ebtoleranceMean_, EcalBarrel, EcalEndcap, EcalTriggerTower, eetoleranceMean_, ecaldqm::MESet::fill(), ecaldqm::MESet::ConstBin::getBinContent(), ecaldqm::MESet::ConstBin::getBinEntries(), ecaldqm::MESet::ConstBin::getBinError(), ecaldqm::DQWorker::getEcalDQMSetupObjects(), ecaldqm::DQWorker::GetElectronicsMap(), ecaldqm::DQWorker::GetTrigTowerMap(), EcalPhiSymFlatTableProducers_cfi::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, ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::mask, SiStripPI::mean, ecaldqm::DQWorker::MEs_, LaserClient_cfi::minChannelEntries, minChannelEntries_, minChannelEntriesFwd_, TimingClient_cfi::minTowerEntries, minTowerEntries_, minTowerEntriesFwd_, EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING, quality, SiStripPI::rms, ecaldqm::scConstituents(), ecaldqm::DQWorkerClient::sources_, mathSSE::sqrt(), ecaldqm::DQWorkerClient::statusManager_, EBDetId::switchZSide(), EEDetId::switchZSide(), tailPopulThreshold_, ecaldqm::DQWorker::timestamp_, toleranceMeanFwd_, toleranceRMS_, toleranceRMSFwd_, ecaldqm::MESet::iterator::toNextChannel(), EBDetId::zside(), and EEDetId::zside().

                                              {
     MESet& meQuality(MEs_.at("Quality"));
     MESet& meMeanSM(MEs_.at("MeanSM"));
     MESet& meMeanAll(MEs_.at("MeanAll"));
     MESet& meFwdBkwdDiff(MEs_.at("FwdBkwdDiff"));
     MESet& meFwdvBkwd(MEs_.at("FwdvBkwd"));
     MESet& meRMSMap(MEs_.at("RMSMap"));
     MESet& meRMSAll(MEs_.at("RMSAll"));
     MESet& meProjEta(MEs_.at("ProjEta"));
     MESet& meProjPhi(MEs_.at("ProjPhi"));
     MESet& meQualitySummary(MEs_.at("QualitySummary"));
 
     MESet const& sTimeAllMap(sources_.at("TimeAllMap"));
     MESet const& sTimeMap(sources_.at("TimeMap"));
     MESet const& sTimeMapByLS(sources_.at("TimeMapByLS"));
     MESet const& sChStatus(sources_.at("ChStatus"));
 
     uint32_t mask(1 << EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING);
 
     MESet::iterator qEnd(meQuality.end(GetElectronicsMap()));
 
     MESet::iterator rItr(GetElectronicsMap(), meRMSMap);
     MESet::const_iterator tItr(GetElectronicsMap(), sTimeMap);
     MESet::const_iterator tLSItr(GetElectronicsMap(), sTimeMapByLS);
 
     float EBentries(0.), EEentries(0.);
     float EBmean(0.), EEmean(0.);
     float EBrms(0.), EErms(0.);
     for (MESet::iterator qItr(meQuality.beginChannel(GetElectronicsMap())); qItr != qEnd;
          qItr.toNextChannel(GetElectronicsMap())) {
       tItr = qItr;
       rItr = qItr;
 
       DetId id(qItr->getId());
 
       int minChannelEntries(minChannelEntries_);
       float meanThresh;
       float rmsThresh(toleranceRMS_);
 
       if (id.subdetId() == EcalBarrel)
         meanThresh = ebtoleranceMean_;
       else
         meanThresh = eetoleranceMean_;
 
       if (isForward(id)) {
         minChannelEntries = minChannelEntriesFwd_;
         meanThresh = toleranceMeanFwd_;
         rmsThresh = toleranceRMSFwd_;
       }
 
       bool doMask(meQuality.maskMatches(id, mask, statusManager_, GetTrigTowerMap()));
 
       float entries(tItr->getBinEntries());
 
       if (entries < minChannelEntries) {
         qItr->setBinContent(doMask ? kMUnknown : kUnknown);
         rItr->setBinContent(-1.);
         continue;
       }
 
       float mean(tItr->getBinContent());
       float rms(tItr->getBinError() * sqrt(entries));
 
       meMeanSM.fill(getEcalDQMSetupObjects(), id, mean);
       meMeanAll.fill(getEcalDQMSetupObjects(), id, mean);
       meProjEta.fill(getEcalDQMSetupObjects(), id, mean);
       meProjPhi.fill(getEcalDQMSetupObjects(), id, mean);
       meRMSAll.fill(getEcalDQMSetupObjects(), id, rms);
       rItr->setBinContent(rms);
 
       bool negative(false);
       float posTime(0.);
 
       if (id.subdetId() == EcalBarrel) {
         EBDetId ebid(id);
         if (ebid.zside() < 0) {
           negative = true;
           EBDetId posId(EBDetId::switchZSide(ebid));
           posTime = sTimeMap.getBinContent(getEcalDQMSetupObjects(), posId);
         }
       } else {
         EEDetId eeid(id);
         if (eeid.zside() < 0) {
           negative = true;
           EEDetId posId(EEDetId::switchZSide(eeid));
           posTime = sTimeMap.getBinContent(getEcalDQMSetupObjects(), posId);
         }
       }
       if (negative) {
         meFwdBkwdDiff.fill(getEcalDQMSetupObjects(), id, posTime - mean);
         meFwdvBkwd.fill(getEcalDQMSetupObjects(), id, mean, posTime);
       }
 
       if (std::abs(mean) > meanThresh || rms > rmsThresh)
         qItr->setBinContent(doMask ? kMBad : kBad);
       else
         qItr->setBinContent(doMask ? kMGood : kGood);
 
       // For Trend plots:
       tLSItr = qItr;
       float entriesLS(tLSItr->getBinEntries());
       float meanLS(tLSItr->getBinContent());
       float rmsLS(tLSItr->getBinError() * sqrt(entriesLS));
       int chStatus = static_cast<int>(sChStatus.getBinContent(getEcalDQMSetupObjects(), id));
 
       if (entriesLS < minChannelEntries)
         continue;
       if (chStatus != EcalChannelStatusCode::kOk)
         continue;  // exclude problematic channels
 
       // Keep running count of timing mean, rms, and N_hits
       if (id.subdetId() == EcalBarrel) {
         EBmean += meanLS;
         EBrms += rmsLS;
         EBentries += entriesLS;
       } else {
         EEmean += meanLS;
         EErms += rmsLS;
         EEentries += entriesLS;
       }
 
     }  // channel loop
 
     // Fill Timing Trend plots at each LS
     MESet& meTrendMean(MEs_.at("TrendMean"));
     MESet& meTrendRMS(MEs_.at("TrendRMS"));
     if (EBentries > 0.) {
       if (std::abs(EBmean) > 0.)
         meTrendMean.fill(getEcalDQMSetupObjects(), EcalBarrel, double(timestamp_.iLumi), EBmean / EBentries);
       if (std::abs(EBrms) > 0.)
         meTrendRMS.fill(getEcalDQMSetupObjects(), EcalBarrel, double(timestamp_.iLumi), EBrms / EBentries);
     }
     if (EEentries > 0.) {
       if (std::abs(EEmean) > 0.)
         meTrendMean.fill(getEcalDQMSetupObjects(), EcalEndcap, double(timestamp_.iLumi), EEmean / EEentries);
       if (std::abs(EErms) > 0.)
         meTrendRMS.fill(getEcalDQMSetupObjects(), EcalEndcap, double(timestamp_.iLumi), EErms / EEentries);
     }
 
     MESet::iterator qsEnd(meQualitySummary.end(GetElectronicsMap()));
 
     for (MESet::iterator qsItr(meQualitySummary.beginChannel(GetElectronicsMap())); qsItr != qsEnd;
          qsItr.toNextChannel(GetElectronicsMap())) {
       DetId tId(qsItr->getId());
 
       std::vector<DetId> ids;
 
       if (tId.subdetId() == EcalTriggerTower)
         ids = GetTrigTowerMap()->constituentsOf(EcalTrigTowerDetId(tId));
       else
         ids = scConstituents(EcalScDetId(tId));
 
       int minTowerEntries(minTowerEntries_);
       float meanThresh;
       float rmsThresh(toleranceRMS_);
 
       if (tId.subdetId() == EcalBarrel)
         meanThresh = ebtoleranceMean_;
       else
         meanThresh = eetoleranceMean_;
 
       if (isForward(tId)) {
         minTowerEntries = minTowerEntriesFwd_;
         meanThresh = toleranceMeanFwd_;
         rmsThresh = toleranceRMSFwd_;
       }
 
       // tower entries != sum(channel entries) because of the difference in timing cut at the source
       float summaryEntries(sTimeAllMap.getBinEntries(getEcalDQMSetupObjects(), tId));
 
       float towerEntries(0.);
       float towerMean(0.);
       float towerMean2(0.);
 
       bool doMask(false);
 
       for (std::vector<DetId>::iterator idItr(ids.begin()); idItr != ids.end(); ++idItr) {
         DetId& id(*idItr);
 
         doMask |= meQuality.maskMatches(id, mask, statusManager_, GetTrigTowerMap());
 
         MESet::const_iterator tmItr(GetElectronicsMap(), sTimeMap, id);
 
         float entries(tmItr->getBinEntries());
         if (entries < 0.)
           continue;
         towerEntries += entries;
         float mean(tmItr->getBinContent());
         towerMean += mean * entries;
         float rms(tmItr->getBinError() * sqrt(entries));
         towerMean2 += (rms * rms + mean * mean) * entries;
       }
 
       double quality(doMask ? kMUnknown : kUnknown);
       if (towerEntries / ids.size() > minTowerEntries / 25.) {
         if (summaryEntries < towerEntries * (1. - tailPopulThreshold_))  // large timing deviation
           quality = doMask ? kMBad : kBad;
         else {
           towerMean /= towerEntries;
           towerMean2 /= towerEntries;
 
           float towerRMS(sqrt(towerMean2 - towerMean * towerMean));
 
           if (std::abs(towerMean) > meanThresh || towerRMS > rmsThresh)
             quality = doMask ? kMBad : kBad;
           else
             quality = doMask ? kMGood : kGood;
         }
       }
       qsItr->setBinContent(quality);
     }
   }

◆ setParams()

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 29 of file TimingClient.cc.

References l1extraParticles_cfi::_params, ebtoleranceMean_, eetoleranceMean_, minChannelEntries_, minChannelEntriesFwd_, minTowerEntries_, minTowerEntriesFwd_, tailPopulThreshold_, toleranceMeanFwd_, toleranceRMS_, and toleranceRMSFwd_.

                                                              {
     ebtoleranceMean_ = _params.getUntrackedParameter<double>("ebtoleranceMean");
     eetoleranceMean_ = _params.getUntrackedParameter<double>("eetoleranceMean");
     toleranceMeanFwd_ = _params.getUntrackedParameter<double>("toleranceMeanFwd");
     toleranceRMS_ = _params.getUntrackedParameter<double>("toleranceRMS");
     toleranceRMSFwd_ = _params.getUntrackedParameter<double>("toleranceRMSFwd");
     minChannelEntries_ = _params.getUntrackedParameter<int>("minChannelEntries");
     minChannelEntriesFwd_ = _params.getUntrackedParameter<int>("minChannelEntriesFwd");
     minTowerEntries_ = _params.getUntrackedParameter<int>("minTowerEntries");
     minTowerEntriesFwd_ = _params.getUntrackedParameter<int>("minChannelEntriesFwd");
     tailPopulThreshold_ = _params.getUntrackedParameter<double>("tailPopulThreshold");
   }

Member Data Documentation

◆ ebtoleranceMean_

float ecaldqm::TimingClient::ebtoleranceMean_

private

Definition at line 18 of file TimingClient.h.

Referenced by producePlots(), and setParams().

◆ eetoleranceMean_

float ecaldqm::TimingClient::eetoleranceMean_

private

Definition at line 19 of file TimingClient.h.

Referenced by producePlots(), and setParams().

◆ minChannelEntries_

int ecaldqm::TimingClient::minChannelEntries_

private

Definition at line 23 of file TimingClient.h.

Referenced by producePlots(), and setParams().

◆ minChannelEntriesFwd_

int ecaldqm::TimingClient::minChannelEntriesFwd_

private

Definition at line 24 of file TimingClient.h.

Referenced by producePlots(), and setParams().

◆ minTowerEntries_

int ecaldqm::TimingClient::minTowerEntries_

private

Definition at line 25 of file TimingClient.h.

Referenced by producePlots(), and setParams().

◆ minTowerEntriesFwd_

int ecaldqm::TimingClient::minTowerEntriesFwd_

private

Definition at line 26 of file TimingClient.h.

Referenced by producePlots(), and setParams().

◆ tailPopulThreshold_

float ecaldqm::TimingClient::tailPopulThreshold_

private

Definition at line 27 of file TimingClient.h.

Referenced by producePlots(), and setParams().

◆ toleranceMeanFwd_

float ecaldqm::TimingClient::toleranceMeanFwd_

private

Definition at line 20 of file TimingClient.h.

Referenced by producePlots(), and setParams().

◆ toleranceRMS_

float ecaldqm::TimingClient::toleranceRMS_

private

Definition at line 21 of file TimingClient.h.

Referenced by producePlots(), and setParams().

◆ toleranceRMSFwd_

float ecaldqm::TimingClient::toleranceRMSFwd_

private

Definition at line 22 of file TimingClient.h.

Referenced by producePlots(), and setParams().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TimingClient()

◆ ~TimingClient()

Member Function Documentation

◆ producePlots()

◆ setParams()

Member Data Documentation

◆ ebtoleranceMean_

◆ eetoleranceMean_

◆ minChannelEntries_

◆ minChannelEntriesFwd_

◆ minTowerEntries_

◆ minTowerEntriesFwd_

◆ tailPopulThreshold_

◆ toleranceMeanFwd_

◆ toleranceRMS_

◆ toleranceRMSFwd_