[Rivet-svn] r3862 - schools/2012-CERN/setups/day2/rivet

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Wed Jul 25 09:31:51 BST 2012
Previous message: [Rivet-svn] r3861 - in schools/2012-CERN/setups/day2/herwig: ME QED shower
Next message: [Rivet-svn] r3863 - schools/2012-CERN/setups/day2/sherpa/QED
Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]
Author: hoeth
Date: Wed Jul 25 09:31:51 2012
New Revision: 3862

Log:
Rivet MC_ZJETS_NOCLUS plugin analysis for day2 tutorial

Added:
   schools/2012-CERN/setups/day2/rivet/
   schools/2012-CERN/setups/day2/rivet/MC_ZJETS_NOCLUS.cc
      - copied, changed from r3842, schools/2010-Lauterbad/setups/day2/rivet/MC_ZJETS_NOCLUS.cc
   schools/2012-CERN/setups/day2/rivet/MC_ZJETS_NOCLUS.info
   schools/2012-CERN/setups/day2/rivet/MC_ZJETS_NOCLUS.plot

Copied and modified: schools/2012-CERN/setups/day2/rivet/MC_ZJETS_NOCLUS.cc (from r3842, schools/2010-Lauterbad/setups/day2/rivet/MC_ZJETS_NOCLUS.cc)
==============================================================================
--- schools/2010-Lauterbad/setups/day2/rivet/MC_ZJETS_NOCLUS.cc	Wed Jul 18 14:40:46 2012	(r3842, copy source)
+++ schools/2012-CERN/setups/day2/rivet/MC_ZJETS_NOCLUS.cc	Wed Jul 25 09:31:51 2012	(r3862)
@@ -15,9 +15,7 @@
     /// Default constructor
     MC_ZJETS_NOCLUS()
       : MC_JetAnalysis("MC_ZJETS_NOCLUS", 4, "Jets")
-    {
-      setNeedsCrossSection(true);
-    }
+    {    }
 
 
     /// @name Analysis methods
