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

[blackhole at projects.hepforge.org]

Author: fsiegert
Date: Tue Jun 22 17:55:33 2010
New Revision: 2503

Log:
OPAL_2001 (jet angles): Implement missing cut on jet energy, although it
probably won't make much difference for ycut=0.008 jets.

Modified:
   trunk/src/Analyses/OPAL_2001_S4553896.cc

Modified: trunk/src/Analyses/OPAL_2001_S4553896.cc
==============================================================================
--- trunk/src/Analyses/OPAL_2001_S4553896.cc	Tue Jun 22 13:43:25 2010	(r2502)
+++ trunk/src/Analyses/OPAL_2001_S4553896.cc	Tue Jun 22 17:55:33 2010	(r2503)
@@ -108,9 +108,12 @@
       }
 
       const FastJets& fastjets = applyProjection<FastJets>(event, "Jets");
-      vector<fastjet::PseudoJet> jets;
       if (fastjets.clusterSeq()) {
-        jets = fastjet::sorted_by_E(fastjets.clusterSeq()->exclusive_jets_ycut(0.008));
+        vector<fastjet::PseudoJet> jets;
+        foreach (const fastjet::PseudoJet& jet,
+                 fastjet::sorted_by_E(fastjets.clusterSeq()->exclusive_jets_ycut(0.008))) {
+          if (jet.E()>3.0*GeV) jets.push_back(jet);
+        }
         if (jets.size() == 4) {
           _h_BZ->fill(fabs(calc_BZ(jets)), weight);
           _h_KSW->fill(calc_KSW(jets), weight);