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& vectorSumPt() const { return _vpt; }

    double scalarSumPt() const { return _spt; }

 


    const Vector3& vectorSumEt() const { return _vet; }

    double scalarSumEt() 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