NoPileUpPFMEtProducer.cc File Reference

#include "RecoMET/METPUSubtraction/plugins/NoPileUpPFMEtProducer.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "RecoMET/METAlgorithms/interface/significanceAlgo.h"
#include "DataFormats/METReco/interface/SigInputObj.h"
#include <cmath>
#include "FWCore/Framework/interface/MakerMacros.h"

Go to the source code of this file.

Functions
void	addToCommonMETData (CommonMETData &metData, const reco::Candidate::LorentzVector &p4)
reco::METCovMatrix	computePFMEtSignificance (const std::vector< metsig::SigInputObj > &metSignObjects)
void	finalizeCommonMETData (CommonMETData &metData)
int	findBestMatchingLepton (const std::vector< reco::Candidate::LorentzVector > &leptons, const reco::Candidate::LorentzVector &p4_ref)
void	initializeCommonMETData (CommonMETData &metData)
void	printCommonMETData (const std::string &label, const CommonMETData &metData)
void	printMVAMEtJetInfo (const std::string &label, int idx, const reco::PUSubMETCandInfo &jet)
void	printMVAMEtPFCandInfo (const std::string &label, int idx, const reco::PUSubMETCandInfo &pfCand)
void	printP4 (const std::string &label_part1, int idx, const std::string &label_part2, const reco::Candidate &candidate)
static const edm::ParameterSetDescriptionFillerPluginFactory::PMaker< edm::ParameterSetDescriptionFiller< NoPileUpPFMEtProducer > >	s_filler__LINE__ ("NoPileUpPFMEtProducer")
static const edm::MakerPluginFactory::PMaker< edm::WorkerMaker< NoPileUpPFMEtProducer > >	s_maker__LINE__ ("NoPileUpPFMEtProducer")
void	scaleAndAddPFMEtSignObjects (std::vector< metsig::SigInputObj > &metSignObjects_scaled, const std::vector< metsig::SigInputObj > &metSignObjects, double sf, double sfMin, double sfMax)

Variables
const double	defaultPFMEtResolutionX = 10.
const double	defaultPFMEtResolutionY = 10.
const double	epsilon = 1.e-9

Function Documentation

void addToCommonMETData	(	CommonMETData &	metData,
		const reco::Candidate::LorentzVector &	p4
	)

Definition at line 94 of file NoPileUpPFMEtProducer.cc.

References CommonMETData::mex, CommonMETData::mey, CommonMETData::mez, and CommonMETData::sumet.

Referenced by NoPileUpPFMEtProducer::produce().

{
  metData.mex   += p4.px();
  metData.mey   += p4.py();
  metData.mez   += p4.pz();
  metData.sumet += p4.pt();
}

reco::METCovMatrix computePFMEtSignificance

(

const std::vector< metsig::SigInputObj > &

metSignObjects

)

Definition at line 145 of file NoPileUpPFMEtProducer.cc.

References funct::abs(), metsig::significanceAlgo::addObjects(), defaultPFMEtResolutionX, defaultPFMEtResolutionY, epsilon, and metsig::significanceAlgo::getSignifMatrix().

Referenced by NoPileUpPFMEtProducer::produce().

{
  reco::METCovMatrix pfMEtCov;
  if ( metSignObjects.size() >= 2 ) {
    metsig::significanceAlgo pfMEtSignAlgorithm;
    pfMEtSignAlgorithm.addObjects(metSignObjects);
    pfMEtCov = pfMEtSignAlgorithm.getSignifMatrix();
  } 
  
  double det = 0;
  pfMEtCov.Det(det);
  if ( std::abs(det) < epsilon ) {
    edm::LogWarning("computePFMEtSignificance") 
      << "Inversion of PFMEt covariance matrix failed, det = " << det
      << " --> replacing covariance matrix by resolution defaults !!";    
    pfMEtCov(0,0) = defaultPFMEtResolutionX*defaultPFMEtResolutionX;
    pfMEtCov(0,1) = 0.;
    pfMEtCov(1,0) = 0.;
    pfMEtCov(1,1) = defaultPFMEtResolutionY*defaultPFMEtResolutionY;
  }

  return pfMEtCov;
}

