fixing issue #19590

cms-sw · Jul 18, 2017 · 4d41ee9 · 4d41ee9
1 parent 9d873f4
commit 4d41ee9
Show file tree

Hide file tree

Showing 3 changed files with 14 additions and 2 deletions.
diff --git a/DQMOffline/L1Trigger/interface/L1TStage2CaloLayer2Offline.h b/DQMOffline/L1Trigger/interface/L1TStage2CaloLayer2Offline.h
@@ -123,6 +123,9 @@ class L1TStage2CaloLayer2Offline: public DQMEDAnalyzer {
   std::vector<double> ettEfficiencyBins_;
   std::vector<double> httEfficiencyBins_;
 
+  double recoHTTMaxEta_;
+  double recoMHTMaxEta_;
+
   // TODO: add turn-on cuts (vectors of doubles)
   // Histograms
   MonitorElement* h_nVertex_;

diff --git a/DQMOffline/L1Trigger/python/L1TStage2CaloLayer2Offline_cfi.py b/DQMOffline/L1Trigger/python/L1TStage2CaloLayer2Offline_cfi.py
@@ -71,6 +71,9 @@
     mhtEfficiencyBins=cms.vdouble(mhtEfficiencyBins),
     ettEfficiencyBins=cms.vdouble(ettEfficiencyBins),
     httEfficiencyBins=cms.vdouble(httEfficiencyBins),
+
+    recoHTTMaxEta=cms.double(2.5),
+    recoMHTMaxEta=cms.double(2.5),
 )
 
 l1tStage2CaloLayer2OfflineDQMEmu = l1tStage2CaloLayer2OfflineDQM.clone(

diff --git a/DQMOffline/L1Trigger/src/L1TStage2CaloLayer2Offline.cc b/DQMOffline/L1Trigger/src/L1TStage2CaloLayer2Offline.cc
@@ -39,6 +39,8 @@ L1TStage2CaloLayer2Offline::L1TStage2CaloLayer2Offline(const edm::ParameterSet&
         mhtEfficiencyBins_(ps.getParameter < std::vector<double> > ("mhtEfficiencyBins")),
         ettEfficiencyBins_(ps.getParameter < std::vector<double> > ("ettEfficiencyBins")),
         httEfficiencyBins_(ps.getParameter < std::vector<double> > ("httEfficiencyBins")),
+        recoHTTMaxEta_(ps.getParameter <double>("recoHTTMaxEta")),
+        recoMHTMaxEta_(ps.getParameter <double>("recoMHTMaxEta")),
         h_controlPlots_()
 {
   edm::LogInfo("L1TStage2CaloLayer2Offline") << "Constructor "
@@ -174,8 +176,12 @@ void L1TStage2CaloLayer2Offline::fillEnergySums(edm::Event const& e, const unsig
       continue;
     }
     TVector2 jetVec(et * cos(jet->phi()), et * sin(jet->phi()));
-    recoHTT += et;
-    mht -= jetVec;
+    if(std::abs(jet->eta()) < recoHTTMaxEta_){
+      recoHTT += et;
+    }
+    if(std::abs(jet->eta()) < recoMHTMaxEta_){
+      mht -= jetVec;
+    }
   }
   recoETT = recoHTT;
   recoMHT = mht.Mod();