L1TauTriggerFilterObjectProducer: added addCollectionTag and minHwIso

Author: kandrosov

RecoTauTag/HLTProducers/src/L1TauTriggerFilterObjectProducer.cc

@@ -9,9 +9,11 @@ class L1TauTriggerFilterObjectProducer : public HLTFilter {
 public:
   explicit L1TauTriggerFilterObjectProducer(const edm::ParameterSet& cfg)
       : HLTFilter(cfg),
-        tausToken_(consumes<l1t::TauBxCollection>(cfg.getParameter<edm::InputTag>("taus"))),
+        tausSrc_(cfg.getParameter<edm::InputTag>("taus")),
+        tausToken_(consumes<l1t::TauBxCollection>(tausSrc_)),
         selectedBx_(cfg.getParameter<std::vector<int>>("selectedBx")),
         minPt_(cfg.getParameter<double>("minPt")),
+        minHwIso_(cfg.getParameter<int>("minHwIso")),
         nExpected_(cfg.getParameter<int>("nExpected")) {}
 
   static void fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
@@ -21,13 +23,16 @@ class L1TauTriggerFilterObjectProducer : public HLTFilter {
     desc.add<std::vector<int>>("selectedBx", std::vector<int>())
         ->setComment("bunch crossings to select, empty means all");
     desc.add<double>("minPt", 0)->setComment("select taus with pt > minPt. minPt=0 means no pt cut");
+    desc.add<int>("minHwIso", 0)->setComment("select taus with hwIso >= minHwIso. minHwIso=0 means no hwIso cut");
     desc.add<int>("nExpected", 0)->setComment("minimal number of taus per event to pass the filter");
     descriptions.addWithDefaultLabel(desc);
   }
 
   bool hltFilter(edm::Event& event,
                  const edm::EventSetup& eventsetup,
                  trigger::TriggerFilterObjectWithRefs& filterproduct) const override {
+    if (saveTags())
+      filterproduct.addCollectionTag(tausSrc_);
     const auto& taus = event.getHandle(tausToken_);
     int nPassed = 0;
     for (int bx_index = taus->getFirstBX(); bx_index <= taus->getLastBX(); ++bx_index) {
@@ -36,7 +41,7 @@ class L1TauTriggerFilterObjectProducer : public HLTFilter {
       const unsigned bx_index_shift = taus->begin(bx_index) - taus->begin();
       unsigned index_in_bx = 0;
       for (auto it = taus->begin(bx_index); it != taus->end(bx_index); ++it, ++index_in_bx) {
-        if (it->pt() <= minPt_)
+        if (it->pt() <= minPt_ || it->hwIso() < minHwIso_)
           continue;
         const l1t::TauRef tauRef(taus, bx_index_shift + index_in_bx);
         filterproduct.addObject(trigger::TriggerL1Tau, tauRef);
@@ -47,9 +52,11 @@ class L1TauTriggerFilterObjectProducer : public HLTFilter {
   }
 
 private:
+  const edm::InputTag tausSrc_;
   const edm::EDGetTokenT<l1t::TauBxCollection> tausToken_;
   const std::vector<int> selectedBx_;
   const double minPt_;
+  const int minHwIso_;
   const int nExpected_;
 };