#include <TimingClient.h>

Inheritance diagram for ecaldqm::TimingClient:

Public Member Functions
void	producePlots (ProcessType) override

	TimingClient ()

	~TimingClient ()

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 ()

bool	retrieveSource (DQMStore::IGetter &, ProcessType)

bool	runsOn (ProcessType _type) const

void	setStatusManager (StatusManager const &_manager)

virtual	~DQWorkerClient ()

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 &)

std::string const &	getName () const

bool	onlineMode () const

void	setEventNumber (edm::EventNumber_t _e)

void	setLumiNumber (edm::LuminosityBlockNumber_t _l)

void	setRunNumber (edm::RunNumber_t _r)

void	setTime (time_t _t)

virtual	~DQWorker ()

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)

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 ( )

inline

Definition at line 11 of file TimingClient.h.

11 {}

Member Function Documentation

void ecaldqm::TimingClient::producePlots ( ProcessType )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 44 of file TimingClient.cc.

References funct::abs(), EcalTrigTowerConstituentsMap::constituentsOf(), EcalBarrel, EcalTriggerTower, python.tagInventory::entries, ecaldqm::MESet::ConstBin::getBinContent(), ecaldqm::MESet::ConstBin::getBinEntries(), ecaldqm::MESet::ConstBin::getBinError(), ecaldqm::getTrigTowerMap(), ecaldqm::isForward(), ecaldqm::DQWorkerClient::kBad, ecaldqm::DQWorkerClient::kGood, ecaldqm::DQWorkerClient::kMBad, ecaldqm::DQWorkerClient::kMGood, ecaldqm::DQWorkerClient::kMUnknown, ecaldqm::DQWorkerClient::kUnknown, timingPdfMaker::mean, ecaldqm::DQWorker::MEs_, LaserClient_cfi::minChannelEntries, minChannelEntries_, minChannelEntriesFwd_, TimingClient_cfi::minTowerEntries, minTowerEntries_, minTowerEntriesFwd_, EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING, HLT_25ns14e33_v1_cff::quality, plotscripts::rms(), ecaldqm::scConstituents(), ecaldqm::DQWorkerClient::sources_, mathSSE::sqrt(), ecaldqm::DQWorkerClient::statusManager_, EBDetId::switchZSide(), EEDetId::switchZSide(), tailPopulThreshold_, 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"));
 
     uint32_t mask(1 << EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING);
 
     MESet::iterator qEnd(meQuality.end());
 
     MESet::iterator rItr(meRMSMap);
     MESet::const_iterator tItr(sTimeMap);
 
     for(MESet::iterator qItr(meQuality.beginChannel()); qItr != qEnd; qItr.toNextChannel()){
 
       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_));
 
       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(id, mean);
       meMeanAll.fill(id, mean);
       meProjEta.fill(id, mean);
       meProjPhi.fill(id, mean);
       meRMSAll.fill(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(posId);
         }
       }
       else{
         EEDetId eeid(id);
         if(eeid.zside() < 0){
           negative = true;
           EEDetId posId(EEDetId::switchZSide(eeid));
           posTime = sTimeMap.getBinContent(posId);
         }
       }
       if(negative){
         meFwdBkwdDiff.fill(id, posTime - mean);
         meFwdvBkwd.fill(id, mean, posTime);
       }
 
       if(abs(mean) > meanThresh || rms > rmsThresh)
         qItr->setBinContent(doMask ? kMBad : kBad);
       else
         qItr->setBinContent(doMask ? kMGood : kGood);
     }
 
     MESet::iterator qsEnd(meQualitySummary.end());
 
     for(MESet::iterator qsItr(meQualitySummary.beginChannel()); qsItr != qsEnd; qsItr.toNextChannel()){
 
       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(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_);
 
         MESet::const_iterator tmItr(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(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 30 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");
   }