ecaldqm::TimingClient Class Reference

#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 &)
	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
int	minChannelEntries_
int	minChannelEntriesFwd_
int	minTowerEntries_
int	minTowerEntriesFwd_
float	tailPopulThreshold_
float	toleranceMean_
float	toleranceMeanFwd_
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 8 of file TimingClient.h.

Constructor & Destructor Documentation

ecaldqm::TimingClient::TimingClient

(

)

Definition at line 13 of file TimingClient.cc.

References ecaldqm::DQWorkerClient::qualitySummaries_.

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

ecaldqm::TimingClient::~TimingClient ( )

inlineoverride

Definition at line 11 of file TimingClient.h.

{}

Member Function Documentation

void ecaldqm::TimingClient::producePlots ( ProcessType  )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 40 of file TimingClient.cc.

References funct::abs(), ecaldqm::MESetCollection::at(), EcalTrigTowerConstituentsMap::constituentsOf(), EcalBarrel, EcalEndcap, EcalTriggerTower, ecaldqm::MESet::fill(), ecaldqm::MESet::ConstBin::getBinContent(), ecaldqm::MESet::ConstBin::getBinEntries(), ecaldqm::MESet::ConstBin::getBinError(), ecaldqm::DQWorker::getEcalDQMSetupObjects(), ecaldqm::DQWorker::GetElectronicsMap(), ecaldqm::DQWorker::GetTrigTowerMap(), gpuClustering::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, 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_, toleranceMean_, 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(toleranceMean_);
      float rmsThresh(toleranceRMS_);

      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));
      float chStatus(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(toleranceMean_);
      float rmsThresh(toleranceRMS_);

      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);
    }
  }

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 28 of file TimingClient.cc.

References edm::ParameterSet::getUntrackedParameter(), minChannelEntries_, minChannelEntriesFwd_, minTowerEntries_, minTowerEntriesFwd_, tailPopulThreshold_, toleranceMean_, toleranceMeanFwd_, toleranceRMS_, and toleranceRMSFwd_.

                                                             {
    toleranceMean_ = _params.getUntrackedParameter<double>("toleranceMean");
    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

int ecaldqm::TimingClient::minChannelEntries_

private

Definition at line 22 of file TimingClient.h.

Referenced by producePlots(), and setParams().

int ecaldqm::TimingClient::minChannelEntriesFwd_

private

Definition at line 23 of file TimingClient.h.

Referenced by producePlots(), and setParams().

int ecaldqm::TimingClient::minTowerEntries_

private

Definition at line 24 of file TimingClient.h.

Referenced by producePlots(), and setParams().

int ecaldqm::TimingClient::minTowerEntriesFwd_

private

Definition at line 25 of file TimingClient.h.

Referenced by producePlots(), and setParams().

float ecaldqm::TimingClient::tailPopulThreshold_

private

Definition at line 26 of file TimingClient.h.

Referenced by producePlots(), and setParams().

float ecaldqm::TimingClient::toleranceMean_

private

Definition at line 18 of file TimingClient.h.

Referenced by producePlots(), and setParams().

float ecaldqm::TimingClient::toleranceMeanFwd_

private

Definition at line 19 of file TimingClient.h.

Referenced by producePlots(), and setParams().

float ecaldqm::TimingClient::toleranceRMS_

private

Definition at line 20 of file TimingClient.h.

Referenced by producePlots(), and setParams().

float ecaldqm::TimingClient::toleranceRMSFwd_

private

Definition at line 21 of file TimingClient.h.

Referenced by producePlots(), and setParams().