Skip to content

Commit

Permalink
Make ML plots available only for online dqm
Browse files Browse the repository at this point in the history
  • Loading branch information
abhih1 committed Nov 8, 2021
1 parent 603d1aa commit 1de34de
Show file tree
Hide file tree
Showing 2 changed files with 36 additions and 12 deletions.
10 changes: 10 additions & 0 deletions DQM/EcalMonitorClient/src/MLClient.cc
Expand Up @@ -22,9 +22,19 @@ namespace ecaldqm {
PUcorr_slope_ = _params.getUntrackedParameter<double>("PUcorr_slope");
PUcorr_intercept_ = _params.getUntrackedParameter<double>("PUcorr_intercept");
avgOcc_ = _params.getUntrackedParameter<std::vector<double>>("avgOcc");
if (!onlineMode_) {
MEs_.erase(std::string("MLQualitySummary"));
MEs_.erase(std::string("EventsperMLImage"));
sources_.erase(std::string("PU"));
sources_.erase(std::string("NumEvents"));
sources_.erase(std::string("DigiAllByLumi"));
sources_.erase(std::string("AELoss"));
}
}

void MLClient::producePlots(ProcessType) {
if (!onlineMode_)
return;
using namespace std;
MESet& meMLQualitySummary(MEs_.at("MLQualitySummary"));
MESet& meEventsperMLImage(MEs_.at("EventsperMLImage"));
Expand Down
38 changes: 26 additions & 12 deletions DQM/EcalMonitorTasks/src/OccupancyTask.cc
Expand Up @@ -15,6 +15,13 @@ namespace ecaldqm {
tpThreshold_ = _params.getUntrackedParameter<double>("tpThreshold");
lumiTag = _params.getParameter<edm::InputTag>("scalers");
lumiCheck_ = _params.getUntrackedParameter<bool>("lumiCheck", false);
if (!onlineMode_) {
MEs_.erase(std::string("PU"));
MEs_.erase(std::string("NEvents"));
MEs_.erase(std::string("TrendEventsperLumi"));
MEs_.erase(std::string("TrendPUperLumi"));
MEs_.erase(std::string("AELoss"));
}
}

void OccupancyTask::setTokens(edm::ConsumesCollector& _collector) {
Expand Down Expand Up @@ -43,10 +50,12 @@ namespace ecaldqm {
MEs_.at("DigiAllByLumi").reset(GetElectronicsMap());
MEs_.at("TPDigiThrAllByLumi").reset(GetElectronicsMap());
MEs_.at("RecHitThrAllByLumi").reset(GetElectronicsMap());
MEs_.at("PU").reset(GetElectronicsMap(), -1);
MEs_.at("NEvents").reset(GetElectronicsMap(), -1);
nEv = 0;
FindPUinLS = true;
if (onlineMode_) {
MEs_.at("PU").reset(GetElectronicsMap(), -1);
MEs_.at("NEvents").reset(GetElectronicsMap(), -1);
FindPUinLS = true;
}
}
nEv++;
MESet& meLaserCorrProjEta(MEs_.at("LaserCorrProjEta"));
Expand Down Expand Up @@ -95,13 +104,15 @@ namespace ecaldqm {
}

void OccupancyTask::endLuminosityBlock(edm::LuminosityBlock const&, edm::EventSetup const&) {
MESet& meNEvents(static_cast<MESet&>(MEs_.at("NEvents")));
MESet& meTrendEventsperLumi(MEs_.at("TrendEventsperLumi"));
MESet& meTrendPUperLumi(MEs_.at("TrendPUperLumi"));

meNEvents.fill(getEcalDQMSetupObjects(), double(nEv));
meTrendEventsperLumi.fill(getEcalDQMSetupObjects(), EcalBarrel, double(timestamp_.iLumi), double(nEv));
meTrendPUperLumi.fill(getEcalDQMSetupObjects(), EcalBarrel, double(timestamp_.iLumi), double(scal_pu));
if (onlineMode_) {
MESet& meNEvents(static_cast<MESet&>(MEs_.at("NEvents")));
MESet& meTrendEventsperLumi(MEs_.at("TrendEventsperLumi"));
MESet& meTrendPUperLumi(MEs_.at("TrendPUperLumi"));

meNEvents.fill(getEcalDQMSetupObjects(), double(nEv));
meTrendEventsperLumi.fill(getEcalDQMSetupObjects(), EcalBarrel, double(timestamp_.iLumi), double(nEv));
meTrendPUperLumi.fill(getEcalDQMSetupObjects(), EcalBarrel, double(timestamp_.iLumi), double(scal_pu));
}
}

template <typename DigiCollection>
Expand All @@ -114,7 +125,9 @@ namespace ecaldqm {
MESet& meDigiDCC(MEs_.at("DigiDCC"));
MESet& meDigi1D(MEs_.at("Digi1D"));
MESet& meTrendNDigi(MEs_.at("TrendNDigi"));
MESet& meAELoss(MEs_.at("AELoss"));
MESet* meAELoss = nullptr;
if (onlineMode_)
meAELoss = &MEs_.at("AELoss");

std::for_each(_digis.begin(), _digis.end(), [&](typename DigiCollection::Digi const& digi) {
DetId id(digi.id());
Expand All @@ -124,7 +137,8 @@ namespace ecaldqm {
meDigiAll.fill(getEcalDQMSetupObjects(), id);
meDigiAllByLumi.fill(getEcalDQMSetupObjects(), id);
meDigiDCC.fill(getEcalDQMSetupObjects(), id);
meAELoss.fill(getEcalDQMSetupObjects(), id);
if (onlineMode_)
meAELoss->fill(getEcalDQMSetupObjects(), id);
});

int iSubdet(_collection == kEBDigi ? EcalBarrel : EcalEndcap);
Expand Down

0 comments on commit 1de34de

Please sign in to comment.