Merge pull request #191 from susy2015/resTopPtCut

Add upper cut option on resolved top Pt
susy2015 · Feb 16, 2020 · 8752555 · 8752555
2 parents 036253a + 5439074
commit 8752555
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/TopTagger/interface/TTMBasicClusterAlgo.h b/TopTagger/interface/TTMBasicClusterAlgo.h
@@ -38,7 +38,7 @@ class TTMBasicClusterAlgo : public TTModule, public TTMConstituentReqs
     bool doDijet_;
 
     //trijet variables
-    double minTopCandMass_, maxTopCandMass_, dRMaxTrijet_, minTrijetAK4JetPt_, midTrijetAK4JetPt_, maxTrijetAK4JetPt_;
+    double minTopCandMass_, maxTopCandMass_, dRMaxTrijet_, minTrijetAK4JetPt_, midTrijetAK4JetPt_, maxTrijetAK4JetPt_, maxTrijetPt_;
     bool doTrijet_;
     int nbSeed_;
 

diff --git a/TopTagger/src/TTMBasicClusterAlgo.cpp b/TopTagger/src/TTMBasicClusterAlgo.cpp
@@ -30,6 +30,7 @@ void TTMBasicClusterAlgo::getParameters(const cfg::CfgDocument* cfgDoc, const st
     minTrijetAK4JetPt_  = cfgDoc->get("minTrijetAK4JetPt", localCxt, -999.0);
     midTrijetAK4JetPt_  = cfgDoc->get("midTrijetAK4JetPt", localCxt, -999.0);
     maxTrijetAK4JetPt_  = cfgDoc->get("maxTrijetAK4JetPt", localCxt, -999.0);
+    maxTrijetPt_        = cfgDoc->get("maxTrijetPt",       localCxt, -999.0);
 
     //get vars for TTMConstituentReqs
     TTMConstituentReqs::getParameters(cfgDoc, localContextName);
@@ -161,8 +162,9 @@ void TTMBasicClusterAlgo::fillTriplet(const Constituent* const c1, const Constit
     //mass window on the top candidate mass
     double m123 = topCand.p().M();
     bool passMassWindow = (minTopCandMass_ < m123) && (m123 < maxTopCandMass_);
+    bool passMaxPt = (maxTrijetPt_ < 0.0) || topCand.p().Pt() < maxTrijetPt_;
 
-    if(topCand.getDRmax() < dRMaxTrijet_ && passMassWindow)
+    if(topCand.getDRmax() < dRMaxTrijet_ && passMassWindow && passMaxPt)
     {
         topCandidates.push_back(topCand);
     }