Add upper cut on jet time to HLTJetTimingFilter

cms-sw · Mar 24, 2023 · f43f616 · f43f616
1 parent 23b415a
commit f43f616
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/HLTrigger/JetMET/plugins/HLTJetTimingFilter.h b/HLTrigger/JetMET/plugins/HLTJetTimingFilter.h
@@ -37,6 +37,7 @@ class HLTJetTimingFilter : public HLTFilter {
   // Thresholds for selection
   const unsigned int minJets_;
   const double jetTimeThresh_;
+  const double jetMaxTimeThresh_;
   const double jetEcalEtForTimingThresh_;
   const unsigned int jetCellsForTimingThresh_;
   const double minPt_;
@@ -54,6 +55,7 @@ HLTJetTimingFilter<T>::HLTJetTimingFilter(const edm::ParameterSet& iConfig)
           consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("jetEcalEtForTiming"))},
       minJets_{iConfig.getParameter<unsigned int>("minJets")},
       jetTimeThresh_{iConfig.getParameter<double>("jetTimeThresh")},
+      jetMaxTimeThresh_{iConfig.getParameter<double>("jetMaxTimeThresh")},
       jetEcalEtForTimingThresh_{iConfig.getParameter<double>("jetEcalEtForTimingThresh")},
       jetCellsForTimingThresh_{iConfig.getParameter<unsigned int>("jetCellsForTimingThresh")},
       minPt_{iConfig.getParameter<double>("minJetPt")} {}
@@ -74,7 +76,7 @@ bool HLTJetTimingFilter<T>::hltFilter(edm::Event& iEvent,
   for (auto iterJet = jets->begin(); iterJet != jets->end(); ++iterJet) {
     edm::Ref<std::vector<T>> const caloJetRef(jets, std::distance(jets->begin(), iterJet));
     if (iterJet->pt() > minPt_ and jetTimes[caloJetRef] > jetTimeThresh_ and
-        jetEcalEtForTiming[caloJetRef] > jetEcalEtForTimingThresh_ and
+        jetTimes[caloJetRef] < jetMaxTimeThresh_ and jetEcalEtForTiming[caloJetRef] > jetEcalEtForTimingThresh_ and
         jetCellsForTiming[caloJetRef] > jetCellsForTimingThresh_) {
       // add caloJetRef to the event
       filterproduct.addObject(trigger::TriggerJet, caloJetRef);
@@ -97,6 +99,7 @@ void HLTJetTimingFilter<T>::fillDescriptions(edm::ConfigurationDescriptions& des
                           edm::InputTag("hltDisplacedHLTCaloJetCollectionProducerMidPtTiming", "jetEcalEtForTiming"));
   desc.add<unsigned int>("minJets", 1);
   desc.add<double>("jetTimeThresh", 1.);
+  desc.add<double>("jetMaxTimeThresh", 999999);
   desc.add<unsigned int>("jetCellsForTimingThresh", 5);
   desc.add<double>("jetEcalEtForTimingThresh", 10.);
   desc.add<double>("minJetPt", 40.);