240 std::cout << __FUNCTION__ <<
":" << __LINE__ <<
":pvHandle handle not found!" << std::endl;
243 const int nvtx = pvHandle->size();
271 met = caloMETs->front();
274 met = pfMETs->front();
277 met = genMETs->front();
280 met = patMET->front();
283 const double SumET =
met.sumEt();
284 const double METSig =
met.mEtSig();
285 const double MET =
met.pt();
286 const double MEx =
met.px();
287 const double MEy =
met.py();
288 const double METPhi =
met.phi();
300 bool isvalidgenmet =
false;
306 isvalidgenmet =
true;
312 isvalidgenmet =
true;
326 else if (
MET > 20 &&
MET < 40)
328 else if (
MET > 40 &&
MET < 60)
330 else if (
MET > 60 &&
MET < 80)
332 else if (
MET > 80 &&
MET < 100)
334 else if (
MET > 100 &&
MET < 150)
336 else if (
MET > 150 &&
MET < 200)
338 else if (
MET > 200 &&
MET < 300)
340 else if (
MET > 300 &&
MET < 400)
342 else if (
MET > 400 &&
MET < 500)
348 edm::LogInfo(
"OutputInfo") <<
" failed to retrieve data required by MET Task: genMetTrue";
363 edm::LogInfo(
"OutputInfo") <<
" failed to retrieve data required by MET Task: genMetCalo";
373 const double caloEmEtFraction = calomet->
emEtFraction();
374 const double caloHadEtInHB = calomet->
hadEtInHB();
375 const double caloHadEtInHO = calomet->
hadEtInHO();
376 const double caloHadEtInHE = calomet->
hadEtInHE();
377 const double caloHadEtInHF = calomet->
hadEtInHF();
378 const double caloEmEtInEB = calomet->
emEtInEB();
379 const double caloEmEtInEE = calomet->
emEtInEE();
380 const double caloEmEtInHF = calomet->
emEtInHF();
401 genmet = &(genMETs->front());
403 const double NeutralEMEtFraction =
genmet->NeutralEMEtFraction();
404 const double NeutralHadEtFraction =
genmet->NeutralHadEtFraction();
405 const double ChargedEMEtFraction =
genmet->ChargedEMEtFraction();
406 const double ChargedHadEtFraction =
genmet->ChargedHadEtFraction();
407 const double MuonEtFraction =
genmet->MuonEtFraction();
408 const double InvisibleEtFraction =
genmet->InvisibleEtFraction();
436 const pat::MET *patmet = &(patMET->front());
Analysis-level MET class.
double Type6EtFraction() const
double maxEtInEmTowers() const
MonitorElement * mChargedEMEtFraction
MonitorElement * mMuonEtFraction
double neutralHadronEt() const
MonitorElement * mMETDifference_GenMETTrue
MonitorElement * mMETUnc_MuonEnUp
bool isPFMET() const
True if this pat::MET was made from a reco::pfMET.
MonitorElement * mMETUnc_ElectronEnDown
double chargedHadronEt() const
double chargedHadronEtFraction() const
MonitorElement * mMETUnc_ElectronEnUp
MonitorElement * mMETDifference_GenMETTrue_MET400to500
MonitorElement * mMETUnc_TauEnDown
double ChargedEMEtFraction() const
MonitorElement * mPFelectronEtFraction
MonitorElement * mCaloMaxEtInEmTowers
MonitorElement * mCaloSETInpHF
genmet
do same things for GEN
MonitorElement * mCaloHadEtInHF
T const * product() const
double neutralHadronEtFraction() const
double etFractionHadronic() const
double HFHadronEtFraction() const
MonitorElement * mPFphotonEtFraction
double HFHadronEt() const
double electronEt() const
MonitorElement * mPFHFEMEt
genMetCalo
____________________________________________________________________________||
double Type7EtFraction() const
MonitorElement * mMETDifference_GenMETTrue_MET40to60
MonitorElement * mCaloEmEtFraction
MonitorElement * mCaloEmEtInEE
MonitorElement * mCaloSETInmHF
MonitorElement * mNeutralEMEtFraction
MonitorElement * mMETDifference_GenMETTrue_MET20to40
MonitorElement * mMETDifference_GenMETTrue_MET60to80
MonitorElement * mCaloMaxEtInHadTowers
MonitorElement * mInvisibleEtFraction
MonitorElement * mCaloEtFractionHadronic
double CaloSETInpHF() const
MonitorElement * mPFmuonEt
MonitorElement * mMETUnc_UnclusteredEnUp
MonitorElement * mPFelectronEt
double HFEMEtFraction() const
double emEtFraction() const
MonitorElement * mPFphotonEt
MonitorElement * mCaloEmEtInEB
MonitorElement * mMETUnc_JetEnUp
MonitorElement * mCaloHadEtInHE
edm::EDGetTokenT< reco::PFMETCollection > pfMETsToken_
edm::EDGetTokenT< reco::CaloMETCollection > caloMETsToken_
MonitorElement * mPFneutralHadronEtFraction
MonitorElement * mMETUnc_JetResUp
MonitorElement * mPFHFHadronEt
double maxEtInHadTowers() const
edm::EDGetTokenT< std::vector< reco::Vertex > > pvToken_
edm::EDGetTokenT< reco::GenMETCollection > genMETsToken_
MonitorElement * mMETUnc_JetEnDown
MonitorElement * mMETDifference_GenMETTrue_MET150to200
MonitorElement * mMETDeltaPhi_GenMETCalo
MonitorElement * mMETDifference_GenMETTrue_MET100to150
MonitorElement * mCaloHadEtInHB
MonitorElement * mMETUnc_PhotonEnDown
MonitorElement * mMETDifference_GenMETTrue_MET300to400
MonitorElement * mMETDifference_GenMETCalo
MonitorElement * mMETUnc_JetResDown
MonitorElement * mCaloHadEtInHO
Log< level::Info, false > LogInfo
MonitorElement * mMETDifference_GenMETTrue_MET500
MonitorElement * mPFchargedHadronEt
MonitorElement * mChargedHadEtFraction
MonitorElement * mMETDeltaPhi_GenMETTrue
MonitorElement * mPFneutralHadronEt
MonitorElement * mMETUnc_TauEnUp
edm::EDGetTokenT< reco::GenMETCollection > genMETsCaloToken_
double NeutralHadEtFraction() const
edm::EDGetTokenT< reco::GenMETCollection > genMETsTrueToken_
MonitorElement * mPFmuonEtFraction
MonitorElement * mMETFine
MonitorElement * mMETUnc_MuonEnDown
double electronEtFraction() const
MonitorElement * mNeutralHadEtFraction
edm::EDGetTokenT< pat::METCollection > patMETToken_
double ChargedHadEtFraction() const
double photonEtFraction() const
MonitorElement * mPFHFEMEtFraction
MonitorElement * mMETUnc_UnclusteredEnDown
genMetTrue
____________________________________________________________________________||
MonitorElement * mMETUnc_PhotonEnUp
MonitorElement * mNvertex
MonitorElement * mMET_Nvtx
MonitorElement * mPFchargedHadronEtFraction
MonitorElement * mMETDifference_GenMETTrue_MET0to20
double NeutralEMFraction() const
double CaloSETInmHF() const
MonitorElement * mMETDifference_GenMETTrue_MET200to300
double muonEtFraction() const
double shiftedPt(METUncertainty shift, METCorrectionLevel level=Type1) const
MonitorElement * mMETDifference_GenMETTrue_MET80to100
MonitorElement * mCaloEmEtInHF
double MuonEtFraction() const
MonitorElement * mPFHFHadronEtFraction