PFMETAlgorithmMVA Class Reference

#include <PFMETAlgorithmMVA.h>

Public Member Functions
void	evaluateMVA ()
double	getCovU1 () const
double	getCovU2 () const
double	getDPhi () const
reco::Candidate::LorentzVector	getMEt () const
const reco::METCovMatrix &	getMEtCov () const
double	getU () const
void	initialize (const edm::EventSetup &)
	PFMETAlgorithmMVA (const edm::ParameterSet &cfg)
void	print (std::ostream &) const
void	setHasPhotons (bool hasPhotons)
void	setInput (const std::vector< reco::PUSubMETCandInfo > &, const std::vector< reco::PUSubMETCandInfo > &, const std::vector< reco::PUSubMETCandInfo > &, const std::vector< reco::Vertex::Point > &)
	~PFMETAlgorithmMVA ()

Private Member Functions
void	evaluateCovU1 ()
void	evaluateCovU2 ()
void	evaluateDPhi ()
void	evaluateU ()
void	setInput (double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double)

Private Attributes
edm::ParameterSet	cfg_
double	chargedSumLeptonPx_
double	chargedSumLeptonPy_
bool	hasPhotons_
Float_t	jet1Eta_
Float_t	jet1Phi_
Float_t	jet1Pt_
Float_t	jet2Eta_
Float_t	jet2Phi_
Float_t	jet2Pt_
bool	loadMVAfromDB_
Float_t *	mvaInputCovU1_
Float_t *	mvaInputCovU2_
Float_t *	mvaInputDPhi_
Float_t *	mvaInputU_
reco::Candidate::LorentzVector	mvaMEt_
reco::METCovMatrix	mvaMEtCov_
std::string	mvaNameCovU1_
std::string	mvaNameCovU2_
std::string	mvaNameDPhi_
std::string	mvaNameU_
Float_t	mvaOutputCovU1_
Float_t	mvaOutputCovU2_
Float_t	mvaOutputDPhi_
Float_t	mvaOutputU_
const GBRForest *	mvaReaderCovU1_
const GBRForest *	mvaReaderCovU2_
const GBRForest *	mvaReaderDPhi_
const GBRForest *	mvaReaderU_
int	mvaType_
Float_t	npuPhi_
Float_t	npuSumEt_
Float_t	npuU_
Float_t	numJets_
Float_t	numJetsPtGt30_
Float_t	numVertices_
Float_t	pfPhi_
Float_t	pfSumEt_
Float_t	pfU_
Float_t	pucPhi_
Float_t	pucSumEt_
Float_t	pucU_
Float_t	puMEt_
Float_t	puPhi_
Float_t	puSumEt_
double	sumLeptonPx_
double	sumLeptonPy_
Float_t	tkPhi_
Float_t	tkSumEt_
Float_t	tkU_
MvaMEtUtilities	utils_

Detailed Description

MVA based algorithm for computing the particle-flow missing Et

Authors

Phil Harris, CERN Christian Veelken, LLR

Definition at line 31 of file PFMETAlgorithmMVA.h.

Constructor & Destructor Documentation

PFMETAlgorithmMVA::PFMETAlgorithmMVA

(

const edm::ParameterSet &

cfg

)

Definition at line 47 of file PFMETAlgorithmMVA.cc.

References edm::ParameterSet::getParameter(), kBaseline, loadMVAfromDB_, mvaInputCovU1_, mvaInputCovU2_, mvaInputDPhi_, mvaInputU_, and mvaType_.

  : utils_(cfg),
    mvaInputU_(nullptr),
    mvaInputDPhi_(nullptr),
    mvaInputCovU1_(nullptr),
    mvaInputCovU2_(nullptr),
    mvaReaderU_(nullptr),
    mvaReaderDPhi_(nullptr),
    mvaReaderCovU1_(nullptr),
    mvaReaderCovU2_(nullptr),
    cfg_(cfg)
{
  mvaType_ = kBaseline;

  loadMVAfromDB_ = cfg.getParameter<bool>("loadMVAfromDB");
  
  mvaInputU_     = new Float_t[25];
  mvaInputDPhi_  = new Float_t[23];
  mvaInputCovU1_ = new Float_t[26];
  mvaInputCovU2_ = new Float_t[26];
}

PFMETAlgorithmMVA::~PFMETAlgorithmMVA

(

)

Definition at line 69 of file PFMETAlgorithmMVA.cc.

References loadMVAfromDB_, mvaInputCovU1_, mvaInputCovU2_, mvaInputDPhi_, mvaInputU_, mvaReaderCovU1_, mvaReaderCovU2_, mvaReaderDPhi_, and mvaReaderU_.

