CMS 3D CMS Logo

List of all members | Classes | Public Types | Public Member Functions | Static Public Member Functions | Protected Types | Protected Member Functions | Protected Attributes
GEMDQMHarvester Class Reference
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, std::map< int, int > > mapIdxToChamberInOcc_
 
std::map< int, int > mapNumChPerChamber_
 
std::map< IdChamber, NumStatusmapNumStatusChambersSummary_
 
std::map< IdVFAT, NumStatusmapNumStatusVFATsSummary_
 
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 nCodeMasked_ = 5
 
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

◆ TableStatusNum

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

Definition at line 144 of file GEMDQMHarvester.cc.

◆ TableStatusOcc

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

Definition at line 143 of file GEMDQMHarvester.cc.

Constructor & Destructor Documentation

◆ GEMDQMHarvester()

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

Definition at line 159 of file GEMDQMHarvester.cc.

References looper::cfg.

159  {
160  fCutErr_ = cfg.getParameter<double>("cutErr");
161  fCutLowErr_ = cfg.getParameter<double>("cutLowErr");
162  fCutWarn_ = cfg.getParameter<double>("cutWarn");
163  nResolutionLumi_ = cfg.getParameter<int>("resolutionLumi");
164  bIsStatusChambersInit_ = false;
165 }
Bool_t bIsStatusChambersInit_

◆ ~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 418 of file GEMDQMHarvester.cc.

419  {
420  if (fNumErr < 0)
421  return nCodeMasked_;
422  if (fNumErr > fCutErr_ * fAll) // The error status criterion
423  return nCodeError_;
424  else if (fNumErr > fCutLowErr_ * fAll) // The low-error status criterion
425  return nCodeLowError_;
426  else if (fNumWarn > fCutWarn_ * fAll) // The warning status criterion
427  return nCodeWarning_;
428  else if (fNumOcc > 0)
429  return nCodeFine_;
430 
431  return 0;
432 }
const Int_t nCodeMasked_
const Int_t nCodeWarning_
const Int_t nCodeError_
const Int_t nCodeLowError_
const Int_t nCodeFine_

◆ copyLabels()

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

Definition at line 330 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().

330  {
331  Int_t nBinX = h2Src->getNbinsX(), nBinY = h2Src->getNbinsY();
332 
333  for (Int_t i = 1; i <= nBinX; i++) {
334  h2Dst->setBinLabel(i, h2Src->getTH2F()->GetXaxis()->GetBinLabel(i), 1);
335  }
336  for (Int_t i = 1; i <= nBinY; i++) {
337  h2Dst->setBinLabel(i, h2Src->getTH2F()->GetYaxis()->GetBinLabel(i), 2);
338  }
339  h2Dst->setTitle(h2Src->getTitle());
340  h2Dst->setXTitle(h2Src->getAxisTitle(1));
341  h2Dst->setYTitle(h2Src->getAxisTitle(2));
342 }
virtual void setTitle(const std::string &title)
set (ie. change) histogram/profile title
virtual std::string getAxisTitle(int axis=1) const
get x-, y- or z-axis title (axis=1, 2, 3 respectively)
virtual TH2F * getTH2F() const
virtual void setXTitle(std::string const &title)
virtual void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
virtual int getNbinsY() const
get # of bins in Y-axis
virtual std::string getTitle() const
get MonitorElement title
virtual void setYTitle(std::string const &title)
virtual int getNbinsX() const
get # of bins in X-axis

◆ createInactiveChannelFracHist()

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

Definition at line 760 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.

762  {
763  std::string strTitle = "The fraction of inactive channels in " + strSuffix.substr(1);
764  MonitorElement *h2InactiveChannel =
765  store->book1D("inactive_frac_chamber" + strSuffix, strTitle, nNumChamber, 0.5, nNumChamber + 0.5);
766  h2InactiveChannel->setXTitle("Chamber");
767  h2InactiveChannel->setYTitle("Fraction of inactive channels");
768  for (Int_t i = 1; i <= nNumChamber; i++) {
769  h2InactiveChannel->setBinLabel(i, Form("%i", i), 1);
770  }
771 
772  std::string (*funcNameCh)(std::string, Int_t) = getNameChamberOccNull;
773 
774  if (strSuffix.find("_GE11") != std::string::npos) {
775  funcNameCh = getNameChamberOccGE11;
776  } else if (strSuffix.find("_GE21") != std::string::npos) {
777  funcNameCh = getNameChamberOccGE21;
778  }
779 
780  for (Int_t nIdxCh = 1; nIdxCh <= nNumChamber; nIdxCh++) {
781  std::string strNameCh = funcNameCh(strSuffix, nIdxCh);
782  MonitorElement *h2SrcChamberOcc = store->get(strNameCh);
783  if (h2SrcChamberOcc == nullptr) {
784  // FIXME: It's about sending a message
785  continue;
786  }
787 
788  Int_t nNumBinX = h2SrcChamberOcc->getNbinsX();
789  Int_t nNumBinY = h2SrcChamberOcc->getNbinsY();
790  Int_t nNumAllChannel = nNumBinX * nNumBinY;
791  auto *histData = h2SrcChamberOcc->getTH2F();
792  auto *pdData = histData->GetArray();
793  Int_t nNumChannelInactive = 0;
794  for (Int_t j = 1; j <= nNumBinY; j++)
795  for (Int_t i = 1; i <= nNumBinX; i++) {
796  if (pdData[j * (nNumBinX + 2) + i] <= 0) {
797  nNumChannelInactive++;
798  }
799  }
800  h2InactiveChannel->setBinContent(nIdxCh, ((Double_t)nNumChannelInactive) / nNumAllChannel);
801  }
802 }
std::string getNameChamberOccGE21(std::string strSuffix, Int_t nIdxCh)
std::string getNameChamberOccNull(std::string strSuffix, Int_t nIdxChamber)
virtual TH2F * getTH2F() const
virtual void setXTitle(std::string const &title)
std::string getNameChamberOccGE11(std::string strSuffix, Int_t nIdxCh)
virtual void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
virtual int getNbinsY() const
get # of bins in Y-axis
virtual void setYTitle(std::string const &title)
virtual void setBinContent(int binx, double content)
set content of bin (1-D)
virtual MonitorElement * get(std::string const &fullpath) const
Definition: DQMStore.cc:712
virtual int getNbinsX() const
get # of bins in X-axis
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98

