BeamThrust.hh

// -*- C++ -*-
#ifndef RIVET_BeamThrust_HH
#define RIVET_BeamThrust_HH
 
#include "Rivet/Projection.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Event.hh"
 
namespace Rivet {
 
 
  class BeamThrust : public Projection {
  public:
 
    BeamThrust() {}
 
    BeamThrust(const FinalState& fsp) {
      setName("BeamThrust");
      declare(fsp, "FS");
    }
 
    RIVET_DEFAULT_PROJ_CLONE(BeamThrust);
 
    using Projection::operator =;
 
 
  protected:
 
    void project(const Event& e) {
      const vector<Particle> ps
        = apply<FinalState>(e, "FS").particles();
      calc(ps);
    }
 
    CmpState compare(const Projection& p) const {
      return mkNamedPCmp(p, "FS");
    }
 
  public:
 
    double beamthrust() const { return _beamthrust; }
 
 
    void calc(const FinalState& fs);
 
    void calc(const vector<Particle>& fsparticles);
 
    void calc(const vector<FourMomentum>& fsmomenta);
 
 
 
  protected:
 
    double _beamthrust;
 
 
  protected:
 
    void _calcBeamThrust(const vector<FourMomentum>& fsmomenta);
 
  };
 
}
 
#endif