Page Menu
Home
HEPForge
Search
Configure Global Search
Log In
Files
F7877917
No One
Temporary
Actions
View File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Flag For Later
Size
150 KB
Subscribers
None
View Options
Index: trunk/code/ALICE_2015_I1343112.cc
===================================================================
--- trunk/code/ALICE_2015_I1343112.cc (revision 461)
+++ trunk/code/ALICE_2015_I1343112.cc (revision 462)
@@ -1,181 +1,181 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/FastJets.hh"
#include "Rivet/Tools/ParticleIdUtils.hh"
#include <boost/lexical_cast.hpp>
#include "Rivet/Tools/Logging.hh"
/// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder...
namespace Rivet {
class ALICE_2015_I1343112 : public Analysis {
public:
/// @name Constructors etc.
//@{
/// Constructor
ALICE_2015_I1343112()
: Analysis("ALICE_2015_I1343112"),
_hadronetamax(1.), _jetradius(0.2), _jetetamax(0.5), _jetptmin(20.), _jetptmax(120.)
{ }
//@}
public:
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void init() {
_centedges += 0.0, 0.1, 0.3;
_biases += 0., 3., 5., 7., 10.;
_ncentbins = 2;
_nbiasbins = 5;
FinalState fs(-_hadronetamax, _hadronetamax, 0.*GeV);
addProjection(fs, "FS");
FastJets fj(fs, FastJets::ANTIKT, _jetradius);
fj.useInvisibles();
addProjection(fj, "Jets");
_sumwtcentbins[0] = 0.;
_sumwtcentbins[1] = 0.;
_sumwtcentbins[2] = 0.;
_h_jetpt[0] = bookHisto1D(1, 1, 1);
_h_jetpt[1] = bookHisto1D(1, 1, 2);
_h_jetpt[2] = bookHisto1D(1, 1, 3);
_h_jetpt_pp_incl = bookHisto1D("jetpt_pp_incl",refData(2, 1, 1));
_h_jetpt_pp_bias = bookHisto1D("jetpt_pp_bias",refData(2, 1, 1));
_s_jetpt_pp_ratio = bookScatter2D(2, 1, 1);
_h_jetpt_AA_bias[0] = bookHisto1D("jetpt_AA_bias_0",refData(3, 1, 1));
_h_jetpt_AA_bias[1] = bookHisto1D("jetpt_AA_bias_3",refData(3, 1, 2));
_h_jetpt_AA_bias[2] = bookHisto1D("jetpt_AA_bias_5",refData(3, 1, 1));
_h_jetpt_AA_bias[3] = bookHisto1D("jetpt_AA_bias_7",refData(3, 1, 3));
_h_jetpt_AA_bias[4] = bookHisto1D("jetpt_AA_bias_10",refData(3, 1, 4));
_s_jetpt_AA_ratios[0] = bookScatter2D(3, 1, 1);
_s_jetpt_AA_ratios[1] = bookScatter2D(3, 1, 2);
_s_jetpt_AA_ratios[2] = bookScatter2D(3, 1, 3);
_s_jetpt_AA_ratios[3] = bookScatter2D(3, 1, 4);
_h_raa[0] = bookHisto1D(4, 1, 1);
_h_raa[1] = bookHisto1D(4, 1, 2);
}
/// Perform the per-event analysis
void analyze(const Event& event) {
const double weight = event.weight();
const double cent = (event.genEvent()->heavy_ion()?event.genEvent()->heavy_ion()->impact_parameter():-1.);
Cut cuts = Cuts::abseta < _jetetamax && Cuts::pT > _jetptmin*GeV && Cuts::pT < _jetptmax*GeV;
const Jets jets = applyProjection<FastJets>(event, "Jets").jetsByPt(cuts);
if (cent < 0.) _sumwtcentbins[0] += weight;
if (cent < 0. || (cent >= _centedges[0] && cent < _centedges[1])) _sumwtcentbins[1] += weight;
if (cent < 0. || (cent >= _centedges[1] && cent < _centedges[2])) _sumwtcentbins[2] += weight;
for (size_t k = 0; k < _nbiasbins; ++k) {
foreach (Jet jet, jets) {
// track bias
foreach (Particle part, jet.particles()) {
if (PID::charge(part.pdgId()) != 0. && part.momentum().pT() > _biases[k]*GeV) {
double jetpt(jet.momentum().pT());
if (cent < 0.) {
if (k == 0) _h_jetpt_pp_incl->fill(jetpt,weight);
if (k == 2) {
_h_jetpt[0]->fill(jetpt,weight);
_h_jetpt_pp_bias->fill(jetpt,weight);
_h_raa[0]->fill(jetpt,weight);
_h_raa[1]->fill(jetpt,weight);
}
}
if (cent >= _centedges[0] && cent < _centedges[1]) {
_h_jetpt_AA_bias[k]->fill(jetpt,weight);
if (k == 2) {
_h_jetpt[1]->fill(jetpt,weight);
_h_raa[0]->fill(jetpt,weight);
}
}
if (cent >= _centedges[0] && cent < _centedges[2]) {
if (k == 2) {
_h_jetpt[2]->fill(jetpt,weight);
_h_raa[1]->fill(jetpt,weight);
}
}
break;
}
}
}
}
}
/// Normalise histograms etc., after the run
void finalize() {
scale(_h_jetpt[0],1./(_sumwtcentbins[0]>0.?_sumwtcentbins[0]:1.));
scale(_h_jetpt[1],1./(_sumwtcentbins[1]>0.?_sumwtcentbins[1]:1.));
scale(_h_jetpt[2],1./(_sumwtcentbins[2]>0.?_sumwtcentbins[2]:1.));
scale(_h_raa[0],crossSection()/(_sumwtcentbins[0]>0.?_sumwtcentbins[0]:1.));
scale(_h_raa[1],crossSection()/(_sumwtcentbins[1]>0.?_sumwtcentbins[1]:1.));
scale(_h_jetpt_pp_incl,crossSection()/(_sumwtcentbins[0]>0.?_sumwtcentbins[0]:1.));
scale(_h_jetpt_pp_bias,crossSection()/(_sumwtcentbins[0]>0.?_sumwtcentbins[0]:1.));
for (size_t k = 0; k < _nbiasbins; ++k) {
scale(_h_jetpt_AA_bias[k],crossSection()/(_sumwtcentbins[0]>0.?_sumwtcentbins[0]:1.));
}
- divide(_h_jetpt_pp_bias, _h_jetpt_pp_bias, _s_jetpt_pp_ratio);
+ divide(_h_jetpt_pp_bias, _h_jetpt_pp_incl, _s_jetpt_pp_ratio);
divide(_h_jetpt_AA_bias[0], _h_jetpt_AA_bias[2], _s_jetpt_AA_ratios[0]);
divide(_h_jetpt_AA_bias[1], _h_jetpt_AA_bias[2], _s_jetpt_AA_ratios[1]);
divide(_h_jetpt_AA_bias[3], _h_jetpt_AA_bias[2], _s_jetpt_AA_ratios[2]);
divide(_h_jetpt_AA_bias[4], _h_jetpt_AA_bias[2], _s_jetpt_AA_ratios[3]);
getLog() << Log::DEBUG << "sumOfWeights() = " << sumOfWeights() << std::endl;
}
//@}
private:
double _hadronetamax, _jetradius, _jetetamax, _jetptmin, _jetptmax;
double _sumwtcentbins[3];
vector<double> _centedges, _biases;
size_t _ncentbins, _nbiasbins;
private:
/// @name Histograms
//@{
Histo1DPtr _h_jetpt[3];
Histo1DPtr _h_jetpt_pp_incl;
Histo1DPtr _h_jetpt_pp_bias;
Histo1DPtr _h_jetpt_AA_bias[5];
Histo1DPtr _h_raa[2];
Scatter2DPtr _s_jetpt_pp_ratio;
Scatter2DPtr _s_jetpt_AA_ratios[4];
//@}
};
// The hook for the plugin system
DECLARE_RIVET_PLUGIN(ALICE_2015_I1343112);
}
Index: trunk/code/ALICE_2014_I1263194.plot
===================================================================
--- trunk/code/ALICE_2014_I1263194.plot (revision 461)
+++ trunk/code/ALICE_2014_I1263194.plot (revision 462)
@@ -1,6 +1,327 @@
-## BEGIN PLOT /ALICE_2014_I1263194/*
-#RatioPlotYMin=0
-#RatioPlotYMin=1
-#RatioPlotYMax=10
-## END PLOT
+# BEGIN PLOT /ALICE_2014_I1263194/d02-x01-y01
+Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d02-x01-y02
+Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d03-x01-y01
+Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d03-x01-y02
+Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d04-x01-y01
+Title=Pb+Pb 30-50, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d04-x01-y02
+Title=Pb+Pb 30-50, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d05-x01-y01
+Title=Pb+Pb 50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d05-x01-y02
+Title=Pb+Pb 50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d06-x01-y01
+Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d06-x01-y02
+Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d07-x01-y01
+Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d07-x01-y02
+Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+
+# BEGIN PLOT /ALICE_2014_I1263194/d08-x01-y01
+Title=Pb+Pb 30-50, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d08-x01-y02
+Title=Pb+Pb 30-50, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d09-x01-y01
+Title=Pb+Pb 50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d09-x01-y02
+Title=Pb+Pb 50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d10-x01-y01
+Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d10-x01-y02
+Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d11-x01-y01
+Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d11-x01-y02
+Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d12-x01-y01
+Title=Pb+Pb 30-50, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d12-x01-y02
+Title=Pb+Pb 30-50, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d13-x01-y01
+Title=Pb+Pb 50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d13-x01-y02
+Title=Pb+Pb 50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d14-x01-y01
+Title=Pb+Pb 0-10, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=N( PT(leading track) > 0.15 GeV ) / N( PT(leading track) > 5 GeV )
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d14-x01-y02
+Title=Pb+Pb 0-10, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=N( PT(leading track) > 0.15 GeV ) / N( PT(leading track) > 5 GeV )
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d15-x01-y01
+Title=Pb+Pb 50-80, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=N( PT(leading track) > 0.15 GeV ) / N( PT(leading track) > 5 GeV )
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d15-x01-y02
+Title=Pb+Pb 50-80, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=N( PT(leading track) > 0.15 GeV ) / N( PT(leading track) > 5 GeV )
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d16-x01-y01
+Title=Pb+Pb 0-10, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=N( PT(leading track) > 10 GeV ) / N( PT(leading track) > 5 GeV )
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d16-x01-y02
+Title=Pb+Pb 0-10, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=N( PT(leading track) > 10 GeV ) / N( PT(leading track) > 5 GeV )
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d17-x01-y01
+Title=Pb+Pb 50-80, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=N( PT(leading track) > 10 GeV ) / N( PT(leading track) > 5 GeV )
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d17-x01-y02
+Title=Pb+Pb 50-80, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=N( PT(leading track) > 10 GeV ) / N( PT(leading track) > 5 GeV )
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d18-x01-y01
+Title=Pb+Pb 00-10/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d18-x01-y02
+Title=Pb+Pb 00-10/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d19-x01-y01
+Title=Pb+Pb 10-30/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d19-x01-y02
+Title=Pb+Pb 10-30/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d20-x01-y01
+Title=Pb+Pb 30-50/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d20-x01-y02
+Title=Pb+Pb 30-50/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d21-x01-y01
+Title=Pb+Pb 00-10/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d21-x01-y02
+Title=Pb+Pb 00-10/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d22-x01-y01
+Title=Pb+Pb 10-30/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d22-x01-y02
+Title=Pb+Pb 10-30/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d23-x01-y01
+Title=Pb+Pb 30-50/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d23-x01-y02
+Title=Pb+Pb 30-50/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d24-x01-y01
+Title=Pb+Pb 00-10/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d24-x01-y02
+Title=Pb+Pb 00-10/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d25-x01-y01
+Title=Pb+Pb 10-30/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d25-x01-y02
+Title=Pb+Pb 10-30/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d26-x01-y01
+Title=Pb+Pb 30-50/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2014_I1263194/d26-x01-y02
+Title=Pb+Pb 30-50/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{cp}$
+LogY=0
+# END PLOT
Index: trunk/code/ALICE_2014_I1263194.cc
===================================================================
--- trunk/code/ALICE_2014_I1263194.cc (revision 461)
+++ trunk/code/ALICE_2014_I1263194.cc (revision 462)
@@ -1,156 +1,192 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/FastJets.hh"
/// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder...
namespace Rivet {
class ALICE_2014_I1263194 : public Analysis {
public:
/// @name Constructors etc.
//@{
/// Constructor
ALICE_2014_I1263194()
: Analysis("ALICE_2014_I1263194")
{ }
//@}
public:
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void init() {
_centedges += 0.0, 0.1, 0.3, 0.5, 0.8;
_jetparams += 0.2, 0.3;
_biases += 0., 5., 10.;
_ncentbins = 4;
_njetbins = 2;
_nbiasbins = 3;
ChargedFinalState cfs(-0.9, 0.9, 0.15*GeV);
addProjection(cfs, "CFS");
for (size_t j = 0; j < _njetbins; ++j) {
stringstream ss; ss << "Jets" << j;
const string pname = ss.str();
_names[j] = pname;
FastJets fj(cfs, FastJets::ANTIKT, _jetparams[j]);
addProjection(fj, pname);
}
for (size_t i = 0; i < _ncentbins; ++i) {
_sumwtcentbins[i] = 0.;
if (i < _ncentbins-1) _rsumwtcentbins[i] = 0.;
for (size_t j = 0; j < _njetbins; ++j) {
for (size_t k = 0; k < _nbiasbins; ++k) {
_h_jetspec[i][j][k] = bookHisto1D(2+4*k+i, 1, j+1);
if (i < _ncentbins-1) {
_h_rcp[i][j][k] = bookHisto1D(18+3*k+i, 1, j+1);
}
}
}
- }
+ }
+
+ for (size_t i = 0; i < 2; i++) {
+ for (size_t j = 0; j < _njetbins; ++j) {
+ for (size_t k = 0; k < _nbiasbins-1; ++k) {
+ _s_rbias[i][j][k] = bookScatter2D(14+2*k+i, 1, j+1);
+ stringstream ss1; ss1 << "num_"<<i<<"_"<<j<<"_"<<k;
+ stringstream ss2; ss2 << "den_"<<i<<"_"<<j<<"_"<<k;
+ _h_rbiasnum[i][j][k] = bookHisto1D(ss1.str(),refData(14+2*k+i, 1, j+1));
+ _h_rbiasden[i][j][k] = bookHisto1D(ss2.str(),refData(14+2*k+i, 1, j+1));
+ }
+ }
+ }
}
/// Perform the per-event analysis
void analyze(const Event& event) {
const double weight = event.weight();
const double cent = (event.genEvent()->heavy_ion()?event.genEvent()->heavy_ion()->impact_parameter():-1.);
for (size_t i = 0; i < _ncentbins; ++i) {
if (cent < 0. || (cent >= _centedges[i] && cent < _centedges[i+1])) {
_sumwtcentbins[i] += weight;
if (i < _ncentbins-1) _rsumwtcentbins[i] += weight;
else {
for (size_t l = 0; l < _ncentbins-1; ++l) {
_rsumwtcentbins[l] += weight;
}
}
}
}
for (size_t j = 0; j < _njetbins; ++j) {
Cut cuts = Cuts::abseta < 0.5 && Cuts::pT > 10*GeV && Cuts::pT < 120.*GeV;
const Jets jets = applyProjection<FastJets>(event, _names[j]).jetsByPt(cuts);
for (size_t i = 0; i < _ncentbins; ++i) {
if (cent < 0. || (cent >= _centedges[i] && cent < _centedges[i+1])) {
// track bias
for (size_t k = 0; k < _nbiasbins; ++k) {
foreach (Jet jet, jets) {
vector< Particle > & particles(jet.particles());
foreach (Particle part, particles) {
if (part.momentum().pT() > _biases[k]*GeV) {
_h_jetspec[i][j][k]->fill(jet.momentum().pT(),weight);
if (i < _ncentbins-1) _h_rcp[i][j][k]->fill(jet.momentum().pT(),weight);
else {
for (size_t l = 0; l < _ncentbins-1; ++l) {
_h_rcp[l][j][k]->fill(jet.momentum().pT(),weight);
}
- }
+ }
+ if (i == 0 && k == 0) _h_rbiasnum[0][j][0]->fill(jet.momentum().pT(),weight);
+ if (i == _ncentbins-1 && k == 0) _h_rbiasnum[1][j][0]->fill(jet.momentum().pT(),weight);
+ if (i == 0 && k == 1) {
+ _h_rbiasden[0][j][0]->fill(jet.momentum().pT(),weight);
+ _h_rbiasden[0][j][1]->fill(jet.momentum().pT(),weight);
+ }
+ if (i == _ncentbins-1 && k == 0) {
+ _h_rbiasden[1][j][0]->fill(jet.momentum().pT(),weight);
+ _h_rbiasden[1][j][1]->fill(jet.momentum().pT(),weight);
+ }
+ if (i == 0 && k == 2) _h_rbiasnum[0][j][1]->fill(jet.momentum().pT(),weight);
+ if (i == _ncentbins-1 && k == 2) _h_rbiasnum[1][j][1]->fill(jet.momentum().pT(),weight);
break;
}
}
}
}
}
}
}
}
/// Normalise histograms etc., after the run
void finalize() {
for (size_t i = 0; i < _ncentbins; ++i) {
for (size_t j = 0; j < _njetbins; ++j) {
for (size_t k = 0; k < _nbiasbins; ++k) {
scale(_h_jetspec[i][j][k],1.e6/(_sumwtcentbins[i]>0.?_sumwtcentbins[i]:1.));
if (i < _ncentbins-1) scale(_h_rcp[i][j][k],crossSection()/(_rsumwtcentbins[i]>0.?_rsumwtcentbins[i]:1.));
}
}
- }
- }
+ }
+
+ for (size_t i = 0; i < 2; ++i) {
+ for (size_t j = 0; j < _njetbins; ++j) {
+ for (size_t k = 0; k < _nbiasbins-1; ++k) {
+ divide(_h_rbiasnum[i][j][k], _h_rbiasden[i][j][k], _s_rbias[i][j][k]);
+ }
+ }
+ }
+
+ }
//@}
private:
vector<double> _jetparams, _centedges, _biases;
double _sumwtcentbins[4],_rsumwtcentbins[3];
size_t _ncentbins, _njetbins, _nbiasbins;
string _names[2];
private:
/// @name Histograms
//@{
Histo1DPtr _h_jetspec[4][2][3];
- Histo1DPtr _h_rcp[3][2][3];
+ Histo1DPtr _h_rcp[3][2][3];
+ Scatter2DPtr _s_rbias[2][2][2];
+ Histo1DPtr _h_rbiasnum[2][2][2];
+ Histo1DPtr _h_rbiasden[2][2][2];
//@}
};
// The hook for the plugin system
DECLARE_RIVET_PLUGIN(ALICE_2014_I1263194);
}
Index: trunk/code/ATLAS_2013_I1240088.plot
===================================================================
--- trunk/code/ATLAS_2013_I1240088.plot (revision 461)
+++ trunk/code/ATLAS_2013_I1240088.plot (revision 462)
@@ -1,13 +1,295 @@
# BEGIN PLOT /ATLAS_2013_I1240088/*
-#Title=[Uncomment and insert title for histogram d01-x01-y01 here]
-#XLabel=[Uncomment and insert x-axis label for histogram d01-x01-y01 here]
-#YLabel=[Uncomment and insert y-axis label for histogram d01-x01-y01 here]
+YMin=0.5
+YMax=0.8
NormalizeToIntegral=1
-LogY=0
# END PLOT
-# BEGIN PLOT /ATLAS_2013_I1240088/*
-YMin=0.5
-YMax=0.8
+# BEGIN PLOT /ATLAS_2013_I1240088/d01-x01-y01
+Title=$5-10$, $45 < p_\perp < 60\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d01-x02-y01
+Title=$5-10$, $45 < p_\perp < 60\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d01-x01-y02
+Title=$5-10$, $60 < p_\perp < 80\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d01-x02-y02
+Title=$5-10$, $60 < p_\perp < 80\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d01-x01-y03
+Title=$5-10$, $80 < p_\perp < 110\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d01-x02-y03
+Title=$5-10$, $80 < p_\perp < 110\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d01-x01-y04
+Title=$5-10$, $110 < p_\perp < 160\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d01-x02-y04
+Title=$5-10$, $110 < p_\perp < 160\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d02-x01-y01
+Title=$10-20$, $45 < p_\perp < 60\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d02-x02-y01
+Title=$10-20$, $45 < p_\perp < 60\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d02-x01-y02
+Title=$10-20$, $60 < p_\perp < 80\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d02-x02-y02
+Title=$10-20$, $60 < p_\perp < 80\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d02-x01-y03
+Title=$10-20$, $80 < p_\perp < 110\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d02-x02-y03
+Title=$10-20$, $80 < p_\perp < 110\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d02-x01-y04
+Title=$10-20$, $110 < p_\perp < 160\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d02-x02-y04
+Title=$10-20$, $110 < p_\perp < 160\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d03-x01-y01
+Title=$20-30$, $45 < p_\perp < 60\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+PolyMarker=*
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d03-x02-y01
+Title=$20-30$, $45 < p_\perp < 60\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d03-x01-y02
+Title=$20-30$, $60 < p_\perp < 80\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d03-x02-y02
+Title=$20-30$, $60 < p_\perp < 80\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d03-x01-y03
+Title=$20-30$, $80 < p_\perp < 110\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d03-x02-y03
+Title=$20-30$, $80 < p_\perp < 110\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d03-x01-y04
+Title=$20-30$, $110 < p_\perp < 160\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d03-x02-y04
+Title=$20-30$, $110 < p_\perp < 160\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d04-x01-y01
+Title=$30-40$, $45 < p_\perp < 60\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d04-x02-y01
+Title=$30-40$, $45 < p_\perp < 60\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d04-x01-y02
+Title=$30-40$, $60 < p_\perp < 80\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d04-x02-y02
+Title=$30-40$, $60 < p_\perp < 80\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d04-x01-y03
+Title=$30-40$, $80 < p_\perp < 110\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d04-x02-y03
+Title=$30-40$, $80 < p_\perp < 110\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d04-x01-y04
+Title=$30-40$, $110 < p_\perp < 160\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d04-x02-y04
+Title=$30-40$, $110 < p_\perp < 160\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d05-x01-y01
+Title=$40-50$, $45 < p_\perp < 60\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d05-x02-y01
+Title=$40-50$, $45 < p_\perp < 60\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d05-x01-y02
+Title=$40-50$, $60 < p_\perp < 80\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d05-x02-y02
+Title=$40-50$, $60 < p_\perp < 80\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d05-x01-y03
+Title=$40-50$, $80 < p_\perp < 110\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d05-x02-y03
+Title=$40-50$, $80 < p_\perp < 110\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d05-x01-y04
+Title=$40-50$, $110 < p_\perp < 160\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d05-x02-y04
+Title=$40-50$, $110 < p_\perp < 160\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d06-x01-y01
+Title=$50-60$, $45 < p_\perp < 60\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d06-x02-y01
+Title=$50-60$, $45 < p_\perp < 60\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d06-x01-y02
+Title=$50-60$, $60 < p_\perp < 80\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d06-x02-y02
+Title=$50-60$, $60 < p_\perp < 80\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d06-x01-y03
+Title=$50-60$, $80 < p_\perp < 110\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d06-x02-y03
+Title=$50-60$, $80 < p_\perp < 110\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d06-x01-y04
+Title=$50-60$, $110 < p_\perp < 160\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1240088/d06-x02-y04
+Title=$50-60$, $110 < p_\perp < 160\,\text{GeV}$
+XLabel=$\Delta \phi$
+YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$
# END PLOT
Index: trunk/code/CMS_2013_I1256590.yoda
===================================================================
--- trunk/code/CMS_2013_I1256590.yoda (revision 461)
+++ trunk/code/CMS_2013_I1256590.yoda (revision 462)
@@ -1,206 +1,206 @@
-# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d01-x01-y03
+BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d01-x01-y03
Path=/REF/CMS_2013_I1256590/d01-x01-y03
Title=
Type=Scatter2D
# xval xerr- xerr+ yval yerr- yerr+
2.500000e-02 2.500000e-02 2.500000e-02 9.996849e-01 1.534981e-02 1.496040e-02
7.500000e-02 2.500000e-02 2.500000e-02 9.914790e-01 4.180134e-02 4.529186e-02
1.250000e-01 2.500000e-02 2.500000e-02 9.037203e-01 4.181550e-02 3.772316e-02
1.750000e-01 2.500000e-02 2.500000e-02 1.016727e+00 5.312961e-02 5.671926e-02
2.250000e-01 2.500000e-02 2.500000e-02 1.273673e+00 1.439611e-01 1.400741e-01
2.750000e-01 2.500000e-02 2.500000e-02 1.382885e+00 1.060964e-01 1.097498e-01
-# END YODA_SCATTER2D
+END YODA_SCATTER2D
-# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d01-x01-y02
+BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d01-x01-y02
Path=/REF/CMS_2013_I1256590/d01-x01-y02
Title=
Type=Scatter2D
# xval xerr- xerr+ yval yerr- yerr+
2.500000e-02 2.500000e-02 2.500000e-02 1.266784e+01 0.000000e+00 0.000000e+00
7.500000e-02 2.500000e-02 2.500000e-02 4.565196e+00 0.000000e+00 0.000000e+00
1.250000e-01 2.500000e-02 2.500000e-02 1.470847e+00 0.000000e+00 0.000000e+00
1.750000e-01 2.500000e-02 2.500000e-02 6.714706e-01 0.000000e+00 0.000000e+00
2.250000e-01 2.500000e-02 2.500000e-02 3.694601e-01 0.000000e+00 0.000000e+00
2.750000e-01 2.500000e-02 2.500000e-02 2.007718e-01 0.000000e+00 0.000000e+00
-# END YODA_SCATTER2D
+END YODA_SCATTER2D
-# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d01-x01-y01
+BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d01-x01-y01
Path=/REF/CMS_2013_I1256590/d01-x01-y01
Title=
Type=Scatter2D
# xval xerr- xerr+ yval yerr- yerr+
2.500000e-02 2.500000e-02 2.500000e-02 1.266784e+01 0.000000e+00 0.000000e+00
7.500000e-02 2.500000e-02 2.500000e-02 4.565196e+00 0.000000e+00 0.000000e+00
1.250000e-01 2.500000e-02 2.500000e-02 1.298714e+00 0.000000e+00 0.000000e+00
1.750000e-01 2.500000e-02 2.500000e-02 6.883952e-01 0.000000e+00 0.000000e+00
2.250000e-01 2.500000e-02 2.500000e-02 4.738880e-01 0.000000e+00 0.000000e+00
2.750000e-01 2.500000e-02 2.500000e-02 2.774874e-01 0.000000e+00 0.000000e+00
-# END YODA_SCATTER2D
+END YODA_SCATTER2D
-# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d03-x01-y01
+BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d03-x01-y01
Path=/REF/CMS_2013_I1256590/d03-x01-y01
Title=
Type=Scatter2D
# xval xerr- xerr+ yval yerr- yerr+
2.500000e-01 2.500000e-01 2.500000e-01 4.589493e-02 2.640026e-02 2.607364e-02
7.500000e-01 2.500000e-01 2.500000e-01 3.302992e-01 1.372295e-01 1.308755e-01
1.250000e+00 2.500000e-01 2.500000e-01 7.553553e-01 2.064520e-01 2.113501e-01
1.750000e+00 2.500000e-01 2.500000e-01 1.219367e+00 1.948854e-01 1.904400e-01
2.250000e+00 2.500000e-01 2.500000e-01 1.560549e+00 1.710532e-01 1.668595e-01
2.750000e+00 2.500000e-01 2.500000e-01 1.804274e+00 2.437255e-01 2.817904e-01
3.250000e+00 2.500000e-01 2.500000e-01 2.086064e+00 4.080724e-01 3.258002e-01
3.750000e+00 2.500000e-01 2.500000e-01 2.447121e+00 4.787017e-01 4.655116e-01
4.250000e+00 2.500000e-01 2.500000e-01 2.556013e+00 3.771229e-01 3.566198e-01
4.750000e+00 2.500000e-01 2.500000e-01 1.273626e+00 4.731263e-01 4.790329e-01
5.250000e+00 2.500000e-01 2.500000e-01 1.213483e-01 5.143852e-02 1.293548e-01
-# END YODA_SCATTER2D
+END YODA_SCATTER2D
-# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d03-x01-y03
+BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d03-x01-y03
Path=/REF/CMS_2013_I1256590/d03-x01-y03
Title=
Type=Scatter2D
# xval xerr- xerr+ yval yerr- yerr+
2.500000e-01 2.500000e-01 2.500000e-01 9.511278e-01 2.506266e-01 2.631579e-01
7.500000e-01 2.500000e-01 2.500000e-01 9.949875e-01 1.817043e-01 1.879699e-01
1.250000e+00 2.500000e-01 2.500000e-01 8.508772e-01 1.190476e-01 1.190476e-01
1.750000e+00 2.500000e-01 2.500000e-01 8.132832e-01 1.002506e-01 9.398496e-02
2.250000e+00 2.500000e-01 2.500000e-01 7.380952e-01 8.145363e-02 8.145363e-02
2.750000e+00 2.500000e-01 2.500000e-01 7.506266e-01 8.145363e-02 8.145363e-02
3.250000e+00 2.500000e-01 2.500000e-01 8.508772e-01 8.771930e-02 8.771930e-02
3.750000e+00 2.500000e-01 2.500000e-01 1.095238e+00 1.190476e-01 1.127820e-01
4.250000e+00 2.500000e-01 2.500000e-01 1.508772e+00 1.879699e-01 1.879699e-01
4.750000e+00 2.500000e-01 2.500000e-01 1.765664e+00 3.320802e-01 3.320802e-01
5.250000e+00 2.500000e-01 2.500000e-01 2.223058e+00 6.892231e-01 6.892231e-01
-# END YODA_SCATTER2D
+END YODA_SCATTER2D
-# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d03-x01-y02
+BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d03-x01-y02
Path=/REF/CMS_2013_I1256590/d03-x01-y02
Title=
Type=Scatter2D
# xval xerr- xerr+ yval yerr- yerr+
2.500000e-01 2.500000e-01 2.500000e-01 4.793716e-02 0.000000e+00 0.000000e+00
7.500000e-01 2.500000e-01 2.500000e-01 3.302992e-01 0.000000e+00 0.000000e+00
1.250000e+00 2.500000e-01 2.500000e-01 8.990452e-01 0.000000e+00 0.000000e+00
1.750000e+00 2.500000e-01 2.500000e-01 1.494066e+00 0.000000e+00 0.000000e+00
2.250000e+00 2.500000e-01 2.500000e-01 2.147498e+00 0.000000e+00 0.000000e+00
2.750000e+00 2.500000e-01 2.500000e-01 2.447121e+00 0.000000e+00 0.000000e+00
3.250000e+00 2.500000e-01 2.500000e-01 2.447121e+00 0.000000e+00 0.000000e+00
3.750000e+00 2.500000e-01 2.500000e-01 2.243057e+00 0.000000e+00 0.000000e+00
4.250000e+00 2.500000e-01 2.500000e-01 1.702521e+00 0.000000e+00 0.000000e+00
4.750000e+00 2.500000e-01 2.500000e-01 7.231755e-01 0.000000e+00 0.000000e+00
5.250000e+00 2.500000e-01 2.500000e-01 5.462544e-02 0.000000e+00 0.000000e+00
-# END YODA_SCATTER2D
+END YODA_SCATTER2D
-# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d04-x01-y02
+BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d04-x01-y02
Path=/REF/CMS_2013_I1256590/d04-x01-y02
Title=
Type=Scatter2D
# xval xerr- xerr+ yval yerr- yerr+
1.097304e+00 9.730405e-02 9.730405e-02 9.843526e-01 0.000000e+00 0.000000e+00
1.310850e+00 1.162410e-01 1.162410e-01 1.032045e+00 0.000000e+00 0.000000e+00
1.565954e+00 1.388630e-01 1.388630e-01 1.032045e+00 0.000000e+00 0.000000e+00
1.870703e+00 1.658865e-01 1.658865e-01 9.689500e-01 0.000000e+00 0.000000e+00
2.234760e+00 1.981700e-01 1.981700e-01 9.537884e-01 0.000000e+00 0.000000e+00
2.669665e+00 2.367350e-01 2.367350e-01 8.407319e-01 0.000000e+00 0.000000e+00
3.189208e+00 2.828065e-01 2.828065e-01 7.648242e-01 0.000000e+00 0.000000e+00
3.809856e+00 3.378430e-01 3.378430e-01 6.132976e-01 0.000000e+00 0.000000e+00
4.551290e+00 4.035905e-01 4.035905e-01 5.238152e-01 0.000000e+00 0.000000e+00
5.437013e+00 4.821335e-01 4.821335e-01 4.473886e-01 0.000000e+00 0.000000e+00
6.495107e+00 5.759600e-01 5.759600e-01 3.644549e-01 0.000000e+00 0.000000e+00
7.759116e+00 6.880480e-01 6.880480e-01 2.922493e-01 0.000000e+00 0.000000e+00
9.269112e+00 8.219480e-01 8.219480e-01 2.418587e-01 0.000000e+00 0.000000e+00
1.107298e+01 9.819200e-01 9.819200e-01 1.764312e-01 0.000000e+00 0.000000e+00
1.322790e+01 1.173000e+00 1.173000e+00 1.370839e-01 0.000000e+00 0.000000e+00
1.580216e+01 1.401255e+00 1.401255e+00 1.015896e-01 0.000000e+00 0.000000e+00
1.887738e+01 1.673970e+00 1.673970e+00 6.957701e-02 0.000000e+00 0.000000e+00
2.255110e+01 1.999745e+00 1.999745e+00 4.917913e-02 0.000000e+00 0.000000e+00
2.693975e+01 2.388910e+00 2.388910e+00 3.368195e-02 0.000000e+00 0.000000e+00
3.218248e+01 2.853820e+00 2.853820e+00 2.306820e-02 0.000000e+00 0.000000e+00
3.844549e+01 3.409190e+00 3.409190e+00 1.287031e-02 0.000000e+00 0.000000e+00
4.592735e+01 4.072655e+00 4.072655e+00 6.532336e-03 0.000000e+00 0.000000e+00
-# END YODA_SCATTER2D
+END YODA_SCATTER2D
-# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d04-x01-y03
+BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d04-x01-y03
Path=/REF/CMS_2013_I1256590/d04-x01-y03
Title=
Type=Scatter2D
# xval xerr- xerr+ yval yerr- yerr+
1.097304e+00 9.730405e-02 9.730405e-02 1.050239e+00 5.598086e-01 5.645933e-01
1.310850e+00 1.162410e-01 1.162410e-01 6.339713e-01 4.354067e-01 4.354067e-01
1.565954e+00 1.388630e-01 1.388630e-01 6.818182e-01 3.636364e-01 3.684211e-01
1.870703e+00 1.658865e-01 1.658865e-01 4.043062e-01 2.918660e-01 2.918660e-01
2.234760e+00 1.981700e-01 1.981700e-01 2.224880e-01 2.344498e-01 2.392344e-01
2.669665e+00 2.367350e-01 2.367350e-01 5.502392e-02 1.818182e-01 1.866029e-01
3.189208e+00 2.828065e-01 2.828065e-01 -5.502392e-02 1.387560e-01 1.435407e-01
3.809856e+00 3.378430e-01 3.378430e-01 -6.937799e-02 1.052632e-01 1.004785e-01
4.551290e+00 4.035905e-01 4.035905e-01 -9.808612e-02 7.177033e-02 7.177033e-02
5.437013e+00 4.821335e-01 4.821335e-01 -1.028708e-01 5.263158e-02 5.263158e-02
6.495107e+00 5.759600e-01 5.759600e-01 -8.851675e-02 3.827751e-02 3.349282e-02
7.759116e+00 6.880480e-01 6.880480e-01 -8.373206e-02 2.870813e-02 2.392344e-02
9.269112e+00 8.219480e-01 8.219480e-01 -6.937799e-02 1.435407e-02 2.392344e-02
1.107298e+01 9.819200e-01 9.819200e-01 -4.066986e-02 1.913876e-02 1.435407e-02
1.322790e+01 1.173000e+00 1.173000e+00 -3.588517e-02 1.435407e-02 1.435407e-02
1.580216e+01 1.401255e+00 1.401255e+00 -2.631579e-02 9.569378e-03 1.435407e-02
1.887738e+01 1.673970e+00 1.673970e+00 -7.177033e-03 1.435407e-02 4.784689e-03
2.255110e+01 1.999745e+00 1.999745e+00 -7.177033e-03 9.569378e-03 4.784689e-03
2.693975e+01 2.388910e+00 2.388910e+00 -2.392344e-03 0.000000e+00 0.000000e+00
3.218248e+01 2.853820e+00 2.853820e+00 2.392344e-03 0.000000e+00 0.000000e+00
3.844549e+01 3.409190e+00 3.409190e+00 2.392344e-03 0.000000e+00 0.000000e+00
4.592735e+01 4.072655e+00 4.072655e+00 -2.392344e-03 0.000000e+00 0.000000e+00
-# END YODA_SCATTER2D
+END YODA_SCATTER2D
-# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d04-x01-y01
+BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d04-x01-y01
Path=/REF/CMS_2013_I1256590/d04-x01-y01
Title=
Type=Scatter2D
# xval xerr- xerr+ yval yerr- yerr+
1.097304e+00 9.730405e-02 9.730405e-02 2.033383e+00 4.782024e-01 4.627072e-01
1.310850e+00 1.162410e-01 1.162410e-01 1.656449e+00 3.489591e-01 3.451167e-01
1.565954e+00 1.388630e-01 1.388630e-01 1.709530e+00 2.947629e-01 3.238530e-01
1.870703e+00 1.658865e-01 1.658865e-01 1.370839e+00 2.541162e-01 2.341778e-01
2.234760e+00 1.981700e-01 1.981700e-01 1.189440e+00 2.204898e-01 1.813987e-01
2.669665e+00 2.367350e-01 2.367350e-01 8.954777e-01 1.774117e-01 1.696391e-01
3.189208e+00 2.828065e-01 2.828065e-01 7.068302e-01 1.400369e-01 1.339018e-01
3.809856e+00 3.378430e-01 3.378430e-01 5.491943e-01 1.088062e-01 1.040393e-01
4.551290e+00 4.035905e-01 4.035905e-01 4.267141e-01 8.454043e-02 8.890473e-02
5.437013e+00 4.821335e-01 4.821335e-01 3.421737e-01 6.779131e-02 5.845289e-02
6.495107e+00 5.759600e-01 5.759600e-01 2.700890e-01 4.656965e-02 4.613880e-02
7.759116e+00 6.880480e-01 6.880480e-01 2.065706e-01 3.561760e-02 3.528807e-02
9.269112e+00 8.219480e-01 8.219480e-01 1.709530e-01 2.947629e-02 2.607165e-02
1.107298e+01 9.819200e-01 9.819200e-01 1.328274e-01 2.115517e-02 2.269066e-02
1.322790e+01 1.173000e+00 1.173000e+00 1.015896e-01 1.883195e-02 1.924512e-02
1.580216e+01 1.401255e+00 1.401255e+00 7.648242e-02 1.611232e-02 1.593491e-02
1.887738e+01 1.673970e+00 1.673970e+00 5.667932e-02 1.332960e-02 1.073732e-02
2.255110e+01 1.999745e+00 1.999745e+00 4.006265e-02 1.083773e-02 1.149923e-02
2.693975e+01 2.388910e+00 2.388910e+00 2.831749e-02 8.615028e-03 9.295896e-03
3.218248e+01 2.853820e+00 2.853820e+00 1.909070e-02 6.620015e-03 6.266981e-03
3.844549e+01 3.409190e+00 3.409190e+00 1.170828e-02 4.750580e-03 4.856209e-03
4.592735e+01 4.072655e+00 4.072655e+00 6.636175e-03 2.933691e-03 3.207351e-03
-# END YODA_SCATTER2D
+END YODA_SCATTER2D
-# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d02-x01-y01
+BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d02-x01-y01
Path=/REF/CMS_2013_I1256590/d02-x01-y01
Title=
Type=Scatter2D
# xval xerr- xerr+ yval yerr- yerr+
2.500000e-02 2.500000e-02 2.500000e-02 6.375449e-01 0.000000e+00 0.000000e+00
7.500000e-02 2.500000e-02 2.500000e-02 8.643001e-01 0.000000e+00 0.000000e+00
1.250000e-01 2.500000e-02 2.500000e-02 9.278651e-01 0.000000e+00 0.000000e+00
1.750000e-01 2.500000e-02 2.500000e-02 9.612714e-01 0.000000e+00 0.000000e+00
2.250000e-01 2.500000e-02 2.500000e-02 9.858084e-01 0.000000e+00 0.000000e+00
2.750000e-01 2.500000e-02 2.500000e-02 9.979320e-01 0.000000e+00 0.000000e+00
-# END YODA_SCATTER2D
+END YODA_SCATTER2D
-# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d02-x01-y02
+BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d02-x01-y02
Path=/REF/CMS_2013_I1256590/d02-x01-y02
Title=
Type=Scatter2D
# xval xerr- xerr+ yval yerr- yerr+
2.500000e-02 2.500000e-02 2.500000e-02 6.357711e-01 0.000000e+00 0.000000e+00
7.500000e-02 2.500000e-02 2.500000e-02 8.678478e-01 0.000000e+00 0.000000e+00
1.250000e-01 2.500000e-02 2.500000e-02 9.366820e-01 0.000000e+00 0.000000e+00
1.750000e-01 2.500000e-02 2.500000e-02 9.700882e-01 0.000000e+00 0.000000e+00
2.250000e-01 2.500000e-02 2.500000e-02 9.892897e-01 0.000000e+00 0.000000e+00
2.750000e-01 2.500000e-02 2.500000e-02 9.978656e-01 0.000000e+00 0.000000e+00
-# END YODA_SCATTER2D
+END YODA_SCATTER2D
Index: trunk/code/jetspectrum.cc
===================================================================
--- trunk/code/jetspectrum.cc (revision 461)
+++ trunk/code/jetspectrum.cc (revision 462)
@@ -1,148 +1,197 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/FastJets.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Tools/ParticleIdUtils.hh"
#include <boost/lexical_cast.hpp>
namespace Rivet {
/// @brief jet pt spectrum
class JETSPECTRUM : public Analysis {
public:
/// Constructor
JETSPECTRUM()
: Analysis("JETSPECTRUM"),
_sumWeightSelected1(0.0),_sumWeightSelected2(0.0)
{
setNeedsCrossSection(true);
}
/// Book projections and histograms
void init() {
ChargedFinalState cfs(-0.9, 0.9, 0.150*GeV);
+ FinalState fs(-3., 3., 0.);
addProjection(cfs, "CFS");
addProjection(FastJets(cfs, FastJets::ANTIKT, 0.2), "Jets1");
addProjection(FastJets(cfs, FastJets::ANTIKT, 0.3), "Jets2");
+ addProjection(FastJets(fs, FastJets::ANTIKT, 0.4), "Jets3");
pt_02 = bookHisto1D("jetspec1", 24, 30.0, 150.0,
"jet spectrum R=0.2",
"$p_\\perp\\ [GeV]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
pt_03 = bookHisto1D("jetspec2", 24, 30.0, 150.0,
"jet spectrum R=0.3",
"$p_\\perp\\ [GeV]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
n_02 = bookHisto1D("njet1", 20, 0.0, 20.0,
"number of jet2 R=0.2",
"$N_\\mathrm{jet}$",
"$\\mathrm{d}\\sigma/\\mathrm{d}N_\\mathrm{jet}$");
n_03 = bookHisto1D("njet2", 20, 0.0, 20.0,
"number of jet2 R=0.3",
"$N_\\mathrm{jet}$",
"$\\mathrm{d}\\sigma/\\mathrm{d}N_\\mathrm{jet}$");
+ _h_jetpt_05 = bookHisto1D("jetpt_eta05",45 , 50.0, 500.0,
+ "jet spectrum R=0.4 |\\eta|<0.5",
+ "$p_\\perp\\ [GeV]$",
+ "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
+ _h_jetpt_1 = bookHisto1D("jetpt_eta1", 45, 50.0, 500.0,
+ "jet spectrum R=0.4 |\\eta|<1",
+ "$p_\\perp\\ [GeV]$",
+ "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
+ _h_jetpt_2 = bookHisto1D("jetpt_eta2", 45, 50.0, 500.0,
+ "jet spectrum R=0.4 |\\eta|<0.5",
+ "$p_\\perp\\ [GeV]$",
+ "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
+ _h_jeteta_50 = bookHisto1D("jeteta_pt50",30 , 0.0, 3.0,
+ "jet eta R=0.4 $p_\\perp > 50$ GeV",
+ "$|\\eta|$",
+ "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$");
+ _h_jeteta_100 = bookHisto1D("jeteta_pt100",30 , 0.0, 3.0,
+ "jet eta R=0.4 $p_\\perp > 100$ GeV",
+ "$|\\eta|$",
+ "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$");
+ _h_jeteta_250 = bookHisto1D("jeteta_pt250",30 , 0.0, 3.0,
+ "jet eta R=0.4 $p_\\perp > 250$ GeV",
+ "$|\\eta|$",
+ "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$");
}
/// Do the analysis
void analyze(const Event& event) {
const double weight = event.weight();
// const double sigma = 15.;
// std::cout<<std::endl<<"new event "<<std::endl;
const FastJets& fastjets1 = applyProjection<FastJets>(event, "Jets1");
/* std::list<Jet> jetlist1;
foreach (Jet jet , fastjets1.jets()){
// std::cout<<"Et before: "<<jet.EtSum()<<" , "<<jet.momentum().Et()<<std::endl;
double r1(1.0*rand()/RAND_MAX), r2(1.0*rand()/RAND_MAX);
double fluc(max(sqrt(-2.*log(r1))*cos(2.*M_PI*r2)*sigma,-jet.totalEnergy()));
const FourMomentum mom(fluc,0,0,0);
jet.addParticle(mom);
jetlist1.push_back(jet);
// std::cout<<"Et after: "<<jet.EtSum()<<" , "<<jet.momentum().Et()<<std::endl;
}
// std::cout<<std::endl<<jetlist.front().EtSum()<<std::endl;
jetlist1.sort(cmpJetsByPt);
// std::cout<<jetlist.front().EtSum()<<std::endl<<std::endl;*/
const Jets alljets1 = fastjets1.jetsByPt();
Jets jets1;
foreach (Jet jet, alljets1) {
// foreach (Jet jet, jetlist1) {
// std::cout<<"Et final: "<<jet.EtSum()<<" , "<<jet.momentum().Et()<<std::endl;
if (fabs(jet.momentum().eta())<0.5 && jet.momentum().Et()>30.)
jets1.push_back(jet);
}
const FastJets& fastjets2 = applyProjection<FastJets>(event, "Jets2");
/* std::list<Jet> jetlist2;
foreach (Jet jet , fastjets2.jets()){
// std::cout<<"Et before: "<<jet.EtSum()<<" , "<<jet.momentum().Et()<<std::endl;
double r1(1.0*rand()/RAND_MAX), r2(1.0*rand()/RAND_MAX);
double fluc(max(sqrt(-2.*log(r1))*cos(2.*M_PI*r2)*sigma,-jet.totalEnergy()));
const FourMomentum mom(fluc,0,0,0);
jet.addParticle(mom);
jetlist2.push_back(jet);
// std::cout<<"Et after: "<<jet.EtSum()<<" , "<<jet.momentum().Et()<<std::endl;
}
// std::cout<<std::endl<<jetlist.front().EtSum()<<std::endl;
jetlist2.sort(cmpJetsByPt);
// std::cout<<jetlist.front().EtSum()<<std::endl<<std::endl;*/
const Jets alljets2 = fastjets2.jetsByPt();
Jets jets2;
foreach (Jet jet, alljets2) {
// foreach (Jet jet, jetlist2) {
// std::cout<<"Et final: "<<jet.EtSum()<<" , "<<jet.momentum().Et()<<std::endl;
if (fabs(jet.momentum().eta())<0.5 && jet.momentum().Et()>30.)
jets2.push_back(jet);
}
n_02->fill(jets1.size(),weight);
if (jets1.size() > 0) {
_sumWeightSelected1 += event.weight();
foreach (Jet jet, jets1){
pt_02->fill(jet.momentum().pT(),weight);
}
}
n_03->fill(jets2.size(),weight);
if (jets2.size() > 0) {
_sumWeightSelected2 += event.weight();
foreach (Jet jet, jets2){
pt_03->fill(jet.momentum().pT(),weight);
}
}
+
+ const FastJets& fastjets3 = applyProjection<FastJets>(event, "Jets3");
+ const Jets alljets3 = fastjets3.jetsByPt();
+ foreach (Jet jet, alljets3) {
+ double pt = jet.momentum().pT();
+ double eta = jet.momentum().abseta();
+ if (eta < 0.5) _h_jetpt_05->fill(pt, weight);
+ if (eta < 1.0) _h_jetpt_1->fill(pt, weight);
+ if (eta < 2.0) _h_jetpt_2->fill(pt, weight);
+ if (pt > 50.) _h_jeteta_50->fill(eta, weight);
+ if (pt > 100.) _h_jeteta_100->fill(eta, weight);
+ if (pt > 250.) _h_jeteta_250->fill(eta, weight);
+
+
+ }
}
/// Finalize
void finalize() {
/* scale(pt_02,1./_sumWeightSelected1);
scale(pt_03,1./_sumWeightSelected2);*/
scale(pt_02,crossSection()/sumOfWeights());
scale(pt_03,crossSection()/sumOfWeights());
scale(n_02,1./sumOfWeights());
scale(n_03,1./sumOfWeights());
+ scale(_h_jetpt_05,crossSection()/sumOfWeights());
+ scale(_h_jetpt_1,crossSection()/sumOfWeights());
+ scale(_h_jetpt_2,crossSection()/sumOfWeights());
+ scale(_h_jeteta_50,crossSection()/sumOfWeights());
+ scale(_h_jeteta_100,crossSection()/sumOfWeights());
+ scale(_h_jeteta_250,crossSection()/sumOfWeights());
getLog() << Log::DEBUG << "sumOfWeights() = " << sumOfWeights() << std::endl;
}
private:
double _sumWeightSelected1;
double _sumWeightSelected2;
Histo1DPtr pt_02;
Histo1DPtr pt_03;
Histo1DPtr n_02;
Histo1DPtr n_03;
+ Histo1DPtr _h_jetpt_05, _h_jetpt_1, _h_jetpt_2;
+ Histo1DPtr _h_jeteta_50, _h_jeteta_100, _h_jeteta_250;
};
// This global object acts as a hook for the plugin system
AnalysisBuilder<JETSPECTRUM> plugin_JETSPECTRUM;
}
Index: trunk/code/ALICE_GIRTH_4MOMSUB.plot
===================================================================
--- trunk/code/ALICE_GIRTH_4MOMSUB.plot (revision 461)
+++ trunk/code/ALICE_GIRTH_4MOMSUB.plot (revision 462)
@@ -1,7 +1,10 @@
# BEGIN PLOT /ALICE_GIRTH_4MOMSUB/d01-x01-y01
LogY=0
+Title=Pb+Pb $(0-10\%)$ (2.76 TeV)
+XLabel=$g$
+YLabel=$(1/N_{\text{jets}}) \ \mathrm{d}\,N/\mathrm{d}\,g$
# END PLOT
# BEGIN PLOT /ALICE_GIRTH_4MOMSUB/pp_g_*
LogY=0
# END PLOT
Index: trunk/code/ATLAS_2013_I1228693.plot
===================================================================
--- trunk/code/ATLAS_2013_I1228693.plot (revision 461)
+++ trunk/code/ATLAS_2013_I1228693.plot (revision 462)
@@ -1,8 +1,84 @@
# BEGIN PLOT /ATLAS_2013_I1228693/d01-x01-y01
-#Title=[Uncomment and insert title for histogram d01-x01-y01 here]
-#XLabel=[Uncomment and insert x-axis label for histogram d01-x01-y01 here]
-#YLabel=[Uncomment and insert y-axis label for histogram d01-x01-y01 here]
-# + any additional plot settings you might like, see make-plots documentation
+Title=$|y| < 0.3$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d02-x01-y01
+Title=$0.3 < |y| < 0.8$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d03-x01-y01
+Title=$0.8 < |y| < 1.2$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d04-x01-y01
+Title=$1.2 < |y| < 2.1$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d05-x01-y01
+Title=$2.1 < |y| < 2.8$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d06-x01-y01
+Title=$2.8 < |y| < 3.6$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d07-x01-y01
+Title=$3.6 < |y| < 4.4$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d08-x01-y01
+Title=$|y| < 0.3$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d09-x01-y01
+Title=$0.3 < |y| < 0.8$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d10-x01-y01
+Title=$0.8 < |y| < 1.2$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d11-x01-y01
+Title=$1.2 < |y| < 2.1$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d12-x01-y01
+Title=$2.1 < |y| < 2.8$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d13-x01-y01
+Title=$2.8 < |y| < 3.6$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2013_I1228693/d14-x01-y01
+Title=$3.6 < |y| < 4.4$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$
# END PLOT
-# ... add more histograms as you need them ...
Index: trunk/code/hadronspectra.cc
===================================================================
--- trunk/code/hadronspectra.cc (revision 461)
+++ trunk/code/hadronspectra.cc (revision 462)
@@ -1,324 +1,354 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/IdentifiedFinalState.hh"
namespace Rivet {
/// @brief hadron spectra in pp
class HADRONSPECTRA : public Analysis {
public:
/// Constructor
HADRONSPECTRA()
: Analysis("HADRONSPECTRA"),
_sumWeightSelected(0.0)
{
setNeedsCrossSection(true);
}
/// Book projections and histograms
void init() {
IdentifiedFinalState pi0fs(-10., 10., 0.*GeV);
pi0fs.acceptId(111);
addProjection(pi0fs, "PI0FS");
IdentifiedFinalState pifs(-10., 10., 0.*GeV);
pifs.acceptId(211);
addProjection(pifs, "PIFS");
IdentifiedFinalState kfs(-10., 10., 0.*GeV);
kfs.acceptId(130);
kfs.acceptId(310);
kfs.acceptId(311);
kfs.acceptId(321);
addProjection(kfs, "KFS");
IdentifiedFinalState pfs(-10., 10., 0.*GeV);
pfs.acceptId(2212);
addProjection(pfs, "PFS");
ChargedFinalState cfs(-10., 10., 0.*GeV);
addProjection(cfs, "CFS");
FinalState fs(-10., 10., 0.*GeV);
addProjection(fs, "FS");
vector<double> ptbins;
ptbins += 12.0, 14.4,19.2, 24.0, 28.8, 35.2, 41.6, 48.0, 60.8,73.6,86.4,103.6,120.8,140,165,250,400;
_h_pT_pi0_035 = bookHisto1D("pi0ptspec", 50, 0.0, 200.0,
"pi0 $p_\\perp$ spectrum $|\\eta|<0.35$",
"$p_\\perp\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
_h_pT_ch_08 = bookHisto1D("chptspec08", 50, 0.0, 200.0,
"charged hadron $p_\\perp$ spectrum $|\\eta|<0.8$",
"$p_\\perp\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
_h_pT_ch_1 = bookHisto1D("chptspec1", 50, 0.0, 200.0,
"charged hadron $p_\\perp$ spectrum $|\\eta|<1.$",
"$p_\\perp\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
_h_pT_ch_2 = bookHisto1D("chptspec2", 50, 0.0, 200.0,
"charged hadron $p_\\perp$ spectrum $|\\eta|<2.$",
"$p_\\perp\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
_h_pT_ch_3 = bookHisto1D("chptspec3", 50, 0.0, 200.0,
"charged hadron $p_\\perp$ spectrum $|\\eta|<3.$",
"$p_\\perp\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
_h_pT_ch_10 = bookHisto1D("chptspec10", 50, 0.0, 200.0,
"charged hadron $p_\\perp$ spectrum $|\\eta|<10.$",
"$p_\\perp\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
_h_E_ch_1 = bookHisto1D("chespec1", 50, 0.0, 1000.0,
"charged hadron energy spectrum $|\\eta|<1.$",
"$E\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}E$");
_h_E_ch_2 = bookHisto1D("chespec2", 50, 0.0, 1000.0,
"charged hadron energy spectrum $|\\eta|<2.$",
"$E\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}E$");
_h_E_ch_3 = bookHisto1D("chespec3", 50, 0.0, 1000.0,
"charged hadron energy spectrum $|\\eta|<3.$",
"$E\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}E$");
_h_E_ch_10 = bookHisto1D("chespec10", 50, 0.0, 1000.0,
"charged hadron energy spectrum $|\\eta|<10.$",
"$E\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}E$");
_h_pT_ch_1h = bookHisto1D("chptspec1high", 40, 0.0, 1000.0,
"charged hadron $p_\\perp$ spectrum $|\\eta|<1.$",
"$p_\\perp\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
_h_pT_ch_1h2 = bookHisto1D("chptspec1high2", 20, 112.2, 1000.0,
"charged hadron $p_\\perp$ spectrum $|\\eta|<1.$",
"$p_\\perp\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
_h_pT_ch_2h = bookHisto1D("chptspec2high", 40, 0.0, 1000.0,
"charged hadron $p_\\perp$ spectrum $|\\eta|<2.$",
"$p_\\perp\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
_h_pT_ch_3h = bookHisto1D("chptspec3high", 40, 0.0, 1000.0,
"charged hadron $p_\\perp$ spectrum $|\\eta|<3.$",
"$p_\\perp\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
_h_pT_ch_10h = bookHisto1D("chptspec10high", 40, 0.0, 1000.0,
"charged hadron $p_\\perp$ spectrum $|\\eta|<10.$",
"$p_\\perp\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$");
_h_E_ch_1h = bookHisto1D("chespec1high", 100, 0.0, 2000.0,
"charged hadron energy spectrum $|\\eta|<1.$",
"$E\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}E$");
_h_E_ch_2h = bookHisto1D("chespec2high", 100, 0.0, 2000.0,
"charged hadron energy spectrum $|\\eta|<2.$",
"$E\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}E$");
_h_E_ch_3h = bookHisto1D("chespec3high", 100, 0.0, 2000.0,
"charged hadron energy spectrum $|\\eta|<3.$",
"$E\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}E$");
_h_E_ch_10h = bookHisto1D("chespec10high", 100, 0.0, 2000.0,
"charged hadron energy spectrum $|\\eta|<10.$",
"$E\\ [\\mathrm{GeV}]$",
"$\\mathrm{d}\\sigma/\\mathrm{d}E$");
+ _h_abseta_pt1 = bookHisto1D("abseta_pt1", 30, 0.0, 3.0,
+ "hadrons $\\p_\\perp > 1\\ \\mathrm{GeV}$",
+ "$|\\eta|$",
+ "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$");
+
+ _h_abseta_pt10 = bookHisto1D("abseta_pt10", 30, 0.0, 3.0,
+ "hadrons $\\p_\\perp > 10\\ \\mathrm{GeV}$",
+ "$|\\eta|$",
+ "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$");
+
+ _h_abseta_pt50 = bookHisto1D("abseta_pt50", 30, 0.0, 3.0,
+ "hadrons $\\p_\\perp > 50\\ \\mathrm{GeV}$",
+ "$|\\eta|$",
+ "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$");
+
+ _h_abseta_pt100 = bookHisto1D("abseta_pt100", 30, 0.0, 3.0,
+ "hadrons $\\p_\\perp > 100\\ \\mathrm{GeV}$",
+ "$|\\eta|$",
+ "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$");
+
_h_chpT1 = bookHisto1D("chpT1", ptbins);
_h_chpT10 = bookHisto1D("chpT10", ptbins);
_h_allpT1 = bookHisto1D("allpT1", ptbins);
_h_allpT10 = bookHisto1D("allpT10", ptbins);
_h_pi0pT1 = bookHisto1D("pi0pT1", ptbins);
_h_pi0pT10 = bookHisto1D("pi0pT10", ptbins);
_h_pipT1 = bookHisto1D("pipT1", ptbins);
_h_pipT10 = bookHisto1D("pipT10", ptbins);
_h_KpT1 = bookHisto1D("KpT1", ptbins);
_h_KpT10 = bookHisto1D("KpT10", ptbins);
_h_ppT1 = bookHisto1D("ppT1", ptbins);
_h_ppT10 = bookHisto1D("ppT10", ptbins);
}
/// Do the analysis
void analyze(const Event& event) {
const double weight = event.weight();
const IdentifiedFinalState& pi0fs = applyProjection<IdentifiedFinalState>(event, "PI0FS");
foreach (const Particle& p, pi0fs.particles()) {
const double pT = p.momentum().pT() / GeV;
if (fabs(p.momentum().eta()) < 0.35) _h_pT_pi0_035->fill(pT, weight);
if (fabs(p.momentum().eta()) < 1.) _h_pi0pT1->fill(pT, weight);
_h_pi0pT10->fill(pT, weight);
}
const IdentifiedFinalState& pifs = applyProjection<IdentifiedFinalState>(event, "PIFS");
foreach (const Particle& p, pifs.particles()) {
const double pT = p.momentum().pT() / GeV;
if (fabs(p.momentum().eta()) < 1.) _h_pipT1->fill(pT, weight);
_h_pipT10->fill(pT, weight);
}
const IdentifiedFinalState& kfs = applyProjection<IdentifiedFinalState>(event, "KFS");
foreach (const Particle& p, kfs.particles()) {
const double pT = p.momentum().pT() / GeV;
if (fabs(p.momentum().eta()) < 1.) _h_KpT1->fill(pT, weight);
_h_KpT10->fill(pT, weight);
}
const IdentifiedFinalState& pfs = applyProjection<IdentifiedFinalState>(event, "PFS");
foreach (const Particle& p, pfs.particles()) {
const double pT = p.momentum().pT() / GeV;
if (fabs(p.momentum().eta()) < 1.) _h_ppT1->fill(pT, weight);
_h_ppT10->fill(pT, weight);
}
const ChargedFinalState& cfs = applyProjection<ChargedFinalState>(event, "CFS");
foreach (const Particle& p, cfs.particles()) {
if (fabs(p.momentum().eta()) < 0.8) {
const double pT = p.momentum().pT() / GeV;
_h_pT_ch_08->fill(pT, weight);
}
if (fabs(p.momentum().eta()) < 1.) {
const double pT = p.momentum().pT() / GeV;
const double E = p.momentum().E() / GeV;
_h_pT_ch_1->fill(pT, weight);
_h_E_ch_1->fill(E, weight);
_h_pT_ch_1h->fill(pT, weight);
_h_pT_ch_1h2->fill(pT, weight/pT);
_h_E_ch_1h->fill(E, weight);
_h_chpT1->fill(pT,weight);
}
if (fabs(p.momentum().eta()) < 2.) {
const double pT = p.momentum().pT() / GeV;
const double E = p.momentum().E() / GeV;
_h_pT_ch_2->fill(pT, weight);
_h_E_ch_2->fill(E, weight);
_h_pT_ch_2h->fill(pT, weight);
_h_E_ch_2h->fill(E, weight);
}
if (fabs(p.momentum().eta()) < 3.) {
const double pT = p.momentum().pT() / GeV;
const double E = p.momentum().E() / GeV;
_h_pT_ch_3->fill(pT, weight);
_h_E_ch_3->fill(E, weight);
_h_pT_ch_3h->fill(pT, weight);
_h_E_ch_3h->fill(E, weight);
}
if (fabs(p.momentum().eta()) < 10.) {
const double pT = p.momentum().pT() / GeV;
const double E = p.momentum().E() / GeV;
_h_pT_ch_10->fill(pT, weight);
_h_E_ch_10->fill(E, weight);
_h_pT_ch_10h->fill(pT, weight);
_h_E_ch_10h->fill(E, weight);
_h_chpT10->fill(pT,weight);
}
}
const FinalState& fs = applyProjection<FinalState>(event, "FS");
foreach (const Particle& p, fs.particles()) {
const double pT = p.momentum().pT() / GeV;
- if (fabs(p.momentum().eta()) < 1.) _h_allpT1->fill(pT,weight);
- if (fabs(p.momentum().eta()) < 10.) _h_allpT10->fill(pT,weight);
+ const double eta = p.momentum().abseta();
+ if (eta < 1.) _h_allpT1->fill(pT,weight);
+ if (eta < 10.) _h_allpT10->fill(pT,weight);
+ if (pT > 1.) _h_abseta_pt1->fill(eta,weight);
+ if (pT > 10.) _h_abseta_pt10->fill(eta,weight);
+ if (pT > 50.) _h_abseta_pt50->fill(eta,weight);
+ if (pT > 100.) _h_abseta_pt100->fill(eta,weight);
}
_sumWeightSelected += event.weight();
}
/// Finalize
void finalize() {
scale(_h_pT_pi0_035,crossSection()/(_sumWeightSelected));
scale(_h_pT_ch_08,crossSection()/(_sumWeightSelected));
scale(_h_pT_ch_1,crossSection()/(_sumWeightSelected));
scale(_h_pT_ch_2,crossSection()/(_sumWeightSelected));
scale(_h_pT_ch_3,crossSection()/(_sumWeightSelected));
scale(_h_pT_ch_10,crossSection()/(_sumWeightSelected));
scale(_h_E_ch_1,crossSection()/(_sumWeightSelected));
scale(_h_E_ch_2,crossSection()/(_sumWeightSelected));
scale(_h_E_ch_3,crossSection()/(_sumWeightSelected));
scale(_h_E_ch_10,crossSection()/(_sumWeightSelected));
scale(_h_pT_ch_1h,crossSection()/(_sumWeightSelected));
scale(_h_pT_ch_1h2,1./(2.*M_PI*_sumWeightSelected));
scale(_h_pT_ch_2h,crossSection()/(_sumWeightSelected));
scale(_h_pT_ch_3h,crossSection()/(_sumWeightSelected));
scale(_h_pT_ch_10h,crossSection()/(_sumWeightSelected));
scale(_h_E_ch_1h,crossSection()/(_sumWeightSelected));
scale(_h_E_ch_2h,crossSection()/(_sumWeightSelected));
scale(_h_E_ch_3h,crossSection()/(_sumWeightSelected));
scale(_h_E_ch_10h,crossSection()/(_sumWeightSelected));
scale(_h_chpT1,crossSection()/(_sumWeightSelected));
scale(_h_chpT10,crossSection()/(_sumWeightSelected));
scale(_h_allpT1,crossSection()/(_sumWeightSelected));
scale(_h_allpT10,crossSection()/(_sumWeightSelected));
scale(_h_pi0pT1,crossSection()/(_sumWeightSelected));
scale(_h_pi0pT10,crossSection()/(_sumWeightSelected));
scale(_h_pipT1,crossSection()/(_sumWeightSelected));
scale(_h_pipT10,crossSection()/(_sumWeightSelected));
scale(_h_KpT1,crossSection()/(_sumWeightSelected));
scale(_h_KpT10,crossSection()/(_sumWeightSelected));
scale(_h_ppT1,crossSection()/(_sumWeightSelected));
scale(_h_ppT10,crossSection()/(_sumWeightSelected));
+ scale(_h_abseta_pt1,crossSection()/(_sumWeightSelected));
+ scale(_h_abseta_pt10,crossSection()/(_sumWeightSelected));
+ scale(_h_abseta_pt50,crossSection()/(_sumWeightSelected));
+ scale(_h_abseta_pt100,crossSection()/(_sumWeightSelected));
getLog() << Log::DEBUG << "sumOfWeights() = " << sumOfWeights() << std::endl;
getLog() << Log::DEBUG << "_sumWeightSelected = " << _sumWeightSelected << std::endl;
}
private:
double _sumWeightSelected;
Histo1DPtr _h_pT_pi0_035;
Histo1DPtr _h_pT_ch_08;
Histo1DPtr _h_pT_ch_1;
Histo1DPtr _h_pT_ch_2;
Histo1DPtr _h_pT_ch_3;
Histo1DPtr _h_pT_ch_10;
Histo1DPtr _h_E_ch_1;
Histo1DPtr _h_E_ch_2;
Histo1DPtr _h_E_ch_3;
Histo1DPtr _h_E_ch_10;
Histo1DPtr _h_pT_ch_1h;
Histo1DPtr _h_pT_ch_1h2;
Histo1DPtr _h_pT_ch_2h;
Histo1DPtr _h_pT_ch_3h;
Histo1DPtr _h_pT_ch_10h;
Histo1DPtr _h_E_ch_1h;
Histo1DPtr _h_E_ch_2h;
Histo1DPtr _h_E_ch_3h;
Histo1DPtr _h_E_ch_10h;
Histo1DPtr _h_chpT1;
Histo1DPtr _h_chpT10;
Histo1DPtr _h_allpT1;
Histo1DPtr _h_allpT10;
Histo1DPtr _h_pi0pT1;
Histo1DPtr _h_pi0pT10;
Histo1DPtr _h_pipT1;
Histo1DPtr _h_pipT10;
Histo1DPtr _h_KpT1;
Histo1DPtr _h_KpT10;
Histo1DPtr _h_ppT1;
Histo1DPtr _h_ppT10;
+ Histo1DPtr _h_abseta_pt1, _h_abseta_pt10, _h_abseta_pt50, _h_abseta_pt100;
};
// This global object acts as a hook for the plugin system
AnalysisBuilder<HADRONSPECTRA> plugin_HADRONSPECTRA;
}
Index: trunk/code/ALICE_2017_I1512107.cc
===================================================================
--- trunk/code/ALICE_2017_I1512107.cc (revision 461)
+++ trunk/code/ALICE_2017_I1512107.cc (revision 462)
@@ -1,515 +1,520 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/FastJets.hh"
#include "fastjet/ClusterSequence.hh"
#include "fastjet/ClusterSequence.hh"
+#ifndef MODE
+#warning MODE not defined, using default
+#define MODE 0
+#endif
+
namespace Rivet {
using namespace fastjet;
#if MODE==0
class ALICE_2017_I1512107_4MOMSUB : public Analysis {
#elif MODE==1
class ALICE_2017_I1512107_GRIDSUB1 : public Analysis {
#elif MODE==2
class ALICE_2017_I1512107_GRIDSUB2 : public Analysis {
#elif MODE==3
class ALICE_2017_I1512107_CHARGED : public Analysis {
#endif
public:
/// Constructor
#if MODE==0
ALICE_2017_I1512107_4MOMSUB() : Analysis("ALICE_2017_I1512107_4MOMSUB")
#elif MODE==1
ALICE_2017_I1512107_GRIDSUB1() : Analysis("ALICE_2017_I1512107_GRIDSUB1")
#elif MODE==2
ALICE_2017_I1512107_GRIDSUB2() : Analysis("ALICE_2017_I1512107_GRIDSUB2")
#elif MODE==3
ALICE_2017_I1512107_CHARGED() : Analysis("ALICE_2017_I1512107_CHARGED")
#endif
{
_jetR = 0.4; _ptmin = 50.; _etamax = 0.5; _etaloose = 0.9;
_nphibins = 120; _netabins = 160;
//_nphibins = 60; _netabins = 80;
_ptscalefac = 0.75; _massscalefac = 0.67;
}
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void init() {
_ptedges += 60.0, 80.0, 100.0, 120.0;
//FinalState fs(-_etaloose, _etaloose, 0.*GeV);
FinalState fs(-3., 3., 0.*GeV);
addProjection(fs,"FS");
ChargedFinalState cfs(-_etaloose, _etaloose, 0.*GeV);
addProjection(cfs,"CFS");
FastJets tfj(cfs, FastJets::ANTIKT, _jetR);
addProjection(tfj, "TrackJets");
FastJets ffj(fs, FastJets::ANTIKT, _jetR);
addProjection(ffj, "FullJets");
//_h_scalefac = bookHisto1D("scalefac", 50, 0., 1.);
_h_jetmass[0] = bookHisto1D(1, 1, 1);
_h_jetmass[1] = bookHisto1D(2, 1, 1);
_h_jetmass[2] = bookHisto1D(3, 1, 1);
_h_jetpt = bookHisto1D("jetpT", 40, 100., 400.);
_h_jetpt_resc = bookHisto1D("jetpT_resc", 40, 100., 400.);
_h_JetMass = bookHisto1D("JetMass", 50, 0., 100.);
_h_JetpT = bookHisto1D("JetpT", 30, 100, 400);
}
PseudoJets extractThermalMomenta(const Event & event) {
PseudoJets thermom;
foreach (const HepMC::GenParticle* p, particles(event.genEvent())) {
FourMomentum mom(p->momentum());
if (p->status() == 3) {
PseudoJet pjet(mom.px(),mom.py(),mom.pz(),mom.E());
thermom.push_back(pjet);
}
}
return thermom;
}
/// 4-momentum subtraction
FourMomentum SubtractJetMom(Jet jet, PseudoJets * thermom, double fac = 1.) {
if (thermom->empty()) return jet.momentum();
FourMomentum sub(0.,0.,0.,0.);
foreach (Particle part, jet.constituents()) {
if (part.E() < 1e-5 && part.E() > 1e-7) {
FourMomentum test(part.momentum());
test *= 10000.;
//cout<<"Found dummy :"<<Get4Mom(pj)<<endl;
double dRmin(10.);
PseudoJet matchthmom;
foreach (PseudoJet tm, *thermom) {
double dR = deltaR(test.eta(), test.phi(ZERO_2PI), Get4Mom(tm).eta(),Get4Mom(tm).phi(ZERO_2PI));
if (dR < dRmin) {
dRmin = dR;
matchthmom = tm;
}
}
if (dRmin < 1e-5) {
//cout<<"Found matching scattering centre :"<<Get4Mom(matchscatcen)<<endl;
sub += fac*Get4Mom(matchthmom);
}
else cout<<"Error: did not find scattering centre matching dummy.\n";
}
}
FourMomentum jetmom(jet.momentum());
//cout<<"Original momentum: "<<jetmom<<endl;
jetmom -= sub;
//cout<<"subtracted momentum: "<<jetmom<<endl;
//cout<<"Subtracted pt: "<<subpt<<" "<<pjet.pt()<<endl;
return jetmom;
}
/// initialise and fill grid for grid subtration version 2
vector<vector<FourMomentum> > fillGrid(const Event & event, std::set<std::pair<size_t,size_t> > & gridevent) {
// initialise grid
vector<vector<FourMomentum> > grid;
grid.resize(_nphibins);
for (size_t i = 0; i < _nphibins; ++i) {
grid[i].resize(_netabins);
for (size_t k = 0; k < _nphibins; ++k) {
grid[i][k] = FourMomentum(0., 0., 0., 0.);
}
}
// fill grid
foreach (const HepMC::GenParticle* p, particles(event.genEvent())) {
FourMomentum mom(p->momentum());
if (fabs(mom.eta()) < 4. && mom.E() > 1e-5) {
size_t phibin, etabin;
phibin = int(mom.phi(ZERO_2PI)*_nphibins/(2.*M_PI));
etabin = int((mom.eta() + 4.)*_netabins/8.);
if (phibin < 0 || phibin > _nphibins-1) {std::cout<<"Error: "<<mom.phi(ZERO_2PI)<<" --> "<<phibin<<std::endl; exit(1);}
if (etabin < 0 || etabin > _netabins-1) {std::cout<<"Error: "<<mom.eta()<<" --> "<<etabin<<std::endl; exit(1);}
/*if (p->status() == 3) {
grid[phibin][etabin] -= mom;
}
else {
grid[phibin][etabin] += mom;
}*/
double cellphi((phibin+0.5)*2*M_PI/_nphibins), celleta((etabin+0.5)*8./_netabins-4.);
double celltheta(2.*atan(exp(-celleta)));
if (p->status() == 3) {
double rho(grid[phibin][etabin].E() - mom.E());
grid[phibin][etabin].setE(rho);
grid[phibin][etabin].setPx(fabs(rho)*sin(celltheta)*cos(cellphi));
grid[phibin][etabin].setPy(fabs(rho)*sin(celltheta)*sin(cellphi));
grid[phibin][etabin].setPz(fabs(rho)*cos(celltheta));
}
else {
double rho(grid[phibin][etabin].E() + mom.E());
grid[phibin][etabin].setE(rho);
grid[phibin][etabin].setPx(rho*sin(celltheta)*cos(cellphi));
grid[phibin][etabin].setPy(rho*sin(celltheta)*sin(cellphi));
grid[phibin][etabin].setPz(rho*cos(celltheta));
}
gridevent.insert(std::pair<size_t,size_t>(phibin,etabin));
}
}
return grid;
}
/// determine fraction of pt carried by charged particles in jet
double getScaleFac(const Jet jet, double ptcut, const FinalState * fs = NULL) {
Particles parts;
if (fs) parts = fs->particles();
else parts = jet.constituents();
double scalefac;
double ntracks(0.), nparts(0.);
FourMomentum chmom(0.,0.,0.,0.), allmom(0.,0.,0.,0.);
foreach (Particle part, parts) {
if (deltaR(part,jet) < _jetR && part.momentum().E() > 1e-5) {
nparts += 1.;
allmom += part.momentum();
if (PID::charge(part.pdgId()) != 0 && part.pT() > ptcut) {
ntracks += 1.;
chmom += part.momentum();
}
}
}
scalefac = (nparts>0.?ntracks/nparts:1);
//scalefac = (allmom.pT()>0.?chmom.pT()/allmom.pT():1.);
return scalefac;
}
/// initialise and fill grid for grid subtration version 1
void fillJetGrid(Jet jet, PseudoJets * thermom, vector<vector<FourMomentum> > & grid, vector<vector<FourMomentum> > & thgrid, std::set<std::pair<size_t,size_t> > & gridevent) {
// fill grid
foreach (Particle part, jet.constituents()) {
FourMomentum mom(part.momentum());
size_t phibin, etabin;
if (mom.E() < 1e-5 && mom.E() > 1e-7) {
/*FourMomentum test(mom);
test *= 10000.;
//cout<<"Found dummy :"<<Get4Mom(mom)<<endl;
double dRmin(10.);
PseudoJet matchthmom;
foreach (PseudoJet tm, *thermom) {
double dR = deltaR(test.eta(), test.phi(ZERO_2PI), Get4Mom(tm).eta(),Get4Mom(tm).phi(ZERO_2PI));
if (dR < dRmin) {
dRmin = dR;
matchthmom = tm;
}
}
if (dRmin < 1e-5) {
mom = Get4Mom(matchthmom);
phibin = int(mom.phi(ZERO_2PI)*_nphibins/(2.*M_PI));
etabin = int((mom.eta() + 4.)*_netabins/8.);
if (phibin < 0 || phibin > _nphibins-1) {std::cout<<"Error: "<<mom.phi(ZERO_2PI)<<" --> "<<phibin<<std::endl; exit(1);}
if (etabin < 0 || etabin > _netabins-1) {std::cout<<"Error: "<<mom.eta()<<" --> "<<etabin<<std::endl; exit(1);}
thgrid[phibin][etabin] += mom;*/
/*double cellphi((phibin+0.5)*2*M_PI/_nphibins), celleta((etabin+0.5)*8./_netabins-4.);
double celltheta(2.*atan(exp(-celleta)));
double rho(thgrid[phibin][etabin].E() + mom.E());
thgrid[phibin][etabin].setE(rho);
thgrid[phibin][etabin].setPx(rho*sin(celltheta)*cos(cellphi));
thgrid[phibin][etabin].setPy(rho*sin(celltheta)*sin(cellphi));
thgrid[phibin][etabin].setPz(rho*cos(celltheta));*/
/* gridevent.insert(std::pair<size_t,size_t>(phibin,etabin));
}
else cout<<"Error: did not find scattering centre matching dummy.\n";*/
}
else {
phibin = int(mom.phi(ZERO_2PI)*_nphibins/(2.*M_PI));
etabin = int((mom.eta() + 4.)*_netabins/8.);
if (phibin < 0 || phibin > _nphibins-1) {std::cout<<"Error: "<<mom.phi(ZERO_2PI)<<" --> "<<phibin<<std::endl; exit(1);}
if (etabin < 0 || etabin > _netabins-1) {std::cout<<"Error: "<<mom.eta()<<" --> "<<etabin<<std::endl; exit(1);}
//grid[phibin][etabin] += mom;
double cellphi((phibin+0.5)*2*M_PI/_nphibins), celleta((etabin+0.5)*8./_netabins-4.);
double celltheta(2.*atan(exp(-celleta)));
double rho(grid[phibin][etabin].E() + mom.E());
grid[phibin][etabin].setE(rho);
grid[phibin][etabin].setPx(rho*sin(celltheta)*cos(cellphi));
grid[phibin][etabin].setPy(rho*sin(celltheta)*sin(cellphi));
grid[phibin][etabin].setPz(rho*cos(celltheta));
gridevent.insert(std::pair<size_t,size_t>(phibin,etabin));
}
}
foreach (PseudoJet pjet, *thermom) {
FourMomentum mom = Get4Mom(pjet);
double dR(deltaR(mom, jet));
if (dR < _jetR) {
size_t phibin, etabin;
phibin = int(mom.phi(ZERO_2PI)*_nphibins/(2.*M_PI));
etabin = int((mom.eta() + 4.)*_netabins/8.);
if (phibin < 0 || phibin > _nphibins-1) {std::cout<<"Error: "<<mom.phi(ZERO_2PI)<<" --> "<<phibin<<std::endl; exit(1);}
if (etabin < 0 || etabin > _netabins-1) {std::cout<<"Error: "<<mom.eta()<<" --> "<<etabin<<std::endl; exit(1);}
//thgrid[phibin][etabin] += mom;
double cellphi((phibin+0.5)*2*M_PI/_nphibins), celleta((etabin+0.5)*8./_netabins-4.);
double celltheta(2.*atan(exp(-celleta)));
double rho(thgrid[phibin][etabin].E() + mom.E());
thgrid[phibin][etabin].setE(rho);
thgrid[phibin][etabin].setPx(rho*sin(celltheta)*cos(cellphi));
thgrid[phibin][etabin].setPy(rho*sin(celltheta)*sin(cellphi));
thgrid[phibin][etabin].setPz(rho*cos(celltheta));
gridevent.insert(std::pair<size_t,size_t>(phibin,etabin));
}
}
return;
}
/// 4-momentum adapter
FourMomentum Get4Mom(PseudoJet pjet){
FourMomentum mom(pjet.E(), pjet.px(), pjet.py(), pjet.pz());
return mom;
}
/// Perform the per-event analysis
void analyze(const Event& event) {
const double weight = event.weight();
const double cent = (event.genEvent()->heavy_ion()?event.genEvent()->heavy_ion()->impact_parameter():-1.);
if (cent > 0.1) vetoEvent;
#if MODE==0
// cout<<"Karo: event number "<<event.genEvent()->event_number()<<endl;
//Cut cuts = Cuts::abseta < _etaloose && Cuts::pT > _ptmin*GeV;
Cut cuts = Cuts::abseta < 2.0 && Cuts::pT > _ptmin*GeV;
const Jets jets = applyProjection<FastJets>(event, "FullJets").jetsByPt(cuts);
PseudoJets thermom = extractThermalMomenta(event);
foreach (Jet jet, jets) {
//double scalefac;
//scalefac = getScaleFac(jet, 0.);
//_h_scalefac->fill(scalefac, weight);
FourMomentum jetsub = SubtractJetMom(jet, &thermom, 1.);
// if (jet.momentum().pT() > 100.) {
// cout<<"Karo: "<<jet.momentum().pT()<<" -> "<<jetsub.pT()<<endl;
// }
if (abs(jetsub.eta()) < _etamax) {
double mass2(jetsub.mass2());
double mass(mass2>0.?sqrt(mass2):-1.);
//mass *= scalefac;
mass *= _massscalefac;
//double jetpt(scalefac*jetsub.pT());
double jetpt(jetsub.pT());
_h_jetpt->fill(jetpt,weight);
jetpt *= _ptscalefac;
_h_jetpt_resc->fill(jetpt,weight);
for (size_t i=0; i<3; ++i) {
if (jetpt > _ptedges[i] && jetpt <= _ptedges[i+1]) _h_jetmass[i]->fill(mass,weight);
}
}
if (jetsub.pT() > 100.*GeV) {
// cout<<"Karo: "<<jetsub.pT()<<endl;
double mass2(jetsub.mass2());
double mass(mass2>0.?sqrt(mass2):-1.);
_h_JetMass->fill(mass,weight);
_h_JetpT->fill(jetsub.pT(),weight);
}
}
#elif MODE==1
//Cut cuts = Cuts::abseta < _etaloose && Cuts::pT > _ptmin*GeV;
Cut cuts = Cuts::abseta < 2.0 && Cuts::pT > _ptmin*GeV;
const Jets jets = applyProjection<FastJets>(event, "FullJets").jetsByPt(cuts);
PseudoJets thermom = extractThermalMomenta(event);
// initialise grid
vector<vector<FourMomentum> > grid, thgrid;
grid.resize(_nphibins);
thgrid.resize(_nphibins);
for (size_t i = 0; i < _nphibins; ++i) {
grid[i].resize(_netabins);
thgrid[i].resize(_netabins);
for (size_t k = 0; k < _nphibins; ++k) {
grid[i][k] = FourMomentum(0., 0., 0., 0.);
thgrid[i][k] = FourMomentum(0., 0., 0., 0.);
}
}
std::set<std::pair<size_t,size_t> > gridevent;
foreach (Jet jet, jets) {
fillJetGrid(jet, &thermom, grid, thgrid, gridevent);
}
double negenergy(0.);
PseudoJets pevent;
for (set<pair<size_t,size_t> >::iterator it=gridevent.begin(); it!=gridevent.end(); ++it) {
size_t i,k;
i = it->first;
k = it->second;
FourMomentum diff(grid[i][k] - thgrid[i][k]);
if (grid[i][k].pT() > thgrid[i][k].pT() && diff.E() > 0.) {
PseudoJet pjet(diff.px(),diff.py(),diff.pz(),diff.E());
pevent.push_back(pjet);
}
else negenergy += diff.E();
}
JetDefinition jet_def(antikt_algorithm, _jetR);
ClusterSequence cs(pevent, jet_def);
PseudoJets pjets = sorted_by_pt(cs.inclusive_jets(_ptmin));
//cout<<"size comparison: "<<jets.size()<<" vs. "<<pjets.size()<<endl;
foreach (PseudoJet pjet, pjets) {
FourMomentum jetsub = Get4Mom(pjet);
if (abs(jetsub.eta()) < _etamax) {
double mass2(jetsub.mass2());
double mass(mass2>0.?sqrt(mass2):-1.);
mass *= _massscalefac;
double jetpt(jetsub.pT());
_h_jetpt->fill(jetpt,weight);
jetpt *= _ptscalefac;
_h_jetpt_resc->fill(jetpt,weight);
for (size_t i=0; i<3; ++i) {
if (jetpt > _ptedges[i] && jetpt <= _ptedges[i+1]){
_h_jetmass[i]->fill(mass,weight);
}
}
}
if (jetsub.pT() > 100.*GeV) {
double mass2(jetsub.mass2());
double mass(mass2>0.?sqrt(mass2):-1.);
_h_JetMass->fill(mass,weight);
_h_JetpT->fill(jetsub.pT(),weight);
}
}
#elif MODE==2
const FinalState parts = applyProjection<FinalState>(event, "FS");
PseudoJets pevent, pjets;
std::set<std::pair<size_t,size_t> > gridevent;
vector<vector<FourMomentum> > grid = fillGrid(event, gridevent);
for (set<pair<size_t,size_t> >::iterator it=gridevent.begin(); it!=gridevent.end(); ++it) {
size_t i,k;
i = it->first;
k = it->second;
if (grid[i][k].E() > 0. ) {
PseudoJet part(grid[i][k].px(),grid[i][k].py(),grid[i][k].pz(),grid[i][k].E());
pevent.push_back(part);
}
}
JetDefinition jet_def(antikt_algorithm, _jetR);
ClusterSequence cs(pevent, jet_def);
pjets = sorted_by_pt(cs.inclusive_jets(_ptmin));
foreach (PseudoJet pjet, pjets) {
FourMomentum jetsub = Get4Mom(pjet);
if (abs(jetsub.eta()) < _etamax) {
//double scalefac;
//scalefac = getScaleFac(Jet(pjet),0.,&parts);
//_h_scalefac->fill(scalefac,weight);
double mass2(jetsub.mass2());
double mass(mass2>0.?sqrt(mass2):-1.);
//mass *= scalefac;
mass *= _massscalefac;
//double jetpt(scalefac*jetsub.pT());
double jetpt(jetsub.pT());
_h_jetpt->fill(jetpt,weight);
jetpt *= _ptscalefac;
_h_jetpt_resc->fill(jetpt,weight);
for (size_t i=0; i<3; ++i) {
if (jetpt > _ptedges[i] && jetpt <= _ptedges[i+1]){
_h_jetmass[i]->fill(mass,weight);
}
}
}
if (jetsub.pT() > 100.*GeV) {
double mass2(jetsub.mass2());
double mass(mass2>0.?sqrt(mass2):-1.);
_h_JetMass->fill(mass,weight);
_h_JetpT->fill(jetsub.pT(),weight);
}
}
#elif MODE==3
const FinalState parts = applyProjection<FinalState>(event, "FS");
Cut cuts = Cuts::abseta < _etaloose && Cuts::pT > _ptmin*GeV;
const Jets jets = applyProjection<FastJets>(event, "TrackJets").jetsByPt(cuts);
foreach (Jet jet, jets) {
//double scalefac;
//scalefac = getScaleFac(jet, 0., &parts);
//_h_scalefac->fill(scalefac, weight);
PseudoJets thermom = extractThermalMomenta(event);
//FourMomentum jetsub = SubtractJetMom(jet, &thermom, scalefac);
FourMomentum jetsub = SubtractJetMom(jet, &thermom, _ptscalefac);
if (abs(jetsub.eta()) < _etamax) {
double mass2(jetsub.mass2());
double mass(mass2>0.?sqrt(mass2):-1.);
double jetpt(jetsub.pT());
_h_jetpt_resc->fill(jetpt,weight);
for (size_t i=0; i<3; ++i) {
if (jetpt > _ptedges[i] && jetpt <= _ptedges[i+1]) _h_jetmass[i]->fill(mass,weight);
}
}
}
#endif
}
/// Normalise histograms etc., after the run
void finalize() {
//normalize(_h_scalefac);
for (size_t i=0; i<3; ++i) {
normalize(_h_jetmass[i]);
}
scale(_h_jetpt, crossSection()/sumOfWeights());
scale(_h_jetpt_resc, crossSection()/sumOfWeights());
normalize(_h_JetMass);
scale(_h_JetpT, crossSection()/picobarn/sumOfWeights());
}
private:
double _jetR, _ptmin, _etamax, _etaloose;
double _massscalefac, _ptscalefac;
size_t _netabins, _nphibins;
vector<double> _ptedges;
//Histo1DPtr _h_scalefac;
Histo1DPtr _h_jetmass[3];
Histo1DPtr _h_jetpt, _h_jetpt_resc;
Histo1DPtr _h_JetMass, _h_JetpT;
};
// The hook for the plugin system
#if MODE==0
DECLARE_RIVET_PLUGIN(ALICE_2017_I1512107_4MOMSUB);
#elif MODE==1
DECLARE_RIVET_PLUGIN(ALICE_2017_I1512107_GRIDSUB1);
#elif MODE==2
DECLARE_RIVET_PLUGIN(ALICE_2017_I1512107_GRIDSUB2);
#elif MODE==3
DECLARE_RIVET_PLUGIN(ALICE_2017_I1512107_CHARGED);
#endif
}
Index: trunk/code/ALICE_2012_I1127497.plot
===================================================================
--- trunk/code/ALICE_2012_I1127497.plot (revision 461)
+++ trunk/code/ALICE_2012_I1127497.plot (revision 462)
@@ -1,20 +1,126 @@
+# BEGIN PLOT /ALICE_2012_I1127497/d01-x01-y01
+Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 00-05
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d02-x01-y01
+Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 05-10
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d03-x01-y01
+Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 10-20
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d04-x01-y01
+Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 20-30
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d05-x01-y01
+Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 30-40
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d06-x01-y01
+Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 40-50
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d07-x01-y01
+Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 50-60
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d08-x01-y01
+Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 60-70
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d09-x01-y01
+Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 70-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$
+# END PLOT
+
# BEGIN PLOT /ALICE_2012_I1127497/d10-x01-y01
Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 00-10
XLabel=$p_\perp\ [\mathrm{GeV}]$
YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$
LogY=1
FullRange=1
NormalizeToIntegral=1
-#DrawOnly=/home/jewel/trunk/ALICE_2012_I1127497.aida/REF/ALICE_2012_I1127497/d10-x01-y01 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-10_wo_1000-1099.aida/ALICE_2012_I1127497/d10-x01-y01 /home/jewel/trunk/data/jewel-2.0.0_vac_fixed2_5-350.aida/ALICE_2012_I1127497/d10-x01-y01
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d16-x01-y01
+Title=charged hadrons ($|\eta|<1$), 00-05/pp
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{AA}$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d17-x01-y01
+Title=charged hadrons ($|\eta|<1$), 05-10/pp
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{AA}$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d18-x01-y01
+Title=charged hadrons ($|\eta|<1$), 10-20/pp
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{AA}$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d19-x01-y01
+Title=charged hadrons ($|\eta|<1$), 20-30/pp
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{AA}$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d20-x01-y01
+Title=charged hadrons ($|\eta|<1$), 30-40/pp
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{AA}$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d21-x01-y01
+Title=charged hadrons ($|\eta|<1$), 40-50/pp
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{AA}$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d22-x01-y01
+Title=charged hadrons ($|\eta|<1$), 50-60/pp
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{AA}$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d23-x01-y01
+Title=charged hadrons ($|\eta|<1$), 60-70/pp
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{AA}$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2012_I1127497/d24-x01-y01
+Title=charged hadrons ($|\eta|<1$), 70-80/pp
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{AA}$
# END PLOT
# BEGIN PLOT /ALICE_2012_I1127497/d25-x01-y01
Title=charged hadrons ($|\eta|<1$), 00-10/pp
XLabel=$p_\perp\ [\mathrm{GeV}]$
YLabel=$R_{AA}$
LogY=0
FullRange=1
YMax=1.
-#DrawOnly=/home/jewel/trunk/ALICE_2012_I1127497.aida/REF/ALICE_2012_I1127497/d25-x01-y01 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-10_yres_wo_1000-1099.aida/ALICE_2012_I1127497/d25-x01-y01
# END PLOT
Index: trunk/code/ATLAS_2012_I1126965.plot
===================================================================
--- trunk/code/ATLAS_2012_I1126965.plot (revision 461)
+++ trunk/code/ATLAS_2012_I1126965.plot (revision 462)
@@ -1,118 +1,154 @@
-# BEGIN PLOT /ATLAS_2012_I1126965/d0*
+# BEGIN PLOT /ATLAS_2012_I1126965/*
LogY=0
FullRange=1
YMax=1.
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y01
Title=R=0.2, 00-10/60-80
XLabel=$p_\perp\ [\mathrm{GeV}]$
YLabel=$R_{CP}$
-LogY=0
-FullRange=1
-YMax=1.
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y02
Title=R=0.3, 00-10/60-80
XLabel=$p_\perp\ [\mathrm{GeV}]$
YLabel=$R_{CP}$
-LogY=0
-FullRange=1
-YMax=1.
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y03
Title=R=0.4, 00-10/60-80
XLabel=$p_\perp\ [\mathrm{GeV}]$
YLabel=$R_{CP}$
-LogY=0
-FullRange=1
-YMax=1.
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y04
Title=R=0.5, 00-10/60-80
XLabel=$p_\perp\ [\mathrm{GeV}]$
YLabel=$R_{CP}$
-LogY=0
-FullRange=1
-YMax=1.
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/jetpt0*
#DrawOnly=
# END PLOT
-# BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y02
-#DrawOnly=
-# END PLOT
-
-# BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y03
-#DrawOnly=
-# END PLOT
-
-# BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y04
-#DrawOnly=
+# BEGIN PLOT /ATLAS_2012_I1126965/d03-x01-y01
+Title=R=0.2, 10-20/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d03-x01-y02
-#DrawOnly=
+Title=R=0.3, 10-20/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d03-x01-y03
-#DrawOnly=
+Title=R=0.4, 10-20/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d03-x01-y04
-#DrawOnly=
+Title=R=0.5, 10-20/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1126965/d04-x01-y01
+Title=R=0.2, 20-30/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d04-x01-y02
-#DrawOnly=
+Title=R=0.3, 20-30/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d04-x01-y03
-#DrawOnly=
+Title=R=0.4, 20-30/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d04-x01-y04
-#DrawOnly=
+Title=R=0.5, 20-30/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1126965/d05-x01-y01
+Title=R=0.2, 30-40/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d05-x01-y02
-#DrawOnly=
+Title=R=0.3, 30-40/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d05-x01-y03
-#DrawOnly=
+Title=R=0.4, 30-40/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d05-x01-y04
-#DrawOnly=
+Title=R=0.5, 30-40/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1126965/d06-x01-y01
+Title=R=0.2, 40-50/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d06-x01-y02
-#DrawOnly=
+Title=R=0.3, 40-50/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d06-x01-y03
-#DrawOnly=
+Title=R=0.4, 40-50/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d06-x01-y04
-#DrawOnly=
+Title=R=0.5, 40-50/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1126965/d07-x01-y01
+Title=R=0.2, 50-60/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d07-x01-y02
-#DrawOnly=
+Title=R=0.3, 50-60/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d07-x01-y03
-#DrawOnly=
+Title=R=0.4, 50-60/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
# BEGIN PLOT /ATLAS_2012_I1126965/d07-x01-y04
-#DrawOnly=
+Title=R=0.5, 50-60/60-80
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_{CP}$
# END PLOT
Index: trunk/code/ATLAS_2013_I1240088.cc
===================================================================
--- trunk/code/ATLAS_2013_I1240088.cc (revision 461)
+++ trunk/code/ATLAS_2013_I1240088.cc (revision 462)
@@ -1,124 +1,125 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/FastJets.hh"
#include <boost/lexical_cast.hpp>
/// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder...
namespace Rivet {
class ATLAS_2013_I1240088 : public Analysis {
public:
/// @name Constructors etc.
//@{
/// Constructor
ATLAS_2013_I1240088()
: Analysis("ATLAS_2013_I1240088")
{
- normalise = false;
}
//@}
public:
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void init() {
FinalState fs(-4.0, 4.0, 0.*GeV);
addProjection(fs, "FS");
FastJets fj02(fs, FastJets::ANTIKT, 0.2);
fj02.useInvisibles();
addProjection(fj02, "Jets02");
- _centedges += 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6;
+ //_centedges += 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6;
+ _centedges += 0.05, 0.1;
_ptedges += 45.0, 60.0, 80.0, 110.0, 160.0;
_ncentbins = _centedges.size()-1;
_nptbins = _ptedges.size()-1;
for (size_t i = 0; i < _ncentbins; ++i) {
for (size_t j = 0; j < _nptbins; ++j) {
+ normalise[i][j] = false;
_histos[i][j] = bookHisto1D(i+1, 2, j+1);
}
}
}
/// Perform the per-event analysis
void analyze(const Event& event) {
const double weight = event.weight();
const double cent = (event.genEvent()->heavy_ion()?event.genEvent()->heavy_ion()->impact_parameter():-1.);
Cut cuts = Cuts::abseta < 2.1 && Cuts::pT > _ptedges.front()*GeV && Cuts::pT < _ptedges.back()*GeV;
Jets alljets = applyProjection<FastJets>(event, "Jets02").jetsByPt(cuts);
for (size_t i = 0; i < _ncentbins; ++i) {
if (cent < 0. || (cent >= _centedges[i] && cent < _centedges[i+1])) {
for (size_t j = 0; j < _nptbins; ++j) {
foreach (Jet jet, alljets) {
double pt(jet.momentum().pT());
if (pt > _ptedges[j]*GeV && pt < _ptedges[j+1]*GeV) {
- normalise = true;
+ normalise[i][j] = true;
double phi(jet.phi());
if (phi > M_PI/2. && phi < M_PI) phi = M_PI - phi;
if (phi > M_PI && phi < 3.*M_PI/2.) phi = phi - M_PI;
if (phi > 3.*M_PI/2.) phi = 2.*M_PI - phi;
_histos[i][j]->fill(phi,weight);
}
}
}
}
}
}
/// Normalise histograms etc., after the run
void finalize() {
- if (!normalise) return;
+ //if (!normalise) return;
for (size_t i = 0; i < _ncentbins; ++i) {
for (size_t j = 0; j < _nptbins; ++j) {
- normalize(_histos[i][j]);
+ if (normalise[i][j]) normalize(_histos[i][j]);
}
}
}
//@}
private:
// Data members like post-cuts event weight counters go here
vector<double> _ptedges;
vector<double> _centedges;
size_t _ncentbins, _nptbins;
- bool normalise;
+ bool normalise[1][4];
private:
/// @name Histograms
//@{
- Histo1DPtr _histos[6][4];
+ Histo1DPtr _histos[1][4];
//@}
};
// The hook for the plugin system
DECLARE_RIVET_PLUGIN(ATLAS_2013_I1240088);
}
Index: trunk/code/ALICE_2015_I1343112.plot
===================================================================
--- trunk/code/ALICE_2015_I1343112.plot (revision 461)
+++ trunk/code/ALICE_2015_I1343112.plot (revision 462)
@@ -1,8 +1,67 @@
# BEGIN PLOT /ALICE_2015_I1343112/d01-x01-y01
-#Title=[Uncomment and insert title for histogram d01-x01-y01 here]
-#XLabel=[Uncomment and insert x-axis label for histogram d01-x01-y01 here]
-#YLabel=[Uncomment and insert y-axis label for histogram d01-x01-y01 here]
-# + any additional plot settings you might like, see make-plots documentation
+Title=p+p, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$1/N_\text{coll}\, 1/N_\text{evt} \mathrm{d}^2N/\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2015_I1343112/d01-x01-y02
+Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$1/N_\text{coll}\, 1/N_\text{evt} \mathrm{d}^2N/\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2015_I1343112/d01-x01-y03
+Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$1/N_\text{coll}\, 1/N_\text{evt} \mathrm{d}^2N/\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ALICE_2015_I1343112/d02-x01-y01
+Title=p+p, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$N_\text{jets}^\text{5 GeV bias}/N_\text{jets}^\text{0 GeV bias}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2015_I1343112/d03-x01-y01
+Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$N_\text{jets}^\text{0 GeV bias}/N_\text{jets}^\text{5 GeV bias}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2015_I1343112/d03-x01-y02
+Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 3\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$N_\text{jets}^\text{3 GeV bias}/N_\text{jets}^\text{5 GeV bias}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2015_I1343112/d03-x01-y03
+Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 7\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$N_\text{jets}^\text{7 GeV bias}/N_\text{jets}^\text{5 GeV bias}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2015_I1343112/d03-x01-y04
+Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$N_\text{jets}^\text{10 GeV bias}/N_\text{jets}^\text{5 GeV bias}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2015_I1343112/d04-x01-y01
+Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{AA}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ALICE_2015_I1343112/d04-x01-y02
+Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$
+XLabel=$p_\perp\ [\mathrm{GeV}]$
+YLabel=$R_\text{AA}$
+LogY=0
# END PLOT
-# ... add more histograms as you need them ...
Index: trunk/code/ALICE_2017_I1512107_4MOMSUB.plot
===================================================================
--- trunk/code/ALICE_2017_I1512107_4MOMSUB.plot (revision 461)
+++ trunk/code/ALICE_2017_I1512107_4MOMSUB.plot (revision 462)
@@ -1,19 +1,28 @@
# BEGIN PLOT /ALICE_2017_I1512107_4MOMSUB/chfrac_rat
LogY=0
# END PLOT
# BEGIN PLOT /ALICE_2017_I1512107_4MOMSUB/jetmass
LogY=0
# END PLOT
# BEGIN PLOT /ALICE_2017_I1512107_4MOMSUB/d01-x01-y01
+XLabel=$M_\text{ch jet}\ [\text{GeV}^2]$
+YLabel=$1/N_\text{jets}\ \mathrm{d}N/\mathrm{d}M_\text{ch jet}\ [\mathrm{GeV}^{-2}]$
+Title=$60\,\mathrm{GeV} < p_{\perp,\,\text{ch jet}} < 80\,\mathrm{GeV}$ - 0-10 Pb+Pb
LogY=0
# END PLOT
# BEGIN PLOT /ALICE_2017_I1512107_4MOMSUB/d02-x01-y01
+XLabel=$M_\text{ch jet}\ [\text{GeV}^2]$
+YLabel=$1/N_\text{jets}\ \mathrm{d}N/\mathrm{d}M_\text{ch jet}\ [\mathrm{GeV}^{-2}]$
+Title=$80\,\mathrm{GeV} < p_{\perp,\,\text{ch jet}} < 100\,\mathrm{GeV}$ - 0-10 Pb+Pb
LogY=0
# END PLOT
# BEGIN PLOT /ALICE_2017_I1512107_4MOMSUB/d03-x01-y01
+XLabel=$M_\text{ch jet}\ [\text{GeV}^2]$
+YLabel=$1/N_\text{jets}\ \mathrm{d}N/\mathrm{d}M_\text{ch jet}\ [\mathrm{GeV}^{-2}]$
+Title=$100\,\mathrm{GeV} < p_{\perp,\,\text{ch jet}} < 120\,\mathrm{GeV}$ - 0-10 Pb+Pb
LogY=0
# END PLOT
Index: trunk/code/ATLAS_CONF_2014_028.plot
===================================================================
--- trunk/code/ATLAS_CONF_2014_028.plot (revision 461)
+++ trunk/code/ATLAS_CONF_2014_028.plot (revision 462)
@@ -1,8 +1,546 @@
# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x01-y01
-#Title=[Uncomment and insert title for histogram d01-x01-y01 here]
-#XLabel=[Uncomment and insert x-axis label for histogram d01-x01-y01 here]
-#YLabel=[Uncomment and insert y-axis label for histogram d01-x01-y01 here]
-# + any additional plot settings you might like, see make-plots documentation
+Title=$R=0.2$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
# END PLOT
-# ... add more histograms as you need them ...
+# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x01-y02
+Title=$R=0.2$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x01-y03
+Title=$R=0.2$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x01-y04
+Title=$R=0.2$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x02-y01
+Title=$R=0.2$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x02-y02
+Title=$R=0.2$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x02-y03
+Title=$R=0.2$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x02-y04
+Title=$R=0.2$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x03-y01
+Title=$R=0.2$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x03-y02
+Title=$R=0.2$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x03-y03
+Title=$R=0.2$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x03-y04
+Title=$R=0.2$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x01-y01
+Title=$R=0.3$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x01-y02
+Title=$R=0.3$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x01-y03
+Title=$R=0.3$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x01-y04
+Title=$R=0.3$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x02-y01
+Title=$R=0.3$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x02-y02
+Title=$R=0.3$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x02-y03
+Title=$R=0.3$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x02-y04
+Title=$R=0.3$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x03-y01
+Title=$R=0.3$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x03-y02
+Title=$R=0.3$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x03-y03
+Title=$R=0.3$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x03-y04
+Title=$R=0.3$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x01-y01
+Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x01-y02
+Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x01-y03
+Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x01-y04
+Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x02-y01
+Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x02-y02
+Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x02-y03
+Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x02-y04
+Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x03-y01
+Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x03-y02
+Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x03-y03
+Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x03-y04
+Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$R_{\Delta R}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x01-y01
+Title=$R=0.2$, $E_\perp^\text{test} > 80\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x01-y02
+Title=$R=0.2$, $E_\perp^\text{test} > 80\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x01-y03
+Title=$R=0.2$, $E_\perp^\text{test} > 80\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x01-y04
+Title=$R=0.2$, $E_\perp^\text{test} > 80\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x02-y01
+Title=$R=0.2$, $E_\perp^\text{test} > 90\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x02-y02
+Title=$R=0.2$, $E_\perp^\text{test} > 90\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x02-y03
+Title=$R=0.2$, $E_\perp^\text{test} > 90\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x02-y04
+Title=$R=0.2$, $E_\perp^\text{test} > 90\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x03-y01
+Title=$R=0.2$, $E_\perp^\text{test} > 110\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x03-y02
+Title=$R=0.2$, $E_\perp^\text{test} > 110\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x03-y03
+Title=$R=0.2$, $E_\perp^\text{test} > 110\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x03-y04
+Title=$R=0.2$, $E_\perp^\text{test} > 110\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x01-y01
+Title=$R=0.3$, $E_\perp^\text{test} > 80\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x01-y02
+Title=$R=0.3$, $E_\perp^\text{test} > 80\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x01-y03
+Title=$R=0.3$, $E_\perp^\text{test} > 80\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x01-y04
+Title=$R=0.3$, $E_\perp^\text{test} > 80\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x02-y01
+Title=$R=0.3$, $E_\perp^\text{test} > 90\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x02-y02
+Title=$R=0.3$, $E_\perp^\text{test} > 90\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x02-y03
+Title=$R=0.3$, $E_\perp^\text{test} > 90\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x02-y04
+Title=$R=0.3$, $E_\perp^\text{test} > 90\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x03-y01
+Title=$R=0.3$, $E_\perp^\text{test} > 110\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x03-y02
+Title=$R=0.3$, $E_\perp^\text{test} > 110\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x03-y03
+Title=$R=0.3$, $E_\perp^\text{test} > 110\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x03-y04
+Title=$R=0.3$, $E_\perp^\text{test} > 110\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x01-y01
+Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x01-y02
+Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x01-y03
+Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x01-y04
+Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x02-y01
+Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x02-y02
+Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x02-y03
+Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x02-y04
+Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x03-y01
+Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 0-10
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x03-y02
+Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 10-20
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x03-y03
+Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 20-40
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x03-y04
+Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x01-y01
+Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 0-10 / 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x01-y02
+Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 10-20 / 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x01-y03
+Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 20-40 / 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x02-y01
+Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 0-10 / 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x02-y02
+Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 10-20 / 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x02-y03
+Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 20-40 / 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x03-y01
+Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 0-10 / 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x03-y02
+Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 10-20 / 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x03-y03
+Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 20-40 / 40-80
+XLabel=$E_\perp^\text{test}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x01-y01
+Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 0-10 / 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x01-y02
+Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 10-20 / 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x01-y03
+Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 20-40 / 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x02-y01
+Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 0-10 / 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x02-y02
+Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 10-20 / 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x02-y03
+Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 20-40 / 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x03-y01
+Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 0-10 / 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x03-y02
+Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 10-20 / 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x03-y03
+Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 20-40 / 40-80
+XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$
+YLabel=$\rho_{R_{\Delta R}}$
+# END PLOT
+
+## BEGIN PLOT /ATLAS_CONF_2014_028/d9*
+#Title=dummy histogram
+## END PLOT
+
+
+
Index: trunk/code/CMS_2012_I1090064.plot
===================================================================
--- trunk/code/CMS_2012_I1090064.plot (revision 461)
+++ trunk/code/CMS_2012_I1090064.plot (revision 462)
@@ -1,155 +1,215 @@
-# BEGIN PLOT /CMS_2012_I1090064/d01-x01-y0*
-#NormalizeToIntegral=1
-Scale=0.10472
-# END PLOT
-
-# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y0*
-#NormalizeToIntegral=1
-Scale=0.06
-# END PLOT
-
-# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y0*
-#NormalizeToIntegral=1
-Scale=0.06
+# BEGIN PLOT /CMS_2012_I1090064/*
+FullRange=1
# END PLOT
-# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y0*
-#NormalizeToIntegral=1
-Scale=0.1
+# BEGIN PLOT /CMS_2012_I1090064/d01-x01-y0*
+Scale=0.10472
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d01-x01-y01
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d01-x01-y01 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d01-x01-y01
+Title=PbPb 0-20, $120\ \mathrm{GeV} < p_{\perp,1} < 150\ \mathrm{GeV}$
+XLabel=$\Delta \phi$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d01-x01-y02
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d01-x01-y02 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d01-x01-y02
+Title=PbPb 0-20, $150\ \mathrm{GeV} < p_{\perp,1} < 180\ \mathrm{GeV}$
+XLabel=$\Delta \phi$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d01-x01-y03
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d01-x01-y03 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d01-x01-y03
+Title=PbPb 0-20, $180\ \mathrm{GeV} < p_{\perp,1} < 220\ \mathrm{GeV}$
+XLabel=$\Delta \phi$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d01-x01-y04
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d01-x01-y04 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d01-x01-y04
+Title=PbPb 0-20, $220\ \mathrm{GeV} < p_{\perp,1} < 260\ \mathrm{GeV}$
+XLabel=$\Delta \phi$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d01-x01-y05
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d01-x01-y05 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d01-x01-y05
-# END PLOT
+#Title=PbPb 0-20, $260\ \mathrm{GeV} < p_{\perp,1} < 300\ \mathrm{GeV}$
+XLabel=$\Delta \phi$
+YLabel=event fraction
+ END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d01-x01-y06
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d01-x01-y06 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d01-x01-y06
+Title=PbPb 0-20, $300\ \mathrm{GeV} < p_{\perp,1} < 500\ \mathrm{GeV}$
+XLabel=$\Delta \phi$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d02-x01-y01
LogY=0
XMax=220.
YMin=0.7
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d02-x01-y01 /home/jewel/trunk/data/jewel-2.0.0_vac_fixed_5-350.aida/CMS_2012_I1090064/d02-x01-y01
+Title=di-jet fraction in pp
+XLabel=$p_{\perp,1}$ [GeV]
+YLabel=$N_\mathrm{dijet}/N_\mathrm{leading jet}$
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d02-x01-y02
LogY=0
XMax=370.
YMin=0.7
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d02-x01-y02 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d02-x01-y02
+Title=di-jet fraction in PbPb 0-20
+XLabel=$p_{\perp,1}$ [GeV]
+YLabel=$N_\mathrm{dijet}/N_\mathrm{leading jet}$
# END PLOT
-# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y01
+# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y0*
+Scale=0.06
LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d03-x01-y01 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d03-x01-y01
+# END PLOT
+
+# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y01
+Title=di-jet asymmetry in PbPb 0-20, $120\ \mathrm{GeV} < p_{\perp,1} < 150\ \mathrm{GeV}$
+XLabel=$A_J$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y02
-LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d03-x01-y02 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d03-x01-y02
+Title=di-jet asymmetry in PbPb 0-20, $150\ \mathrm{GeV} < p_{\perp,1} < 180\ \mathrm{GeV}$
+XLabel=$A_J$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y03
-LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d03-x01-y03 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d03-x01-y03
+Title=di-jet asymmetry in PbPb 0-20, $180\ \mathrm{GeV} < p_{\perp,1} < 220\ \mathrm{GeV}$
+XLabel=$A_J$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y04
-LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d03-x01-y04 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d03-x01-y04
+Title=di-jet asymmetry in PbPb 0-20, $220\ \mathrm{GeV} < p_{\perp,1} < 260\ \mathrm{GeV}$
+XLabel=$A_J$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y05
-LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d03-x01-y05 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d03-x01-y05
+Title=di-jet asymmetry in PbPb 0-20, $260\ \mathrm{GeV} < p_{\perp,1} < 300\ \mathrm{GeV}$
+XLabel=$A_J$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y06
-LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d03-x01-y06 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d03-x01-y06
+Title=di-jet asymmetry in PbPb 0-20, $300\ \mathrm{GeV} < p_{\perp,1} < 500\ \mathrm{GeV}$
+XLabel=$A_J$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y0*
+Scale=0.06
LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d04-x01-y01 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-10_wo_1000-1199.aida/CMS_2012_I1090064/d04-x01-y01
# END PLOT
-# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y01
+# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y01
+Title=di-jet asymmetry in PbPb 0-10, $p_{\perp,1} > 120\ \mathrm{GeV}$
+XLabel=$A_J$
+YLabel=event fraction
+# END PLOT
+
+# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y02
+Title=di-jet asymmetry in PbPb 10-20, $p_{\perp,1} > 120\ \mathrm{GeV}$
+XLabel=$A_J$
+YLabel=event fraction
+# END PLOT
+
+# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y03
+Title=di-jet asymmetry in PbPb 20-30, $p_{\perp,1} > 120\ \mathrm{GeV}$
+XLabel=$A_J$
+YLabel=event fraction
+# END PLOT
+
+# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y04
+Title=di-jet asymmetry in PbPb 30-50, $p_{\perp,1} > 120\ \mathrm{GeV}$
+XLabel=$A_J$
+YLabel=event fraction
+# END PLOT
+
+# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y05
+Title=di-jet asymmetry in PbPb 50-70, $p_{\perp,1} > 120\ \mathrm{GeV}$
+XLabel=$A_J$
+YLabel=event fraction
+# END PLOT
+
+# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y06
+Title=di-jet asymmetry in PbPb 70-100, $p_{\perp,1} > 120\ \mathrm{GeV}$
+XLabel=$A_J$
+YLabel=event fraction
+# END PLOT
+
+# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y0*
+Scale=0.1
LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d05-x01-y01 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d05-x01-y01
+# END PLOT
+
+# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y01
+Title=di-jet $p_\perp$ ratio in PbPb 0-20, $120\ \mathrm{GeV} < p_{\perp,1} < 150\ \mathrm{GeV}$
+XLabel=$p_{\perp,2}/p_{\perp,1}$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y02
-LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d05-x01-y02 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d05-x01-y02
+Title=di-jet $p_\perp$ ratio in PbPb 0-20, $150\ \mathrm{GeV} < p_{\perp,1} < 180\ \mathrm{GeV}$
+XLabel=$p_{\perp,2}/p_{\perp,1}$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y03
-LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d05-x01-y03 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d05-x01-y03
+Title=di-jet $p_\perp$ ratio in PbPb 0-20, $180\ \mathrm{GeV} < p_{\perp,1} < 220\ \mathrm{GeV}$
+XLabel=$p_{\perp,2}/p_{\perp,1}$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y04
-LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d05-x01-y04 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d05-x01-y04
+Title=di-jet $p_\perp$ ratio in PbPb 0-20, $220\ \mathrm{GeV} < p_{\perp,1} < 260\ \mathrm{GeV}$
+XLabel=$p_{\perp,2}/p_{\perp,1}$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y05
-LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d05-x01-y05 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d05-x01-y05
+Title=di-jet $p_\perp$ ratio in PbPb 0-20, $260\ \mathrm{GeV} < p_{\perp,1} < 300\ \mathrm{GeV}$
+XLabel=$p_{\perp,2}/p_{\perp,1}$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y06
-LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d05-x01-y06 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d05-x01-y06
+Title=di-jet $p_\perp$ ratio in PbPb 0-20, $300\ \mathrm{GeV} < p_{\perp,1} < 500\ \mathrm{GeV}$
+XLabel=$p_{\perp,2}/p_{\perp,1}$
+YLabel=event fraction
# END PLOT
# BEGIN PLOT /CMS_2012_I1090064/d06-x01-y0*
YMin=0.5
XMax=300.
LogY=0
-FullRange=1
-#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d06-x01-y01 /home/jewel/trunk/data/jewel-2.0.0_vac_fixed_5-350.aida/CMS_2012_I1090064/d06-x01-y01
# END PLOT
+# BEGIN PLOT /CMS_2012_I1090064/d06-x01-y01
+Title=mean di-jet $p_\perp$ ratio in pp
+XLabel=$p_{\perp,1}$ [GeV]
+YLabel=$\langle p_{\perp,2}/p_{\perp,1}\rangle$
+# END PLOT
+# BEGIN PLOT /CMS_2012_I1090064/d06-x01-y02
+Title=mean di-jet $p_\perp$ ratio in PbPb 0-20
+XLabel=$p_{\perp,1}$ [GeV]
+YLabel=$\langle p_{\perp,2}/p_{\perp,1}\rangle$
+# END PLOT
+
+# BEGIN PLOT /CMS_2012_I1090064/d06-x01-y03
+Title=mean di-jet $p_\perp$ ratio in PbPb 20-50
+XLabel=$p_{\perp,1}$ [GeV]
+YLabel=$\langle p_{\perp,2}/p_{\perp,1}\rangle$
+# END PLOT
+# BEGIN PLOT /CMS_2012_I1090064/d06-x01-y04
+Title=mean di-jet $p_\perp$ ratio in PbPb 50-100
+XLabel=$p_{\perp,1}$ [GeV]
+YLabel=$\langle p_{\perp,2}/p_{\perp,1}\rangle$
+# END PLOT
Index: trunk/code/CMS_HIN_12_004.cc
===================================================================
--- trunk/code/CMS_HIN_12_004.cc (revision 461)
+++ trunk/code/CMS_HIN_12_004.cc (revision 462)
@@ -1,109 +1,109 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/FastJets.hh"
/// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder...
namespace Rivet {
class CMS_HIN_12_004 : public Analysis {
public:
/// @name Constructors etc.
//@{
/// Constructor
CMS_HIN_12_004()
: Analysis("CMS_HIN_12_004")
{
- setNeedsCrossSection(true);
+ //setNeedsCrossSection(true);
}
//@}
public:
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void init() {
_centedges += 0.0, 0.05, 0.1, 0.3, 0.5, 0.7, 0.9;
_ncentbins = _centedges.size()-1;
FinalState fs(-5.0, 5.0, 0.*GeV);
addProjection(fs, "FS");
FastJets fj(fs, FastJets::ANTIKT, 0.3);
fj.useInvisibles();
addProjection(fj, "jets03");
for (size_t i = 0; i < _ncentbins; ++i) {
_sumwtcentbins[i]=0.;
stringstream ss; ss << "centbin_" << i;
const string pname = ss.str();
_histos[i] = bookHisto1D(pname,refData(1, 1, i+1));
_ratios[i] = bookHisto1D(1, 1, i+1);
}
}
/// Perform the per-event analysis
void analyze(const Event& event) {
const double weight = event.weight();
const double cent = (event.genEvent()->heavy_ion()?event.genEvent()->heavy_ion()->impact_parameter():-1.);
Cut cuts = Cuts::abseta < 2.0 && Cuts::pT > 100.*GeV && Cuts::pT < 300.*GeV;
const Jets jets = applyProjection<FastJets>(event, "jets03").jetsByPt(cuts);
for (size_t i = 0; i < _ncentbins; ++i) {
if (cent < 0. || (cent >= _centedges[i] && cent < _centedges[i+1])) {
_sumwtcentbins[i]+=weight;
// std::cout<<i<<" "<<cent<<" "<<jets.size()<<std::endl;
foreach (Jet jet, jets){
// std::cout<<jet.momentum().pT()<<std::endl;
_histos[i]->fill(jet.momentum().pT(),weight);
_ratios[i]->fill(jet.momentum().pT(),weight);
}
}
}
}
/// Normalise histograms etc., after the run
void finalize() {
for (size_t i = 0; i < _ncentbins; ++i) {
scale(_histos[i],crossSection()/(_sumwtcentbins[i]>0.?_sumwtcentbins[i]:1.));
scale(_ratios[i],crossSection()/(_sumwtcentbins[i]>0.?_sumwtcentbins[i]:1.));
}
}
//@}
private:
vector<double> _centedges;
double _sumwtcentbins[6];
size_t _ncentbins;
private:
Histo1DPtr _histos[6];
Histo1DPtr _ratios[6];
};
// The hook for the plugin system
DECLARE_RIVET_PLUGIN(CMS_HIN_12_004);
}
File Metadata
Details
Attached
Mime Type
text/x-diff
Expires
Tue, Nov 19, 4:43 PM (1 d, 9 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3805221
Default Alt Text
(150 KB)
Attached To
rJEWELSVN jewelsvn
Event Timeline
Log In to Comment