[Rivet-svn] r1888 - trunk/src/Analyses

[blackhole at projects.hepforge.org]

Author: buckley
Date: Thu Oct  8 11:29:05 2009
New Revision: 1888

Log:
Adding SUSY m_ll lepton pT cuts

Modified:
   trunk/src/Analyses/MC_LHC_SUSY.cc

Modified: trunk/src/Analyses/MC_LHC_SUSY.cc
==============================================================================
--- trunk/src/Analyses/MC_LHC_SUSY.cc	Wed Oct  7 19:00:11 2009	(r1887)
+++ trunk/src/Analyses/MC_LHC_SUSY.cc	Thu Oct  8 11:29:05 2009	(r1888)
@@ -146,7 +146,10 @@
         _hist_phi_e->fill(mapAngleMPiToPi(p.phi()), weight);
         _hist_eta_e->fill(p.eta(), weight);
         _hist_pt_e->fill(p.pT()/GeV, weight);
-        if (PID::threeCharge(e.pdgId()) > 0) epluses += p; else eminuses += p;
+        // Add sufficiently hard leptons to collections for m_ll histo
+        if (p.pT()/GeV > 20) {
+          if (PID::threeCharge(e.pdgId()) > 0) epluses += p; else eminuses += p;
+        }
       }
 
       /// @todo Resum photons around muons
@@ -160,7 +163,10 @@
         _hist_phi_mu->fill(mapAngleMPiToPi(p.phi()), weight);
         _hist_eta_mu->fill(p.eta(), weight);
         _hist_pt_mu->fill(p.pT()/GeV, weight);
-        if (PID::threeCharge(mu.pdgId()) > 0) mupluses += p; else muminuses += p;
+        // Add sufficiently hard leptons to collections for m_ll histo
+        if (p.pT()/GeV > 20) {
+          if (PID::threeCharge(mu.pdgId()) > 0) mupluses += p; else muminuses += p;
+        }
       }
 
       // Fill final state non-isolated photon histos