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

[blackhole at projects.hepforge.org]

Author: holsch
Date: Tue Oct 27 22:32:43 2009
New Revision: 1965

Log:
Adding first bits of MC analysis for ttbar studies, this commit includes the Pseudorapidity distribution only.

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

Added: trunk/data/anainfo/MC_LHC_TTBAR.info
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/anainfo/MC_LHC_TTBAR.info	Tue Oct 27 22:32:43 2009	(r1965)
@@ -0,0 +1,17 @@
+Name: MC_LHC_TTBAR
+Year: 
+Summary: MC analysis for ttbar studies at the LHC
+Experiment: 
+Collider: LHC
+SpiresID: 
+Status: UNVALIDATED
+Authors:
+ - Holger Schulz <hschulz at physik.hu-berlin.de>
+References:
+RunInfo:
+  For Pythia6, set MSEL=6
+  For fortran Herwig/Jimmy select IPROC=1706
+NumEvents: 10000
+PtCuts: 
+Description:
+  This is a pure Monte Carlo study for t-tbar production at the LHC.

Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am	Tue Oct 27 14:18:14 2009	(r1964)
+++ trunk/data/anainfo/Makefile.am	Tue Oct 27 22:32:43 2009	(r1965)
@@ -56,6 +56,7 @@
   MC_LHC_DIJET.info \
   MC_LHC_PHOTONJETUE.info \
   MC_LHC_SUSY.info \
+  MC_LHC_TTBAR.info \
   MC_LHC_WANALYSIS.info \
   MC_LHC_ZANALYSIS.info \
   MC_TVT1960_PHOTONJETS.info \

Added: trunk/data/plotinfo/MC_LHC_TTBAR.plot
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/plotinfo/MC_LHC_TTBAR.plot	Tue Oct 27 22:32:43 2009	(r1965)
@@ -0,0 +1,8 @@
+# BEGIN PLOT /MC_LHC_TTBAR/d01-x01-y01
+Title=[Insert title for histogram d01-x01-y01 here]
+XLabel=[Insert x-axis label for histogram d01-x01-y01 here]
+YLabel=[Insert y-axis label for histogram d01-x01-y01 here]
+# + any additional plot settings you might like, see make-plots documentation
+# END PLOT
+
+# ... add more histograms as you need them ...

Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am	Tue Oct 27 14:18:14 2009	(r1964)
+++ trunk/data/plotinfo/Makefile.am	Tue Oct 27 22:32:43 2009	(r1965)
@@ -56,6 +56,7 @@
   MC_LHC_ZANALYSIS.plot \
   MC_LHC_WANALYSIS.plot \
   MC_LHC_SUSY.plot \
+  MC_LHC_TTBAR.plot \
   MC_TVT1960_PHOTONJETS.plot \
   MC_TVT1960_ZJETS.plot \
   OPAL_1998_S3780481.plot \

Added: trunk/src/Analyses/MC_LHC_TTBAR.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/src/Analyses/MC_LHC_TTBAR.cc	Tue Oct 27 22:32:43 2009	(r1965)
@@ -0,0 +1,85 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Projections/ChargedFinalState.hh"
+/// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder...
+
+namespace Rivet {
+
+
+  class MC_LHC_TTBAR : public Analysis {
+  public:
+
+    /// @name Constructors etc.
+    //@{
+
+    /// Constructor
+    MC_LHC_TTBAR()
+      : Analysis("MC_LHC_TTBAR") 
+    {
+      /// @todo Set approriate for your analysis
+      setBeams(PROTON, PROTON);
+      
+      /// @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() {
+
+      const ChargedFinalState cfs(-5.0, 5.0);
+      addProjection(cfs, "CFS");
+
+      /// @todo Book histograms here, e.g.:
+       _histPseudorapidity = bookHistogram1D("pseudorap", 20, -5.0, 5.0);
+
+    }
+
+
+    /// Perform the per-event analysis
+    void analyze(const Event& event) {
+      const double weight = event.weight();
+      const ChargedFinalState& cfs = applyProjection<ChargedFinalState>(event, "CFS");
+      
+      foreach (const Particle& p, cfs.particles()) {
+        double eta = p.momentum().pseudorapidity();
+        _histPseudorapidity->fill(eta, weight);
+     }
+
+    }
+
+
+    /// Normalise histograms etc., after the run
+    void finalize() {
+    scale(_histPseudorapidity, 1.0/sumOfWeights());
+      
+    }
+
+
+  private:
+
+    /// @name Histograms
+    //@{
+
+    AIDA::IHistogram1D *_histPseudorapidity;
+    //@}
+
+  };
+
+
+
+  // This global object acts as a hook for the plugin system
+  AnalysisBuilder<MC_LHC_TTBAR> plugin_MC_LHC_TTBAR;
+
+
+}

Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am	Tue Oct 27 14:18:14 2009	(r1964)
+++ trunk/src/Analyses/Makefile.am	Tue Oct 27 22:32:43 2009	(r1965)
@@ -112,4 +112,5 @@
     MC_LHC_PHOTONJETUE.cc \
     MC_LHC_WANALYSIS.cc \
     MC_LHC_ZANALYSIS.cc \
-    MC_LHC_SUSY.cc
+    MC_LHC_SUSY.cc \
+    MC_LHC_TTBAR.cc