◆ createLumiFuncHist()

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

Definition at line 600 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().

603  {
604  auto &nNumCh = mapNumChPerChamber_[nIdxLayer];
605 
606  MonitorElement *h2Summary;
607 
608  //Int_t nLumiCurrLowRes = ( ( nLumiCurr - 1 ) / nResolutionLumi_ ) * nResolutionLumi_;
609  Int_t nNumBinLumi = ((nLumiCurr - 1) / nResolutionLumi_) + 1;
610  Int_t nMaxBin = 0;
611 
612  // Creating or Summoning the corresponding histogram
613  if (mapHistLumiFunc_.find(nIdxLayer) == mapHistLumiFunc_.end()) {
615  h2Summary = store->book2S("chamberStatus_inLumi" + strSuffix,
616  "Chamber status on lumi-block " + strSuffix.substr(1),
618  1.0,
619  (Float_t)(nMaxLumi_ + 1),
620  //nNumBinLumi, 1.0, (Float_t)( nLumiCurr + 1 ),
621  nNumCh,
622  0.5,
623  nNumCh + 0.5);
624  mapHistLumiFunc_[nIdxLayer] = h2Summary;
625 
626  h2Summary->setXTitle("Luminosity block");
627  h2Summary->setYTitle("Chamber");
628  for (Int_t i = 1; i <= nNumCh; i++) {
629  h2Summary->setBinLabel(i, Form("%i", i), 2);
630  }
631  } else {
632  h2Summary = mapHistLumiFunc_[nIdxLayer];
633  }
634 
635  for (Int_t nIdxCh = 1; nIdxCh <= nNumCh; nIdxCh++) {
636  auto &listStatus = mapStatusChambersSummary_[{nIdxLayer, nIdxCh}];
637 
638  Int_t nIdxStatus = 0;
639  for (Int_t nIdxLumi = 0; nIdxLumi < nNumBinLumi; nIdxLumi++) {
640  // Lumis covered by these values (nLumiStart <=, <= nLumiEnd) are counted for the current bin
641  Int_t nLumiStart = 1 + nIdxLumi * nResolutionLumi_;
642  Int_t nLumiEnd = (1 + nIdxLumi) * nResolutionLumi_;
643  if (nLumiEnd > nLumiCurr)
644  nLumiEnd = nLumiCurr;
645 
646  Int_t nStatusSum = 0;
647  while (true) { // No worries, nIdxStatus must increase and reach at listStatus.size()
648  // True: It was too past so that
649  // the lumi range of listStatus[ nIdxStatus ] is out of the coverage of the current bin
650  if (listStatus[nIdxStatus].nLumiEnd < nLumiStart) {
651  nIdxStatus++;
652  if (nIdxStatus >= (int)listStatus.size()) {
653  break; // For safety
654  }
655  continue;
656  }
657 
658  nStatusSum = listStatus[nIdxStatus].nStatus;
659 
660  // True: This is the last item of listStatus which is covered by the current bin
661  if (nIdxStatus + 1 >= (int)listStatus.size() || listStatus[nIdxStatus].nLumiEnd >= nLumiEnd) {
662  break;
663  }
664 
665  nIdxStatus++;
666  if (nIdxStatus >= (int)listStatus.size()) {
667  break; // For safety
668  }
669  }
670 
671  h2Summary->setBinContent(nIdxLumi + 1, nIdxCh, nStatusSum);
672  if (nMaxBin < nIdxLumi + 1)
673  nMaxBin = nIdxLumi + 1;
674  }
675  }
676 
677  for (Int_t nX = 1; nX <= nMaxBin; nX++) {
678  h2Summary->setBinContent(nX, 0, 1);
679  }
680 }
MonitorElement * book2S(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:254
const std::string strDirSummary_
void setCurrentFolder(std::string const &fullpath) override
Definition: DQMStore.h:647
std::map< IdChamber, std::vector< StatusInfo > > mapStatusChambersSummary_
virtual void setXTitle(std::string const &title)
std::map< int, MonitorElement * > mapHistLumiFunc_
virtual void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
const Int_t nMaxLumi_
virtual void setYTitle(std::string const &title)
virtual void setBinContent(int binx, double content)
set content of bin (1-D)
std::map< int, int > mapNumChPerChamber_

◆ createSummaryHist()

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

Definition at line 391 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().

391  {
392  //store->setCurrentFolder(strDirSummary_);
393 
394  Int_t nBinX = h2Src->getNbinsX(), nBinY = h2Src->getNbinsY();
395  h2Sum = store->book2D("reportSummaryMap", "", nBinX, 0.5, nBinX + 0.5, nBinY, 0.5, nBinY + 0.5);
396  h2Sum->setTitle("Summary plot");
397  h2Sum->setXTitle("Chamber");
398  h2Sum->setYTitle("Layer");
399 
400  for (Int_t i = 1; i <= nBinX; i++)
401  h2Sum->setBinLabel(i, h2Src->getTH2F()->GetXaxis()->GetBinLabel(i), 1);
402  for (Int_t i = 1; i <= nBinY; i++)
403  h2Sum->setBinLabel(i, listLayer_[i].substr(1), 2);
404 }
virtual void setTitle(const std::string &title)
set (ie. change) histogram/profile title
std::vector< std::string > listLayer_
virtual TH2F * getTH2F() const
virtual void setXTitle(std::string const &title)
virtual void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
virtual int getNbinsY() const
get # of bins in Y-axis
virtual void setYTitle(std::string const &title)
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:212
virtual int getNbinsX() const
get # of bins in X-axis

◆ createSummaryVFAT()

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

Definition at line 406 of file GEMDQMHarvester.cc.

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

409  {
410  //store->setCurrentFolder(strDirStatus_);
411  //store->setCurrentFolder(strDirSummary_);
412 
413  Int_t nBinX = h2Src->getNbinsX(), nBinY = h2Src->getNbinsY();
414  h2Sum = store->book2D("vfat_statusSummary" + strSuffix, "", nBinX, 0.5, nBinX + 0.5, nBinY, -0.5, nBinY - 0.5);
415  copyLabels(h2Src, h2Sum);
416 }
void copyLabels(MonitorElement *h2Src, MonitorElement *h2Dst)
virtual int getNbinsY() const
get # of bins in Y-axis
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:212
virtual int getNbinsX() const
get # of bins in X-axis

◆ createTableWatchingSummary()

void GEMDQMHarvester::createTableWatchingSummary ( )
protected

Definition at line 312 of file GEMDQMHarvester.cc.

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

312  {
314  return;
315 
316  for (const auto &[nIdxLayer, nNumCh] : mapNumChPerChamber_) {
317  for (Int_t i = 1; i <= nNumCh; i++) {
318  mapStatusChambersSummary_[{nIdxLayer, i}] = std::vector<StatusInfo>();
319  mapNumStatusChambersSummary_[{nIdxLayer, i}] = NumStatus();
320  for (Int_t j = 1; j <= nNumVFATs_; j++) {
321  mapStatusVFATsSummary_[{nIdxLayer, i, j}] = std::vector<StatusInfo>();
322  mapNumStatusVFATsSummary_[{nIdxLayer, i, j}] = NumStatus();
323  }
324  }
325  }
326 
327  bIsStatusChambersInit_ = true;
328 }
std::map< IdChamber, std::vector< StatusInfo > > mapStatusChambersSummary_
std::map< IdVFAT, NumStatus > mapNumStatusVFATsSummary_
std::map< IdVFAT, std::vector< StatusInfo > > mapStatusVFATsSummary_
Bool_t bIsStatusChambersInit_
const Int_t nNumVFATs_
std::map< IdChamber, NumStatus > mapNumStatusChambersSummary_
std::map< int, int > mapNumChPerChamber_

◆ 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 176 of file GEMDQMHarvester.cc.

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

179  {
181  Int_t nLumiCurr = iLumi.id().luminosityBlock();
182  drawSummaryHistogram(store, nLumiCurr);
183 }
void drawSummaryHistogram(edm::Service< DQMStore > &store, Int_t nLumiCurr)

◆ drawSummaryHistogram()

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

Definition at line 185 of file GEMDQMHarvester.cc.

References dqm::implementation::IBooker::bookFloat(), dqm::impl::MonitorElement::Fill(), dqm::implementation::IGetter::get(), dqm::impl::MonitorElement::getAxisTitle(), dqm::impl::MonitorElement::getNbinsX(), dqm::legacy::MonitorElement::getTH2F(), dqm::impl::MonitorElement::getTitle(), mps_fire::i, dqm::implementation::DQMStore::setCurrentFolder(), dqm::impl::MonitorElement::setTitle(), dqm::impl::MonitorElement::setXTitle(), dqm::impl::MonitorElement::setYTitle(), contentValuesCheck::ss, and AlCaHLTBitMon_QueryRunRegistry::string.

185  {
186  Float_t fReportSummary = -1.0;
187 
188  std::string strSrcStatusA = "GEM/DAQStatus/chamberAllStatus";
189  std::string strSrcStatusE = "GEM/DAQStatus/chamberErrors";
190  std::string strSrcStatusW = "GEM/DAQStatus/chamberWarnings";
191  std::string strSrcStatusEVFAT = "GEM/DAQStatus/chamberVFATErrors";
192  std::string strSrcStatusWVFAT = "GEM/DAQStatus/chamberVFATWarnings";
193  std::string strSrcStatusEOH = "GEM/DAQStatus/chamberOHErrors";
194  std::string strSrcStatusWOH = "GEM/DAQStatus/chamberOHWarnings";
195  std::string strSrcStatusEAMC = "GEM/DAQStatus/chamberAMCErrors";
196  std::string strSrcStatusWAMC = "GEM/DAQStatus/chamberAMCWarnings";
197  std::string strSrcStatusEAMC13 = "GEM/DAQStatus/chamberAMC13Errors";
198 
199  std::string strSrcVFATOcc = "GEM/Digis/occ";
200  std::string strSrcVFATStatusW = "GEM/DAQStatus/vfat_statusWarnSum";
201  std::string strSrcVFATStatusE = "GEM/DAQStatus/vfat_statusErrSum";
202 
204 
205  MonitorElement *h2SrcStatusA = store->get(strSrcStatusA);
206  MonitorElement *h2SrcStatusE = store->get(strSrcStatusE);
207  MonitorElement *h2SrcStatusW = store->get(strSrcStatusW);
208  MonitorElement *h2SrcStatusEVFAT = store->get(strSrcStatusEVFAT);
209  MonitorElement *h2SrcStatusWVFAT = store->get(strSrcStatusWVFAT);
210  MonitorElement *h2SrcStatusEOH = store->get(strSrcStatusEOH);
211  MonitorElement *h2SrcStatusWOH = store->get(strSrcStatusWOH);
212  MonitorElement *h2SrcStatusEAMC = store->get(strSrcStatusEAMC);
213  MonitorElement *h2SrcStatusWAMC = store->get(strSrcStatusWAMC);
214  MonitorElement *h2SrcStatusEAMC13 = store->get(strSrcStatusEAMC13);
215 
216  std::string strTitleSummary = "summary";
217 
218  getGeometryInfo(store, h2SrcStatusE);
219 
220  if (h2SrcStatusA != nullptr && h2SrcStatusE != nullptr && h2SrcStatusW != nullptr) {
221  MonitorElement *h2Sum = nullptr;
222  createSummaryHist(store, h2SrcStatusE, h2Sum);
224 
225  std::vector<MonitorElement *> listOccPlots(listLayer_.size() + 1); // The index starts at 1
226  for (const auto &strSuffix : listLayer_) {
227  if (mapIdxLayer_.find(strSuffix) == mapIdxLayer_.end())
228  continue;
229  auto nIdxLayer = mapIdxLayer_[strSuffix];
230  MonitorElement *h2SrcVFATOcc = store->get(strSrcVFATOcc + strSuffix);
231  if (h2SrcVFATOcc == nullptr)
232  continue;
233  listOccPlots[nIdxLayer] = h2SrcVFATOcc;
234  mapIdxToChamberInOcc_[nIdxLayer] = {};
235  // Obtaining the bin indices of chambers from their labels
236  for (Int_t i = 1; i <= h2SrcVFATOcc->getNbinsX(); i++) {
237  std::string strLabel = h2SrcVFATOcc->getTH2F()->GetXaxis()->GetBinLabel(i);
238  std::string strSrc = strLabel;
239  std::vector<int> listNumExtract;
240 
241  while (!strSrc.empty()) {
242  auto nPosDigit = strSrc.find_first_of("0123456789");
243  if (nPosDigit == std::string::npos)
244  break;
245  std::stringstream ss;
246  ss << strSrc.substr(nPosDigit);
247  Int_t nExtract;
248  ss >> nExtract;
249  if (!ss.eof()) {
250  ss >> strSrc;
251  } else {
252  strSrc = "";
253  }
254  listNumExtract.push_back(nExtract);
255  }
256 
257  if (listNumExtract.empty()) { // Errneous case; but the job should not be dead
258  edm::LogError("GEMDQMHarvester") << "Error: Wrong label of GEM VFAT occupancy plot: " << strLabel;
259  } else {
260  mapIdxToChamberInOcc_[nIdxLayer][listNumExtract[0]] = i;
261  }
262  }
263  }
264 
265  fReportSummary = refineSummaryHistogram(strTitleSummary,
266  h2Sum,
267  listOccPlots,
268  h2SrcStatusA,
269  h2SrcStatusE,
270  h2SrcStatusW,
271  h2SrcStatusEVFAT,
272  h2SrcStatusWVFAT,
273  h2SrcStatusEOH,
274  h2SrcStatusWOH,
275  h2SrcStatusEAMC,
276  h2SrcStatusWAMC,
277  h2SrcStatusEAMC13,
278  nLumiCurr);
279 
280  for (const auto &strSuffix : listLayer_) {
281  if (mapIdxLayer_.find(strSuffix) == mapIdxLayer_.end())
282  continue;
283  auto nIdxLayer = mapIdxLayer_[strSuffix];
284  MonitorElement *h2SrcVFATOcc = store->get(strSrcVFATOcc + strSuffix);
285  MonitorElement *h2SrcVFATStatusW = store->get(strSrcVFATStatusW + strSuffix);
286  MonitorElement *h2SrcVFATStatusE = store->get(strSrcVFATStatusE + strSuffix);
287  if (h2SrcVFATOcc == nullptr || h2SrcVFATStatusW == nullptr || h2SrcVFATStatusE == nullptr)
288  continue;
289 
290  MonitorElement *h2SumVFAT = nullptr;
291  createSummaryVFAT(store, h2SrcVFATStatusE, strSuffix, h2SumVFAT);
292  refineSummaryVFAT(strSuffix, h2SumVFAT, h2SrcVFATOcc, h2SrcVFATStatusE, h2SrcVFATStatusW, nLumiCurr, nIdxLayer);
293  TString strNewTitle = h2SrcVFATStatusE->getTitle();
294  h2SumVFAT->setTitle((const char *)strNewTitle.ReplaceAll("errors", "errors/warnings"));
295  h2SumVFAT->setXTitle(h2SrcVFATStatusE->getAxisTitle(1));
296  h2SumVFAT->setYTitle(h2SrcVFATStatusE->getAxisTitle(2));
297 
298  createLumiFuncHist(store, strSuffix, nIdxLayer, nLumiCurr);
299  }
300  }
301 
302  for (const auto &strSuffix : listLayer_) {
303  if (mapIdxLayer_.find(strSuffix) == mapIdxLayer_.end())
304  continue;
305  Int_t nIdxLayer = mapIdxLayer_[strSuffix];
306  createInactiveChannelFracHist(store, strSuffix, mapNumChPerChamber_[nIdxLayer]);
307  }
308 
309  store->bookFloat("reportSummary")->Fill(fReportSummary);
310 }
virtual void setTitle(const std::string &title)
set (ie. change) histogram/profile title
std::map< std::string, int > mapIdxLayer_
MonitorElement * bookFloat(TString const &name, FUNC onbooking=NOOP())
Definition: DQMStore.h:80
void createSummaryVFAT(edm::Service< DQMStore > &store, MonitorElement *h2Src, std::string strSuffix, MonitorElement *&h2Sum)
void createInactiveChannelFracHist(edm::Service< DQMStore > &store, std::string strSuffix, Int_t nNumChamber)
const std::string strDirSummary_
void setCurrentFolder(std::string const &fullpath) override
Definition: DQMStore.h:647
virtual std::string getAxisTitle(int axis=1) const
get x-, y- or z-axis title (axis=1, 2, 3 respectively)
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)
Log< level::Error, false > LogError
void createTableWatchingSummary()
std::vector< std::string > listLayer_
void Fill(long long x)
virtual TH2F * getTH2F() const
virtual void setXTitle(std::string const &title)
void getGeometryInfo(edm::Service< DQMStore > &store, MonitorElement *h2Src)
virtual std::string getTitle() const
get MonitorElement title
virtual void setYTitle(std::string const &title)
std::map< int, std::map< int, int > > mapIdxToChamberInOcc_
virtual MonitorElement * get(std::string const &fullpath) const
Definition: DQMStore.cc:712
virtual int getNbinsX() const
get # of bins in X-axis
void createSummaryHist(edm::Service< DQMStore > &store, MonitorElement *h2Src, MonitorElement *&h2Sum)
Int_t refineSummaryVFAT(std::string strName, MonitorElement *h2Sum, MonitorElement *h2SrcOcc, MonitorElement *h2SrcStatusE, MonitorElement *h2SrcStatusW, Int_t nLumiCurr, Int_t nIdxLayer)
std::map< int, int > mapNumChPerChamber_
void createLumiFuncHist(edm::Service< DQMStore > &store, std::string strSuffix, Int_t nIdxLayer, Int_t nLumiCurr)