{
  delete mvaInputU_;
  delete mvaInputDPhi_;
  delete mvaInputCovU1_;
  delete mvaInputCovU2_;

  if ( !loadMVAfromDB_ ) {
    delete mvaReaderU_;
    delete mvaReaderDPhi_;
    delete mvaReaderCovU1_;
    delete mvaReaderCovU2_;
  }
}

Member Function Documentation

void PFMETAlgorithmMVA::evaluateCovU1 ( )

private

Definition at line 314 of file PFMETAlgorithmMVA.cc.

References GBRForest::GetResponse(), jet1Eta_, jet1Phi_, jet1Pt_, jet2Eta_, jet2Phi_, jet2Pt_, mvaInputCovU1_, mvaOutputCovU1_, mvaOutputDPhi_, mvaOutputU_, mvaReaderCovU1_, npuPhi_, npuSumEt_, npuU_, numJets_, numJetsPtGt30_, numVertices_, pfPhi_, pfSumEt_, pfU_, pucPhi_, pucSumEt_, pucU_, puMEt_, puPhi_, puSumEt_, tkPhi_, tkSumEt_, and tkU_.

Referenced by evaluateMVA().

{ 
  mvaInputCovU1_[0]  = pfSumEt_; // PH: helps flattens response vs. Nvtx
  mvaInputCovU1_[1]  = numVertices_;
  mvaInputCovU1_[2]  = pfU_;
  mvaInputCovU1_[3]  = pfPhi_;
  mvaInputCovU1_[4]  = tkSumEt_;
  mvaInputCovU1_[5]  = tkU_;
  mvaInputCovU1_[6]  = tkPhi_;
  mvaInputCovU1_[7]  = npuSumEt_;
  mvaInputCovU1_[8]  = npuU_;
  mvaInputCovU1_[9]  = npuPhi_;
  mvaInputCovU1_[10] = puSumEt_;
  mvaInputCovU1_[11] = puMEt_;
  mvaInputCovU1_[12] = puPhi_;
  mvaInputCovU1_[13] = pucSumEt_;
  mvaInputCovU1_[14] = pucU_;
  mvaInputCovU1_[15] = pucPhi_;
  mvaInputCovU1_[16] = jet1Pt_;
  mvaInputCovU1_[17] = jet1Eta_;
  mvaInputCovU1_[18] = jet1Phi_;
  mvaInputCovU1_[19] = jet2Pt_;
  mvaInputCovU1_[20] = jet2Eta_;
  mvaInputCovU1_[21] = jet2Phi_;
  mvaInputCovU1_[22] = numJets_;
  mvaInputCovU1_[23] = numJetsPtGt30_;
  mvaInputCovU1_[24] = pfPhi_ + mvaOutputDPhi_;
  mvaInputCovU1_[25] = mvaOutputU_*pfU_;
  mvaOutputCovU1_    = mvaReaderCovU1_->GetResponse(mvaInputCovU1_)*mvaOutputU_*pfU_;
}

void PFMETAlgorithmMVA::evaluateCovU2 ( )

private

Definition at line 345 of file PFMETAlgorithmMVA.cc.

References GBRForest::GetResponse(), jet1Eta_, jet1Phi_, jet1Pt_, jet2Eta_, jet2Phi_, jet2Pt_, mvaInputCovU2_, mvaOutputCovU2_, mvaOutputDPhi_, mvaOutputU_, mvaReaderCovU2_, npuPhi_, npuSumEt_, npuU_, numJets_, numJetsPtGt30_, numVertices_, pfPhi_, pfSumEt_, pfU_, pucPhi_, pucSumEt_, pucU_, puMEt_, puPhi_, puSumEt_, tkPhi_, tkSumEt_, and tkU_.

Referenced by evaluateMVA().

{ 
  mvaInputCovU2_[0]  = pfSumEt_; // PH: helps flattens response vs. Nvtx
  mvaInputCovU2_[1]  = numVertices_;
  mvaInputCovU2_[2]  = pfU_;
  mvaInputCovU2_[3]  = pfPhi_;
  mvaInputCovU2_[4]  = tkSumEt_;
  mvaInputCovU2_[5]  = tkU_;
  mvaInputCovU2_[6]  = tkPhi_;
  mvaInputCovU2_[7]  = npuSumEt_;
  mvaInputCovU2_[8]  = npuU_;
  mvaInputCovU2_[9]  = npuPhi_;
  mvaInputCovU2_[10] = puSumEt_;
  mvaInputCovU2_[11] = puMEt_;
  mvaInputCovU2_[12] = puPhi_;
  mvaInputCovU2_[13] = pucSumEt_;
  mvaInputCovU2_[14] = pucU_;
  mvaInputCovU2_[15] = pucPhi_;
  mvaInputCovU2_[16] = jet1Pt_;
  mvaInputCovU2_[17] = jet1Eta_;
  mvaInputCovU2_[18] = jet1Phi_;
  mvaInputCovU2_[19] = jet2Pt_;
  mvaInputCovU2_[20] = jet2Eta_;
  mvaInputCovU2_[21] = jet2Phi_;
  mvaInputCovU2_[22] = numJets_;
  mvaInputCovU2_[23] = numJetsPtGt30_;
  mvaInputCovU2_[24] = pfPhi_ + mvaOutputDPhi_;
  mvaInputCovU2_[25] = mvaOutputU_*pfU_;
  mvaOutputCovU2_    = mvaReaderCovU2_->GetResponse(mvaInputCovU2_)*mvaOutputU_*pfU_;
}

