2 #ifndef RIVET_JetShape_HH
3 #define RIVET_JetShape_HH
5 #include "Rivet/Rivet.hh"
6 #include "Rivet/Projection.hh"
7 #include "Rivet/Projections/JetAlg.hh"
8 #include "Rivet/Particle.hh"
9 #include "Rivet/Event.hh"
10 #include "Rivet/Tools/Utils.hh"
53 double rmin,
double rmax,
size_t nbins,
54 double ptmin=0,
double ptmax=MAXDOUBLE,
55 double absrapmin=-MAXDOUBLE,
double absrapmax=-MAXDOUBLE,
60 double ptmin=0,
double ptmax=MAXDOUBLE,
61 double absrapmin=-MAXDOUBLE,
double absrapmax=-MAXDOUBLE,
85 return _binedges.size() - 1;
90 return _diffjetshapes.size();
95 return _binedges.front();
100 return _binedges.back();
105 return _ptcuts.first;
110 return _ptcuts.second;
116 return _binedges[rbin];
122 return _binedges[rbin+1];
129 return (_binedges[rbin] + _binedges[rbin+1])/2.0;
136 return _diffjetshapes[ijet][rbin];
144 for (
size_t i = 0; i <= rbin; ++i) {
145 rtn += _diffjetshapes[ijet][i];
174 vector<double> _binedges;
177 pair<double, double> _ptcuts;
180 pair<double, double> _rapcuts;
192 vector< vector<double> > _diffjetshapes;