ecaldqm::TrigPrimClient Class Reference

#include <TrigPrimClient.h>

Inheritance diagram for ecaldqm::TrigPrimClient:

Public Member Functions
void	producePlots (ProcessType) override
	TrigPrimClient ()
	~TrigPrimClient () 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 std::string &	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 () noexcept(false)

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

Private Attributes
float	errorFractionThreshold_
int	minEntries_
bool	sourceFromEmul_
float	TTF4MaskingAlarmThreshold_

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

Constructor & Destructor Documentation

TrigPrimClient()

ecaldqm::TrigPrimClient::TrigPrimClient

(

)

Definition at line 13 of file TrigPrimClient.cc.

      : DQWorkerClient(), minEntries_(0), errorFractionThreshold_(0.), TTF4MaskingAlarmThreshold_(0.) {
    qualitySummaries_.insert("EmulQualitySummary");
  }

References ecaldqm::DQWorkerClient::qualitySummaries_.

~TrigPrimClient()

ecaldqm::TrigPrimClient::~TrigPrimClient ( )

inlineoverride

Definition at line 10 of file TrigPrimClient.h.

{}

Member Function Documentation

producePlots()

void ecaldqm::TrigPrimClient::producePlots ( ProcessType  )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 31 of file TrigPrimClient.cc.

                                               {
    MESet* meNonSingleSummary = nullptr;
    MESet* meTimingSummary = nullptr;
    MESet* sEtEmulError = nullptr;
    MESet* sMatchedIndex = nullptr;
 
    MESet& meEmulQualitySummary(MEs_.at("EmulQualitySummary"));
    MESet& meTrendTTF4Flags(MEs_.at("TrendTTF4Flags"));
 
    MESet const& sTPDigiThrAll(sources_.at("TPDigiThrAll"));
    MESetNonObject const& sLHCStatusByLumi(static_cast<MESetNonObject&>(sources_.at("LHCStatusByLumi")));
 
    uint32_t mask(1 << EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING);
 
    // Store # of entries for Occupancy analysis
    std::vector<float> Nentries(nDCC, 0.);
 
    double currentLHCStatus = sLHCStatusByLumi.getFloatValue();
    bool statsCheckEnabled =
        ((currentLHCStatus > 10.5 && currentLHCStatus < 11.5) ||
         currentLHCStatus < 0);  // currentLHCStatus = -1 is the default when no beam info is available
 
    for (unsigned iTT(0); iTT < EcalTrigTowerDetId::kSizeForDenseIndexing; iTT++) {
      EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(iTT));
 
      bool doMask(meEmulQualitySummary.maskMatches(ttid, mask, statusManager_));
 
      if (sourceFromEmul_) {
        sEtEmulError = &sources_.at("EtEmulError");
        sMatchedIndex = &sources_.at("MatchedIndex");
        meNonSingleSummary = &MEs_.at("NonSingleSummary");
        meTimingSummary = &MEs_.at("TimingSummary");
        float towerEntries(0.);
        float tMax(0.5);
        float nMax(0.);
        for (int iBin(0); iBin < 6; iBin++) {
          float entries(sMatchedIndex->getBinContent(ttid, iBin + 1));
          towerEntries += entries;
 
          if (entries > nMax) {
            nMax = entries;
            tMax = iBin == 0 ? -0.5 : iBin + 0.5;  // historical reasons.. much clearer to say "no entry = -0.5"
          }
        }
        meTimingSummary->setBinContent(ttid, tMax);
        if (towerEntries < minEntries_) {
          meEmulQualitySummary.setBinContent(ttid, doMask ? kMUnknown : kUnknown);
          continue;
        }
 
        float nonsingleFraction(1. - nMax / towerEntries);
 
        if (nonsingleFraction > 0.) {
          meNonSingleSummary->setBinContent(ttid, nonsingleFraction);
        }
 
        if (sEtEmulError->getBinContent(ttid) / towerEntries > errorFractionThreshold_) {
          meEmulQualitySummary.setBinContent(ttid, doMask ? kMBad : kBad);
        } else {
          meEmulQualitySummary.setBinContent(ttid, doMask ? kMGood : kGood);
        }
      }
 
      // Keep count for Occupancy analysis
      unsigned iDCC(dccId(ttid) - 1);
      Nentries[iDCC] += sTPDigiThrAll.getBinContent(ttid);
    }
 
    // Fill TTF4 v Masking ME
    // NOT an occupancy plot: only tells you if non-zero TTF4 occupancy was seen
    // without giving info about how many were seen
    MESet& meTTF4vMask(MEs_.at("TTF4vMask"));
    MESet& meTTF4vMaskByLumi(MEs_.at("TTF4vMaskByLumi"));
    MESet const& sTTFlags4(sources_.at("TTFlags4"));
    MESet const& sTTFlags4ByLumi(sources_.at("TTFlags4ByLumi"));
    MESet const& sTTMaskMapAll(sources_.at("TTMaskMapAll"));
 
    std::vector<float> nWithTTF4(nDCC,
                                 0.);  // counters to keep track of number of towers in a DCC that have TTF4 flag set
    int nWithTTF4_EE = 0;              // total number of towers in EE with TTF4
    int nWithTTF4_EB = 0;              // total number of towers in EB with TTF4
    // Loop over all TTs
    for (unsigned iTT(0); iTT < EcalTrigTowerDetId::kSizeForDenseIndexing; iTT++) {
      EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(iTT));
      unsigned iDCC(dccId(ttid) - 1);
      bool isMasked(sTTMaskMapAll.getBinContent(ttid) > 0.);
      bool hasTTF4(sTTFlags4.getBinContent(ttid) > 0.);
      bool hasTTF4InThisLumiSection(sTTFlags4ByLumi.getBinContent(ttid) > 0.);
      if (hasTTF4InThisLumiSection) {
        nWithTTF4[iDCC]++;
        if (ttid.subDet() == EcalBarrel)
          nWithTTF4_EB++;
        else if (ttid.subDet() == EcalEndcap)
          nWithTTF4_EE++;
      }
      if (isMasked) {
        if (hasTTF4) {
          meTTF4vMask.setBinContent(ttid, 12);  // Masked, has TTF4
        } else {
          meTTF4vMask.setBinContent(ttid, 11);  // Masked, no TTF4
        }
        if (hasTTF4InThisLumiSection) {
          meTTF4vMaskByLumi.setBinContent(ttid, 12);  // Masked, has TTF4
        } else {
          meTTF4vMaskByLumi.setBinContent(ttid, 11);  // Masked, no TTF4
        }
      } else {
        if (hasTTF4)
          meTTF4vMask.setBinContent(ttid, 13);  // not Masked, has TTF4
        if (hasTTF4InThisLumiSection)
          meTTF4vMaskByLumi.setBinContent(ttid, 13);  // not Masked, has TTF4
      }
    }  // TT loop
 
    // Fill trend plots for number of TTs with TTF4 flag set
    meTrendTTF4Flags.fill(EcalBarrel, double(timestamp_.iLumi), nWithTTF4_EB);
    meTrendTTF4Flags.fill(EcalEndcap, double(timestamp_.iLumi), nWithTTF4_EE);
 
    // Quality check: set an entire FED to BAD if a more than 80% of the TTs in that FED show any DCC-SRP flag mismatch errors
    // Fill flag mismatch statistics
    std::vector<float> nTTs(nDCC, 0.);
    std::vector<float> nTTFMismath(nDCC, 0.);
    MESet const& sTTFMismatch(sources_.at("TTFMismatch"));
    for (unsigned iTT(0); iTT < EcalTrigTowerDetId::kSizeForDenseIndexing; iTT++) {
      EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(iTT));
      unsigned iDCC(dccId(ttid) - 1);
      if (sTTFMismatch.getBinContent(ttid) > 0.)
        nTTFMismath[iDCC]++;
      nTTs[iDCC]++;
    }
    // Analyze flag mismatch statistics and TTF4 fraction statistics
    for (unsigned iTT(0); iTT < EcalTrigTowerDetId::kSizeForDenseIndexing; iTT++) {
      EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(iTT));
      unsigned iDCC(dccId(ttid) - 1);
      if (nTTFMismath[iDCC] > 0.8 * nTTs[iDCC] || nWithTTF4[iDCC] > TTF4MaskingAlarmThreshold_ * nTTs[iDCC]) {
        meEmulQualitySummary.setBinContent(ttid,
                                           meEmulQualitySummary.maskMatches(ttid, mask, statusManager_) ? kMBad : kBad);
      }
    }
 
    // Quality check: set entire FED to BAD if its occupancy begins to vanish
    // Fill FED statistics from TP digi 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 += Nentries[iDCC];
        rmsFEDEB += Nentries[iDCC] * Nentries[iDCC];
        nFEDEB++;
      } else {
        meanFEDEE += Nentries[iDCC];
        rmsFEDEE += Nentries[iDCC] * Nentries[iDCC];
        nFEDEE++;
      }
    }
    meanFEDEB /= float(nFEDEB);
    rmsFEDEB /= float(nFEDEB);
    meanFEDEE /= float(nFEDEE);
    rmsFEDEE /= float(nFEDEE);
    rmsFEDEB = sqrt(std::abs(rmsFEDEB - meanFEDEB * meanFEDEB));
    rmsFEDEE = sqrt(std::abs(rmsFEDEE - meanFEDEE * meanFEDEE));
    // Analyze FED statistics
    float meanFED(0.), rmsFED(0.), nRMS(5.);
    for (unsigned iTT(0); iTT < EcalTrigTowerDetId::kSizeForDenseIndexing; iTT++) {
      EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(iTT));
      unsigned iDCC(dccId(ttid) - 1);
      if (iDCC >= kEBmLow && iDCC <= kEBpHigh) {
        meanFED = meanFEDEB;
        rmsFED = rmsFEDEB;
      } else {
        meanFED = meanFEDEE;
        rmsFED = rmsFEDEE;
      }
      float threshold(meanFED < nRMS * rmsFED ? minEntries_ : meanFED - nRMS * rmsFED);
      if ((meanFED > 100. && Nentries[iDCC] < threshold) && statsCheckEnabled)
        meEmulQualitySummary.setBinContent(ttid,
                                           meEmulQualitySummary.maskMatches(ttid, mask, statusManager_) ? kMBad : kBad);
    }
 
  }  // producePlots()