@@ -25,7 +23,8 @@
 
     /// Book histograms
     void init() {
-      ZFinder zfinder(-3.5, 3.5, 25.0*GeV, ELECTRON, 65.0*GeV, 115.0*GeV, 0.2, false, true);
+      FinalState fs;
+      ZFinder zfinder(fs, -3.5, 3.5, 25.0*GeV, ELECTRON, 65.0*GeV, 115.0*GeV, 0.2, false, true);
       addProjection(zfinder, "ZFinder");
       FastJets jetpro(zfinder.remainingFinalState(), FastJets::KT, 0.7);
       addProjection(jetpro, "Jets");

Added: schools/2012-CERN/setups/day2/rivet/MC_ZJETS_NOCLUS.info
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2012-CERN/setups/day2/rivet/MC_ZJETS_NOCLUS.info	Wed Jul 25 09:31:51 2012	(r3862)
@@ -0,0 +1,17 @@
+Name: MC_ZJETS_NOCLUS
+Summary: Monte Carlo validation observables for $Z[e^+ \, e^-]$ + jets production
+Status: VALIDATED
+Authors:
+ - Frank Siegert <frank.siegert at cern.ch>
+References:
+RunInfo:
+  $e^+ e^-$ + jets analysis. Needs mass cut on lepton pair to avoid
+  photon singularity, e.g. a min range of $66 < m_{ee} < 116$ GeV
+NumEvents: 1000000
+PtCuts: [0]
+Description:
+  Available observables are Z mass, pT of jets 1-4, jet multiplicity,
+  $\Delta\eta(Z, \text{jet1})$, $\Delta R(\text{jet2}, \text{jet3})$,
+  differential jet rates 0->1, 1->2, 2->3, 3->4, integrated 0--4 jet
+  rates.
+NeedCrossSection: yes

Added: schools/2012-CERN/setups/day2/rivet/MC_ZJETS_NOCLUS.plot
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2012-CERN/setups/day2/rivet/MC_ZJETS_NOCLUS.plot	Wed Jul 25 09:31:51 2012	(r3862)
@@ -0,0 +1,402 @@
+# BEGIN PLOT /MC_ZJETS_NOCLUS/Z_jet1_dR
+Title=Separation between Z boson and leading jet
+XLabel=$\Delta R(\text{Z, 1st jet})$
+YLabel=$\text{d}\sigma/\text{d}\Delta R(\text{Z, 1st jet})$ [pb]
+LegendXPos=0.35
+LegendYPos=0.5
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/Z_jet1_deta
+Title=
+XLabel=$\Delta{\eta}(\text{Z, 1st jet})$
+YLabel=$\text{d}\sigma/\text{d}\Delta{\eta}(\text{Z, 1st jet})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/Z_mass
+Title=Z mass
+XLabel=$m_{\text{Z}}$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}m_{\text{Z}}$ [pb/GeV]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/Z_pT$
+Title=Z $p_\perp$
+XLabel=$p_\perp^{\text{Z}}$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{Z}}$ [pb/GeV]
+LogX=1
+LegendXPos=0.05
+LegendYPos=0.5
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/Z_pT_peak
+Title=Z $p_\perp$ in peak region
+XLabel=$p_\perp^{\text{Z}}$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{Z}}$ [pb/GeV]
+LogY=0
+LegendYPos=0.5
+LegendXPos=0.30
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/Z_y
+Title=Z rapidity
+XLabel=$y_{\text{Z}}$
+YLabel=$\text{d}\sigma/\text{d}y_{\text{Z}}$ [pb]
+LegendXPos=0.35
+LegendYPos=0.5
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/Z_phi
+Title=Z azimuthal angle
+XLabel=$\phi_{\text{Z}}$
+YLabel=$\text{d}\sigma/\text{d}\phi_{\text{Z}}$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/lepton_pT
+Title=Lepton $p_\perp$
+XLabel=$p_\perp^{\text{l}}$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{l}}$ [pb/GeV]
+LogX=1
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/lepton_eta
+Title=Lepton $\eta$
+XLabel=$\eta_{\text{l}}$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}\eta_{\text{l}}$ [pb/GeV]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jets_dR_
+Title=$\Delta R$ separation between jets
+LegendXPos=0.10
+LegendYPos=0.5
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jets_dR_12
+XLabel=$\Delta{R}(\text{jet 1, jet 2})$
+YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 1, jet 2})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jets_dR_13
+XLabel=$\Delta{R}(\text{jet 1, jet 3})$
+YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 1, jet 3})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jets_dR_23
+XLabel=$\Delta{R}(\text{jet 2, jet 3})$
+YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 2, jet 3})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jets_deta_
+Title=Pseudorapidity separation between jets
+LegendYPos=0.5
+LegendXPos=0.30
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jets_deta_12
+XLabel=$\Delta\eta(\text{jet 1, jet 2})$
+YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 1, jet 2})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jets_deta_13
+XLabel=$\Delta\eta(\text{jet 1, jet 3})$
+YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 1, jet 3})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jets_deta_23
+XLabel=$\Delta\eta(\text{jet 2, jet 3})$
+YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 2, jet 3})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jets_dphi_
+Title=Pseudorapidity separation between jets
+LegendXPos=0.1
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jets_dphi_12
+XLabel=$\Delta\phi(\text{jet 1, jet 2})$
+YLabel=$\text{d}\sigma/\text{d}\Delta\phi(\text{jet 1, jet 2})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jets_dphi_13
+XLabel=$\Delta\phi(\text{jet 1, jet 3})$
+YLabel=$\text{d}\sigma/\text{d}\Delta\phi(\text{jet 1, jet 3})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jets_dphi_23
+XLabel=$\Delta\phi(\text{jet 2, jet 3})$
+YLabel=$\text{d}\sigma/\text{d}\Delta\phi(\text{jet 2, jet 3})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_eta_1
+Title=Pseudorapidity of leading jet
+XLabel=$\eta(\text{jet 1})$
+YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 1})$ [pb]
+LegendYPos=0.5
+LegendXPos=0.30
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_eta_2
+Title=Pseudorapidity of second jet
+XLabel=$\eta(\text{jet 2})$
+YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 2})$ [pb]
+LegendYPos=0.5
+LegendXPos=0.30
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_eta_3
+Title=Pseudorapidity of third jet
+XLabel=$\eta(\text{jet 3})$
+YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 3})$ [pb]
+LegendYPos=0.5
+LegendXPos=0.30
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_eta_4
+Title=Pseudorapidity of fourth jet
+XLabel=$\eta(\text{jet 4})$
+YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 4})$ [pb]
+LegendYPos=0.5
+LegendXPos=0.30
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_eta_pmratio_1
+Title=Pseudorapidity $+/-$ ratio of first jet
+XLabel=$\eta(\text{jet 1})_+/\eta(\text{jet 1})_-$
+YLabel=$|\eta(\text{jet 1}|$
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_eta_pmratio_2
+Title=Pseudorapidity $+/-$ ratio of second jet
+XLabel=$\eta(\text{jet 2})_+/\eta(\text{jet 2})_-$
+YLabel=$|\eta(\text{jet 2}|$
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_eta_pmratio_3
+Title=Pseudorapidity $+/-$ ratio of third jet
+XLabel=$\eta(\text{jet 3})_+/\eta(\text{jet 3})_-$
+YLabel=$|\eta(\text{jet 3}|$
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_eta_pmratio_4
+Title=Pseudorapidity $+/-$ ratio of fourth jet
+XLabel=$\eta(\text{jet 4})_+/\eta(\text{jet 4})_-$
+YLabel=$|\eta(\text{jet 4}|$
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_y_1
+Title=Rapidity of first jet
+XLabel=$y(\text{jet 1})$
+YLabel=$\text{d}\sigma/\text{d}y(\text{jet 1})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_y_2
+Title=Rapidity of second jet
+XLabel=$y(\text{jet 2})$
+YLabel=$\text{d}\sigma/\text{d}y(\text{jet 2})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_y_3
+Title=Rapidity of third jet
+XLabel=$y(\text{jet 3})$
+YLabel=$\text{d}\sigma/\text{d}y(\text{jet 3})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_y_4
+Title=Rapidity of fourth jet
+XLabel=$y(\text{jet 4})$
+YLabel=$\text{d}\sigma/\text{d}y(\text{jet 4})$ [pb]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_y_pmratio_1
+Title=Rapidity $+/-$ ratio of first jet
+XLabel=$y(\text{jet 1})_+/y(\text{jet 1})_-$
+YLabel=$|y(\text{jet 1}|$
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_y_pmratio_2
+Title=Rapidity $+/-$ ratio of second jet
+XLabel=$y(\text{jet 2})_+/y(\text{jet 2})_-$
+YLabel=$|y(\text{jet 2}|$
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_y_pmratio_3
+Title=Rapidity $+/-$ ratio of third jet
+XLabel=$y(\text{jet 3})_+/y(\text{jet 3})_-$
+YLabel=$|y(\text{jet 3}|$
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_y_pmratio_4
+Title=Rapidity $+/-$ ratio of fourth jet
+XLabel=$y(\text{jet 4})_+/y(\text{jet 4})_-$
+YLabel=$|y(\text{jet 4}|$
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_mass_1
+Title=Mass of first jet
+XLabel=$m(\text{jet 1})$
+YLabel=$\text{d}\sigma/\text{d}m(\text{jet 1})$ [pb/GeV]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_mass_2
+Title=Mass of second jet
+XLabel=$m(\text{jet 2})$
+YLabel=$\text{d}\sigma/\text{d}m(\text{jet 2})$ [pb/GeV]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_mass_3
+Title=Mass of third jet
+XLabel=$m(\text{jet 3})$
+YLabel=$\text{d}\sigma/\text{d}m(\text{jet 3})$ [pb/GeV]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_mass_4
+Title=Mass of fourth jet
+XLabel=$m(\text{jet 4})$
+YLabel=$\text{d}\sigma/\text{d}m(\text{jet 4})$ [pb/GeV]
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_multi_exclusive
+Title=Exclusive jet multiplicity
+XLabel=$N_{\text{jet}}$
+YLabel=$\sigma(N_{\text{jet}})$ [pb]
+XMajorTickMarks=10
+XMinorTickMarks=0
+ErrorBands=1
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_multi_inclusive
+Title=Inclusive jet multiplicity
+XLabel=$N_{\text{jet}}$
+YLabel=$\sigma(\geq N_{\text{jet}})$ [pb]
+XMajorTickMarks=10
+XMinorTickMarks=0
+ErrorBands=1
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/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_ZJETS_NOCLUS/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=0.75
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/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_ZJETS_NOCLUS/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_ZJETS_NOCLUS/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_ZJETS_NOCLUS/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_ZJETS_NOCLUS/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}/\text{GeV})$ [pb]
+LegendXPos=0.15
+LegendYPos=0.5
+Rebin=2
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/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}/\text{GeV})$ [pb]
+LegendXPos=0.15
+LegendYPos=0.5
+Rebin=2
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/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}/\text{GeV})$ [pb]
+LegendXPos=0.15
+LegendYPos=0.5
+Rebin=2
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/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}/\text{GeV})$ [pb]
+LegendXPos=0.15
+LegendYPos=0.5
+Rebin=2
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_pT_1
+Title=Transverse momentum of leading jet
+XLabel=$p_\perp(\text{jet 1})$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 1})$ [pb/GeV]
+LogX=1
+LegendXPos=0.05
+LegendYPos=0.5
+XMin=20.0
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_pT_2
+Title=Transverse momentum of second jet
+XLabel=$p_\perp(\text{jet 2})$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 2})$ [pb/GeV]
+LogX=1
+LegendXPos=0.05
+LegendYPos=0.5
+XMin=20.0
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_pT_3
+Title=Transverse momentum of third jet
+XLabel=$p_\perp(\text{jet 3})$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 3})$ [pb/GeV]
+LogX=1
+LegendXPos=0.05
+LegendYPos=0.5
+XMin=20.0
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/jet_pT_4
+Title=Transverse momentum of fourth jet
+XLabel=$p_\perp(\text{jet 4})$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 4})$ [pb/GeV]
+LogX=1
+LegendXPos=0.05
+LegendYPos=0.5
+XMin=20.0
+# END PLOT
+
+# BEGIN PLOT /MC_ZJETS_NOCLUS/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
Previous message: [Rivet-svn] r3861 - in schools/2012-CERN/setups/day2/herwig: ME QED shower
Next message: [Rivet-svn] r3863 - schools/2012-CERN/setups/day2/sherpa/QED
Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]
More information about the Rivet-svn mailing list