#include <TestPulseClient.h>

Inheritance diagram for ecaldqm::TestPulseClient:

Public Member Functions
void	producePlots (ProcessType) override

	TestPulseClient ()

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

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

const EcalDQMSetupObjects	getEcalDQMSetupObjects ()

const EcalElectronicsMapping *	GetElectronicsMap ()

const CaloGeometry *	GetGeometry ()

const std::string &	getName () const

const CaloTopology *	GetTopology ()

const EcalTrigTowerConstituentsMap *	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	setTime (time_t _t)

virtual	~DQWorker () noexcept(false)

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

Private Attributes
std::vector< float >	amplitudeThreshold_

std::map< int, unsigned >	gainToME_

int	minChannelEntries_

std::vector< float >	PNAmplitudeThreshold_

std::map< int, unsigned >	pnGainToME_

std::vector< float >	tolerancePNRMS_

std::vector< float >	toleranceRMS_

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 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 TestPulseClient.h.

Constructor & Destructor Documentation

◆ TestPulseClient()

ecaldqm::TestPulseClient::TestPulseClient ( )

Definition at line 15 of file TestPulseClient.cc.

       : DQWorkerClient(),
         gainToME_(),
         pnGainToME_(),
         minChannelEntries_(0),
         amplitudeThreshold_(0),
         toleranceRMS_(0),
         PNAmplitudeThreshold_(0),
         tolerancePNRMS_(0) {}

◆ ~TestPulseClient()

ecaldqm::TestPulseClient::~TestPulseClient ( )

inlineoverride

Definition at line 10 of file TestPulseClient.h.

10 {}

Member Function Documentation

◆ producePlots()

