[Rivet-svn] r1659 - in trunk: include/Rivet include/Rivet/Analyses src/Analyses

[blackhole at projects.hepforge.org]

Author: fsiegert
Date: Mon Jul  6 23:47:38 2009
New Revision: 1659

Log:
Add preliminary version of CDF_2008_S8093652. No data on HepData yet =>
not tested.

Added:
   trunk/include/Rivet/Analyses/CDF_2008_S8093652.hh
   trunk/src/Analyses/CDF_2008_S8093652.cc
Modified:
   trunk/include/Rivet/Makefile.am
   trunk/src/Analyses/Makefile.am
   trunk/src/Analyses/StdAnalyses.cc

Added: trunk/include/Rivet/Analyses/CDF_2008_S8093652.hh
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/include/Rivet/Analyses/CDF_2008_S8093652.hh	Mon Jul  6 23:47:38 2009	(r1659)
@@ -0,0 +1,103 @@
+// -*- C++ -*-
+#ifndef RIVET_CDF_2008_S8093652_HH
+#define RIVET_CDF_2008_S8093652_HH
+
+#include "Rivet/Analysis.hh"
+#include "Rivet/Tools/BinnedHistogram.hh"
+
+namespace Rivet {
+
+
+  class CDF_2008_S8093652 : public Analysis {
+
+  public:
+
+    /// @name Construction
+    //@{
+    /// Constructor
+    CDF_2008_S8093652();
+
+    /// Factory method 
+    static Analysis* create() {
+      return new CDF_2008_S8093652();
+    }
+    //@}
+
+
+    /// @name Publication metadata
+    //@{
+    /// A short description of the analysis. 
+    string spiresId() const {
+      return "8093652";
+    }
+    /// A short description of the analysis.
+    string summary() const {
+      return "Dijet mass spectrum";
+    }    
+    /// Experiment which performed and published this analysis. 
+    string experiment() const {
+      return "CDF";
+    }
+    /// Collider on which the experiment ran.
+    string collider() const {
+      return "Tevatron Run 2";
+    }
+    /// When published (preprint year according to SPIRES). 
+    string year() const {
+      return "2008";
+    }
+    /// Names & emails of analysis authors.
+    vector<string> authors() const {
+      vector<string> ret;
+      ret += "Frank Siegert <frank.siegert at durham.ac.uk>";
+      return ret;
+    }
+    /// A full description of the analysis.
+    string description() const {
+      ostringstream os;
+      os << "Dijet mass spectrum  "
+         << "from 0.2 TeV to 1.4 TeV "
+         << "in ppbar collisions at $\\sqrt{s} = 1.96$ TeV, based on "
+         << "an integrated luminosity of $1.13 \\text{fb}^{-1}$.";
+      return os.str();
+    }
+    /// Information about the events needed as input for this analysis.
+    string runInfo() const {
+      ostringstream os;
+      os << "Tevatron Run 2 conditions:\n" << endl << endl
+         << "* ppbar -> jets at 1960 GeV";
+      return os.str();
+    }
+    string status() const {
+      return "UNVALIDATED";
+    }
+    /// Publication references.
+    vector<string> references() const {
+      vector<string> ret;
+      ret.push_back("arXiv:0812.4036 [hep-ex]");
+      return ret;
+    }
+    //@}
+
+
+    /// @name Analysis methods
+    //@{ 
+    void init();
+    void analyze(const Event& event);
+    void finalize();
+    //@}
+
+
+  private:
+
+    /// @name Histograms
+    //@{
+    AIDA::IHistogram1D* _h_m_dijet;
+    //@}
+    
+  };
+
+
+}
+
+#endif

Modified: trunk/include/Rivet/Makefile.am
==============================================================================
--- trunk/include/Rivet/Makefile.am	Mon Jul  6 23:45:57 2009	(r1658)
+++ trunk/include/Rivet/Makefile.am	Mon Jul  6 23:47:38 2009	(r1659)
@@ -53,6 +53,7 @@
   Analyses/CDF_2008_LEADINGJETS.hh \
   Analyses/CDF_2008_NOTE_9351.hh \
   Analyses/CDF_2008_S7828950.hh \
+  Analyses/CDF_2008_S8093652.hh \
   Analyses/CDF_2008_S8095620.hh \
   Analyses/CDF_2009_S8233977.hh \
   Analyses/D0_2001_S4674421.hh  \

