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

[blackhole at projects.hepforge.org]

Author: holsch
Date: Sun Jan 17 18:55:19 2010
New Revision: 2218

Log:
Add D0 measurement of WpT from 1998. Plots can be found here:http://users.hepforge.org/~holsch/Rivet/Validation/D0_1998_S3711838-check/

Added:
   trunk/src/Analyses/D0_1998_S3711838.cc
Modified:
   trunk/src/Analyses/Makefile.am

Added: trunk/src/Analyses/D0_1998_S3711838.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/src/Analyses/D0_1998_S3711838.cc	Sun Jan 17 18:55:19 2010	(r2218)
@@ -0,0 +1,112 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/WFinder.hh"
+#include "Rivet/Tools/ParticleIdUtils.hh"
+/// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder...
+
+namespace Rivet {
+
+
+  class D0_1998_S3711838 : public Analysis {
+  public:
+
+    /// @name Constructors etc.
+    //@{
+
+    /// Constructor
+    D0_1998_S3711838()
+      : Analysis("D0_1998_S3711838") 
+    {
+      /// @todo Set approriate for your analysis
+      setBeams(PROTON, ANTIPROTON);
+      
+      /// @todo Set whether your finalize method needs the generator cross section
+      setNeedsCrossSection(false);
+    }
+
+    //@}
+
+
+  public:
+
+    /// @name Analysis methods
+    //@{
+
+    /// Book histograms and initialise projections before the run
+    void init() {
+      WFinder wfe(-5, 5, 0.0*GeV, ELECTRON, 60.0*GeV, 100.0*GeV, 0.2);
+      addProjection(wfe, "WFe");
+
+      // Book histogram
+      _hist_wpt = bookHistogram1D(1, 1, 1);
+
+    }
+
+
+    /// Perform the per-event analysis
+    void analyze(const Event& event) {
+      const double weight = event.weight();
+      const WFinder& wf = applyProjection<WFinder>(event, "WFe");
+      if (wf.size() == 0) {
+        getLog() << Log::DEBUG << "No W candidates found: vetoing" << endl;
+        vetoEvent;
+      }
+
+      // Require the electron to have ET > 25 GeV  and |eta| < 1.1
+      FourMomentum p_e;
+      int chg_e = 0;
+
+      foreach (const Particle& l, wf.constituentsFinalState().particles()) {
+        const FourMomentum pl = l.momentum();
+        if (abs(l.pdgId()) == ELECTRON) {
+          chg_e = PID::threeCharge(l.pdgId());
+          p_e = pl;
+          const double eta_e = fabs(p_e.pseudorapidity());
+          if ( (pl.Et()/GeV < 25.0) || (eta_e > 1.1) ) {
+            getLog() << Log::DEBUG << l.pdgId() << " ET,pT,eta: " << pl.Et()/GeV << "," << pl.pT()/GeV << "," << eta_e << " fails electron cut" << endl;
+            vetoEvent;
+          }
+        }
+        // Require ETmiss to be > 25 GeV
+        else if (abs(l.pdgId()) == NU_E) {
+          FourMomentum p_nu = l.momentum();
+          if (p_nu.Et()/GeV < 25.0) {
+            getLog() << Log::DEBUG << l.pdgId() << " ET(miss): " << p_nu.Et() << "fails ETmiss cut" << endl;
+            vetoEvent;
+          }
+        }
+      }
+      assert(chg_e != 0);
+
+      FourMomentum pW = wf.particles()[0].momentum();
+      getLog() << Log::DEBUG << "Dilepton mass = " << pW.mass()/GeV << " GeV"  << endl;
+      getLog() << Log::DEBUG << "Dilepton pT   = " << pW.pT()/GeV << " GeV" << endl;
+      _hist_wpt->fill(pW.pT()/GeV, weight);
+
+    }
+
+
+    /// Normalise histograms etc., after the run
+    void finalize() {
+      
+    }
+
+    //@}
+
+  private:
+
+    /// @name Histograms
+    AIDA::IHistogram1D *_hist_wpt;
+
+  };
+
+
+
+  // This global object acts as a hook for the plugin system
+  AnalysisBuilder<D0_1998_S3711838> plugin_D0_1998_S3711838;
+
+
+}

Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am	Sun Jan 17 18:27:07 2010	(r2217)
+++ trunk/src/Analyses/Makefile.am	Sun Jan 17 18:55:19 2010	(r2218)
@@ -26,6 +26,7 @@
     CDF_2008_S8095620.cc \
     D0_1996_S3214044.cc \
     D0_1996_S3324664.cc \
+    D0_1998_S3711838.cc \
     DELPHI_2003_WUD_03_11.cc \
     E735_1998_S3905616.cc \
     H1_1995_S3167097.cc \