void ecaldqm::TestPulseClient::producePlots ( ProcessType )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 108 of file TestPulseClient.cc.

                                                 {
     using namespace std;
  
     MESetMulti& meQuality(static_cast<MESetMulti&>(MEs_.at("Quality")));
     MESetMulti& meAmplitudeRMS(static_cast<MESetMulti&>(MEs_.at("AmplitudeRMS")));
     MESetMulti& meQualitySummary(static_cast<MESetMulti&>(MEs_.at("QualitySummary")));
     MESetMulti& mePNQualitySummary(static_cast<MESetMulti&>(MEs_.at("PNQualitySummary")));
  
     MESetMulti const& sAmplitude(static_cast<MESetMulti const&>(sources_.at("Amplitude")));
     MESetMulti const& sPNAmplitude(static_cast<MESetMulti const&>(sources_.at("PNAmplitude")));
  
     for (map<int, unsigned>::iterator gainItr(gainToME_.begin()); gainItr != gainToME_.end(); ++gainItr) {
       meQuality.use(gainItr->second);
       meQualitySummary.use(gainItr->second);
       meAmplitudeRMS.use(gainItr->second);
  
       sAmplitude.use(gainItr->second);
  
       uint32_t mask(0);
       switch (gainItr->first) {
         case 1:
           mask |= (1 << EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_MEAN_ERROR |
                    1 << EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_RMS_ERROR);
           break;
         case 6:
           mask |= (1 << EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_MEAN_ERROR |
                    1 << EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_RMS_ERROR);
           break;
         case 12:
           mask |= (1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR |
                    1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_RMS_ERROR);
           break;
         default:
           break;
       }
  
       MESet::iterator qEnd(meQuality.end(GetElectronicsMap()));
       MESet::iterator rItr(GetElectronicsMap(), meAmplitudeRMS);
       MESet::const_iterator aItr(GetElectronicsMap(), sAmplitude);
       for (MESet::iterator qItr(meQuality.beginChannel(GetElectronicsMap())); qItr != qEnd;
            qItr.toNextChannel(GetElectronicsMap())) {
         DetId id(qItr->getId());
  
         bool doMask(meQuality.maskMatches(id, mask, statusManager_, GetTrigTowerMap()));
  
         aItr = qItr;
         rItr = qItr;
  
         float entries(aItr->getBinEntries());
  
         if (entries < minChannelEntries_) {
           qItr->setBinContent(doMask ? kMUnknown : kUnknown);
           continue;
         }
  
         float amp(aItr->getBinContent());
         float rms(aItr->getBinError() * sqrt(entries));
  
         rItr->setBinContent(rms);
  
         if (amp < amplitudeThreshold_[gainItr->second] || rms > toleranceRMS_[gainItr->second])
           qItr->setBinContent(doMask ? kMBad : kBad);
         else
           qItr->setBinContent(doMask ? kMGood : kGood);
       }
  
       towerAverage_(meQualitySummary, meQuality, 0.2);
     }
  
     for (map<int, unsigned>::iterator gainItr(pnGainToME_.begin()); gainItr != pnGainToME_.end(); ++gainItr) {
       mePNQualitySummary.use(gainItr->second);
  
       sPNAmplitude.use(gainItr->second);
  
       uint32_t mask(0);
       switch (gainItr->first) {
         case 1:
           mask |= (1 << EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_MEAN_ERROR |
                    1 << EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_RMS_ERROR);
           break;
         case 16:
           mask |= (1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR |
                    1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_RMS_ERROR);
           break;
         default:
           break;
       }
  
       for (unsigned iDCC(0); iDCC < nDCC; ++iDCC) {
         if (memDCCIndex(iDCC + 1) == unsigned(-1))
           continue;
  
         for (unsigned iPN(0); iPN < 10; ++iPN) {
           int subdet(0);
           if (iDCC >= kEBmLow && iDCC <= kEBpHigh)
             subdet = EcalBarrel;
           else
             subdet = EcalEndcap;
  
           EcalPnDiodeDetId id(subdet, iDCC + 1, iPN + 1);
  
           bool doMask(mePNQualitySummary.maskMatches(id, mask, statusManager_, GetTrigTowerMap()));
  
           float amp(sPNAmplitude.getBinContent(getEcalDQMSetupObjects(), id));
           float entries(sPNAmplitude.getBinEntries(getEcalDQMSetupObjects(), id));
           float rms(sPNAmplitude.getBinError(getEcalDQMSetupObjects(), id) * sqrt(entries));
  
           if (entries < minChannelEntries_) {
             mePNQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, doMask ? kMUnknown : kUnknown);
             continue;
           }
  
           if (amp < PNAmplitudeThreshold_[gainItr->second] || rms > tolerancePNRMS_[gainItr->second])
             mePNQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, doMask ? kMBad : kBad);
           else
             mePNQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, doMask ? kMGood : kGood);
         }
       }
     }
   }

References amplitudeThreshold_, ecaldqm::MESetCollection::at(), ecaldqm::MESetMulti::beginChannel(), EcalBarrel, EcalEndcap, ecaldqm::MESetMulti::end(), gainToME_, ecaldqm::MESetMulti::getBinContent(), ecaldqm::MESetMulti::getBinEntries(), ecaldqm::MESetMulti::getBinError(), ecaldqm::DQWorker::getEcalDQMSetupObjects(), ecaldqm::DQWorker::GetElectronicsMap(), ecaldqm::DQWorker::GetTrigTowerMap(), triggerObjects_cff::id, ecaldqm::DQWorkerClient::kBad, ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::DQWorkerClient::kGood, ecaldqm::DQWorkerClient::kMBad, ecaldqm::DQWorkerClient::kMGood, ecaldqm::DQWorkerClient::kMUnknown, ecaldqm::DQWorkerClient::kUnknown, ecaldqm::MESetMulti::maskMatches(), ecaldqm::memDCCIndex(), ecaldqm::DQWorker::MEs_, minChannelEntries_, ecaldqm::nDCC, PNAmplitudeThreshold_, pnGainToME_, SiStripPI::rms, ecaldqm::MESetMulti::setBinContent(), ecaldqm::DQWorkerClient::sources_, mathSSE::sqrt(), ecaldqm::DQWorkerClient::statusManager_, EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR, EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_RMS_ERROR, EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_MEAN_ERROR, EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_RMS_ERROR, EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_MEAN_ERROR, EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_RMS_ERROR, tolerancePNRMS_, toleranceRMS_, ecaldqm::MESet::iterator::toNextChannel(), ecaldqm::DQWorkerClient::towerAverage_(), and ecaldqm::MESetMulti::use().

