From 30607cfa2bc3ede9f8f1327aec59b3af18f59d04 Mon Sep 17 00:00:00 2001
From: Marco Link <marco.link@student.kit.edu>
Date: Wed, 19 Jun 2019 10:31:25 +0200
Subject: [PATCH] changed tHW cuts
---
src/MVAvarsJABDTthw.cpp | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/src/MVAvarsJABDTthw.cpp b/src/MVAvarsJABDTthw.cpp
index 9a63855..f9867ab 100644
--- a/src/MVAvarsJABDTthw.cpp
+++ b/src/MVAvarsJABDTthw.cpp
@@ -252,14 +252,16 @@ bool MVAvarsJABDTthw::SkipEvent(const std::vector<TLorentzVector> &selectedJetP4
const std::vector<int> &jets_idx)
{
int btop = jets_idx.at(tHWIndexes::tHW_btop_idx);
-// int whaddau1 = jets_idx.at(tHWIndexes::tHW_whaddau1_idx);
-// int whaddau2 = jets_idx.at(tHWIndexes::tHW_whaddau2_idx);
+ int whaddau1 = jets_idx.at(tHWIndexes::tHW_whaddau1_idx);
+ int whaddau2 = jets_idx.at(tHWIndexes::tHW_whaddau2_idx);
int hdau1 = jets_idx.at(tHWIndexes::tHW_hdau1_idx);
int hdau2 = jets_idx.at(tHWIndexes::tHW_hdau2_idx);
if (!JetIsCentral(selectedJetP4[btop])) return true;
if (!JetIsCentral(selectedJetP4[hdau1])) return true;
if (!JetIsCentral(selectedJetP4[hdau2])) return true;
+ if (JetIsTagged(selectedJetP4[whaddau1], selectedJetCSV.at(whaddau1))) return true;
+ if (JetIsTagged(selectedJetP4[whaddau2], selectedJetCSV.at(whaddau2))) return true;
int btags = 0;
if(JetIsTagged(selectedJetP4.at(btop), selectedJetCSV.at(btop))) btags++;
--
GitLab