void PFMETAlgorithmMVA::evaluateDPhi ( )

private

Definition at line 286 of file PFMETAlgorithmMVA.cc.

References GBRForest::GetResponse(), jet1Eta_, jet1Phi_, jet1Pt_, jet2Eta_, jet2Phi_, jet2Pt_, mvaInputDPhi_, mvaOutputDPhi_, mvaReaderDPhi_, npuPhi_, npuSumEt_, npuU_, numJets_, numJetsPtGt30_, numVertices_, pfPhi_, pfU_, pucPhi_, pucSumEt_, pucU_, puMEt_, puPhi_, puSumEt_, tkPhi_, tkSumEt_, and tkU_.

Referenced by evaluateMVA().

{ 
  mvaInputDPhi_[0]  = numVertices_;
  mvaInputDPhi_[1]  = pfU_;
  mvaInputDPhi_[2]  = pfPhi_;
  mvaInputDPhi_[3]  = tkSumEt_;
  mvaInputDPhi_[4]  = tkU_;
  mvaInputDPhi_[5]  = tkPhi_;
  mvaInputDPhi_[6]  = npuSumEt_;
  mvaInputDPhi_[7]  = npuU_;
  mvaInputDPhi_[8]  = npuPhi_;
  mvaInputDPhi_[9]  = puSumEt_;
  mvaInputDPhi_[10] = puMEt_;
  mvaInputDPhi_[11] = puPhi_;
  mvaInputDPhi_[12] = pucSumEt_;
  mvaInputDPhi_[13] = pucU_;
  mvaInputDPhi_[14] = pucPhi_;
  mvaInputDPhi_[15] = jet1Pt_;
  mvaInputDPhi_[16] = jet1Eta_;
  mvaInputDPhi_[17] = jet1Phi_;
  mvaInputDPhi_[18] = jet2Pt_;
  mvaInputDPhi_[19] = jet2Eta_;
  mvaInputDPhi_[20] = jet2Phi_;
  mvaInputDPhi_[21] = numJets_;
  mvaInputDPhi_[22] = numJetsPtGt30_;
  mvaOutputDPhi_    = mvaReaderDPhi_->GetResponse(mvaInputDPhi_);
}

void PFMETAlgorithmMVA::evaluateMVA

(

)

Definition at line 220 of file PFMETAlgorithmMVA.cc.

References funct::cos(), evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), hasPhotons_, bookConverter::max, mvaMEt_, mvaMEtCov_, mvaOutputCovU1_, mvaOutputCovU2_, mvaOutputDPhi_, mvaOutputU_, npuU_, pfPhi_, pfU_, colinearityKinematic::Phi, Pi, funct::sin(), mathSSE::sqrt(), sumLeptonPx_, sumLeptonPy_, and tkU_.

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

{
  // CV: MVAs needs to be evaluated in order { DPhi, U1, CovU1, CovU2 }
  //     as MVA for U1 (CovU1, CovU2) uses output of DPhi (DPhi and U1) MVA
  evaluateDPhi();
  evaluateU();
  evaluateCovU1();
  evaluateCovU2();

  // compute MET(Photon check)
  if(hasPhotons_) { 
    //Fix events with unphysical properties
    double sumLeptonPt = std::max(sqrt(sumLeptonPx_*sumLeptonPx_+sumLeptonPy_*sumLeptonPy_),1.);
    if(tkU_/sumLeptonPt < 0.1 || npuU_/sumLeptonPt <  0.1 ) mvaOutputU_      = 1.;
    if(tkU_/sumLeptonPt < 0.1 || npuU_/sumLeptonPt <  0.1 ) mvaOutputDPhi_   = 0.;
  }
  double U      = pfU_*mvaOutputU_;
  double Phi    = pfPhi_ + mvaOutputDPhi_;
  if ( U < 0. ) Phi += Pi;
  double cosPhi = cos(Phi);
  double sinPhi = sin(Phi);
  double metPx  = U*cosPhi - sumLeptonPx_;
  double metPy  = U*sinPhi - sumLeptonPy_;
  double metPt  = sqrt(metPx*metPx + metPy*metPy);
  mvaMEt_.SetCoordinates(metPx, metPy, 0., metPt);

  // compute MET uncertainties in dirrections parallel and perpendicular to hadronic recoil
  // (neglecting uncertainties on lepton momenta)
  mvaMEtCov_(0, 0) =  mvaOutputCovU1_*cosPhi*cosPhi + mvaOutputCovU2_*sinPhi*sinPhi;
  mvaMEtCov_(0, 1) = -mvaOutputCovU1_*sinPhi*cosPhi + mvaOutputCovU2_*sinPhi*cosPhi;
  mvaMEtCov_(1, 0) = mvaMEtCov_(0, 1);
  mvaMEtCov_(1, 1) =  mvaOutputCovU1_*sinPhi*sinPhi + mvaOutputCovU2_*cosPhi*cosPhi;
}

