METAlgo Class Reference

#include <RecoMET/METAlgorithms/interface/METAlgo.h>

Public Member Functions
	METAlgo ()
CommonMETData	run (const edm::View< reco::Candidate > &candidates, double globalThreshold=0.0)
virtual	~METAlgo ()

Detailed Description

Description: Calculates MET for given input

Implementation: [Notes on implementation]

Definition at line 28 of file METAlgo.h.

Constructor & Destructor Documentation

METAlgo::METAlgo ( )

inline

Definition at line 31 of file METAlgo.h.

{}

virtual METAlgo::~METAlgo ( )

inlinevirtual

Definition at line 32 of file METAlgo.h.

References objects.IsoTrackAnalyzer::candidates, met_cff::globalThreshold, and run().

{}

Member Function Documentation

CommonMETData METAlgo::run	(	const edm::View< reco::Candidate > &	candidates,
		double	globalThreshold = 0.0
	)

Definition at line 16 of file METAlgo.cc.

References edm::View< T >::begin(), edm::View< T >::end(), stringResolutionProvider_cfi::et, met_cff::globalThreshold, CommonMETData::met, RazorAnalyzer::met, CommonMETData::mex, CommonMETData::mey, CommonMETData::mez, p4, and CommonMETData::sumet.

Referenced by cms::ElseMETProducer::produce(), cms::PFClusterMETProducer::produce(), cms::CaloMETProducer::produce(), reco::PFMETProducerMVA::produce(), cms::PFMETProducer::produce(), and ~METAlgo().

{
  math::XYZTLorentzVector p4;
  for(auto cand = candidates.begin(); cand != candidates.end(); ++cand)
    {
      if( !(cand->et() > globalThreshold) ) continue;
      p4 += cand->p4();
    }
  math::XYZTLorentzVector met = -p4;


  CommonMETData ret;
  ret.mex   = met.Px();
  ret.mey   = met.Py();

  ret.mez   = met.Pz(); // included here since it might be useful
                        // for Data Quality Monitering as it should be
                        // symmetrically distributed about the origin

  ret.met   = met.Pt();

  double et = 0.0;
  for(auto cand = candidates.begin(); cand != candidates.end(); ++cand)
    {
      if( !(cand->et() > globalThreshold) ) continue;
      et += cand->et();
    }

  ret.sumet = et;

  return ret;
}