Moved module declaration to SealModule, added checking for empty region or empty CICADAScore

Author: Max Zhao

DQM/L1TMonitor/plugins/SealModule.cc

@@ -27,6 +27,9 @@ DEFINE_FWK_MODULE(L1TGMT);
 #include "DQM/L1TMonitor/interface/L1TStage2CaloLayer1.h"
 DEFINE_FWK_MODULE(L1TStage2CaloLayer1);
 
+#include "DQM/L1TMonitor/interface/L1TCaloLayer1Summary.h"
+DEFINE_FWK_MODULE(L1TCaloLayer1Summary);
+
 #include "DQM/L1TMonitor/interface/L1TStage2CaloLayer2.h"
 DEFINE_FWK_MODULE(L1TStage2CaloLayer2);
 

DQM/L1TMonitor/src/L1TCaloLayer1Summary.cc

@@ -9,7 +9,8 @@ L1TCaloLayer1Summary::L1TCaloLayer1Summary(const edm::ParameterSet& iConfig)
           consumes<L1CaloRegionCollection>(iConfig.getParameter<edm::InputTag>("caloLayer1Regions"))),
       simRegionsToken_(consumes<L1CaloRegionCollection>(iConfig.getParameter<edm::InputTag>("simRegions"))),
       fedRawData_(consumes<FEDRawDataCollection>(iConfig.getParameter<edm::InputTag>("fedRawDataLabel"))),
-      histFolder_(iConfig.getParameter<std::string>("histFolder")) {}
+      histFolder_(iConfig.getParameter<std::string>("histFolder")) {
+}
 
 // ------------ method called for each event  ------------
 void L1TCaloLayer1Summary::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
@@ -77,15 +78,21 @@ void L1TCaloLayer1Summary::analyze(const edm::Event& iEvent, const edm::EventSet
     }
   }
 
-  float caloCICADAScore = iEvent.get(caloLayer1CICADAScoreToken_)[0];
+  auto caloCICADAScores = iEvent.get(caloLayer1CICADAScoreToken_);
   auto gtCICADAScores = iEvent.get(gtCICADAScoreToken_);
-  float simCICADAScore = iEvent.get(simCICADAScoreToken_)[0];
+  auto simCICADAScores = iEvent.get(simCICADAScoreToken_);
 
-  histoSimCICADAScore->Fill(simCICADAScore);
-  histoCaloMinusSim->Fill(caloCICADAScore - simCICADAScore);
-  histoCaloLayer1CICADAScore->Fill(caloCICADAScore);
-  histoGtCICADAScore->Fill(gtCICADAScores.at(0, 0));
-  histoCaloMinusGt->Fill(gtCICADAScores.at(0, 0) - caloCICADAScore);
+  if (caloCICADAScores.size() > 0) {
+    histoCaloLayer1CICADAScore->Fill(caloCICADAScores[0]);
+    if (gtCICADAScores.size() > 0) {
+      histoGtCICADAScore->Fill(gtCICADAScores.at(0, 0));
+      histoCaloMinusGt->Fill(caloCICADAScores[0] - gtCICADAScores.at(0, 0));
+    }
+    if (simCICADAScores.size() > 0) {
+      histoSimCICADAScore->Fill(simCICADAScores[0]);
+      histoCaloMinusSim->Fill(caloCICADAScores[0] - simCICADAScores[0]);
+    }
+  }
 }
 
 void L1TCaloLayer1Summary::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const&, edm::EventSetup const&) {
@@ -119,7 +126,4 @@ void L1TCaloLayer1Summary::fillDescriptions(edm::ConfigurationDescriptions& desc
   desc.add<edm::InputTag>("fedRawDataLabel", edm::InputTag("rawDataCollector"));
   desc.add<std::string>("histFolder", "L1T/L1TCaloLayer1Summary");
   descriptions.add("l1tCaloLayer1Summary", desc);
-}
-
-//define this as a plug-in
-DEFINE_FWK_MODULE(L1TCaloLayer1Summary);
+}

L1Trigger/L1TCaloLayer1/plugins/L1TCaloSummary.cc

@@ -107,6 +107,7 @@ class L1TCaloSummary : public edm::stream::EDProducer<> {
   int fwVersion;
 
   edm::EDGetTokenT<L1CaloRegionCollection> regionToken;
+  edm::EDGetTokenT<L1CaloRegionCollection> backupRegionToken;
 
   UCTLayer1* layer1;
 
@@ -142,6 +143,7 @@ L1TCaloSummary<INPUT, OUTPUT>::L1TCaloSummary(const edm::ParameterSet& iConfig)
       verbose(iConfig.getParameter<bool>("verbose")),
       fwVersion(iConfig.getParameter<int>("firmwareVersion")),
       regionToken(consumes<L1CaloRegionCollection>(iConfig.getParameter<edm::InputTag>("caloLayer1Regions"))),
+      backupRegionToken(consumes<L1CaloRegionCollection>(edm::InputTag("simCaloStage2Layer1Digis"))),
       loader(hls4mlEmulator::ModelLoader(iConfig.getParameter<string>("CICADAModelVersion"))),
       overwriteWithTestPatterns(iConfig.getParameter<bool>("useTestPatterns")),
       testPatterns(iConfig.getParameter<std::vector<edm::ParameterSet>>("testPatterns")) {
@@ -198,6 +200,12 @@ void L1TCaloSummary<INPUT, OUTPUT>::produce(edm::Event& iEvent, const edm::Event
   if (!iEvent.getByToken(regionToken, regionCollection))
     edm::LogError("L1TCaloSummary") << "UCT: Failed to get regions from region collection!";
   iEvent.getByToken(regionToken, regionCollection);
+
+  if (regionCollection->size() == 0) {
+    iEvent.getByToken(backupRegionToken, regionCollection);
+    edm::LogWarning("L1TCaloSummary") << "Switched to emulated regions since data regions was empty.\n";
+  }
+
   //Model input
   //This is done as a flat vector input, but future versions may involve 2D input
   //This will have to be handled later