void PFMETAlgorithmMVA::evaluateU ( )

private

Definition at line 256 of file PFMETAlgorithmMVA.cc.

References GBRForest::GetResponse(), jet1Eta_, jet1Phi_, jet1Pt_, jet2Eta_, jet2Phi_, jet2Pt_, mvaInputU_, mvaOutputDPhi_, mvaOutputU_, mvaReaderU_, npuPhi_, npuSumEt_, npuU_, numJets_, numJetsPtGt30_, numVertices_, pfPhi_, pfSumEt_, pfU_, pucPhi_, pucSumEt_, pucU_, puMEt_, puPhi_, puSumEt_, tkPhi_, tkSumEt_, and tkU_.

Referenced by evaluateMVA().

{ 
  mvaInputU_[0]  = pfSumEt_; // PH: helps flattens response vs. Nvtx
  mvaInputU_[1]  = numVertices_;
  mvaInputU_[2]  = pfU_;
  mvaInputU_[3]  = pfPhi_;
  mvaInputU_[4]  = tkSumEt_;
  mvaInputU_[5]  = tkU_;
  mvaInputU_[6]  = tkPhi_;
  mvaInputU_[7]  = npuSumEt_;
  mvaInputU_[8]  = npuU_;
  mvaInputU_[9]  = npuPhi_;
  mvaInputU_[10] = puSumEt_;
  mvaInputU_[11] = puMEt_;
  mvaInputU_[12] = puPhi_;
  mvaInputU_[13] = pucSumEt_;
  mvaInputU_[14] = pucU_;
  mvaInputU_[15] = pucPhi_;
  mvaInputU_[16] = jet1Pt_;
  mvaInputU_[17] = jet1Eta_;
  mvaInputU_[18] = jet1Phi_;
  mvaInputU_[19] = jet2Pt_;
  mvaInputU_[20] = jet2Eta_;
  mvaInputU_[21] = jet2Phi_;
  mvaInputU_[22] = numJets_;
  mvaInputU_[23] = numJetsPtGt30_;
  mvaInputU_[24] = pfPhi_ + mvaOutputDPhi_;
  mvaOutputU_    = mvaReaderU_->GetResponse(mvaInputU_);
}

double PFMETAlgorithmMVA::getCovU1 ( ) const

inline

Definition at line 55 of file PFMETAlgorithmMVA.h.

References mvaOutputCovU1_.

{ return mvaOutputCovU1_; }

double PFMETAlgorithmMVA::getCovU2 ( ) const

inline

Definition at line 56 of file PFMETAlgorithmMVA.h.

References mvaOutputCovU2_.

{ return mvaOutputCovU2_; }

double PFMETAlgorithmMVA::getDPhi ( ) const

inline

Definition at line 54 of file PFMETAlgorithmMVA.h.

References mvaOutputDPhi_.

{ return mvaOutputDPhi_;  }

reco::Candidate::LorentzVector PFMETAlgorithmMVA::getMEt ( ) const

inline

Definition at line 50 of file PFMETAlgorithmMVA.h.

References mvaMEt_.

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

{ return mvaMEt_;    }

const reco::METCovMatrix& PFMETAlgorithmMVA::getMEtCov ( ) const

inline

Definition at line 51 of file PFMETAlgorithmMVA.h.

References mvaMEtCov_.

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

{ return mvaMEtCov_; }

double PFMETAlgorithmMVA::getU ( ) const

inline

Definition at line 53 of file PFMETAlgorithmMVA.h.

References mvaOutputU_.

{ return mvaOutputU_;    }

void PFMETAlgorithmMVA::initialize

(

const edm::EventSetup &

es

)

Definition at line 84 of file PFMETAlgorithmMVA.cc.

