FParameter.hh

// -*- C++ -*-
#ifndef RIVET_FParameter_HH
#define RIVET_FParameter_HH
 
#include "Rivet/Projection.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Event.hh"
 
namespace Rivet {
 
 
  class FParameter : public Projection {
  public:
 
 
    FParameter(const FinalState& fsp);
 
    RIVET_DEFAULT_PROJ_CLONE(FParameter);
 
 
    using Projection::operator =;
 
 
  protected:
 
    void project(const Event& e);
 
    //CmpState compare(const Projection& p) const;
    // Taken from Thrust.hh
    CmpState compare(const Projection& p) const {
      return mkNamedPCmp(p, "FS");
    }
 
 
  public:
 
    void clear();
 
    double F() const { return lambda1() >= lambda2() ? lambda2()/lambda1() : lambda1()/lambda2(); }
 
    double lambda1() const { return _lambdas[0]; }
    double lambda2() const { return _lambdas[1]; }
 
 
 
    void calc(const FinalState& fs);
 
    void calc(const vector<Particle>& fsparticles);
 
    void calc(const vector<FourMomentum>& fsmomenta);
 
    void calc(const vector<Vector3>& fsmomenta);
 
 
 
  protected:
 
    vector<double> _lambdas;
 
 
  protected:
 
    void _calcFParameter(const vector<Vector3>& fsmomenta);
 
  };
 
 
}
 
#endif