Added: trunk/src/Analyses/CDF_2008_S8093652.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/src/Analyses/CDF_2008_S8093652.cc	Mon Jul  6 23:47:38 2009	(r1659)
@@ -0,0 +1,66 @@
+// -*- C++ -*-
+#include "Rivet/Analyses/CDF_2008_S8093652.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/ChargedFinalState.hh"
+#include "Rivet/Projections/FastJets.hh"
+#include "Rivet/RivetAIDA.hh"
+
+namespace Rivet {
+
+
+  CDF_2008_S8093652::CDF_2008_S8093652()
+  {
+    setBeams(PROTON, ANTIPROTON);
+    setNeedsCrossSection(true);
+    
+    ChargedFinalState cfs;
+    FastJets conefinder(cfs, FastJets::CDFMIDPOINT, 0.7);
+    addProjection(conefinder, "ConeFinder");
+  } 
+
+
+  // Book histograms
+  void CDF_2008_S8093652::init() {
+    _h_m_dijet = bookHistogram1D(
+        1, 1, 1, "Dijet mass spectrum",
+        "$m_{\\text{jj}} [GeV/c^2]$",
+        "$d\\sigma/dm_{\\text{jj}} [pb/(GeV/c^2)]$");
+  }
+
+
+
+  // Do the analysis 
+  void CDF_2008_S8093652::analyze(const Event & e) {
+    double weight = e.weight();
+
+    const JetAlg& jetpro = applyProjection<JetAlg>(e, "ConeFinder");
+    const Jets& jets = jetpro.jetsByPt();
+    
+    if(jets.size()<2) {
+      vetoEvent;
+    }
+    
+    FourMomentum j0(jets[0].momentum());
+    FourMomentum j1(jets[1].momentum());
+    
+    if (fabs(j1.rapidity())>1.0 || fabs(j0.rapidity())>1.0) {
+      vetoEvent;
+    }
+    
+    double mjj = FourMomentum(j0+j1).mass();
+    
+    _h_m_dijet->fill(mjj, weight);
+  }
+
+
+
+  // Finalize
+  void CDF_2008_S8093652::finalize() {
+    /// Scale by L_eff = sig_MC * L_exp / num_MC
+    const double lumi_mc = sumOfWeights() / crossSection();
+    const double scalefactor =  1 / lumi_mc;
+    scale(_h_m_dijet, scalefactor);
+  }
+
+}

Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am	Mon Jul  6 23:45:57 2009	(r1658)
+++ trunk/src/Analyses/Makefile.am	Mon Jul  6 23:47:38 2009	(r1659)
@@ -20,6 +20,7 @@
     CDF_2008_S7540469.cc \
     CDF_2008_S7782535.cc \
     CDF_2008_S7828950.cc \
+    CDF_2008_S8093652.cc \
     CDF_2008_S8095620.cc \
     CDF_2009_S8233977.cc \
     D0_2001_S4674421.cc \

Modified: trunk/src/Analyses/StdAnalyses.cc
==============================================================================
--- trunk/src/Analyses/StdAnalyses.cc	Mon Jul  6 23:45:57 2009	(r1658)
+++ trunk/src/Analyses/StdAnalyses.cc	Mon Jul  6 23:47:38 2009	(r1659)
@@ -41,6 +41,7 @@
 #include "Rivet/Analyses/CDF_2008_LEADINGJETS.hh"
 #include "Rivet/Analyses/CDF_2008_S7782535.hh"
 #include "Rivet/Analyses/CDF_2008_S7828950.hh"
+#include "Rivet/Analyses/CDF_2008_S8093652.hh"
 #include "Rivet/Analyses/CDF_2008_S8095620.hh"
 #include "Rivet/Analyses/CDF_2009_S8233977.hh"
 #include "Rivet/Analyses/D0_2001_S4674421.hh"
@@ -116,6 +117,7 @@
     fns["CDF_2008_NOTE_9351"] = Rivet::CDF_2008_NOTE_9351::create;
     fns["CDF_2008_LEADINGJETS"] = Rivet::CDF_2008_LEADINGJETS::create;
     fns["CDF_2008_S7828950"] = Rivet::CDF_2008_S7828950::create;
+    fns["CDF_2008_S8093652"] = Rivet::CDF_2008_S8093652::create;
     fns["CDF_2008_S8095620"] = Rivet::CDF_2008_S8095620::create;
     fns["CDF_2009_S8233977"] = Rivet::CDF_2009_S8233977::create;
     fns["D0_2001_S4674421"] = Rivet::D0_2001_S4674421::create;