Inheritance diagram for GEMDQMHarvester:

Classes
class	NumStatus

struct	PreStatusInfo

Public Types
typedef std::tuple< int, int >	IdChamber

typedef std::tuple< int, int, int >	IdVFAT

typedef struct GEMDQMHarvester::PreStatusInfo	StatusInfo

Public Member Functions
	GEMDQMHarvester (const edm::ParameterSet &)

	~GEMDQMHarvester () override

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Protected Types
typedef std::vector< std::vector< Int_t > >	TableStatusNum

typedef std::vector< std::vector< Float_t > >	TableStatusOcc

Protected Member Functions
Int_t	assessOneBin (std::string strName, Int_t nIdxX, Int_t nIdxY, Float_t fAll, Float_t fNumOcc, Float_t fNumErr, Float_t fNumWarn)

void	copyLabels (MonitorElement h2Src, MonitorElement h2Dst)

void	createInactiveChannelFracHist (edm::Service< DQMStore > &store, std::string strSuffix, Int_t nNumChamber)

void	createLumiFuncHist (edm::Service< DQMStore > &store, std::string strSuffix, Int_t nIdxLayer, Int_t nLumiCurr)

void	createSummaryHist (edm::Service< DQMStore > &store, MonitorElement h2Src, MonitorElement &h2Sum)

void	createSummaryVFAT (edm::Service< DQMStore > &store, MonitorElement h2Src, std::string strSuffix, MonitorElement &h2Sum)

void	createTableWatchingSummary ()

void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override

void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &iLumi, edm::EventSetup const &) override

void	drawSummaryHistogram (edm::Service< DQMStore > &store, Int_t nLumiCurr)

void	getGeometryInfo (edm::Service< DQMStore > &store, MonitorElement *h2Src)

Float_t	refineSummaryHistogram (std::string strName, MonitorElement h2Sum, std::vector< MonitorElement > &listOccPlots, MonitorElement h2SrcStatusA, MonitorElement h2SrcStatusE, MonitorElement h2SrcStatusW, MonitorElement h2SrcStatusEVFAT, MonitorElement h2SrcStatusWVFAT, MonitorElement h2SrcStatusEOH, MonitorElement h2SrcStatusWOH, MonitorElement h2SrcStatusEAMC, MonitorElement h2SrcStatusWAMC, MonitorElement h2SrcStatusEAMC13, Int_t nLumiCurr)

Int_t	refineSummaryVFAT (std::string strName, MonitorElement h2Sum, MonitorElement h2SrcOcc, MonitorElement h2SrcStatusE, MonitorElement h2SrcStatusW, Int_t nLumiCurr, Int_t nIdxLayer)

Int_t	UpdateStatusChamber (Int_t nIdxLayer, Int_t nIdxCh, Int_t nLumiCurr, NumStatus numStatus)

Int_t	UpdateStatusChamber (Int_t nIdxLayer, Int_t nIdxCh, Int_t nIdxVFAT, Int_t nLumiCurr, NumStatus numStatus)

Int_t	UpdateStatusChamber (std::vector< StatusInfo > &listStatus, NumStatus &numStatus, Int_t nLumiCurr, NumStatus numStatusNew)

Protected Attributes
Bool_t	bIsStatusChambersInit_

Float_t	fCutErr_

Float_t	fCutLowErr_

Float_t	fCutWarn_

std::vector< std::string >	listLayer_

std::map< int, MonitorElement * >	mapHistLumiFunc_

std::map< std::string, int >	mapIdxLayer_

std::map< int, int >	mapNumChPerChamber_

std::map< IdChamber, NumStatus >	mapNumStatusChambersSummary_

std::map< IdVFAT, NumStatus >	mapNumStatusVFATsSummary_

std::map< IdChamber, std::vector< StatusInfo > >	mapStatusChambersSummary_

std::map< IdVFAT, std::vector< StatusInfo > >	mapStatusVFATsSummary_

const Int_t	nBitErrAMC13_ = 1

const Int_t	nBitErrAMC_ = 2

const Int_t	nBitErrOH_ = 4

const Int_t	nBitErrVFAT_ = 6

const Int_t	nBitOcc_ = 0

const Int_t	nBitWarnAMC_ = 3

const Int_t	nBitWarnOH_ = 5

const Int_t	nBitWarnVFAT_ = 7

const Int_t	nCodeError_ = 2

const Int_t	nCodeFine_ = 1

const Int_t	nCodeLowError_ = 4

const Int_t	nCodeWarning_ = 3

const Int_t	nMaxLumi_ = 6000

const Int_t	nNumVFATs_ = 24

Int_t	nResolutionLumi_

const std::string	strDirRecHit_ = "GEM/RecHits"

const std::string	strDirStatus_ = "GEM/DAQStatus"

const std::string	strDirSummary_ = "GEM/EventInfo"

Detailed Description

Definition at line 22 of file GEMDQMHarvester.cc.

Member Typedef Documentation

◆ IdChamber

typedef std::tuple<int, int> GEMDQMHarvester::IdChamber

Definition at line 28 of file GEMDQMHarvester.cc.

◆ IdVFAT

typedef std::tuple<int, int, int> GEMDQMHarvester::IdVFAT

Definition at line 29 of file GEMDQMHarvester.cc.

◆ StatusInfo

typedef struct GEMDQMHarvester::PreStatusInfo GEMDQMHarvester::StatusInfo

◆ TableStatusNum

typedef std::vector<std::vector<Int_t> > GEMDQMHarvester::TableStatusNum

protected

Definition at line 143 of file GEMDQMHarvester.cc.

◆ TableStatusOcc

typedef std::vector<std::vector<Float_t> > GEMDQMHarvester::TableStatusOcc

protected

Definition at line 142 of file GEMDQMHarvester.cc.

Constructor & Destructor Documentation

◆ GEMDQMHarvester()

GEMDQMHarvester::GEMDQMHarvester ( const edm::ParameterSet & cfg )

Definition at line 157 of file GEMDQMHarvester.cc.

