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

[blackhole at projects.hepforge.org]

Author: fsiegert
Date: Wed Nov  4 11:57:36 2009
New Revision: 1995

Log:
Add chargedtrackmulti >= 2 requirement to photon fragmentation analysis.
This is a standard hadronic Z selection cut at LEP.
It becomes necessary only in the stupid ;-) case where one includes diphoton
events "e e -> gamma gamma" in the Monte-Carlo

Modified:
   trunk/src/Analyses/ALEPH_1996_S3196992.cc

Modified: trunk/src/Analyses/ALEPH_1996_S3196992.cc
==============================================================================
--- trunk/src/Analyses/ALEPH_1996_S3196992.cc	Wed Nov  4 10:58:19 2009	(r1994)
+++ trunk/src/Analyses/ALEPH_1996_S3196992.cc	Wed Nov  4 11:57:36 2009	(r1995)
@@ -3,6 +3,7 @@
 #include "Rivet/RivetAIDA.hh"
 #include "Rivet/Tools/Logging.hh"
 #include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/ChargedFinalState.hh"
 #include "Rivet/Projections/IdentifiedFinalState.hh"
 #include "Rivet/Projections/FastJets.hh"
 #include "Rivet/Projections/Thrust.hh"
@@ -31,6 +32,7 @@
       ifs.acceptId(PHOTON);
       addProjection(ifs, "Photons");
       addProjection(Thrust(fs), "Thrust");
+      addProjection(ChargedFinalState(), "CFS");
 
       // Book histograms
       _h_z_2jet_001 = bookHistogram1D(1, 1, 1);
@@ -47,6 +49,10 @@
     /// Perform the per-event analysis
     void analyze(const Event& event) {
       const double weight = event.weight();
+
+      if (applyProjection<FinalState>(event, "CFS").particles().size()<2) {
+        vetoEvent;
+      }
       
       const ParticleVector allphotons = applyProjection<IdentifiedFinalState>(event, "Photons").particles();
       ParticleVector photons;