◆ fillDescriptions()

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

Definition at line 167 of file GEMDQMHarvester.cc.

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

167  {
169  desc.add<double>("cutErr", 0.05);
170  desc.add<double>("cutLowErr", 0.00);
171  desc.add<double>("cutWarn", 0.05);
172  desc.add<int>("resolutionLumi", 1);
173  descriptions.add("GEMDQMHarvester", desc);
174 }
void add(std::string const &label, ParameterSetDescription const &psetDescription)

◆ getGeometryInfo()

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

Definition at line 344 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.

344  {
345  listLayer_.clear();
346  mapIdxLayer_.clear();
347  mapNumChPerChamber_.clear();
348 
349  if (h2Src != nullptr) { // For online and offline
350  Int_t nBinY = h2Src->getNbinsY();
351  listLayer_.push_back("");
352 
353  for (Int_t i = 1; i <= nBinY; i++) {
354  std::string strLabelFull = h2Src->getTH2F()->GetYaxis()->GetBinLabel(i);
355  auto nPos = strLabelFull.find(';');
356  auto strLayer = strLabelFull.substr(nPos + 1);
357  Int_t nBinXActual = (Int_t)(h2Src->getBinContent(0, i) + 0.5);
358  if (nBinXActual > 108) { // When the number seems wrong
359  if (strLayer.find("GE11") != std::string::npos) {
360  nBinXActual = 36;
361  } else if (strLayer.find("GE21") != std::string::npos) {
362  nBinXActual = 18;
363  } else if (strLayer.find("GE01") != std::string::npos) {
364  nBinXActual = 36;
365  }
366  }
367  listLayer_.push_back(strLayer);
368  mapIdxLayer_[strLayer] = i;
369  mapNumChPerChamber_[i] = nBinXActual;
370  }
371  } else { // For others (validation and...?)
372  listLayer_.push_back("");
373  if (store->get("GEM/Digis/occupancy_GE11-M-L1/occ_GE11-M-01L1-S") != nullptr) {
374  listLayer_.push_back("_GE11-P-L2");
375  listLayer_.push_back("_GE11-P-L1");
376  listLayer_.push_back("_GE11-M-L1");
377  listLayer_.push_back("_GE11-M-L2");
378  mapIdxLayer_["_GE11-P-L2"] = 1;
379  mapIdxLayer_["_GE11-P-L1"] = 2;
380  mapIdxLayer_["_GE11-M-L1"] = 3;
381  mapIdxLayer_["_GE11-M-L2"] = 4;
382  mapNumChPerChamber_[1] = 36;
383  mapNumChPerChamber_[2] = 36;
384  mapNumChPerChamber_[3] = 36;
385  mapNumChPerChamber_[4] = 36;
386  }
387  // FIXME: How about GE21 and ME0?
388  }
389 }
std::map< std::string, int > mapIdxLayer_
std::vector< std::string > listLayer_
virtual TH2F * getTH2F() const
virtual int getNbinsY() const
get # of bins in Y-axis
virtual MonitorElement * get(std::string const &fullpath) const
Definition: DQMStore.cc:712
std::map< int, int > mapNumChPerChamber_
virtual double getBinContent(int binx) const
get content of bin (1-D)

