Rivet 3.1.9
Thrust.hh
1// -*- C++ -*-
2#ifndef RIVET_Thrust_HH
3#define RIVET_Thrust_HH
4
5#include "Rivet/Projection.hh"
6#include "Rivet/Projections/AxesDefinition.hh"
7#include "Rivet/Projections/FinalState.hh"
8#include "Rivet/Event.hh"
9
10namespace Rivet {
11
12
42 class Thrust : public AxesDefinition {
43 public:
44
46 Thrust() {}
47
48 Thrust(const FinalState& fsp) {
49 setName("Thrust");
50 declare(fsp, "FS");
51 }
52
55
57 using Projection::operator =;
58
59
60 protected:
61
63 void project(const Event& e) {
64 const vector<Particle> ps
65 = applyProjection<FinalState>(e, "FS").particles();
66 calc(ps);
67 }
68
70 CmpState compare(const Projection& p) const {
71 return mkNamedPCmp(p, "FS");
72 }
73
74
75 public:
76
79 double thrust() const { return _thrusts[0]; }
81 double thrustMajor() const { return _thrusts[1]; }
83 double thrustMinor() const { return _thrusts[2]; }
85 double oblateness() const { return _thrusts[1] - _thrusts[2]; }
87
90 const Vector3& thrustAxis() const { return _thrustAxes[0]; }
92 const Vector3& thrustMajorAxis() const { return _thrustAxes[1]; }
94 const Vector3& thrustMinorAxis() const { return _thrustAxes[2]; }
96
98 const Vector3& axis1() const { return thrustAxis(); }
99 const Vector3& axis2() const { return thrustMajorAxis(); }
100 const Vector3& axis3() const { return thrustMinorAxis(); }
102
103
104 public:
105
108
109
111 void calc(const FinalState& fs);
112
114 void calc(const vector<Particle>& fsparticles);
115
117 void calc(const vector<FourMomentum>& fsmomenta);
118
120 void calc(const vector<Vector3>& threeMomenta);
121
123
124
125 protected:
126
128 vector<double> _thrusts;
129
131 vector<Vector3> _thrustAxes;
132
133
134 protected:
135
137 void _calcThrust(const vector<Vector3>& fsmomenta);
138
139 };
140
141}
142
143#endif
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
const PROJ & declare(const PROJ &proj, const std::string &name)
Register a contained projection (user-facing version)
Definition ProjectionApplier.hh:170
Base class for all Rivet projections.
Definition Projection.hh:29
void setName(const std::string &name)
Used by derived classes to set their name.
Definition Projection.hh:142
Cmp< Projection > mkNamedPCmp(const Projection &otherparent, const std::string &pname) const
Get the e+ e- thrust basis and the thrust, thrust major and thrust minor scalars.
Definition Thrust.hh:42
DEFAULT_RIVET_PROJ_CLONE(Thrust)
Clone on the heap.
const Vector3 & thrustAxis() const
Definition Thrust.hh:90
double oblateness() const
The oblateness, .
Definition Thrust.hh:85
const Vector3 & thrustMinorAxis() const
The thrust minor axis (axis perpendicular to thrust and thrust major).
Definition Thrust.hh:94
const Vector3 & axis1() const
AxesDefinition axis accessors.
Definition Thrust.hh:98
void calc(const vector< Vector3 > &threeMomenta)
Manually calculate the thrust, without engaging the caching system.
void calc(const vector< FourMomentum > &fsmomenta)
Manually calculate the thrust, without engaging the caching system.
double thrustMajor() const
The thrust major scalar, , (thrust along thrust major axis).
Definition Thrust.hh:81
double thrust() const
Definition Thrust.hh:79
void project(const Event &e)
Perform the projection on the Event.
Definition Thrust.hh:63
CmpState compare(const Projection &p) const
Compare projections.
Definition Thrust.hh:70
double thrustMinor() const
The thrust minor scalar, , (thrust along thrust minor axis).
Definition Thrust.hh:83
void calc(const vector< Particle > &fsparticles)
Manually calculate the thrust, without engaging the caching system.
Thrust()
Constructor.
Definition Thrust.hh:46
const Vector3 & axis2() const
The 2nd most significant ("major") axis.
Definition Thrust.hh:99
void calc(const FinalState &fs)
Manually calculate the thrust, without engaging the caching system.
const Vector3 & thrustMajorAxis() const
The thrust major axis (axis of max thrust perpendicular to thrust axis).
Definition Thrust.hh:92
const Vector3 & axis3() const
The least significant ("minor") axis.
Definition Thrust.hh:100
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