void finalizeCommonMETData

(

CommonMETData &

metData

)

Definition at line 102 of file NoPileUpPFMEtProducer.cc.

References CommonMETData::met, CommonMETData::mex, CommonMETData::mey, CommonMETData::phi, and mathSSE::sqrt().

Referenced by NoPileUpPFMEtProducer::produce().

{
  metData.met = sqrt(metData.mex*metData.mex + metData.mey*metData.mey);
  metData.phi = atan2(metData.mey, metData.mex);
}

int findBestMatchingLepton	(	const std::vector< reco::Candidate::LorentzVector > &	leptons,
		const reco::Candidate::LorentzVector &	p4_ref
	)

Definition at line 108 of file NoPileUpPFMEtProducer.cc.

References reco::deltaR2().

Referenced by NoPileUpPFMEtProducer::produce().

{
  int leptonIdx_dR2min = -1;
  double dR2min = 1.e+3;
  int leptonIdx = 0;
  for ( std::vector<reco::Candidate::LorentzVector>::const_iterator lepton = leptons.begin();
    lepton != leptons.end(); ++lepton ) {
    double dR2 = deltaR2(*lepton, p4_ref);
    if ( leptonIdx_dR2min == -1 || dR2 < dR2min ) {
      leptonIdx_dR2min = leptonIdx;
      dR2min = dR2;
    }
    ++leptonIdx;
  }
  assert(leptonIdx_dR2min >= 0 && leptonIdx_dR2min < (int)leptons.size());
  return leptonIdx_dR2min;
}

void initializeCommonMETData

(

CommonMETData &

metData

)

Definition at line 84 of file NoPileUpPFMEtProducer.cc.

References CommonMETData::met, CommonMETData::mex, CommonMETData::mey, CommonMETData::mez, CommonMETData::phi, and CommonMETData::sumet.

Referenced by NoPileUpPFMEtProducer::produce().

{
  metData.met   = 0.;
  metData.mex   = 0.;
  metData.mey   = 0.;
  metData.mez   = 0.;
  metData.sumet = 0.;
  metData.phi   = 0.;
}

void printCommonMETData	(	const std::string &	label,
		const CommonMETData &	metData
	)

Definition at line 175 of file NoPileUpPFMEtProducer.cc.

References gather_cfg::cout, CommonMETData::mex, CommonMETData::mey, and CommonMETData::sumet.

{
  std::cout << label << ": Px = " << metData.mex << ", Py = " << metData.mey << ", sumEt = " << metData.sumet << std::endl;
}

void printMVAMEtJetInfo	(	const std::string &	label,
		int	idx,
		const reco::PUSubMETCandInfo &	jet
	)

Definition at line 180 of file NoPileUpPFMEtProducer.cc.

References gather_cfg::cout, reco::PUSubMETCandInfo::kHS, reco::PUSubMETCandInfo::kPU, reco::PUSubMETCandInfo::p4(), reco::PUSubMETCandInfo::passesLooseJetId(), and reco::PUSubMETCandInfo::type().

{
  std::cout << label << " #" << idx << " (";
  if      ( jet.type() == reco::PUSubMETCandInfo::kHS ) std::cout << "no-PU";
  else if ( jet.type() == reco::PUSubMETCandInfo::kPU   ) std::cout << "PU";
  std::cout << "): Pt = " << jet.p4().pt() << ", eta = " << jet.p4().eta() << ", phi = " << jet.p4().phi();
  std::cout << " id. flags: anti-noise = " << jet.passesLooseJetId() << std::endl;
  std::cout << std::endl;
}

void printMVAMEtPFCandInfo	(	const std::string &	label,
		int	idx,
		const reco::PUSubMETCandInfo &	pfCand
	)

Definition at line 190 of file NoPileUpPFMEtProducer.cc.

