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, edm::ValueMap< float > const *weights=nullptr)
virtual	~METAlgo ()

Detailed Description

Description: Calculates MET for given input

Implementation: [Notes on implementation]

Definition at line 29 of file METAlgo.h.

Constructor & Destructor Documentation

METAlgo()

METAlgo::METAlgo ( )

inline

Definition at line 31 of file METAlgo.h.

{}

~METAlgo()

virtual METAlgo::~METAlgo ( )

inlinevirtual

Definition at line 32 of file METAlgo.h.

{}

Member Function Documentation

run()

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

Definition at line 16 of file METAlgo.cc.

References HLT_2023v12_cff::candidates, EgHLTOffHistBins_cfi::et, HLT_2023v12_cff::globalThreshold, BTaggingMonitor_cfi::met, runTheMatrix::ret, mps_merge::weight, and hltDeepSecondaryVertexTagInfosPFPuppi_cfi::weights.

Referenced by reco::PFMETProducerMVA::produce().

                                                              {
  math::XYZTLorentzVector p4;
  for (auto const& candPtr : candidates.ptrs()) {
    const reco::Candidate* cand = candPtr.get();
    float weight = (weights != nullptr) ? (*weights)[candPtr] : 1.0;
    if (!(cand->et() * weight > globalThreshold))
      continue;
    p4 += cand->p4() * weight;
  }
  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 const& candPtr : candidates.ptrs()) {
    const reco::Candidate* cand = candPtr.get();
    float weight = (weights != nullptr) ? (*weights)[candPtr] : 1.0;
    if (!(cand->et() * weight > globalThreshold))
      continue;
    et += cand->et() * weight;
  }

  ret.sumet = et;

  return ret;
}