|
|
[Rivet-svn] r3593 - in trunk: data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Thu Mar 1 09:59:22 GMT 2012
Author: buckley
Date: Thu Mar 1 09:59:22 2012
New Revision: 3593
Log:
A couple of DOS to Unix newline consistency fixes
Modified:
trunk/data/plotinfo/LHCB_2011_I917009.plot
trunk/src/Analyses/CMS_2011_S9215166.cc
Modified: trunk/data/plotinfo/LHCB_2011_I917009.plot
==============================================================================
--- trunk/data/plotinfo/LHCB_2011_I917009.plot Wed Feb 29 19:50:30 2012 (r3592)
+++ trunk/data/plotinfo/LHCB_2011_I917009.plot Thu Mar 1 09:59:22 2012 (r3593)
@@ -1,236 +1,236 @@
-# BEGIN PLOT /LHCB_2011_I917009/d01-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 0.65$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-ErrorBars=1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d01-x01-y02
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.65 < p_\perp < 1.00$ Gev/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-LegendXPos=0.05
-LegendYPos=0.4
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d01-x01-y03
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-LegendXPos=0.05
-LegendYPos=0.35
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d02-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 0.65$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendXPos=0.05
-LegendYPos=0.35
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d02-x01-y02
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.65 < p_\perp < 1.00$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendXPos=0.05
-LegendYPos=0.33
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d02-x01-y03
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-YMax=0.7
-FullRange=1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d03-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ Gev/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-LegendXPos=0.05
-LegendYPos=0.55
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d04-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendXPos=0.05
-LegendYPos=0.35
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d05-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($2.0 < y < 4.0$)
-XLabel=$p_\perp$ [GeV/$c$]
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-LegendYPos=0.48
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d06-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($2.0 < y < 4.0$)
-XLabel=$p_\perp$ [GeV/$c$]
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendYPos=0.3
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d07-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity loss ($\Delta y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-LegendYPos=0.5
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d08-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ Gev/$c$)
-XLabel=Rapidity loss ($\Delta y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendYPos=0.3
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d09-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 0.65$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-YMin=0.65
-YMax=1.2
-FullRange=1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d09-x01-y02
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.65 < p_\perp < 1.00$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-YMin=0.65
-LegendYPos=0.3
-LegendXPos=0.1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d09-x01-y03
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-YMin=0.65
-LegendYPos=0.3
-LegendXPos=0.1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d10-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 0.65$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-YMin=0.1
-LegendXPos=0.4
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d10-x01-y02
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.65 < p_\perp < 1.00$ GeV/$c$)
-XLabel=Rapidity (y)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-YMin=0.1
-YMax=0.48
-FullRange=1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d10-x01-y03
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-YMin=0.1
-YMax=0.7
-FullRange=1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d11-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-YMin=0.65
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d12-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.5$ GeV/$c$)
-XLabel=Rapidity (y)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-YMin=0.1
-YMax=0.43
-FullRange=1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d13-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($2.0 < y < 4.5$)
-XLabel=$p_\perp$ [GeV/$c$]
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-YMin=0.65
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d14-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($2.0 < y < 4.5$)
-XLabel=$p_\perp$ [GeV/$c$]
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendYPos=0.33
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d15-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity loss ($\Delta y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-YMin=0.25
-LegendYPos=0.65
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d16-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity loss ($\Delta y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendYPos=0.33
-RatioPlot=0
-# END PLOT
+# BEGIN PLOT /LHCB_2011_I917009/d01-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 0.65$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+ErrorBars=1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d01-x01-y02
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.65 < p_\perp < 1.00$ Gev/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.4
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d01-x01-y03
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.35
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d02-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 0.65$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.35
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d02-x01-y02
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.65 < p_\perp < 1.00$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.33
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d02-x01-y03
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+YMax=0.7
+FullRange=1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d03-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ Gev/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.55
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d04-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.35
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d05-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($2.0 < y < 4.0$)
+XLabel=$p_\perp$ [GeV/$c$]
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+LegendYPos=0.48
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d06-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($2.0 < y < 4.0$)
+XLabel=$p_\perp$ [GeV/$c$]
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendYPos=0.3
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d07-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity loss ($\Delta y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+LegendYPos=0.5
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d08-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ Gev/$c$)
+XLabel=Rapidity loss ($\Delta y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendYPos=0.3
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d09-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 0.65$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+YMin=0.65
+YMax=1.2
+FullRange=1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d09-x01-y02
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.65 < p_\perp < 1.00$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+YMin=0.65
+LegendYPos=0.3
+LegendXPos=0.1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d09-x01-y03
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+YMin=0.65
+LegendYPos=0.3
+LegendXPos=0.1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d10-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 0.65$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+YMin=0.1
+LegendXPos=0.4
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d10-x01-y02
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.65 < p_\perp < 1.00$ GeV/$c$)
+XLabel=Rapidity (y)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+YMin=0.1
+YMax=0.48
+FullRange=1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d10-x01-y03
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+YMin=0.1
+YMax=0.7
+FullRange=1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d11-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+YMin=0.65
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d12-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.5$ GeV/$c$)
+XLabel=Rapidity (y)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+YMin=0.1
+YMax=0.43
+FullRange=1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d13-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($2.0 < y < 4.5$)
+XLabel=$p_\perp$ [GeV/$c$]
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+YMin=0.65
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d14-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($2.0 < y < 4.5$)
+XLabel=$p_\perp$ [GeV/$c$]
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendYPos=0.33
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d15-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity loss ($\Delta y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+YMin=0.25
+LegendYPos=0.65
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d16-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity loss ($\Delta y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendYPos=0.33
+RatioPlot=0
+# END PLOT
Modified: trunk/src/Analyses/CMS_2011_S9215166.cc
==============================================================================
--- trunk/src/Analyses/CMS_2011_S9215166.cc Wed Feb 29 19:50:30 2012 (r3592)
+++ trunk/src/Analyses/CMS_2011_S9215166.cc Thu Mar 1 09:59:22 2012 (r3593)
@@ -1,113 +1,113 @@
-// -*- C++ -*-
-#include "Rivet/Analysis.hh"
-#include "Rivet/RivetAIDA.hh"
-#include "Rivet/Tools/Logging.hh"
-#include "Rivet/Projections/FinalState.hh"
-#include "Rivet/Projections/ChargedFinalState.hh"
-#include "Rivet/Projections/FastJets.hh"
-#include "Rivet/Projections/VetoedFinalState.hh"
-
-namespace Rivet {
-
-
- class CMS_2011_S9215166 : public Analysis {
- public:
-
- /// Constructor
- CMS_2011_S9215166() : Analysis("CMS_2011_S9215166"), _weightMB(0.), _weightDiJet(0.) { }
-
-
- void init() {
- const FinalState fs(-6.0, 6.0, 0.0*GeV);
- addProjection(fs, "FS");
- addProjection(FastJets(fs, FastJets::ANTIKT, 0.5), "Jets");
-
- VetoedFinalState fsv(fs);
- fsv.vetoNeutrinos();
- fsv.addVetoPairDetail(MUON, 0.0*GeV, 99999.9*GeV);
- addProjection(fsv, "fsv");
-
- // For the MB ND selection
- const ChargedFinalState fschrgd(-6.0,6.0,0.0*GeV);
- addProjection(fschrgd, "fschrgd");
- VetoedFinalState fschrgdv(fschrgd);
- fschrgdv.vetoNeutrinos();
- addProjection(fschrgdv, "fschrgdv");
-
- if (fuzzyEquals(sqrtS()/GeV, 900, 1E-3)) {
- _hist_mb = bookHistogram1D(1, 1, 1); // energy flow in MB, 0.9 TeV
- _hist_dijet = bookHistogram1D(2, 1, 1); // energy flow in dijet events, 0.9 TeV
- } else if (fuzzyEquals(sqrtS()/GeV, 7000, 1E-3)) {
- _hist_mb = bookHistogram1D(3, 1, 1); // energy flow in MB, 7 TeV
- _hist_dijet = bookHistogram1D(4, 1, 1); // energy flow in dijet events, 7 TeV
- }
- }
-
-
- void analyze(const Event& event) {
- const double weight = event.weight();
-
- // Skip if the event is empty
- const FinalState& fsv = applyProjection<FinalState>(event, "fsv");
- if (fsv.empty()) vetoEvent;
-
- // Veto diffractive topologies according to defined hadron level
- double count_chrg_forward = 0;
- double count_chrg_backward = 0;
- const FinalState& fschrgdv = applyProjection<FinalState>(event, "fschrgdv");
- foreach (const Particle& p, fschrgdv.particles()) {
- if (3.9 < p.momentum().eta() && p.momentum().eta() < 4.4) count_chrg_forward++;
- if (-4.4 < p.momentum().eta() && p.momentum().eta() < -3.9) count_chrg_backward++;
- }
- if (count_chrg_forward == 0 || count_chrg_backward == 0) vetoEvent;
- /// @todo "Diffractive" veto should really also veto dijet events?
-
-
- // MINIMUM BIAS EVENTS
- _weightMB += weight;
- foreach (const Particle& p, fsv.particles()) {
- _hist_mb->fill(fabs(p.momentum().eta()), weight*p.momentum().E()/GeV);
- }
-
-
- // DIJET EVENTS
- double PTCUT = -1.0;
- if (fuzzyEquals(sqrtS()/GeV, 900, 1E-3)) PTCUT = 8.0*GeV;
- else if (fuzzyEquals(sqrtS()/GeV, 7000, 1E-3)) PTCUT = 20.0*GeV;
- const FastJets& jetpro = applyProjection<FastJets>(event, "Jets");
- const Jets jets = jetpro.jetsByPt(PTCUT);
- if (jets.size() >= 2) {
- // eta cut for the central jets
- if (fabs(jets[0].momentum().eta()) < 2.5 && fabs(jets[1].momentum().eta()) < 2.5) {
- // Back to back condition of the jets
- const double diffphi = deltaPhi(jets[1].momentum().phi(), jets[0].momentum().phi());
- if (diffphi-PI < 1.0) {
- _weightDiJet += weight;
- foreach (const Particle& p, fsv.particles()) {
- _hist_dijet->fill(fabs(p.momentum().eta()), weight*p.momentum().E()/GeV);
- }
- }
- }
- }
-
- }
-
-
- void finalize() {
- scale(_hist_mb , 0.5/_weightMB );
- scale(_hist_dijet, 0.5/_weightDiJet);
- }
-
-
- private:
-
- AIDA::IHistogram1D *_hist_mb, *_hist_dijet;
- double _weightMB,_weightDiJet;
-
- };
-
-
- // Hook for the plugin system
- DECLARE_RIVET_PLUGIN(CMS_2011_S9215166);
-
-}
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/ChargedFinalState.hh"
+#include "Rivet/Projections/FastJets.hh"
+#include "Rivet/Projections/VetoedFinalState.hh"
+
+namespace Rivet {
+
+
+ class CMS_2011_S9215166 : public Analysis {
+ public:
+
+ /// Constructor
+ CMS_2011_S9215166() : Analysis("CMS_2011_S9215166"), _weightMB(0.), _weightDiJet(0.) { }
+
+
+ void init() {
+ const FinalState fs(-6.0, 6.0, 0.0*GeV);
+ addProjection(fs, "FS");
+ addProjection(FastJets(fs, FastJets::ANTIKT, 0.5), "Jets");
+
+ VetoedFinalState fsv(fs);
+ fsv.vetoNeutrinos();
+ fsv.addVetoPairDetail(MUON, 0.0*GeV, 99999.9*GeV);
+ addProjection(fsv, "fsv");
+
+ // For the MB ND selection
+ const ChargedFinalState fschrgd(-6.0,6.0,0.0*GeV);
+ addProjection(fschrgd, "fschrgd");
+ VetoedFinalState fschrgdv(fschrgd);
+ fschrgdv.vetoNeutrinos();
+ addProjection(fschrgdv, "fschrgdv");
+
+ if (fuzzyEquals(sqrtS()/GeV, 900, 1E-3)) {
+ _hist_mb = bookHistogram1D(1, 1, 1); // energy flow in MB, 0.9 TeV
+ _hist_dijet = bookHistogram1D(2, 1, 1); // energy flow in dijet events, 0.9 TeV
+ } else if (fuzzyEquals(sqrtS()/GeV, 7000, 1E-3)) {
+ _hist_mb = bookHistogram1D(3, 1, 1); // energy flow in MB, 7 TeV
+ _hist_dijet = bookHistogram1D(4, 1, 1); // energy flow in dijet events, 7 TeV
+ }
+ }
+
+
+ void analyze(const Event& event) {
+ const double weight = event.weight();
+
+ // Skip if the event is empty
+ const FinalState& fsv = applyProjection<FinalState>(event, "fsv");
+ if (fsv.empty()) vetoEvent;
+
+ // Veto diffractive topologies according to defined hadron level
+ double count_chrg_forward = 0;
+ double count_chrg_backward = 0;
+ const FinalState& fschrgdv = applyProjection<FinalState>(event, "fschrgdv");
+ foreach (const Particle& p, fschrgdv.particles()) {
+ if (3.9 < p.momentum().eta() && p.momentum().eta() < 4.4) count_chrg_forward++;
+ if (-4.4 < p.momentum().eta() && p.momentum().eta() < -3.9) count_chrg_backward++;
+ }
+ if (count_chrg_forward == 0 || count_chrg_backward == 0) vetoEvent;
+ /// @todo "Diffractive" veto should really also veto dijet events?
+
+
+ // MINIMUM BIAS EVENTS
+ _weightMB += weight;
+ foreach (const Particle& p, fsv.particles()) {
+ _hist_mb->fill(fabs(p.momentum().eta()), weight*p.momentum().E()/GeV);
+ }
+
+
+ // DIJET EVENTS
+ double PTCUT = -1.0;
+ if (fuzzyEquals(sqrtS()/GeV, 900, 1E-3)) PTCUT = 8.0*GeV;
+ else if (fuzzyEquals(sqrtS()/GeV, 7000, 1E-3)) PTCUT = 20.0*GeV;
+ const FastJets& jetpro = applyProjection<FastJets>(event, "Jets");
+ const Jets jets = jetpro.jetsByPt(PTCUT);
+ if (jets.size() >= 2) {
+ // eta cut for the central jets
+ if (fabs(jets[0].momentum().eta()) < 2.5 && fabs(jets[1].momentum().eta()) < 2.5) {
+ // Back to back condition of the jets
+ const double diffphi = deltaPhi(jets[1].momentum().phi(), jets[0].momentum().phi());
+ if (diffphi-PI < 1.0) {
+ _weightDiJet += weight;
+ foreach (const Particle& p, fsv.particles()) {
+ _hist_dijet->fill(fabs(p.momentum().eta()), weight*p.momentum().E()/GeV);
+ }
+ }
+ }
+ }
+
+ }
+
+
+ void finalize() {
+ scale(_hist_mb , 0.5/_weightMB );
+ scale(_hist_dijet, 0.5/_weightDiJet);
+ }
+
+
+ private:
+
+ AIDA::IHistogram1D *_hist_mb, *_hist_dijet;
+ double _weightMB,_weightDiJet;
+
+ };
+
+
+ // Hook for the plugin system
+ DECLARE_RIVET_PLUGIN(CMS_2011_S9215166);
+
+}
More information about the Rivet-svn
mailing list
|