[Rivet-svn] r3070 - in trunk: data/plotinfo include/Rivet/Analyses src/Analyses

[blackhole at projects.hepforge.org]

Author: fsiegert
Date: Fri Apr 29 16:21:04 2011
New Revision: 3070

Log:
Add HT(jets) as observable to MC_JetAnalysis.

Modified:
   trunk/data/plotinfo/MC_HJETS.plot
   trunk/data/plotinfo/MC_JETS.plot
   trunk/data/plotinfo/MC_PHOTONJETS.plot
   trunk/data/plotinfo/MC_WJETS.plot
   trunk/data/plotinfo/MC_WWJETS.plot
   trunk/data/plotinfo/MC_ZJETS.plot
   trunk/data/plotinfo/MC_ZZJETS.plot
   trunk/include/Rivet/Analyses/MC_JetAnalysis.hh
   trunk/src/Analyses/MC_JetAnalysis.cc

Modified: trunk/data/plotinfo/MC_HJETS.plot
==============================================================================
--- trunk/data/plotinfo/MC_HJETS.plot	Fri Apr 29 15:54:21 2011	(r3069)
+++ trunk/data/plotinfo/MC_HJETS.plot	Fri Apr 29 16:21:04 2011	(r3070)
@@ -415,3 +415,10 @@
 LegendYPos=0.5
 XMin=20.0
 # END PLOT
+
+# BEGIN PLOT /MC_HJETS/jet_HT
+Title=Scalar sum of jet transverse momenta (HT)
+XLabel=$HT$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}HT$ [pb/GeV]
+LogX=1
+# END PLOT

Modified: trunk/data/plotinfo/MC_JETS.plot
==============================================================================
--- trunk/data/plotinfo/MC_JETS.plot	Fri Apr 29 15:54:21 2011	(r3069)
+++ trunk/data/plotinfo/MC_JETS.plot	Fri Apr 29 16:21:04 2011	(r3070)
@@ -352,3 +352,10 @@
 LegendYPos=0.5
 XMin=20.0
 # END PLOT
+
+# BEGIN PLOT /MC_JETS/jet_HT
+Title=Scalar sum of jet transverse momenta (HT)
+XLabel=$HT$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}HT$ [pb/GeV]
+LogX=1
+# END PLOT

Modified: trunk/data/plotinfo/MC_PHOTONJETS.plot
==============================================================================
--- trunk/data/plotinfo/MC_PHOTONJETS.plot	Fri Apr 29 15:54:21 2011	(r3069)
+++ trunk/data/plotinfo/MC_PHOTONJETS.plot	Fri Apr 29 16:21:04 2011	(r3070)
@@ -383,3 +383,10 @@
 LegendYPos=0.5
 XMin=20.0
 # END PLOT
+
+# BEGIN PLOT /MC_PHOTONJETS/jet_HT
+Title=Scalar sum of jet transverse momenta (HT)
+XLabel=$HT$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}HT$ [pb/GeV]
+LogX=1
+# END PLOT

Modified: trunk/data/plotinfo/MC_WJETS.plot
==============================================================================
--- trunk/data/plotinfo/MC_WJETS.plot	Fri Apr 29 15:54:21 2011	(r3069)
+++ trunk/data/plotinfo/MC_WJETS.plot	Fri Apr 29 16:21:04 2011	(r3070)
@@ -455,3 +455,10 @@
 LegendYPos=0.5
 XMin=20.0
 # END PLOT
+
+# BEGIN PLOT /MC_WJETS/jet_HT
+Title=Scalar sum of jet transverse momenta (HT)
+XLabel=$HT$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}HT$ [pb/GeV]
+LogX=1
+# END PLOT

Modified: trunk/data/plotinfo/MC_WWJETS.plot
==============================================================================
--- trunk/data/plotinfo/MC_WWJETS.plot	Fri Apr 29 15:54:21 2011	(r3069)
+++ trunk/data/plotinfo/MC_WWJETS.plot	Fri Apr 29 16:21:04 2011	(r3070)
@@ -532,3 +532,10 @@
 LegendYPos=0.5
 XMin=20.0
 # END PLOT