◆ setParams()

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 25 of file TestPulseClient.cc.

                                                                 {
     minChannelEntries_ = _params.getUntrackedParameter<int>("minChannelEntries");
  
     std::vector<int> MGPAGains(_params.getUntrackedParameter<std::vector<int> >("MGPAGains"));
     std::vector<int> MGPAGainsPN(_params.getUntrackedParameter<std::vector<int> >("MGPAGainsPN"));
  
     MESet::PathReplacements repl;
  
     MESetMulti const& amplitude(static_cast<MESetMulti const&>(sources_.at("Amplitude")));
     unsigned nG(MGPAGains.size());
     for (unsigned iG(0); iG != nG; ++iG) {
       int gain(MGPAGains[iG]);
       if (gain != 1 && gain != 6 && gain != 12)
         throw cms::Exception("InvalidConfiguration") << "MGPA gain";
       repl["gain"] = std::to_string(gain);
       gainToME_[gain] = amplitude.getIndex(repl);
     }
  
     repl.clear();
  
     MESetMulti const& pnAmplitude(static_cast<MESetMulti const&>(sources_.at("PNAmplitude")));
     unsigned nGPN(MGPAGainsPN.size());
     for (unsigned iG(0); iG != nGPN; ++iG) {
       int gain(MGPAGainsPN[iG]);
       if (gain != 1 && gain != 16)
         throw cms::Exception("InvalidConfiguration") << "PN MGPA gain";
       repl["pngain"] = std::to_string(gain);
       pnGainToME_[gain] = pnAmplitude.getIndex(repl);
     }
  
     amplitudeThreshold_.resize(nG);
     toleranceRMS_.resize(nG);
  
     std::vector<double> inAmplitudeThreshold(_params.getUntrackedParameter<std::vector<double> >("amplitudeThreshold"));
     std::vector<double> inToleranceRMS(_params.getUntrackedParameter<std::vector<double> >("toleranceRMS"));
  
     for (std::map<int, unsigned>::iterator gainItr(gainToME_.begin()); gainItr != gainToME_.end(); ++gainItr) {
       unsigned iME(gainItr->second);
       unsigned iGain(0);
       switch (gainItr->first) {
         case 1:
           iGain = 0;
           break;
         case 6:
           iGain = 1;
           break;
         case 12:
           iGain = 2;
           break;
       }
  
       amplitudeThreshold_[iME] = inAmplitudeThreshold[iGain];
       toleranceRMS_[iME] = inToleranceRMS[iGain];
     }
  
     PNAmplitudeThreshold_.resize(nGPN);
     tolerancePNRMS_.resize(nGPN);
  
     std::vector<double> inPNAmplitudeThreshold(
         _params.getUntrackedParameter<std::vector<double> >("PNAmplitudeThreshold"));
     std::vector<double> inTolerancePNRMS(_params.getUntrackedParameter<std::vector<double> >("tolerancePNRMS"));
  
     for (std::map<int, unsigned>::iterator gainItr(pnGainToME_.begin()); gainItr != pnGainToME_.end(); ++gainItr) {
       unsigned iME(gainItr->second);
       unsigned iGain(0);
       switch (gainItr->first) {
         case 1:
           iGain = 0;
           break;
         case 16:
           iGain = 1;
           break;
       }
  
       PNAmplitudeThreshold_[iME] = inPNAmplitudeThreshold[iGain];
       tolerancePNRMS_[iME] = inTolerancePNRMS[iGain];
     }
  
     qualitySummaries_.insert("Quality");
     qualitySummaries_.insert("QualitySummary");
     qualitySummaries_.insert("PNQualitySummary");
   }