◆ 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 435 of file GEMDQMHarvester.cc.

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

448  {
449  Int_t nBinY = h2Sum->getNbinsY();
450  Int_t nAllBin = 0, nFineBin = 0;
451  for (Int_t j = 1; j <= nBinY; j++) {
452  Int_t nBinX = mapNumChPerChamber_[j];
453  auto h2SrcOcc = listOccPlots[j];
454  auto &mapIdxOccChamber = mapIdxToChamberInOcc_[j];
455  Int_t nBinYOcc = 0;
456  if (h2SrcOcc != nullptr) {
457  nBinYOcc = h2SrcOcc->getNbinsY();
458  }
459 
460  h2Sum->setBinContent(0, j, nBinX);
461  for (Int_t i = 1; i <= nBinX; i++) {
462  Float_t fOcc = 0;
463  Int_t nIdxChOcc = mapIdxOccChamber[i];
464  for (Int_t r = 1; r <= nBinYOcc; r++) {
465  fOcc += h2SrcOcc->getBinContent(nIdxChOcc, r);
466  }
467 
468  Float_t fStatusAll = h2SrcStatusA->getBinContent(i, j);
469  Float_t fStatusErr = h2SrcStatusE->getBinContent(i, j);
470  Float_t fStatusWarn = h2SrcStatusW->getBinContent(i, j);
471  Float_t fStatusErrVFAT = h2SrcStatusEVFAT != nullptr ? h2SrcStatusEVFAT->getBinContent(i, j) : 0;
472  Float_t fStatusWarnVFAT = h2SrcStatusWVFAT != nullptr ? h2SrcStatusWVFAT->getBinContent(i, j) : 0;
473  Float_t fStatusErrOH = h2SrcStatusEOH != nullptr ? h2SrcStatusEOH->getBinContent(i, j) : 0;
474  Float_t fStatusWarnOH = h2SrcStatusWOH != nullptr ? h2SrcStatusWOH->getBinContent(i, j) : 0;
475  Float_t fStatusErrAMC = h2SrcStatusEAMC != nullptr ? h2SrcStatusEAMC->getBinContent(i, j) : 0;
476  Float_t fStatusWarnAMC = h2SrcStatusWAMC != nullptr ? h2SrcStatusWAMC->getBinContent(i, j) : 0;
477  Float_t fStatusErrAMC13 = h2SrcStatusEAMC13 != nullptr ? h2SrcStatusEAMC13->getBinContent(i, j) : 0;
478  NumStatus numStatus(fStatusAll,
479  fOcc,
480  fStatusErrVFAT,
481  fStatusWarnVFAT,
482  fStatusErrOH,
483  fStatusWarnOH,
484  fStatusErrAMC,
485  fStatusWarnAMC,
486  fStatusErrAMC13);
487  UpdateStatusChamber(j, i, nLumiCurr, numStatus);
488 
489  Int_t nRes = assessOneBin(strName, i, j, fStatusAll, fOcc, fStatusErr, fStatusWarn);
490  if (nRes == 1)
491  nFineBin++;
492 
493  h2Sum->setBinContent(i, j, (Float_t)nRes);
494  nAllBin++;
495  }
496  }
497 
498  return ((Float_t)nFineBin) / nAllBin;
499 }
Int_t assessOneBin(std::string strName, Int_t nIdxX, Int_t nIdxY, Float_t fAll, Float_t fNumOcc, Float_t fNumErr, Float_t fNumWarn)
virtual int getNbinsY() const
get # of bins in Y-axis
virtual void setBinContent(int binx, double content)
set content of bin (1-D)
std::map< int, std::map< int, int > > mapIdxToChamberInOcc_
Int_t UpdateStatusChamber(Int_t nIdxLayer, Int_t nIdxCh, Int_t nLumiCurr, NumStatus numStatus)
std::map< int, int > mapNumChPerChamber_
virtual double getBinContent(int binx) const
get content of bin (1-D)

