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< mvaMEtUtilities::leptonInfo > &, const std::vector< mvaMEtUtilities::JetInfo > &, const std::vector< mvaMEtUtilities::pfCandInfo > &, const std::vector< reco::Vertex::Point > &)
	~PFMETAlgorithmMVA ()

Protected 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)

Protected Attributes
edm::ParameterSet	cfg_
double	chargedSumLeptonPx_
double	chargedSumLeptonPy_
double	dZcut_
bool	hasPhotons_
bool	is42_
bool	isOld42_
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 dZcut_, edm::ParameterSet::getParameter(), is42_, isOld42_, 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;

  dZcut_ = cfg.getParameter<double>("dZcut");
  isOld42_ = cfg.getParameter<bool>("useOld42");
  
  loadMVAfromDB_ = cfg.getParameter<bool>("loadMVAfromDB");
  
  is42_ = cfg.getParameter<bool>("is42");
  
  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 74 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 ( )

protected

Definition at line 326 of file PFMETAlgorithmMVA.cc.

References GBRForest::GetResponse(), isOld42_, 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_;
  if ( !isOld42_ ) mvaOutputCovU1_ *= mvaOutputCovU1_; // PH: Training is not on the square anymore
}

void PFMETAlgorithmMVA::evaluateCovU2 ( )

protected

Definition at line 358 of file PFMETAlgorithmMVA.cc.

References GBRForest::GetResponse(), isOld42_, 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_;
  if ( !isOld42_ ) mvaOutputCovU2_ *= mvaOutputCovU2_; // PH: Training is not on the square anymore
}

void PFMETAlgorithmMVA::evaluateDPhi ( )

protected

Definition at line 298 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 232 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_; // CV: U is actually minus the hadronic recoil in the event
  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 ( )

protected

