2#ifndef RIVET_Sphericity_HH
3#define RIVET_Sphericity_HH
5#include "Rivet/Projection.hh"
6#include "Rivet/Projections/AxesDefinition.hh"
7#include "Rivet/Projections/FinalState.hh"
8#include "Rivet/Event.hh"
9#include "Rivet/Jet.fhh"
69 using Projection::operator =;
89 double sphericity()
const {
return 3.0 / 2.0 * (lambda2() + lambda3()); }
91 double transSphericity()
const {
return 2.0 * lambda2() / ( lambda1() + lambda2() ); }
95 double aplanarity()
const {
return 3 / 2.0 * lambda3(); }
120 double lambda1()
const {
return _lambdas[0]; }
121 double lambda2()
const {
return _lambdas[1]; }
122 double lambda3()
const {
return _lambdas[2]; }
125 Vector3 mkEigenVector(Matrix3 A,
const double &lambda);
141 void calc(
const vector<FourMomentum>& momenta);
146 void calc(
const vector<Vector3>& momenta);
154 vector<double> _lambdas;
157 vector<Vector3> _sphAxes;
160 const double _regparam;
Base class for projections which define a spatial basis.
Definition AxesDefinition.hh:19
Representation of a HepMC event, and enabler of Projection caching.
Definition Event.hh:22
Project out all final-state particles in an event. Probably the most important projection in Rivet!
Definition FinalState.hh:12
Specialised vector of Jet objects.
Definition Jet.hh:23
Specialised vector of Particle objects.
Definition Particle.hh:25
Base class for all Rivet projections.
Definition Projection.hh:29
Calculate the sphericity event shape.
Definition Sphericity.hh:52
const Vector3 & axis3() const
The least significant ("minor") axis.
Definition Sphericity.hh:114
const Vector3 & axis1() const
Definition Sphericity.hh:112
void project(const Event &e)
Perform the projection on the Event.
double transSphericity() const
Transverse sphericity.
Definition Sphericity.hh:91
void calc(const Particles &particles)
Manually calculate the sphericity, without engaging the caching system.
CmpState compare(const Projection &p) const
Compare with other projections.
Sphericity(double rparam=2.0)
Constructor.
Definition Sphericity.hh:59
double planarity() const
Planarity.
Definition Sphericity.hh:93
void calc(const FinalState &fs)
Manually calculate the sphericity, without engaging the caching system.
void clear()
Reset the projection.
void calc(const vector< FourMomentum > &momenta)
Manually calculate the sphericity, without engaging the caching system.
const Vector3 & sphericityAxis() const
Sphericity axis.
Definition Sphericity.hh:102
double aplanarity() const
Aplanarity.
Definition Sphericity.hh:95
void calc(const Jets &jets)
Manually calculate the sphericity, without engaging the caching system.
const Vector3 & sphericityMinorAxis() const
Sphericity minor axis.
Definition Sphericity.hh:106
const Vector3 & axis2() const
The 2nd most significant ("major") axis.
Definition Sphericity.hh:113
const Vector3 & sphericityMajorAxis() const
Sphericity major axis.
Definition Sphericity.hh:104
double sphericity() const
Sphericity.
Definition Sphericity.hh:89
void calc(const vector< Vector3 > &momenta)
Manually calculate the sphericity, without engaging the caching system.
DEFAULT_RIVET_PROJ_CLONE(Sphericity)
Clone on the heap.
Three-dimensional specialisation of Vector.
Definition Vector3.hh:40
double p(const ParticleBase &p)
Unbound function access to p.
Definition ParticleBaseUtils.hh:684
Definition MC_Cent_pPb.hh:10