Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Lumi-based status plots on GEM onlineDQM, backport to CMSSW_12_2_X #37204

Merged
merged 1 commit into from Mar 14, 2022
Merged

Lumi-based status plots on GEM onlineDQM, backport to CMSSW_12_2_X #37204

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
30 changes: 30 additions & 0 deletions DQM/GEM/interface/GEMDAQStatusSource.h
View file
Open in desktop
Expand Up @@ -16,6 +16,9 @@

#include "Validation/MuonGEMHits/interface/GEMValidationUtils.h"

#include "CondFormats/DataRecord/interface/GEMeMapRcd.h"
#include "CondFormats/GEMObjects/interface/GEMeMap.h"
#include "CondFormats/GEMObjects/interface/GEMROMapping.h"
#include "DataFormats/GEMDigi/interface/GEMDigiCollection.h"
#include "DataFormats/GEMDigi/interface/GEMVFATStatusCollection.h"
#include "DataFormats/GEMDigi/interface/GEMOHStatusCollection.h"
Expand All @@ -35,6 +38,8 @@ class GEMDAQStatusSource : public GEMDQMBase {
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions);

protected:
void LoadROMap(edm::EventSetup const &iSetup);

void dqmBeginRun(edm::Run const &, edm::EventSetup const &) override{};
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
void analyze(edm::Event const &e, edm::EventSetup const &eSetup) override;
Expand All @@ -44,6 +49,18 @@ class GEMDAQStatusSource : public GEMDQMBase {
bErr = true;
};

void FillStatusSummaryPlot(std::map<ME4IdsKey, bool> &mapChamber,
MonitorElement *h2Plot,
std::map<ME4IdsKey, bool> *pmapSummary = nullptr) {
for (auto const &[key4, bFlag] : mapChamber) { // bFlag is not used
ME3IdsKey key3 = key4Tokey3(key4);
Int_t nChamber = keyToChamber(key4);
h2Plot->Fill(nChamber, mapStationToIdx_[key3]);
if (pmapSummary != nullptr)
(*pmapSummary)[key4] = true;
}
};

