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

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Wed Mar 3 17:04:54 GMT 2010
Previous message: [Rivet-svn] r2306 - trunk/bin
Next message: [Rivet-svn] r2308 - trunk/src/Analyses
Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]
Author: fsiegert
Date: Wed Mar  3 17:04:54 2010
New Revision: 2307

Log:
Add MC_JETS analysis for generic jet production MC validation.

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

Added: trunk/data/anainfo/MC_JETS.info
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/anainfo/MC_JETS.info	Wed Mar  3 17:04:54 2010	(r2307)
@@ -0,0 +1,13 @@
+Name: MC_JETS
+Summary: Monte Carlo validation observables for jet production
+Status: VALIDATED
+Authors:
+ - Frank Siegert <frank.siegert at durham.ac.uk>
+References:
+RunInfo:
+  Pure QCD jet production events at an arbitrary collider.
+NumEvents: 1000000
+PtCuts: [15]
+Description:
+  Jets with $p_\perp>20$ GeV are constructed with a $k_\perp$ jet finder with
+  $D=0.7$ and projected onto many different observables.

Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am	Wed Mar  3 16:11:51 2010	(r2306)
+++ trunk/data/anainfo/Makefile.am	Wed Mar  3 17:04:54 2010	(r2307)
@@ -69,6 +69,7 @@
   MC_PHOTONJETUE.info \
   MC_SUSY.info \
   MC_TTBAR.info \
+  MC_JETS.info \
   MC_PHOTONJETS.info \
   MC_WJETS.info \
   MC_ZJETS.info \

Added: trunk/data/plotinfo/MC_JETS.plot
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/plotinfo/MC_JETS.plot	Wed Mar  3 17:04:54 2010	(r2307)
@@ -0,0 +1,221 @@
+# BEGIN PLOT /MC_JETS/jets_dR_12
+Title=
+XLabel=$\Delta{R}(\text{jet 1, jet 2})$
+YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 1, jet 2})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jets_dR_13
+Title=
+XLabel=$\Delta{R}(\text{jet 1, jet 3})$
+YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 1, jet 3})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jets_dR_14
+Title=
+XLabel=$\Delta{R}(\text{jet 1, jet 4})$
+YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 1, jet 4})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jets_dR_23
+Title=
+XLabel=$\Delta{R}(\text{jet 2, jet 3})$
+YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 2, jet 3})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jets_dR_24
+Title=
+XLabel=$\Delta{R}(\text{jet 2, jet 4})$
+YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 2, jet 4})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jets_dR_34
+Title=
+XLabel=$\Delta{R}(\text{jet 3, jet 4})$
+YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 3, jet 4})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jets_deta_12
+Title=
+XLabel=$\Delta\eta(\text{jet 1, jet 2})$
+YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 1, jet 2})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jets_deta_13
+Title=
+XLabel=$\Delta\eta(\text{jet 1, jet 3})$
+YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 1, jet 3})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jets_deta_14
+Title=
+XLabel=$\Delta\eta(\text{jet 1, jet 4})$
+YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 1, jet 4})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jets_deta_23
+Title=
+XLabel=$\Delta\eta(\text{jet 2, jet 3})$
+YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 2, jet 3})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jets_deta_24
+Title=
+XLabel=$\Delta\eta(\text{jet 2, jet 4})$
+YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 2, jet 4})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jets_deta_34
+Title=
+XLabel=$\Delta\eta(\text{jet 3, jet 4})$
+YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 3, jet 4})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jet_eta_1
+Title=
+XLabel=$\eta(\text{jet 1})$
+YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 1})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jet_eta_2
+Title=
+XLabel=$\eta(\text{jet 2})$
+YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 2})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jet_eta_3
+Title=
+XLabel=$\eta(\text{jet 3})$
+YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 3})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jet_eta_4
+Title=
+XLabel=$\eta(\text{jet 4})$
+YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 4})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jet_multi_exclusive
+Title=Exclusive jet multiplicity
+XLabel=$N_{\text{jet}}$
+YLabel=$\sigma(N_{\text{jet}})$
+XMajorTickMarks=10
+XMinorTickMarks=0
+ErrorBands=1
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jet_multi_inclusive
+Title=Inclusive jet multiplicity
+XLabel=$N_{\text{jet}}$
+YLabel=$\sigma(\geq N_{\text{jet}})$
+XMajorTickMarks=10
+XMinorTickMarks=0
+ErrorBands=1
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jet_multi_ratio
+Title=Ratio of jet multiplicity
+XLabel=$N_{\text{jet}}$
+YLabel=$\sigma(\geq N_{\text{jet}})/\sigma(\geq N_{\text{jet}}-1)$
+XMajorTickMarks=10
+XMinorTickMarks=0
+LogY=0
+ErrorBands=1
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/log10_R_0
+Title=$\log_{10}$(Integrated $0$ jet rate in $k_\perp$ [GeV])
+XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
+YLabel=$R_{0}$
+Rebin=2
+LegendYPos=0.8
+LegendXPos=1.2
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/log10_R_1
+Title=$\log_{10}$(Integrated $1$ jet rate in $k_\perp$ [GeV])
+XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
+YLabel=$R_{1}$
+Rebin=2
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/log10_R_2
+Title=$\log_{10}$(Integrated $2$ jet rate in $k_\perp$ [GeV])
+XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
+YLabel=$R_{2}$
+Rebin=2
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/log10_R_3
+Title=$\log_{10}$(Integrated $3$ jet rate in $k_\perp$ [GeV])
+XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
+YLabel=$R_{3}$
+Rebin=2
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/log10_R_4
+Title=$\log_{10}$(Integrated $4$ jet rate in $k_\perp$ [GeV])
+XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
+YLabel=$R_{\geq4}$
+Rebin=2
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/log10_d_01
+Title=$\log_{10}$($k_\perp$ jet resolution $0 \to 1$ [GeV])
+XLabel=$\log_{10}(d_{01}/\text{GeV})$
+YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{01})$
+LegendXPos=0.5
+LegendYPos=0.5
+Rebin=2
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/log10_d_12
+Title=$\log_{10}$($k_\perp$ jet resolution $1 \to 2$ [GeV])
+XLabel=$\log_{10}(d_{12}/\text{GeV})$
+YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{12})$
+LegendXPos=0.5
+LegendYPos=0.5
+Rebin=2
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/log10_d_23
+Title=$\log_{10}$($k_\perp$ jet resolution $2 \to 3$ [GeV])
+XLabel=$\log_{10}(d_{23}/\text{GeV})$
+YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{23})$
+LegendXPos=0.5
+LegendYPos=0.5
+Rebin=2
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/log10_d_34
+Title=$\log_{10}$($k_\perp$ jet resolution $3 \to 4$ [GeV])
+XLabel=$\log_{10}(d_{34}/\text{GeV})$
+YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{34})$
+LegendXPos=0.5
+LegendYPos=0.5
+Rebin=2
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jet_pT_1
+Title=
+XLabel=$p_\perp(\text{jet 1})$
+YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 1})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jet_pT_2
+Title=
+XLabel=$p_\perp(\text{jet 2})$
+YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 2})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jet_pT_3
+Title=
+XLabel=$p_\perp(\text{jet 3})$
+YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 3})$
+# END PLOT
+
+# BEGIN PLOT /MC_JETS/jet_pT_4
+Title=
+XLabel=$p_\perp(\text{jet 4})$
+YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 4})$
+# END PLOT
+

Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am	Wed Mar  3 16:11:51 2010	(r2306)
+++ trunk/data/plotinfo/Makefile.am	Wed Mar  3 17:04:54 2010	(r2307)
@@ -67,6 +67,7 @@
   MC_PHOTONJETUE.plot \
   MC_SUSY.plot \
   MC_TTBAR.plot \
