METFinder.hh

// -*- C++ -*-
#ifndef RIVET_METFinder_HH
#define RIVET_METFinder_HH
 
#include "Rivet/Projection.hh"
 
namespace Rivet {
 

  class METFinder : public Projection {
  public:

    using Projection::operator=;


    virtual const Vector3& vectorSumPt() const = 0;
    virtual const Vector3& vectorPt() const {
      return vectorSumPt();
    }

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

    double missingPt() const {
      return vectorPt().mod();
    }

    double pTmiss() const {
      return missingPt();
    }

    double scalarPtMiss() const {
      return missingPt();
    }

    double mpt() const {
      return missingPt();
    }

    virtual double scalarSumPt() const = 0;
    double scalarPt() const {
      return scalarSumPt();
    }

    double spt() const {
      return scalarSumPt();
    }

 


    virtual const Vector3& vectorSumEt() const = 0;
    virtual const Vector3& vectorEt() const {
      return vectorSumEt();
    }

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

    double missingEt() const {
      return vectorMissingEt().mod();
    }

    double eTmiss() const {
      return missingEt();
    }

    double scalarEtMiss() const {
      return missingEt();
    }

    double met() const {
      return missingEt();
    }

    virtual double scalarSumEt() const = 0;
    double scalarEt() const {
      return scalarSumEt();
    }

    double set() const {
      return scalarSumEt();
    }

 

    virtual void reset() { }
  };
 
 
}
 
#endif