+
+# BEGIN PLOT /MC_WWJETS/jet_HT
+Title=Scalar sum of jet transverse momenta (HT)
+XLabel=$HT$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}HT$ [pb/GeV]
+LogX=1
+# END PLOT

Modified: trunk/data/plotinfo/MC_ZJETS.plot
==============================================================================
--- trunk/data/plotinfo/MC_ZJETS.plot	Fri Apr 29 15:54:21 2011	(r3069)
+++ trunk/data/plotinfo/MC_ZJETS.plot	Fri Apr 29 16:21:04 2011	(r3070)
@@ -417,3 +417,10 @@
 LegendYPos=0.5
 XMin=20.0
 # END PLOT
+
+# BEGIN PLOT /MC_ZJETS/jet_HT
+Title=Scalar sum of jet transverse momenta (HT)
+XLabel=$HT$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}HT$ [pb/GeV]
+LogX=1
+# END PLOT

Modified: trunk/data/plotinfo/MC_ZZJETS.plot
==============================================================================
--- trunk/data/plotinfo/MC_ZZJETS.plot	Fri Apr 29 15:54:21 2011	(r3069)
+++ trunk/data/plotinfo/MC_ZZJETS.plot	Fri Apr 29 16:21:04 2011	(r3070)
@@ -517,3 +517,10 @@
 LegendYPos=0.5
 XMin=20.0
 # END PLOT
+
+# BEGIN PLOT /MC_ZZJETS/jet_HT
+Title=Scalar sum of jet transverse momenta (HT)
+XLabel=$HT$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}HT$ [pb/GeV]
+LogX=1
+# END PLOT

Modified: trunk/include/Rivet/Analyses/MC_JetAnalysis.hh
==============================================================================
--- trunk/include/Rivet/Analyses/MC_JetAnalysis.hh	Fri Apr 29 15:54:21 2011	(r3069)
+++ trunk/include/Rivet/Analyses/MC_JetAnalysis.hh	Fri Apr 29 16:21:04 2011	(r3070)
@@ -57,6 +57,7 @@
     AIDA::IHistogram1D * _h_jet_multi_exclusive;
     AIDA::IHistogram1D * _h_jet_multi_inclusive;
     AIDA::IDataPointSet * _h_jet_multi_ratio;
+    AIDA::IHistogram1D * _h_jet_HT;
     //@}
 
   };

Modified: trunk/src/Analyses/MC_JetAnalysis.cc
==============================================================================
--- trunk/src/Analyses/MC_JetAnalysis.cc	Fri Apr 29 15:54:21 2011	(r3069)
+++ trunk/src/Analyses/MC_JetAnalysis.cc	Fri Apr 29 16:21:04 2011	(r3070)
@@ -78,6 +78,7 @@
     _h_jet_multi_exclusive = bookHistogram1D("jet_multi_exclusive", m_njet+3, -0.5, m_njet+3-0.5);
     _h_jet_multi_inclusive = bookHistogram1D("jet_multi_inclusive", m_njet+3, -0.5, m_njet+3-0.5);
     _h_jet_multi_ratio = bookDataPointSet("jet_multi_ratio", m_njet+2, 0.5, m_njet+3-0.5);
+    _h_jet_HT = bookHistogram1D("jet_HT", logBinEdges(50, m_jetptcut, sqrtS()/GeV/2.0));
   }
 
 
@@ -173,6 +174,12 @@
         _h_jet_multi_inclusive->fill(i, weight);
       }
     }
+
+    double HT=0.0;
+    foreach (const Jet& jet, jets) {
+      HT += jet.momentum().pT();
+    }
+    _h_jet_HT->fill(HT, weight);
   }
 
 
@@ -229,6 +236,7 @@
 
     scale(_h_jet_multi_exclusive, crossSection()/sumOfWeights());
     scale(_h_jet_multi_inclusive, crossSection()/sumOfWeights());
+    scale(_h_jet_HT, crossSection()/sumOfWeights());
   }