MissingMomentum.hh

// -*- C++ -*-
#ifndef RIVET_MissingMomentum_HH
#define RIVET_MissingMomentum_HH
 
#include "Rivet/Config/RivetCommon.hh"
#include "Rivet/Projection.hh"
#include "Rivet/Projections/METFinder.hh"
#include "Rivet/Projections/VisibleFinalState.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Event.hh"
 
namespace Rivet {
 
 
  class MissingMomentum : public METFinder {
  public:
 
    using METFinder::operator=;
 
    MissingMomentum(const FinalState& fs) {
      setName("MissingMomentum");
      declare(fs, "FS");
      declare(VisibleFinalState(fs), "VisibleFS");
    }
 
    MissingMomentum(const Cut& c=Cuts::open())
      : MissingMomentum(FinalState(c))
    {    }
 
 
    RIVET_DEFAULT_PROJ_CLONE(MissingMomentum);
 
    using Projection::operator =;
 
 
 
    const FourMomentum visibleMomentum(double mass=0*GeV) const;
    const FourMomentum visibleMom(double mass=0*GeV) const { return visibleMomentum(mass); }
 
    const FourMomentum missingMomentum(double mass=0*GeV) const { return visibleMomentum(mass).reverse(); }
    const FourMomentum missingMom(double mass=0*GeV) const { return missingMomentum(mass); }
 
 
 
 
    const ThreeMomentum& vectorPt() const { return _vpt; }
 
    const ThreeMomentum vectorPtMiss() const { return -vectorPt(); }
 
    double scalarPtMiss() const { return vectorPtMiss().pT(); }
 
    double scalarPt() const { return _spt; }
    // /// Alias for scalarPt
    // double spt() const { return scalarPt(); }
 
 
 
 
    const Vector3& vectorEt() const { return _vet; }
 
    const Vector3 vectorEtMiss() const { return -vectorEt(); }
 
    double scalarEtMiss() const { return vectorEtMiss().perp(); }
 
    double scalarEt() const { return _set; }
 
 
 
    void clear();
 
 
  protected:
 
    void project(const Event& e);
 
    CmpState compare(const Projection& p) const;
 
 
  protected:
 
    FourMomentum _momentum;
 
    double _set, _spt;
 
    Vector3 _vet;
    ThreeMomentum _vpt;
 
  };
 
 
 
  using MissingMom = MissingMomentum;
 
 
}
 
#endif