[Rivet-svn] r1843 - in trunk: data/anainfo data/plotinfo src/Analyses

[blackhole at projects.hepforge.org]

Author: fsiegert
Date: Thu Oct  1 17:04:43 2009
New Revision: 1843

Log:
Add CDF_2009_S8383952.
Only reference data is still missing, waiting for Mike to put it into
HepData.

Added:
   trunk/data/anainfo/CDF_2009_S8383952.info
   trunk/data/plotinfo/CDF_2009_S8383952.plot
   trunk/src/Analyses/CDF_2009_S8383952.cc
Modified:
   trunk/data/anainfo/Makefile.am
   trunk/data/plotinfo/Makefile.am
   trunk/src/Analyses/Makefile.am

Added: trunk/data/anainfo/CDF_2009_S8383952.info
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/anainfo/CDF_2009_S8383952.info	Thu Oct  1 17:04:43 2009	(r1843)
@@ -0,0 +1,21 @@
+Name: CDF_2009_S8383952
+Year: 2009
+Summary: Z rapidity measurement
+Experiment: CDF
+Collider: Tevatron Run 2
+SpiresID: 8383952
+Status: UNFINISHED
+Authors:
+ - Frank Siegert <frank.siegert at durham.ac.uk>
+References:
+ - arXiv:0908.3914 [hep-ex]
+RunInfo:
+  "ppbar -> e+ e- + jets at 1960 GeV. Needs mass cut on lepton pair to 
+  avoid photon singularity: $m$ in range $66 < m_ee < 116$ GeV"
+NumEvents: 1000000
+PtCuts: [0]
+Description:
+  CDF measurement of the total cross section and rapidity distribution,
+  $d\sigma/dy$, for $q\bar{q}\to \gamma^{*}/Z\to e^{+}e^{-}$ events in the $Z$
+  boson mass region ($66<M_{ee}<116$ GeV/c$^2$) produced in $p\bar{p}$
+  collisions at $\sqrt{s}=1.96$ TeV with 2.1 fb$^{-1}$ of integrated luminosity.

Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am	Wed Sep 30 10:33:46 2009	(r1842)
+++ trunk/data/anainfo/Makefile.am	Thu Oct  1 17:04:43 2009	(r1843)
@@ -23,6 +23,7 @@
   CDF_2008_S8093652.info \
   CDF_2008_S8095620.info \
   CDF_2009_S8233977.info \
+  CDF_2009_S8383952.info \
   D0_1996_S3214044.info \
   D0_1996_S3324664.info \
   D0_2001_S4674421.info \

Added: trunk/data/plotinfo/CDF_2009_S8383952.plot
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/plotinfo/CDF_2009_S8383952.plot	Thu Oct  1 17:04:43 2009	(r1843)
@@ -0,0 +1,5 @@
+# BEGIN PLOT /CDF_2009_S8383952/d01-x01-y01
+Title=$e^+ e^-$ pair rapidity
+XLabel=$|y_{e^+ e^-}|$
+YLabel=$\mathrm{d}\sigma/\mathrm{d}y$ [pb]
+# END PLOT

Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am	Wed Sep 30 10:33:46 2009	(r1842)
+++ trunk/data/plotinfo/Makefile.am	Thu Oct  1 17:04:43 2009	(r1843)
@@ -23,6 +23,7 @@
   CDF_2008_S8093652.plot \
   CDF_2008_S8095620.plot \
   CDF_2009_S8233977.plot \
+  CDF_2009_S8383952.plot \
   D0_1996_S3214044.plot \
   D0_1996_S3324664.plot \
   D0_2001_S4674421.plot \

Added: trunk/src/Analyses/CDF_2009_S8383952.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/src/Analyses/CDF_2009_S8383952.cc	Thu Oct  1 17:04:43 2009	(r1843)
@@ -0,0 +1,88 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/ZFinder.hh"
+
+namespace Rivet {
+
+
+  class CDF_2009_S8383952 : public Analysis {
+  public:
+
+    /// @name Constructors etc.
+    //@{
+
+    /// Constructor
+    CDF_2009_S8383952()
+      : Analysis("CDF_2009_S8383952") 
+    {
+      setBeams(PROTON, ANTIPROTON);
+      setNeedsCrossSection(true);
+    }
+
+    //@}
+
+
+  public:
+
+    /// @name Analysis methods
+    //@{
+
+    /// Book histograms and initialise projections before the run
+    void init() {
+
+      /// Initialise and register projections here
+      ZFinder zfinder(-2.8, 2.8, 20.0*GeV, ELECTRON,
+                      66.0*GeV, 116.0*GeV, 0.2);
+      addProjection(zfinder, "ZFinder");
+
+
+      /// Book histograms here, e.g.:
+      _h_yZ = bookHistogram1D(1, 1, 1);
+
+    }
+
+
+    /// Perform the per-event analysis
+    void analyze(const Event& event) {
+      const double weight = event.weight();
+
+      const ZFinder& zfinder = applyProjection<ZFinder>(event, "ZFinder");
+      if (zfinder.particles().size() == 1) {
+        double yZ = fabs(zfinder.particles()[0].momentum().rapidity());
+        _h_yZ->fill(yZ, weight);
+      }
+      else {
+        getLog() << Log::DEBUG << "no unique lepton pair found." << endl;
+      }
+
+    }
+
+
+    /// Normalise histograms etc., after the run
+    void finalize() {
+      scale(_h_yZ, crossSection()/sumOfWeights());
+    }
+
+    //@}
+
+
+  private:
+
+    /// @name Histograms
+    //@{
+
+    AIDA::IHistogram1D *_h_yZ;
+    //@}
+
+  };
+
+
+
+  // This global object acts as a hook for the plugin system
+  AnalysisBuilder<CDF_2009_S8383952> plugin_CDF_2009_S8383952;
+
+
+}

Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am	Wed Sep 30 10:33:46 2009	(r1842)
+++ trunk/src/Analyses/Makefile.am	Thu Oct  1 17:04:43 2009	(r1843)
@@ -32,7 +32,8 @@
     CDF_2008_S7828950.cc \
     CDF_2008_S8093652.cc \
     CDF_2008_S8095620.cc \
-    CDF_2009_S8233977.cc
+    CDF_2009_S8233977.cc \
+    CDF_2009_S8383952.cc
 
 
 lib_LTLIBRARIES += RivetD0Analyses.la