[Rivet-svn] r4079 - in branches/2012-06-aidarivet: . src/Analyses

[blackhole at projects.hepforge.org]

Author: buckley
Date: Wed Dec 12 15:08:58 2012
New Revision: 4079

Log:
Correction to jet acceptance in CMS_2011_S9120041, from Sercan Sen: thanks!

Modified:
   branches/2012-06-aidarivet/ChangeLog
   branches/2012-06-aidarivet/src/Analyses/CMS_2011_S9120041.cc

Modified: branches/2012-06-aidarivet/ChangeLog
==============================================================================
--- branches/2012-06-aidarivet/ChangeLog	Wed Dec 12 13:26:31 2012	(r4078)
+++ branches/2012-06-aidarivet/ChangeLog	Wed Dec 12 15:08:58 2012	(r4079)
@@ -1,3 +1,7 @@
+2012-12-12  Andy Buckley  <andy.buckley at cern.ch>
+
+	* Correction to jet acceptance in CMS_2011_S9120041, from Sercan Sen: thanks!
+
 2012-12-12  Hendrik Hoeth <hendrik.hoeth at cern.ch>
 
 	* Added CMS_2012_PAS_QCD_11_010 analysis

Modified: branches/2012-06-aidarivet/src/Analyses/CMS_2011_S9120041.cc
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/CMS_2011_S9120041.cc	Wed Dec 12 13:26:31 2012	(r4078)
+++ branches/2012-06-aidarivet/src/Analyses/CMS_2011_S9120041.cc	Wed Dec 12 15:08:58 2012	(r4079)
@@ -11,6 +11,7 @@
 
 namespace Rivet {
 
+
   // UE charged particles vs. leading jet
   class CMS_2011_S9120041 : public Analysis {
   public:
@@ -24,8 +25,6 @@
       addProjection(cfs, "CFS");
 
       const ChargedFinalState cfsforjet(-2.5, 2.5, 500*MeV);
-      addProjection(cfsforjet, "CFSforjet");
-
       const FastJets jetpro(cfsforjet, FastJets::SISCONE, 0.5);
       addProjection(jetpro, "Jets");
 
@@ -60,10 +59,15 @@
     void analyze(const Event& event) {
       const double weight = event.weight();
 
-      Jets jets = applyProjection<FastJets>(event, "Jets").jetsByPt(1.0*GeV);
-      if (jets.size() < 1) vetoEvent;
-
-      FourMomentum p_lead = jets[0].momentum();
+      // Find the lead jet, applying a restriction that the jets must be within |eta| < 2.
+      FourMomentum p_lead;
+      foreach (const Jet& j, applyProjection<FastJets>(event, "Jets").jetsByPt(1.0*GeV)) {
+        if (fabs(j.momentum().eta()) < 2.0) {
+          p_lead = j.momentum();
+          break;
+        }
+      }
+      if (p_lead.isZero()) vetoEvent;
       const double philead = p_lead.phi();
       const double pTlead  = p_lead.pT();