MissingMomentum.hh

// -*- C++ -*-
#ifndef RIVET_MissingMomentum_HH
#define RIVET_MissingMomentum_HH

#include "Rivet/Config/RivetCommon.hh"
#include "Rivet/Projection.hh"
#include "Rivet/Projections/VisibleFinalState.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Event.hh"

namespace Rivet {


  class MissingMomentum : public Projection {
  public:

    MissingMomentum(const Cut& c=Cuts::open()) {
      setName("MissingMomentum");
      FinalState fs(c);
      addProjection(fs, "FS");
      addProjection(VisibleFinalState(fs), "VisibleFS");
    }


    MissingMomentum(const FinalState& fs) {
      setName("MissingMomentum");
      addProjection(fs, "FS");
      addProjection(VisibleFinalState(fs), "VisibleFS");
    }


    DEFAULT_RIVET_PROJ_CLONE(MissingMomentum);




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

    const Vector3 vectorMissingPt() const { return -vectorPt(); }
    // Alias
    const Vector3 vectorMPT() const { return vectorMissingPt(); }

    double missingPt() const { return vectorPt().mod(); }
    // /// Alias for missingPt
    // double mpt() const { return missingPt(); }

    double scalarPt() const { return _spt; }
    // /// Alias for scalarPt
    // double spt() const { return scalarPt(); }





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

    const Vector3 vectorMissingEt() const { return -vectorEt(); }
    // Alias
    const Vector3 vectorMET() const { return vectorMissingEt(); }

    double missingEt() const { return vectorEt().mod(); }
    double met() const { return missingEt(); }

    double scalarEt() const { return _set; }
    double set() const { return scalarEt(); }



  public:

    void clear();


  protected:

    void project(const Event& e);

    int compare(const Projection& p) const;


  private:

    FourMomentum _momentum;

    double _set, _spt;

    Vector3 _vet, _vpt;

  };


}

#endif