References cfg_, edm::ParameterSet::getParameter(), mvaPFMET_cff::loadMVAfromDB, loadMVAfromDB_, mvaNameCovU1_, mvaNameCovU2_, mvaNameDPhi_, mvaNameU_, mvaReaderCovU1_, mvaReaderCovU2_, mvaReaderDPhi_, mvaReaderU_, and AlCaHLTBitMon_QueryRunRegistry::string.

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

{
  if ( loadMVAfromDB_ ) {
    edm::ParameterSet cfgInputRecords = cfg_.getParameter<edm::ParameterSet>("inputRecords");
    mvaNameU_       = cfgInputRecords.getParameter<std::string>("U");
    mvaReaderU_     = loadMVAfromDB(es, mvaNameU_);
    mvaNameDPhi_    = cfgInputRecords.getParameter<std::string>("DPhi");
    mvaReaderDPhi_  = loadMVAfromDB(es, mvaNameDPhi_);
    mvaNameCovU1_   = cfgInputRecords.getParameter<std::string>("CovU1");
    mvaReaderCovU1_ = loadMVAfromDB(es, mvaNameCovU1_);
    mvaNameCovU2_   = cfgInputRecords.getParameter<std::string>("CovU2");
    mvaReaderCovU2_ = loadMVAfromDB(es, mvaNameCovU2_);
  } else {
    edm::ParameterSet cfgInputFileNames = cfg_.getParameter<edm::ParameterSet>("inputFileNames");
    
    mvaNameU_      = "U1Correction";
    mvaNameDPhi_   = "PhiCorrection";
    mvaNameCovU1_  = "CovU1";
    mvaNameCovU2_  = "CovU2";
    
    edm::FileInPath inputFileNameU = cfgInputFileNames.getParameter<edm::FileInPath>("U");
    mvaReaderU_     = loadMVAfromFile(inputFileNameU, mvaNameU_);
    edm::FileInPath inputFileNameDPhi = cfgInputFileNames.getParameter<edm::FileInPath>("DPhi");
    mvaReaderDPhi_  = loadMVAfromFile(inputFileNameDPhi, mvaNameDPhi_);
    edm::FileInPath inputFileNameCovU1 = cfgInputFileNames.getParameter<edm::FileInPath>("CovU1");
    mvaReaderCovU1_ = loadMVAfromFile(inputFileNameCovU1, mvaNameCovU1_);
    edm::FileInPath inputFileNameCovU2 = cfgInputFileNames.getParameter<edm::FileInPath>("CovU2");
    mvaReaderCovU2_ = loadMVAfromFile(inputFileNameCovU2, mvaNameCovU2_);
  }
}

void PFMETAlgorithmMVA::print

(

std::ostream &

stream

)

const

Definition at line 375 of file PFMETAlgorithmMVA.cc.

References jet1Eta_, jet1Phi_, jet1Pt_, jet2Eta_, jet2Phi_, jet2Pt_, mvaOutputCovU1_, mvaOutputCovU2_, mvaOutputDPhi_, mvaOutputU_, npuPhi_, npuSumEt_, npuU_, numJets_, numJetsPtGt30_, numVertices_, pfPhi_, pfSumEt_, pfU_, pucPhi_, pucSumEt_, pucU_, puMEt_, puPhi_, puSumEt_, mathSSE::sqrt(), sumLeptonPx_, sumLeptonPy_, tkPhi_, tkSumEt_, and tkU_.

{
  stream << "<PFMETAlgorithmMVA::print>:" << std::endl;
  stream << " PF: sumEt = " << pfSumEt_ << ", U = " << pfU_ << ", phi = " << pfPhi_ << std::endl;
  stream << " TK: sumEt = " << tkSumEt_ << ", U = " << tkU_ << ", phi = " << tkPhi_ << std::endl;
  stream << " NPU: sumEt = " << npuSumEt_ << ", U = " << npuU_ << ", phi = " << npuPhi_ << std::endl;
  stream << " PU: sumEt = " << puSumEt_ << ", MEt = " << puMEt_ << ", phi = " << puPhi_ << std::endl;
  stream << " PUC: sumEt = " << pucSumEt_ << ", U = " << pucU_ << ", phi = " << pucPhi_ << std::endl;
  stream << " jet1: Pt = " << jet1Pt_ << ", eta = " << jet1Eta_ << ", phi = " << jet1Phi_ << std::endl;
  stream << " jet2: Pt = " << jet2Pt_ << ", eta = " << jet2Eta_ << ", phi = " << jet2Phi_ << std::endl;
  stream << " num. jets = " << numJets_ << " (" << numJetsPtGt30_ << " with Pt > 30 GeV)" << std::endl;
  stream << " num. vertices = " << numVertices_ << std::endl;
  stream << " MVA output: U = " << mvaOutputU_ << ", dPhi = " << mvaOutputDPhi_ << "," 
     << " covU1 = " << mvaOutputCovU1_ << ", covU2 = " << mvaOutputCovU2_ << std::endl;
  stream << " sum(leptons): Pt = " << sqrt(sumLeptonPx_*sumLeptonPx_ + sumLeptonPy_*sumLeptonPy_) << ","
     << " phi = " << atan2(sumLeptonPy_, sumLeptonPx_) << " "
     << "(Px = " << sumLeptonPx_ << ", Py = " << sumLeptonPy_ << ")" << std::endl;
}

