#include <TestPulseClient.h>

Inheritance diagram for ecaldqm::TestPulseClient:

Public Member Functions
void	producePlots (ProcessType) override

	TestPulseClient ()

	~TestPulseClient ()

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

Constructor & Destructor Documentation

ecaldqm::TestPulseClient::TestPulseClient ( )

Definition at line 16 of file TestPulseClient.cc.

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

ecaldqm::TestPulseClient::~TestPulseClient ( )

inline

Definition at line 11 of file TestPulseClient.h.

References producePlots(), and setParams().

11 {}

Member Function Documentation

void ecaldqm::TestPulseClient::producePlots ( ProcessType )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 106 of file TestPulseClient.cc.

References amplitudeThreshold_, ecaldqm::MESetMulti::beginChannel(), DEFINE_ECALDQM_WORKER, EcalBarrel, EcalEndcap, ecaldqm::MESetMulti::end(), gainToME_, ecaldqm::MESetMulti::getBinContent(), ecaldqm::MESetMulti::getBinEntries(), ecaldqm::MESetMulti::getBinError(), hcalTTPDigis_cfi::id, ecaldqm::DQWorkerClient::kBad, ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::DQWorkerClient::kGood, ecaldqm::DQWorkerClient::kMBad, ecaldqm::DQWorkerClient::kMGood, ecaldqm::DQWorkerClient::kMUnknown, ecaldqm::DQWorkerClient::kUnknown, RecoTauDiscriminantConfiguration::mask, ecaldqm::MESetMulti::maskMatches(), ecaldqm::memDCCIndex(), ecaldqm::DQWorker::MEs_, minChannelEntries_, ecaldqm::nDCC, PNAmplitudeThreshold_, pnGainToME_, RecoTauDiscriminantConfiguration::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().

Referenced by ~TestPulseClient().

   {
     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());
       MESet::iterator rItr(meAmplitudeRMS);
       MESet::const_iterator aItr(sAmplitude);
       for(MESet::iterator qItr(meQuality.beginChannel()); qItr != qEnd; qItr.toNextChannel()){
 
         DetId id(qItr->getId());
 
         bool doMask(meQuality.maskMatches(id, mask, statusManager_));
 
         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_));
 
           float amp(sPNAmplitude.getBinContent(id));
           float entries(sPNAmplitude.getBinEntries(id));
           float rms(sPNAmplitude.getBinError(id) * sqrt(entries));
 
           if(entries < minChannelEntries_){
             mePNQualitySummary.setBinContent(id, doMask ? kMUnknown : kUnknown);
             continue;
           }
 
           if(amp < PNAmplitudeThreshold_[gainItr->second] || rms > tolerancePNRMS_[gainItr->second])
             mePNQualitySummary.setBinContent(id, doMask ? kMBad : kBad);
           else
             mePNQualitySummary.setBinContent(id, doMask ? kMGood : kGood);
         }
       }
     }
   }

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 29 of file TestPulseClient.cc.

References CustomPhysics_cfi::amplitude, amplitudeThreshold_, Exception, muonCSCDigis_cfi::gain, gainToME_, ecaldqm::MESetMulti::getIndex(), edm::ParameterSet::getUntrackedParameter(), minChannelEntries_, PNAmplitudeThreshold_, pnGainToME_, ecaldqm::DQWorkerClient::qualitySummaries_, ecaldqm::DQWorkerClient::sources_, tolerancePNRMS_, and toleranceRMS_.

Referenced by ~TestPulseClient().

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