46 m_bitAlgTechTrig_ = -1;
48 miniaodfilterdec = -1;
58 triggerResultsToken_ = consumes<edm::TriggerResults>(
edm::InputTag(triggerResultsLabel_));
63 isMiniAODMet_ = (
std::string(
"miniaod") == MetType_);
73 jetID_ValueMapToken_ = consumes<edm::ValueMap<reco::JetID> >(inputJetIDValueMap);
87 ptThreshold_ =
parameters.getParameter<
double>(
"ptThreshold");
90 pfMetToken_ = consumes<reco::PFMETCollection>(
edm::InputTag(metCollectionLabel_));
93 caloMetToken_ = consumes<reco::CaloMETCollection>(
edm::InputTag(metCollectionLabel_));
96 patMetToken_ = consumes<pat::METCollection>(
edm::InputTag(metCollectionLabel_));
102 fill_met_high_level_histo =
parameters.getParameter<
bool>(
"fillMetHighLevel");
103 fillCandidateMap_histos =
parameters.getParameter<
bool>(
"fillCandidateMaps");
105 hTriggerLabelsIsSet_ =
false;
109 diagnosticsParameters_ = pSet.
getParameter<std::vector<edm::ParameterSet> >(
"METDiagonisticsParameters");
116 bypassAllPVChecks_ = cleaningParameters_.getParameter<
bool>(
"bypassAllPVChecks");
117 bypassAllDCSChecks_ = cleaningParameters_.getParameter<
bool>(
"bypassAllDCSChecks");
118 runcosmics_ =
parameters.getUntrackedParameter<
bool>(
"runcosmics");
119 onlyCleaned_ =
parameters.getUntrackedParameter<
bool>(
"onlyCleaned");
120 vertexTag_ = cleaningParameters_.getParameter<
edm::InputTag>(
"vertexCollection");
121 vertexToken_ = consumes<std::vector<reco::Vertex> >(
edm::InputTag(vertexTag_));
124 gtTag_ = cleaningParameters_.getParameter<
edm::InputTag>(
"gtLabel");
125 gtToken_ = consumes<L1GlobalTriggerReadoutRecord>(
edm::InputTag(gtTag_));
130 caloJetsToken_ = consumes<reco::CaloJetCollection>(jetCollectionLabel_);
133 pfJetsToken_ = consumes<reco::PFJetCollection>(jetCollectionLabel_);
135 patJetsToken_ = consumes<pat::JetCollection>(jetCollectionLabel_);
141 hbheNoiseFilterResultToken_ = consumes<bool>(hbheNoiseFilterResultTag_);
143 hbheIsoNoiseFilterResultToken_ = consumes<bool>(hbheNoiseIsoFilterResultTag_);
145 CSCHaloResultToken_ = consumes<bool>(CSCHaloResultTag_);
147 CSCHalo2015ResultToken_ = consumes<bool>(CSCHalo2015ResultTag_);
149 EcalDeadCellTriggerToken_ = consumes<bool>(EcalDeadCellTriggerTag_);
151 EcalDeadCellBoundaryToken_ = consumes<bool>(EcalDeadCellBoundaryTag_);
153 eeBadScFilterToken_ = consumes<bool>(eeBadScFilterTag_);
155 HcalStripHaloToken_ = consumes<bool>(HcalStripHaloTag_);
159 METFilterMiniAODToken_ = consumes<edm::TriggerResults>(METFilterMiniAODLabel_);
162 METFilterMiniAODToken2_ = consumes<edm::TriggerResults>(METFilterMiniAODLabel2_);
166 nbinsPV_ =
parameters.getParameter<
int>(
"pVBin");
167 nPVMin_ =
parameters.getParameter<
double>(
"pVMin");
168 nPVMax_ =
parameters.getParameter<
double>(
"pVMax");
171 for (edm::VParameterSet::const_iterator it = triggerSelectedSubFolders_.begin();
172 it != triggerSelectedSubFolders_.end();
175 triggerFolderExpr_.push_back(it->getParameter<std::vector<std::string> >(
"hltPaths"));
176 triggerFolderLabels_.push_back(it->getParameter<
std::string>(
"label"));
188 for (std::vector<GenericTriggerEventFlag*>::const_iterator it = triggerFolderEventFlag_.begin();
189 it != triggerFolderEventFlag_.end();
200 if (!folderNames_.empty()) {
201 folderNames_.clear();
204 folderNames_.push_back(
"Uncleaned");
207 folderNames_.push_back(
"Uncleaned");
209 folderNames_.push_back(
"Cleaned");
210 folderNames_.push_back(
"DiJet");
211 if (!isMiniAODMet_) {
212 folderNames_.push_back(
"ZJets");
215 for (std::vector<std::string>::const_iterator ic = folderNames_.begin(); ic != folderNames_.end(); ic++) {
216 bookMESet(DirName +
"/" + *ic, ibooker, map_dijet_MEs);
223 std::map<std::string, MonitorElement*>& map_of_MEs) {
224 bool bLumiSecPlot = fill_met_high_level_histo;
226 bool fillPFCandidatePlots =
false;
227 bool fillZPlots =
false;
229 if (DirName.find(
"Cleaned") != std::string::npos) {
230 fillPFCandidatePlots =
true;
231 bookMonitorElement(DirName, ibooker, map_of_MEs, bLumiSecPlot, fillPFCandidatePlots, fillZPlots);
238 }
else if (DirName.find(
"ZJets") != std::string::npos) {
239 fillPFCandidatePlots =
false;
241 bookMonitorElement(DirName, ibooker, map_of_MEs, bLumiSecPlot, fillPFCandidatePlots, fillZPlots);
243 bookMonitorElement(DirName, ibooker, map_of_MEs, bLumiSecPlot, fillPFCandidatePlots, fillZPlots);
250 std::map<std::string, MonitorElement*>& map_of_MEs,
251 bool bLumiSecPlot =
false,
252 bool fillPFCandPlots =
false,
253 bool fillZPlots =
false) {
255 std::cout <<
"bookMonitorElement " << DirName << std::endl;
260 meZJets_u_par = ibooker.
book1D(
"u_parallel_Z_inc",
"u_parallel_Z_inc", 50, -1000., 75);
262 meZJets_u_par = ibooker.
book1D(
"u_parallel_Z_inc",
"u_parallel_Z_inc", 50, -800., 75);
264 meZJets_u_par_ZPt_0_15 = ibooker.
book1D(
"u_parallel_ZPt_0_15",
"u_parallel_ZPt_0_15", 50, -100, 75);
265 meZJets_u_par_ZPt_15_30 = ibooker.
book1D(
"u_parallel_ZPt_15_30",
"u_parallel_ZPt_15_30", 50, -100, 50);
266 meZJets_u_par_ZPt_30_55 = ibooker.
book1D(
"u_parallel_ZPt_30_55",
"u_parallel_ZPt_30_55", 50, -175, 50);
267 meZJets_u_par_ZPt_55_75 = ibooker.
book1D(
"u_parallel_ZPt_55_75",
"u_parallel_ZPt_55_75", 50, -175, 0);
268 meZJets_u_par_ZPt_75_150 = ibooker.
book1D(
"u_parallel_ZPt_75_150",
"u_parallel_ZPt_75_150", 50, -300, 0);
270 meZJets_u_par_ZPt_150_290 = ibooker.
book1D(
"u_parallel_ZPt_150_290",
"u_parallel_ZPt_150_290", 50, -750, -100);
272 meZJets_u_par_ZPt_150_290 = ibooker.
book1D(
"u_parallel_ZPt_150_290",
"u_parallel_ZPt_150_290", 50, -450, -50);
275 meZJets_u_par_ZPt_290 = ibooker.
book1D(
"u_parallel_ZPt_290",
"u_parallel_ZPt_290", 50, -1000., -350.);
277 meZJets_u_par_ZPt_290 = ibooker.
book1D(
"u_parallel_ZPt_290",
"u_parallel_ZPt_290", 50, -750., -150.);
279 meZJets_u_perp = ibooker.
book1D(
"u_perp_Z_inc",
"u_perp_Z_inc", 50, -85., 85.);
280 meZJets_u_perp_ZPt_0_15 = ibooker.
book1D(
"u_perp_ZPt_0_15",
"u_perp_ZPt_0_15", 50, -85., 85.);
281 meZJets_u_perp_ZPt_15_30 = ibooker.
book1D(
"u_perp_ZPt_15_30",
"u_perp_ZPt_15_30", 50, -85., 85.);
282 meZJets_u_perp_ZPt_30_55 = ibooker.
book1D(
"u_perp_ZPt_30_55",
"u_perp_ZPt_30_55", 50, -85., 85.);
283 meZJets_u_perp_ZPt_55_75 = ibooker.
book1D(
"u_perp_ZPt_55_75",
"u_perp_ZPt_55_75", 50, -85., 85.);
284 meZJets_u_perp_ZPt_75_150 = ibooker.
book1D(
"u_perp_ZPt_75_150",
"u_perp_ZPt_75_150", 50, -85., 85.);
285 meZJets_u_perp_ZPt_150_290 = ibooker.
book1D(
"u_perp_ZPt_150_290",
"u_perp_ZPt_150_290", 50, -85., 85.);
286 meZJets_u_perp_ZPt_290 = ibooker.
book1D(
"u_perp_ZPt_290",
"u_perp_ZPt_290", 50, -85., 85.);
288 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_parallel_Z_inc", meZJets_u_par));
290 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_parallel_ZPt_0_15", meZJets_u_par_ZPt_0_15));
292 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_parallel_ZPt_15_30", meZJets_u_par_ZPt_15_30));
294 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_parallel_ZPt_30_55", meZJets_u_par_ZPt_30_55));
296 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_parallel_ZPt_55_75", meZJets_u_par_ZPt_55_75));
298 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_parallel_ZPt_75_150", meZJets_u_par_ZPt_75_150));
300 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_parallel_ZPt_150_290", meZJets_u_par_ZPt_150_290));
302 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_parallel_ZPt_290", meZJets_u_par_ZPt_290));
304 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_perp_Z_inc", meZJets_u_perp));
306 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_perp_ZPt_0_15", meZJets_u_perp_ZPt_0_15));
308 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_perp_ZPt_15_30", meZJets_u_perp_ZPt_15_30));
310 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_perp_ZPt_30_55", meZJets_u_perp_ZPt_30_55));
312 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_perp_ZPt_55_75", meZJets_u_perp_ZPt_55_75));
314 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_perp_ZPt_75_150", meZJets_u_perp_ZPt_75_150));
316 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_perp_ZPt_150_290", meZJets_u_perp_ZPt_150_290));
318 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"u_perp_ZPt_290", meZJets_u_perp_ZPt_290));
322 hTrigger = ibooker.
book1D(
"triggerResults",
"triggerResults", 500, 0, 500);
323 for (
unsigned int i = 0;
i < allTriggerNames_.size();
i++) {
324 if (
i < (
unsigned int)hTrigger->getNbinsX()) {
325 if (!hTriggerLabelsIsSet_) {
330 hTriggerLabelsIsSet_ =
true;
332 hMEx = ibooker.
book1D(
"MEx",
"MEx", 200, -500, 500);
333 hMEy = ibooker.
book1D(
"MEy",
"MEy", 200, -500, 500);
334 hMET = ibooker.
book1D(
"MET",
"MET", 200, 0, 1000);
335 hMET_2 = ibooker.
book1D(
"MET_2",
"MET Range 2", 200, 0, 2000);
336 hSumET = ibooker.
book1D(
"SumET",
"SumET", 400, 0, 4000);
337 hMETSig = ibooker.
book1D(
"METSig",
"METSig", 51, 0, 51);
340 hMET_logx = ibooker.
book1D(
"MET_logx",
"MET_logx", 40, -1, 9);
341 hSumET_logx = ibooker.
book1D(
"SumET_logx",
"SumET_logx", 40, -1, 9);
344 hMEy->setAxisTitle(
"MEy [GeV]", 1);
345 hMET->setAxisTitle(
"MET [GeV]", 1);
346 hMET_2->setAxisTitle(
"MET [GeV]", 1);
347 hSumET->setAxisTitle(
"SumET [GeV]", 1);
348 hMETSig->setAxisTitle(
"METSig", 1);
349 hMETPhi->setAxisTitle(
"METPhi [rad]", 1);
350 hMET_logx->setAxisTitle(
"log(MET) [GeV]", 1);
351 hSumET_logx->setAxisTitle(
"log(SumET) [GeV]", 1);
353 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"triggerResults", hTrigger));
354 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MEx", hMEx));
355 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MEy", hMEy));
356 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET", hMET));
357 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_2", hMET_2));
358 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"SumET", hSumET));
359 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METSig", hMETSig));
360 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhi", hMETPhi));
361 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_logx", hMET_logx));
362 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"SumET_logx", hSumET_logx));
364 hMET_HBHENoiseFilter = ibooker.
book1D(
"MET_HBHENoiseFilter",
"MET_HBHENoiseFiltered", 200, 0, 1000);
365 hMET_CSCTightHaloFilter = ibooker.
book1D(
"MET_CSCTightHaloFilter",
"MET_CSCTightHaloFiltered", 200, 0, 1000);
366 hMET_eeBadScFilter = ibooker.
book1D(
"MET_eeBadScFilter",
"MET_eeBadScFiltered", 200, 0, 1000);
367 hMET_HBHEIsoNoiseFilter = ibooker.
book1D(
"MET_HBHEIsoNoiseFilter",
"MET_HBHEIsoNoiseFiltered", 200, 0, 1000);
368 hMET_CSCTightHalo2015Filter =
369 ibooker.
book1D(
"MET_CSCTightHalo2015Filter",
"MET_CSCTightHalo2015Filtered", 200, 0, 1000);
370 hMET_EcalDeadCellTriggerFilter =
371 ibooker.
book1D(
"MET_EcalDeadCellTriggerFilter",
"MET_EcalDeadCellTriggerFiltered", 200, 0, 1000);
372 hMET_EcalDeadCellBoundaryFilter =
373 ibooker.
book1D(
"MET_EcalDeadCellBoundaryFilter",
"MET_EcalDeadCellBoundaryFiltered", 200, 0, 1000);
374 hMET_HcalStripHaloFilter = ibooker.
book1D(
"MET_HcalStripHaloFilter",
"MET_HcalStripHaloFiltered", 200, 0, 1000);
377 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_HBHENoiseFilter", hMET_HBHENoiseFilter));
379 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_CSCTightHaloFilter", hMET_CSCTightHaloFilter));
380 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_eeBadScFilter", hMET_eeBadScFilter));
382 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_HBHEIsoNoiseFilter", hMET_HBHEIsoNoiseFilter));
383 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_CSCTightHalo2015Filter",
384 hMET_CSCTightHalo2015Filter));
385 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_EcalDeadCellTriggerFilter",
386 hMET_EcalDeadCellTriggerFilter));
387 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_EcalDeadCellBoundaryFilter",
388 hMET_EcalDeadCellBoundaryFilter));
390 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_HcalStripHaloFilter", hMET_HcalStripHaloFilter));
394 meMEx_profile = ibooker.
bookProfile(
"MEx_profile",
"met.px()", nbinsPV_, nPVMin_, nPVMax_, 200, -500, 500);
395 meMEy_profile = ibooker.
bookProfile(
"MEy_profile",
"met.py()", nbinsPV_, nPVMin_, nPVMax_, 200, -500, 500);
396 meMET_profile = ibooker.
bookProfile(
"MET_profile",
"met.pt()", nbinsPV_, nPVMin_, nPVMax_, 200, 0, 1000);
397 meSumET_profile = ibooker.
bookProfile(
"SumET_profile",
"met.sumEt()", nbinsPV_, nPVMin_, nPVMax_, 400, 0, 4000);
401 meMEy_profile->setAxisTitle(
"nvtx", 1);
402 meMET_profile->setAxisTitle(
"nvtx", 1);
403 meSumET_profile->setAxisTitle(
"nvtx", 1);
405 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MEx_profile", meMEx_profile));
406 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MEy_profile", meMEy_profile));
407 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_profile", meMET_profile));
408 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"SumET_profile", meSumET_profile));
411 hCaloHadEtInHB = ibooker.
book1D(
"CaloHadEtInHB",
"CaloHadEtInHB", 50, 0, 2000);
413 hCaloHadEtInHO = ibooker.
book1D(
"CaloHadEtInHO",
"CaloHadEtInHO", 25, 0, 500);
415 hCaloHadEtInHE = ibooker.
book1D(
"CaloHadEtInHE",
"CaloHadEtInHE", 50, 0, 2000);
417 hCaloHadEtInHF = ibooker.
book1D(
"CaloHadEtInHF",
"CaloHadEtInHF", 50, 0, 1000);
419 hCaloEmEtInHF = ibooker.
book1D(
"CaloEmEtInHF",
"CaloEmEtInHF", 25, 0, 500);
421 hCaloEmEtInEE = ibooker.
book1D(
"CaloEmEtInEE",
"CaloEmEtInEE", 50, 0, 1000);
423 hCaloEmEtInEB = ibooker.
book1D(
"CaloEmEtInEB",
"CaloEmEtInEB", 50, 0, 2000);
426 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"CaloHadEtInHO", hCaloHadEtInHO));
427 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"CaloHadEtInHF", hCaloHadEtInHF));
428 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"CaloHadEtInHE", hCaloHadEtInHE));
429 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"CaloHadEtInHB", hCaloHadEtInHB));
430 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"CaloEmEtInHF", hCaloEmEtInHF));
431 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"CaloEmEtInEE", hCaloEmEtInEE));
432 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"CaloEmEtInEB", hCaloEmEtInEB));
434 hCaloMETPhi020 = ibooker.
book1D(
"CaloMETPhi020",
"CaloMETPhi020", 60, -
M_PI,
M_PI);
435 hCaloMETPhi020->
setAxisTitle(
"METPhi [rad] (MET>20 GeV)", 1);
437 hCaloEtFractionHadronic = ibooker.
book1D(
"CaloEtFractionHadronic",
"CaloEtFractionHadronic", 50, 0, 1);
438 hCaloEtFractionHadronic->
setAxisTitle(
"Hadronic Et Fraction", 1);
439 hCaloEmEtFraction = ibooker.
book1D(
"CaloEmEtFraction",
"CaloEmEtFraction", 50, 0, 1);
442 hCaloEmEtFraction020 = ibooker.
book1D(
"CaloEmEtFraction020",
"CaloEmEtFraction020", 50, 0, 1);
443 hCaloEmEtFraction020->
setAxisTitle(
"EM Et Fraction (MET>20 GeV)", 1);
445 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"CaloMETPhi020", hCaloMETPhi020));
447 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"CaloEtFractionHadronic", hCaloEtFractionHadronic));
448 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"CaloEmEtFraction", hCaloEmEtFraction));
450 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"CaloEmEtFraction020", hCaloEmEtFraction020));
454 if (fillPFCandPlots &&
455 fillCandidateMap_histos) {
457 meCHF_Barrel = ibooker.
book1D(
"PfChargedHadronEtFractionBarrel",
"chargedHadronEtFractionBarrel", 50, 0, 1);
459 ibooker.
book1D(
"PfChargedHadronEtFractionEndcapPlus",
"chargedHadronEtFractionEndcapPlus", 50, 0, 1);
461 ibooker.
book1D(
"PfChargedHadronEtFractionEndcapMinus",
"chargedHadronEtFractionEndcapMinus", 50, 0, 1);
462 meCHF_Barrel_BXm1Empty = ibooker.
book1D(
463 "PfChargedHadronEtFractionBarrel_BXm1Empty",
"chargedHadronEtFractionBarrel prev empty bunch", 50, 0, 1);
464 meCHF_EndcapPlus_BXm1Empty = ibooker.
book1D(
"PfChargedHadronEtFractionEndcapPlus_BXm1Empty",
465 "chargedHadronEtFractionEndcapPlus prev empty bunch",
469 meCHF_EndcapMinus_BXm1Empty = ibooker.
book1D(
"PfChargedHadronEtFractionEndcapMinus_BXm1Empty",
470 "chargedHadronEtFractionEndcapMinus prev empty bunch",
474 meCHF_Barrel_BXm1Filled = ibooker.
book1D(
"PfChargedHadronEtFractionBarrel_BXm1Filled",
475 "chargedHadronEtFractionBarrel prev filled 2 bunches",
479 meCHF_EndcapPlus_BXm1Filled = ibooker.
book1D(
"PfChargedHadronEtFractionEndcapPlus_BXm1Filled",
480 "chargedHadronEtFractionEndcapPlus prev filled bunch",
484 meCHF_EndcapMinus_BXm1Filled = ibooker.
book1D(
"PfChargedHadronEtFractionEndcapMinus_BXm1Filled",
485 "chargedHadronEtFractionEndcapMinus prev filled bunch",
491 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfChargedHadronEtFractionBarrel", meCHF_Barrel));
492 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfChargedHadronEtFractionEndcapPlus",
494 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
495 DirName +
"/" +
"PfChargedHadronEtFractionEndcapMinus", meCHF_EndcapMinus));
496 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
497 DirName +
"/" +
"PfChargedHadronEtFractionBarrel_BXm1Empty", meCHF_Barrel_BXm1Empty));
498 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
499 DirName +
"/" +
"PfChargedHadronEtFractionEndcapPlus_BXm1Empty", meCHF_EndcapPlus_BXm1Empty));
500 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
501 DirName +
"/" +
"PfChargedHadronEtFractionEndcapMinus_BXm1Empty", meCHF_EndcapMinus_BXm1Empty));
505 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
506 DirName +
"/" +
"PfChargedHadronEtFractionBarrel_BXm1Filled", meCHF_Barrel_BXm1Filled));
507 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
508 DirName +
"/" +
"PfChargedHadronEtFractionEndcapPlus_BXm1Filled", meCHF_EndcapPlus_BXm1Filled));
509 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
510 DirName +
"/" +
"PfChargedHadronEtFractionEndcapMinus_BXm1Filled", meCHF_EndcapMinus_BXm1Filled));
515 meNHF_Barrel = ibooker.
book1D(
"PfNeutralHadronEtFractionBarrel",
"neutralHadronEtFractionBarrel", 50, 0, 1);
517 ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapPlus",
"neutralHadronEtFractionEndcapPlus", 50, 0, 1);
519 ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapMinus",
"neutralHadronEtFractionEndcapMinus", 50, 0, 1);
520 meNHF_Barrel_BXm1Empty = ibooker.
book1D(
521 "PfNeutralHadronEtFractionBarrel_BXm1Empty",
"neutralHadronEtFractionBarrel prev empty bunch", 50, 0, 1);
522 meNHF_EndcapPlus_BXm1Empty = ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapPlus_BXm1Empty",
523 "neutralHadronEtFractionEndcapPlus prev empty bunch",
527 meNHF_EndcapMinus_BXm1Empty = ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapMinus_BXm1Empty",
528 "neutralHadronEtFractionEndcapMinus prev empty bunch",
535 meNHF_Barrel_BXm1Filled = ibooker.
book1D(
"PfNeutralHadronEtFractionBarrel_BXm1Filled",
536 "neutralHadronEtFractionBarrel prev filled 2 bunches",
540 meNHF_EndcapPlus_BXm1Filled = ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapPlus_BXm1Filled",
541 "neutralHadronEtFractionEndcapPlus prev filled bunch",
545 meNHF_EndcapMinus_BXm1Filled = ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapMinus_BXm1Filled",
546 "neutralHadronEtFractionEndcapMinus prev filled bunch",
555 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfNeutralHadronEtFractionBarrel", meNHF_Barrel));
556 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfNeutralHadronEtFractionEndcapPlus",
558 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
559 DirName +
"/" +
"PfNeutralHadronEtFractionEndcapMinus", meNHF_EndcapMinus));
560 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
561 DirName +
"/" +
"PfNeutralHadronEtFractionBarrel_BXm1Empty", meNHF_Barrel_BXm1Empty));
562 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
563 DirName +
"/" +
"PfNeutralHadronEtFractionEndcapPlus_BXm1Empty", meNHF_EndcapPlus_BXm1Empty));
564 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
565 DirName +
"/" +
"PfNeutralHadronEtFractionEndcapMinus_BXm1Empty", meNHF_EndcapMinus_BXm1Empty));
569 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
570 DirName +
"/" +
"PfNeutralHadronEtFractionBarrel_BXm1Filled", meNHF_Barrel_BXm1Filled));
571 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
572 DirName +
"/" +
"PfNeutralHadronEtFractionEndcapPlus_BXm1Filled", meNHF_EndcapPlus_BXm1Filled));
573 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
574 DirName +
"/" +
"PfNeutralHadronEtFractionEndcapMinus_BXm1Filled", meNHF_EndcapMinus_BXm1Filled));
579 mePhF_Barrel = ibooker.
book1D(
"PfPhotonEtFractionBarrel",
"photonEtFractionBarrel", 50, 0, 1);
580 mePhF_EndcapPlus = ibooker.
book1D(
"PfPhotonEtFractionEndcapPlus",
"photonEtFractionEndcapPlus", 50, 0, 1);
581 mePhF_EndcapMinus = ibooker.
book1D(
"PfPhotonEtFractionEndcapMinus",
"photonEtFractionEndcapMinus", 50, 0, 1);
582 mePhF_Barrel_BXm1Empty =
583 ibooker.
book1D(
"PfPhotonEtFractionBarrel_BXm1Empty",
"photonEtFractionBarrel prev empty bunch", 50, 0, 1);
584 mePhF_EndcapPlus_BXm1Empty = ibooker.
book1D(
585 "PfPhotonEtFractionEndcapPlus_BXm1Empty",
"photonEtFractionEndcapPlus prev empty bunch", 50, 0, 1);
586 mePhF_EndcapMinus_BXm1Empty = ibooker.
book1D(
587 "PfPhotonEtFractionEndcapMinus_BXm1Empty",
"photonEtFractionEndcapMinus prev empty bunch", 50, 0, 1);
591 mePhF_Barrel_BXm1Filled = ibooker.
book1D(
592 "PfPhotonEtFractionBarrel_BXm1Filled",
"photonEtFractionBarrel prev filled 2 bunches", 50, 0, 1);
593 mePhF_EndcapPlus_BXm1Filled = ibooker.
book1D(
594 "PfPhotonEtFractionEndcapPlus_BXm1Filled",
"photonEtFractionEndcapPlus prev filled bunch", 50, 0, 1);
595 mePhF_EndcapMinus_BXm1Filled = ibooker.
book1D(
596 "PfPhotonEtFractionEndcapMinus_BXm1Filled",
"photonEtFractionEndcapMinus prev filled bunch", 50, 0, 1);
602 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfPhotonEtFractionBarrel", mePhF_Barrel));
604 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfPhotonEtFractionEndcapPlus", mePhF_EndcapPlus));
605 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfPhotonEtFractionEndcapMinus",
607 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfPhotonEtFractionBarrel_BXm1Empty",
608 mePhF_Barrel_BXm1Empty));
609 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
610 DirName +
"/" +
"PfPhotonEtFractionEndcapPlus_BXm1Empty", mePhF_EndcapPlus_BXm1Empty));
611 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
612 DirName +
"/" +
"PfPhotonEtFractionEndcapMinus_BXm1Empty", mePhF_EndcapMinus_BXm1Empty));
616 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfPhotonEtFractionBarrel_BXm1Filled",
617 mePhF_Barrel_BXm1Filled));
618 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
619 DirName +
"/" +
"PfPhotonEtFractionEndcapPlus_BXm1Filled", mePhF_EndcapPlus_BXm1Filled));
620 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
621 DirName +
"/" +
"PfPhotonEtFractionEndcapMinus_BXm1Filled", mePhF_EndcapMinus_BXm1Filled));
626 meHFHadF_Plus = ibooker.
book1D(
"PfHFHadronEtFractionPlus",
"HFHadronEtFractionPlus", 50, 0, 1);
627 meHFHadF_Minus = ibooker.
book1D(
"PfHFHadronEtFractionMinus",
"HFHadronEtFractionMinus", 50, 0, 1);
628 meHFHadF_Plus_BXm1Empty =
629 ibooker.
book1D(
"PfHFHadronEtFractionPlus_BXm1Empty",
"HFHadronEtFractionPlus prev empty bunch", 50, 0, 1);
630 meHFHadF_Minus_BXm1Empty =
631 ibooker.
book1D(
"PfHFHadronEtFractionMinus_BXm1Empty",
"HFHadronEtFractionMinus prev empty bunch", 50, 0, 1);
634 meHFHadF_Plus_BXm1Filled =
635 ibooker.
book1D(
"PfHFHadronEtFractionPlus_BXm1Filled",
"HFHadronEtFractionPlus prev filled bunch", 50, 0, 1);
636 meHFHadF_Minus_BXm1Filled = ibooker.
book1D(
637 "PfHFHadronEtFractionMinus_BXm1Filled",
"HFHadronEtFractionMinus prev filled bunch", 50, 0, 1);
642 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFHadronEtFractionPlus", meHFHadF_Plus));
644 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFHadronEtFractionMinus", meHFHadF_Minus));
645 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFHadronEtFractionPlus_BXm1Empty",
646 meHFHadF_Plus_BXm1Empty));
647 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFHadronEtFractionMinus_BXm1Empty",
648 meHFHadF_Minus_BXm1Empty));
651 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFHadronEtFractionPlus_BXm1Filled",
652 meHFHadF_Plus_BXm1Filled));
653 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
654 DirName +
"/" +
"PfHFHadronEtFractionMinus_BXm1Filled", meHFHadF_Minus_BXm1Filled));
658 meHFEMF_Plus = ibooker.
book1D(
"PfHFEMEtFractionPlus",
"HFEMEtFractionPlus", 50, 0, 1);
659 meHFEMF_Minus = ibooker.
book1D(
"PfHFEMEtFractionMinus",
"HFEMEtFractionMinus", 50, 0, 1);
660 meHFEMF_Plus_BXm1Empty =
661 ibooker.
book1D(
"PfHFEMEtFractionPlus_BXm1Empty",
"HFEMEtFractionPlus prev empty bunch", 50, 0, 1);
662 meHFEMF_Minus_BXm1Empty =
663 ibooker.
book1D(
"PfHFEMEtFractionMinus_BXm1Empty",
"HFEMEtFractionMinus prev empty bunch", 50, 0, 1);
666 meHFEMF_Plus_BXm1Filled =
667 ibooker.
book1D(
"PfHFEMEtFractionPlus_BXm1Filled",
"HFEMEtFractionPlus prev filled bunch", 50, 0, 1);
668 meHFEMF_Minus_BXm1Filled =
669 ibooker.
book1D(
"PfHFEMEtFractionMinus_BXm1Filled",
"HFEMEtFractionMinus prev filled bunch", 50, 0, 1);
674 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFEMEtFractionPlus", meHFEMF_Plus));
676 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFEMEtFractionMinus", meHFEMF_Minus));
677 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFEMEtFractionPlus_BXm1Empty",
678 meHFEMF_Plus_BXm1Empty));
679 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFEMEtFractionMinus_BXm1Empty",
680 meHFEMF_Minus_BXm1Empty));
683 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFEMEtFractionPlus_BXm1Filled",
684 meHFEMF_Plus_BXm1Filled));
685 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFEMEtFractionMinus_BXm1Filled",
686 meHFEMF_Minus_BXm1Filled));
747 mePhotonEtFraction_BXm1Empty =
748 ibooker.
book1D(
"PfPhotonEtFraction_BXm1Empty",
"photonEtFraction() prev empty bunch", 50, 0, 1);
749 mePhotonEtFraction_BXm1Filled =
750 ibooker.
book1D(
"PfPhotonEtFraction_BXm1Filled",
"photonEtFraction() prev filled bunch", 50, 0, 1);
751 meNeutralHadronEtFraction_BXm1Empty = ibooker.
book1D(
752 "PfNeutralHadronEtFraction_BXm1Empty",
"neutralHadronEtFraction() prev empty bunch", 50, 0, 1);
753 meNeutralHadronEtFraction_BXm1Filled = ibooker.
book1D(
754 "PfNeutralHadronEtFraction_BXm1Filled",
"neutralHadronEtFraction() prev filled bunch", 50, 0, 1);
755 meChargedHadronEtFraction_BXm1Empty = ibooker.
book1D(
756 "PfChargedHadronEtFraction_BXm1Empty",
"chargedHadronEtFraction() prev empty bunch", 50, 0, 1);
757 meChargedHadronEtFraction_BXm1Filled = ibooker.
book1D(
758 "PfChargedHadronEtFraction_BXm1Filled",
"chargedHadronEtFraction() prev filled bunch", 50, 0, 1);
759 meMET_BXm1Empty = ibooker.
book1D(
"MET_BXm1Empty",
"MET prev empty bunch", 200, 0, 1000);
760 meMET_BXm1Filled = ibooker.
book1D(
"MET_BXm1Filled",
"MET prev filled bunch", 200, 0, 1000);
761 meSumET_BXm1Empty = ibooker.
book1D(
"SumET_BXm1Empty",
"SumET prev empty bunch", 400, 0, 4000);
762 meSumET_BXm1Filled = ibooker.
book1D(
"SumET_BXm1Filled",
"SumET prev filled bunch", 400, 0, 4000);
764 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfPhotonEtFraction_BXm1Empty",
765 mePhotonEtFraction_BXm1Empty));
766 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfPhotonEtFraction_BXm1Filled",
767 mePhotonEtFraction_BXm1Filled));
768 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfNeutralHadronEtFraction_BXm1Empty",
769 meNeutralHadronEtFraction_BXm1Empty));
770 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfNeutralEtFraction_BXm1Filled",
771 meNeutralHadronEtFraction_BXm1Filled));
772 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfChargedHadronEtFraction_BXm1Empty",
773 meChargedHadronEtFraction_BXm1Empty));
774 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfChargedEtFraction_BXm1Filled",
775 meChargedHadronEtFraction_BXm1Filled));
776 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_BXm1Empty", meMET_BXm1Empty));
777 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MET_BXm1Filled", meMET_BXm1Filled));
779 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"SumET_BXm1Empty", meSumET_BXm1Empty));
781 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"SumET_BXm1Filled", meSumET_BXm1Filled));
783 meMETPhiChargedHadronsBarrel_BXm1Filled = ibooker.
book1D(
784 "METPhiChargedHadronsBarrel_BXm1Filled",
"METPhi_PFChargedHadronsBarrel prev bunch filled", 50, -
M_PI,
M_PI);
785 meMETPhiChargedHadronsEndcapPlus_BXm1Filled =
786 ibooker.
book1D(
"METPhiChargedHadronsEndcapPlus_BXm1Filled",
787 "METPhi_PFChargedHadronsEndcapPlus prev bunch filled",
791 meMETPhiChargedHadronsEndcapMinus_BXm1Filled =
792 ibooker.
book1D(
"METPhiChargedHadronsEndcapMinus_BXm1Filled",
793 "METPhi_PFChargedHadronsEndcapMinus prev bunch filled",
797 meMETPhiNeutralHadronsBarrel_BXm1Filled = ibooker.
book1D(
798 "METPhiNeutralHadronsBarrel_BXm1Filled",
"METPhi_PFNeutralHadronsBarrel prev bunch filled", 50, -
M_PI,
M_PI);
799 meMETPhiNeutralHadronsEndcapPlus_BXm1Filled =
800 ibooker.
book1D(
"METPhiNeutralHadronsEndcapPlus_BXm1Filled",
801 "METPhi_PFNeutralHadronsEndcapPlus prev bunch filled",
805 meMETPhiNeutralHadronsEndcapMinus_BXm1Filled =
806 ibooker.
book1D(
"METPhiNeutralHadronsEndcapMinus_BXm1Filled",
807 "METPhi_PFNeutralHadronsEndcapMinus prev bunch filled",
811 meMETPhiPhotonsBarrel_BXm1Filled = ibooker.
book1D(
812 "METPhiPhotonsBarrel_BXm1Filled",
"METPhi_PFPhotonsBarrel prev bunch filled", 50, -
M_PI,
M_PI);
813 meMETPhiPhotonsEndcapPlus_BXm1Filled = ibooker.
book1D(
814 "METPhiPhotonsEndcapPlus_BXm1Filled",
"METPhi_PFPhotonsEndcapPlus prev bunch filled", 50, -
M_PI,
M_PI);
815 meMETPhiPhotonsEndcapMinus_BXm1Filled = ibooker.
book1D(
816 "METPhiPhotonsEndcapMinus_BXm1Filled",
"METPhi_PFPhotonsEndcapMinus prev bunch filled", 50, -
M_PI,
M_PI);
817 meMETPhiHFHadronsPlus_BXm1Filled = ibooker.
book1D(
818 "METPhiHFHadronsPlus_BXm1Filled",
"METPhi_PFHFHadronsPlus prev bunch filled", 50, -
M_PI,
M_PI);
819 meMETPhiHFHadronsMinus_BXm1Filled = ibooker.
book1D(
820 "METPhiHFHadronsMinus_BXm1Filled",
"METPhi_PFHFHadronsMinus prev bunch filled", 50, -
M_PI,
M_PI);
821 meMETPhiHFEGammasPlus_BXm1Filled = ibooker.
book1D(
822 "METPhiHFEGammasPlus_BXm1Filled",
"METPhi_PFHFEGammasPlus prev bunch filled", 50, -
M_PI,
M_PI);
823 meMETPhiHFEGammasMinus_BXm1Filled = ibooker.
book1D(
824 "METPhiHFEGammasMinus_BXm1Filled",
"METPhi_PFHFEGammasMinus prev bunch filled", 50, -
M_PI,
M_PI);
826 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
827 DirName +
"/" +
"METPhiChargedHadronsBarrel_BXm1Filled", meMETPhiChargedHadronsBarrel_BXm1Filled));
828 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
829 DirName +
"/" +
"METPhiChargedHadronsEndcapPlus_BXm1Filled", meMETPhiChargedHadronsEndcapPlus_BXm1Filled));
831 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiChargedHadronsEndcapMinus_BXm1Filled",
832 meMETPhiChargedHadronsEndcapMinus_BXm1Filled));
833 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
834 DirName +
"/" +
"METPhiNeutralHadronsBarrel_BXm1Filled", meMETPhiNeutralHadronsBarrel_BXm1Filled));
835 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
836 DirName +
"/" +
"METPhiNeutralHadronsEndcapPlus_BXm1Filled", meMETPhiNeutralHadronsEndcapPlus_BXm1Filled));
838 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiNeutralHadronsEndcapMinus_BXm1Filled",
839 meMETPhiNeutralHadronsEndcapMinus_BXm1Filled));
840 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiPhotonsBarrel_BXm1Filled",
841 meMETPhiPhotonsBarrel_BXm1Filled));
842 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiPhotonsEndcapPlus_BXm1Filled",
843 meMETPhiPhotonsEndcapPlus_BXm1Filled));
844 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiPhotonsEndcapMinus_BXm1Filled",
845 meMETPhiPhotonsEndcapMinus_BXm1Filled));
846 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiHFHadronsPlus_BXm1Filled",
847 meMETPhiHFHadronsPlus_BXm1Filled));
848 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiHFHadronsMinus_BXm1Filled",
849 meMETPhiHFHadronsMinus_BXm1Filled));
850 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiHFEGammasPlus_BXm1Filled",
851 meMETPhiHFEGammasPlus_BXm1Filled));
852 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiHFEGammasMinus_BXm1Filled",
853 meMETPhiHFEGammasMinus_BXm1Filled));
855 meMETPhiChargedHadronsBarrel_BXm1Empty = ibooker.
book1D(
856 "METPhiChargedHadronsBarrel_BXm1Empty",
"METPhi_PFChargedHadronsBarrel prev bunch empty", 50, -
M_PI,
M_PI);
857 meMETPhiChargedHadronsEndcapPlus_BXm1Empty =
858 ibooker.
book1D(
"METPhiChargedHadronsEndcapPlus_BXm1Empty",
859 "METPhi_PFChargedHadronsEndcapPlus prev bunch empty",
863 meMETPhiChargedHadronsEndcapMinus_BXm1Empty =
864 ibooker.
book1D(
"METPhiChargedHadronsEndcapMinus_BXm1Empty",
865 "METPhi_PFChargedHadronsEndcapMinus prev bunch empty",
869 meMETPhiNeutralHadronsBarrel_BXm1Empty = ibooker.
book1D(
870 "METPhiNeutralHadronsBarrel_BXm1Empty",
"METPhi_PFNeutralHadronsBarrel prev bunch empty", 50, -
M_PI,
M_PI);
871 meMETPhiNeutralHadronsEndcapPlus_BXm1Empty =
872 ibooker.
book1D(
"METPhiNeutralHadronsEndcapPlus_BXm1Empty",
873 "METPhi_PFNeutralHadronsEndcapPlus prev bunch empty",
877 meMETPhiNeutralHadronsEndcapMinus_BXm1Empty =
878 ibooker.
book1D(
"METPhiNeutralHadronsEndcapMinus_BXm1Empty",
879 "METPhi_PFNeutralHadronsEndcapMinus prev bunch empty",
883 meMETPhiPhotonsBarrel_BXm1Empty =
884 ibooker.
book1D(
"METPhiPhotonsBarrel_BXm1Empty",
"METPhi_PFPhotonsBarrel prev bunch empty", 50, -
M_PI,
M_PI);
885 meMETPhiPhotonsEndcapPlus_BXm1Empty = ibooker.
book1D(
886 "METPhiPhotonsEndcapPlus_BXm1Empty",
"METPhi_PFPhotonsEndcapPlus prev bunch empty", 50, -
M_PI,
M_PI);
887 meMETPhiPhotonsEndcapMinus_BXm1Empty = ibooker.
book1D(
888 "METPhiPhotonsEndcapMinus_BXm1Empty",
"METPhi_PFPhotonsEndcapMinus prev bunch empty", 50, -
M_PI,
M_PI);
889 meMETPhiHFHadronsPlus_BXm1Empty =
890 ibooker.
book1D(
"METPhiHFHadronsPlus_BXm1Empty",
"METPhi_PFHFHadronsPlus prev bunch empty", 50, -
M_PI,
M_PI);
891 meMETPhiHFHadronsMinus_BXm1Empty = ibooker.
book1D(
892 "METPhiHFHadronsMinus_BXm1Empty",
"METPhi_PFHFHadronsMinus prev bunch empty", 50, -
M_PI,
M_PI);
893 meMETPhiHFEGammasPlus_BXm1Empty =
894 ibooker.
book1D(
"METPhiHFEGammasPlus_BXm1Empty",
"METPhi_PFHFEGammasPlus prev bunch empty", 50, -
M_PI,
M_PI);
895 meMETPhiHFEGammasMinus_BXm1Empty = ibooker.
book1D(
896 "METPhiHFEGammasMinus_BXm1Empty",
"METPhi_PFHFEGammasMinus prev bunch empty", 50, -
M_PI,
M_PI);
898 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
899 DirName +
"/" +
"METPhiChargedHadronsBarrel_BXm1Empty", meMETPhiChargedHadronsBarrel_BXm1Empty));
900 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
901 DirName +
"/" +
"METPhiChargedHadronsEndcapPlus_BXm1Empty", meMETPhiChargedHadronsEndcapPlus_BXm1Empty));
902 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
903 DirName +
"/" +
"METPhiChargedHadronsEndcapMinus_BXm1Empty", meMETPhiChargedHadronsEndcapMinus_BXm1Empty));
904 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
905 DirName +
"/" +
"METPhiNeutralHadronsBarrel_BXm1Empty", meMETPhiNeutralHadronsBarrel_BXm1Empty));
906 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
907 DirName +
"/" +
"METPhiNeutralHadronsEndcapPlus_BXm1Empty", meMETPhiNeutralHadronsEndcapPlus_BXm1Empty));
908 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
909 DirName +
"/" +
"METPhiNeutralHadronsEndcapMinus_BXm1Empty", meMETPhiNeutralHadronsEndcapMinus_BXm1Empty));
910 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiPhotonsBarrel_BXm1Empty",
911 meMETPhiPhotonsBarrel_BXm1Empty));
912 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiPhotonsEndcapPlus_BXm1Empty",
913 meMETPhiPhotonsEndcapPlus_BXm1Empty));
914 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiPhotonsEndcapMinus_BXm1Empty",
915 meMETPhiPhotonsEndcapMinus_BXm1Empty));
916 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiHFHadronsPlus_BXm1Empty",
917 meMETPhiHFHadronsPlus_BXm1Empty));
918 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiHFHadronsMinus_BXm1Empty",
919 meMETPhiHFHadronsMinus_BXm1Empty));
920 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiHFEGammasPlus_BXm1Empty",
921 meMETPhiHFEGammasPlus_BXm1Empty));
922 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiHFEGammasMinus_BXm1Empty",
923 meMETPhiHFEGammasMinus_BXm1Empty));
925 meMETPhiChargedHadronsBarrel =
926 ibooker.
book1D(
"METPhiChargedHadronsBarrel",
"METPhi_PFChargedHadronsBarrel", 50, -
M_PI,
M_PI);
927 meMETPhiChargedHadronsEndcapPlus =
928 ibooker.
book1D(
"METPhiChargedHadronsEndcapPlus",
"METPhi_PFChargedHadronsEndcapPlus", 50, -
M_PI,
M_PI);
929 meMETPhiChargedHadronsEndcapMinus =
930 ibooker.
book1D(
"METPhiChargedHadronsEndcapMinus",
"METPhi_PFChargedHadronsEndcapMinus", 50, -
M_PI,
M_PI);
931 meMETPhiNeutralHadronsBarrel =
932 ibooker.
book1D(
"METPhiNeutralHadronsBarrel",
"METPhi_PFNeutralHadronsBarrel", 50, -
M_PI,
M_PI);
933 meMETPhiNeutralHadronsEndcapPlus =
934 ibooker.
book1D(
"METPhiNeutralHadronsEndcapPlus",
"METPhi_PFNeutralHadronsEndcapPlus", 50, -
M_PI,
M_PI);
935 meMETPhiNeutralHadronsEndcapMinus =
936 ibooker.
book1D(
"METPhiNeutralHadronsEndcapMinus",
"METPhi_PFNeutralHadronsEndcapMinus", 50, -
M_PI,
M_PI);
937 meMETPhiPhotonsBarrel = ibooker.
book1D(
"METPhiPhotonsBarrel",
"METPhi_PFPhotonsBarrel", 50, -
M_PI,
M_PI);
938 meMETPhiPhotonsEndcapPlus =
939 ibooker.
book1D(
"METPhiPhotonsEndcapPlus",
"METPhi_PFPhotonsEndcapPlus", 50, -
M_PI,
M_PI);
940 meMETPhiPhotonsEndcapMinus =
941 ibooker.
book1D(
"METPhiPhotonsEndcapMinus",
"METPhi_PFPhotonsEndcapMinus", 50, -
M_PI,
M_PI);
942 meMETPhiHFHadronsPlus = ibooker.
book1D(
"METPhiHFHadronsPlus",
"METPhi_PFHFHadronsPlus", 50, -
M_PI,
M_PI);
943 meMETPhiHFHadronsMinus = ibooker.
book1D(
"METPhiHFHadronsMinus",
"METPhi_PFHFHadronsMinus", 50, -
M_PI,
M_PI);
944 meMETPhiHFEGammasPlus = ibooker.
book1D(
"METPhiHFEGammasPlus",
"METPhi_PFHFEGammasPlus", 50, -
M_PI,
M_PI);
945 meMETPhiHFEGammasMinus = ibooker.
book1D(
"METPhiHFEGammasMinus",
"METPhi_PFHFEGammasMinus", 50, -
M_PI,
M_PI);
947 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiChargedHadronsBarrel",
948 meMETPhiChargedHadronsBarrel));
949 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiChargedHadronsEndcapPlus",
950 meMETPhiChargedHadronsEndcapPlus));
951 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiChargedHadronsEndcapMinus",
952 meMETPhiChargedHadronsEndcapMinus));
953 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiNeutralHadronsBarrel",
954 meMETPhiNeutralHadronsBarrel));
955 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiNeutralHadronsEndcapPlus",
956 meMETPhiNeutralHadronsEndcapPlus));
957 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiNeutralHadronsEndcapMinus",
958 meMETPhiNeutralHadronsEndcapMinus));
960 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiPhotonsBarrel", meMETPhiPhotonsBarrel));
961 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiPhotonsEndcapPlus",
962 meMETPhiPhotonsEndcapPlus));
963 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiPhotonsEndcapMinus",
964 meMETPhiPhotonsEndcapMinus));
966 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiHFHadronsPlus", meMETPhiHFHadronsPlus));
968 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiHFHadronsMinus", meMETPhiHFHadronsMinus));
970 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiHFEGammasPlus", meMETPhiHFEGammasPlus));
972 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METPhiHFEGammasMinus", meMETPhiHFEGammasMinus));
975 if (fillPFCandPlots && fillCandidateMap_histos) {
976 if (!profilePFCand_x_.empty()) {
977 etaMinPFCand_.clear();
978 etaMaxPFCand_.clear();
980 countsPFCand_.clear();
983 profilePFCand_x_.clear();
984 profilePFCand_y_.clear();
985 profilePFCand_x_name_.clear();
986 profilePFCand_y_name_.clear();
988 for (std::vector<edm::ParameterSet>::const_iterator
v = diagnosticsParameters_.begin();
989 v != diagnosticsParameters_.end();
991 double etaMinPFCand =
v->getParameter<
double>(
"etaMin");
992 double etaMaxPFCand =
v->getParameter<
double>(
"etaMax");
993 int nMinPFCand =
v->getParameter<
int>(
"nMin");
994 int nMaxPFCand =
v->getParameter<
int>(
"nMax");
995 int nbinsPFCand =
v->getParameter<
double>(
"nbins");
998 etaMinPFCand_.push_back(etaMinPFCand);
999 etaMaxPFCand_.push_back(etaMaxPFCand);
1000 typePFCand_.push_back(
v->getParameter<
int>(
"type"));
1001 countsPFCand_.push_back(0);
1002 MExPFCand_.push_back(0.);
1003 MEyPFCand_.push_back(0.);
1005 profilePFCand_x_.push_back(
1014 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
1015 DirName +
"/" + profilePFCand_x_name_[profilePFCand_x_name_.size() - 1],
1016 profilePFCand_x_[profilePFCand_x_.size() - 1]));
1017 profilePFCand_y_.push_back(
1026 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
1027 DirName +
"/" + profilePFCand_y_name_[profilePFCand_y_name_.size() - 1],
1028 profilePFCand_y_[profilePFCand_y_.size() - 1]));
1032 if (isMiniAODMet_) {
1033 if (fillPFCandPlots &&
1034 fillCandidateMap_histos) {
1035 if (!profilePFCand_x_.empty()) {
1036 etaMinPFCand_.clear();
1037 etaMaxPFCand_.clear();
1038 typePFCand_.clear();
1039 countsPFCand_.clear();
1040 profilePFCand_x_.clear();
1041 profilePFCand_y_.clear();
1043 for (std::vector<edm::ParameterSet>::const_iterator
v = diagnosticsParameters_.begin();
1044 v != diagnosticsParameters_.end();
1046 double etaMinPFCand =
v->getParameter<
double>(
"etaMin");
1047 double etaMaxPFCand =
v->getParameter<
double>(
"etaMax");
1049 etaMinPFCand_.push_back(etaMinPFCand);
1050 etaMaxPFCand_.push_back(etaMaxPFCand);
1051 typePFCand_.push_back(
v->getParameter<
int>(
"type"));
1052 countsPFCand_.push_back(0);
1053 MExPFCand_.push_back(0.);
1054 MEyPFCand_.push_back(0.);
1059 if (isPFMet_ || isMiniAODMet_) {
1060 mePhotonEtFraction = ibooker.
book1D(
"PfPhotonEtFraction",
"photonEtFraction()", 50, 0, 1);
1061 meNeutralHadronEtFraction = ibooker.
book1D(
"PfNeutralHadronEtFraction",
"neutralHadronEtFraction()", 50, 0, 1);
1062 meChargedHadronEtFraction = ibooker.
book1D(
"PfChargedHadronEtFraction",
"chargedHadronEtFraction()", 50, 0, 1);
1063 meHFHadronEtFraction = ibooker.
book1D(
"PfHFHadronEtFraction",
"HFHadronEtFraction()", 50, 0, 1);
1064 meHFEMEtFraction = ibooker.
book1D(
"PfHFEMEtFraction",
"HFEMEtFraction()", 50, 0, 1);
1067 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfPhotonEtFraction", mePhotonEtFraction));
1068 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfNeutralHadronEtFraction",
1069 meNeutralHadronEtFraction));
1070 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfChargedHadronEtFraction",
1071 meChargedHadronEtFraction));
1073 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFHadronEtFraction", meHFHadronEtFraction));
1074 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFEMEtFraction", meHFEMEtFraction));
1076 mePhotonEtFraction_profile =
1077 ibooker.
bookProfile(
"PfPhotonEtFraction_profile",
"photonEtFraction()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1);
1078 meNeutralHadronEtFraction_profile = ibooker.
bookProfile(
1079 "PfNeutralHadronEtFraction_profile",
"neutralHadronEtFraction()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1);
1080 meChargedHadronEtFraction_profile = ibooker.
bookProfile(
1081 "PfChargedHadronEtFraction_profile",
"chargedHadronEtFraction()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1);
1082 meHFHadronEtFraction_profile = ibooker.
bookProfile(
1083 "PfHFHadronEtFraction_profile",
"HFHadronEtFraction()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1);
1084 meHFEMEtFraction_profile =
1085 ibooker.
bookProfile(
"PfHFEMEtFraction_profile",
"HFEMEtFraction()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1);
1087 meNeutralHadronEtFraction_profile->setAxisTitle(
"nvtx", 1);
1088 meChargedHadronEtFraction_profile->setAxisTitle(
"nvtx", 1);
1089 meHFHadronEtFraction_profile->setAxisTitle(
"nvtx", 1);
1090 meHFEMEtFraction_profile->setAxisTitle(
"nvtx", 1);
1092 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfPhotonEtFraction_profile",
1093 mePhotonEtFraction_profile));
1094 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfNeutralHadronEtFraction_profile",
1095 meNeutralHadronEtFraction_profile));
1096 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfChargedHadronEtFraction_profile",
1097 meChargedHadronEtFraction_profile));
1098 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFHadronEtFraction_profile",
1099 meHFHadronEtFraction_profile));
1100 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFEMEtFraction_profile",
1101 meHFEMEtFraction_profile));
1103 mePhotonEt = ibooker.
book1D(
"PfPhotonEt",
"photonEt()", 50, 0, 1000);
1104 meNeutralHadronEt = ibooker.
book1D(
"PfNeutralHadronEt",
"neutralHadronEt()", 50, 0, 1000);
1105 meElectronEt = ibooker.
book1D(
"PfElectronEt",
"electronEt()", 50, 0, 100);
1106 meChargedHadronEt = ibooker.
book1D(
"PfChargedHadronEt",
"chargedHadronEt()", 50, 0, 2000);
1107 meMuonEt = ibooker.
book1D(
"PfMuonEt",
"muonEt()", 50, 0, 100);
1108 meHFHadronEt = ibooker.
book1D(
"PfHFHadronEt",
"HFHadronEt()", 50, 0, 2000);
1109 meHFEMEt = ibooker.
book1D(
"PfHFEMEt",
"HFEMEt()", 50, 0, 1000);
1111 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfPhotonEt", mePhotonEt));
1113 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfNeutralHadronEt", meNeutralHadronEt));
1114 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfElectronEt", meElectronEt));
1116 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfChargedHadronEt", meChargedHadronEt));
1117 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfMuonEt", meMuonEt));
1118 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFHadronEt", meHFHadronEt));
1119 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFEMEt", meHFEMEt));
1121 mePhotonEt_profile =
1122 ibooker.
bookProfile(
"PfPhotonEt_profile",
"photonEt()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1000);
1124 "PfNeutralHadronEt_profile",
"neutralHadronEt()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1000);
1126 "PfChargedHadronEt_profile",
"chargedHadronEt()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1000);
1127 meHFHadronEt_profile =
1128 ibooker.
bookProfile(
"PfHFHadronEt_profile",
"HFHadronEt()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1000);
1129 meHFEMEt_profile = ibooker.
bookProfile(
"PfHFEMEt_profile",
"HFEMEt()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1000);
1132 meNeutralHadronEt_profile->setAxisTitle(
"nvtx", 1);
1133 meChargedHadronEt_profile->setAxisTitle(
"nvtx", 1);
1134 meHFHadronEt_profile->setAxisTitle(
"nvtx", 1);
1135 meHFEMEt_profile->setAxisTitle(
"nvtx", 1);
1138 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfPhotonEt_profile", mePhotonEt_profile));
1139 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfNeutralHadronEt_profile",
1140 meNeutralHadronEt_profile));
1141 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfChargedHadronEt_profile",
1142 meChargedHadronEt_profile));
1144 std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFHadronEt_profile", meHFHadronEt_profile));
1145 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"PfHFEMEt_profile", meHFEMEt_profile));
1149 if (fill_met_high_level_histo) {
1150 hMExLS = ibooker.
book2D(
"MExLS",
"MEx_LS", 200, -200, 200, 250, 0., 2500.);
1152 hMExLS->setAxisTitle(
"Lumi Section", 2);
1153 hMExLS->setOption(
"colz");
1154 hMEyLS = ibooker.
book2D(
"MEyLS",
"MEy_LS", 200, -200, 200, 250, 0., 2500.);
1156 hMEyLS->setAxisTitle(
"Lumi Section", 2);
1157 hMEyLS->setOption(
"colz");
1158 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MExLS", hMExLS));
1159 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"MEyLS", hMEyLS));
1163 hMETRate = ibooker.
book1D(
"METRate",
"METRate", 200, 0, 1000);
1164 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName +
"/" +
"METRate", hMETRate));
1167 lumisecME = ibooker.
book1D(
"lumisec",
"lumisec", 2501, -1., 2500.);
1168 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
"JetMET/lumisec", lumisecME));
1179 if ((techTrig->second).algoName() == m_l1algoname_) {
1180 m_bitAlgTechTrig_ = (techTrig->second).algoBitNumber();
1188 if (hltConfig_.init(iRun, iSetup, triggerResultsLabel_.process(), changed)) {
1200 std::cout <<
"HLTEventAnalyzerAOD::analyze:" 1201 <<
" config extraction failure with process name " << triggerResultsLabel_.process() << std::endl;
1204 allTriggerNames_.clear();
1205 for (
unsigned int i = 0;
i < hltConfig_.size();
i++) {
1206 allTriggerNames_.push_back(hltConfig_.triggerName(
i));
1211 for (std::vector<GenericTriggerEventFlag*>::const_iterator it = triggerFolderEventFlag_.begin();
1212 it != triggerFolderEventFlag_.end();
1214 int pos = it - triggerFolderEventFlag_.begin();
1216 (*it)->initRun(iRun, iSetup);
1217 if (triggerSelectedSubFolders_[pos].exists(
std::string(
"hltDBKey"))) {
1218 if ((*it)->expressionsFromDB((*it)->hltDBKey(), iSetup)[0] !=
"CONFIG_ERROR")
1219 triggerFolderExpr_[pos] = (*it)->expressionsFromDB((*it)->hltDBKey(), iSetup);
1223 if (isMiniAODMet_) {
1224 bool changed_filter =
true;
1225 std::vector<int> initializeFilter(8, -1);
1226 miniaodFilterIndex_ = initializeFilter;
1227 if (FilterhltConfig_.init(iRun, iSetup, METFilterMiniAODLabel_.process(), changed_filter)) {
1228 miniaodfilterdec = 0;
1229 for (
unsigned int i = 0;
i < FilterhltConfig_.size();
i++) {
1233 HBHENoiseStringMiniAOD;
1234 std::size_t
found = search2.find(search);
1235 if (found != std::string::npos) {
1236 miniaodFilterIndex_[0] =
i;
1238 search2 = CSCHaloResultTag_.label().substr(0, CSCHaloResultTag_.label().size() - 3);
1239 found = search2.find(search);
1240 if (found != std::string::npos) {
1241 miniaodFilterIndex_[1] =
i;
1243 search2 = eeBadScFilterTag_.label().substr(0, eeBadScFilterTag_.label().size() - 3);
1244 found = search2.find(search);
1245 if (found != std::string::npos) {
1246 miniaodFilterIndex_[2] =
i;
1248 search2 = HBHEIsoNoiseStringMiniAOD;
1249 found = search2.find(search);
1250 if (found != std::string::npos) {
1251 miniaodFilterIndex_[3] =
i;
1253 search2 = CSCHalo2015ResultTag_.label().substr(0, CSCHalo2015ResultTag_.label().size() - 3);
1254 found = search2.find(search);
1255 if (found != std::string::npos) {
1256 miniaodFilterIndex_[4] =
i;
1258 search2 = EcalDeadCellTriggerTag_.label().substr(0, EcalDeadCellTriggerTag_.label().size() - 3);
1259 found = search2.find(search);
1260 if (found != std::string::npos) {
1261 miniaodFilterIndex_[5] =
i;
1263 search2 = EcalDeadCellBoundaryTag_.label().substr(0, EcalDeadCellBoundaryTag_.label().size() - 3);
1264 found = search2.find(search);
1265 if (found != std::string::npos) {
1266 miniaodFilterIndex_[6] =
i;
1268 search2 = HcalStripHaloTag_.label().substr(0, HcalStripHaloTag_.label().size() - 3);
1269 found = search2.find(search);
1270 if (found != std::string::npos) {
1271 miniaodFilterIndex_[7] =
i;
1274 }
else if (FilterhltConfig_.init(iRun, iSetup, METFilterMiniAODLabel2_.process(), changed_filter)) {
1275 miniaodfilterdec = 1;
1276 for (
unsigned int i = 0;
i < FilterhltConfig_.size();
i++) {
1280 HBHENoiseStringMiniAOD;
1281 std::size_t
found = search2.find(search);
1282 if (found != std::string::npos) {
1283 miniaodFilterIndex_[0] =
i;
1285 search2 = CSCHaloResultTag_.label().substr(0, CSCHaloResultTag_.label().size() - 3);
1286 found = search2.find(search);
1287 if (found != std::string::npos) {
1288 miniaodFilterIndex_[1] =
i;
1290 search2 = eeBadScFilterTag_.label().substr(0, eeBadScFilterTag_.label().size() - 3);
1291 found = search2.find(search);
1292 if (found != std::string::npos) {
1293 miniaodFilterIndex_[2] =
i;
1295 search2 = HBHEIsoNoiseStringMiniAOD;
1296 found = search2.find(search);
1297 if (found != std::string::npos) {
1298 miniaodFilterIndex_[3] =
i;
1300 search2 = CSCHalo2015ResultTag_.label().substr(0, CSCHalo2015ResultTag_.label().size() - 3);
1301 found = search2.find(search);
1302 if (found != std::string::npos) {
1303 miniaodFilterIndex_[4] =
i;
1305 search2 = EcalDeadCellTriggerTag_.label().substr(0, EcalDeadCellTriggerTag_.label().size() - 3);
1306 found = search2.find(search);
1307 if (found != std::string::npos) {
1308 miniaodFilterIndex_[5] =
i;
1310 search2 = EcalDeadCellBoundaryTag_.label().substr(0, EcalDeadCellBoundaryTag_.label().size() - 3);
1311 found = search2.find(search);
1312 if (found != std::string::npos) {
1313 miniaodFilterIndex_[6] =
i;
1315 search2 = HcalStripHaloTag_.label().substr(0, HcalStripHaloTag_.label().size() - 3);
1316 found = search2.find(search);
1317 if (found != std::string::npos) {
1318 miniaodFilterIndex_[7] =
i;
1323 <<
"nothing found with both RECO and reRECO label" << std::endl;
1336 meLumiSec = map_dijet_MEs[
"JetMET/lumisec"];
1339 int totlssecsum = 0;
1340 double totltime = 0.;
1342 tlumisec = meLumiSec->
getTH1F();
1347 for (
int i = 0;
i < (tlumisec->GetNbinsX());
i++) {
1348 if (tlumisec->GetBinContent(
i) != 0) {
1350 totlssecsum += tlumisec->GetBinContent(
i);
1353 int num_per_ls = (double)totlssecsum / (
double)totlsec;
1354 totlsec = totlsec + tlumisec->GetBinContent(tlumisec->GetNbinsX() + 1) / (
double)num_per_ls;
1355 totltime = double(totlsec * 90);
1366 for (std::vector<std::string>::const_iterator ic = folderNames_.begin(); ic != folderNames_.end(); ic++) {
1368 DirName = dirName + *ic;
1369 makeRatePlot(DirName, totltime);
1370 for (std::vector<GenericTriggerEventFlag*>::const_iterator it = triggerFolderEventFlag_.begin();
1371 it != triggerFolderEventFlag_.end();
1373 int pos = it - triggerFolderEventFlag_.begin();
1375 makeRatePlot(DirName +
"/" + triggerFolderLabels_[pos], totltime);
1385 MonitorElement* mMETRate = map_dijet_MEs[DirName +
"/" +
"METRate"];
1390 if (meMET && mMETRate) {
1395 tMETRate = (TH1F*)tMET->Clone(
"METRateHist");
1396 for (
int i = tMETRate->GetNbinsX() - 1;
i >= 0;
i--) {
1397 mMETRate->
setBinContent(
i + 1, tMETRate->GetBinContent(
i + 2) + tMET->GetBinContent(
i + 1));
1399 for (
int i = 0;
i < tMETRate->GetNbinsX();
i++) {
1400 mMETRate->
setBinContent(
i + 1, tMETRate->GetBinContent(
i + 1) / double(totltime));
1409 int myLuminosityBlock;
1411 if (fill_met_high_level_histo) {
1412 lumisecME = map_dijet_MEs[
"JetMET/lumisec"];
1413 if (lumisecME && lumisecME->getRootObject())
1414 lumisecME->Fill(myLuminosityBlock);
1417 if (myLuminosityBlock < LSBegin_)
1419 if (myLuminosityBlock > LSEnd_ && LSEnd_ > 0)
1423 std::cout <<
"METAnalyzer analyze" << std::endl;
1439 std::vector<int> triggerFolderDecisions;
1440 triggerFolderDecisions_ = std::vector<int>(triggerFolderEventFlag_.size(), 0);
1443 iEvent.
getByToken(triggerResultsToken_, triggerResults);
1445 if (triggerResults.isValid()) {
1448 int ntrigs = (*triggerResults).size();
1450 std::cout <<
"ntrigs=" << ntrigs << std::endl;
1452 for (std::vector<GenericTriggerEventFlag*>::const_iterator it = triggerFolderEventFlag_.begin();
1453 it != triggerFolderEventFlag_.end();
1455 unsigned int pos = it - triggerFolderEventFlag_.begin();
1456 bool fd = (*it)->accept(iEvent, iSetup);
1457 triggerFolderDecisions_[pos] =
fd;
1459 allTriggerDecisions_.clear();
1460 for (
unsigned int i = 0;
i < allTriggerNames_.size(); ++
i) {
1461 allTriggerDecisions_.push_back((*triggerResults).accept(
i));
1480 iEvent.
getByToken(caloMetToken_, calometcoll);
1489 if (isMiniAODMet_) {
1495 const MET*
met =
nullptr;
1497 const PFMET* pfmet =
nullptr;
1498 const CaloMET* calomet =
nullptr;
1503 assert(!pfmetcoll->empty());
1504 met = &(pfmetcoll->front());
1505 pfmet = &(pfmetcoll->front());
1508 assert(!calometcoll->empty());
1509 met = &(calometcoll->front());
1510 calomet = &(calometcoll->front());
1512 if (isMiniAODMet_) {
1513 assert(!patmetcoll->empty());
1514 met = &(patmetcoll->front());
1515 patmet = &(patmetcoll->front());
1518 LogTrace(
"METAnalyzer") <<
"[METAnalyzer] Call to the MET analyzer";
1540 bool bJetID =
false;
1541 bool bDiJetID =
false;
1555 LogDebug(
"") <<
"METAnalyzer: Could not find calojet product" << std::endl;
1557 std::cout <<
"METAnalyzer: Could not find calojet product" << std::endl;
1559 collsize = caloJets->size();
1574 iEvent.
getByToken(jetID_ValueMapToken_, jetID_ValueMap_Handle);
1578 if (isMiniAODMet_) {
1581 LogDebug(
"") <<
"METAnalyzer: Could not find patjet product" << std::endl;
1583 std::cout <<
"METAnalyzer: Could not find patjet product" << std::endl;
1585 collsize = patJets->size();
1591 LogDebug(
"") <<
"METAnalyzer: Could not find pfjet product" << std::endl;
1593 std::cout <<
"METAnalyzer: Could not find pfjet product" << std::endl;
1595 collsize = pfJets->size();
1598 unsigned int ind1 = -1;
1600 bool pass_jetID1 =
false;
1601 unsigned int ind2 = -1;
1603 bool pass_jetID2 =
false;
1606 bool pass_correction_flag =
false;
1607 if (!isMiniAODMet_) {
1608 iEvent.
getByToken(jetCorrectorToken_, jetCorr);
1610 pass_correction_flag =
true;
1613 pass_correction_flag =
true;
1616 for (
int ijet = 0; ijet < collsize; ijet++) {
1617 double pt_jet = -10;
1619 bool iscleaned =
false;
1620 if (pass_correction_flag) {
1622 scale = jetCorr->
correction((*caloJets)[ijet]);
1628 scale = jetCorr->
correction((*pfJets)[ijet]);
1632 pt_jet = scale * (*caloJets)[ijet].pt();
1633 if (pt_jet > ptThreshold_) {
1637 iscleaned = jetIDFunctorLoose((*caloJets)[ijet], jetID);
1659 pt_jet = scale * (*pfJets)[ijet].pt();
1660 if (pt_jet > ptThreshold_) {
1661 iscleaned = pfjetIDFunctorLoose((*pfJets)[ijet]);
1664 if (isMiniAODMet_) {
1665 pt_jet = (*patJets)[ijet].pt();
1666 if (pt_jet > ptThreshold_) {
1667 pat::strbitset stringbitset = pfjetIDFunctorLoose.getBitTemplate();
1668 iscleaned = pfjetIDFunctorLoose((*patJets)[ijet], stringbitset);
1677 pass_jetID2 = pass_jetID1;
1680 pass_jetID1 = iscleaned;
1681 }
else if (pt_jet > pt2) {
1684 pass_jetID2 = iscleaned;
1687 if (pass_jetID1 && pass_jetID2) {
1690 dphi = fabs((*caloJets)[ind1].phi() - (*caloJets)[ind2].phi());
1696 dphi = fabs((*pfJets)[ind1].phi() - (*pfJets)[ind2].phi());
1698 if (isMiniAODMet_) {
1699 dphi = fabs((*patJets)[0].phi() - (*patJets)[1].phi());
1701 if (dphi > acos(-1.)) {
1702 dphi = 2 * acos(-1.) - dphi;
1712 iEvent.
getByToken(vertexToken_, vertexHandle);
1714 if (!vertexHandle.
isValid()) {
1715 LogDebug(
"") <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1717 std::cout <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1722 numPV_ = vertexCollection.size();
1724 bool bPrimaryVertex = (bypassAllPVChecks_ || (numPV_ > 0));
1726 bool bZJets =
false;
1733 double pt_muon0 = -1;
1734 double pt_muon1 = -1;
1738 if (Muons.
isValid() && Muons->size() > 1) {
1739 for (
unsigned int i = 0;
i < Muons->size();
i++) {
1740 bool pass_muon_id =
false;
1741 bool pass_muon_iso =
false;
1742 double dxy = fabs((*Muons)[
i].muonBestTrack()->
dxy());
1743 double dz = fabs((*Muons)[
i].muonBestTrack()->
dz());
1745 dxy = fabs((*Muons)[
i].muonBestTrack()->
dxy((*vertexHandle)[0].
position()));
1746 dz = fabs((*Muons)[
i].muonBestTrack()->
dz((*vertexHandle)[0].
position()));
1748 if ((*Muons)[
i].
pt() > 20 && fabs((*Muons)[
i].
eta()) < 2.3) {
1749 if ((*Muons)[
i].isGlobalMuon() && (*Muons)[
i].isPFMuon() &&
1750 (*Muons)[
i].
globalTrack()->hitPattern().numberOfValidMuonHits() > 0 &&
1751 (*Muons)[
i].numberOfMatchedStations() > 1 && dxy < 0.2 && (*Muons)[
i].numberOfMatchedStations() > 1 &&
1753 (*Muons)[
i].
innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5) {
1754 pass_muon_id =
true;
1758 ((*Muons)[
i].pfIsolationR04().sumChargedHadronPt +
1763 if (muonIsoPFdb < 0.12) {
1764 pass_muon_iso =
true;
1767 if (pass_muon_id && pass_muon_iso) {
1768 if ((*Muons)[
i].
pt() > pt_muon0) {
1769 mu_index1 = mu_index0;
1770 pt_muon1 = pt_muon0;
1772 pt_muon0 = (*Muons)[
i].pt();
1773 }
else if ((*Muons)[
i].
pt() > pt_muon1) {
1775 pt_muon1 = (*Muons)[
i].pt();
1780 if (mu_index0 >= 0 && mu_index1 >= 0) {
1781 if ((*Muons)[mu_index0].
charge() * (*Muons)[mu_index1].
charge() < 0) {
1782 zCand = (*Muons)[mu_index0].polarP4() + (*Muons)[mu_index1].polarP4();
1783 if (fabs(zCand.M() - 91.) < 20) {
1793 iEvent.
getByToken(gtToken_, gtReadoutRecord);
1795 bool techTriggerResultBxM2 =
false;
1796 bool techTriggerResultBxM1 =
false;
1797 bool techTriggerResultBx0 =
false;
1799 if (!gtReadoutRecord.
isValid()) {
1800 LogDebug(
"") <<
"METAnalyzer: Could not find GT readout record" << std::endl;
1802 std::cout <<
"METAnalyzer: Could not find GT readout record product" << std::endl;
1810 if (m_bitAlgTechTrig_ > -1 && !technicalTriggerWordBeforeMaskBx0.empty()) {
1811 techTriggerResultBx0 = technicalTriggerWordBeforeMaskBx0.at(m_bitAlgTechTrig_);
1812 if (techTriggerResultBx0 != 0) {
1813 techTriggerResultBxM2 = technicalTriggerWordBeforeMaskBxM2.at(m_bitAlgTechTrig_);
1814 techTriggerResultBxM1 = technicalTriggerWordBeforeMaskBxM1.at(m_bitAlgTechTrig_);
1819 std::vector<bool> trigger_flag(4,
false);
1820 if (techTriggerResultBx0 && techTriggerResultBxM2 &&
1821 techTriggerResultBxM1) {
1822 trigger_flag[0] =
true;
1824 if (techTriggerResultBx0 && techTriggerResultBxM1) {
1825 trigger_flag[1] =
true;
1827 if (techTriggerResultBx0 && !techTriggerResultBxM1) {
1828 trigger_flag[2] =
true;
1830 if (techTriggerResultBx0 && !techTriggerResultBxM2 &&
1831 !techTriggerResultBxM1) {
1832 trigger_flag[3] =
true;
1834 std::vector<bool> filter_decisions(
1836 if (!isMiniAODMet_ &&
1839 iEvent.
getByToken(hbheNoiseFilterResultToken_, HBHENoiseFilterResultHandle);
1840 if (!HBHENoiseFilterResultHandle.
isValid()) {
1841 LogDebug(
"") <<
"METAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
1843 std::cout <<
"METAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
1845 filter_decisions[0] = *HBHENoiseFilterResultHandle;
1847 iEvent.
getByToken(CSCHaloResultToken_, CSCTightHaloFilterResultHandle);
1848 if (!CSCTightHaloFilterResultHandle.
isValid()) {
1849 LogDebug(
"") <<
"METAnalyzer: Could not find CSCTightHaloFilterResultHandle" << std::endl;
1851 std::cout <<
"METAnalyzer: CSCTightHaloFilterResultHandle" << std::endl;
1853 filter_decisions[1] = *CSCTightHaloFilterResultHandle;
1855 iEvent.
getByToken(eeBadScFilterToken_, eeBadScFilterResultHandle);
1856 if (!eeBadScFilterResultHandle.
isValid()) {
1857 LogDebug(
"") <<
"METAnalyzer: Could not find eeBadScFilterResultHandle" << std::endl;
1859 std::cout <<
"METAnalyzer: eeBadScFilterResultHandle" << std::endl;
1861 filter_decisions[2] = *eeBadScFilterResultHandle;
1863 iEvent.
getByToken(hbheIsoNoiseFilterResultToken_, HBHENoiseIsoFilterResultHandle);
1864 if (!HBHENoiseIsoFilterResultHandle.
isValid()) {
1865 LogDebug(
"") <<
"METAnalyzer: Could not find HBHENoiseIsoFilterResult" << std::endl;
1867 std::cout <<
"METAnalyzer: Could not find HBHENoiseIsoFilterResult" << std::endl;
1869 filter_decisions[3] = *HBHENoiseIsoFilterResultHandle;
1871 iEvent.
getByToken(CSCHalo2015ResultToken_, CSCTightHalo2015FilterResultHandle);
1872 if (!CSCTightHalo2015FilterResultHandle.
isValid()) {
1873 LogDebug(
"") <<
"METAnalyzer: Could not find CSCTightHalo2015FilterResultHandle" << std::endl;
1875 std::cout <<
"METAnalyzer: CSCTightHalo2015FilterResultHandle" << std::endl;
1877 filter_decisions[4] = *CSCTightHalo2015FilterResultHandle;
1879 iEvent.
getByToken(EcalDeadCellTriggerToken_, EcalDeadCellTriggerFilterResultHandle);
1880 if (!EcalDeadCellTriggerFilterResultHandle.
isValid()) {
1881 LogDebug(
"") <<
"METAnalyzer: Could not find EcalDeadCellTriggerFilterResultHandle" << std::endl;
1883 std::cout <<
"METAnalyzer: EcalDeadCellTriggerFilterResultHandle" << std::endl;
1885 filter_decisions[5] = *EcalDeadCellTriggerFilterResultHandle;
1887 iEvent.
getByToken(EcalDeadCellBoundaryToken_, EcalDeadCellBoundaryHandle);
1888 if (!EcalDeadCellBoundaryHandle.
isValid()) {
1889 LogDebug(
"") <<
"METAnalyzer: Could not find EcalDeadCellBoundaryHandle" << std::endl;
1891 std::cout <<
"METAnalyzer: EcalDeadCellBoundaryHandle" << std::endl;
1893 filter_decisions[6] = *EcalDeadCellBoundaryHandle;
1895 iEvent.
getByToken(HcalStripHaloToken_, HcalStripHaloFilterHandle);
1896 if (!HcalStripHaloFilterHandle.
isValid()) {
1897 LogDebug(
"") <<
"METAnalyzer: Could not find CSCTightHalo2015FilterResultHandle" << std::endl;
1899 std::cout <<
"METAnalyzer: CSCTightHalo2015FilterResultHandle" << std::endl;
1901 filter_decisions[7] = *HcalStripHaloFilterHandle;
1902 }
else if (isMiniAODMet_) {
1905 iEvent.
getByToken(METFilterMiniAODToken_, metFilterResults);
1906 if (metFilterResults.
isValid()) {
1907 if (miniaodFilterIndex_[0] != -1) {
1908 filter_decisions[0] = metFilterResults->
accept(miniaodFilterIndex_[0]);
1910 if (miniaodFilterIndex_[1] != -1) {
1911 filter_decisions[1] = metFilterResults->
accept(miniaodFilterIndex_[1]);
1913 if (miniaodFilterIndex_[2] != -1) {
1914 filter_decisions[2] = metFilterResults->
accept(miniaodFilterIndex_[2]);
1916 if (miniaodFilterIndex_[3] != -1) {
1917 filter_decisions[3] = metFilterResults->
accept(miniaodFilterIndex_[3]);
1919 if (miniaodFilterIndex_[4] != -1) {
1920 filter_decisions[4] = metFilterResults->
accept(miniaodFilterIndex_[4]);
1922 if (miniaodFilterIndex_[5] != -1) {
1923 filter_decisions[5] = metFilterResults->
accept(miniaodFilterIndex_[5]);
1925 if (miniaodFilterIndex_[6] != -1) {
1926 filter_decisions[6] = metFilterResults->
accept(miniaodFilterIndex_[6]);
1928 if (miniaodFilterIndex_[7] != -1) {
1929 filter_decisions[7] = metFilterResults->
accept(miniaodFilterIndex_[7]);
1932 iEvent.
getByToken(METFilterMiniAODToken2_, metFilterResults);
1933 if (metFilterResults.
isValid()) {
1934 if (miniaodFilterIndex_[0] != -1) {
1935 filter_decisions[0] = metFilterResults->
accept(miniaodFilterIndex_[0]);
1937 if (miniaodFilterIndex_[1] != -1) {
1938 filter_decisions[1] = metFilterResults->
accept(miniaodFilterIndex_[1]);
1940 if (miniaodFilterIndex_[2] != -1) {
1941 filter_decisions[2] = metFilterResults->
accept(miniaodFilterIndex_[2]);
1943 if (miniaodFilterIndex_[3] != -1) {
1944 filter_decisions[3] = metFilterResults->
accept(miniaodFilterIndex_[3]);
1946 if (miniaodFilterIndex_[4] != -1) {
1947 filter_decisions[4] = metFilterResults->
accept(miniaodFilterIndex_[4]);
1949 if (miniaodFilterIndex_[5] != -1) {
1950 filter_decisions[5] = metFilterResults->
accept(miniaodFilterIndex_[5]);
1952 if (miniaodFilterIndex_[6] != -1) {
1953 filter_decisions[6] = metFilterResults->
accept(miniaodFilterIndex_[6]);
1955 if (miniaodFilterIndex_[7] != -1) {
1956 filter_decisions[7] = metFilterResults->
accept(miniaodFilterIndex_[7]);
1961 bool HBHENoiseFilterResultFlag = filter_decisions[0];
1964 bool bHBHENoiseFilter = HBHENoiseFilterResultFlag;
1967 bool bDCSFilter = (bypassAllDCSChecks_ || DCSFilter_->filter(iEvent, iSetup));
1971 for (std::vector<std::string>::const_iterator ic = folderNames_.begin(); ic != folderNames_.end(); ic++) {
1972 bool pass_selection =
false;
1973 if ((*ic ==
"Uncleaned") && (isCaloMet_ || bPrimaryVertex)) {
1975 DirName_old +
"/" + *ic,
1984 pass_selection =
true;
1987 if ((*ic ==
"Cleaned") && bDCSFilter && bHBHENoiseFilter && bPrimaryVertex && bJetID) {
1989 DirName_old +
"/" + *ic,
1998 pass_selection =
true;
2000 if ((*ic ==
"DiJet") && bDCSFilter && bHBHENoiseFilter && bPrimaryVertex && bDiJetID) {
2002 DirName_old +
"/" + *ic,
2011 pass_selection =
true;
2013 if ((*ic ==
"ZJets") && bDCSFilter && bHBHENoiseFilter && bPrimaryVertex && bZJets) {
2015 DirName_old +
"/" + *ic,
2024 pass_selection =
true;
2026 if (pass_selection && isPFMet_) {
2027 DirName = DirName_old +
"/" + *ic;
2040 std::map<std::string, MonitorElement*>& map_of_MEs,
2041 std::vector<bool> techTriggerCase,
2042 std::vector<bool> METFilterDecision) {
2043 bool bLumiSecPlot = fill_met_high_level_histo;
2044 bool fillPFCandidatePlots =
false;
2045 if (DirName.find(
"Cleaned") != std::string::npos) {
2046 fillPFCandidatePlots =
true;
2047 fillMonitorElement(iEvent,
2057 fillPFCandidatePlots,
2060 for (
unsigned int i = 0;
i < triggerFolderLabels_.size();
i++) {
2061 fillPFCandidatePlots =
false;
2062 if (triggerFolderDecisions_[
i]) {
2063 fillMonitorElement(iEvent,
2065 triggerFolderLabels_[i],
2073 fillPFCandidatePlots,
2078 }
else if (DirName.find(
"DiJet") != std::string::npos) {
2079 fillMonitorElement(iEvent,
2089 fillPFCandidatePlots,
2092 for (
unsigned int i = 0;
i < triggerFolderLabels_.size();
i++) {
2093 if (triggerFolderDecisions_[
i])
2094 fillMonitorElement(iEvent,
2096 triggerFolderLabels_[i],
2104 fillPFCandidatePlots,
2108 }
else if (DirName.find(
"ZJets") != std::string::npos) {
2109 fillMonitorElement(iEvent,
2119 fillPFCandidatePlots,
2123 fillMonitorElement(iEvent,
2133 fillPFCandidatePlots,
2148 std::map<std::string, MonitorElement*>& map_of_MEs,
2150 bool fillPFCandidatePlots,
2151 std::vector<bool> techTriggerCase,
2152 std::vector<bool> METFilterDecision) {
2153 bool do_only_Z_histograms =
false;
2154 if (DirName.find(
"ZJets") != std::string::npos) {
2155 do_only_Z_histograms =
true;
2158 double u_x = -met.
px() - zCand.Px();
2159 double u_y = -met.
py() - zCand.Py();
2163 double u_perp =
sqrt(u_x * u_x + u_y * u_y);
2166 if (zCand.Pt() != 0) {
2167 e_Z_x = zCand.Px() / zCand.Pt();
2168 e_Z_y = zCand.Py() / zCand.Pt();
2170 u_par = u_x * e_Z_x + u_y * e_Z_y;
2171 u_perp = -e_Z_y * u_x + e_Z_x * u_y;
2173 meZJets_u_par = map_of_MEs[DirName +
"/" +
"u_parallel_Z_inc"];
2174 if (meZJets_u_par && meZJets_u_par->getRootObject())
2175 meZJets_u_par->Fill(u_par);
2176 if (zCand.Pt() < 15) {
2177 meZJets_u_par_ZPt_0_15 = map_of_MEs[DirName +
"/" +
"u_parallel_ZPt_0_15"];
2178 if (meZJets_u_par_ZPt_0_15 && meZJets_u_par_ZPt_0_15->getRootObject())
2179 meZJets_u_par_ZPt_0_15->Fill(u_par);
2180 }
else if (zCand.Pt() < 30) {
2181 meZJets_u_par_ZPt_15_30 = map_of_MEs[DirName +
"/" +
"u_parallel_ZPt_15_30"];
2182 if (meZJets_u_par_ZPt_15_30 && meZJets_u_par_ZPt_15_30->getRootObject())
2183 meZJets_u_par_ZPt_15_30->Fill(u_par);
2184 }
else if (zCand.Pt() < 55) {
2185 meZJets_u_par_ZPt_30_55 = map_of_MEs[DirName +
"/" +
"u_parallel_ZPt_30_55"];
2186 if (meZJets_u_par_ZPt_30_55 && meZJets_u_par_ZPt_30_55->getRootObject())
2187 meZJets_u_par_ZPt_30_55->Fill(u_par);
2188 }
else if (zCand.Pt() < 75) {
2189 meZJets_u_par_ZPt_55_75 = map_of_MEs[DirName +
"/" +
"u_parallel_ZPt_55_75"];
2190 if (meZJets_u_par_ZPt_55_75 && meZJets_u_par_ZPt_55_75->getRootObject())
2191 meZJets_u_par_ZPt_55_75->Fill(u_par);
2192 }
else if (zCand.Pt() < 150) {
2193 meZJets_u_par_ZPt_75_150 = map_of_MEs[DirName +
"/" +
"u_parallel_ZPt_75_150"];
2194 if (meZJets_u_par_ZPt_75_150 && meZJets_u_par_ZPt_75_150->getRootObject())
2195 meZJets_u_par_ZPt_75_150->Fill(u_par);
2196 }
else if (zCand.Pt() < 290) {
2197 meZJets_u_par_ZPt_150_290 = map_of_MEs[DirName +
"/" +
"u_parallel_ZPt_150_290"];
2198 if (meZJets_u_par_ZPt_150_290 && meZJets_u_par_ZPt_150_290->getRootObject())
2199 meZJets_u_par_ZPt_150_290->Fill(u_par);
2201 meZJets_u_par_ZPt_290 = map_of_MEs[DirName +
"/" +
"u_parallel_ZPt_290"];
2202 if (meZJets_u_par_ZPt_290 && meZJets_u_par_ZPt_290->getRootObject())
2203 meZJets_u_par_ZPt_290->Fill(u_par);
2206 meZJets_u_perp = map_of_MEs[DirName +
"/" +
"u_perp_Z_inc"];
2207 if (meZJets_u_perp && meZJets_u_perp->getRootObject())
2208 meZJets_u_perp->Fill(u_perp);
2209 if (zCand.Pt() < 15) {
2210 meZJets_u_perp_ZPt_0_15 = map_of_MEs[DirName +
"/" +
"u_perp_ZPt_0_15"];
2211 if (meZJets_u_perp_ZPt_0_15 && meZJets_u_perp_ZPt_0_15->getRootObject())
2212 meZJets_u_perp_ZPt_0_15->Fill(u_perp);
2213 }
else if (zCand.Pt() < 30) {
2214 meZJets_u_perp_ZPt_15_30 = map_of_MEs[DirName +
"/" +
"u_perp_ZPt_15_30"];
2215 if (meZJets_u_perp_ZPt_15_30 && meZJets_u_perp_ZPt_15_30->getRootObject())
2216 meZJets_u_perp_ZPt_15_30->Fill(u_perp);
2217 }
else if (zCand.Pt() < 55) {
2218 meZJets_u_perp_ZPt_30_55 = map_of_MEs[DirName +
"/" +
"u_perp_ZPt_30_55"];
2219 if (meZJets_u_perp_ZPt_30_55 && meZJets_u_perp_ZPt_30_55->getRootObject())
2220 meZJets_u_perp_ZPt_30_55->Fill(u_perp);
2221 }
else if (zCand.Pt() < 75) {
2222 meZJets_u_perp_ZPt_55_75 = map_of_MEs[DirName +
"/" +
"u_perp_ZPt_55_75"];
2223 if (meZJets_u_perp_ZPt_55_75 && meZJets_u_perp_ZPt_55_75->getRootObject())
2224 meZJets_u_perp_ZPt_55_75->Fill(u_perp);
2225 }
else if (zCand.Pt() < 150) {
2226 meZJets_u_perp_ZPt_75_150 = map_of_MEs[DirName +
"/" +
"u_perp_ZPt_75_150"];
2227 if (meZJets_u_perp_ZPt_75_150 && meZJets_u_perp_ZPt_75_150->getRootObject())
2228 meZJets_u_perp_ZPt_75_150->Fill(u_perp);
2229 }
else if (zCand.Pt() < 290) {
2230 meZJets_u_perp_ZPt_150_290 = map_of_MEs[DirName +
"/" +
"u_perp_ZPt_150_290"];
2231 if (meZJets_u_perp_ZPt_150_290 && meZJets_u_perp_ZPt_150_290->getRootObject())
2232 meZJets_u_perp_ZPt_150_290->Fill(u_perp);
2234 meZJets_u_perp_ZPt_290 = map_of_MEs[DirName +
"/" +
"u_perp_ZPt_290"];
2235 if (meZJets_u_perp_ZPt_290 && meZJets_u_perp_ZPt_290->getRootObject())
2236 meZJets_u_perp_ZPt_290->Fill(u_perp);
2239 if (!do_only_Z_histograms) {
2241 double SumET = met.
sumEt();
2242 double METSig = met.
mEtSig();
2244 double MET = met.
pt();
2245 double MEx = met.
px();
2246 double MEy = met.
py();
2247 double METPhi = met.
phi();
2249 int myLuminosityBlock;
2253 if (!subFolderName.empty()) {
2254 DirName = DirName +
"/" + subFolderName;
2257 hTrigger = map_of_MEs[DirName +
"/triggerResults"];
2258 if (hTrigger && hTrigger->getRootObject()) {
2259 for (
unsigned int i = 0;
i < allTriggerDecisions_.size();
i++) {
2260 if (
i < (
unsigned int)hTrigger->getNbinsX()) {
2261 hTrigger->Fill(
i + .5, allTriggerDecisions_[
i]);
2266 hMEx = map_of_MEs[DirName +
"/" +
"MEx"];
2267 if (hMEx && hMEx->getRootObject())
2269 hMEy = map_of_MEs[DirName +
"/" +
"MEy"];
2270 if (hMEy && hMEy->getRootObject())
2272 hMET = map_of_MEs[DirName +
"/" +
"MET"];
2273 if (hMET && hMET->getRootObject())
2275 hMET_2 = map_of_MEs[DirName +
"/" +
"MET_2"];
2276 if (hMET_2 && hMET_2->getRootObject())
2288 bool HBHENoiseFilterResult =
false;
2289 bool CSCTightHaloFilterResult =
false;
2290 bool eeBadScFilterResult =
false;
2291 bool HBHEIsoNoiseFilterResult =
false;
2292 bool CSCTightHalo2015FilterResult =
false;
2293 bool EcalDeadCellTriggerFilterResult =
false;
2294 bool EcalDeadCellBoundaryFilterResult =
false;
2295 bool HcalStripHaloFilterResult =
false;
2296 HBHENoiseFilterResult = METFilterDecision[0];
2297 if (HBHENoiseFilterResult) {
2298 hMET_HBHENoiseFilter = map_of_MEs[DirName +
"/" +
"MET_HBHENoiseFilter"];
2299 if (hMET_HBHENoiseFilter && hMET_HBHENoiseFilter->getRootObject())
2300 hMET_HBHENoiseFilter->Fill(MET);
2302 CSCTightHaloFilterResult = METFilterDecision[1];
2303 if (CSCTightHaloFilterResult) {
2304 hMET_CSCTightHaloFilter = map_of_MEs[DirName +
"/" +
"MET_CSCTightHaloFilter"];
2305 if (hMET_CSCTightHaloFilter && hMET_CSCTightHaloFilter->getRootObject())
2306 hMET_CSCTightHaloFilter->Fill(MET);
2308 eeBadScFilterResult = METFilterDecision[2];
2309 if (eeBadScFilterResult) {
2310 hMET_eeBadScFilter = map_of_MEs[DirName +
"/" +
"MET_eeBadScFilter"];
2311 if (hMET_eeBadScFilter && hMET_eeBadScFilter->getRootObject())
2312 hMET_eeBadScFilter->Fill(MET);
2314 HBHEIsoNoiseFilterResult = METFilterDecision[3];
2315 if (HBHEIsoNoiseFilterResult) {
2316 hMET_HBHEIsoNoiseFilter = map_of_MEs[DirName +
"/" +
"MET_HBHEIsoNoiseFilter"];
2317 if (hMET_HBHEIsoNoiseFilter && hMET_HBHEIsoNoiseFilter->getRootObject())
2318 hMET_HBHEIsoNoiseFilter->Fill(MET);
2320 CSCTightHalo2015FilterResult = METFilterDecision[4];
2321 if (CSCTightHalo2015FilterResult) {
2322 hMET_CSCTightHalo2015Filter = map_of_MEs[DirName +
"/" +
"MET_CSCTightHalo2015Filter"];
2323 if (hMET_CSCTightHalo2015Filter && hMET_CSCTightHalo2015Filter->getRootObject())
2324 hMET_CSCTightHalo2015Filter->Fill(MET);
2326 EcalDeadCellTriggerFilterResult = METFilterDecision[5];
2327 if (EcalDeadCellTriggerFilterResult) {
2328 hMET_EcalDeadCellTriggerFilter = map_of_MEs[DirName +
"/" +
"MET_EcalDeadCellTriggerFilter"];
2329 if (hMET_EcalDeadCellTriggerFilter && hMET_EcalDeadCellTriggerFilter->getRootObject())
2330 hMET_EcalDeadCellTriggerFilter->Fill(MET);
2332 EcalDeadCellBoundaryFilterResult = METFilterDecision[6];
2333 if (EcalDeadCellBoundaryFilterResult) {
2334 hMET_EcalDeadCellBoundaryFilter = map_of_MEs[DirName +
"/" +
"MET_EcalDeadCellBoundaryFilter"];
2335 if (hMET_EcalDeadCellBoundaryFilter && hMET_EcalDeadCellBoundaryFilter->getRootObject())
2336 hMET_EcalDeadCellBoundaryFilter->Fill(MET);
2338 HcalStripHaloFilterResult = METFilterDecision[7];
2339 if (HcalStripHaloFilterResult) {
2340 hMET_HcalStripHaloFilter = map_of_MEs[DirName +
"/" +
"MET_HcalStripHaloFilter"];
2341 if (hMET_HcalStripHaloFilter && hMET_HcalStripHaloFilter->getRootObject())
2342 hMET_HcalStripHaloFilter->Fill(MET);
2344 hMETPhi = map_of_MEs[DirName +
"/" +
"METPhi"];
2345 if (hMETPhi && hMETPhi->getRootObject())
2346 hMETPhi->Fill(METPhi);
2347 hSumET = map_of_MEs[DirName +
"/" +
"SumET"];
2348 if (hSumET && hSumET->getRootObject())
2349 hSumET->Fill(SumET);
2350 hMETSig = map_of_MEs[DirName +
"/" +
"METSig"];
2351 if (hMETSig && hMETSig->getRootObject())
2352 hMETSig->Fill(METSig);
2353 hMET_logx = map_of_MEs[DirName +
"/" +
"MET_logx"];
2354 if (hMET_logx && hMET_logx->getRootObject())
2355 hMET_logx->Fill(log10(MET));
2356 hSumET_logx = map_of_MEs[DirName +
"/" +
"SumET_logx"];
2357 if (hSumET_logx && hSumET_logx->getRootObject())
2358 hSumET_logx->Fill(log10(SumET));
2362 meMEx_profile = map_of_MEs[DirName +
"/MEx_profile"];
2363 meMEy_profile = map_of_MEs[DirName +
"/MEy_profile"];
2364 meMET_profile = map_of_MEs[DirName +
"/MET_profile"];
2365 meSumET_profile = map_of_MEs[DirName +
"/SumET_profile"];
2367 if (meMEx_profile && meMEx_profile->getRootObject())
2368 meMEx_profile->Fill(numPV_, MEx);
2369 if (meMEy_profile && meMEy_profile->getRootObject())
2370 meMEy_profile->Fill(numPV_, MEy);
2371 if (meMET_profile && meMET_profile->getRootObject())
2372 meMET_profile->Fill(numPV_, MET);
2373 if (meSumET_profile && meSumET_profile->getRootObject())
2374 meSumET_profile->Fill(numPV_, SumET);
2384 double caloHadEtInHB = calomet->
hadEtInHB();
2385 double caloHadEtInHO = calomet->
hadEtInHO();
2386 double caloHadEtInHE = calomet->
hadEtInHE();
2387 double caloHadEtInHF = calomet->
hadEtInHF();
2388 double caloEmEtInEB = calomet->
emEtInEB();
2389 double caloEmEtInEE = calomet->
emEtInEE();
2390 double caloEmEtInHF = calomet->
emEtInHF();
2392 hCaloHadEtInHB = map_of_MEs[DirName +
"/" +
"CaloHadEtInHB"];
2393 if (hCaloHadEtInHB && hCaloHadEtInHB->getRootObject())
2394 hCaloHadEtInHB->Fill(caloHadEtInHB);
2395 hCaloHadEtInHO = map_of_MEs[DirName +
"/" +
"CaloHadEtInHO"];
2396 if (hCaloHadEtInHO && hCaloHadEtInHO->getRootObject())
2397 hCaloHadEtInHO->Fill(caloHadEtInHO);
2398 hCaloHadEtInHE = map_of_MEs[DirName +
"/" +
"CaloHadEtInHE"];
2399 if (hCaloHadEtInHE && hCaloHadEtInHE->getRootObject())
2400 hCaloHadEtInHE->Fill(caloHadEtInHE);
2401 hCaloHadEtInHF = map_of_MEs[DirName +
"/" +
"CaloHadEtInHF"];
2402 if (hCaloHadEtInHF && hCaloHadEtInHF->getRootObject())
2403 hCaloHadEtInHF->Fill(caloHadEtInHF);
2404 hCaloEmEtInEB = map_of_MEs[DirName +
"/" +
"CaloEmEtInEB"];
2405 if (hCaloEmEtInEB && hCaloEmEtInEB->getRootObject())
2406 hCaloEmEtInEB->Fill(caloEmEtInEB);
2407 hCaloEmEtInEE = map_of_MEs[DirName +
"/" +
"CaloEmEtInEE"];
2408 if (hCaloEmEtInEE && hCaloEmEtInEE->getRootObject())
2409 hCaloEmEtInEE->Fill(caloEmEtInEE);
2410 hCaloEmEtInHF = map_of_MEs[DirName +
"/" +
"CaloEmEtInHF"];
2411 if (hCaloEmEtInHF && hCaloEmEtInHF->getRootObject())
2412 hCaloEmEtInHF->Fill(caloEmEtInHF);
2414 hCaloMETPhi020 = map_of_MEs[DirName +
"/" +
"CaloMETPhi020"];
2415 if (MET > 20. && hCaloMETPhi020 && hCaloMETPhi020->getRootObject()) {
2416 hCaloMETPhi020->Fill(METPhi);
2419 hCaloEtFractionHadronic = map_of_MEs[DirName +
"/" +
"CaloEtFractionHadronic"];
2420 if (hCaloEtFractionHadronic && hCaloEtFractionHadronic->getRootObject())
2421 hCaloEtFractionHadronic->Fill(caloEtFractionHadronic);
2422 hCaloEmEtFraction = map_of_MEs[DirName +
"/" +
"CaloEmEtFraction"];
2423 if (hCaloEmEtFraction && hCaloEmEtFraction->getRootObject())
2424 hCaloEmEtFraction->Fill(caloEmEtFraction);
2425 hCaloEmEtFraction020 = map_of_MEs[DirName +
"/" +
"CaloEmEtFraction020"];
2426 if (MET > 20. && hCaloEmEtFraction020 && hCaloEmEtFraction020->getRootObject())
2427 hCaloEmEtFraction020->Fill(caloEmEtFraction);
2430 if (fillPFCandidatePlots && fillCandidateMap_histos) {
2431 for (
unsigned int i = 0;
i < countsPFCand_.size();
i++) {
2432 countsPFCand_[
i] = 0;
2439 iEvent.
getByToken(pflowToken_, particleFlow);
2441 float pt_sum_CHF_Barrel = 0;
2442 float pt_sum_CHF_Endcap_plus = 0;
2443 float pt_sum_CHF_Endcap_minus = 0;
2444 float pt_sum_NHF_Barrel = 0;
2445 float pt_sum_NHF_Endcap_plus = 0;
2446 float pt_sum_NHF_Endcap_minus = 0;
2447 float pt_sum_PhF_Barrel = 0;
2448 float pt_sum_PhF_Endcap_plus = 0;
2449 float pt_sum_PhF_Endcap_minus = 0;
2450 float pt_sum_HFH_plus = 0;
2451 float pt_sum_HFH_minus = 0;
2452 float pt_sum_HFE_plus = 0;
2453 float pt_sum_HFE_minus = 0;
2455 float px_chargedHadronsBarrel = 0;
2456 float py_chargedHadronsBarrel = 0;
2457 float px_chargedHadronsEndcapPlus = 0;
2458 float py_chargedHadronsEndcapPlus = 0;
2459 float px_chargedHadronsEndcapMinus = 0;
2460 float py_chargedHadronsEndcapMinus = 0;
2461 float px_neutralHadronsBarrel = 0;
2462 float py_neutralHadronsBarrel = 0;
2463 float px_neutralHadronsEndcapPlus = 0;
2464 float py_neutralHadronsEndcapPlus = 0;
2465 float px_neutralHadronsEndcapMinus = 0;
2466 float py_neutralHadronsEndcapMinus = 0;
2467 float px_PhotonsBarrel = 0;
2468 float py_PhotonsBarrel = 0;
2469 float px_PhotonsEndcapPlus = 0;
2470 float py_PhotonsEndcapPlus = 0;
2471 float px_PhotonsEndcapMinus = 0;
2472 float py_PhotonsEndcapMinus = 0;
2473 float px_HFHadronsPlus = 0;
2474 float py_HFHadronsPlus = 0;
2475 float px_HFHadronsMinus = 0;
2476 float py_HFHadronsMinus = 0;
2477 float px_HFEGammasPlus = 0;
2478 float py_HFEGammasPlus = 0;
2479 float px_HFEGammasMinus = 0;
2480 float py_HFEGammasMinus = 0;
2481 for (
unsigned int i = 0;
i < particleFlow->size();
i++) {
2485 if (c.
eta() > (-3.0) && c.
eta() < (-1.392)) {
2486 px_chargedHadronsEndcapMinus -= c.
px();
2487 py_chargedHadronsEndcapMinus -= c.
py();
2488 pt_sum_CHF_Endcap_minus += c.
et();
2489 }
else if (c.
eta() >= (-1.392) && c.
eta() <= 1.392) {
2490 px_chargedHadronsBarrel -= c.
px();
2491 py_chargedHadronsBarrel -= c.
py();
2492 pt_sum_CHF_Barrel += c.
et();
2493 }
else if (c.
eta() > 1.392 && c.
eta() < 3.0) {
2494 px_chargedHadronsEndcapPlus -= c.
px();
2495 py_chargedHadronsEndcapPlus -= c.
py();
2496 pt_sum_CHF_Endcap_plus += c.
et();
2501 if (c.
eta() > (-3.0) && c.
eta() < (-1.392)) {
2502 px_neutralHadronsEndcapMinus -= c.
px();
2503 py_neutralHadronsEndcapMinus -= c.
py();
2504 pt_sum_NHF_Endcap_minus += c.
et();
2505 }
else if (c.
eta() >= (-1.392) && c.
eta() <= 1.392) {
2506 px_neutralHadronsBarrel -= c.
px();
2507 py_neutralHadronsBarrel -= c.
py();
2508 pt_sum_NHF_Barrel += c.
et();
2509 }
else if (c.
eta() > 1.392 && c.
eta() < 3.0) {
2510 px_neutralHadronsEndcapPlus -= c.
px();
2511 py_neutralHadronsEndcapPlus -= c.
py();
2512 pt_sum_NHF_Endcap_plus += c.
et();
2517 if (c.
eta() > (-3.0) && c.
eta() < (-1.479)) {
2518 px_PhotonsEndcapMinus -= c.
px();
2519 py_PhotonsEndcapMinus -= c.
py();
2520 pt_sum_PhF_Endcap_minus += c.
et();
2521 }
else if (c.
eta() >= (-1.479) && c.
eta() <= 1.479) {
2522 px_PhotonsBarrel -= c.
px();
2523 py_PhotonsBarrel -= c.
py();
2524 pt_sum_PhF_Barrel += c.
et();
2525 }
else if (c.
eta() > 1.479 && c.
eta() < 3.0) {
2526 px_PhotonsEndcapPlus -= c.
px();
2527 py_PhotonsEndcapPlus -= c.
py();
2528 pt_sum_PhF_Endcap_plus += c.
et();
2533 if (c.
eta() > (-5.20) && c.
eta() < -2.901376) {
2534 pt_sum_HFH_minus += c.
et();
2535 px_HFHadronsMinus -= c.
px();
2536 py_HFHadronsMinus -= c.
py();
2537 }
else if (c.
eta() > 2.901376 && c.
eta() < 5.20) {
2538 px_HFHadronsPlus -= c.
px();
2539 py_HFHadronsPlus -= c.
py();
2540 pt_sum_HFH_plus += c.
et();
2545 if (c.
eta() > (-5.20) && c.
eta() < -2.901376) {
2546 pt_sum_HFE_minus += c.
et();
2547 px_HFEGammasMinus -= c.
px();
2548 py_HFEGammasMinus -= c.
py();
2549 }
else if (c.
eta() > 2.901376 && c.
eta() < 5.20) {
2550 px_HFEGammasPlus -= c.
px();
2551 py_HFEGammasPlus -= c.
py();
2552 pt_sum_HFE_plus += c.
et();
2555 for (
unsigned int j = 0;
j < typePFCand_.size();
j++) {
2558 if (((c.
eta() > etaMinPFCand_[
j]) && (c.
eta() < etaMaxPFCand_[
j])) ||
2559 ((c.
eta() > (-etaMaxPFCand_[
j])) && (c.
eta() < (-etaMinPFCand_[
j])))) {
2560 countsPFCand_[
j] += 1;
2561 MExPFCand_[
j] -= c.
px();
2562 MEyPFCand_[
j] -= c.
py();
2568 for (
unsigned int j = 0;
j < countsPFCand_.size();
j++) {
2569 profilePFCand_x_[
j] = map_of_MEs[DirName +
"/" + profilePFCand_x_name_[
j]];
2570 if (profilePFCand_x_[
j] && profilePFCand_x_[
j]->getRootObject())
2571 profilePFCand_x_[
j]->
Fill(countsPFCand_[
j], MExPFCand_[j]);
2572 profilePFCand_y_[
j] = map_of_MEs[DirName +
"/" + profilePFCand_y_name_[
j]];
2573 if (profilePFCand_y_[j] && profilePFCand_y_[j]->getRootObject())
2574 profilePFCand_y_[j]->
Fill(countsPFCand_[j], MEyPFCand_[j]);
2576 meCHF_Barrel = map_of_MEs[DirName +
"/" +
"PfChargedHadronEtFractionBarrel"];
2577 if (meCHF_Barrel && meCHF_Barrel->getRootObject())
2578 meCHF_Barrel->Fill(pt_sum_CHF_Barrel / pfmet->
sumEt());
2579 meCHF_EndcapPlus = map_of_MEs[DirName +
"/" +
"PfChargedHadronEtFractionEndcapPlus"];
2580 if (meCHF_EndcapPlus && meCHF_EndcapPlus->getRootObject())
2581 meCHF_EndcapPlus->Fill(pt_sum_CHF_Endcap_plus / pfmet->
sumEt());
2582 meCHF_EndcapMinus = map_of_MEs[DirName +
"/" +
"PfChargedHadronEtFractionEndcapMinus"];
2583 if (meCHF_EndcapMinus && meCHF_EndcapMinus->getRootObject())
2584 meCHF_EndcapMinus->Fill(pt_sum_CHF_Endcap_minus / pfmet->
sumEt());
2585 meNHF_Barrel = map_of_MEs[DirName +
"/" +
"PfNeutralHadronEtFractionBarrel"];
2586 if (meNHF_Barrel && meNHF_Barrel->getRootObject())
2587 meNHF_Barrel->Fill(pt_sum_NHF_Barrel / pfmet->
sumEt());
2588 meNHF_EndcapPlus = map_of_MEs[DirName +
"/" +
"PfNeutralHadronEtFractionEndcapPlus"];
2589 if (meNHF_EndcapPlus && meNHF_EndcapPlus->getRootObject())
2590 meNHF_EndcapPlus->Fill(pt_sum_NHF_Endcap_plus / pfmet->
sumEt());
2591 meNHF_EndcapMinus = map_of_MEs[DirName +
"/" +
"PfNeutralHadronEtFractionEndcapMinus"];
2592 if (meNHF_EndcapMinus && meNHF_EndcapMinus->getRootObject())
2593 meNHF_EndcapMinus->Fill(pt_sum_NHF_Endcap_minus / pfmet->
sumEt());
2594 mePhF_Barrel = map_of_MEs[DirName +
"/" +
"PfPhotonEtFractionBarrel"];
2595 if (mePhF_Barrel && mePhF_Barrel->getRootObject())
2596 mePhF_Barrel->Fill(pt_sum_PhF_Barrel / pfmet->
sumEt());
2597 mePhF_EndcapPlus = map_of_MEs[DirName +
"/" +
"PfPhotonEtFractionEndcapPlus"];
2598 if (mePhF_EndcapPlus && mePhF_EndcapPlus->getRootObject())
2599 mePhF_EndcapPlus->Fill(pt_sum_PhF_Endcap_plus / pfmet->
sumEt());
2600 mePhF_EndcapMinus = map_of_MEs[DirName +
"/" +
"PfPhotonEtFractionEndcapMinus"];
2601 if (mePhF_EndcapMinus && mePhF_EndcapMinus->getRootObject())
2602 mePhF_EndcapMinus->Fill(pt_sum_PhF_Endcap_minus / pfmet->
sumEt());
2603 meHFHadF_Plus = map_of_MEs[DirName +
"/" +
"PfHFHadronEtFractionPlus"];
2604 if (meHFHadF_Plus && meHFHadF_Plus->getRootObject())
2605 meHFHadF_Plus->Fill(pt_sum_HFH_plus / pfmet->
sumEt());
2606 meHFHadF_Minus = map_of_MEs[DirName +
"/" +
"PfHFHadronEtFractionMinus"];
2607 if (meHFHadF_Minus && meHFHadF_Minus->getRootObject())
2608 meHFHadF_Minus->Fill(pt_sum_HFH_minus / pfmet->
sumEt());
2609 meHFEMF_Plus = map_of_MEs[DirName +
"/" +
"PfHFEMEtFractionPlus"];
2610 if (meHFEMF_Plus && meHFEMF_Plus->getRootObject())
2611 meHFEMF_Plus->Fill(pt_sum_HFE_plus / pfmet->
sumEt());
2612 meHFEMF_Minus = map_of_MEs[DirName +
"/" +
"PfHFEMEtFractionMinus"];
2613 if (meHFEMF_Minus && meHFEMF_Minus->getRootObject())
2614 meHFEMF_Minus->Fill(pt_sum_HFE_minus / pfmet->
sumEt());
2618 if (pt_sum_CHF_Barrel) {
2619 meMETPhiChargedHadronsBarrel = map_of_MEs[DirName +
"/" +
"METPhiChargedHadronsBarrel"];
2620 if (meMETPhiChargedHadronsBarrel && meMETPhiChargedHadronsBarrel->getRootObject())
2621 meMETPhiChargedHadronsBarrel->Fill(atan2(py_chargedHadronsBarrel, px_chargedHadronsBarrel));
2623 if (pt_sum_CHF_Endcap_plus) {
2624 meMETPhiChargedHadronsEndcapPlus = map_of_MEs[DirName +
"/" +
"METPhiChargedHadronsEndcapPlus"];
2625 if (meMETPhiChargedHadronsEndcapPlus && meMETPhiChargedHadronsEndcapPlus->getRootObject())
2626 meMETPhiChargedHadronsEndcapPlus->Fill(atan2(py_chargedHadronsEndcapPlus, px_chargedHadronsEndcapPlus));
2628 if (pt_sum_CHF_Endcap_minus) {
2629 meMETPhiChargedHadronsEndcapMinus = map_of_MEs[DirName +
"/" +
"METPhiChargedHadronsEndcapMinus"];
2630 if (meMETPhiChargedHadronsEndcapMinus && meMETPhiChargedHadronsEndcapMinus->getRootObject())
2631 meMETPhiChargedHadronsEndcapMinus->Fill(atan2(py_chargedHadronsEndcapMinus, px_chargedHadronsEndcapMinus));
2633 if (pt_sum_NHF_Barrel) {
2634 meMETPhiNeutralHadronsBarrel = map_of_MEs[DirName +
"/" +
"METPhiNeutralHadronsBarrel"];
2635 if (meMETPhiNeutralHadronsBarrel && meMETPhiNeutralHadronsBarrel->getRootObject())
2636 meMETPhiNeutralHadronsBarrel->Fill(atan2(py_neutralHadronsBarrel, px_neutralHadronsBarrel));
2638 if (pt_sum_NHF_Endcap_plus) {
2639 meMETPhiNeutralHadronsEndcapPlus = map_of_MEs[DirName +
"/" +
"METPhiNeutralHadronsEndcapPlus"];
2640 if (meMETPhiNeutralHadronsEndcapPlus && meMETPhiNeutralHadronsEndcapPlus->getRootObject())
2641 meMETPhiNeutralHadronsEndcapPlus->Fill(atan2(py_neutralHadronsEndcapPlus, px_neutralHadronsEndcapPlus));
2643 if (pt_sum_NHF_Endcap_minus) {
2644 meMETPhiNeutralHadronsEndcapMinus = map_of_MEs[DirName +
"/" +
"METPhiNeutralHadronsEndcapMinus"];
2645 if (meMETPhiNeutralHadronsEndcapMinus && meMETPhiNeutralHadronsEndcapMinus->getRootObject())
2646 meMETPhiNeutralHadronsEndcapMinus->Fill(atan2(py_neutralHadronsEndcapMinus, px_neutralHadronsEndcapMinus));
2648 if (pt_sum_PhF_Barrel) {
2649 meMETPhiPhotonsBarrel = map_of_MEs[DirName +
"/" +
"METPhiPhotonsBarrel"];
2650 if (meMETPhiPhotonsBarrel && meMETPhiPhotonsBarrel->getRootObject())
2651 meMETPhiPhotonsBarrel->Fill(atan2(py_PhotonsBarrel, px_PhotonsBarrel));
2653 if (pt_sum_PhF_Endcap_plus) {
2654 meMETPhiPhotonsEndcapPlus = map_of_MEs[DirName +
"/" +
"METPhiPhotonsEndcapPlus"];
2655 if (meMETPhiPhotonsEndcapPlus && meMETPhiPhotonsEndcapPlus->getRootObject())
2656 meMETPhiPhotonsEndcapPlus->Fill(atan2(py_PhotonsEndcapPlus, px_PhotonsEndcapPlus));
2658 if (pt_sum_PhF_Endcap_minus) {
2659 meMETPhiPhotonsEndcapMinus = map_of_MEs[DirName +
"/" +
"METPhiPhotonsEndcapMinus"];
2660 if (meMETPhiPhotonsEndcapMinus && meMETPhiPhotonsEndcapMinus->getRootObject())
2661 meMETPhiPhotonsEndcapMinus->Fill(atan2(py_PhotonsEndcapMinus, px_PhotonsEndcapMinus));
2663 if (pt_sum_HFH_plus) {
2664 meMETPhiHFHadronsPlus = map_of_MEs[DirName +
"/" +
"METPhiHFHadronsPlus"];
2665 if (meMETPhiHFHadronsPlus && meMETPhiHFHadronsPlus->getRootObject())
2666 meMETPhiHFHadronsPlus->Fill(atan2(py_HFHadronsPlus, px_HFHadronsPlus));
2668 if (pt_sum_HFH_minus) {
2669 meMETPhiHFHadronsMinus = map_of_MEs[DirName +
"/" +
"METPhiHFHadronsMinus"];
2670 if (meMETPhiHFHadronsMinus && meMETPhiHFHadronsMinus->getRootObject())
2671 meMETPhiHFHadronsMinus->Fill(atan2(py_HFHadronsMinus, px_HFHadronsMinus));
2673 if (pt_sum_HFE_plus) {
2674 meMETPhiHFEGammasPlus = map_of_MEs[DirName +
"/" +
"METPhiHFEGammasPlus"];
2675 if (meMETPhiHFEGammasPlus && meMETPhiHFEGammasPlus->getRootObject())
2676 meMETPhiHFEGammasPlus->Fill(atan2(py_HFEGammasPlus, px_HFEGammasPlus));
2678 if (pt_sum_HFE_minus) {
2679 meMETPhiHFEGammasMinus = map_of_MEs[DirName +
"/" +
"METPhiHFEGammasMinus"];
2680 if (meMETPhiHFEGammasMinus && meMETPhiHFEGammasMinus->getRootObject())
2681 meMETPhiHFEGammasMinus->Fill(atan2(py_HFEGammasMinus, px_HFEGammasMinus));
2743 if (techTriggerCase[1]) {
2744 meCHF_Barrel_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfChargedHadronEtFractionBarrel_BXm1Filled"];
2745 if (meCHF_Barrel_BXm1Filled && meCHF_Barrel_BXm1Filled->getRootObject())
2746 meCHF_Barrel_BXm1Filled->Fill(pt_sum_CHF_Barrel / pfmet->
sumEt());
2747 meCHF_EndcapPlus_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfChargedHadronEtFractionEndcapPlus_BXm1Filled"];
2748 if (meCHF_EndcapPlus_BXm1Filled && meCHF_EndcapPlus_BXm1Filled->getRootObject())
2749 meCHF_EndcapPlus_BXm1Filled->Fill(pt_sum_CHF_Endcap_plus / pfmet->
sumEt());
2750 meCHF_EndcapMinus_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfChargedHadronEtFractionEndcapMinus_BXm1Filled"];
2751 if (meCHF_EndcapMinus_BXm1Filled && meCHF_EndcapMinus_BXm1Filled->getRootObject())
2752 meCHF_EndcapMinus_BXm1Filled->Fill(pt_sum_CHF_Endcap_minus / pfmet->
sumEt());
2753 meNHF_Barrel_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfNeutralHadronEtFractionBarrel_BXm1Filled"];
2754 if (meNHF_Barrel_BXm1Filled && meNHF_Barrel_BXm1Filled->getRootObject())
2755 meNHF_Barrel_BXm1Filled->Fill(pt_sum_NHF_Barrel / pfmet->
sumEt());
2756 meNHF_EndcapPlus_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfNeutralHadronEtFractionEndcapPlus_BXm1Filled"];
2757 if (meNHF_EndcapPlus_BXm1Filled && meNHF_EndcapPlus_BXm1Filled->getRootObject())
2758 meNHF_EndcapPlus_BXm1Filled->Fill(pt_sum_NHF_Endcap_plus / pfmet->
sumEt());
2759 meNHF_EndcapMinus_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfNeutralHadronEtFractionEndcapMinus_BXm1Filled"];
2760 if (meNHF_EndcapMinus_BXm1Filled && meNHF_EndcapMinus_BXm1Filled->getRootObject())
2761 meNHF_EndcapMinus_BXm1Filled->Fill(pt_sum_NHF_Endcap_minus / pfmet->
sumEt());
2762 mePhF_Barrel_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfPhotonEtFractionBarrel_BXm1Filled"];
2763 if (mePhF_Barrel_BXm1Filled && mePhF_Barrel_BXm1Filled->getRootObject())
2764 mePhF_Barrel_BXm1Filled->Fill(pt_sum_PhF_Barrel / pfmet->
sumEt());
2765 mePhF_EndcapPlus_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfPhotonEtFractionEndcapPlus_BXm1Filled"];
2766 if (mePhF_EndcapPlus_BXm1Filled && mePhF_EndcapPlus_BXm1Filled->getRootObject())
2767 mePhF_EndcapPlus_BXm1Filled->Fill(pt_sum_PhF_Endcap_plus / pfmet->
sumEt());
2768 mePhF_EndcapMinus_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfPhotonEtFractionEndcapMinus_BXm1Filled"];
2769 if (mePhF_EndcapMinus_BXm1Filled && mePhF_EndcapMinus_BXm1Filled->getRootObject())
2770 mePhF_EndcapMinus_BXm1Filled->Fill(pt_sum_PhF_Endcap_minus / pfmet->
sumEt());
2771 meHFHadF_Plus_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfHFHadronEtFractionPlus_BXm1Filled"];
2772 if (meHFHadF_Plus_BXm1Filled && meHFHadF_Plus_BXm1Filled->getRootObject())
2773 meHFHadF_Plus_BXm1Filled->Fill(pt_sum_HFH_plus / pfmet->
sumEt());
2774 meHFHadF_Minus_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfHFHadronEtFractionMinus_BXm1Filled"];
2775 if (meHFHadF_Minus_BXm1Filled && meHFHadF_Minus_BXm1Filled->getRootObject())
2776 meHFHadF_Minus_BXm1Filled->Fill(pt_sum_HFH_minus / pfmet->
sumEt());
2777 meHFEMF_Plus_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfHFEMEtFractionPlus_BXm1Filled"];
2778 if (meHFEMF_Plus_BXm1Filled && meHFEMF_Plus_BXm1Filled->getRootObject())
2779 meHFEMF_Plus_BXm1Filled->Fill(pt_sum_HFE_plus / pfmet->
sumEt());
2780 meHFEMF_Minus_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfHFEMEtFractionMinus_BXm1Filled"];
2781 if (meHFEMF_Minus_BXm1Filled && meHFEMF_Minus_BXm1Filled->getRootObject())
2782 meHFEMF_Minus_BXm1Filled->Fill(pt_sum_HFE_minus / pfmet->
sumEt());
2783 mePhotonEtFraction_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfPhotonEtFraction_BXm1Filled"];
2784 if (mePhotonEtFraction_BXm1Filled && mePhotonEtFraction_BXm1Filled->getRootObject())
2786 meNeutralHadronEtFraction_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfNeutralHadronEtFraction_BXm1Filled"];
2787 if (meNeutralHadronEtFraction_BXm1Filled && meNeutralHadronEtFraction_BXm1Filled->getRootObject())
2789 meChargedHadronEtFraction_BXm1Filled = map_of_MEs[DirName +
"/" +
"PfChargedHadronEtFraction_BXm1Filled"];
2790 if (meChargedHadronEtFraction_BXm1Filled && meChargedHadronEtFraction_BXm1Filled->getRootObject())
2792 meMET_BXm1Filled = map_of_MEs[DirName +
"/" +
"MET_BXm1Filled"];
2793 if (meMET_BXm1Filled && meMET_BXm1Filled->getRootObject())
2794 meMET_BXm1Filled->Fill(pfmet->
pt());
2795 meSumET_BXm1Filled = map_of_MEs[DirName +
"/" +
"SumET_BXm1Filled"];
2796 if (meSumET_BXm1Filled && meSumET_BXm1Filled->getRootObject())
2797 meSumET_BXm1Filled->Fill(pfmet->
sumEt());
2798 if (pt_sum_CHF_Barrel) {
2799 meMETPhiChargedHadronsBarrel_BXm1Filled =
2800 map_of_MEs[DirName +
"/" +
"METPhiChargedHadronsBarrel_BXm1Filled"];
2801 if (meMETPhiChargedHadronsBarrel_BXm1Filled && meMETPhiChargedHadronsBarrel_BXm1Filled->getRootObject())
2802 meMETPhiChargedHadronsBarrel_BXm1Filled->Fill(atan2(py_chargedHadronsBarrel, px_chargedHadronsBarrel));
2804 if (pt_sum_CHF_Endcap_plus) {
2805 meMETPhiChargedHadronsEndcapPlus_BXm1Filled =
2806 map_of_MEs[DirName +
"/" +
"METPhiChargedHadronsEndcapPlus_BXm1Filled"];
2807 if (meMETPhiChargedHadronsEndcapPlus_BXm1Filled &&
2808 meMETPhiChargedHadronsEndcapPlus_BXm1Filled->getRootObject())
2809 meMETPhiChargedHadronsEndcapPlus_BXm1Filled->Fill(
2810 atan2(py_chargedHadronsEndcapPlus, px_chargedHadronsEndcapPlus));
2812 if (pt_sum_CHF_Endcap_minus) {
2813 meMETPhiChargedHadronsEndcapMinus_BXm1Filled =
2814 map_of_MEs[DirName +
"/" +
"METPhiChargedHadronsEndcapMinus_BXm1Filled"];
2815 if (meMETPhiChargedHadronsEndcapMinus_BXm1Filled &&
2816 meMETPhiChargedHadronsEndcapMinus_BXm1Filled->getRootObject())
2817 meMETPhiChargedHadronsEndcapMinus_BXm1Filled->Fill(
2818 atan2(py_chargedHadronsEndcapMinus, px_chargedHadronsEndcapMinus));
2820 if (pt_sum_NHF_Barrel) {
2821 meMETPhiNeutralHadronsBarrel_BXm1Filled =
2822 map_of_MEs[DirName +
"/" +
"METPhiNeutralHadronsBarrel_BXm1Filled"];
2823 if (meMETPhiNeutralHadronsBarrel_BXm1Filled && meMETPhiNeutralHadronsBarrel_BXm1Filled->getRootObject())
2824 meMETPhiNeutralHadronsBarrel_BXm1Filled->Fill(atan2(py_neutralHadronsBarrel, px_neutralHadronsBarrel));
2826 if (pt_sum_NHF_Endcap_plus) {
2827 meMETPhiNeutralHadronsEndcapPlus_BXm1Filled =
2828 map_of_MEs[DirName +
"/" +
"METPhiNeutralHadronsEndcapPlus_BXm1Filled"];
2829 if (meMETPhiNeutralHadronsEndcapPlus_BXm1Filled &&
2830 meMETPhiNeutralHadronsEndcapPlus_BXm1Filled->getRootObject())
2831 meMETPhiNeutralHadronsEndcapPlus_BXm1Filled->Fill(
2832 atan2(py_neutralHadronsEndcapPlus, px_neutralHadronsEndcapPlus));
2834 if (pt_sum_NHF_Endcap_minus) {
2835 meMETPhiNeutralHadronsEndcapMinus_BXm1Filled =
2836 map_of_MEs[DirName +
"/" +
"METPhiNeutralHadronsEndcapMinus_BXm1Filled"];
2837 if (meMETPhiNeutralHadronsEndcapMinus_BXm1Filled &&
2838 meMETPhiNeutralHadronsEndcapMinus_BXm1Filled->getRootObject())
2839 meMETPhiNeutralHadronsEndcapMinus_BXm1Filled->Fill(
2840 atan2(py_neutralHadronsEndcapMinus, px_neutralHadronsEndcapMinus));
2842 if (pt_sum_PhF_Barrel) {
2843 meMETPhiPhotonsBarrel_BXm1Filled = map_of_MEs[DirName +
"/" +
"METPhiPhotonsBarrel_BXm1Filled"];
2844 if (meMETPhiPhotonsBarrel_BXm1Filled && meMETPhiPhotonsBarrel_BXm1Filled->getRootObject())
2845 meMETPhiPhotonsBarrel_BXm1Filled->Fill(atan2(py_PhotonsBarrel, px_PhotonsBarrel));
2847 if (pt_sum_PhF_Endcap_plus) {
2848 meMETPhiPhotonsEndcapPlus_BXm1Filled = map_of_MEs[DirName +
"/" +
"METPhiPhotonsEndcapPlus_BXm1Filled"];
2849 if (meMETPhiPhotonsEndcapPlus_BXm1Filled && meMETPhiPhotonsEndcapPlus_BXm1Filled->getRootObject())
2850 meMETPhiPhotonsEndcapPlus_BXm1Filled->Fill(atan2(py_PhotonsEndcapPlus, px_PhotonsEndcapPlus));
2852 if (pt_sum_PhF_Endcap_minus) {
2853 meMETPhiPhotonsEndcapMinus_BXm1Filled = map_of_MEs[DirName +
"/" +
"METPhiPhotonsEndcapMinus_BXm1Filled"];
2854 if (meMETPhiPhotonsEndcapMinus_BXm1Filled && meMETPhiPhotonsEndcapMinus_BXm1Filled->getRootObject())
2855 meMETPhiPhotonsEndcapMinus_BXm1Filled->Fill(atan2(py_PhotonsEndcapMinus, px_PhotonsEndcapMinus));
2857 if (pt_sum_HFH_plus) {
2858 meMETPhiHFHadronsPlus_BXm1Filled = map_of_MEs[DirName +
"/" +
"METPhiHFHadronsPlus_BXm1Filled"];
2859 if (meMETPhiHFHadronsPlus_BXm1Filled && meMETPhiHFHadronsPlus_BXm1Filled->getRootObject())
2860 meMETPhiHFHadronsPlus_BXm1Filled->Fill(atan2(py_HFHadronsPlus, px_HFHadronsPlus));
2862 if (pt_sum_HFH_minus) {
2863 meMETPhiHFHadronsMinus_BXm1Filled = map_of_MEs[DirName +
"/" +
"METPhiHFHadronsMinus_BXm1Filled"];
2864 if (meMETPhiHFHadronsMinus_BXm1Filled && meMETPhiHFHadronsMinus_BXm1Filled->getRootObject())
2865 meMETPhiHFHadronsMinus_BXm1Filled->Fill(atan2(py_HFHadronsMinus, px_HFHadronsMinus));
2867 if (pt_sum_HFE_plus) {
2868 meMETPhiHFEGammasPlus_BXm1Filled = map_of_MEs[DirName +
"/" +
"METPhiHFEGammasPlus_BXm1Filled"];
2869 if (meMETPhiHFEGammasPlus_BXm1Filled && meMETPhiHFEGammasPlus_BXm1Filled->getRootObject())
2870 meMETPhiHFEGammasPlus_BXm1Filled->Fill(atan2(py_HFEGammasPlus, px_HFEGammasPlus));
2872 if (pt_sum_HFE_minus) {
2873 meMETPhiHFEGammasMinus_BXm1Filled = map_of_MEs[DirName +
"/" +
"METPhiHFEGammasMinus_BXm1Filled"];
2874 if (meMETPhiHFEGammasMinus_BXm1Filled && meMETPhiHFEGammasMinus_BXm1Filled->getRootObject())
2875 meMETPhiHFEGammasMinus_BXm1Filled->Fill(atan2(py_HFEGammasMinus, px_HFEGammasMinus));
2937 if (techTriggerCase[2]) {
2938 meCHF_Barrel_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfChargedHadronEtFractionBarrel_BXm1Empty"];
2939 if (meCHF_Barrel_BXm1Empty && meCHF_Barrel_BXm1Empty->getRootObject())
2940 meCHF_Barrel_BXm1Empty->Fill(pt_sum_CHF_Barrel / pfmet->
sumEt());
2941 meCHF_EndcapPlus_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfChargedHadronEtFractionEndcapPlus_BXm1Empty"];
2942 if (meCHF_EndcapPlus_BXm1Empty && meCHF_EndcapPlus_BXm1Empty->getRootObject())
2943 meCHF_EndcapPlus_BXm1Empty->Fill(pt_sum_CHF_Endcap_plus / pfmet->
sumEt());
2944 meCHF_EndcapMinus_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfChargedHadronEtFractionEndcapMinus_BXm1Empty"];
2945 if (meCHF_EndcapMinus_BXm1Empty && meCHF_EndcapMinus_BXm1Empty->getRootObject())
2946 meCHF_EndcapMinus_BXm1Empty->Fill(pt_sum_CHF_Endcap_minus / pfmet->
sumEt());
2947 meNHF_Barrel_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfNeutralHadronEtFractionBarrel_BXm1Empty"];
2948 if (meNHF_Barrel_BXm1Empty && meNHF_Barrel_BXm1Empty->getRootObject())
2949 meNHF_Barrel_BXm1Empty->Fill(pt_sum_NHF_Barrel / pfmet->
sumEt());
2950 meNHF_EndcapPlus_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfNeutralHadronEtFractionEndcapPlus_BXm1Empty"];
2951 if (meNHF_EndcapPlus_BXm1Empty && meNHF_EndcapPlus_BXm1Empty->getRootObject())
2952 meNHF_EndcapPlus_BXm1Empty->Fill(pt_sum_NHF_Endcap_plus / pfmet->
sumEt());
2953 meNHF_EndcapMinus_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfNeutralHadronEtFractionEndcapMinus_BXm1Empty"];
2954 if (meNHF_EndcapMinus_BXm1Empty && meNHF_EndcapMinus_BXm1Empty->getRootObject())
2955 meNHF_EndcapMinus_BXm1Empty->Fill(pt_sum_NHF_Endcap_minus / pfmet->
sumEt());
2956 mePhF_Barrel_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfPhotonEtFractionBarrel_BXm1Empty"];
2957 if (mePhF_Barrel_BXm1Empty && mePhF_Barrel_BXm1Empty->getRootObject())
2958 mePhF_Barrel_BXm1Empty->Fill(pt_sum_PhF_Barrel / pfmet->
sumEt());
2959 mePhF_EndcapPlus_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfPhotonEtFractionEndcapPlus_BXm1Empty"];
2960 if (mePhF_EndcapPlus_BXm1Empty && mePhF_EndcapPlus_BXm1Empty->getRootObject())
2961 mePhF_EndcapPlus_BXm1Empty->Fill(pt_sum_PhF_Endcap_plus / pfmet->
sumEt());
2962 mePhF_EndcapMinus_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfPhotonEtFractionEndcapMinus_BXm1Empty"];
2963 if (mePhF_EndcapMinus_BXm1Empty && mePhF_EndcapMinus_BXm1Empty->getRootObject())
2964 mePhF_EndcapMinus_BXm1Empty->Fill(pt_sum_PhF_Endcap_minus / pfmet->
sumEt());
2965 meHFHadF_Plus_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfHFHadronEtFractionPlus_BXm1Empty"];
2966 if (meHFHadF_Plus_BXm1Empty && meHFHadF_Plus_BXm1Empty->getRootObject())
2967 meHFHadF_Plus_BXm1Empty->Fill(pt_sum_HFH_plus / pfmet->
sumEt());
2968 meHFHadF_Minus_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfHFHadronEtFractionMinus_BXm1Empty"];
2969 if (meHFHadF_Minus_BXm1Empty && meHFHadF_Minus_BXm1Empty->getRootObject())
2970 meHFHadF_Minus_BXm1Empty->Fill(pt_sum_HFH_minus / pfmet->
sumEt());
2971 meHFEMF_Plus_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfHFEMEtFractionPlus_BXm1Empty"];
2972 if (meHFEMF_Plus_BXm1Empty && meHFEMF_Plus_BXm1Empty->getRootObject())
2973 meHFEMF_Plus_BXm1Empty->Fill(pt_sum_HFE_plus / pfmet->
sumEt());
2974 meHFEMF_Minus_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfHFEMEtFractionMinus_BXm1Empty"];
2975 if (meHFEMF_Minus_BXm1Empty && meHFEMF_Minus_BXm1Empty->getRootObject())
2976 meHFEMF_Minus_BXm1Empty->Fill(pt_sum_HFE_minus / pfmet->
sumEt());
2977 mePhotonEtFraction_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfPhotonEtFraction_BXm1Empty"];
2978 if (mePhotonEtFraction_BXm1Empty && mePhotonEtFraction_BXm1Empty->getRootObject())
2980 meNeutralHadronEtFraction_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfNeutralHadronEtFraction_BXm1Empty"];
2981 if (meNeutralHadronEtFraction_BXm1Empty && meNeutralHadronEtFraction_BXm1Empty->getRootObject())
2983 meChargedHadronEtFraction_BXm1Empty = map_of_MEs[DirName +
"/" +
"PfChargedHadronEtFraction_BXm1Empty"];
2984 if (meChargedHadronEtFraction_BXm1Empty && meChargedHadronEtFraction_BXm1Empty->getRootObject())
2986 meMET_BXm1Empty = map_of_MEs[DirName +
"/" +
"MET_BXm1Empty"];
2987 if (meMET_BXm1Empty && meMET_BXm1Empty->getRootObject())
2988 meMET_BXm1Empty->Fill(pfmet->
pt());
2989 meSumET_BXm1Empty = map_of_MEs[DirName +
"/" +
"SumET_BXm1Empty"];
2990 if (meSumET_BXm1Empty && meSumET_BXm1Empty->getRootObject())
2991 meSumET_BXm1Empty->Fill(pfmet->
sumEt());
2992 if (pt_sum_CHF_Barrel) {
2993 meMETPhiChargedHadronsBarrel_BXm1Empty = map_of_MEs[DirName +
"/" +
"METPhiChargedHadronsBarrel_BXm1Empty"];
2994 if (meMETPhiChargedHadronsBarrel_BXm1Empty && meMETPhiChargedHadronsBarrel_BXm1Empty->getRootObject())
2995 meMETPhiChargedHadronsBarrel_BXm1Empty->Fill(atan2(py_chargedHadronsBarrel, px_chargedHadronsBarrel));
2997 if (pt_sum_CHF_Endcap_plus) {
2998 meMETPhiChargedHadronsEndcapPlus_BXm1Empty =
2999 map_of_MEs[DirName +
"/" +
"METPhiChargedHadronsEndcapPlus_BXm1Empty"];
3000 if (meMETPhiChargedHadronsEndcapPlus_BXm1Empty &&
3001 meMETPhiChargedHadronsEndcapPlus_BXm1Empty->getRootObject())
3002 meMETPhiChargedHadronsEndcapPlus_BXm1Empty->Fill(
3003 atan2(py_chargedHadronsEndcapPlus, px_chargedHadronsEndcapPlus));
3005 if (pt_sum_CHF_Endcap_minus) {
3006 meMETPhiChargedHadronsEndcapMinus_BXm1Empty =
3007 map_of_MEs[DirName +
"/" +
"METPhiChargedHadronsEndcapMinus_BXm1Empty"];
3008 if (meMETPhiChargedHadronsEndcapMinus_BXm1Empty &&
3009 meMETPhiChargedHadronsEndcapMinus_BXm1Empty->getRootObject())
3010 meMETPhiChargedHadronsEndcapMinus_BXm1Empty->Fill(
3011 atan2(py_chargedHadronsEndcapMinus, px_chargedHadronsEndcapMinus));
3013 if (pt_sum_NHF_Barrel) {
3014 meMETPhiNeutralHadronsBarrel_BXm1Empty = map_of_MEs[DirName +
"/" +
"METPhiNeutralHadronsBarrel_BXm1Empty"];
3015 if (meMETPhiNeutralHadronsBarrel_BXm1Empty && meMETPhiNeutralHadronsBarrel_BXm1Empty->getRootObject())
3016 meMETPhiNeutralHadronsBarrel_BXm1Empty->Fill(atan2(py_neutralHadronsBarrel, px_neutralHadronsBarrel));
3018 if (pt_sum_NHF_Endcap_plus) {
3019 meMETPhiNeutralHadronsEndcapPlus_BXm1Empty =
3020 map_of_MEs[DirName +
"/" +
"METPhiNeutralHadronsEndcapPlus_BXm1Empty"];
3021 if (meMETPhiNeutralHadronsEndcapPlus_BXm1Empty &&
3022 meMETPhiNeutralHadronsEndcapPlus_BXm1Empty->getRootObject())
3023 meMETPhiNeutralHadronsEndcapPlus_BXm1Empty->Fill(
3024 atan2(py_neutralHadronsEndcapPlus, px_neutralHadronsEndcapPlus));
3026 if (pt_sum_NHF_Endcap_minus) {
3027 meMETPhiNeutralHadronsEndcapMinus_BXm1Empty =
3028 map_of_MEs[DirName +
"/" +
"METPhiNeutralHadronsEndcapMinus_BXm1Empty"];
3029 if (meMETPhiNeutralHadronsEndcapMinus_BXm1Empty &&
3030 meMETPhiNeutralHadronsEndcapMinus_BXm1Empty->getRootObject())
3031 meMETPhiNeutralHadronsEndcapMinus_BXm1Empty->Fill(
3032 atan2(py_neutralHadronsEndcapMinus, px_neutralHadronsEndcapMinus));
3034 if (pt_sum_PhF_Barrel) {
3035 meMETPhiPhotonsBarrel_BXm1Empty = map_of_MEs[DirName +
"/" +
"METPhiPhotonsBarrel_BXm1Empty"];
3036 if (meMETPhiPhotonsBarrel_BXm1Empty && meMETPhiPhotonsBarrel_BXm1Empty->getRootObject())
3037 meMETPhiPhotonsBarrel_BXm1Empty->Fill(atan2(py_PhotonsBarrel, px_PhotonsBarrel));
3039 if (pt_sum_PhF_Endcap_plus) {
3040 meMETPhiPhotonsEndcapPlus_BXm1Empty = map_of_MEs[DirName +
"/" +
"METPhiPhotonsEndcapPlus_BXm1Empty"];
3041 if (meMETPhiPhotonsEndcapPlus_BXm1Empty && meMETPhiPhotonsEndcapPlus_BXm1Empty->getRootObject())
3042 meMETPhiPhotonsEndcapPlus_BXm1Empty->Fill(atan2(py_PhotonsEndcapPlus, px_PhotonsEndcapPlus));
3044 if (pt_sum_PhF_Endcap_minus) {
3045 meMETPhiPhotonsEndcapMinus_BXm1Empty = map_of_MEs[DirName +
"/" +
"METPhiPhotonsEndcapMinus_BXm1Empty"];
3046 if (meMETPhiPhotonsEndcapMinus_BXm1Empty && meMETPhiPhotonsEndcapMinus_BXm1Empty->getRootObject())
3047 meMETPhiPhotonsEndcapMinus_BXm1Empty->Fill(atan2(py_PhotonsEndcapMinus, px_PhotonsEndcapMinus));
3049 if (pt_sum_HFH_plus) {
3050 meMETPhiHFHadronsPlus_BXm1Empty = map_of_MEs[DirName +
"/" +
"METPhiHFHadronsPlus_BXm1Empty"];
3051 if (meMETPhiHFHadronsPlus_BXm1Empty && meMETPhiHFHadronsPlus_BXm1Empty->getRootObject())
3052 meMETPhiHFHadronsPlus_BXm1Empty->Fill(atan2(py_HFHadronsPlus, px_HFHadronsPlus));
3054 if (pt_sum_HFH_minus) {
3055 meMETPhiHFHadronsMinus_BXm1Empty = map_of_MEs[DirName +
"/" +
"METPhiHFHadronsMinus_BXm1Empty"];
3056 if (meMETPhiHFHadronsMinus_BXm1Empty && meMETPhiHFHadronsMinus_BXm1Empty->getRootObject())
3057 meMETPhiHFHadronsMinus_BXm1Empty->Fill(atan2(py_HFHadronsMinus, px_HFHadronsMinus));
3059 if (pt_sum_HFE_plus) {
3060 meMETPhiHFEGammasPlus_BXm1Empty = map_of_MEs[DirName +
"/" +
"METPhiHFEGammasPlus_BXm1Empty"];
3061 if (meMETPhiHFEGammasPlus_BXm1Empty && meMETPhiHFEGammasPlus_BXm1Empty->getRootObject())
3062 meMETPhiHFEGammasPlus_BXm1Empty->Fill(atan2(py_HFEGammasPlus, px_HFEGammasPlus));
3064 if (pt_sum_HFE_minus) {
3065 meMETPhiHFEGammasMinus_BXm1Empty = map_of_MEs[DirName +
"/" +
"METPhiHFEGammasMinus_BXm1Empty"];
3066 if (meMETPhiHFEGammasMinus_BXm1Empty && meMETPhiHFEGammasMinus_BXm1Empty->getRootObject())
3067 meMETPhiHFEGammasMinus_BXm1Empty->Fill(atan2(py_HFEGammasMinus, px_HFEGammasMinus));
3076 double pfPhotonEt = pfmet->
photonEt();
3082 double pfMuonEt = pfmet->
muonEt();
3086 double pfHFEMEt = pfmet->
HFEMEt();
3087 mePhotonEtFraction = map_of_MEs[DirName +
"/PfPhotonEtFraction"];
3088 mePhotonEt = map_of_MEs[DirName +
"/PfPhotonEt"];
3089 meNeutralHadronEtFraction = map_of_MEs[DirName +
"/PfNeutralHadronEtFraction"];
3090 meNeutralHadronEt = map_of_MEs[DirName +
"/PfNeutralHadronEt"];
3091 meElectronEt = map_of_MEs[DirName +
"/PfElectronEt"];
3092 meChargedHadronEtFraction = map_of_MEs[DirName +
"/PfChargedHadronEtFraction"];
3093 meChargedHadronEt = map_of_MEs[DirName +
"/PfChargedHadronEt"];
3094 meMuonEt = map_of_MEs[DirName +
"/PfMuonEt"];
3095 meHFHadronEtFraction = map_of_MEs[DirName +
"/PfHFHadronEtFraction"];
3096 meHFHadronEt = map_of_MEs[DirName +
"/PfHFHadronEt"];
3097 meHFEMEtFraction = map_of_MEs[DirName +
"/PfHFEMEtFraction"];
3098 meHFEMEt = map_of_MEs[DirName +
"/PfHFEMEt"];
3100 if (mePhotonEtFraction && mePhotonEtFraction->getRootObject())
3101 mePhotonEtFraction->Fill(pfPhotonEtFraction);
3102 if (mePhotonEt && mePhotonEt->getRootObject())
3103 mePhotonEt->Fill(pfPhotonEt);
3104 if (meNeutralHadronEtFraction && meNeutralHadronEtFraction->getRootObject())
3105 meNeutralHadronEtFraction->Fill(pfNeutralHadronEtFraction);
3106 if (meNeutralHadronEt && meNeutralHadronEt->getRootObject())
3107 meNeutralHadronEt->Fill(pfNeutralHadronEt);
3108 if (meElectronEt && meElectronEt->getRootObject())
3109 meElectronEt->Fill(pfElectronEt);
3110 if (meChargedHadronEtFraction && meChargedHadronEtFraction->getRootObject())
3111 meChargedHadronEtFraction->Fill(pfChargedHadronEtFraction);
3112 if (meChargedHadronEt && meChargedHadronEt->getRootObject())
3113 meChargedHadronEt->Fill(pfChargedHadronEt);
3114 if (meMuonEt && meMuonEt->getRootObject())
3115 meMuonEt->Fill(pfMuonEt);
3116 if (meHFHadronEtFraction && meHFHadronEtFraction->getRootObject())
3117 meHFHadronEtFraction->Fill(pfHFHadronEtFraction);
3118 if (meHFHadronEt && meHFHadronEt->getRootObject())
3119 meHFHadronEt->Fill(pfHFHadronEt);
3120 if (meHFEMEtFraction && meHFEMEtFraction->getRootObject())
3121 meHFEMEtFraction->Fill(pfHFEMEtFraction);
3122 if (meHFEMEt && meHFEMEt->getRootObject())
3123 meHFEMEt->Fill(pfHFEMEt);
3127 mePhotonEtFraction_profile = map_of_MEs[DirName +
"/PfPhotonEtFraction_profile"];
3128 mePhotonEt_profile = map_of_MEs[DirName +
"/PfPhotonEt_profile"];
3129 meNeutralHadronEtFraction_profile = map_of_MEs[DirName +
"/PfNeutralHadronEtFraction_profile"];
3130 meNeutralHadronEt_profile = map_of_MEs[DirName +
"/PfNeutralHadronEt_profile"];
3131 meChargedHadronEtFraction_profile = map_of_MEs[DirName +
"/PfChargedHadronEtFraction_profile"];
3132 meChargedHadronEt_profile = map_of_MEs[DirName +
"/PfChargedHadronEt_profile"];
3133 meHFHadronEtFraction_profile = map_of_MEs[DirName +
"/PfHFHadronEtFraction_profile"];
3134 meHFHadronEt_profile = map_of_MEs[DirName +
"/PfHFHadronEt_profile"];
3135 meHFEMEtFraction_profile = map_of_MEs[DirName +
"/PfHFEMEtFraction_profile"];
3136 meHFEMEt_profile = map_of_MEs[DirName +
"/PfHFEMEt_profile"];
3138 if (mePhotonEtFraction_profile && mePhotonEtFraction_profile->getRootObject())
3139 mePhotonEtFraction_profile->Fill(numPV_, pfPhotonEtFraction);
3140 if (mePhotonEt_profile && mePhotonEt_profile->getRootObject())
3141 mePhotonEt_profile->Fill(numPV_, pfPhotonEt);
3142 if (meNeutralHadronEtFraction_profile && meNeutralHadronEtFraction_profile->getRootObject())
3143 meNeutralHadronEtFraction_profile->Fill(numPV_, pfNeutralHadronEtFraction);
3144 if (meNeutralHadronEt_profile && meNeutralHadronEt_profile->getRootObject())
3145 meNeutralHadronEt_profile->Fill(numPV_, pfNeutralHadronEt);
3146 if (meChargedHadronEtFraction_profile && meChargedHadronEtFraction_profile->getRootObject())
3147 meChargedHadronEtFraction_profile->Fill(numPV_, pfChargedHadronEtFraction);
3148 if (meChargedHadronEt_profile && meChargedHadronEt_profile->getRootObject())
3149 meChargedHadronEt_profile->Fill(numPV_, pfChargedHadronEt);
3150 if (meHFHadronEtFraction_profile && meHFHadronEtFraction_profile->getRootObject())
3151 meHFHadronEtFraction_profile->Fill(numPV_, pfHFHadronEtFraction);
3152 if (meHFHadronEt_profile && meHFHadronEt_profile->getRootObject())
3153 meHFHadronEt_profile->Fill(numPV_, pfHFHadronEt);
3154 if (meHFEMEtFraction_profile && meHFEMEtFraction_profile->getRootObject())
3155 meHFEMEtFraction_profile->Fill(numPV_, pfHFEMEtFraction);
3156 if (meHFEMEt_profile && meHFEMEt_profile->getRootObject())
3157 meHFEMEt_profile->Fill(numPV_, pfHFEMEt);
3160 if (isMiniAODMet_) {
3161 mePhotonEtFraction = map_of_MEs[DirName +
"/PfPhotonEtFraction"];
3162 meNeutralHadronEtFraction = map_of_MEs[DirName +
"/PfNeutralHadronEtFraction"];
3163 meChargedHadronEtFraction = map_of_MEs[DirName +
"/PfChargedHadronEtFraction"];
3164 meHFHadronEtFraction = map_of_MEs[DirName +
"/PfHFHadronEtFraction"];
3165 meHFEMEtFraction = map_of_MEs[DirName +
"/PfHFEMEtFraction"];
3167 if (mePhotonEtFraction && mePhotonEtFraction->getRootObject())
3169 if (meNeutralHadronEtFraction && meNeutralHadronEtFraction->getRootObject())
3171 if (meChargedHadronEtFraction && meChargedHadronEtFraction->getRootObject())
3173 if (meHFHadronEtFraction && meHFHadronEtFraction->getRootObject())
3175 if (meHFEMEtFraction && meHFEMEtFraction->getRootObject())
3179 mePhotonEtFraction_profile = map_of_MEs[DirName +
"/PfPhotonEtFraction_profile"];
3180 meNeutralHadronEtFraction_profile = map_of_MEs[DirName +
"/PfNeutralHadronEtFraction_profile"];
3181 meChargedHadronEtFraction_profile = map_of_MEs[DirName +
"/PfChargedHadronEtFraction_profile"];
3182 meHFHadronEtFraction_profile = map_of_MEs[DirName +
"/PfHFHadronEtFraction_profile"];
3183 meHFEMEtFraction_profile = map_of_MEs[DirName +
"/PfHFEMEtFraction_profile"];
3185 if (mePhotonEtFraction_profile && mePhotonEtFraction_profile->getRootObject())
3187 if (meNeutralHadronEtFraction_profile && meNeutralHadronEtFraction_profile->getRootObject())
3189 if (meChargedHadronEtFraction_profile && meChargedHadronEtFraction_profile->getRootObject())
3191 if (meHFHadronEtFraction_profile && meHFHadronEtFraction_profile->getRootObject())
3192 meHFHadronEtFraction_profile->Fill(numPV_, patmet->
Type6EtFraction());
3193 if (meHFEMEtFraction_profile && meHFEMEtFraction_profile->getRootObject())
3196 mePhotonEt = map_of_MEs[DirName +
"/PfPhotonEt"];
3197 meNeutralHadronEt = map_of_MEs[DirName +
"/PfNeutralHadronEt"];
3198 meChargedHadronEt = map_of_MEs[DirName +
"/PfChargedHadronEt"];
3199 meHFHadronEt = map_of_MEs[DirName +
"/PfHFHadronEt"];
3200 meHFEMEt = map_of_MEs[DirName +
"/PfHFEMEt"];
3201 meMuonEt = map_of_MEs[DirName +
"/PfMuonEt"];
3202 meElectronEt = map_of_MEs[DirName +
"/PfElectronEt"];
3204 if (mePhotonEt && mePhotonEt->getRootObject())
3206 if (meNeutralHadronEt && meNeutralHadronEt->getRootObject())
3208 if (meChargedHadronEt && meChargedHadronEt->getRootObject())
3210 if (meHFHadronEt && meHFHadronEt->getRootObject())
3212 if (meHFEMEt && meHFEMEt->getRootObject())
3214 if (meMuonEt && meMuonEt->getRootObject())
3216 if (meElectronEt && meElectronEt->getRootObject())
3220 mePhotonEt_profile = map_of_MEs[DirName +
"/PfPhotonEt_profile"];
3221 meNeutralHadronEt_profile = map_of_MEs[DirName +
"/PfNeutralHadronEt_profile"];
3222 meChargedHadronEt_profile = map_of_MEs[DirName +
"/PfChargedHadronEt_profile"];
3223 meHFHadronEt_profile = map_of_MEs[DirName +
"/PfHFHadronEt_profile"];
3224 meHFEMEt_profile = map_of_MEs[DirName +
"/PfHFEMEt_profile"];
3226 if (mePhotonEt_profile && mePhotonEt_profile->getRootObject())
3228 if (meNeutralHadronEt_profile && meNeutralHadronEt_profile->getRootObject())
3230 if (meChargedHadronEt_profile && meChargedHadronEt_profile->getRootObject())
3232 if (meHFHadronEt_profile && meHFHadronEt_profile->getRootObject())
3234 if (meHFEMEt_profile && meHFEMEt_profile->getRootObject())
3240 if (fill_met_high_level_histo) {
3241 hMExLS = map_of_MEs[DirName +
"/" +
"MExLS"];
3242 if (hMExLS && hMExLS->getRootObject())
3243 hMExLS->Fill(MEx, myLuminosityBlock);
3244 hMEyLS = map_of_MEs[DirName +
"/" +
"MEyLS"];
3245 if (hMEyLS && hMEyLS->getRootObject())
3246 hMEyLS->Fill(MEy, myLuminosityBlock);
Analysis-level MET class.
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX)
T getParameter(std::string const &) const
double eta() const final
momentum pseudorapidity
double HFEMEtFraction() const
const TechnicalTriggerWord & technicalTriggerWord(int bxInEventValue) const
double ChargedHadEtFraction() const
double MuonEtFraction() const
bool getByToken(EDGetToken token, Handle< PROD > &result) const
std::vector< T >::const_iterator search(const cond::Time_t &val, const std::vector< T > &container)
void setCurrentFolder(std::string const &fullpath)
std::vector< ParameterSet > VParameterSet
double neutralHadronEtFraction() const
double px() const final
x coordinate of momentum vector
bool accept() const
Has at least one path accepted the event?
void setLumiFlag()
this ME is meant to be stored for each luminosity section
edm::LuminosityBlockNumber_t luminosityBlock() const
double pt() const final
transverse momentum
std::vector< Vertex > VertexCollection
collection of Vertex objects
double NeutralEMFraction() const
double correction(const LorentzVector &fJet) const
get correction using Jet information only
Provides a code based selection for trigger and DCS information in order to have no failing filters i...
double ChargedEMEtFraction() const
double Type6EtFraction() const
double Type7EtFraction() const
double NeutralHadEtFraction() const
double HFHadronEtFraction() const
double HFHadronEt() const
double et() const final
transverse energy
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
double photonEtFraction() const
MonitorElement * bookProfile(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, char const *option="s")
std::vector< bool > TechnicalTriggerWord
technical trigger bits (64 bits)
numberOfValidPixelHits
at most 999 lost hits
virtual void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
PF Jet selector for pat::Jets.
double electronEt() const
Jet selector for pat::Jets and for CaloJets.
double etFractionHadronic() const
T const * product() const
virtual void setBinContent(int binx, double content)
set content of bin (1-D)
static std::string const triggerResults("TriggerResults")
double py() const final
y coordinate of momentum vector
double chargedHadronEtFraction() const
Particle reconstructed by the particle flow algorithm.
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
static int position[264][3]
double neutralHadronEt() const
virtual ParticleType particleId() const
virtual TObject const * getRootObject() const
double phi() const final
momentum azimuthal angle
T const * product() const
double emEtFraction() const
double chargedHadronEt() const
virtual void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
math::PtEtaPhiMLorentzVector PolarLorentzVector
Lorentz vector.