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

[blackhole at projects.hepforge.org]

Author: fsiegert
Date: Tue Nov  3 15:38:41 2009
New Revision: 1987

Log:
Add cut on thrust axis angle to ALEPH photon fragmentation analysis.

Modified:
   trunk/src/Analyses/ALEPH_1996_S3196992.cc

Modified: trunk/src/Analyses/ALEPH_1996_S3196992.cc
==============================================================================
--- trunk/src/Analyses/ALEPH_1996_S3196992.cc	Tue Nov  3 14:34:58 2009	(r1986)
+++ trunk/src/Analyses/ALEPH_1996_S3196992.cc	Tue Nov  3 15:38:41 2009	(r1987)
@@ -5,6 +5,7 @@
 #include "Rivet/Projections/FinalState.hh"
 #include "Rivet/Projections/IdentifiedFinalState.hh"
 #include "Rivet/Projections/FastJets.hh"
+#include "Rivet/Projections/Thrust.hh"
 
 namespace Rivet {
 
@@ -29,6 +30,7 @@
       IdentifiedFinalState ifs(-MAXRAPIDITY, +MAXRAPIDITY, 5.0*GeV);
       ifs.acceptId(PHOTON);
       addProjection(ifs, "Photons");
+      addProjection(Thrust(fs), "Thrust");
 
       // Book histograms
       _h_z_2jet_001 = bookHistogram1D(1, 1, 1);
@@ -56,6 +58,11 @@
       if (photons.size()<1) {
         vetoEvent;
       }
+      
+      const Thrust& thrust = applyProjection<Thrust>(event, "Thrust");
+      if (fabs(cos(thrust.thrustAxis().theta()))>0.95) {
+        vetoEvent;
+      }
 
       const FastJets& durjet = applyProjection<FastJets>(event, "DurhamJets");