Page MenuHomeHEPForge
Files F7877174

No OneTemporary
Actions

View Options

diff --git a/analyses/pluginCMS/CMS_2016_I1459051.cc b/analyses/pluginCMS/CMS_2016_I1459051.cc
--- a/analyses/pluginCMS/CMS_2016_I1459051.cc
+++ b/analyses/pluginCMS/CMS_2016_I1459051.cc
@@ -1,95 +1,95 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/FastJets.hh"
#include "Rivet/Tools/BinnedHistogram.hh"
namespace Rivet {
/// Inclusive jet pT at 13 TeV
class CMS_2016_I1459051 : public Analysis {
public:
/// Constructor
DEFAULT_RIVET_ANALYSIS_CTOR(CMS_2016_I1459051);
/// Book histograms and initialize projections:
void init() {
// Initialize the projections
const FinalState fs;
declare(FastJets(fs, FastJets::ANTIKT, 0.4), "JetsAK4");
declare(FastJets(fs, FastJets::ANTIKT, 0.7), "JetsAK7");
// Book sets of histograms, binned in absolute rapidity
// AK7
_hist_sigmaAK7.addHistogram(0.0, 0.5, bookHisto1D(1, 1, 1));
_hist_sigmaAK7.addHistogram(0.5, 1.0, bookHisto1D(2, 1, 1));
_hist_sigmaAK7.addHistogram(1.0, 1.5, bookHisto1D(3, 1, 1));
_hist_sigmaAK7.addHistogram(1.5, 2.0, bookHisto1D(4, 1, 1));
_hist_sigmaAK7.addHistogram(2.0, 2.5, bookHisto1D(5, 1, 1));
_hist_sigmaAK7.addHistogram(2.5, 3.0, bookHisto1D(6, 1, 1));
_hist_sigmaAK7Forward = bookHisto1D(7, 1, 1);
// AK4
_hist_sigmaAK4.addHistogram(0.0, 0.5, bookHisto1D(8, 1, 1));
_hist_sigmaAK4.addHistogram(0.5, 1.0, bookHisto1D(9, 1, 1));
_hist_sigmaAK4.addHistogram(1.0, 1.5, bookHisto1D(10, 1, 1));
_hist_sigmaAK4.addHistogram(1.5, 2.0, bookHisto1D(11, 1, 1));
_hist_sigmaAK4.addHistogram(2.0, 2.5, bookHisto1D(12, 1, 1));
_hist_sigmaAK4.addHistogram(2.5, 3.0, bookHisto1D(13, 1, 1));
_hist_sigmaAK4Forward = bookHisto1D(14, 1, 1);
}
/// Per-event analysis
void analyze(const Event &event) {
const double weight = event.weight();
// AK4 jets
const FastJets& fjAK4 = applyProjection<FastJets>(event, "JetsAK4");
const Jets& jetsAK4 = fjAK4.jets(Cuts::ptIn(114*GeV, 2200.0*GeV) && Cuts::absrap < 4.7);
for (const Jet& j : jetsAK4) {
_hist_sigmaAK4.fill(j.absrap(), j.pT(), weight);
if (inRange(j.absrap(), 3.2, 4.7)) _hist_sigmaAK4Forward->fill(j.pT(), weight);
}
// AK7 jets
const FastJets& fjAK7 = applyProjection<FastJets>(event, "JetsAK7");
const Jets& jetsAK7 = fjAK7.jets(Cuts::ptIn(114*GeV, 2200.0*GeV) && Cuts::absrap < 4.7);
for (const Jet& j : jetsAK7) {
_hist_sigmaAK7.fill(j.absrap(), j.pT(), weight);
if (inRange(j.absrap(), 3.2, 4.7)) _hist_sigmaAK7Forward->fill(j.pT(), weight);
}
}
// Finalize
void finalize() {
- /// @todo What is the cross-section unit?
- _hist_sigmaAK4.scale(crossSection()/sumOfWeights()/2.0, this);
- _hist_sigmaAK7.scale(crossSection()/sumOfWeights()/2.0, this);
- scale(_hist_sigmaAK4Forward,crossSection()/sumOfWeights()/3.0);
- scale(_hist_sigmaAK7Forward,crossSection()/sumOfWeights()/3.0);
+ /// @note Extra division factor is the *signed* dy, i.e. 2*d|y|
+ _hist_sigmaAK4.scale(crossSection()/picobarn/sumOfWeights()/2.0, this);
+ _hist_sigmaAK7.scale(crossSection()/picobarn/sumOfWeights()/2.0, this);
+ scale(_hist_sigmaAK4Forward,crossSection()/picobarn/sumOfWeights()/3.0);
+ scale(_hist_sigmaAK7Forward,crossSection()/picobarn/sumOfWeights()/3.0);
}
/// @name Histograms
//@{
BinnedHistogram<double> _hist_sigmaAK4;
BinnedHistogram<double> _hist_sigmaAK7;
Histo1DPtr _hist_sigmaAK4Forward;
Histo1DPtr _hist_sigmaAK7Forward;
//@}
};
// This global object acts as a hook for the plugin system.
DECLARE_RIVET_PLUGIN(CMS_2016_I1459051);
}
diff --git a/analyses/pluginCMS/CMS_2016_I1459051.plot b/analyses/pluginCMS/CMS_2016_I1459051.plot
--- a/analyses/pluginCMS/CMS_2016_I1459051.plot
+++ b/analyses/pluginCMS/CMS_2016_I1459051.plot
@@ -1,155 +1,167 @@
# BEGIN PLOT /CMS_2016_I1459051/d01-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK7 jets, $0.0 < |y| < 0.5$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d02-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK7 jets, $0.5 < |y| < 1.0$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d03-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK7 jets, $1.0 < |y| < 1.5$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d04-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK7 jets, $1.5 < |y| < 2.0$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d05-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK7 jets, $2.0 < |y| < 2.5$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d06-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK7 jets, $2.5 < |y| < 3.0$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d07-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK7 jets, $3.2 < |y| < 4.7$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d08-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK4 jets, $0.0 < |y| < 0.5$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d09-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK4 jets, $0.5 < |y| < 1.0$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d10-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK4 jets, $1.0 < |y| < 1.5$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d11-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK4 jets, $1.5 < |y| < 2.0$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d12-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK4 jets, $2.0 < |y| < 2.5$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d13-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK4 jets, $2.5 < |y| < 3.0$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
+
# BEGIN PLOT /CMS_2016_I1459051/d14-x01-y01
Title=CMS, $\sqrt(s)$ = 13 TeV, Inclusive AK4 jets, $3.2 < |y| < 4.7$
XLabel=$p_{T}$ [GeV]
YLabel=d$^2\sigma$/d$p_{T}$d$y$ [pb/GeV]
FullRange=1
LogY=1
LogX=1
LegendYPos=0.9
RatioPlotYMax=2
RatioPlotYMin=0.1
# END PLOT
-

File Metadata

Mime Type
text/x-diff
Expires
Tue, Nov 19, 2:53 PM (1 d, 12 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3804853
Default Alt Text
(7 KB)

Event Timeline