References looper::cfg.

                                                            {
   fCutErr_ = cfg.getParameter<double>("cutErr");
   fCutLowErr_ = cfg.getParameter<double>("cutLowErr");
   fCutWarn_ = cfg.getParameter<double>("cutWarn");
   nResolutionLumi_ = cfg.getParameter<int>("resolutionLumi");
   bIsStatusChambersInit_ = false;
 }

◆ ~GEMDQMHarvester()

GEMDQMHarvester::~GEMDQMHarvester ( )

inlineoverride

Definition at line 25 of file GEMDQMHarvester.cc.

25 {};

Member Function Documentation

◆ assessOneBin()

Int_t GEMDQMHarvester::assessOneBin	(	std::string	strName,
		Int_t	nIdxX,
		Int_t	nIdxY,
		Float_t	fAll,
		Float_t	fNumOcc,
		Float_t	fNumErr,
		Float_t	fNumWarn
	)

protected

Definition at line 380 of file GEMDQMHarvester.cc.

                                                                                                                    {
   if (fNumErr > fCutErr_ * fAll)  // The error status criterion
     return nCodeError_;
   else if (fNumErr > fCutLowErr_ * fAll)  // The low-error status criterion
     return nCodeLowError_;
   else if (fNumWarn > fCutWarn_ * fAll)  // The warning status criterion
     return nCodeWarning_;
   else if (fNumOcc > 0)
     return nCodeFine_;
 
   return 0;
 }

◆ copyLabels()

void GEMDQMHarvester::copyLabels	(	MonitorElement *	h2Src,
		MonitorElement *	h2Dst
	)

protected

Definition at line 301 of file GEMDQMHarvester.cc.

References dqm::impl::MonitorElement::getAxisTitle(), dqm::impl::MonitorElement::getNbinsX(), dqm::impl::MonitorElement::getNbinsY(), dqm::legacy::MonitorElement::getTH2F(), dqm::impl::MonitorElement::getTitle(), mps_fire::i, anotherprimaryvertexanalyzer_cfi::nBinX, anotherprimaryvertexanalyzer_cfi::nBinY, dqm::impl::MonitorElement::setBinLabel(), dqm::impl::MonitorElement::setTitle(), dqm::impl::MonitorElement::setXTitle(), and dqm::impl::MonitorElement::setYTitle().

                                                                              {
   Int_t nBinX = h2Src->getNbinsX(), nBinY = h2Src->getNbinsY();
 
   for (Int_t i = 1; i <= nBinX; i++) {
     h2Dst->setBinLabel(i, h2Src->getTH2F()->GetXaxis()->GetBinLabel(i), 1);
   }
   for (Int_t i = 1; i <= nBinY; i++) {
     h2Dst->setBinLabel(i, h2Src->getTH2F()->GetYaxis()->GetBinLabel(i), 2);
   }
   h2Dst->setTitle(h2Src->getTitle());
   h2Dst->setXTitle(h2Src->getAxisTitle(1));
   h2Dst->setYTitle(h2Src->getAxisTitle(2));
 }

◆ createInactiveChannelFracHist()

void GEMDQMHarvester::createInactiveChannelFracHist	(	edm::Service< DQMStore > &	store,
		std::string	strSuffix,
		Int_t	nNumChamber
	)

protected