References funct::abs(), ecaldqm::MESetCollection::at(), ecaldqm::dccId(), EcalTrigTowerDetId::detIdFromDenseIndex(), EcalBarrel, EcalEndcap, errorFractionThreshold_, dqmMemoryStats::float, ecaldqm::MESet::getBinContent(), ecaldqm::DQWorker::Timestamp::iLumi, ecaldqm::DQWorkerClient::kBad, ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::DQWorkerClient::kGood, ecaldqm::DQWorkerClient::kMBad, ecaldqm::DQWorkerClient::kMGood, ecaldqm::DQWorkerClient::kMUnknown, EcalTrigTowerDetId::kSizeForDenseIndexing, ecaldqm::DQWorkerClient::kUnknown, ecaldqm::DQWorker::MEs_, minEntries_, ecaldqm::nDCC, metDiagnosticParameterSet_cfi::nMax, EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING, ecaldqm::MESet::setBinContent(), sourceFromEmul_, ecaldqm::DQWorkerClient::sources_, mathSSE::sqrt(), ecaldqm::DQWorkerClient::statusManager_, EcalTrigTowerDetId::subDet(), remoteMonitoring_LED_IterMethod_cfg::threshold, ecaldqm::DQWorker::timestamp_, and TTF4MaskingAlarmThreshold_.