References gather_cfg::cout, reco::PUSubMETCandInfo::isWithinJet(), reco::PUSubMETCandInfo::kChHS, reco::PUSubMETCandInfo::kChPU, reco::PUSubMETCandInfo::kNeutral, reco::PUSubMETCandInfo::p4(), reco::PUSubMETCandInfo::passesLooseJetId(), AlCaHLTBitMon_QueryRunRegistry::string, and reco::PUSubMETCandInfo::type().

{
  std::cout << label << " #" << idx << " (";
  if      ( pfCand.type() == reco::PUSubMETCandInfo::kChHS ) std::cout << "no-PU charged";
  else if ( pfCand.type() == reco::PUSubMETCandInfo::kChPU   ) std::cout << "PU charged";
  else if ( pfCand.type() == reco::PUSubMETCandInfo::kNeutral         ) std::cout << "neutral";
  std::cout << "): Pt = " << pfCand.p4().pt() << ", eta = " << pfCand.p4().eta() << ", phi = " << pfCand.p4().phi();
  std::string isWithinJet_string;
  if ( pfCand.isWithinJet() ) isWithinJet_string = "true";
  else isWithinJet_string = "false";
  std::cout << " (isWithinJet = " << isWithinJet_string << ")";
  if ( pfCand.isWithinJet() ) std::cout << " Jet id. flags: anti-noise = " << pfCand.passesLooseJetId() << std::endl;
  std::cout << std::endl;
}

void printP4	(	const std::string &	label_part1,
		int	idx,
		const std::string &	label_part2,
		const reco::Candidate &	candidate
	)

Definition at line 169 of file NoPileUpPFMEtProducer.cc.

References reco::Candidate::charge(), gather_cfg::cout, reco::Candidate::eta(), reco::Candidate::phi(), and reco::Candidate::pt().

{
  std::cout << label_part1 << " #" << idx << label_part2 << ": Pt = " << candidate.pt() << ", eta = " << candidate.eta() << ", phi = " << candidate.phi() 
        << " (charge = " << candidate.charge() << ")" << std::endl;
}

static const edm::ParameterSetDescriptionFillerPluginFactory::PMaker<edm::ParameterSetDescriptionFiller< NoPileUpPFMEtProducer > > s_filler__LINE__ ( "NoPileUpPFMEtProducer"  )

static

static const edm::MakerPluginFactory ::PMaker< edm::WorkerMaker< NoPileUpPFMEtProducer > > s_maker__LINE__ ( "NoPileUpPFMEtProducer"  )

static

void scaleAndAddPFMEtSignObjects	(	std::vector< metsig::SigInputObj > &	metSignObjects_scaled,
		const std::vector< metsig::SigInputObj > &	metSignObjects,
		double	sf,
		double	sfMin,
		double	sfMax
	)

Definition at line 126 of file NoPileUpPFMEtProducer.cc.

References metsig::SigInputObj::set().

Referenced by NoPileUpPFMEtProducer::produce().

{
  double sf_value = sf;
  if ( sf_value > sfMax ) sf_value = sfMax;
  if ( sf_value < sfMin ) sf_value = sfMin;
  for ( std::vector<metsig::SigInputObj>::const_iterator metSignObject = metSignObjects.begin();
    metSignObject != metSignObjects.end(); ++metSignObject ) {
    metsig::SigInputObj metSignObject_scaled;
    metSignObject_scaled.set(
      metSignObject->get_type(), 
      sf_value*metSignObject->get_energy(), 
      metSignObject->get_phi(),
      sf_value*metSignObject->get_sigma_e(), 
      metSignObject->get_sigma_tan());
    metSignObjects_scaled.push_back(metSignObject_scaled);
  }
}

Variable Documentation

const double defaultPFMEtResolutionX = 10.

Definition at line 11 of file NoPileUpPFMEtProducer.cc.

Referenced by computePFMEtSignificance().

const double defaultPFMEtResolutionY = 10.

Definition at line 12 of file NoPileUpPFMEtProducer.cc.

Referenced by computePFMEtSignificance().

const double epsilon = 1.e-9

Definition at line 14 of file NoPileUpPFMEtProducer.cc.

Referenced by computePFMEtSignificance().