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

[blackhole at projects.hepforge.org]

Author: buckley
Date: Sat Jun 26 10:44:41 2010
New Revision: 2533

Log:
Fix D0 2004 analysis: the (mostly pointless) scalar ET cut should actually have been a vector missing ET

Modified:
   trunk/src/Analyses/D0_2004_S5992206.cc

Modified: trunk/src/Analyses/D0_2004_S5992206.cc
==============================================================================
--- trunk/src/Analyses/D0_2004_S5992206.cc	Fri Jun 25 20:40:26 2010	(r2532)
+++ trunk/src/Analyses/D0_2004_S5992206.cc	Sat Jun 26 10:44:41 2010	(r2533)
@@ -63,7 +63,7 @@
 
 
     /// Do the analysis
-    void analyze(const Event & event) {
+    void analyze(const Event& event) {
 
       // Analyse and print some info
       const JetAlg& jetpro = applyProjection<JetAlg>(event, "Jets");
@@ -84,8 +84,11 @@
       const double pT1 = jets[0].momentum().pT();
 
       const MissingMomentum& caloMissEt = applyProjection<MissingMomentum>(event, "CalMET");
-      getLog() << Log::DEBUG << "Missing scalar Et = " << caloMissEt.scalarET()/GeV << " GeV" << endl;
-      if (caloMissEt.scalarET() > 0.7*pT1) {
+      getLog() << Log::DEBUG << "Missing vector Et = " << caloMissEt.vectorET()/GeV << " GeV" << endl;
+      if (caloMissEt.vectorET() > 0.7*pT1) {
+        MSG_DEBUG("Vetoing event with too much missing Et: "
+                  << caloMissEt.vectorET()/GeV << " GeV > "
+                  << 0.7*pT1/GeV << " GeV");
         vetoEvent;
       }