+  MC_JETS.plot \
   MC_PHOTONJETS.plot \
   MC_WJETS.plot \
   MC_ZJETS.plot \

Added: trunk/src/Analyses/MC_JETS.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/src/Analyses/MC_JETS.cc	Wed Mar  3 17:04:54 2010	(r2307)
@@ -0,0 +1,45 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Analyses/MC_JetAnalysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/FastJets.hh"
+
+namespace Rivet {
+
+
+  class MC_JETS : public MC_JetAnalysis {
+  public:
+
+    MC_JETS() : MC_JetAnalysis("MC_JETS", 4, "Jets") 
+    {
+      setNeedsCrossSection(true);
+    }
+
+
+  public:
+
+    void init() {
+      FinalState fs;
+      FastJets jetpro(fs, FastJets::KT, 0.7, 20.0*GeV);
+      addProjection(jetpro, "Jets");
+      
+      MC_JetAnalysis::init();
+    }
+
+
+    void analyze(const Event& event) {
+      MC_JetAnalysis::analyze(event);
+    }
+
+
+    void finalize() {
+      MC_JetAnalysis::finalize();
+    }
+
+  };
+
+  AnalysisBuilder<MC_JETS> plugin_MC_JETS;
+
+}

Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am	Wed Mar  3 16:11:51 2010	(r2306)
+++ trunk/src/Analyses/Makefile.am	Wed Mar  3 17:04:54 2010	(r2307)
@@ -127,6 +127,7 @@
 RivetMCAnalyses_la_SOURCES = \
     ExampleAnalysis.cc \
     MC_JetAnalysis.cc \
+    MC_JETS.cc \
     MC_PHOTONJETS.cc \
     MC_WJETS.cc \
     MC_ZJETS.cc \
Previous message: [Rivet-svn] r2306 - trunk/bin
Next message: [Rivet-svn] r2308 - trunk/src/Analyses
Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]
More information about the Rivet-svn mailing list