private:
int ProcessWithMEMap3(BookingHelper &bh, ME3IdsKey key) override;
int ProcessWithMEMap3WithChamber(BookingHelper &bh, ME4IdsKey key) override;
Expand All @@ -53,6 +70,10 @@ class GEMDAQStatusSource : public GEMDQMBase {
void SetLabelOHStatus(MonitorElement *h2Status);
void SetLabelVFATStatus(MonitorElement *h2Status);

edm::ESGetToken<GEMeMap, GEMeMapRcd> gemEMapToken_;
//std::shared_ptr<GEMROMapping> gemROMap_;
const GEMeMap *gemEMap_;

edm::EDGetToken tagVFAT_;
edm::EDGetToken tagOH_;
edm::EDGetToken tagAMC_;
Expand All @@ -72,10 +93,19 @@ class GEMDAQStatusSource : public GEMDQMBase {
MonitorElement *h2SummaryStatusAll;
MonitorElement *h2SummaryStatusWarning;
MonitorElement *h2SummaryStatusError;
MonitorElement *h2SummaryStatusVFATWarning;
MonitorElement *h2SummaryStatusVFATError;
MonitorElement *h2SummaryStatusOHWarning;
MonitorElement *h2SummaryStatusOHError;
MonitorElement *h2SummaryStatusAMCWarning;
MonitorElement *h2SummaryStatusAMCError;
MonitorElement *h2SummaryStatusAMC13Error;

Int_t nBXMin_, nBXMax_;

std::map<UInt_t, int> mapFEDIdToRe_;
std::map<int, std::vector<GEMDetId>> mapAMC13ToListChamber_;
std::map<std::tuple<int, int>, std::vector<GEMDetId>> mapAMCToListChamber_;
Int_t nAMCSlots_;

int nBitAMC13_ = 10;
Expand Down
2 changes: 2 additions & 0 deletions DQM/GEM/plugins/BuildFile.xml
View file
Open in desktop
Expand Up @@ -3,5 +3,7 @@
<use name="FWCore/Framework"/>
<use name="FWCore/ParameterSet"/>
<use name="DQM/GEM"/>
<use name="CondFormats/DataRecord"/>
<use name="CondFormats/GEMObjects"/>
<use name="DQMServices/Core"/>
</library>
130 changes: 102 additions & 28 deletions DQM/GEM/plugins/GEMDAQStatusSource.cc
View file
Open in desktop
Expand Up @@ -10,6 +10,7 @@ GEMDAQStatusSource::GEMDAQStatusSource(const edm::ParameterSet &cfg) : GEMDQMBas
tagAMC13_ = consumes<GEMAMC13StatusCollection>(cfg.getParameter<edm::InputTag>("AMC13InputLabel"));

nAMCSlots_ = cfg.getParameter<Int_t>("AMCSlots");
gemEMapToken_ = esConsumes<GEMeMap, GEMeMapRcd, edm::Transition::BeginRun>();
}

void GEMDAQStatusSource::fillDescriptions(edm::ConfigurationDescriptions &descriptions) {
Expand All @@ -25,6 +26,40 @@ void GEMDAQStatusSource::fillDescriptions(edm::ConfigurationDescriptions &descri
descriptions.add("GEMDAQStatusSource", desc);
}

void GEMDAQStatusSource::LoadROMap(edm::EventSetup const &iSetup) {
auto gemROMap = std::make_shared<GEMROMapping>();
//if (useDBEMap_)
if (true) {
const auto &eMap = iSetup.getData(gemEMapToken_);
auto gemEMap = std::make_unique<GEMeMap>(eMap);
gemEMap->convert(*gemROMap);

for (auto imap : gemEMap->theChamberMap_) {
int nNumChamber = (int)imap.fedId.size();
for (int i = 0; i < nNumChamber; i++) {
unsigned int fedId = imap.fedId[i];
uint8_t amcNum = imap.amcNum[i];
uint8_t gebId = imap.gebId[i];
GEMROMapping::chamEC geb_ec{fedId, amcNum, gebId};
GEMROMapping::chamDC geb_dc = gemROMap->chamberPos(geb_ec);
GEMDetId gemChId = geb_dc.detId;

mapFEDIdToRe_[fedId] = gemChId.region();
mapAMC13ToListChamber_[fedId].push_back(gemChId);
mapAMCToListChamber_[{fedId, amcNum}].push_back(gemChId);
}
}

gemEMap.reset();
} else {
// no EMap in DB, using dummy
// FIXME: How to add mapFEDIdToRe_ and mapDetIdToAMC_??
auto gemEMap = std::make_unique<GEMeMap>();
gemEMap->convertDummy(*gemROMap);
gemEMap.reset();
}
}

void GEMDAQStatusSource::SetLabelAMC13Status(MonitorElement *h2Status) {
unsigned int unBinPos = 1;
h2Status->setBinLabel(unBinPos++, "Good", 2);
Expand Down Expand Up @@ -91,6 +126,9 @@ void GEMDAQStatusSource::SetLabelVFATStatus(MonitorElement *h2Status) {
}

void GEMDAQStatusSource::bookHistograms(DQMStore::IBooker &ibooker, edm::Run const &, edm::EventSetup const &iSetup) {
LoadROMap(iSetup);
if (mapAMC13ToListChamber_.empty() || mapAMCToListChamber_.empty())
return;
initGeometry(iSetup);
if (GEMGeometry_ == nullptr)
return;
Expand Down Expand Up @@ -145,10 +183,24 @@ void GEMDAQStatusSource::bookHistograms(DQMStore::IBooker &ibooker, edm::Run con
h2SummaryStatusAll = CreateSummaryHist(ibooker, "chamberAllStatus");
h2SummaryStatusWarning = CreateSummaryHist(ibooker, "chamberWarnings");
h2SummaryStatusError = CreateSummaryHist(ibooker, "chamberErrors");
h2SummaryStatusVFATWarning = CreateSummaryHist(ibooker, "chamberVFATWarnings");
h2SummaryStatusVFATError = CreateSummaryHist(ibooker, "chamberVFATErrors");
h2SummaryStatusOHWarning = CreateSummaryHist(ibooker, "chamberOHWarnings");
h2SummaryStatusOHError = CreateSummaryHist(ibooker, "chamberOHErrors");
h2SummaryStatusAMCWarning = CreateSummaryHist(ibooker, "chamberAMCWarnings");
h2SummaryStatusAMCError = CreateSummaryHist(ibooker, "chamberAMCErrors");
h2SummaryStatusAMC13Error = CreateSummaryHist(ibooker, "chamberAMC13Errors");

h2SummaryStatusAll->setTitle("Summary of all number of OH or VFAT status of each chambers");
h2SummaryStatusWarning->setTitle("Summary of OH or VFAT warnings of each chambers");
h2SummaryStatusError->setTitle("Summary of OH or VFAT errors of each chambers");
h2SummaryStatusWarning->setTitle("Summary of all warnings of each chambers");
h2SummaryStatusError->setTitle("Summary of all errors of each chambers");
h2SummaryStatusVFATWarning->setTitle("Summary of VFAT warnings of each chambers");
h2SummaryStatusVFATError->setTitle("Summary of VFAT errors of each chambers");
h2SummaryStatusOHWarning->setTitle("Summary of OH warnings of each chambers");
h2SummaryStatusOHError->setTitle("Summary of OH errors of each chambers");
h2SummaryStatusAMCWarning->setTitle("Summary of AMC warnings of each chambers");
h2SummaryStatusAMCError->setTitle("Summary of AMC errors of each chambers");
h2SummaryStatusAMC13Error->setTitle("Summary of AMC13 errors of each chambers");
}

int GEMDAQStatusSource::ProcessWithMEMap3(BookingHelper &bh, ME3IdsKey key) {
Expand Down Expand Up @@ -204,6 +256,17 @@ void GEMDAQStatusSource::analyze(edm::Event const &event, edm::EventSetup const
event.getByToken(tagAMC_, gemAMC);
event.getByToken(tagAMC13_, gemAMC13);

std::map<ME4IdsKey, bool> mapChamberAll;
std::map<ME4IdsKey, bool> mapChamberWarning;
std::map<ME4IdsKey, bool> mapChamberError;
std::map<ME4IdsKey, bool> mapChamberVFATWarning;
std::map<ME4IdsKey, bool> mapChamberVFATError;
std::map<ME4IdsKey, bool> mapChamberOHWarning;
std::map<ME4IdsKey, bool> mapChamberOHError;
std::map<ME4IdsKey, bool> mapChamberAMCWarning;
std::map<ME4IdsKey, bool> mapChamberAMCError;
std::map<ME4IdsKey, bool> mapChamberAMC13Error;

for (auto amc13It = gemAMC13->begin(); amc13It != gemAMC13->end(); ++amc13It) {
int fedId = (*amc13It).first;
if (mapFEDIdToRe_.find(fedId) == mapFEDIdToRe_.end())
Expand Down Expand Up @@ -245,8 +308,16 @@ void GEMDAQStatusSource::analyze(edm::Event const &event, edm::EventSetup const
if (errors.wrongFedId)
FillWithRiseErr(h2AMC13Status_, nXBin, 10, bErr);

if (!bWarn && !bErr)
if (!bWarn && !bErr) {
h2AMC13Status_->Fill(nXBin, 1);
} else {
auto &listChamber = mapAMC13ToListChamber_[fedId];
for (auto gid : listChamber) {
ME4IdsKey key4Ch{gid.region(), gid.station(), gid.layer(), gid.chamber()};
if (bErr)
mapChamberAMC13Error[key4Ch] = false;
}
}
}
}

Expand Down Expand Up @@ -298,16 +369,21 @@ void GEMDAQStatusSource::analyze(edm::Event const &event, edm::EventSetup const
if (errors.BC0locked)
FillWithRiseErr(h2AMCStatus, nAMCNum, 12, bErr);

if (!bWarn && !bErr)
if (!bWarn && !bErr) {
h2AMCStatus->Fill(nAMCNum, 1);
} else {
auto &listChamber = mapAMCToListChamber_[{fedId, nAMCNum}];
for (auto gid : listChamber) {
ME4IdsKey key4Ch{gid.region(), gid.station(), gid.layer(), gid.chamber()};
if (bErr)
mapChamberAMCError[key4Ch] = false;
if (bWarn)
mapChamberAMCWarning[key4Ch] = false;
}
}
}
}

// WARNING: ME4IdsKey for region, station, layer, chamber (not iEta)
std::map<ME4IdsKey, bool> mapChamberAll;
std::map<ME4IdsKey, bool> mapChamberWarning;
std::map<ME4IdsKey, bool> mapChamberError;

for (auto ohIt = gemOH->begin(); ohIt != gemOH->end(); ++ohIt) {
GEMDetId gid = (*ohIt).first;
ME3IdsKey key3{gid.region(), gid.station(), gid.layer()};
Expand Down Expand Up @@ -356,9 +432,9 @@ void GEMDAQStatusSource::analyze(edm::Event const &event, edm::EventSetup const
if (!bWarn && !bErr)
mapStatusOH_.Fill(key3, gid.chamber(), 1);
if (bWarn)
mapChamberWarning[key4] = false;
mapChamberOHWarning[key4] = false;
if (bErr)
mapChamberError[key4] = false;
mapChamberOHError[key4] = false;
mapChamberAll[key4] = true;
}
}
Expand Down Expand Up @@ -407,9 +483,9 @@ void GEMDAQStatusSource::analyze(edm::Event const &event, edm::EventSetup const
if (!bWarn && !bErr)
mapStatusVFATPerCh_.Fill(key4Ch, nIdxVFAT, 1);
if (bWarn)
mapChamberWarning[key4Ch] = false;
mapChamberVFATWarning[key4Ch] = false;
if (bErr)
mapChamberError[key4Ch] = false;
mapChamberVFATError[key4Ch] = false;
if (bWarn)
mapStatusWarnVFATPerLayer_.Fill(key3, gid.chamber(), nIdxVFAT);
if (bErr)
Expand All @@ -425,21 +501,19 @@ void GEMDAQStatusSource::analyze(edm::Event const &event, edm::EventSetup const
h2SummaryStatusAll->Fill(nChamber, mapStationToIdx_[key3]);
}

// Summarizing the warning occupancy
for (auto const &[key4, bWarning] : mapChamberWarning) {
if (mapChamberError.find(key4) != mapChamberError.end()) // Avoiding any double-counting
continue;
ME3IdsKey key3 = key4Tokey3(key4);
Int_t nChamber = keyToChamber(key4);
h2SummaryStatusWarning->Fill(nChamber, mapStationToIdx_[key3]);
}

// Summarizing the error occupancy
for (auto const &[key4, bErr] : mapChamberError) {
ME3IdsKey key3 = key4Tokey3(key4);
Int_t nChamber = keyToChamber(key4);
h2SummaryStatusError->Fill(nChamber, mapStationToIdx_[key3]);
}
// Summarizing all presence of status of each chamber
FillStatusSummaryPlot(mapChamberAll, h2SummaryStatusAll);
// Summarizing all the error and warning occupancy
FillStatusSummaryPlot(mapChamberVFATWarning, h2SummaryStatusVFATWarning, &mapChamberWarning);
FillStatusSummaryPlot(mapChamberVFATError, h2SummaryStatusVFATError, &mapChamberError);
FillStatusSummaryPlot(mapChamberOHWarning, h2SummaryStatusOHWarning, &mapChamberWarning);
FillStatusSummaryPlot(mapChamberOHError, h2SummaryStatusOHError, &mapChamberError);
FillStatusSummaryPlot(mapChamberAMCWarning, h2SummaryStatusAMCWarning, &mapChamberWarning);
FillStatusSummaryPlot(mapChamberAMCError, h2SummaryStatusAMCError, &mapChamberError);
FillStatusSummaryPlot(mapChamberAMC13Error, h2SummaryStatusAMC13Error, &mapChamberError);

FillStatusSummaryPlot(mapChamberWarning, h2SummaryStatusWarning);
FillStatusSummaryPlot(mapChamberError, h2SummaryStatusError);
}

DEFINE_FWK_MODULE(GEMDAQStatusSource);