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