void PFMETAlgorithmMVA::setHasPhotons ( bool  hasPhotons )

inline

Definition at line 41 of file PFMETAlgorithmMVA.h.

References hasPhotons_.

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

{ hasPhotons_ = hasPhotons; }

void PFMETAlgorithmMVA::setInput	(	const std::vector< reco::PUSubMETCandInfo > &	leptons,
		const std::vector< reco::PUSubMETCandInfo > &	jets,
		const std::vector< reco::PUSubMETCandInfo > &	pfCandidates,
		const std::vector< reco::Vertex::Point > &	vertices
	)

Definition at line 116 of file PFMETAlgorithmMVA.cc.

References chargedSumLeptonPx_, chargedSumLeptonPy_, MvaMEtUtilities::computeAllSums(), MvaMEtUtilities::computeRecoil(), MvaMEtUtilities::getCleanedJets(), MvaMEtUtilities::getLeptonsChSumMEX(), MvaMEtUtilities::getLeptonsChSumMEY(), MvaMEtUtilities::getLeptonsSumMEX(), MvaMEtUtilities::getLeptonsSumMEY(), MvaMEtUtilities::kChHS, MvaMEtUtilities::kHS, MvaMEtUtilities::kHSMinusNeutralPU, MvaMEtUtilities::kPF, MvaMEtUtilities::kPU, MvaMEtUtilities::leadJetP4(), CommonMETData::met, MvaMEtUtilities::numJetsAboveThreshold(), CommonMETData::phi, MvaMEtUtilities::subleadJetP4(), CommonMETData::sumet, sumLeptonPx_, sumLeptonPy_, and utils_.

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

{

  
  utils_.computeAllSums( jets, leptons, pfCandidates);
  
  sumLeptonPx_        = utils_.getLeptonsSumMEX();
  sumLeptonPy_        = utils_.getLeptonsSumMEY();

  chargedSumLeptonPx_ = utils_.getLeptonsChSumMEX();
  chargedSumLeptonPy_ = utils_.getLeptonsChSumMEY();

  const std::vector<reco::PUSubMETCandInfo> jets_cleaned = utils_.getCleanedJets();

  CommonMETData pfRecoil_data  = utils_.computeRecoil( MvaMEtUtilities::kPF );
  CommonMETData chHSRecoil_data  = utils_.computeRecoil( MvaMEtUtilities::kChHS );
  CommonMETData hsRecoil_data = utils_.computeRecoil( MvaMEtUtilities::kHS );
  CommonMETData puRecoil_data = utils_.computeRecoil( MvaMEtUtilities::kPU );
  CommonMETData hsMinusNeutralPUMEt_data = utils_.computeRecoil( MvaMEtUtilities::kHSMinusNeutralPU );

  reco::Candidate::LorentzVector jet1P4 = utils_.leadJetP4(jets_cleaned);
  reco::Candidate::LorentzVector jet2P4 = utils_.subleadJetP4(jets_cleaned);

  double pfSumEt       = pfRecoil_data.sumet;
  double pfU           = pfRecoil_data.met;
  double pfPhi         = pfRecoil_data.phi;
  double tkSumEt       = chHSRecoil_data.sumet;
  double tkU           = chHSRecoil_data.met;
  double tkPhi         = chHSRecoil_data.phi;
  double npuSumEt      = hsRecoil_data.sumet;
  double npuU          = hsRecoil_data.met; 
  double npuPhi        = hsRecoil_data.phi;
  double puSumEt       = puRecoil_data.sumet;
  double puMEt         = puRecoil_data.met;
  double puPhi         = puRecoil_data.phi;
  double pucSumEt      = hsMinusNeutralPUMEt_data.sumet; 
  double pucU          = hsMinusNeutralPUMEt_data.met; 
  double pucPhi        = hsMinusNeutralPUMEt_data.phi;
  double jet1Pt        = jet1P4.pt();
  double jet1Eta       = jet1P4.eta();
  double jet1Phi       = jet1P4.phi();
  double jet2Pt        = jet2P4.pt();
  double jet2Eta       = jet2P4.eta();
  double jet2Phi       = jet2P4.phi();
  double numJetsPtGt30 = utils_.numJetsAboveThreshold(jets_cleaned, 30.);
  double numJets       = jets_cleaned.size();
  double numVertices   = vertices.size();

  setInput(pfSumEt, pfU, pfPhi,
       tkSumEt, tkU, tkPhi,
       npuSumEt, npuU, npuPhi,
       puSumEt, puMEt, puPhi,
       pucSumEt, pucU, pucPhi,
       jet1Pt, jet1Eta, jet1Phi,
       jet2Pt, jet2Eta, jet2Phi,
       numJetsPtGt30, numJets, 
       numVertices);

}

