MET.h

Go to the documentation of this file.

// -*- C++ -*-

// Package:    METReco
// Class:      MET
//
// Original authors: Michael Schmitt, Richard Cavanaugh The University of Florida
// changes by: Freya Blekman, Cornell University
//

//____________________________________________________________________________||
#ifndef METRECO_MET_H
#define METRECO_MET_H

//____________________________________________________________________________||
#include "DataFormats/RecoCandidate/interface/RecoCandidate.h"
#include "DataFormats/METReco/interface/CorrMETData.h"
#include <cmath>
#include <vector>
#include <cstring>
#include <Math/SMatrix.h>
#include <Math/SVector.h>

//____________________________________________________________________________||
namespace reco
{
  typedef ROOT::Math::SMatrix<double,2> METCovMatrix;
  
  class MET : public RecoCandidate
  {
  public:

    MET();
    MET( const LorentzVector& p4_, const Point& vtx_ );
    MET( double sumet_, const LorentzVector& p4_, const Point& vtx_ );
    MET( double sumet_, const std::vector<CorrMETData>& corr_,
     const LorentzVector& p4_, const Point& vtx_ );

    MET * clone() const override;

    //________________________________________________________________________||
    //scalar sum of transverse energy over all objects
    double sumEt() const { return sumet; }
    //MET Significance = MET / std::sqrt(SumET)
    double mEtSig() const { return ( sumet ? (this->et() / std::sqrt(sumet)) : (0.0) ); }
    //real MET significance
    double significance() const;
    //longitudinal component of the vector sum of energy over all object
    //(useful for data quality monitoring)
    double e_longitudinal() const {return elongit; }

    //________________________________________________________________________||
    //Define different methods for the corrections to individual MET elements
    std::vector<double> dmEx() const ;
    std::vector<double> dmEy() const ;
    std::vector<double> dsumEt() const ;
    std::vector<CorrMETData> mEtCorr() const { return corr; }

    //________________________________________________________________________||
    void setSignificanceMatrix(const reco::METCovMatrix& matrix);
    reco::METCovMatrix getSignificanceMatrix(void) const;

  private:
    bool overlap( const Candidate & ) const override;
    double sumet;
    double elongit;
    // bookkeeping for the significance
    double signif_dxx;
    double signif_dyy;
    double signif_dyx;
    double signif_dxy;
    std::vector<CorrMETData> corr;
  };
}

//____________________________________________________________________________||
#endif // METRECO_MET_H