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"
87 double sphericity()
const {
return 3.0 / 2.0 * (lambda2() + lambda3()); }
89 double transSphericity()
const {
return 2.0 * lambda2() / ( lambda1() + lambda2() ); }
93 double aplanarity()
const {
return 3 / 2.0 * lambda3(); }
99 const Vector3& sphericityAxis()
const {
return _sphAxes[0]; }
115 double lambda1()
const {
return _lambdas[0]; }
116 double lambda2()
const {
return _lambdas[1]; }
117 double lambda3()
const {
return _lambdas[2]; }
126 void calc(
const FinalState& fs);
129 void calc(
const vector<Particle>& fsparticles);
132 void calc(
const vector<FourMomentum>& fsmomenta);
135 void calc(
const vector<Vector3>& fsmomenta);
143 vector<double> _lambdas;
146 vector<Vector3> _sphAxes;
149 const double _regparam;
154 void _calcSphericity(
const vector<Vector3>& fsmomenta);