void PFMETAlgorithmMVA::setInput ( double  pfSumEt, double  pfU, double  pfPhi, double  tkSumEt, double  tkU, double  tkPhi, double  npuSumEt, double  npuU, double  npuPhi, double  puSumEt, double  puMEt, double  puPhi, double  pucSumEt, double  pucU, double  pucPhi, double  jet1Pt, double  jet1Eta, double  jet1Phi, double  jet2Pt, double  jet2Eta, double  jet2Phi, double  numJetsPtGt30, double  numJets, double  numVertices  )

private

Definition at line 179 of file PFMETAlgorithmMVA.cc.

References jet1Eta_, jet1Phi_, jet1Pt_, jet2Eta_, jet2Phi_, jet2Pt_, npuPhi_, npuSumEt_, npuU_, numJets_, numJetsPtGt30_, numVertices_, pfPhi_, pfSumEt_, pfU_, pucPhi_, pucSumEt_, pucU_, puMEt_, puPhi_, puSumEt_, tkPhi_, tkSumEt_, and tkU_.

{
  // protection against "empty events"
  if ( pfSumEt < 1. ) pfSumEt = 1.;
  
  pfSumEt_       = pfSumEt;
  pfU_           = pfU;
  pfPhi_         = pfPhi;
  tkSumEt_       = tkSumEt/pfSumEt_;
  tkU_           = tkU;
  tkPhi_         = tkPhi;
  npuSumEt_      = npuSumEt/pfSumEt_;
  npuU_          = npuU;
  npuPhi_        = npuPhi;
  puSumEt_       = puSumEt/pfSumEt_;
  puMEt_         = puMEt;
  puPhi_         = puPhi;
  pucSumEt_      = pucSumEt/pfSumEt_;
  pucU_          = pucU;
  pucPhi_        = pucPhi;
  jet1Pt_        = jet1Pt;
  jet1Eta_       = jet1Eta;
  jet1Phi_       = jet1Phi;
  jet2Pt_        = jet2Pt;
  jet2Eta_       = jet2Eta;
  jet2Phi_       = jet2Phi;
  numJetsPtGt30_ = numJetsPtGt30;
  numJets_       = numJets;
  numVertices_   = numVertices;
}

Member Data Documentation

edm::ParameterSet PFMETAlgorithmMVA::cfg_

private

Definition at line 138 of file PFMETAlgorithmMVA.h.

Referenced by initialize().

double PFMETAlgorithmMVA::chargedSumLeptonPx_

private

Definition at line 124 of file PFMETAlgorithmMVA.h.

Referenced by setInput().

double PFMETAlgorithmMVA::chargedSumLeptonPy_

private

Definition at line 125 of file PFMETAlgorithmMVA.h.

Referenced by setInput().

bool PFMETAlgorithmMVA::hasPhotons_

private

Definition at line 85 of file PFMETAlgorithmMVA.h.

Referenced by evaluateMVA(), and setHasPhotons().

Float_t PFMETAlgorithmMVA::jet1Eta_

private

Definition at line 103 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::jet1Phi_

private

Definition at line 104 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::jet1Pt_

private

Definition at line 102 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::jet2Eta_

private

Definition at line 106 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::jet2Phi_

private

Definition at line 107 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::jet2Pt_

private

Definition at line 105 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

bool PFMETAlgorithmMVA::loadMVAfromDB_

private

Definition at line 136 of file PFMETAlgorithmMVA.h.

Referenced by initialize(), PFMETAlgorithmMVA(), and ~PFMETAlgorithmMVA().

Float_t* PFMETAlgorithmMVA::mvaInputCovU1_

private

Definition at line 114 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), PFMETAlgorithmMVA(), and ~PFMETAlgorithmMVA().

Float_t* PFMETAlgorithmMVA::mvaInputCovU2_

private

Definition at line 115 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU2(), PFMETAlgorithmMVA(), and ~PFMETAlgorithmMVA().

Float_t* PFMETAlgorithmMVA::mvaInputDPhi_

private

Definition at line 113 of file PFMETAlgorithmMVA.h.

Referenced by evaluateDPhi(), PFMETAlgorithmMVA(), and ~PFMETAlgorithmMVA().

Float_t* PFMETAlgorithmMVA::mvaInputU_