Definition at line 672 of file GEMDQMHarvester.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IGetter::get(), getNameChamberOccGE11(), getNameChamberOccGE21(), getNameChamberOccNull(), dqm::impl::MonitorElement::getNbinsX(), dqm::impl::MonitorElement::getNbinsY(), dqm::legacy::MonitorElement::getTH2F(), mps_fire::i, dqmiolumiharvest::j, dqm::impl::MonitorElement::setBinContent(), dqm::impl::MonitorElement::setBinLabel(), dqm::impl::MonitorElement::setXTitle(), dqm::impl::MonitorElement::setYTitle(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                        {
   std::string strTitle = "The fraction of inactive channels in " + strSuffix.substr(1);
   MonitorElement *h2InactiveChannel =
       store->book1D("inactive_frac_chamber" + strSuffix, strTitle, nNumChamber, 0.5, nNumChamber + 0.5);
   h2InactiveChannel->setXTitle("Chamber");
   h2InactiveChannel->setYTitle("Fraction of inactive channels");
   for (Int_t i = 1; i <= nNumChamber; i++) {
     h2InactiveChannel->setBinLabel(i, Form("%i", i), 1);
   }
 
   std::string (*funcNameCh)(std::string, Int_t) = getNameChamberOccNull;
 
   if (strSuffix.find("_GE11") != std::string::npos) {
     funcNameCh = getNameChamberOccGE11;
   } else if (strSuffix.find("_GE21") != std::string::npos) {
     funcNameCh = getNameChamberOccGE21;
   }
 
   for (Int_t nIdxCh = 1; nIdxCh <= nNumChamber; nIdxCh++) {
     std::string strNameCh = funcNameCh(strSuffix, nIdxCh);
     MonitorElement *h2SrcChamberOcc = store->get(strNameCh);
     if (h2SrcChamberOcc == nullptr) {
       // FIXME: It's about sending a message
       continue;
     }
 
     Int_t nNumBinX = h2SrcChamberOcc->getNbinsX();
     Int_t nNumBinY = h2SrcChamberOcc->getNbinsY();
     Int_t nNumAllChannel = nNumBinX * nNumBinY;
     auto *histData = h2SrcChamberOcc->getTH2F();
     auto *pdData = histData->GetArray();
     Int_t nNumChannelInactive = 0;
     for (Int_t j = 1; j <= nNumBinY; j++)
       for (Int_t i = 1; i <= nNumBinX; i++) {
         if (pdData[j * (nNumBinX + 2) + i] <= 0) {
           nNumChannelInactive++;
         }
       }
     h2InactiveChannel->setBinContent(nIdxCh, ((Double_t)nNumChannelInactive) / nNumAllChannel);
   }
 }

◆ createLumiFuncHist()

void GEMDQMHarvester::createLumiFuncHist	(	edm::Service< DQMStore > &	store,
		std::string	strSuffix,
		Int_t	nIdxLayer,
		Int_t	nLumiCurr
	)

protected

Definition at line 558 of file GEMDQMHarvester.cc.

References dqm::implementation::IBooker::book2S(), mps_fire::i, dqm::impl::MonitorElement::setBinContent(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::DQMStore::setCurrentFolder(), dqm::impl::MonitorElement::setXTitle(), and dqm::impl::MonitorElement::setYTitle().

                                                           {
   auto &nNumCh = mapNumChPerChamber_[nIdxLayer];
 
   MonitorElement *h2Summary;
 
   //Int_t nLumiCurrLowRes = ( ( nLumiCurr - 1 ) / nResolutionLumi_ ) * nResolutionLumi_;
   Int_t nNumBinLumi = ((nLumiCurr - 1) / nResolutionLumi_) + 1;
   Int_t nMaxBin = 0;
 
   // Creating or Summoning the corresponding histogram
   if (mapHistLumiFunc_.find(nIdxLayer) == mapHistLumiFunc_.end()) {
     store->setCurrentFolder(strDirSummary_);
     h2Summary = store->book2S("chamberStatus_inLumi" + strSuffix,
                               "Chamber status on lumi-block " + strSuffix.substr(1),
                               nMaxLumi_ / nResolutionLumi_,
                               1.0,
                               (Float_t)(nMaxLumi_ + 1),
                               //nNumBinLumi, 1.0, (Float_t)( nLumiCurr + 1 ),
                               nNumCh,
                               0.5,
                               nNumCh + 0.5);
     mapHistLumiFunc_[nIdxLayer] = h2Summary;
 
     h2Summary->setXTitle("Luminosity block");
     h2Summary->setYTitle("Chamber");
     for (Int_t i = 1; i <= nNumCh; i++) {
       h2Summary->setBinLabel(i, Form("%i", i), 2);
     }
   } else {
     h2Summary = mapHistLumiFunc_[nIdxLayer];
   }
 
   for (Int_t nIdxCh = 1; nIdxCh <= nNumCh; nIdxCh++) {
     auto &listStatus = mapStatusChambersSummary_[{nIdxLayer, nIdxCh}];
 
     Int_t nIdxStatus = 0;
     for (Int_t nIdxLumi = 0; nIdxLumi < nNumBinLumi; nIdxLumi++) {
       // Lumis covered by these values (nLumiStart <=, <= nLumiEnd) are counted for the current bin
       Int_t nLumiStart = 1 + nIdxLumi * nResolutionLumi_;
       Int_t nLumiEnd = (1 + nIdxLumi) * nResolutionLumi_;
       if (nLumiEnd > nLumiCurr)
         nLumiEnd = nLumiCurr;
 
       Int_t nStatusSum = 0;
       while (true) {  // No worries, nIdxStatus must increase and reach at listStatus.size()
         // True: It was too past so that
         //       the lumi range of listStatus[ nIdxStatus ] is out of the coverage of the current bin
         if (listStatus[nIdxStatus].nLumiEnd < nLumiStart) {
           nIdxStatus++;
           if (nIdxStatus >= (int)listStatus.size()) {
             break;  // For safety
           }
           continue;
         }
 
         nStatusSum = listStatus[nIdxStatus].nStatus;
 
         // True: This is the last item of listStatus which is covered by the current bin
         if (nIdxStatus + 1 >= (int)listStatus.size() || listStatus[nIdxStatus].nLumiEnd >= nLumiEnd) {
           break;
         }
 
         nIdxStatus++;
         if (nIdxStatus >= (int)listStatus.size()) {
           break;  // For safety
         }
       }
 
       nStatusSum &= ~(1 << nBitOcc_);  // No need of displaying the digi occupancy
       h2Summary->setBinContent(nIdxLumi + 1, nIdxCh, nStatusSum);
       if (nMaxBin < nIdxLumi + 1)
         nMaxBin = nIdxLumi + 1;
     }
   }
 
   for (Int_t nX = 1; nX <= nMaxBin; nX++) {
     h2Summary->setBinContent(nX, 0, 1);
   }
 }

◆ createSummaryHist()

void GEMDQMHarvester::createSummaryHist	(	edm::Service< DQMStore > &	store,
		MonitorElement *	h2Src,
		MonitorElement *&	h2Sum
	)

protected

Definition at line 353 of file GEMDQMHarvester.cc.

References dqm::implementation::IBooker::book2D(), dqm::impl::MonitorElement::getNbinsX(), dqm::impl::MonitorElement::getNbinsY(), dqm::legacy::MonitorElement::getTH2F(), mps_fire::i, anotherprimaryvertexanalyzer_cfi::nBinX, anotherprimaryvertexanalyzer_cfi::nBinY, dqm::impl::MonitorElement::setBinLabel(), dqm::impl::MonitorElement::setTitle(), dqm::impl::MonitorElement::setXTitle(), and dqm::impl::MonitorElement::setYTitle().

                                                                                                                   {
   //store->setCurrentFolder(strDirSummary_);
 
   Int_t nBinX = h2Src->getNbinsX(), nBinY = h2Src->getNbinsY();
   h2Sum = store->book2D("reportSummaryMap", "", nBinX, 0.5, nBinX + 0.5, nBinY, 0.5, nBinY + 0.5);
   h2Sum->setTitle("Summary plot");
   h2Sum->setXTitle("Chamber");
   h2Sum->setYTitle("Layer");
 
   for (Int_t i = 1; i <= nBinX; i++)
     h2Sum->setBinLabel(i, h2Src->getTH2F()->GetXaxis()->GetBinLabel(i), 1);
   for (Int_t i = 1; i <= nBinY; i++)
     h2Sum->setBinLabel(i, listLayer_[i].substr(1), 2);
 }

◆ createSummaryVFAT()

void GEMDQMHarvester::createSummaryVFAT	(	edm::Service< DQMStore > &	store,
		MonitorElement *	h2Src,
		std::string	strSuffix,
		MonitorElement *&	h2Sum
	)

protected

Definition at line 368 of file GEMDQMHarvester.cc.

References dqm::implementation::IBooker::book2D(), dqm::impl::MonitorElement::getNbinsX(), dqm::impl::MonitorElement::getNbinsY(), anotherprimaryvertexanalyzer_cfi::nBinX, and anotherprimaryvertexanalyzer_cfi::nBinY.

                                                                 {
   //store->setCurrentFolder(strDirStatus_);
   //store->setCurrentFolder(strDirSummary_);
 
   Int_t nBinX = h2Src->getNbinsX(), nBinY = h2Src->getNbinsY();
   h2Sum = store->book2D("vfat_statusSummary" + strSuffix, "", nBinX, 0.5, nBinX + 0.5, nBinY, -0.5, nBinY - 0.5);
   copyLabels(h2Src, h2Sum);
 }

◆ createTableWatchingSummary()

void GEMDQMHarvester::createTableWatchingSummary ( )

protected

Definition at line 283 of file GEMDQMHarvester.cc.

References mps_fire::i, and dqmiolumiharvest::j.

                                                  {
   if (bIsStatusChambersInit_)
     return;
 
   for (const auto &[nIdxLayer, nNumCh] : mapNumChPerChamber_) {
     for (Int_t i = 1; i <= nNumCh; i++) {
       mapStatusChambersSummary_[{nIdxLayer, i}] = std::vector<StatusInfo>();
       mapNumStatusChambersSummary_[{nIdxLayer, i}] = NumStatus();
       for (Int_t j = 1; j <= nNumVFATs_; j++) {
         mapStatusVFATsSummary_[{nIdxLayer, i, j}] = std::vector<StatusInfo>();
         mapNumStatusVFATsSummary_[{nIdxLayer, i, j}] = NumStatus();
       }
     }
   }
 
   bIsStatusChambersInit_ = true;
 }

◆ dqmEndJob()

void GEMDQMHarvester::dqmEndJob	(	DQMStore::IBooker &	,
		DQMStore::IGetter &
	)

inlineoverrideprotected

Definition at line 72 of file GEMDQMHarvester.cc.

72 {}; // Cannot use; it is called after dqmSaver

◆ dqmEndLuminosityBlock()

void GEMDQMHarvester::dqmEndLuminosityBlock	(	DQMStore::IBooker &	,
		DQMStore::IGetter &	,
		edm::LuminosityBlock const &	iLumi,
		edm::EventSetup const &
	)

overrideprotected

Definition at line 174 of file GEMDQMHarvester.cc.

References edm::LuminosityBlockBase::id(), and edm::LuminosityBlockID::luminosityBlock().

                                                                    {
   edm::Service<DQMStore> store;
   Int_t nLumiCurr = iLumi.id().luminosityBlock();
   drawSummaryHistogram(store, nLumiCurr);
 }

◆ drawSummaryHistogram()

void GEMDQMHarvester::drawSummaryHistogram	(	edm::Service< DQMStore > &	store,
		Int_t	nLumiCurr
	)

protected

Definition at line 183 of file GEMDQMHarvester.cc.

References dqm::implementation::IBooker::bookFloat(), dqm::impl::MonitorElement::Fill(), dqm::implementation::IGetter::get(), dqm::impl::MonitorElement::getAxisTitle(), dqm::impl::MonitorElement::getTitle(), dqm::implementation::DQMStore::setCurrentFolder(), dqm::impl::MonitorElement::setTitle(), dqm::impl::MonitorElement::setXTitle(), dqm::impl::MonitorElement::setYTitle(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                        {
   Float_t fReportSummary = -1.0;
 
   std::string strSrcStatusA = "GEM/DAQStatus/chamberAllStatus";
   std::string strSrcStatusE = "GEM/DAQStatus/chamberErrors";
   std::string strSrcStatusW = "GEM/DAQStatus/chamberWarnings";
   std::string strSrcStatusEVFAT = "GEM/DAQStatus/chamberVFATErrors";
   std::string strSrcStatusWVFAT = "GEM/DAQStatus/chamberVFATWarnings";
   std::string strSrcStatusEOH = "GEM/DAQStatus/chamberOHErrors";
   std::string strSrcStatusWOH = "GEM/DAQStatus/chamberOHWarnings";
   std::string strSrcStatusEAMC = "GEM/DAQStatus/chamberAMCErrors";
   std::string strSrcStatusWAMC = "GEM/DAQStatus/chamberAMCWarnings";
   std::string strSrcStatusEAMC13 = "GEM/DAQStatus/chamberAMC13Errors";
 
   std::string strSrcVFATOcc = "GEM/Digis/occ";
   std::string strSrcVFATStatusW = "GEM/DAQStatus/vfat_statusWarnSum";
   std::string strSrcVFATStatusE = "GEM/DAQStatus/vfat_statusErrSum";
 
   store->setCurrentFolder(strDirSummary_);
 
   MonitorElement *h2SrcStatusA = store->get(strSrcStatusA);
   MonitorElement *h2SrcStatusE = store->get(strSrcStatusE);
   MonitorElement *h2SrcStatusW = store->get(strSrcStatusW);
   MonitorElement *h2SrcStatusEVFAT = store->get(strSrcStatusEVFAT);
   MonitorElement *h2SrcStatusWVFAT = store->get(strSrcStatusWVFAT);
   MonitorElement *h2SrcStatusEOH = store->get(strSrcStatusEOH);
   MonitorElement *h2SrcStatusWOH = store->get(strSrcStatusWOH);
   MonitorElement *h2SrcStatusEAMC = store->get(strSrcStatusEAMC);
   MonitorElement *h2SrcStatusWAMC = store->get(strSrcStatusWAMC);
   MonitorElement *h2SrcStatusEAMC13 = store->get(strSrcStatusEAMC13);
 
   std::string strTitleSummary = "summary";
 
   getGeometryInfo(store, h2SrcStatusEOH);
 
   if (h2SrcStatusA != nullptr && h2SrcStatusE != nullptr && h2SrcStatusW != nullptr && h2SrcStatusEVFAT != nullptr &&
       h2SrcStatusWVFAT != nullptr && h2SrcStatusEOH != nullptr && h2SrcStatusWOH != nullptr &&
       h2SrcStatusEAMC != nullptr && h2SrcStatusWAMC != nullptr && h2SrcStatusEAMC13 != nullptr) {
     MonitorElement *h2Sum = nullptr;
     createSummaryHist(store, h2SrcStatusEOH, h2Sum);
     createTableWatchingSummary();
 
     std::vector<MonitorElement *> listOccPlots(listLayer_.size() + 1);  // The index starts at 1
     for (const auto &strSuffix : listLayer_) {
       if (mapIdxLayer_.find(strSuffix) == mapIdxLayer_.end())
         continue;
       auto nIdxLayer = mapIdxLayer_[strSuffix];
       MonitorElement *h2SrcVFATOcc = store->get(strSrcVFATOcc + strSuffix);
       if (h2SrcVFATOcc == nullptr)
         continue;
       listOccPlots[nIdxLayer] = h2SrcVFATOcc;
     }
 
     fReportSummary = refineSummaryHistogram(strTitleSummary,
                                             h2Sum,
                                             listOccPlots,
                                             h2SrcStatusA,
                                             h2SrcStatusE,
                                             h2SrcStatusW,
                                             h2SrcStatusEVFAT,
                                             h2SrcStatusWVFAT,
                                             h2SrcStatusEOH,
                                             h2SrcStatusWOH,
                                             h2SrcStatusEAMC,
                                             h2SrcStatusWAMC,
                                             h2SrcStatusEAMC13,
                                             nLumiCurr);
 
     for (const auto &strSuffix : listLayer_) {
       if (mapIdxLayer_.find(strSuffix) == mapIdxLayer_.end())
         continue;
       auto nIdxLayer = mapIdxLayer_[strSuffix];
       MonitorElement *h2SrcVFATOcc = store->get(strSrcVFATOcc + strSuffix);
       MonitorElement *h2SrcVFATStatusW = store->get(strSrcVFATStatusW + strSuffix);
       MonitorElement *h2SrcVFATStatusE = store->get(strSrcVFATStatusE + strSuffix);
       if (h2SrcVFATOcc == nullptr || h2SrcVFATStatusW == nullptr || h2SrcVFATStatusE == nullptr)
         continue;
 
       MonitorElement *h2SumVFAT = nullptr;
       createSummaryVFAT(store, h2SrcVFATStatusE, strSuffix, h2SumVFAT);
       refineSummaryVFAT(strSuffix, h2SumVFAT, h2SrcVFATOcc, h2SrcVFATStatusE, h2SrcVFATStatusW, nLumiCurr, nIdxLayer);
       TString strNewTitle = h2SrcVFATStatusE->getTitle();
       h2SumVFAT->setTitle((const char *)strNewTitle.ReplaceAll("errors", "errors/warnings"));
       h2SumVFAT->setXTitle(h2SrcVFATStatusE->getAxisTitle(1));
       h2SumVFAT->setYTitle(h2SrcVFATStatusE->getAxisTitle(2));
 
       createLumiFuncHist(store, strSuffix, nIdxLayer, nLumiCurr);
     }
   }
 
   for (const auto &strSuffix : listLayer_) {
     if (mapIdxLayer_.find(strSuffix) == mapIdxLayer_.end())
       continue;
     auto nNumChamber = mapNumChPerChamber_[mapIdxLayer_[strSuffix]];
     createInactiveChannelFracHist(store, strSuffix, nNumChamber);
   }
 
   store->bookFloat("reportSummary")->Fill(fReportSummary);
 }

◆ fillDescriptions()

void GEMDQMHarvester::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 165 of file GEMDQMHarvester.cc.

References edm::ConfigurationDescriptions::add(), and submitPVResolutionJobs::desc.

                                                                                  {
   edm::ParameterSetDescription desc;
   desc.add<double>("cutErr", 0.05);
   desc.add<double>("cutLowErr", 0.00);
   desc.add<double>("cutWarn", 0.05);
   desc.add<int>("resolutionLumi", 1);
   descriptions.add("GEMDQMHarvester", desc);
 }

◆ getGeometryInfo()

void GEMDQMHarvester::getGeometryInfo	(	edm::Service< DQMStore > &	store,
		MonitorElement *	h2Src
	)

protected

Definition at line 315 of file GEMDQMHarvester.cc.

References dqm::implementation::IGetter::get(), dqm::impl::MonitorElement::getBinContent(), dqm::impl::MonitorElement::getNbinsY(), dqm::legacy::MonitorElement::getTH2F(), mps_fire::i, anotherprimaryvertexanalyzer_cfi::nBinY, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                         {
   listLayer_.clear();
   mapIdxLayer_.clear();
   mapNumChPerChamber_.clear();
 
   if (h2Src != nullptr) {  // For online and offline
     Int_t nBinY = h2Src->getNbinsY();
     listLayer_.push_back("");
 
     for (Int_t i = 1; i <= nBinY; i++) {
       std::string strLabelFull = h2Src->getTH2F()->GetYaxis()->GetBinLabel(i);
       Int_t nBinXActual = (Int_t)(h2Src->getBinContent(0, i) + 0.5);
       auto nPos = strLabelFull.find(';');
       auto strLayer = strLabelFull.substr(nPos + 1);
       listLayer_.push_back(strLayer);
       mapIdxLayer_[strLayer] = i;
       mapNumChPerChamber_[i] = nBinXActual;
     }
   } else {  // For others (validation and...?)
     listLayer_.push_back("");
     if (store->get("GEM/Digis/occupancy_GE11-M-L1/occ_GE11-M-01L1-S") != nullptr) {
       listLayer_.push_back("_GE11-P-L2");
       listLayer_.push_back("_GE11-P-L1");
       listLayer_.push_back("_GE11-M-L1");
       listLayer_.push_back("_GE11-M-L2");
       mapIdxLayer_["_GE11-P-L2"] = 1;
       mapIdxLayer_["_GE11-P-L1"] = 2;
       mapIdxLayer_["_GE11-M-L1"] = 3;
       mapIdxLayer_["_GE11-M-L2"] = 4;
       mapNumChPerChamber_[1] = 36;
       mapNumChPerChamber_[2] = 36;
       mapNumChPerChamber_[3] = 36;
       mapNumChPerChamber_[4] = 36;
     }
     // FIXME: How about GE21 and ME0?
   }
 }

◆ refineSummaryHistogram()

Float_t GEMDQMHarvester::refineSummaryHistogram	(	std::string	strName,
		MonitorElement *	h2Sum,
		std::vector< MonitorElement *> &	listOccPlots,
		MonitorElement *	h2SrcStatusA,
		MonitorElement *	h2SrcStatusE,
		MonitorElement *	h2SrcStatusW,
		MonitorElement *	h2SrcStatusEVFAT,
		MonitorElement *	h2SrcStatusWVFAT,
		MonitorElement *	h2SrcStatusEOH,
		MonitorElement *	h2SrcStatusWOH,
		MonitorElement *	h2SrcStatusEAMC,
		MonitorElement *	h2SrcStatusWAMC,
		MonitorElement *	h2SrcStatusEAMC13,
		Int_t	nLumiCurr
	)

protected

Definition at line 395 of file GEMDQMHarvester.cc.

References dqm::impl::MonitorElement::getBinContent(), dqm::impl::MonitorElement::getNbinsY(), mps_fire::i, dqmiolumiharvest::j, anotherprimaryvertexanalyzer_cfi::nBinX, anotherprimaryvertexanalyzer_cfi::nBinY, alignCSCRings::r, and dqm::impl::MonitorElement::setBinContent().

                                                                  {
   Int_t nBinY = h2Sum->getNbinsY();
   Int_t nAllBin = 0, nFineBin = 0;
   for (Int_t j = 1; j <= nBinY; j++) {
     Int_t nBinX = (Int_t)(h2SrcStatusE->getBinContent(0, j) + 0.5);
     auto h2SrcOcc = listOccPlots[j];
     Int_t nBinYOcc = 0;
     if (h2SrcOcc != nullptr) {
       nBinYOcc = h2SrcOcc->getNbinsY();
     }
 
     h2Sum->setBinContent(0, j, nBinX);
     for (Int_t i = 1; i <= nBinX; i++) {
       Float_t fOcc = 0;
       for (Int_t r = 1; r <= nBinYOcc; r++) {
         fOcc += h2SrcOcc->getBinContent(i, r);
       }
 
       Float_t fStatusAll = h2SrcStatusA->getBinContent(i, j);
       Float_t fStatusErr = h2SrcStatusE->getBinContent(i, j);
       Float_t fStatusWarn = h2SrcStatusW->getBinContent(i, j);
       Float_t fStatusErrVFAT = h2SrcStatusEVFAT->getBinContent(i, j);
       Float_t fStatusWarnVFAT = h2SrcStatusWVFAT->getBinContent(i, j);
       Float_t fStatusErrOH = h2SrcStatusEOH->getBinContent(i, j);
       Float_t fStatusWarnOH = h2SrcStatusWOH->getBinContent(i, j);
       Float_t fStatusErrAMC = h2SrcStatusEAMC->getBinContent(i, j);
       Float_t fStatusWarnAMC = h2SrcStatusWAMC->getBinContent(i, j);
       Float_t fStatusErrAMC13 = h2SrcStatusEAMC13->getBinContent(i, j);
       NumStatus numStatus(fStatusAll,
                           fOcc,
                           fStatusErrVFAT,
                           fStatusWarnVFAT,
                           fStatusErrOH,
                           fStatusWarnOH,
                           fStatusErrAMC,
                           fStatusWarnAMC,
                           fStatusErrAMC13);
       UpdateStatusChamber(j, i, nLumiCurr, numStatus);
 
       Int_t nRes = assessOneBin(strName, i, j, fStatusAll, fOcc, fStatusErr, fStatusWarn);
       if (nRes == 1)
         nFineBin++;
 
       h2Sum->setBinContent(i, j, (Float_t)nRes);
       nAllBin++;
     }
   }
 
   return ((Float_t)nFineBin) / nAllBin;
 }

◆ refineSummaryVFAT()

Int_t GEMDQMHarvester::refineSummaryVFAT	(	std::string	strName,
		MonitorElement *	h2Sum,
		MonitorElement *	h2SrcOcc,
		MonitorElement *	h2SrcStatusE,
		MonitorElement *	h2SrcStatusW,
		Int_t	nLumiCurr,
		Int_t	nIdxLayer
	)

protected

Definition at line 459 of file GEMDQMHarvester.cc.

References dqm::impl::MonitorElement::getBinContent(), dqm::impl::MonitorElement::getNbinsX(), dqm::impl::MonitorElement::getNbinsY(), mps_fire::i, dqmiolumiharvest::j, anotherprimaryvertexanalyzer_cfi::nBinX, anotherprimaryvertexanalyzer_cfi::nBinY, and dqm::impl::MonitorElement::setBinContent().

                                                           {
   Int_t nBinY = h2Sum->getNbinsY();
   for (Int_t j = 1; j <= nBinY; j++) {
     Int_t nBinX = h2Sum->getNbinsX();
     for (Int_t i = 1; i <= nBinX; i++) {
       Float_t fOcc = h2SrcOcc->getBinContent(i, j);
       Float_t fStatusErr = h2SrcStatusE->getBinContent(i, j);
       Float_t fStatusWarn = h2SrcStatusW->getBinContent(i, j);
       Float_t fStatusAll = fOcc + fStatusErr + fStatusWarn;
       NumStatus numStatus(fStatusAll, fOcc, fStatusErr, fStatusWarn, 0, 0, 0, 0, 0);
       UpdateStatusChamber(nIdxLayer, i, j, nLumiCurr, numStatus);
 
       Int_t nRes = assessOneBin(strName, i, j, fStatusAll, fOcc, fStatusErr, fStatusWarn);
       h2Sum->setBinContent(i, j, (Float_t)nRes);
     }
   }
 
   return 0;
 }

◆ UpdateStatusChamber() [1/3]

Int_t GEMDQMHarvester::UpdateStatusChamber	(	Int_t	nIdxLayer,
		Int_t	nIdxCh,
		Int_t	nLumiCurr,
		NumStatus	numStatus
	)

protected

Definition at line 485 of file GEMDQMHarvester.cc.

                                                                                                               {
   if (!bIsStatusChambersInit_)
     return 0;
   if (0 >= nIdxCh || nIdxCh > mapNumChPerChamber_[nIdxLayer])
     return 0;
   auto &listStatus = mapStatusChambersSummary_[{nIdxLayer, nIdxCh}];
   auto &numStatusPrev = mapNumStatusChambersSummary_[{nIdxLayer, nIdxCh}];
   return UpdateStatusChamber(listStatus, numStatusPrev, nLumiCurr, numStatus);
 }

◆ UpdateStatusChamber() [2/3]

Int_t GEMDQMHarvester::UpdateStatusChamber	(	Int_t	nIdxLayer,
		Int_t	nIdxCh,
		Int_t	nIdxVFAT,
		Int_t	nLumiCurr,
		NumStatus	numStatus
	)

protected

Definition at line 495 of file GEMDQMHarvester.cc.

                                                                                          {
   if (!bIsStatusChambersInit_)
     return 0;
   if (0 >= nIdxCh || nIdxCh > mapNumChPerChamber_[nIdxLayer])
     return 0;
   if (0 >= nIdxVFAT || nIdxVFAT > nNumVFATs_)
     return 0;
   auto &listStatus = mapStatusVFATsSummary_[{nIdxLayer, nIdxCh, nIdxVFAT}];
   auto &numStatusPrev = mapNumStatusVFATsSummary_[{nIdxLayer, nIdxCh, nIdxVFAT}];
   return UpdateStatusChamber(listStatus, numStatusPrev, nLumiCurr, numStatus);
 }

◆ UpdateStatusChamber() [3/3]

Int_t GEMDQMHarvester::UpdateStatusChamber	(	std::vector< StatusInfo > &	listStatus,
		NumStatus &	numStatus,
		Int_t	nLumiCurr,
		NumStatus	numStatusNew
	)

protected

Definition at line 508 of file GEMDQMHarvester.cc.

                                                                    {
   // First of all, the current lumi section will be assessed, of which the result will be stored in nStatus
   Int_t nStatus = 0;
 
   Float_t fNumAddErrVFAT = numStatusNew.fNumErrVFAT_ - numStatus.fNumErrVFAT_;
   Float_t fNumAddWarnVFAT = numStatusNew.fNumWarnVFAT_ - numStatus.fNumWarnVFAT_;
   Float_t fNumAddErrOH = numStatusNew.fNumErrOH_ - numStatus.fNumErrOH_;
   Float_t fNumAddWarnOH = numStatusNew.fNumWarnOH_ - numStatus.fNumWarnOH_;
   Float_t fNumAddErrAMC = numStatusNew.fNumErrAMC_ - numStatus.fNumErrAMC_;
   Float_t fNumAddWarnAMC = numStatusNew.fNumWarnAMC_ - numStatus.fNumWarnAMC_;
   Float_t fNumAddErrAMC13 = numStatusNew.fNumErrAMC13_ - numStatus.fNumErrAMC13_;
 
   numStatus.fNumTotal_ = numStatusNew.fNumTotal_;
   numStatus.fNumOcc_ = numStatusNew.fNumOcc_;
   numStatus.fNumErrVFAT_ = numStatusNew.fNumErrVFAT_;
   numStatus.fNumWarnVFAT_ = numStatusNew.fNumWarnVFAT_;
   numStatus.fNumErrOH_ = numStatusNew.fNumErrOH_;
   numStatus.fNumWarnOH_ = numStatusNew.fNumWarnOH_;
   numStatus.fNumErrAMC_ = numStatusNew.fNumErrAMC_;
   numStatus.fNumWarnAMC_ = numStatusNew.fNumWarnAMC_;
   numStatus.fNumErrAMC13_ = numStatusNew.fNumErrAMC13_;
 
   nStatus = (numStatusNew.fNumOcc_ > 0 ? 1 << nBitOcc_ : 0) | (fNumAddErrAMC13 > 0 ? 1 << nBitErrAMC13_ : 0) |
             (fNumAddErrAMC > 0 ? 1 << nBitErrAMC_ : 0) | (fNumAddWarnAMC > 0 ? 1 << nBitWarnAMC_ : 0) |
             (fNumAddErrOH > 0 ? 1 << nBitErrOH_ : 0) | (fNumAddWarnOH > 0 ? 1 << nBitWarnOH_ : 0) |
             (fNumAddErrVFAT > 0 ? 1 << nBitErrVFAT_ : 0) | (fNumAddWarnVFAT > 0 ? 1 << nBitWarnVFAT_ : 0);
 
   // Only used in the next if statement; See statusLast
   StatusInfo statusNew;
   statusNew.nLumiStart = nLumiCurr;
   statusNew.nLumiEnd = nLumiCurr;
   statusNew.nStatus = nStatus;
 
   if (listStatus.empty()) {
     listStatus.push_back(statusNew);
   } else {
     auto &statusLastPre = listStatus.back();
     if (statusLastPre.nStatus == nStatus) {
       statusLastPre.nLumiEnd = nLumiCurr;
     } else {
       listStatus.push_back(statusNew);
     }
   }
 
   return 0;
 }

Member Data Documentation

◆ bIsStatusChambersInit_

Bool_t GEMDQMHarvester::bIsStatusChambersInit_

protected

Definition at line 154 of file GEMDQMHarvester.cc.

◆ fCutErr_

Float_t GEMDQMHarvester::fCutErr_

protected

Definition at line 116 of file GEMDQMHarvester.cc.

◆ fCutLowErr_

Float_t GEMDQMHarvester::fCutLowErr_

protected

Definition at line 116 of file GEMDQMHarvester.cc.

◆ fCutWarn_

Float_t GEMDQMHarvester::fCutWarn_

protected

Definition at line 116 of file GEMDQMHarvester.cc.

◆ listLayer_

std::vector<std::string> GEMDQMHarvester::listLayer_

protected

Definition at line 150 of file GEMDQMHarvester.cc.

◆ mapHistLumiFunc_

std::map<int, MonitorElement *> GEMDQMHarvester::mapHistLumiFunc_

protected

Definition at line 153 of file GEMDQMHarvester.cc.

◆ mapIdxLayer_

std::map<std::string, int> GEMDQMHarvester::mapIdxLayer_

protected

Definition at line 151 of file GEMDQMHarvester.cc.

◆ mapNumChPerChamber_

std::map<int, int> GEMDQMHarvester::mapNumChPerChamber_

protected

Definition at line 152 of file GEMDQMHarvester.cc.

◆ mapNumStatusChambersSummary_

std::map<IdChamber, NumStatus> GEMDQMHarvester::mapNumStatusChambersSummary_

protected

Definition at line 147 of file GEMDQMHarvester.cc.

◆ mapNumStatusVFATsSummary_

std::map<IdVFAT, NumStatus> GEMDQMHarvester::mapNumStatusVFATsSummary_

protected

Definition at line 148 of file GEMDQMHarvester.cc.

◆ mapStatusChambersSummary_

std::map<IdChamber, std::vector<StatusInfo> > GEMDQMHarvester::mapStatusChambersSummary_

protected

Definition at line 145 of file GEMDQMHarvester.cc.

◆ mapStatusVFATsSummary_

std::map<IdVFAT, std::vector<StatusInfo> > GEMDQMHarvester::mapStatusVFATsSummary_

protected

Definition at line 146 of file GEMDQMHarvester.cc.

◆ nBitErrAMC13_

const Int_t GEMDQMHarvester::nBitErrAMC13_ = 1

protected

Definition at line 133 of file GEMDQMHarvester.cc.

◆ nBitErrAMC_

const Int_t GEMDQMHarvester::nBitErrAMC_ = 2

protected

Definition at line 132 of file GEMDQMHarvester.cc.

◆ nBitErrOH_

const Int_t GEMDQMHarvester::nBitErrOH_ = 4

protected

Definition at line 130 of file GEMDQMHarvester.cc.

◆ nBitErrVFAT_

const Int_t GEMDQMHarvester::nBitErrVFAT_ = 6

protected

Definition at line 128 of file GEMDQMHarvester.cc.

◆ nBitOcc_

const Int_t GEMDQMHarvester::nBitOcc_ = 0

protected

Definition at line 134 of file GEMDQMHarvester.cc.

◆ nBitWarnAMC_

const Int_t GEMDQMHarvester::nBitWarnAMC_ = 3

protected

Definition at line 131 of file GEMDQMHarvester.cc.

◆ nBitWarnOH_

const Int_t GEMDQMHarvester::nBitWarnOH_ = 5

protected

Definition at line 129 of file GEMDQMHarvester.cc.

◆ nBitWarnVFAT_

const Int_t GEMDQMHarvester::nBitWarnVFAT_ = 7

protected

Definition at line 127 of file GEMDQMHarvester.cc.

◆ nCodeError_

const Int_t GEMDQMHarvester::nCodeError_ = 2

protected

Definition at line 123 of file GEMDQMHarvester.cc.

◆ nCodeFine_

const Int_t GEMDQMHarvester::nCodeFine_ = 1

protected

Definition at line 122 of file GEMDQMHarvester.cc.

◆ nCodeLowError_

const Int_t GEMDQMHarvester::nCodeLowError_ = 4

protected

Definition at line 125 of file GEMDQMHarvester.cc.

◆ nCodeWarning_

const Int_t GEMDQMHarvester::nCodeWarning_ = 3

protected

Definition at line 124 of file GEMDQMHarvester.cc.

◆ nMaxLumi_

const Int_t GEMDQMHarvester::nMaxLumi_ = 6000

protected

Definition at line 138 of file GEMDQMHarvester.cc.

◆ nNumVFATs_

const Int_t GEMDQMHarvester::nNumVFATs_ = 24

protected

Definition at line 136 of file GEMDQMHarvester.cc.

◆ nResolutionLumi_

Int_t GEMDQMHarvester::nResolutionLumi_

protected

Definition at line 140 of file GEMDQMHarvester.cc.

◆ strDirRecHit_

const std::string GEMDQMHarvester::strDirRecHit_ = "GEM/RecHits"

protected

Definition at line 119 of file GEMDQMHarvester.cc.

◆ strDirStatus_

const std::string GEMDQMHarvester::strDirStatus_ = "GEM/DAQStatus"

protected

Definition at line 120 of file GEMDQMHarvester.cc.

◆ strDirSummary_

const std::string GEMDQMHarvester::strDirSummary_ = "GEM/EventInfo"

protected

Definition at line 118 of file GEMDQMHarvester.cc.

Classes

Public Types

Public Member Functions

Static Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Detailed Description

Member Typedef Documentation

◆ IdChamber

◆ IdVFAT

◆ StatusInfo

◆ TableStatusNum

◆ TableStatusOcc

Constructor & Destructor Documentation

◆ GEMDQMHarvester()

◆ ~GEMDQMHarvester()

Member Function Documentation

◆ assessOneBin()

◆ copyLabels()

◆ createInactiveChannelFracHist()

◆ createLumiFuncHist()

◆ createSummaryHist()

◆ createSummaryVFAT()

◆ createTableWatchingSummary()

◆ dqmEndJob()

◆ dqmEndLuminosityBlock()

◆ drawSummaryHistogram()

◆ fillDescriptions()

◆ getGeometryInfo()

◆ refineSummaryHistogram()

◆ refineSummaryVFAT()

◆ UpdateStatusChamber() [1/3]

◆ UpdateStatusChamber() [2/3]

◆ UpdateStatusChamber() [3/3]

Member Data Documentation

◆ bIsStatusChambersInit_

◆ fCutErr_

◆ fCutLowErr_

◆ fCutWarn_

◆ listLayer_

◆ mapHistLumiFunc_

◆ mapIdxLayer_

◆ mapNumChPerChamber_

◆ mapNumStatusChambersSummary_

◆ mapNumStatusVFATsSummary_

◆ mapStatusChambersSummary_

◆ mapStatusVFATsSummary_

◆ nBitErrAMC13_

◆ nBitErrAMC_

◆ nBitErrOH_

◆ nBitErrVFAT_

◆ nBitOcc_

◆ nBitWarnAMC_

◆ nBitWarnOH_

◆ nBitWarnVFAT_

◆ nCodeError_

◆ nCodeFine_

◆ nCodeLowError_

◆ nCodeWarning_

◆ nMaxLumi_

◆ nNumVFATs_

◆ nResolutionLumi_

◆ strDirRecHit_

◆ strDirStatus_

◆ strDirSummary_