#include <OccupancyClient.h>

Inheritance diagram for ecaldqm::OccupancyClient:

Public Member Functions
	OccupancyClient ()

void	producePlots (ProcessType) override

	~OccupancyClient ()

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

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
float	deviationThreshold_

int	minHits_

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

Constructor & Destructor Documentation

ecaldqm::OccupancyClient::OccupancyClient ( )

Definition at line 13 of file OccupancyClient.cc.

References ecaldqm::DQWorkerClient::qualitySummaries_.

                                    :
     DQWorkerClient(),
     minHits_(0),
     deviationThreshold_(0.)
   {
     qualitySummaries_.insert("QualitySummary");
   }

ecaldqm::OccupancyClient::~OccupancyClient ( )

inline

Definition at line 11 of file OccupancyClient.h.

References producePlots(), and setParams().

11 {}

Member Function Documentation

void ecaldqm::OccupancyClient::producePlots ( ProcessType )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 29 of file OccupancyClient.cc.

References funct::abs(), ecaldqm::dccId(), DEFINE_ECALDQM_WORKER, EcalTrigTowerDetId::detIdFromDenseIndex(), deviationThreshold_, EcalTriggerTower, objects.autophobj::float, ecaldqm::getTrigTowerMap(), triggerObjects_cff::id, photons_cff::ids, EcalTrigTowerDetId::ieta(), ecaldqm::DQWorkerClient::kBad, ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::DQWorkerClient::kGood, ecaldqm::DQWorkerClient::kMBad, ecaldqm::DQWorkerClient::kMGood, EcalTrigTowerDetId::kSizeForDenseIndexing, RecoTauDiscriminantConfiguration::mask, ecaldqm::DQWorker::MEs_, minHits_, ecaldqm::nDCC, EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR, EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_ERROR, EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING, jets_cff::quality, ecaldqm::scConstituents(), ecaldqm::DQWorkerClient::sources_, mathSSE::sqrt(), ecaldqm::DQWorkerClient::statusManager_, electronIdCutBased_cfi::threshold, ecaldqm::MESet::const_iterator::toNextChannel(), ecaldqm::MESet::iterator::toNextChannel(), and EcalTrigTowerConstituentsMap::towerOf().