◆ 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 501 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().

507  {
508  Int_t nBinY = h2Sum->getNbinsY();
509  for (Int_t j = 1; j <= nBinY; j++) {
510  Int_t nBinX = h2Sum->getNbinsX();
511  for (Int_t i = 1; i <= nBinX; i++) {
512  Float_t fOcc = h2SrcOcc->getBinContent(i, j);
513  Float_t fStatusErr = h2SrcStatusE->getBinContent(i, j);
514  Float_t fStatusWarn = h2SrcStatusW->getBinContent(i, j);
515  Float_t fStatusAll = fOcc + fStatusErr + fStatusWarn;
516  NumStatus numStatus(fStatusAll, fOcc, fStatusErr, fStatusWarn, 0, 0, 0, 0, 0);
517  UpdateStatusChamber(nIdxLayer, i, j, nLumiCurr, numStatus);
518 
519  Int_t nRes = assessOneBin(strName, i, j, fStatusAll, fOcc, fStatusErr, fStatusWarn);
520  h2Sum->setBinContent(i, j, (Float_t)nRes);
521  }
522  }
523 
524  return 0;
525 }
Int_t assessOneBin(std::string strName, Int_t nIdxX, Int_t nIdxY, Float_t fAll, Float_t fNumOcc, Float_t fNumErr, Float_t fNumWarn)
virtual int getNbinsY() const
get # of bins in Y-axis
virtual void setBinContent(int binx, double content)
set content of bin (1-D)
Int_t UpdateStatusChamber(Int_t nIdxLayer, Int_t nIdxCh, Int_t nLumiCurr, NumStatus numStatus)
virtual int getNbinsX() const
get # of bins in X-axis
virtual double getBinContent(int binx) const
get content of bin (1-D)

