Rivet  1.8.3
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 
10 namespace Rivet {
11 
12 
44  class Thrust : public AxesDefinition {
45  public:
46 
48  Thrust() {}
49 
50  Thrust(const FinalState& fsp) {
51  setName("Thrust");
52  addProjection(fsp, "FS");
53  }
54 
56  virtual const Projection* clone() const {
57  return new Thrust(*this);
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  int 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  private:
126 
128  vector<double> _thrusts;
129 
131  vector<Vector3> _thrustAxes;
132 
133  private:
134 
136  void _calcThrust(const vector<Vector3>& fsmomenta);
137 
138  };
139 
140 }
141 
142 #endif
const Vector3 & axis3() const
The least significant (&quot;minor&quot;) axis.
Definition: Thrust.hh:100
double thrustMajor() const
The thrust major scalar, , (thrust along thrust major axis).
Definition: Thrust.hh:81
void setName(const std::string &name)
Used by derived classes to set their name.
Definition: Projection.hh:120
const Vector3 & thrustMajorAxis() const
The thrust major axis (axis of max thrust perpendicular to thrust axis).
Definition: Thrust.hh:92
const Vector3 & axis1() const
AxesDefinition axis accessors.
Definition: Thrust.hh:98
const Vector3 & thrustMinorAxis() const
The thrust minor axis (axis perpendicular to thrust and thrust major).
Definition: Thrust.hh:94
Base class for projections which define a spatial basis.
Definition: AxesDefinition.hh:20
virtual const Projection * clone() const
Clone on the heap.
Definition: Thrust.hh:56
Thrust()
Constructor.
Definition: Thrust.hh:48
double thrustMinor() const
The thrust minor scalar, , (thrust along thrust minor axis).
Definition: Thrust.hh:83
double oblateness() const
The oblateness, .
Definition: Thrust.hh:85
double thrust() const
Definition: Thrust.hh:79
Definition: Event.hh:22
void calc(const FinalState &fs)
Manually calculate the thrust, without engaging the caching system.
Definition: Thrust.cc:10
void project(const Event &e)
Perform the projection on the Event.
Definition: Thrust.hh:63
int compare(const Projection &p) const
Compare projections.
Definition: Thrust.hh:70
Get the e+ e- thrust basis and the thrust, thrust major and thrust minor scalars. ...
Definition: Thrust.hh:44
Project out all final-state particles in an event. Probably the most important projection in Rivet! ...
Definition: FinalState.hh:14
const PROJ & addProjection(const PROJ &proj, const std::string &name)
Definition: ProjectionApplier.hh:113
Cmp< Projection > mkNamedPCmp(const Projection &otherparent, const std::string &pname) const
Definition: Projection.cc:51
const Vector3 & axis2() const
The 2nd most significant (&quot;major&quot;) axis.
Definition: Thrust.hh:99
Base class for all Rivet projections.
Definition: Projection.hh:28
Three-dimensional specialisation of Vector.
Definition: Vector3.hh:26
const Vector3 & thrustAxis() const
Definition: Thrust.hh:90