setParams()

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 18 of file TrigPrimClient.cc.

                                                               {
    minEntries_ = _params.getUntrackedParameter<int>("minEntries");
    errorFractionThreshold_ = _params.getUntrackedParameter<double>("errorFractionThreshold");
    TTF4MaskingAlarmThreshold_ = _params.getUntrackedParameter<double>("TTF4MaskingAlarmThreshold");
    sourceFromEmul_ = _params.getUntrackedParameter<bool>("sourceFromEmul");
    if (!sourceFromEmul_) {
      MEs_.erase(std::string("NonSingleSummary"));
      MEs_.erase(std::string("TimingSummary"));
      sources_.erase(std::string("EtEmulError"));
      sources_.erase(std::string("MatchedIndex"));
    }
  }

References l1extraParticles_cfi::_params, ecaldqm::MESetCollection::erase(), errorFractionThreshold_, ecaldqm::DQWorker::MEs_, minEntries_, sourceFromEmul_, ecaldqm::DQWorkerClient::sources_, AlCaHLTBitMon_QueryRunRegistry::string, and TTF4MaskingAlarmThreshold_.

Member Data Documentation

errorFractionThreshold_

float ecaldqm::TrigPrimClient::errorFractionThreshold_

private

Definition at line 18 of file TrigPrimClient.h.

Referenced by producePlots(), and setParams().

minEntries_

int ecaldqm::TrigPrimClient::minEntries_

private

Definition at line 17 of file TrigPrimClient.h.

Referenced by producePlots(), and setParams().

sourceFromEmul_

bool ecaldqm::TrigPrimClient::sourceFromEmul_

private

Definition at line 21 of file TrigPrimClient.h.

Referenced by producePlots(), and setParams().

TTF4MaskingAlarmThreshold_

float ecaldqm::TrigPrimClient::TTF4MaskingAlarmThreshold_

private

Definition at line 19 of file TrigPrimClient.h.

Referenced by producePlots(), and setParams().