◆ UpdateStatusChamber() [1/3]

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

Definition at line 527 of file GEMDQMHarvester.cc.

527  {
529  return 0;
530  if (0 >= nIdxCh || nIdxCh > mapNumChPerChamber_[nIdxLayer])
531  return 0;
532  auto &listStatus = mapStatusChambersSummary_[{nIdxLayer, nIdxCh}];
533  auto &numStatusPrev = mapNumStatusChambersSummary_[{nIdxLayer, nIdxCh}];
534  return UpdateStatusChamber(listStatus, numStatusPrev, nLumiCurr, numStatus);
535 }
std::map< IdChamber, std::vector< StatusInfo > > mapStatusChambersSummary_
Bool_t bIsStatusChambersInit_
Int_t UpdateStatusChamber(Int_t nIdxLayer, Int_t nIdxCh, Int_t nLumiCurr, NumStatus numStatus)
std::map< IdChamber, NumStatus > mapNumStatusChambersSummary_
std::map< int, int > mapNumChPerChamber_

◆ 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 537 of file GEMDQMHarvester.cc.

538  {
540  return 0;
541  if (0 >= nIdxCh || nIdxCh > mapNumChPerChamber_[nIdxLayer])
542  return 0;
543  if (0 >= nIdxVFAT || nIdxVFAT > nNumVFATs_)
544  return 0;
545  auto &listStatus = mapStatusVFATsSummary_[{nIdxLayer, nIdxCh, nIdxVFAT}];
546  auto &numStatusPrev = mapNumStatusVFATsSummary_[{nIdxLayer, nIdxCh, nIdxVFAT}];
547  return UpdateStatusChamber(listStatus, numStatusPrev, nLumiCurr, numStatus);
548 }
std::map< IdVFAT, NumStatus > mapNumStatusVFATsSummary_
std::map< IdVFAT, std::vector< StatusInfo > > mapStatusVFATsSummary_
Bool_t bIsStatusChambersInit_
const Int_t nNumVFATs_
Int_t UpdateStatusChamber(Int_t nIdxLayer, Int_t nIdxCh, Int_t nLumiCurr, NumStatus numStatus)
std::map< int, int > mapNumChPerChamber_

