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 90 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 143 of file NoPileUpPFMEtProducer.cc.

References funct::abs(), metsig::significanceAlgo::addObjects(), defaultPFMEtResolutionX, defaultPFMEtResolutionY, geometryDiff::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 97 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 102 of file NoPileUpPFMEtProducer.cc.

References cms::cuda::assert(), and 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 81 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 170 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 175 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 186 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 165 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 120 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(), and PFMEtSignInterfaceBase::operator()().

const double defaultPFMEtResolutionY = 10.

Definition at line 12 of file NoPileUpPFMEtProducer.cc.

Referenced by computePFMEtSignificance(), and PFMEtSignInterfaceBase::operator()().

const double epsilon = 1.e-9

Definition at line 14 of file NoPileUpPFMEtProducer.cc.