2 #ifndef RIVET_Spherocity_HH
3 #define RIVET_Spherocity_HH
5 #include "Rivet/Projection.hh"
6 #include "Rivet/Projections/AxesDefinition.hh"
7 #include "Rivet/Projections/FinalState.hh"
8 #include "Rivet/Event.hh"
41 : _calculatedSpherocity(false)
57 const vector<Particle> ps
58 = applyProjection<FinalState>(e,
"FS").particles();
104 void calc(
const vector<Particle>& fsparticles);
107 void calc(
const vector<FourMomentum>& fsmomenta);
110 void calc(
const vector<Vector3>& threeMomenta);
118 vector<double> _spherocities;
121 vector<Vector3> _spherocityAxes;
124 bool _calculatedSpherocity;
130 void _calcSpherocity(
const vector<Vector3>& fsmomenta);