Add Eta cut to HLT Egamma objects.

Will squash soon. Co-authored-by: Marino Missiroli <m.missiroli@cern.ch>
cms-sw · Oct 20, 2022 · 4877938 · 4877938
1 parent 0fb364f
commit 4877938
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 3 deletions.
diff --git a/HLTrigger/Egamma/plugins/HLTEgammaEtFilter.cc b/HLTrigger/Egamma/plugins/HLTEgammaEtFilter.cc
@@ -22,6 +22,8 @@ HLTEgammaEtFilter::HLTEgammaEtFilter(const edm::ParameterSet& iConfig) : HLTFilt
   inputTag_ = iConfig.getParameter<edm::InputTag>("inputTag");
   etcutEB_ = iConfig.getParameter<double>("etcutEB");
   etcutEE_ = iConfig.getParameter<double>("etcutEE");
+  minEtaCut_ = iConfig.getParameter<double>("minEtaCut");
+  maxEtaCut_ = iConfig.getParameter<double>("maxEtaCut");
   ncandcut_ = iConfig.getParameter<int>("ncandcut");
   l1EGTag_ = iConfig.getParameter<edm::InputTag>("l1EGCand");
   inputToken_ = consumes<trigger::TriggerFilterObjectWithRefs>(inputTag_);
@@ -34,6 +36,8 @@ void HLTEgammaEtFilter::fillDescriptions(edm::ConfigurationDescriptions& descrip
   desc.add<edm::InputTag>("l1EGCand", edm::InputTag("hltL1IsoRecoEcalCandidate"));
   desc.add<double>("etcutEB", 1.0);
   desc.add<double>("etcutEE", 1.0);
+  desc.add<double>("minEtaCut", -9999.0);
+  desc.add<double>("maxEtaCut", 9999.0);
   desc.add<int>("ncandcut", 1);
   descriptions.add("hltEgammaEtFilter", desc);
 }
@@ -72,6 +76,9 @@ bool HLTEgammaEtFilter::hltFilter(edm::Event& iEvent,
   for (auto& recoecalcand : recoecalcands) {
     ref = recoecalcand;
 
+    if ((ref->eta() < minEtaCut_) or (ref->eta() > maxEtaCut_))
+      continue;
+
     if ((fabs(ref->eta()) < 1.479 && ref->et() >= etcutEB_) || (fabs(ref->eta()) >= 1.479 && ref->et() >= etcutEE_)) {
       n++;
       // std::cout << "Passed eta: " << ref->eta() << std::endl;

diff --git a/HLTrigger/Egamma/plugins/HLTEgammaEtFilter.h b/HLTrigger/Egamma/plugins/HLTEgammaEtFilter.h
@@ -31,9 +31,11 @@ class HLTEgammaEtFilter : public HLTFilter {
 private:
   edm::InputTag inputTag_;  // input tag identifying product contains egammas
   edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> inputToken_;
-  double etcutEB_;  // Barrel Et threshold in GeV
-  double etcutEE_;  // Endcap Et threshold in GeV
-  int ncandcut_;    // number of egammas required
+  double etcutEB_;    // Barrel Et threshold in GeV
+  double etcutEE_;    // Endcap Et threshold in GeV
+  double minEtaCut_;  // Min pseudorapidity cut
+  double maxEtaCut_;  // Max pseudorapidity cut
+  int ncandcut_;      // number of egammas required
 
   edm::InputTag l1EGTag_;
 };