◆ UpdateStatusChamber() [3/3]

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

Definition at line 550 of file GEMDQMHarvester.cc.

References GEMDQMHarvester::NumStatus::fNumErrAMC13_, GEMDQMHarvester::NumStatus::fNumErrAMC_, GEMDQMHarvester::NumStatus::fNumErrOH_, GEMDQMHarvester::NumStatus::fNumErrVFAT_, GEMDQMHarvester::NumStatus::fNumOcc_, GEMDQMHarvester::NumStatus::fNumTotal_, GEMDQMHarvester::NumStatus::fNumWarnAMC_, GEMDQMHarvester::NumStatus::fNumWarnOH_, GEMDQMHarvester::NumStatus::fNumWarnVFAT_, GEMDQMHarvester::PreStatusInfo::nLumiEnd, GEMDQMHarvester::PreStatusInfo::nLumiStart, and GEMDQMHarvester::PreStatusInfo::nStatus.

553  {
554  // First of all, the current lumi section will be assessed, of which the result will be stored in nStatus
555  Int_t nStatus = 0;
556 
557  Float_t fNumAddErrVFAT = numStatusNew.fNumErrVFAT_ - numStatus.fNumErrVFAT_;
558  Float_t fNumAddWarnVFAT = numStatusNew.fNumWarnVFAT_ - numStatus.fNumWarnVFAT_;
559  Float_t fNumAddErrOH = numStatusNew.fNumErrOH_ - numStatus.fNumErrOH_;
560  Float_t fNumAddWarnOH = numStatusNew.fNumWarnOH_ - numStatus.fNumWarnOH_;
561  Float_t fNumAddErrAMC = numStatusNew.fNumErrAMC_ - numStatus.fNumErrAMC_;
562  Float_t fNumAddWarnAMC = numStatusNew.fNumWarnAMC_ - numStatus.fNumWarnAMC_;
563  Float_t fNumAddErrAMC13 = numStatusNew.fNumErrAMC13_ - numStatus.fNumErrAMC13_;
564 
565  numStatus.fNumTotal_ = numStatusNew.fNumTotal_;
566  numStatus.fNumOcc_ = numStatusNew.fNumOcc_;
567  numStatus.fNumErrVFAT_ = numStatusNew.fNumErrVFAT_;
568  numStatus.fNumWarnVFAT_ = numStatusNew.fNumWarnVFAT_;
569  numStatus.fNumErrOH_ = numStatusNew.fNumErrOH_;
570  numStatus.fNumWarnOH_ = numStatusNew.fNumWarnOH_;
571  numStatus.fNumErrAMC_ = numStatusNew.fNumErrAMC_;
572  numStatus.fNumWarnAMC_ = numStatusNew.fNumWarnAMC_;
573  numStatus.fNumErrAMC13_ = numStatusNew.fNumErrAMC13_;
574 
575  nStatus = (numStatusNew.fNumOcc_ > 0 ? 1 << nBitOcc_ : 0) | (fNumAddErrAMC13 > 0 ? 1 << nBitErrAMC13_ : 0) |
576  (fNumAddErrAMC > 0 ? 1 << nBitErrAMC_ : 0) | (fNumAddWarnAMC > 0 ? 1 << nBitWarnAMC_ : 0) |
577  (fNumAddErrOH > 0 ? 1 << nBitErrOH_ : 0) | (fNumAddWarnOH > 0 ? 1 << nBitWarnOH_ : 0) |
578  (fNumAddErrVFAT > 0 ? 1 << nBitErrVFAT_ : 0) | (fNumAddWarnVFAT > 0 ? 1 << nBitWarnVFAT_ : 0);
579 
580  // Only used in the next if statement; See statusLast
581  StatusInfo statusNew;
582  statusNew.nLumiStart = nLumiCurr;
583  statusNew.nLumiEnd = nLumiCurr;
584  statusNew.nStatus = nStatus;
585 
586  if (listStatus.empty()) {
587  listStatus.push_back(statusNew);
588  } else {
589  auto &statusLastPre = listStatus.back();
590  if (statusLastPre.nStatus == nStatus) {
591  statusLastPre.nLumiEnd = nLumiCurr;
592  } else {
593  listStatus.push_back(statusNew);
594  }
595  }
596 
597  return 0;
598 }
const Int_t nBitErrVFAT_
const Int_t nBitErrOH_
struct GEMDQMHarvester::PreStatusInfo StatusInfo
const Int_t nBitOcc_
const Int_t nBitErrAMC13_
const Int_t nBitErrAMC_
const Int_t nBitWarnVFAT_
const Int_t nBitWarnAMC_
const Int_t nBitWarnOH_