private

Definition at line 112 of file PFMETAlgorithmMVA.h.

Referenced by evaluateU(), PFMETAlgorithmMVA(), and ~PFMETAlgorithmMVA().

reco::Candidate::LorentzVector PFMETAlgorithmMVA::mvaMEt_

private

Definition at line 127 of file PFMETAlgorithmMVA.h.

Referenced by evaluateMVA(), and getMEt().

reco::METCovMatrix PFMETAlgorithmMVA::mvaMEtCov_

private

Definition at line 129 of file PFMETAlgorithmMVA.h.

Referenced by evaluateMVA(), and getMEtCov().

std::string PFMETAlgorithmMVA::mvaNameCovU1_

private

Definition at line 81 of file PFMETAlgorithmMVA.h.

Referenced by initialize().

std::string PFMETAlgorithmMVA::mvaNameCovU2_

private

Definition at line 82 of file PFMETAlgorithmMVA.h.

Referenced by initialize().

std::string PFMETAlgorithmMVA::mvaNameDPhi_

private

Definition at line 80 of file PFMETAlgorithmMVA.h.

Referenced by initialize().

std::string PFMETAlgorithmMVA::mvaNameU_

private

Definition at line 79 of file PFMETAlgorithmMVA.h.

Referenced by initialize().

Float_t PFMETAlgorithmMVA::mvaOutputCovU1_

private

Definition at line 119 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateMVA(), getCovU1(), and print().

Float_t PFMETAlgorithmMVA::mvaOutputCovU2_

private

Definition at line 120 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU2(), evaluateMVA(), getCovU2(), and print().

Float_t PFMETAlgorithmMVA::mvaOutputDPhi_

private

Definition at line 118 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateMVA(), evaluateU(), getDPhi(), and print().

Float_t PFMETAlgorithmMVA::mvaOutputU_

private

Definition at line 117 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateMVA(), evaluateU(), getU(), and print().

const GBRForest* PFMETAlgorithmMVA::mvaReaderCovU1_

private

Definition at line 133 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), initialize(), and ~PFMETAlgorithmMVA().

const GBRForest* PFMETAlgorithmMVA::mvaReaderCovU2_

private

Definition at line 134 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU2(), initialize(), and ~PFMETAlgorithmMVA().

const GBRForest* PFMETAlgorithmMVA::mvaReaderDPhi_

private

Definition at line 132 of file PFMETAlgorithmMVA.h.

Referenced by evaluateDPhi(), initialize(), and ~PFMETAlgorithmMVA().

const GBRForest* PFMETAlgorithmMVA::mvaReaderU_

private

Definition at line 131 of file PFMETAlgorithmMVA.h.

Referenced by evaluateU(), initialize(), and ~PFMETAlgorithmMVA().

int PFMETAlgorithmMVA::mvaType_

private

Definition at line 84 of file PFMETAlgorithmMVA.h.

Referenced by PFMETAlgorithmMVA().

Float_t PFMETAlgorithmMVA::npuPhi_

private

Definition at line 95 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::npuSumEt_

private

Definition at line 93 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::npuU_

private

Definition at line 94 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateMVA(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::numJets_

private

Definition at line 109 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::numJetsPtGt30_

private

Definition at line 108 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::numVertices_

private

Definition at line 110 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::pfPhi_

private

Definition at line 89 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateMVA(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::pfSumEt_

private

Definition at line 87 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::pfU_

private

Definition at line 88 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateMVA(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::pucPhi_

private

Definition at line 101 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::pucSumEt_

private

Definition at line 99 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::pucU_

private

Definition at line 100 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::puMEt_

private

Definition at line 97 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::puPhi_

private

Definition at line 98 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::puSumEt_

private

Definition at line 96 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

double PFMETAlgorithmMVA::sumLeptonPx_

private

Definition at line 122 of file PFMETAlgorithmMVA.h.

Referenced by evaluateMVA(), print(), and setInput().

double PFMETAlgorithmMVA::sumLeptonPy_

private

Definition at line 123 of file PFMETAlgorithmMVA.h.

Referenced by evaluateMVA(), print(), and setInput().

Float_t PFMETAlgorithmMVA::tkPhi_

private

Definition at line 92 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::tkSumEt_

private

Definition at line 90 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateU(), print(), and setInput().

Float_t PFMETAlgorithmMVA::tkU_

private

Definition at line 91 of file PFMETAlgorithmMVA.h.

Referenced by evaluateCovU1(), evaluateCovU2(), evaluateDPhi(), evaluateMVA(), evaluateU(), print(), and setInput().

MvaMEtUtilities PFMETAlgorithmMVA::utils_

private

Definition at line 77 of file PFMETAlgorithmMVA.h.

Referenced by setInput().