#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 <math.h>
#include "FWCore/Framework/interface/MakerMacros.h"

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::MVAMEtJetInfo &jet)

void	printMVAMEtPFCandInfo (const std::string &label, int idx, const reco::MVAMEtPFCandInfo &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 96 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 157 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 104 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 110 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 86 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 187 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::MVAMEtJetInfo &	jet
	)

Definition at line 192 of file NoPileUpPFMEtProducer.cc.

References gather_cfg::cout, reco::MVAMEtJetInfo::kNoPileUp, reco::MVAMEtJetInfo::kPileUp, reco::MVAMEtJetInfo::p4_, reco::MVAMEtJetInfo::passesLooseJetId_, and reco::MVAMEtJetInfo::type_.

 {
   std::cout << label << " #" << idx << " (";
   if      ( jet.type_ == reco::MVAMEtJetInfo::kNoPileUp ) std::cout << "no-PU";
   else if ( jet.type_ == reco::MVAMEtJetInfo::kPileUp   ) 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::MVAMEtPFCandInfo &	pfCand
	)

Definition at line 202 of file NoPileUpPFMEtProducer.cc.

References gather_cfg::cout, reco::MVAMEtPFCandInfo::isWithinJet_, reco::MVAMEtPFCandInfo::kNeutral, reco::MVAMEtPFCandInfo::kNoPileUpCharged, reco::MVAMEtPFCandInfo::kPileUpCharged, reco::MVAMEtPFCandInfo::p4_, reco::MVAMEtPFCandInfo::passesLooseJetId_, AlCaHLTBitMon_QueryRunRegistry::string, and reco::MVAMEtPFCandInfo::type_.

 {
   std::cout << label << " #" << idx << " (";
   if      ( pfCand.type_ == reco::MVAMEtPFCandInfo::kNoPileUpCharged ) std::cout << "no-PU charged";
   else if ( pfCand.type_ == reco::MVAMEtPFCandInfo::kPileUpCharged   ) std::cout << "PU charged";
   else if ( pfCand.type_ == reco::MVAMEtPFCandInfo::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 181 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 128 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 13 of file NoPileUpPFMEtProducer.cc.

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

const double defaultPFMEtResolutionY = 10.

Definition at line 14 of file NoPileUpPFMEtProducer.cc.

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

const double epsilon = 1.e-9

Definition at line 16 of file NoPileUpPFMEtProducer.cc.

Functions

Variables

Function Documentation

Variable Documentation