MissingMomentum.hh

// -*- C++ -*-
#ifndef RIVET_MissingMomentum_HH
#define RIVET_MissingMomentum_HH
 
#include "Rivet/Event.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Projection.hh"
#include "Rivet/Config/RivetCommon.hh"
#include "Rivet/Projections/METFinder.hh"
#include "Rivet/Projections/VisibleFinalState.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