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

[blackhole at projects.hepforge.org]

Author: holsch
Date: Thu Oct 15 19:04:42 2009
New Revision: 1917

Log:
Remove constraint for primary vertex since it made Jimmy crash and should not be an issue in generators.

Modified:
   trunk/src/Analyses/D0_2004_S5992206.cc

Modified: trunk/src/Analyses/D0_2004_S5992206.cc
==============================================================================
--- trunk/src/Analyses/D0_2004_S5992206.cc	Thu Oct 15 18:37:23 2009	(r1916)
+++ trunk/src/Analyses/D0_2004_S5992206.cc	Thu Oct 15 19:04:42 2009	(r1917)
@@ -3,7 +3,6 @@
 #include "Rivet/RivetAIDA.hh"
 #include "Rivet/Tools/Logging.hh"
 #include "Rivet/Projections/FastJets.hh"
-#include "Rivet/Projections/PVertex.hh"
 #include "Rivet/Projections/TotalVisibleMomentum.hh"
 
 namespace Rivet {
@@ -48,7 +47,6 @@
       addProjection(fs, "FS");
       addProjection(FastJets(fs, FastJets::D0ILCONE, 0.7, 6*GeV), "Jets");
       addProjection(TotalVisibleMomentum(fs), "CalMET");
-      addProjection(PVertex(), "PV");
       
       // Veto neutrinos, and muons with pT above 1.0 GeV
       VetoedFinalState vfs(fs);
@@ -71,10 +69,6 @@
       const JetAlg& jetpro = applyProjection<JetAlg>(event, "Jets");
       getLog() << Log::DEBUG << "Jet multiplicity before any pT cut = " << jetpro.size() << endl;
       
-      // Find vertex and check  that its z-component is < 50 cm from the nominal IP
-      const PVertex& pv = applyProjection<PVertex>(event, "PV");
-      if (fabs(pv.position().z())/cm > 50.0) vetoEvent;
-      
       const Jets jets  = jetpro.jetsByPt(40.0*GeV);
       if (jets.size() >= 2) {
         getLog() << Log::DEBUG << "Jet multiplicity after pT > 40 GeV cut = " << jets.size() << endl;