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;
Definition: MC_JetAnalysis.hh:9
double rBinMin(size_t rbin) const
Central value for bin rbin.
Definition: JetShape.hh:114
double rBinMid(size_t rbin) const
Central value for bin rbin.
Definition: JetShape.hh:126
double ptMax() const
value.
Definition: JetShape.hh:109
int compare(const Projection &p) const
Compare projections.
Definition: JetShape.cc:39
bool inRange(NUM value, NUM low, NUM high, RangeBoundary lowbound=CLOSED, RangeBoundary highbound=OPEN)
Determine if value is in the range low to high, for floating point numbers.
Definition: MathUtils.hh:109
size_t numJets() const
Number of jets which passed cuts.
Definition: JetShape.hh:89
void project(const Event &e)
Apply the projection to the event.
Definition: JetShape.cc:105
virtual const Projection * clone() const
Clone on the heap.
Definition: JetShape.hh:65
double intJetShape(size_t ijet, size_t rbin) const
Return value of integrated jet shape profile histo bin.
Definition: JetShape.hh:140
JetShape(const JetAlg &jetalg, double rmin, double rmax, size_t nbins, double ptmin=0, double ptmax=MAXDOUBLE, double absrapmin=-MAXDOUBLE, double absrapmax=-MAXDOUBLE, RapScheme rapscheme=RAPIDITY)
Constructor from histo range and number of bins.
Definition: JetShape.cc:9
std::vector< Jet > Jets
Typedef for a collection of Jet objects.
Definition: Jet.hh:177
double diffJetShape(size_t ijet, size_t rbin) const
Return value of differential jet shape profile histo bin.
Definition: JetShape.hh:133
void calc(const Jets &jets)
Do the calculation directly on a supplied collection of Jet objects.
Definition: JetShape.cc:62
size_t numBins() const
Number of equidistant radius bins.
Definition: JetShape.hh:84
double rMax() const
value.
Definition: JetShape.hh:99
double rMin() const
value.
Definition: JetShape.hh:94
RapScheme
Enum for rapidity variable to be used in calculating , applying rapidity cuts, etc.
Definition: MathHeader.hh:60
double ptMin() const
value.
Definition: JetShape.hh:104
double rBinMax(size_t rbin) const
Central value for bin rbin.
Definition: JetShape.hh:120
void clear()
Reset projection between events.
Definition: JetShape.cc:57
Base class for all Rivet projections.
Definition: Projection.hh:28
Abstract base class for projections which can return a set of Jets.
Definition: JetAlg.hh:65
Calculate the jet shape.
Definition: JetShape.hh:45