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

[blackhole at projects.hepforge.org]

Author: buckley
Date: Wed Oct 12 16:18:13 2011
New Revision: 3437

Log:
Adding a jet multiplicity plot to MC_TTBAR, and making it a validated analysis

Modified:
   trunk/data/anainfo/MC_TTBAR.info
   trunk/src/Analyses/MC_TTBAR.cc
   trunk/src/Analyses/Makefile.am

Modified: trunk/data/anainfo/MC_TTBAR.info
==============================================================================
--- trunk/data/anainfo/MC_TTBAR.info	Wed Oct 12 16:17:15 2011	(r3436)
+++ trunk/data/anainfo/MC_TTBAR.info	Wed Oct 12 16:18:13 2011	(r3437)
@@ -1,6 +1,6 @@
 Name: MC_TTBAR
 Summary: MC analysis for ttbar studies
-Status: UNVALIDATED
+Status: VALIDATED
 Authors:
  - Hendrik Hoeth <hendrik.hoeth at cern.ch>
  - Andy Buckley <andy.buckley at cern.ch>

Modified: trunk/src/Analyses/MC_TTBAR.cc
==============================================================================
--- trunk/src/Analyses/MC_TTBAR.cc	Wed Oct 12 16:17:15 2011	(r3436)
+++ trunk/src/Analyses/MC_TTBAR.cc	Wed Oct 12 16:18:13 2011	(r3437)
@@ -39,6 +39,8 @@
       addProjection(MissingMomentum(fs), "MissingET");
 
       // Booking of histograms
+      _h_njets = bookHistogram1D("jet_mult", 11, -0.5, 10.5);
+      //
       _h_jet_1_pT = bookHistogram1D("jet_1_pT", 50, 0, 500);
       _h_jet_2_pT = bookHistogram1D("jet_2_pT", 50, 0, 400);
       _h_jet_3_pT = bookHistogram1D("jet_3_pT", 50, 0, 300);
@@ -103,6 +105,7 @@
       // Insist that the hardest 4 jets pass pT hardness cuts. If we don't find
       // at least 4 such jets, we abandon this event.
       const Jets jets = jetpro.jetsByPt(30*GeV);
+      _h_njets->fill(jets.size(), weight);
       double ht = 0.0;
       foreach (const Jet& j, jets) { ht += j.momentum().pT(); }
       _h_jet_HT->fill(ht/GeV, weight);
@@ -194,6 +197,7 @@
 
 
     void finalize() {
+      scale(_h_njets, 1/_sumwPassedLepJetMET);
       scale(_h_jet_1_pT, 1/_sumwPassedLepJetMET);
       scale(_h_jet_2_pT, 1/_sumwPassedLepJetMET);
       scale(_h_jet_3_pT, 1/_sumwPassedLepJetMET);
@@ -219,6 +223,7 @@
     // @name Histogram data members
     //@{
 
+    AIDA::IHistogram1D *_h_njets;
     AIDA::IHistogram1D *_h_jet_1_pT, *_h_jet_2_pT, *_h_jet_3_pT, *_h_jet_4_pT;
     AIDA::IHistogram1D *_h_jet_HT;
     AIDA::IHistogram1D *_h_bjet_1_pT, *_h_bjet_2_pT;

Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am	Wed Oct 12 16:17:15 2011	(r3436)
+++ trunk/src/Analyses/Makefile.am	Wed Oct 12 16:18:13 2011	(r3437)
@@ -256,11 +256,11 @@
     MC_ZZJETS.cc \
     MC_LEADINGJETS.cc \
     MC_DIPHOTON.cc \
-    MC_SUSY.cc
+    MC_SUSY.cc \
+    MC_TTBAR.cc \
+    MC_XS.cc
 if ENABLE_UNVALIDATED
 RivetMCAnalyses_la_SOURCES += \
     MC_DIJET.cc \
-    MC_PHOTONJETUE.cc \
-    MC_TTBAR.cc \
-    MC_XS.cc
+    MC_PHOTONJETUE.cc
 endif