Member Data Documentation

◆ bIsStatusChambersInit_

Bool_t GEMDQMHarvester::bIsStatusChambersInit_
protected

Definition at line 156 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 151 of file GEMDQMHarvester.cc.

◆ mapHistLumiFunc_

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

Definition at line 155 of file GEMDQMHarvester.cc.

◆ mapIdxLayer_

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

Definition at line 152 of file GEMDQMHarvester.cc.

◆ mapIdxToChamberInOcc_

std::map<int, std::map<int, int> > GEMDQMHarvester::mapIdxToChamberInOcc_
protected

Definition at line 154 of file GEMDQMHarvester.cc.

◆ mapNumChPerChamber_

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

Definition at line 153 of file GEMDQMHarvester.cc.

◆ mapNumStatusChambersSummary_

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

Definition at line 148 of file GEMDQMHarvester.cc.

◆ mapNumStatusVFATsSummary_

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

Definition at line 149 of file GEMDQMHarvester.cc.

◆ mapStatusChambersSummary_

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

Definition at line 146 of file GEMDQMHarvester.cc.

◆ mapStatusVFATsSummary_

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

Definition at line 147 of file GEMDQMHarvester.cc.

◆ nBitErrAMC13_

const Int_t GEMDQMHarvester::nBitErrAMC13_ = 1
protected

Definition at line 134 of file GEMDQMHarvester.cc.

◆ nBitErrAMC_

const Int_t GEMDQMHarvester::nBitErrAMC_ = 2
protected

Definition at line 133 of file GEMDQMHarvester.cc.

◆ nBitErrOH_

const Int_t GEMDQMHarvester::nBitErrOH_ = 4
protected

Definition at line 131 of file GEMDQMHarvester.cc.

◆ nBitErrVFAT_

const Int_t GEMDQMHarvester::nBitErrVFAT_ = 6
protected

Definition at line 129 of file GEMDQMHarvester.cc.

◆ nBitOcc_

const Int_t GEMDQMHarvester::nBitOcc_ = 0
protected

Definition at line 135 of file GEMDQMHarvester.cc.

◆ nBitWarnAMC_

const Int_t GEMDQMHarvester::nBitWarnAMC_ = 3
protected

Definition at line 132 of file GEMDQMHarvester.cc.

◆ nBitWarnOH_

const Int_t GEMDQMHarvester::nBitWarnOH_ = 5
protected

Definition at line 130 of file GEMDQMHarvester.cc.

◆ nBitWarnVFAT_

const Int_t GEMDQMHarvester::nBitWarnVFAT_ = 7
protected

Definition at line 128 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.

◆ nCodeMasked_

const Int_t GEMDQMHarvester::nCodeMasked_ = 5
protected

Definition at line 126 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 139 of file GEMDQMHarvester.cc.

◆ nNumVFATs_

const Int_t GEMDQMHarvester::nNumVFATs_ = 24
protected

Definition at line 137 of file GEMDQMHarvester.cc.

◆ nResolutionLumi_

Int_t GEMDQMHarvester::nResolutionLumi_
protected

Definition at line 141 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.