Referenced by ~OccupancyClient().

   {
     using namespace std;
 
     // number of allowed ieta indices
     // EE-: -28 to -1 with -27, -25 empty
     // EE+: 1 to 28 with 26, 28 empty
     unsigned const nPhiRings(56);
 
     MESet& meQualitySummary(MEs_.at("QualitySummary"));
 //     MESet& meHotDigi(MEs_.at("HotDigi"));
 //     MESet& meHotRecHitThr(MEs_.at("HotRecHitThr"));
 //     MESet& meHotTPDigiThr(MEs_.at("HotTPDigiThr"));
 
     MESet const& sDigi(sources_.at("DigiAll"));
     MESet const& sRecHitThr(sources_.at("RecHitThrAll"));
     MESet const& sTPDigiThr(sources_.at("TPDigiThrAll"));
 
     uint32_t mask(1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR |
                   1 << EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING |
                   1 << EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_ERROR);
 
     double digiPhiRingMean[nPhiRings];
     std::fill_n(digiPhiRingMean, nPhiRings, 0.);
     double rechitPhiRingMean[nPhiRings];
     std::fill_n(rechitPhiRingMean, nPhiRings, 0.);
     int numCrystals[nPhiRings]; // this is static, but is easier to count now
     std::fill_n(numCrystals, nPhiRings, 0);
 
     MESet::const_iterator dEnd(sDigi.end());
     MESet::const_iterator rItr(sRecHitThr);
     for(MESet::const_iterator dItr(sDigi.beginChannel()); dItr != dEnd; dItr.toNextChannel()){
       rItr = dItr;
 
       float entries(dItr->getBinContent());
       float rhentries(rItr->getBinContent());
 
       DetId id(dItr->getId());
       int ieta(0);
       if(id.subdetId() == EcalTriggerTower) // barrel
         ieta = EcalTrigTowerDetId(id).ieta();
       else{
         std::vector<DetId> ids(scConstituents(EcalScDetId(id)));
         if(ids.empty()) continue;
         ieta = getTrigTowerMap()->towerOf(ids[0]).ieta();
       }
 
       unsigned index(ieta < 0 ? ieta + 28 : ieta + 27);
 
       digiPhiRingMean[index] += entries;
       rechitPhiRingMean[index] += rhentries;
       numCrystals[index] += 1;
     }
 
     for(unsigned ie(0); ie < nPhiRings; ie++){
       digiPhiRingMean[ie] /= numCrystals[ie];
       rechitPhiRingMean[ie] /= numCrystals[ie];
     }
 
     // Store # of entries for Occupancy analysis
     std::vector<float> Nentries(nDCC,0.);   // digis
     std::vector<float> Nrhentries(nDCC,0.); // (filtered) rechits
 
     // second round to find hot towers
     for(MESet::const_iterator dItr(sDigi.beginChannel()); dItr != dEnd; dItr.toNextChannel()){
       DetId id(dItr->getId());
 
       bool doMask(meQualitySummary.maskMatches(id, mask, statusManager_));
 
       rItr = dItr;
 
       float entries(dItr->getBinContent());
       float rhentries(rItr->getBinContent());
 
       int ieta(0);
       if(id.subdetId() == EcalTriggerTower) // barrel
         ieta = EcalTrigTowerDetId(id).ieta();
       else{
         std::vector<DetId> ids(scConstituents(EcalScDetId(id)));
         if(ids.empty()) continue;
         ieta = getTrigTowerMap()->towerOf(ids[0]).ieta();
       }
 
       unsigned index(ieta < 0 ? ieta + 28 : ieta + 27);
 
       int quality(doMask ? kMGood : kGood);
 
       if(entries > minHits_ && entries > digiPhiRingMean[index] * deviationThreshold_){
         //        meHotDigi->fill(id);
         quality = doMask ? kMBad : kBad;
       }
       if(rhentries > minHits_ && rhentries > rechitPhiRingMean[index] * deviationThreshold_){
         //        meHotRecHitThr->fill(id);
         quality = doMask ? kMBad : kBad;
       }
 
       meQualitySummary.setBinContent(id, double(quality));
 
       // Keep count of digis & rechits for Occupancy analysis
       unsigned iDCC( dccId(id)-1 );
       if ( entries   > minHits_ ) Nentries[iDCC]   += entries;
       if ( rhentries > minHits_ ) Nrhentries[iDCC] += rhentries;
 
     }
 
     double tpdigiPhiRingMean[nPhiRings];
     std::fill_n(tpdigiPhiRingMean, nPhiRings, 0.);
 
     for(unsigned iTT(0); iTT < EcalTrigTowerDetId::kSizeForDenseIndexing; ++iTT){
       EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(iTT));
       float entries(sTPDigiThr.getBinContent(ttid));
 
       unsigned index(ttid.ieta() < 0 ? ttid.ieta() + 28 : ttid.ieta() + 27);
 
       tpdigiPhiRingMean[index] += entries;
     }
 
     for(int ie(0); ie < 28; ie++){
       float denom(-1.);
       if(ie < 27) denom = 72.;
       else denom = 36.;
       tpdigiPhiRingMean[ie] /= denom;
       tpdigiPhiRingMean[55 - ie] /= denom;
     }
 
     for(unsigned iTT(0); iTT < EcalTrigTowerDetId::kSizeForDenseIndexing; ++iTT){
       EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(iTT));
 
       float entries(sTPDigiThr.getBinContent(ttid));
 
       unsigned index(ttid.ieta() < 0 ? ttid.ieta() + 28 : ttid.ieta() + 27);
 
       int quality(kGood);
 
       if(entries > minHits_ && entries > tpdigiPhiRingMean[index] * deviationThreshold_){
         //        meHotTPDigiThr.fill(ttid);
         quality = kBad;
       }
 
       if(quality != kBad) continue;
 
       std::vector<DetId> ids(getTrigTowerMap()->constituentsOf(ttid));
       for(unsigned iD(0); iD < ids.size(); ++iD){
         DetId& id(ids[iD]);
 
         int quality(meQualitySummary.getBinContent(id));
         if(quality == kMBad || quality == kBad) continue;
 
         meQualitySummary.setBinContent(id, meQualitySummary.maskMatches(id, mask, statusManager_) ? kMBad : kBad);
       }
     }
 
     // Quality check: set entire FED to BAD if its occupancy begins to vanish
     // Fill FED statistics from (filtered) RecHit Occupancy 
     float meanFEDEB(0), meanFEDEE(0), rmsFEDEB(0), rmsFEDEE(0);
     unsigned int nFEDEB(0), nFEDEE(0);
     for ( unsigned iDCC(0); iDCC < nDCC; iDCC++ ) {
       if ( iDCC >=kEBmLow && iDCC <= kEBpHigh) {
         meanFEDEB += Nrhentries[iDCC];
         rmsFEDEB  += Nrhentries[iDCC]*Nrhentries[iDCC];
         nFEDEB++;
       }
       else {
         meanFEDEE += Nrhentries[iDCC];
         rmsFEDEE  += Nrhentries[iDCC]*Nrhentries[iDCC];
         nFEDEE++;
       }
     }
     meanFEDEB /= float( nFEDEB ); rmsFEDEB /= float( nFEDEB );
     meanFEDEE /= float( nFEDEE ); rmsFEDEE /= float( nFEDEE );
     rmsFEDEB   = sqrt( abs(rmsFEDEB - meanFEDEB*meanFEDEB) );
     rmsFEDEE   = sqrt( abs(rmsFEDEE - meanFEDEE*meanFEDEE) );
     // Analyze FED statistics
     float meanFED(0.), rmsFED(0.), nRMS(5.);
     for ( MESet::iterator qsItr(meQualitySummary.beginChannel()); qsItr != meQualitySummary.end(); qsItr.toNextChannel() ) {
       DetId id( qsItr->getId() );
       unsigned iDCC( dccId(id)-1 );
       if ( iDCC >= kEBmLow && iDCC <= kEBpHigh ) {
         meanFED = meanFEDEB;
         rmsFED  = rmsFEDEB;
       }
       else {
         meanFED = meanFEDEE;
         rmsFED  = rmsFEDEE;
       }
       float threshold( meanFED < nRMS*rmsFED ? minHits_ : meanFED - nRMS*rmsFED );
       if ( meanFED > 1000. && Nrhentries[iDCC] < threshold )
         meQualitySummary.setBinContent( id, meQualitySummary.maskMatches(id, mask, statusManager_) ? kMBad : kBad );
     }
 
   } // producePlots()

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 22 of file OccupancyClient.cc.

References deviationThreshold_, edm::ParameterSet::getUntrackedParameter(), and minHits_.

Referenced by ~OccupancyClient().

   {
     minHits_ = _params.getUntrackedParameter<int>("minHits");
     deviationThreshold_ = _params.getUntrackedParameter<double>("deviationThreshold");
   }

Member Data Documentation

float ecaldqm::OccupancyClient::deviationThreshold_

private

Definition at line 19 of file OccupancyClient.h.

Referenced by producePlots(), and setParams().

int ecaldqm::OccupancyClient::minHits_

private

Definition at line 18 of file OccupancyClient.h.

Referenced by producePlots(), and setParams().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation