Merge pull request #114 from blinkseb/fix/electron_id

Fix HWW electron id
cp3-llbb · Nov 23, 2016 · f37d854 · f37d854
2 parents d217142 + 389080e
commit f37d854
Show file tree

Hide file tree

Showing 3 changed files with 9 additions and 26 deletions.
diff --git a/interface/HHAnalyzer.h b/interface/HHAnalyzer.h
@@ -39,6 +39,7 @@ class HHAnalyzer: public Framework::Analyzer {
             m_electron_loose_wp_name = config.getUntrackedParameter<std::string>("electrons_loose_wp_name");
             m_electron_medium_wp_name = config.getUntrackedParameter<std::string>("electrons_medium_wp_name");
             m_electron_tight_wp_name = config.getUntrackedParameter<std::string>("electrons_tight_wp_name");
+            m_electron_hlt_safe_wp_name = config.getUntrackedParameter<std::string>("electrons_hlt_safe_wp_name");
 
             m_jetEtaCut = config.getUntrackedParameter<double>("jetEtaCut", 2.4);
             m_jetPtCut = config.getUntrackedParameter<double>("jetPtCut", 20);
@@ -230,6 +231,7 @@ class HHAnalyzer: public Framework::Analyzer {
         std::string m_electron_loose_wp_name;
         std::string m_electron_medium_wp_name;
         std::string m_electron_tight_wp_name;
+        std::string m_electron_hlt_safe_wp_name;
         bool m_applyBJetRegression;
         std::unordered_map<std::string, std::unique_ptr<BinnedValues>> m_hlt_efficiencies;
 };

diff --git a/plugins/HHAnalyzer.cc b/plugins/HHAnalyzer.cc
@@ -65,45 +65,25 @@ void HHAnalyzer::analyze(const edm::Event& event, const edm::EventSetup&, const
     leptons.clear();
     ll.clear();
 
-    static auto electron_pass_HWW_id = [&allelectrons](size_t index) {
-        bool result = true;
+    static auto electron_pass_HWW_id = [&allelectrons, this](size_t index) {
         auto electron = allelectrons.products[index];
 
-        static auto dEtaInSeed = [&electron]() {
-            return electron->superCluster().isNonnull() && electron->superCluster()->seed().isNonnull() ? electron->deltaEtaSuperClusterTrackAtVtx() - electron->superCluster()->eta() + electron->superCluster()->seed()->eta() : std::numeric_limits<float>::max();
-        };
+        // Use POG HLT-safe id, which is the same id used by HWW
+
+        bool result = allelectrons.ids[index][m_electron_hlt_safe_wp_name];
 
-        // Cuts described at https://twiki.cern.ch/twiki/pub/CMS/HWW2016TriggerAndIdIsoScaleFactorsResults/AN-16-172_temp.pdf
+        // Add dxy and d0 cuts described at https://twiki.cern.ch/twiki/pub/CMS/HWW2016TriggerAndIdIsoScaleFactorsResults/AN-16-172_temp.pdf
         // page 24
         if (electron->isEB()) {
-            result &= electron->ecalPFClusterIso() < 0.160;
-            result &= electron->full5x5_sigmaIetaIeta() < 0.011;
-            result &= std::abs(dEtaInSeed()) < 0.004;
-            result &= std::abs(electron->deltaPhiSuperClusterTrackAtVtx()) < 0.020;
             result &= std::abs(allelectrons.dz[index]) < 0.373;
             result &= std::abs(allelectrons.dxy[index]) < 0.1;
         } else {
-            result &= electron->ecalPFClusterIso() < 0.120;
-            result &= electron->full5x5_sigmaIetaIeta() < 0.031;
             result &= std::abs(allelectrons.dz[index]) < 0.602;
             result &= std::abs(allelectrons.dxy[index]) < 0.2;
         }
 
-        // Common cuts to EE & EB
-        result &= electron->hcalPFClusterIso() < 0.120;
-        result &= electron->trackIso() < 0.08;
-        result &= electron->hadronicOverEm() < 0.060;
-
-        // Code from https://github.com/ikrav/cmssw/blob/egm_id_80X_v1/RecoEgamma/ElectronIdentification/plugins/cuts/GsfEleEInverseMinusPInverseCut.cc#L43-L45
-        const float ecal_energy_inverse = 1.0 / electron->ecalEnergy();
-        const float eSCoverP = electron->eSuperClusterOverP();
-        const float ooEmooP = std::abs(1.0 - eSCoverP) * ecal_energy_inverse;
-        result &= ooEmooP < 0.013;
-
         result &= (electron->gsfTrack()->hitPattern().numberOfHits(reco::HitPattern::MISSING_INNER_HITS)) < 1;
 
-        // Pass conversion veto is already part of all electron IDs
-
         return result;
     };
 
@@ -128,7 +108,7 @@ void HHAnalyzer::analyze(const edm::Event& event, const edm::EventSetup&, const
 
             // For electrons, isolation requirement is already included in ID
             ele.iso_L = ele.id_L;
-            ele.iso_T = ele.iso_T;
+            ele.iso_T = ele.id_T;
             ele.iso_HWW = ele.iso_T;
 
             ele.gen_matched = allelectrons.matched[ielectron];

diff --git a/test/HHConfiguration.py b/test/HHConfiguration.py
@@ -53,6 +53,7 @@
             electrons_loose_wp_name = cms.untracked.string("cutBasedElectronID-Summer16-80X-V1-loose"),
             electrons_medium_wp_name = cms.untracked.string("cutBasedElectronID-Summer16-80X-V1-medium"),
             electrons_tight_wp_name = cms.untracked.string("cutBasedElectronID-Summer16-80X-V1-tight"),
+            electrons_hlt_safe_wp_name = cms.untracked.string("cutBasedElectronHLTPreselection-Summer16-V1"),
             jetEtaCut = cms.untracked.double(2.4),
             jetPtCut = cms.untracked.double(20),
             # BTAG INFO