Definition at line 268 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 89 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_ ) {
    //std::cout << "<PFMETAlgorithmMVA::initialize>: loading BDTs from Database." << std::endl;
    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 {
    //std::cout << "<PFMETAlgorithmMVA::initialize>: loading BDTs from ROOT files." << std::endl;
    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 389 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< mvaMEtUtilities::leptonInfo > &	leptons,
		const std::vector< mvaMEtUtilities::JetInfo > &	jets,
		const std::vector< mvaMEtUtilities::pfCandInfo > &	pfCandidates,
		const std::vector< reco::Vertex::Point > &	vertices
	)

CommonMETData sumLeptons = utils_.computePFCandSum(pfCandidates_leptons, dZcut_, 2);

Definition at line 123 of file PFMETAlgorithmMVA.cc.

References chargedSumLeptonPx_, chargedSumLeptonPy_, mvaMEtUtilities::cleanJets(), mvaMEtUtilities::computeNegNoPURecoil(), mvaMEtUtilities::computeNegPFRecoil(), mvaMEtUtilities::computeNegPUCRecoil(), mvaMEtUtilities::computeNegTrackRecoil(), mvaMEtUtilities::computePUMEt(), mvaMEtUtilities::computeSumLeptons(), dZcut_, is42_, mvaMEtUtilities::leadJetP4(), CommonMETData::met, CommonMETData::mex, CommonMETData::mey, mvaMEtUtilities::numJetsAboveThreshold(), CommonMETData::phi, mvaMEtUtilities::subleadJetP4(), CommonMETData::sumet, sumLeptonPx_, sumLeptonPy_, and utils_.

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

{
  //std::vector<mvaMEtUtilities::pfCandInfo> pfCandidates_leptons = utils_.cleanPFCands(pfCandidates, leptons, 0.3, true);
  //const std::vector<mvaMEtUtilities::pfCandInfo> pfCandidates_cleaned = utils_.cleanPFCands(pfCandidates, leptons, 0.3, false);

  CommonMETData        sumLeptons = utils_.computeSumLeptons(leptons, false);
  CommonMETData chargedSumLeptons = utils_.computeSumLeptons(leptons, true);
  
  sumLeptonPx_        = sumLeptons.mex;
  sumLeptonPy_        = sumLeptons.mey;

  chargedSumLeptonPx_ = chargedSumLeptons.mex;
  chargedSumLeptonPy_ = chargedSumLeptons.mey;

  double ptThreshold = -1000.;
  if ( is42_ ) ptThreshold = 1.;  //PH: For 42 training added a pT cut of 1 GeV on corrected Jets
  std::vector<mvaMEtUtilities::JetInfo> jets_cleaned = utils_.cleanJets(jets, leptons, ptThreshold, 0.5);

  CommonMETData pfRecoil_data  = utils_.computeNegPFRecoil   (sumLeptons       , pfCandidates,               dZcut_);
  CommonMETData tkRecoil_data  = utils_.computeNegTrackRecoil(chargedSumLeptons, pfCandidates,               dZcut_);
  CommonMETData npuRecoil_data = utils_.computeNegNoPURecoil (chargedSumLeptons, pfCandidates, jets_cleaned, dZcut_);
  CommonMETData pucRecoil_data = utils_.computeNegPUCRecoil  (sumLeptons       , pfCandidates, jets_cleaned, dZcut_);
  CommonMETData puMEt_data     = utils_.computePUMEt         (                   pfCandidates, jets_cleaned, 0.2); //dZCut bug

  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       = tkRecoil_data.sumet;
  double tkU           = tkRecoil_data.met;
  double tkPhi         = tkRecoil_data.phi;
  double npuSumEt      = npuRecoil_data.sumet;
  double npuU          = npuRecoil_data.met; 
  double npuPhi        = npuRecoil_data.phi;
  double puSumEt       = puMEt_data.sumet;
  double puMEt         = puMEt_data.met;
  double puPhi         = puMEt_data.phi;
  double pucSumEt      = pucRecoil_data.sumet; 
  double pucU          = pucRecoil_data.met; 
  double pucPhi        = pucRecoil_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  )

protected

Definition at line 191 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_.

{
  // CV: add 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_

protected

Definition at line 141 of file PFMETAlgorithmMVA.h.

Referenced by initialize().

double PFMETAlgorithmMVA::chargedSumLeptonPx_

protected

Definition at line 127 of file PFMETAlgorithmMVA.h.

Referenced by setInput().

double PFMETAlgorithmMVA::chargedSumLeptonPy_

protected

Definition at line 128 of file PFMETAlgorithmMVA.h.

Referenced by setInput().

double PFMETAlgorithmMVA::dZcut_

protected

Definition at line 88 of file PFMETAlgorithmMVA.h.

Referenced by PFMETAlgorithmMVA(), and setInput().

bool PFMETAlgorithmMVA::hasPhotons_

protected

Definition at line 87 of file PFMETAlgorithmMVA.h.

Referenced by evaluateMVA(), and setHasPhotons().

bool PFMETAlgorithmMVA::is42_

protected

Definition at line 85 of file PFMETAlgorithmMVA.h.

Referenced by PFMETAlgorithmMVA(), and setInput().

bool PFMETAlgorithmMVA::isOld42_

protected

Definition at line 86 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::jet1Eta_

protected

Definition at line 106 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::jet1Phi_

protected

Definition at line 107 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::jet1Pt_

protected

Definition at line 105 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::jet2Eta_

protected

Definition at line 109 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::jet2Phi_

protected

Definition at line 110 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::jet2Pt_

protected

Definition at line 108 of file PFMETAlgorithmMVA.h.

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

bool PFMETAlgorithmMVA::loadMVAfromDB_

protected

Definition at line 139 of file PFMETAlgorithmMVA.h.

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

Float_t* PFMETAlgorithmMVA::mvaInputCovU1_

protected

Definition at line 117 of file PFMETAlgorithmMVA.h.

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

Float_t* PFMETAlgorithmMVA::mvaInputCovU2_

protected

Definition at line 118 of file PFMETAlgorithmMVA.h.

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

Float_t* PFMETAlgorithmMVA::mvaInputDPhi_

protected

Definition at line 116 of file PFMETAlgorithmMVA.h.

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

Float_t* PFMETAlgorithmMVA::mvaInputU_

protected

Definition at line 115 of file PFMETAlgorithmMVA.h.

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

reco::Candidate::LorentzVector PFMETAlgorithmMVA::mvaMEt_

protected

Definition at line 130 of file PFMETAlgorithmMVA.h.

Referenced by evaluateMVA(), and getMEt().

reco::METCovMatrix PFMETAlgorithmMVA::mvaMEtCov_

protected

Definition at line 132 of file PFMETAlgorithmMVA.h.

Referenced by evaluateMVA(), and getMEtCov().

std::string PFMETAlgorithmMVA::mvaNameCovU1_

protected

Definition at line 81 of file PFMETAlgorithmMVA.h.

Referenced by initialize().

std::string PFMETAlgorithmMVA::mvaNameCovU2_

protected

Definition at line 82 of file PFMETAlgorithmMVA.h.

Referenced by initialize().

std::string PFMETAlgorithmMVA::mvaNameDPhi_

protected

Definition at line 80 of file PFMETAlgorithmMVA.h.

Referenced by initialize().

std::string PFMETAlgorithmMVA::mvaNameU_

protected

Definition at line 79 of file PFMETAlgorithmMVA.h.

Referenced by initialize().

Float_t PFMETAlgorithmMVA::mvaOutputCovU1_

protected

Definition at line 122 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::mvaOutputCovU2_

protected

Definition at line 123 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::mvaOutputDPhi_

protected

Definition at line 121 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::mvaOutputU_

protected

Definition at line 120 of file PFMETAlgorithmMVA.h.

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

const GBRForest* PFMETAlgorithmMVA::mvaReaderCovU1_

protected

Definition at line 136 of file PFMETAlgorithmMVA.h.

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

const GBRForest* PFMETAlgorithmMVA::mvaReaderCovU2_

protected

Definition at line 137 of file PFMETAlgorithmMVA.h.

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

const GBRForest* PFMETAlgorithmMVA::mvaReaderDPhi_

protected

Definition at line 135 of file PFMETAlgorithmMVA.h.

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

const GBRForest* PFMETAlgorithmMVA::mvaReaderU_

protected

Definition at line 134 of file PFMETAlgorithmMVA.h.

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

int PFMETAlgorithmMVA::mvaType_

protected

Definition at line 84 of file PFMETAlgorithmMVA.h.

Referenced by PFMETAlgorithmMVA().

Float_t PFMETAlgorithmMVA::npuPhi_

protected

Definition at line 98 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::npuSumEt_

protected

Definition at line 96 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::npuU_

protected

Definition at line 97 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::numJets_

protected

Definition at line 112 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::numJetsPtGt30_

protected

Definition at line 111 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::numVertices_

protected

Definition at line 113 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::pfPhi_

protected

Definition at line 92 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::pfSumEt_

protected

Definition at line 90 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::pfU_

protected

Definition at line 91 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::pucPhi_

protected

Definition at line 104 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::pucSumEt_

protected

Definition at line 102 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::pucU_

protected

Definition at line 103 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::puMEt_

protected

Definition at line 100 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::puPhi_

protected

Definition at line 101 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::puSumEt_

protected

Definition at line 99 of file PFMETAlgorithmMVA.h.

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

double PFMETAlgorithmMVA::sumLeptonPx_

protected

Definition at line 125 of file PFMETAlgorithmMVA.h.

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

double PFMETAlgorithmMVA::sumLeptonPy_

protected

Definition at line 126 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::tkPhi_

protected

Definition at line 95 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::tkSumEt_

protected

Definition at line 93 of file PFMETAlgorithmMVA.h.

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

Float_t PFMETAlgorithmMVA::tkU_

protected

Definition at line 94 of file PFMETAlgorithmMVA.h.

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

mvaMEtUtilities PFMETAlgorithmMVA::utils_

protected

Definition at line 77 of file PFMETAlgorithmMVA.h.

Referenced by setInput().