67 isMiniAODMet_ = (
std::string(
"miniaod") ==MetType_);
77 jetID_ValueMapToken_= consumes< edm::ValueMap<reco::JetID> >(inputJetIDValueMap);
91 ptThreshold_ =
parameters.getParameter<
double>(
"ptThreshold");
95 pfMetToken_= consumes<reco::PFMETCollection>(
edm::InputTag(metCollectionLabel_));
98 caloMetToken_= consumes<reco::CaloMETCollection>(
edm::InputTag(metCollectionLabel_));
101 patMetToken_= consumes<pat::METCollection>(
edm::InputTag(metCollectionLabel_));
107 fill_met_high_level_histo =
parameters.getParameter<
bool>(
"fillMetHighLevel");
108 fillCandidateMap_histos =
parameters.getParameter<
bool>(
"fillCandidateMaps");
110 hTriggerLabelsIsSet_ =
false;
114 diagnosticsParameters_ = pSet.
getParameter<std::vector<edm::ParameterSet> >(
"METDiagonisticsParameters");
121 bypassAllPVChecks_ = cleaningParameters_.getParameter<
bool>(
"bypassAllPVChecks");
122 bypassAllDCSChecks_ = cleaningParameters_.getParameter<
bool>(
"bypassAllDCSChecks");
123 runcosmics_ =
parameters.getUntrackedParameter<
bool>(
"runcosmics");
124 onlyCleaned_ =
parameters.getUntrackedParameter<
bool>(
"onlyCleaned");
125 vertexTag_ = cleaningParameters_.getParameter<
edm::InputTag>(
"vertexCollection");
126 vertexToken_ = consumes<std::vector<reco::Vertex> >(
edm::InputTag(vertexTag_));
129 gtTag_ = cleaningParameters_.getParameter<
edm::InputTag>(
"gtLabel");
130 gtToken_= consumes<L1GlobalTriggerReadoutRecord>(
edm::InputTag(gtTag_));
134 if (isCaloMet_) caloJetsToken_ = consumes<reco::CaloJetCollection>(jetCollectionLabel_);
136 if (isPFMet_) pfJetsToken_ = consumes<reco::PFJetCollection>(jetCollectionLabel_);
137 if (isMiniAODMet_) patJetsToken_ = consumes<pat::JetCollection>(jetCollectionLabel_);
143 hbheNoiseFilterResultToken_=consumes<bool>(hbheNoiseFilterResultTag_);
145 hbheIsoNoiseFilterResultToken_ =consumes<bool> (hbheNoiseFilterResultTag_);
147 CSCHaloResultToken_=consumes<bool>(CSCHaloResultTag_);
149 CSCHalo2015ResultToken_=consumes<bool> (CSCHalo2015ResultTag_);
151 EcalDeadCellTriggerToken_=consumes<bool>(EcalDeadCellTriggerTag_);
153 EcalDeadCellBoundaryToken_=consumes<bool>(EcalDeadCellBoundaryTag_);
155 eeBadScFilterToken_=consumes<bool>(eeBadScFilterTag_);
157 HcalStripHaloToken_=consumes<bool>(HcalStripHaloTag_);
160 METFilterMiniAODToken_=consumes<edm::TriggerResults>(METFilterMiniAODLabel_);
163 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(); it!= triggerSelectedSubFolders_.end(); it++) {
173 triggerFolderExpr_.push_back(it->getParameter<std::vector<std::string> >(
"hltPaths"));
174 triggerFolderLabels_.push_back(it->getParameter<
std::string>(
"label"));
179 verbose_ =
parameters.getParameter<
int>(
"verbose");
183 nCh = std::vector<unsigned int>(10,
static_cast<unsigned int>(0));
184 nEv = std::vector<unsigned int>(2,
static_cast<unsigned int>(0));
190 for (std::vector<GenericTriggerEventFlag *>::const_iterator it = triggerFolderEventFlag_.begin(); it!= triggerFolderEventFlag_.end(); it++) {
203 if(!folderNames_.empty()){
204 folderNames_.clear();
207 folderNames_.push_back(
"Uncleaned");
210 folderNames_.push_back(
"Uncleaned");
212 folderNames_.push_back(
"Cleaned");
213 folderNames_.push_back(
"DiJet");
215 folderNames_.push_back(
"ZJets");
218 for (std::vector<std::string>::const_iterator ic = folderNames_.begin();
219 ic != folderNames_.end(); ic++){
220 bookMESet(DirName+
"/"+*ic, ibooker,map_dijet_MEs);
228 bool bLumiSecPlot=fill_met_high_level_histo;
230 bool fillPFCandidatePlots=
false;
231 bool fillZPlots=
false;
233 if (DirName.find(
"Cleaned")!=std::string::npos) {
234 fillPFCandidatePlots=
true;
235 bookMonitorElement(DirName,ibooker,map_of_MEs,bLumiSecPlot,fillPFCandidatePlots,fillZPlots);
242 }
else if (DirName.find(
"ZJets")!=std::string::npos) {
243 fillPFCandidatePlots=
false;
245 bookMonitorElement(DirName,ibooker,map_of_MEs,bLumiSecPlot,fillPFCandidatePlots,fillZPlots);
247 bookMonitorElement(DirName,ibooker,map_of_MEs,bLumiSecPlot,fillPFCandidatePlots,fillZPlots);
254 if (verbose_)
std::cout <<
"bookMonitorElement " << DirName << std::endl;
259 meZJets_u_par = ibooker.
book1D(
"u_parallel_Z_inc",
"u_parallel_Z_inc", 50, -1000., 75);
261 meZJets_u_par = ibooker.
book1D(
"u_parallel_Z_inc",
"u_parallel_Z_inc", 50, -800., 75);
263 meZJets_u_par_ZPt_0_15 = ibooker.
book1D(
"u_parallel_ZPt_0_15",
"u_parallel_ZPt_0_15", 50, -100, 75);
264 meZJets_u_par_ZPt_15_30 = ibooker.
book1D(
"u_parallel_ZPt_15_30",
"u_parallel_ZPt_15_30", 50, -100, 50);
265 meZJets_u_par_ZPt_30_55 = ibooker.
book1D(
"u_parallel_ZPt_30_55",
"u_parallel_ZPt_30_55", 50, -175, 50);
266 meZJets_u_par_ZPt_55_75 = ibooker.
book1D(
"u_parallel_ZPt_55_75",
"u_parallel_ZPt_55_75", 50, -175, 0);
267 meZJets_u_par_ZPt_75_150 = ibooker.
book1D(
"u_parallel_ZPt_75_150",
"u_parallel_ZPt_75_150", 50, -300, 0);
269 meZJets_u_par_ZPt_150_290 = ibooker.
book1D(
"u_parallel_ZPt_150_290",
"u_parallel_ZPt_150_290", 50, -750, -100);
271 meZJets_u_par_ZPt_150_290 = ibooker.
book1D(
"u_parallel_ZPt_150_290",
"u_parallel_ZPt_150_290", 50, -450, -50);
274 meZJets_u_par_ZPt_290 = ibooker.
book1D(
"u_parallel_ZPt_290",
"u_parallel_ZPt_290", 50, -1000., -350.);
276 meZJets_u_par_ZPt_290 = ibooker.
book1D(
"u_parallel_ZPt_290",
"u_parallel_ZPt_290", 50, -750., -150.);
278 meZJets_u_perp = ibooker.
book1D(
"u_perp_Z_inc",
"u_perp_Z_inc", 50, -85., 85.);
279 meZJets_u_perp_ZPt_0_15 = ibooker.
book1D(
"u_perp_ZPt_0_15",
"u_perp_ZPt_0_15", 50, -85., 85.);
280 meZJets_u_perp_ZPt_15_30 = ibooker.
book1D(
"u_perp_ZPt_15_30",
"u_perp_ZPt_15_30", 50, -85., 85.);
281 meZJets_u_perp_ZPt_30_55 = ibooker.
book1D(
"u_perp_ZPt_30_55",
"u_perp_ZPt_30_55", 50, -85., 85.);
282 meZJets_u_perp_ZPt_55_75 = ibooker.
book1D(
"u_perp_ZPt_55_75",
"u_perp_ZPt_55_75", 50, -85., 85.);
283 meZJets_u_perp_ZPt_75_150 = ibooker.
book1D(
"u_perp_ZPt_75_150",
"u_perp_ZPt_75_150", 50, -85., 85.);
284 meZJets_u_perp_ZPt_150_290 = ibooker.
book1D(
"u_perp_ZPt_150_290",
"u_perp_ZPt_150_290", 50, -85., 85.);
285 meZJets_u_perp_ZPt_290 = ibooker.
book1D(
"u_perp_ZPt_290",
"u_perp_ZPt_290", 50, -85., 85.);
287 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_parallel_Z_inc",meZJets_u_par));
288 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_parallel_ZPt_0_15",meZJets_u_par_ZPt_0_15));
289 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_parallel_ZPt_15_30",meZJets_u_par_ZPt_15_30));
290 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_parallel_ZPt_30_55",meZJets_u_par_ZPt_30_55));
291 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_parallel_ZPt_55_75",meZJets_u_par_ZPt_55_75));
292 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_parallel_ZPt_75_150",meZJets_u_par_ZPt_75_150));
293 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_parallel_ZPt_150_290",meZJets_u_par_ZPt_150_290));
294 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_parallel_ZPt_290",meZJets_u_par_ZPt_290));
296 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_perp_Z_inc",meZJets_u_perp));
297 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_perp_ZPt_0_15",meZJets_u_perp_ZPt_0_15));
298 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_perp_ZPt_15_30",meZJets_u_perp_ZPt_15_30));
299 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_perp_ZPt_30_55",meZJets_u_perp_ZPt_30_55));
300 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_perp_ZPt_55_75",meZJets_u_perp_ZPt_55_75));
301 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_perp_ZPt_75_150",meZJets_u_perp_ZPt_75_150));
302 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_perp_ZPt_150_290",meZJets_u_perp_ZPt_150_290));
303 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"u_perp_ZPt_290",meZJets_u_perp_ZPt_290));
307 hTrigger = ibooker.
book1D(
"triggerResults",
"triggerResults", 500, 0, 500);
308 for (
unsigned int i = 0;
i<allTriggerNames_.size();
i++){
309 if(
i<(
unsigned int)hTrigger->getNbinsX()){
310 if (!hTriggerLabelsIsSet_) {
315 if (!hTriggerLabelsIsSet_)
for (
int i = allTriggerNames_.size();
i<hTrigger->getNbinsX();
i++){
316 hTrigger->setBinLabel(
i+1,
"");
318 hTriggerLabelsIsSet_ =
true;
320 hMEx = ibooker.
book1D(
"MEx",
"MEx", 200, -500, 500);
321 hMEy = ibooker.
book1D(
"MEy",
"MEy", 200, -500, 500);
322 hMET = ibooker.
book1D(
"MET",
"MET", 200, 0, 1000);
323 hMET_2 = ibooker.
book1D(
"MET_2",
"MET Range 2",200, 0, 2000);
324 hSumET = ibooker.
book1D(
"SumET",
"SumET", 400, 0, 4000);
325 hMETSig = ibooker.
book1D(
"METSig",
"METSig", 51, 0, 51);
327 hMET_logx = ibooker.
book1D(
"MET_logx",
"MET_logx", 40, -1, 9);
328 hSumET_logx = ibooker.
book1D(
"SumET_logx",
"SumET_logx", 40, -1, 9);
331 hMEy ->setAxisTitle(
"MEy [GeV]", 1);
332 hMET ->setAxisTitle(
"MET [GeV]", 1);
333 hMET_2 ->setAxisTitle(
"MET [GeV]", 1);
334 hSumET ->setAxisTitle(
"SumET [GeV]", 1);
335 hMETSig ->setAxisTitle(
"METSig", 1);
336 hMETPhi ->setAxisTitle(
"METPhi [rad]", 1);
337 hMET_logx ->setAxisTitle(
"log(MET) [GeV]", 1);
338 hSumET_logx->setAxisTitle(
"log(SumET) [GeV]", 1);
340 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"triggerResults",hTrigger));
341 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MEx",hMEx));
342 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MEy",hMEy));
343 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET",hMET));
344 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_2",hMET_2));
345 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"SumET",hSumET));
346 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METSig",hMETSig));
347 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhi",hMETPhi));
348 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_logx",hMET_logx));
349 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"SumET_logx",hSumET_logx));
351 hMET_HBHENoiseFilter = ibooker.
book1D(
"MET_HBHENoiseFilter",
"MET_HBHENoiseFiltered", 200, 0, 1000);
352 hMET_CSCTightHaloFilter = ibooker.
book1D(
"MET_CSCTightHaloFilter",
"MET_CSCTightHaloFiltered", 200, 0, 1000);
353 hMET_eeBadScFilter = ibooker.
book1D(
"MET_eeBadScFilter",
"MET_eeBadScFiltered", 200, 0, 1000);
354 hMET_HBHEIsoNoiseFilter = ibooker.
book1D(
"MET_HBHEIsoNoiseFilter",
"MET_HBHEIsoNoiseFiltered", 200, 0, 1000);
355 hMET_CSCTightHalo2015Filter = ibooker.
book1D(
"MET_CSCTightHalo2015Filter",
"MET_CSCTightHalo2015Filtered", 200, 0, 1000);
356 hMET_EcalDeadCellTriggerFilter = ibooker.
book1D(
"MET_EcalDeadCellTriggerFilter",
"MET_EcalDeadCellTriggerFiltered", 200, 0, 1000);
357 hMET_EcalDeadCellBoundaryFilter = ibooker.
book1D(
"MET_EcalDeadCellBoundaryFilter",
"MET_EcalDeadCellBoundaryFiltered", 200, 0, 1000);
358 hMET_HcalStripHaloFilter = ibooker.
book1D(
"MET_HcalStripHaloFilter",
"MET_HcalStripHaloFiltered", 200, 0, 1000);
360 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_HBHENoiseFilter",hMET_HBHENoiseFilter));
361 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_CSCTightHaloFilter",hMET_CSCTightHaloFilter));
362 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_eeBadScFilter",hMET_eeBadScFilter));
363 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_HBHEIsoNoiseFilter",hMET_HBHEIsoNoiseFilter));
364 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_CSCTightHalo2015Filter",hMET_CSCTightHalo2015Filter));
365 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_EcalDeadCellTriggerFilter",hMET_EcalDeadCellTriggerFilter));
366 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_EcalDeadCellBoundaryFilter",hMET_EcalDeadCellBoundaryFilter));
367 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_HcalStripHaloFilter",hMET_HcalStripHaloFilter));
372 meMEx_profile = ibooker.
bookProfile(
"MEx_profile",
"met.px()", nbinsPV_, nPVMin_, nPVMax_, 200, -500, 500);
373 meMEy_profile = ibooker.
bookProfile(
"MEy_profile",
"met.py()", nbinsPV_, nPVMin_, nPVMax_, 200, -500, 500);
374 meMET_profile = ibooker.
bookProfile(
"MET_profile",
"met.pt()", nbinsPV_, nPVMin_, nPVMax_, 200, 0, 1000);
375 meSumET_profile = ibooker.
bookProfile(
"SumET_profile",
"met.sumEt()", nbinsPV_, nPVMin_, nPVMax_, 400, 0, 4000);
379 meMEy_profile ->setAxisTitle(
"nvtx", 1);
380 meMET_profile ->setAxisTitle(
"nvtx", 1);
381 meSumET_profile->setAxisTitle(
"nvtx", 1);
383 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MEx_profile",meMEx_profile));
384 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MEy_profile",meMEy_profile));
385 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_profile",meMET_profile));
386 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"SumET_profile",meSumET_profile));
390 hCaloHadEtInHB = ibooker.
book1D(
"CaloHadEtInHB",
"CaloHadEtInHB",50,0,2000);
392 hCaloHadEtInHO = ibooker.
book1D(
"CaloHadEtInHO",
"CaloHadEtInHO",25,0,500);
394 hCaloHadEtInHE = ibooker.
book1D(
"CaloHadEtInHE",
"CaloHadEtInHE",50,0,2000);
396 hCaloHadEtInHF = ibooker.
book1D(
"CaloHadEtInHF",
"CaloHadEtInHF",50,0,1000);
398 hCaloEmEtInHF = ibooker.
book1D(
"CaloEmEtInHF" ,
"CaloEmEtInHF" ,25,0,500);
400 hCaloEmEtInEE = ibooker.
book1D(
"CaloEmEtInEE" ,
"CaloEmEtInEE" ,50,0,1000);
402 hCaloEmEtInEB = ibooker.
book1D(
"CaloEmEtInEB" ,
"CaloEmEtInEB" ,50,0,2000);
405 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CaloHadEtInHO",hCaloHadEtInHO));
406 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CaloHadEtInHF",hCaloHadEtInHF));
407 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CaloHadEtInHE",hCaloHadEtInHE));
408 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CaloHadEtInHB",hCaloHadEtInHB));
409 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CaloEmEtInHF",hCaloEmEtInHF));
410 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CaloEmEtInEE",hCaloEmEtInEE));
411 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CaloEmEtInEB",hCaloEmEtInEB));
413 hCaloMETPhi020 = ibooker.
book1D(
"CaloMETPhi020",
"CaloMETPhi020", 60, -
M_PI,
M_PI);
414 hCaloMETPhi020 ->
setAxisTitle(
"METPhi [rad] (MET>20 GeV)", 1);
416 hCaloEtFractionHadronic = ibooker.
book1D(
"CaloEtFractionHadronic",
"CaloEtFractionHadronic",50,0,1);
417 hCaloEtFractionHadronic->
setAxisTitle(
"Hadronic Et Fraction",1);
418 hCaloEmEtFraction = ibooker.
book1D(
"CaloEmEtFraction",
"CaloEmEtFraction" ,50,0,1);
421 hCaloEmEtFraction020 = ibooker.
book1D(
"CaloEmEtFraction020",
"CaloEmEtFraction020" ,50,0,1);
422 hCaloEmEtFraction020->
setAxisTitle(
"EM Et Fraction (MET>20 GeV)",1);
425 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CaloMETPhi020",hCaloMETPhi020));
426 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CaloEtFractionHadronic",hCaloEtFractionHadronic));
427 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CaloEmEtFraction", hCaloEmEtFraction));
428 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CaloEmEtFraction020",hCaloEmEtFraction020));
433 if(fillPFCandPlots && fillCandidateMap_histos){
458 meCHF_Barrel = ibooker.
book1D(
"PfChargedHadronEtFractionBarrel",
"chargedHadronEtFractionBarrel", 50, 0, 1);
459 meCHF_EndcapPlus = ibooker.
book1D(
"PfChargedHadronEtFractionEndcapPlus",
"chargedHadronEtFractionEndcapPlus", 50, 0, 1);
460 meCHF_EndcapMinus = ibooker.
book1D(
"PfChargedHadronEtFractionEndcapMinus",
"chargedHadronEtFractionEndcapMinus", 50, 0, 1);
461 meCHF_Barrel_BXm1Empty = ibooker.
book1D(
"PfChargedHadronEtFractionBarrel_BXm1Empty",
"chargedHadronEtFractionBarrel prev empty bunch", 50, 0, 1);
462 meCHF_EndcapPlus_BXm1Empty = ibooker.
book1D(
"PfChargedHadronEtFractionEndcapPlus_BXm1Empty",
"chargedHadronEtFractionEndcapPlus prev empty bunch", 50, 0, 1);
463 meCHF_EndcapMinus_BXm1Empty = ibooker.
book1D(
"PfChargedHadronEtFractionEndcapMinus_BXm1Empty",
"chargedHadronEtFractionEndcapMinus prev empty bunch", 50, 0, 1);
467 meCHF_Barrel_BXm1Filled = ibooker.
book1D(
"PfChargedHadronEtFractionBarrel_BXm1Filled",
"chargedHadronEtFractionBarrel prev filled 2 bunches", 50, 0, 1);
468 meCHF_EndcapPlus_BXm1Filled = ibooker.
book1D(
"PfChargedHadronEtFractionEndcapPlus_BXm1Filled",
"chargedHadronEtFractionEndcapPlus prev filled bunch", 50, 0, 1);
469 meCHF_EndcapMinus_BXm1Filled = ibooker.
book1D(
"PfChargedHadronEtFractionEndcapMinus_BXm1Filled",
"chargedHadronEtFractionEndcapMinus prev filled bunch", 50, 0, 1);
474 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEtFractionBarrel", meCHF_Barrel));
475 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEtFractionEndcapPlus", meCHF_EndcapPlus));
476 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEtFractionEndcapMinus", meCHF_EndcapMinus));
477 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEtFractionBarrel_BXm1Empty", meCHF_Barrel_BXm1Empty));
478 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEtFractionEndcapPlus_BXm1Empty", meCHF_EndcapPlus_BXm1Empty));
479 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEtFractionEndcapMinus_BXm1Empty", meCHF_EndcapMinus_BXm1Empty));
483 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEtFractionBarrel_BXm1Filled", meCHF_Barrel_BXm1Filled));
484 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEtFractionEndcapPlus_BXm1Filled", meCHF_EndcapPlus_BXm1Filled));
485 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEtFractionEndcapMinus_BXm1Filled", meCHF_EndcapMinus_BXm1Filled));
490 meNHF_Barrel = ibooker.
book1D(
"PfNeutralHadronEtFractionBarrel",
"neutralHadronEtFractionBarrel", 50, 0, 1);
491 meNHF_EndcapPlus = ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapPlus",
"neutralHadronEtFractionEndcapPlus", 50, 0, 1);
492 meNHF_EndcapMinus = ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapMinus",
"neutralHadronEtFractionEndcapMinus", 50, 0, 1);
493 meNHF_Barrel_BXm1Empty = ibooker.
book1D(
"PfNeutralHadronEtFractionBarrel_BXm1Empty",
"neutralHadronEtFractionBarrel prev empty bunch", 50, 0, 1);
494 meNHF_EndcapPlus_BXm1Empty = ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapPlus_BXm1Empty",
"neutralHadronEtFractionEndcapPlus prev empty bunch", 50, 0, 1);
495 meNHF_EndcapMinus_BXm1Empty = ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapMinus_BXm1Empty",
"neutralHadronEtFractionEndcapMinus prev empty bunch", 50, 0, 1);
499 meNHF_Barrel_BXm1Filled = ibooker.
book1D(
"PfNeutralHadronEtFractionBarrel_BXm1Filled",
"neutralHadronEtFractionBarrel prev filled 2 bunches", 50, 0, 1);
500 meNHF_EndcapPlus_BXm1Filled = ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapPlus_BXm1Filled",
"neutralHadronEtFractionEndcapPlus prev filled bunch", 50, 0, 1);
501 meNHF_EndcapMinus_BXm1Filled = ibooker.
book1D(
"PfNeutralHadronEtFractionEndcapMinus_BXm1Filled",
"neutralHadronEtFractionEndcapMinus prev filled bunch", 50, 0, 1);
506 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEtFractionBarrel", meNHF_Barrel));
507 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEtFractionEndcapPlus", meNHF_EndcapPlus));
508 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEtFractionEndcapMinus", meNHF_EndcapMinus));
509 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEtFractionBarrel_BXm1Empty", meNHF_Barrel_BXm1Empty));
510 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEtFractionEndcapPlus_BXm1Empty", meNHF_EndcapPlus_BXm1Empty));
511 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEtFractionEndcapMinus_BXm1Empty", meNHF_EndcapMinus_BXm1Empty));
515 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEtFractionBarrel_BXm1Filled", meNHF_Barrel_BXm1Filled));
516 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEtFractionEndcapPlus_BXm1Filled", meNHF_EndcapPlus_BXm1Filled));
517 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEtFractionEndcapMinus_BXm1Filled", meNHF_EndcapMinus_BXm1Filled));
522 mePhF_Barrel = ibooker.
book1D(
"PfPhotonEtFractionBarrel",
"photonEtFractionBarrel", 50, 0, 1);
523 mePhF_EndcapPlus = ibooker.
book1D(
"PfPhotonEtFractionEndcapPlus",
"photonEtFractionEndcapPlus", 50, 0, 1);
524 mePhF_EndcapMinus = ibooker.
book1D(
"PfPhotonEtFractionEndcapMinus",
"photonEtFractionEndcapMinus", 50, 0, 1);
525 mePhF_Barrel_BXm1Empty = ibooker.
book1D(
"PfPhotonEtFractionBarrel_BXm1Empty",
"photonEtFractionBarrel prev empty bunch", 50, 0, 1);
526 mePhF_EndcapPlus_BXm1Empty = ibooker.
book1D(
"PfPhotonEtFractionEndcapPlus_BXm1Empty",
"photonEtFractionEndcapPlus prev empty bunch", 50, 0, 1);
527 mePhF_EndcapMinus_BXm1Empty = ibooker.
book1D(
"PfPhotonEtFractionEndcapMinus_BXm1Empty",
"photonEtFractionEndcapMinus prev empty bunch", 50, 0, 1);
531 mePhF_Barrel_BXm1Filled = ibooker.
book1D(
"PfPhotonEtFractionBarrel_BXm1Filled",
"photonEtFractionBarrel prev filled 2 bunches", 50, 0, 1);
532 mePhF_EndcapPlus_BXm1Filled = ibooker.
book1D(
"PfPhotonEtFractionEndcapPlus_BXm1Filled",
"photonEtFractionEndcapPlus prev filled bunch", 50, 0, 1);
533 mePhF_EndcapMinus_BXm1Filled = ibooker.
book1D(
"PfPhotonEtFractionEndcapMinus_BXm1Filled",
"photonEtFractionEndcapMinus prev filled bunch", 50, 0, 1);
538 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFractionBarrel", mePhF_Barrel));
539 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFractionEndcapPlus", mePhF_EndcapPlus));
540 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFractionEndcapMinus", mePhF_EndcapMinus));
541 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFractionBarrel_BXm1Empty", mePhF_Barrel_BXm1Empty));
542 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFractionEndcapPlus_BXm1Empty", mePhF_EndcapPlus_BXm1Empty));
543 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFractionEndcapMinus_BXm1Empty", mePhF_EndcapMinus_BXm1Empty));
547 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFractionBarrel_BXm1Filled", mePhF_Barrel_BXm1Filled));
548 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFractionEndcapPlus_BXm1Filled", mePhF_EndcapPlus_BXm1Filled));
549 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFractionEndcapMinus_BXm1Filled", mePhF_EndcapMinus_BXm1Filled));
554 meHFHadF_Plus = ibooker.
book1D(
"PfHFHadronEtFractionPlus",
"HFHadronEtFractionPlus", 50, 0, 1);
555 meHFHadF_Minus = ibooker.
book1D(
"PfHFHadronEtFractionMinus",
"HFHadronEtFractionMinus", 50, 0, 1);
556 meHFHadF_Plus_BXm1Empty = ibooker.
book1D(
"PfHFHadronEtFractionPlus_BXm1Empty",
"HFHadronEtFractionPlus prev empty bunch", 50, 0, 1);
557 meHFHadF_Minus_BXm1Empty = ibooker.
book1D(
"PfHFHadronEtFractionMinus_BXm1Empty",
"HFHadronEtFractionMinus prev empty bunch", 50, 0, 1);
560 meHFHadF_Plus_BXm1Filled = ibooker.
book1D(
"PfHFHadronEtFractionPlus_BXm1Filled",
"HFHadronEtFractionPlus prev filled bunch", 50, 0, 1);
561 meHFHadF_Minus_BXm1Filled = ibooker.
book1D(
"PfHFHadronEtFractionMinus_BXm1Filled",
"HFHadronEtFractionMinus prev filled bunch", 50, 0, 1);
565 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFHadronEtFractionPlus", meHFHadF_Plus));
566 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFHadronEtFractionMinus", meHFHadF_Minus));
567 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFHadronEtFractionPlus_BXm1Empty", meHFHadF_Plus_BXm1Empty));
568 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFHadronEtFractionMinus_BXm1Empty", meHFHadF_Minus_BXm1Empty));
571 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFHadronEtFractionPlus_BXm1Filled", meHFHadF_Plus_BXm1Filled));
572 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFHadronEtFractionMinus_BXm1Filled", meHFHadF_Minus_BXm1Filled));
576 meHFEMF_Plus = ibooker.
book1D(
"PfHFEMEtFractionPlus",
"HFEMEtFractionPlus", 50, 0, 1);
577 meHFEMF_Minus = ibooker.
book1D(
"PfHFEMEtFractionMinus",
"HFEMEtFractionMinus", 50, 0, 1);
578 meHFEMF_Plus_BXm1Empty = ibooker.
book1D(
"PfHFEMEtFractionPlus_BXm1Empty",
"HFEMEtFractionPlus prev empty bunch", 50, 0, 1);
579 meHFEMF_Minus_BXm1Empty = ibooker.
book1D(
"PfHFEMEtFractionMinus_BXm1Empty",
"HFEMEtFractionMinus prev empty bunch", 50, 0, 1);
582 meHFEMF_Plus_BXm1Filled = ibooker.
book1D(
"PfHFEMEtFractionPlus_BXm1Filled",
"HFEMEtFractionPlus prev filled bunch", 50, 0, 1);
583 meHFEMF_Minus_BXm1Filled = ibooker.
book1D(
"PfHFEMEtFractionMinus_BXm1Filled",
"HFEMEtFractionMinus prev filled bunch", 50, 0, 1);
587 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFEMEtFractionPlus", meHFEMF_Plus));
588 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFEMEtFractionMinus", meHFEMF_Minus));
589 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFEMEtFractionPlus_BXm1Empty", meHFEMF_Plus_BXm1Empty));
590 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFEMEtFractionMinus_BXm1Empty", meHFEMF_Minus_BXm1Empty));
593 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFEMEtFractionPlus_BXm1Filled", meHFEMF_Plus_BXm1Filled));
594 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFEMEtFractionMinus_BXm1Filled", meHFEMF_Minus_BXm1Filled));
655 mePhotonEtFraction_BXm1Empty = ibooker.
book1D(
"PfPhotonEtFraction_BXm1Empty",
"photonEtFraction() prev empty bunch", 50, 0, 1);
656 mePhotonEtFraction_BXm1Filled = ibooker.
book1D(
"PfPhotonEtFraction_BXm1Filled",
"photonEtFraction() prev filled bunch", 50, 0, 1);
657 meNeutralHadronEtFraction_BXm1Empty = ibooker.
book1D(
"PfNeutralHadronEtFraction_BXm1Empty",
"neutralHadronEtFraction() prev empty bunch", 50, 0, 1);
658 meNeutralHadronEtFraction_BXm1Filled = ibooker.
book1D(
"PfNeutralHadronEtFraction_BXm1Filled",
"neutralHadronEtFraction() prev filled bunch", 50, 0, 1);
659 meChargedHadronEtFraction_BXm1Empty = ibooker.
book1D(
"PfChargedHadronEtFraction_BXm1Empty",
"chargedHadronEtFraction() prev empty bunch", 50, 0, 1);
660 meChargedHadronEtFraction_BXm1Filled = ibooker.
book1D(
"PfChargedHadronEtFraction_BXm1Filled",
"chargedHadronEtFraction() prev filled bunch", 50, 0, 1);
661 meMET_BXm1Empty = ibooker.
book1D(
"MET_BXm1Empty",
"MET prev empty bunch", 200, 0, 1000);
662 meMET_BXm1Filled = ibooker.
book1D(
"MET_BXm1Filled",
"MET prev filled bunch", 200, 0, 1000);
663 meSumET_BXm1Empty = ibooker.
book1D(
"SumET_BXm1Empty",
"SumET prev empty bunch", 400, 0, 4000);
664 meSumET_BXm1Filled = ibooker.
book1D(
"SumET_BXm1Filled",
"SumET prev filled bunch", 400, 0, 4000);
666 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFraction_BXm1Empty" ,mePhotonEtFraction_BXm1Empty));
667 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFraction_BXm1Filled" ,mePhotonEtFraction_BXm1Filled));
668 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEtFraction_BXm1Empty" ,meNeutralHadronEtFraction_BXm1Empty));
669 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralEtFraction_BXm1Filled" ,meNeutralHadronEtFraction_BXm1Filled));
670 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEtFraction_BXm1Empty" ,meChargedHadronEtFraction_BXm1Empty));
671 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedEtFraction_BXm1Filled" ,meChargedHadronEtFraction_BXm1Filled));
672 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_BXm1Empty" ,meMET_BXm1Empty));
673 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MET_BXm1Filled" ,meMET_BXm1Filled));
674 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"SumET_BXm1Empty" ,meSumET_BXm1Empty));
675 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"SumET_BXm1Filled",meSumET_BXm1Filled));
677 meMETPhiChargedHadronsBarrel_BXm1Filled = ibooker.
book1D(
"METPhiChargedHadronsBarrel_BXm1Filled",
"METPhi_PFChargedHadronsBarrel prev bunch filled", 50, -
M_PI,
M_PI);
678 meMETPhiChargedHadronsEndcapPlus_BXm1Filled = ibooker.
book1D(
"METPhiChargedHadronsEndcapPlus_BXm1Filled",
"METPhi_PFChargedHadronsEndcapPlus prev bunch filled", 50, -
M_PI,
M_PI);
679 meMETPhiChargedHadronsEndcapMinus_BXm1Filled = ibooker.
book1D(
"METPhiChargedHadronsEndcapMinus_BXm1Filled",
"METPhi_PFChargedHadronsEndcapMinus prev bunch filled", 50, -
M_PI,
M_PI);
680 meMETPhiNeutralHadronsBarrel_BXm1Filled = ibooker.
book1D(
"METPhiNeutralHadronsBarrel_BXm1Filled",
"METPhi_PFNeutralHadronsBarrel prev bunch filled", 50, -
M_PI,
M_PI);
681 meMETPhiNeutralHadronsEndcapPlus_BXm1Filled = ibooker.
book1D(
"METPhiNeutralHadronsEndcapPlus_BXm1Filled",
"METPhi_PFNeutralHadronsEndcapPlus prev bunch filled", 50, -
M_PI,
M_PI);
682 meMETPhiNeutralHadronsEndcapMinus_BXm1Filled = ibooker.
book1D(
"METPhiNeutralHadronsEndcapMinus_BXm1Filled",
"METPhi_PFNeutralHadronsEndcapMinus prev bunch filled", 50, -
M_PI,
M_PI);
683 meMETPhiPhotonsBarrel_BXm1Filled = ibooker.
book1D(
"METPhiPhotonsBarrel_BXm1Filled",
"METPhi_PFPhotonsBarrel prev bunch filled", 50, -
M_PI,
M_PI);
684 meMETPhiPhotonsEndcapPlus_BXm1Filled = ibooker.
book1D(
"METPhiPhotonsEndcapPlus_BXm1Filled",
"METPhi_PFPhotonsEndcapPlus prev bunch filled", 50, -
M_PI,
M_PI);
685 meMETPhiPhotonsEndcapMinus_BXm1Filled = ibooker.
book1D(
"METPhiPhotonsEndcapMinus_BXm1Filled",
"METPhi_PFPhotonsEndcapMinus prev bunch filled", 50, -
M_PI,
M_PI);
686 meMETPhiHFHadronsPlus_BXm1Filled = ibooker.
book1D(
"METPhiHFHadronsPlus_BXm1Filled",
"METPhi_PFHFHadronsPlus prev bunch filled", 50, -
M_PI,
M_PI);
687 meMETPhiHFHadronsMinus_BXm1Filled = ibooker.
book1D(
"METPhiHFHadronsMinus_BXm1Filled",
"METPhi_PFHFHadronsMinus prev bunch filled", 50, -
M_PI,
M_PI);
688 meMETPhiHFEGammasPlus_BXm1Filled = ibooker.
book1D(
"METPhiHFEGammasPlus_BXm1Filled",
"METPhi_PFHFEGammasPlus prev bunch filled", 50, -
M_PI,
M_PI);
689 meMETPhiHFEGammasMinus_BXm1Filled = ibooker.
book1D(
"METPhiHFEGammasMinus_BXm1Filled",
"METPhi_PFHFEGammasMinus prev bunch filled", 50, -
M_PI,
M_PI);
691 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiChargedHadronsBarrel_BXm1Filled" ,meMETPhiChargedHadronsBarrel_BXm1Filled));
692 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiChargedHadronsEndcapPlus_BXm1Filled" ,meMETPhiChargedHadronsEndcapPlus_BXm1Filled));
693 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiChargedHadronsEndcapMinus_BXm1Filled" ,meMETPhiChargedHadronsEndcapMinus_BXm1Filled));
694 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiNeutralHadronsBarrel_BXm1Filled" ,meMETPhiNeutralHadronsBarrel_BXm1Filled));
695 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiNeutralHadronsEndcapPlus_BXm1Filled" ,meMETPhiNeutralHadronsEndcapPlus_BXm1Filled));
696 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiNeutralHadronsEndcapMinus_BXm1Filled" ,meMETPhiNeutralHadronsEndcapMinus_BXm1Filled));
697 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiPhotonsBarrel_BXm1Filled" ,meMETPhiPhotonsBarrel_BXm1Filled));
698 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiPhotonsEndcapPlus_BXm1Filled" ,meMETPhiPhotonsEndcapPlus_BXm1Filled));
699 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiPhotonsEndcapMinus_BXm1Filled" ,meMETPhiPhotonsEndcapMinus_BXm1Filled));
700 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiHFHadronsPlus_BXm1Filled" ,meMETPhiHFHadronsPlus_BXm1Filled));
701 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiHFHadronsMinus_BXm1Filled" ,meMETPhiHFHadronsMinus_BXm1Filled));
702 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiHFEGammasPlus_BXm1Filled" ,meMETPhiHFEGammasPlus_BXm1Filled));
703 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiHFEGammasMinus_BXm1Filled" ,meMETPhiHFEGammasMinus_BXm1Filled));
705 meMETPhiChargedHadronsBarrel_BXm1Empty = ibooker.
book1D(
"METPhiChargedHadronsBarrel_BXm1Empty",
"METPhi_PFChargedHadronsBarrel prev bunch empty", 50, -
M_PI,
M_PI);
706 meMETPhiChargedHadronsEndcapPlus_BXm1Empty = ibooker.
book1D(
"METPhiChargedHadronsEndcapPlus_BXm1Empty",
"METPhi_PFChargedHadronsEndcapPlus prev bunch empty", 50, -
M_PI,
M_PI);
707 meMETPhiChargedHadronsEndcapMinus_BXm1Empty = ibooker.
book1D(
"METPhiChargedHadronsEndcapMinus_BXm1Empty",
"METPhi_PFChargedHadronsEndcapMinus prev bunch empty", 50, -
M_PI,
M_PI);
708 meMETPhiNeutralHadronsBarrel_BXm1Empty = ibooker.
book1D(
"METPhiNeutralHadronsBarrel_BXm1Empty",
"METPhi_PFNeutralHadronsBarrel prev bunch empty", 50, -
M_PI,
M_PI);
709 meMETPhiNeutralHadronsEndcapPlus_BXm1Empty = ibooker.
book1D(
"METPhiNeutralHadronsEndcapPlus_BXm1Empty",
"METPhi_PFNeutralHadronsEndcapPlus prev bunch empty", 50, -
M_PI,
M_PI);
710 meMETPhiNeutralHadronsEndcapMinus_BXm1Empty = ibooker.
book1D(
"METPhiNeutralHadronsEndcapMinus_BXm1Empty",
"METPhi_PFNeutralHadronsEndcapMinus prev bunch empty", 50, -
M_PI,
M_PI);
711 meMETPhiPhotonsBarrel_BXm1Empty = ibooker.
book1D(
"METPhiPhotonsBarrel_BXm1Empty",
"METPhi_PFPhotonsBarrel prev bunch empty", 50, -
M_PI,
M_PI);
712 meMETPhiPhotonsEndcapPlus_BXm1Empty = ibooker.
book1D(
"METPhiPhotonsEndcapPlus_BXm1Empty",
"METPhi_PFPhotonsEndcapPlus prev bunch empty", 50, -
M_PI,
M_PI);
713 meMETPhiPhotonsEndcapMinus_BXm1Empty = ibooker.
book1D(
"METPhiPhotonsEndcapMinus_BXm1Empty",
"METPhi_PFPhotonsEndcapMinus prev bunch empty", 50, -
M_PI,
M_PI);
714 meMETPhiHFHadronsPlus_BXm1Empty = ibooker.
book1D(
"METPhiHFHadronsPlus_BXm1Empty",
"METPhi_PFHFHadronsPlus prev bunch empty", 50, -
M_PI,
M_PI);
715 meMETPhiHFHadronsMinus_BXm1Empty = ibooker.
book1D(
"METPhiHFHadronsMinus_BXm1Empty",
"METPhi_PFHFHadronsMinus prev bunch empty", 50, -
M_PI,
M_PI);
716 meMETPhiHFEGammasPlus_BXm1Empty = ibooker.
book1D(
"METPhiHFEGammasPlus_BXm1Empty",
"METPhi_PFHFEGammasPlus prev bunch empty", 50, -
M_PI,
M_PI);
717 meMETPhiHFEGammasMinus_BXm1Empty = ibooker.
book1D(
"METPhiHFEGammasMinus_BXm1Empty",
"METPhi_PFHFEGammasMinus prev bunch empty", 50, -
M_PI,
M_PI);
719 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiChargedHadronsBarrel_BXm1Empty" ,meMETPhiChargedHadronsBarrel_BXm1Empty));
720 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiChargedHadronsEndcapPlus_BXm1Empty" ,meMETPhiChargedHadronsEndcapPlus_BXm1Empty));
721 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiChargedHadronsEndcapMinus_BXm1Empty" ,meMETPhiChargedHadronsEndcapMinus_BXm1Empty));
722 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiNeutralHadronsBarrel_BXm1Empty" ,meMETPhiNeutralHadronsBarrel_BXm1Empty));
723 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiNeutralHadronsEndcapPlus_BXm1Empty" ,meMETPhiNeutralHadronsEndcapPlus_BXm1Empty));
724 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiNeutralHadronsEndcapMinus_BXm1Empty" ,meMETPhiNeutralHadronsEndcapMinus_BXm1Empty));
725 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiPhotonsBarrel_BXm1Empty" ,meMETPhiPhotonsBarrel_BXm1Empty));
726 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiPhotonsEndcapPlus_BXm1Empty" ,meMETPhiPhotonsEndcapPlus_BXm1Empty));
727 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiPhotonsEndcapMinus_BXm1Empty" ,meMETPhiPhotonsEndcapMinus_BXm1Empty));
728 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiHFHadronsPlus_BXm1Empty" ,meMETPhiHFHadronsPlus_BXm1Empty));
729 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiHFHadronsMinus_BXm1Empty" ,meMETPhiHFHadronsMinus_BXm1Empty));
730 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiHFEGammasPlus_BXm1Empty" ,meMETPhiHFEGammasPlus_BXm1Empty));
731 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiHFEGammasMinus_BXm1Empty" ,meMETPhiHFEGammasMinus_BXm1Empty));
733 meMETPhiChargedHadronsBarrel = ibooker.
book1D(
"METPhiChargedHadronsBarrel",
"METPhi_PFChargedHadronsBarrel", 50, -
M_PI,
M_PI);
734 meMETPhiChargedHadronsEndcapPlus = ibooker.
book1D(
"METPhiChargedHadronsEndcapPlus",
"METPhi_PFChargedHadronsEndcapPlus", 50, -
M_PI,
M_PI);
735 meMETPhiChargedHadronsEndcapMinus = ibooker.
book1D(
"METPhiChargedHadronsEndcapMinus",
"METPhi_PFChargedHadronsEndcapMinus", 50, -
M_PI,
M_PI);
736 meMETPhiNeutralHadronsBarrel = ibooker.
book1D(
"METPhiNeutralHadronsBarrel",
"METPhi_PFNeutralHadronsBarrel", 50, -
M_PI,
M_PI);
737 meMETPhiNeutralHadronsEndcapPlus = ibooker.
book1D(
"METPhiNeutralHadronsEndcapPlus",
"METPhi_PFNeutralHadronsEndcapPlus", 50, -
M_PI,
M_PI);
738 meMETPhiNeutralHadronsEndcapMinus = ibooker.
book1D(
"METPhiNeutralHadronsEndcapMinus",
"METPhi_PFNeutralHadronsEndcapMinus", 50, -
M_PI,
M_PI);
739 meMETPhiPhotonsBarrel = ibooker.
book1D(
"METPhiPhotonsBarrel",
"METPhi_PFPhotonsBarrel", 50, -
M_PI,
M_PI);
740 meMETPhiPhotonsEndcapPlus = ibooker.
book1D(
"METPhiPhotonsEndcapPlus",
"METPhi_PFPhotonsEndcapPlus", 50, -
M_PI,
M_PI);
741 meMETPhiPhotonsEndcapMinus = ibooker.
book1D(
"METPhiPhotonsEndcapMinus",
"METPhi_PFPhotonsEndcapMinus", 50, -
M_PI,
M_PI);
742 meMETPhiHFHadronsPlus = ibooker.
book1D(
"METPhiHFHadronsPlus",
"METPhi_PFHFHadronsPlus", 50, -
M_PI,
M_PI);
743 meMETPhiHFHadronsMinus = ibooker.
book1D(
"METPhiHFHadronsMinus",
"METPhi_PFHFHadronsMinus", 50, -
M_PI,
M_PI);
744 meMETPhiHFEGammasPlus = ibooker.
book1D(
"METPhiHFEGammasPlus",
"METPhi_PFHFEGammasPlus", 50, -
M_PI,
M_PI);
745 meMETPhiHFEGammasMinus = ibooker.
book1D(
"METPhiHFEGammasMinus",
"METPhi_PFHFEGammasMinus", 50, -
M_PI,
M_PI);
747 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiChargedHadronsBarrel" ,meMETPhiChargedHadronsBarrel));
748 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiChargedHadronsEndcapPlus" ,meMETPhiChargedHadronsEndcapPlus));
749 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiChargedHadronsEndcapMinus" ,meMETPhiChargedHadronsEndcapMinus));
750 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiNeutralHadronsBarrel" ,meMETPhiNeutralHadronsBarrel));
751 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiNeutralHadronsEndcapPlus" ,meMETPhiNeutralHadronsEndcapPlus));
752 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiNeutralHadronsEndcapMinus" ,meMETPhiNeutralHadronsEndcapMinus));
753 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiPhotonsBarrel" ,meMETPhiPhotonsBarrel));
754 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiPhotonsEndcapPlus" ,meMETPhiPhotonsEndcapPlus));
755 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiPhotonsEndcapMinus" ,meMETPhiPhotonsEndcapMinus));
756 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiHFHadronsPlus" ,meMETPhiHFHadronsPlus));
757 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiHFHadronsMinus" ,meMETPhiHFHadronsMinus));
758 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiHFEGammasPlus" ,meMETPhiHFEGammasPlus));
759 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METPhiHFEGammasMinus" ,meMETPhiHFEGammasMinus));
762 if(fillPFCandPlots && fillCandidateMap_histos){
763 if(!etaMinPFCand_.empty()){
764 etaMinPFCand_.clear();
765 etaMaxPFCand_.clear();
767 nbinsPFCand_.clear();
768 countsPFCand_.clear();
771 profilePFCand_x_.clear();
772 profilePFCand_y_.clear();
773 occupancyPFCand_.clear();
775 multiplicityPFCand_.clear();
776 profilePFCand_x_name_.clear();
777 profilePFCand_y_name_.clear();
778 occupancyPFCand_name_.clear();
779 ptPFCand_name_.clear();
780 multiplicityPFCand_name_.clear();
782 for (std::vector<edm::ParameterSet>::const_iterator
v = diagnosticsParameters_.begin();
v!=diagnosticsParameters_.end();
v++) {
783 int etaNBinsPFCand =
v->getParameter<
int>(
"etaNBins");
784 double etaMinPFCand =
v->getParameter<
double>(
"etaMin");
785 double etaMaxPFCand =
v->getParameter<
double>(
"etaMax");
786 int phiNBinsPFCand =
v->getParameter<
int>(
"phiNBins");
787 double phiMinPFCand =
v->getParameter<
double>(
"phiMin");
788 double phiMaxPFCand =
v->getParameter<
double>(
"phiMax");
789 int nMinPFCand =
v->getParameter<
int>(
"nMin");
790 int nMaxPFCand =
v->getParameter<
int>(
"nMax");
791 int nbinsPFCand =
v->getParameter<
double>(
"nbins");
794 etaMinPFCand_.push_back(etaMinPFCand);
795 etaMaxPFCand_.push_back(etaMaxPFCand);
796 nbinsPFCand_.push_back(nbinsPFCand);
797 typePFCand_.push_back(
v->getParameter<
int>(
"type"));
798 countsPFCand_.push_back(0);
799 MExPFCand_.push_back(0.);
800 MEyPFCand_.push_back(0.);
804 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+profilePFCand_x_name_[profilePFCand_x_name_.size()-1], profilePFCand_x_[profilePFCand_x_.size()-1]));
807 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+profilePFCand_y_name_[profilePFCand_y_name_.size()-1], profilePFCand_y_[profilePFCand_y_.size()-1]));
809 multiplicityPFCand_name_.push_back(
std::string(
v->getParameter<
std::string>(
"name")).append(
"_multiplicity_").c_str());
810 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+ multiplicityPFCand_name_[multiplicityPFCand_name_.size()-1], multiplicityPFCand_[multiplicityPFCand_.size()-1]));
813 occupancyPFCand_name_.push_back(
std::string(
v->getParameter<
std::string>(
"name")).append(
"_occupancy_").c_str());
819 if(etaMinPFCand*etaMaxPFCand<0){
820 occupancyPFCand_.push_back(ibooker.
book2D(
std::string(
v->getParameter<
std::string>(
"name")).append(
"_occupancy_").c_str(),
std::string(
v->getParameter<
std::string>(
"name"))+
"occupancy", etaNBinsPFCand, etaMinPFCand, etaMaxPFCand, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand));
823 const int nbins_eta_endcap=2*(etaNBinsPFCand+1);
824 double eta_limits_endcap[nbins_eta_endcap];
825 for(
int i=0;
i<nbins_eta_endcap;
i++){
826 if(
i<(etaNBinsPFCand+1)){
827 eta_limits_endcap[
i]=etaMinPFCand+
i*(etaMaxPFCand-etaMinPFCand)/(
double)etaNBinsPFCand;
829 eta_limits_endcap[
i]= -etaMaxPFCand +(
i- (etaNBinsPFCand+1) )*(etaMaxPFCand-etaMinPFCand)/(
double)etaNBinsPFCand;
832 TH2F* hist_temp_occup =
new TH2F((occupancyPFCand_name_[occupancyPFCand_name_.size()-1]).c_str(),
"occupancy",nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
833 occupancyPFCand_.push_back(ibooker.
book2D(occupancyPFCand_name_[occupancyPFCand_name_.size()-1],hist_temp_occup));
834 TH2F* hist_temp_pt =
new TH2F((ptPFCand_name_[ptPFCand_name_.size()-1]).c_str(),
"pt",nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
835 ptPFCand_.push_back(ibooker.
book2D(ptPFCand_name_[ptPFCand_name_.size()-1], hist_temp_pt));
838 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+ occupancyPFCand_name_[occupancyPFCand_name_.size()-1], occupancyPFCand_[occupancyPFCand_.size()-1]));
839 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+ ptPFCand_name_[ptPFCand_name_.size()-1], ptPFCand_[ptPFCand_.size()-1]));
842 mProfileIsoPFChHad_TrackOccupancy=ibooker.
book2D(
"IsoPfChHad_Track_profile",
"Isolated PFChHadron Tracker_occupancy", 108, -2.7, 2.7, 160, -
M_PI,
M_PI);
843 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"IsoPfChHad_Track_profile" ,mProfileIsoPFChHad_TrackOccupancy));
844 mProfileIsoPFChHad_TrackPt=ibooker.
book2D(
"IsoPfChHad_TrackPt",
"Isolated PFChHadron TrackPt", 108, -2.7, 2.7, 160, -
M_PI,
M_PI);
845 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"IsoPfChHad_TrackPt" ,mProfileIsoPFChHad_TrackPt));
848 mProfileIsoPFChHad_EcalOccupancyCentral = ibooker.
book2D(
"IsoPfChHad_ECAL_profile_central",
"IsolatedPFChHa ECAL occupancy (Barrel)", 180, -1.479, 1.479, 125, -
M_PI,
M_PI);
849 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"IsoPfChHad_ECAL_profile_central" ,mProfileIsoPFChHad_EcalOccupancyCentral));
850 mProfileIsoPFChHad_EMPtCentral=ibooker.
book2D(
"IsoPfChHad_EMPt_central",
"Isolated PFChHadron HadPt_central", 180, -1.479, 1.479, 360, -
M_PI,
M_PI);
851 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"IsoPfChHad_EMPt_central" ,mProfileIsoPFChHad_EMPtCentral));
853 mProfileIsoPFChHad_EcalOccupancyEndcap = ibooker.
book2D(
"IsoPfChHad_ECAL_profile_endcap",
"IsolatedPFChHa ECAL occupancy (Endcap)", 110, -2.75, 2.75, 125, -
M_PI,
M_PI);
854 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"IsoPfChHad_ECAL_profile_endcap" ,mProfileIsoPFChHad_EcalOccupancyEndcap));
855 mProfileIsoPFChHad_EMPtEndcap=ibooker.
book2D(
"IsoPfChHad_EMPt_endcap",
"Isolated PFChHadron EMPt_endcap", 110, -2.75, 2.75, 125, -
M_PI,
M_PI);
856 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"IsoPfChHad_EMPt_endcap" ,mProfileIsoPFChHad_EMPtEndcap));
858 const int nbins_eta=16;
860 double eta_limits[nbins_eta]={-2.650,-2.500,-2.322,-2.172,-2.043,-1.930,-1.830,-1.740,1.740,1.830,1.930,2.043,2.172,2.3122,2.500,2.650};
862 TH2F* hist_temp_HCAL =
new TH2F(
"IsoPfChHad_HCAL_profile_endcap",
"IsolatedPFChHa HCAL occupancy (outer endcap)",nbins_eta-1,eta_limits, 36, -
M_PI,
M_PI);
863 TH2F* hist_tempPt_HCAL=(TH2F*)hist_temp_HCAL->Clone(
"Isolated PFCHHadron HadPt (outer endcap)");
865 mProfileIsoPFChHad_HcalOccupancyCentral = ibooker.
book2D(
"IsoPfChHad_HCAL_profile_central",
"IsolatedPFChHa HCAL occupancy (Central Part)", 40, -1.740, 1.740, 72, -
M_PI,
M_PI);
866 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"IsoPfChHad_HCAL_profile_central" ,mProfileIsoPFChHad_HcalOccupancyCentral));
867 mProfileIsoPFChHad_HadPtCentral=ibooker.
book2D(
"IsoPfChHad_HadPt_central",
"Isolated PFChHadron HadPt_central", 40, -1.740, 1.740, 72, -
M_PI,
M_PI);
868 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"IsoPfChHad_HadPt_central" ,mProfileIsoPFChHad_HadPtCentral));
870 mProfileIsoPFChHad_HcalOccupancyEndcap = ibooker.
book2D(
"IsoPfChHad_HCAL_profile_endcap",hist_temp_HCAL);
871 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"IsoPfChHad_HCAL_profile_endcap" ,mProfileIsoPFChHad_HcalOccupancyEndcap));
872 mProfileIsoPFChHad_HadPtEndcap=ibooker.
book2D(
"IsoPfChHad_HadPt_endcap",hist_tempPt_HCAL);
873 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"IsoPfChHad_HadPt_endcap" ,mProfileIsoPFChHad_HadPtEndcap));
877 if(fillPFCandPlots && fillCandidateMap_histos){
878 if(!etaMinPFCand_.empty()){
879 etaMinPFCand_.clear();
880 etaMaxPFCand_.clear();
882 nbinsPFCand_.clear();
883 countsPFCand_.clear();
884 profilePFCand_x_.clear();
885 profilePFCand_y_.clear();
886 occupancyPFCand_.clear();
888 occupancyPFCand_.clear();
890 occupancyPFCand_name_puppiNolepWeight_.clear();
891 ptPFCand_name_puppiNolepWeight_.clear();
892 occupancyPFCand_name_puppiNolepWeight_.clear();
893 ptPFCand_name_puppiNolepWeight_.clear();
894 multiplicityPFCand_name_.clear();
896 for (std::vector<edm::ParameterSet>::const_iterator
v = diagnosticsParameters_.begin();
v!=diagnosticsParameters_.end();
v++) {
897 int etaNBinsPFCand =
v->getParameter<
int>(
"etaNBins");
898 double etaMinPFCand =
v->getParameter<
double>(
"etaMin");
899 double etaMaxPFCand =
v->getParameter<
double>(
"etaMax");
900 int phiNBinsPFCand =
v->getParameter<
int>(
"phiNBins");
901 double phiMinPFCand =
v->getParameter<
double>(
"phiMin");
902 double phiMaxPFCand =
v->getParameter<
double>(
"phiMax");
903 int nbinsPFCand =
v->getParameter<
double>(
"nbins");
906 etaMinPFCand_.push_back(etaMinPFCand);
907 etaMaxPFCand_.push_back(etaMaxPFCand);
908 nbinsPFCand_.push_back(nbinsPFCand);
909 typePFCand_.push_back(
v->getParameter<
int>(
"type"));
910 countsPFCand_.push_back(0);
911 MExPFCand_.push_back(0.);
912 MEyPFCand_.push_back(0.);
915 occupancyPFCand_name_.push_back(
std::string(
v->getParameter<
std::string>(
"name")).append(
"_occupancy_").c_str());
918 occupancyPFCand_name_puppiNolepWeight_.push_back(
std::string(
v->getParameter<
std::string>(
"name")).append(
"_occupancy_puppiNolepWeight_").c_str());
919 ptPFCand_name_puppiNolepWeight_.push_back(
std::string(
v->getParameter<
std::string>(
"name")).append(
"_pt_puppiNolepWeight_").c_str());
923 if(etaMinPFCand*etaMaxPFCand<0){
924 occupancyPFCand_.push_back(ibooker.
book2D(
std::string(
v->getParameter<
std::string>(
"name")).append(
"_occupancy_").c_str(),
std::string(
v->getParameter<
std::string>(
"name"))+
"occupancy", etaNBinsPFCand, etaMinPFCand, etaMaxPFCand, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand));
926 occupancyPFCand_puppiNolepWeight_.push_back(ibooker.
book2D(
std::string(
v->getParameter<
std::string>(
"name")).append(
"_occupancy_puppiNolepWeight_").c_str(),
std::string(
v->getParameter<
std::string>(
"name"))+
"occupancy_puppiNolepWeight", etaNBinsPFCand, etaMinPFCand, etaMaxPFCand, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand));
927 ptPFCand_puppiNolepWeight_.push_back(ibooker.
book2D(
std::string(
v->getParameter<
std::string>(
"name")).append(
"_pt_puppiNolepWeight_").c_str(),
std::string(
v->getParameter<
std::string>(
"name"))+
"pt_puppiNolepWeight", etaNBinsPFCand, etaMinPFCand, etaMaxPFCand, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand));
929 const int nbins_eta_endcap=2*(etaNBinsPFCand+1);
930 double eta_limits_endcap[nbins_eta_endcap];
931 for(
int i=0;
i<nbins_eta_endcap;
i++){
932 if(
i<(etaNBinsPFCand+1)){
933 eta_limits_endcap[
i]=etaMinPFCand+
i*(etaMaxPFCand-etaMinPFCand)/(
double)etaNBinsPFCand;
935 eta_limits_endcap[
i]= -etaMaxPFCand +(
i- (etaNBinsPFCand+1) )*(etaMaxPFCand-etaMinPFCand)/(
double)etaNBinsPFCand;
938 TH2F* hist_temp_occup =
new TH2F((occupancyPFCand_name_[occupancyPFCand_name_.size()-1]).c_str(),
"occupancy",nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
939 occupancyPFCand_.push_back(ibooker.
book2D(occupancyPFCand_name_[occupancyPFCand_name_.size()-1],hist_temp_occup));
940 TH2F* hist_temp_pt =
new TH2F((ptPFCand_name_[ptPFCand_name_.size()-1]).c_str(),
"pt",nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
941 ptPFCand_.push_back(ibooker.
book2D(ptPFCand_name_[ptPFCand_name_.size()-1], hist_temp_pt));
942 TH2F* hist_temp_occup_puppiNolepWeight =
new TH2F((occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size()-1]).c_str(),
"occupancy_puppiNolepWeight",nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
943 occupancyPFCand_puppiNolepWeight_.push_back(ibooker.
book2D(occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size()-1],hist_temp_occup_puppiNolepWeight));
944 TH2F* hist_temp_pt_puppiNolepWeight =
new TH2F((ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size()-1]).c_str(),
"pt_puppiNolepWeight",nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
945 ptPFCand_puppiNolepWeight_.push_back(ibooker.
book2D(ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size()-1], hist_temp_pt_puppiNolepWeight));
947 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+ occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size()-1], occupancyPFCand_puppiNolepWeight_[occupancyPFCand_puppiNolepWeight_.size()-1]));
948 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+ ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size()-1], ptPFCand_puppiNolepWeight_[ptPFCand_puppiNolepWeight_.size()-1]));
949 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+ occupancyPFCand_name_[occupancyPFCand_name_.size()-1], occupancyPFCand_[occupancyPFCand_.size()-1]));
950 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+ ptPFCand_name_[ptPFCand_name_.size()-1], ptPFCand_[ptPFCand_.size()-1]));
955 if(isPFMet_ || isMiniAODMet_){
956 mePhotonEtFraction = ibooker.
book1D(
"PfPhotonEtFraction",
"photonEtFraction()", 50, 0, 1);
957 meNeutralHadronEtFraction = ibooker.
book1D(
"PfNeutralHadronEtFraction",
"neutralHadronEtFraction()", 50, 0, 1);
958 meChargedHadronEtFraction = ibooker.
book1D(
"PfChargedHadronEtFraction",
"chargedHadronEtFraction()", 50, 0, 1);
959 meHFHadronEtFraction = ibooker.
book1D(
"PfHFHadronEtFraction",
"HFHadronEtFraction()", 50, 0, 1);
960 meHFEMEtFraction = ibooker.
book1D(
"PfHFEMEtFraction",
"HFEMEtFraction()", 50, 0, 1);
962 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFraction" ,mePhotonEtFraction));
963 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEtFraction",meNeutralHadronEtFraction));
964 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEtFraction",meChargedHadronEtFraction));
965 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFHadronEtFraction" ,meHFHadronEtFraction));
966 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFEMEtFraction" ,meHFEMEtFraction));
968 mePhotonEtFraction_profile = ibooker.
bookProfile(
"PfPhotonEtFraction_profile",
"photonEtFraction()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1);
969 meNeutralHadronEtFraction_profile = ibooker.
bookProfile(
"PfNeutralHadronEtFraction_profile",
"neutralHadronEtFraction()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1);
970 meChargedHadronEtFraction_profile = ibooker.
bookProfile(
"PfChargedHadronEtFraction_profile",
"chargedHadronEtFraction()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1);
971 meHFHadronEtFraction_profile = ibooker.
bookProfile(
"PfHFHadronEtFraction_profile",
"HFHadronEtFraction()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1);
972 meHFEMEtFraction_profile = ibooker.
bookProfile(
"PfHFEMEtFraction_profile",
"HFEMEtFraction()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1);
974 meNeutralHadronEtFraction_profile->setAxisTitle(
"nvtx", 1);
975 meChargedHadronEtFraction_profile->setAxisTitle(
"nvtx", 1);
976 meHFHadronEtFraction_profile ->setAxisTitle(
"nvtx", 1);
977 meHFEMEtFraction_profile ->setAxisTitle(
"nvtx", 1);
979 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEtFraction_profile" ,mePhotonEtFraction_profile));
980 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEtFraction_profile" ,meNeutralHadronEtFraction_profile));
981 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEtFraction_profile" ,meChargedHadronEtFraction_profile));
982 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFHadronEtFraction_profile" ,meHFHadronEtFraction_profile));
983 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFEMEtFraction_profile" ,meHFEMEtFraction_profile));
985 mePhotonEt = ibooker.
book1D(
"PfPhotonEt",
"photonEt()", 50, 0, 1000);
986 meNeutralHadronEt = ibooker.
book1D(
"PfNeutralHadronEt",
"neutralHadronEt()", 50, 0, 1000);
987 meElectronEt = ibooker.
book1D(
"PfElectronEt",
"electronEt()", 50, 0, 100);
988 meChargedHadronEt = ibooker.
book1D(
"PfChargedHadronEt",
"chargedHadronEt()", 50, 0, 2000);
989 meMuonEt = ibooker.
book1D(
"PfMuonEt",
"muonEt()", 50, 0, 100);
990 meHFHadronEt = ibooker.
book1D(
"PfHFHadronEt",
"HFHadronEt()", 50, 0, 2000);
991 meHFEMEt = ibooker.
book1D(
"PfHFEMEt",
"HFEMEt()", 50, 0, 1000);
993 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEt" ,mePhotonEt));
994 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEt" ,meNeutralHadronEt));
995 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfElectronEt" ,meElectronEt));
996 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEt" ,meChargedHadronEt));
997 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfMuonEt" ,meMuonEt));
998 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFHadronEt" ,meHFHadronEt));
999 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFEMEt" ,meHFEMEt));
1001 mePhotonEt_profile = ibooker.
bookProfile(
"PfPhotonEt_profile",
"photonEt()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1000);
1002 meNeutralHadronEt_profile = ibooker.
bookProfile(
"PfNeutralHadronEt_profile",
"neutralHadronEt()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1000);
1003 meChargedHadronEt_profile = ibooker.
bookProfile(
"PfChargedHadronEt_profile",
"chargedHadronEt()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1000);
1004 meHFHadronEt_profile = ibooker.
bookProfile(
"PfHFHadronEt_profile",
"HFHadronEt()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1000);
1005 meHFEMEt_profile = ibooker.
bookProfile(
"PfHFEMEt_profile",
"HFEMEt()", nbinsPV_, nPVMin_, nPVMax_, 50, 0, 1000);
1008 meNeutralHadronEt_profile ->setAxisTitle(
"nvtx", 1);
1009 meChargedHadronEt_profile ->setAxisTitle(
"nvtx", 1);
1010 meHFHadronEt_profile ->setAxisTitle(
"nvtx", 1);
1011 meHFEMEt_profile ->setAxisTitle(
"nvtx", 1);
1013 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfPhotonEt_profile" ,mePhotonEt_profile));
1014 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfNeutralHadronEt_profile" ,meNeutralHadronEt_profile));
1015 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfChargedHadronEt_profile" ,meChargedHadronEt_profile));
1016 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFHadronEt_profile" ,meHFHadronEt_profile));
1017 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PfHFEMEt_profile" ,meHFEMEt_profile));
1021 if (fill_met_high_level_histo){
1022 hMExLS = ibooker.
book2D(
"MExLS",
"MEx_LS",200,-200,200,250,0.,2500.);
1024 hMExLS->setAxisTitle(
"Lumi Section",2);
1025 hMExLS->getTH2F()->SetOption(
"colz");
1026 hMEyLS = ibooker.
book2D(
"MEyLS",
"MEy_LS",200,-200,200,250,0.,2500.);
1028 hMEyLS->setAxisTitle(
"Lumi Section",2);
1029 hMEyLS->getTH2F()->SetOption(
"colz");
1030 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MExLS",hMExLS));
1031 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MEyLS",hMEyLS));
1035 hMETRate = ibooker.
book1D(
"METRate",
"METRate", 200, 0, 1000);
1036 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"METRate",hMETRate));
1040 lumisecME = ibooker.
book1D(
"lumisec",
"lumisec", 2501, -1., 2500.);
1041 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
"JetMET/lumisec",lumisecME));
1052 if ((techTrig->second).algoName() == m_l1algoname_) {
1053 m_bitAlgTechTrig_=(techTrig->second).algoBitNumber();
1061 if (hltConfig_.init(iRun,iSetup,triggerResultsLabel_.process(),changed)) {
1072 if (verbose_)
std::cout <<
"HLTEventAnalyzerAOD::analyze:"
1073 <<
" config extraction failure with process name "
1074 << triggerResultsLabel_.process() << std::endl;
1077 allTriggerNames_.clear();
1078 for (
unsigned int i = 0;
i<hltConfig_.size();
i++) {
1079 allTriggerNames_.push_back(hltConfig_.triggerName(
i));
1084 for ( std::vector<GenericTriggerEventFlag *>::const_iterator it = triggerFolderEventFlag_.begin(); it!= triggerFolderEventFlag_.end(); it++) {
1085 int pos = it - triggerFolderEventFlag_.begin();
1087 (*it)->initRun( iRun, iSetup );
1088 if (triggerSelectedSubFolders_[pos].exists(
std::string(
"hltDBKey"))) {
1089 if ((*it)->expressionsFromDB((*it)->hltDBKey(), iSetup)[0] !=
"CONFIG_ERROR")
1090 triggerFolderExpr_[pos] = (*it)->expressionsFromDB((*it)->hltDBKey(), iSetup);
1095 bool changed_filter=
true;
1096 std::vector<int>initializeFilter(8,-1);
1097 miniaodFilterIndex_=initializeFilter;
1098 if (FilterhltConfig_.init(iRun,iSetup,METFilterMiniAODLabel_.process(),changed_filter)){
1100 for(
unsigned int i=0;
i<FilterhltConfig_.size();
i++){
1103 std::size_t
found=search2.find(search);
1104 if(found!=std::string::npos){
1105 miniaodFilterIndex_[0]=
i;
1107 search2=CSCHaloResultTag_.label().substr(0,CSCHaloResultTag_.label().size()-3);
1108 found=search2.find(search);
1109 if(found!=std::string::npos){
1110 miniaodFilterIndex_[1]=
i;
1112 search2=eeBadScFilterTag_.label().substr(0,eeBadScFilterTag_.label().size()-3);
1113 found=search2.find(search);
1114 if(found!=std::string::npos){
1115 miniaodFilterIndex_[2]=
i;
1117 search2=HBHEIsoNoiseStringMiniAOD;
1118 found=search2.find(search);
1119 if(found!=std::string::npos){
1120 miniaodFilterIndex_[3]=
i;
1122 search2=CSCHalo2015ResultTag_.label().substr(0,CSCHalo2015ResultTag_.label().size()-3);
1123 found=search2.find(search);
1124 if(found!=std::string::npos){
1125 miniaodFilterIndex_[4]=
i;
1127 search2=EcalDeadCellTriggerTag_.label().substr(0,EcalDeadCellTriggerTag_.label().size()-3);
1128 found=search2.find(search);
1129 if(found!=std::string::npos){
1130 miniaodFilterIndex_[5]=
i;
1132 search2=EcalDeadCellBoundaryTag_.label().substr(0,EcalDeadCellBoundaryTag_.label().size()-3);
1133 found=search2.find(search);
1134 if(found!=std::string::npos){
1135 miniaodFilterIndex_[6]=
i;
1137 search2=HcalStripHaloTag_.label().substr(0,HcalStripHaloTag_.label().size()-3);
1138 found=search2.find(search);
1139 if(found!=std::string::npos){
1140 miniaodFilterIndex_[7]=
i;
1143 }
else if(FilterhltConfig_.init(iRun,iSetup,METFilterMiniAODLabel2_.process(),changed_filter)){
1145 for(
unsigned int i=0;
i<FilterhltConfig_.size();
i++){
1148 std::size_t
found=search2.find(search);
1149 if(found!=std::string::npos){
1150 miniaodFilterIndex_[0]=
i;
1152 search2=CSCHaloResultTag_.label().substr(0,CSCHaloResultTag_.label().size()-3);
1153 found=search2.find(search);
1154 if(found!=std::string::npos){
1155 miniaodFilterIndex_[1]=
i;
1157 search2=eeBadScFilterTag_.label().substr(0,eeBadScFilterTag_.label().size()-3);
1158 found=search2.find(search);
1159 if(found!=std::string::npos){
1160 miniaodFilterIndex_[2]=
i;
1162 search2=HBHEIsoNoiseStringMiniAOD;
1163 found=search2.find(search);
1164 if(found!=std::string::npos){
1165 miniaodFilterIndex_[3]=
i;
1167 search2=CSCHalo2015ResultTag_.label().substr(0,CSCHalo2015ResultTag_.label().size()-3);
1168 found=search2.find(search);
1169 if(found!=std::string::npos){
1170 miniaodFilterIndex_[4]=
i;
1172 search2=EcalDeadCellTriggerTag_.label().substr(0,EcalDeadCellTriggerTag_.label().size()-3);
1173 found=search2.find(search);
1174 if(found!=std::string::npos){
1175 miniaodFilterIndex_[5]=
i;
1177 search2=EcalDeadCellBoundaryTag_.label().substr(0,EcalDeadCellBoundaryTag_.label().size()-3);
1178 found=search2.find(search);
1179 if(found!=std::string::npos){
1180 miniaodFilterIndex_[6]=
i;
1182 search2=HcalStripHaloTag_.label().substr(0,HcalStripHaloTag_.label().size()-3);
1183 found=search2.find(search);
1184 if(found!=std::string::npos){
1185 miniaodFilterIndex_[7]=
i;
1189 std::cout<<
"nothing found with both RECO and reRECO label"<<std::endl;
1205 meLumiSec = map_dijet_MEs[
"JetMET/lumisec"];
1211 tlumisec = meLumiSec->
getTH1F();
1216 for (
int i=0;
i< (tlumisec->GetNbinsX());
i++){
1217 if (tlumisec->GetBinContent(
i)!=0){
1219 totlssecsum+=tlumisec->GetBinContent(
i);
1222 int num_per_ls=(double)totlssecsum/(
double)totlsec;
1223 totlsec=totlsec+tlumisec->GetBinContent(tlumisec->GetNbinsX()+1)/(
double)num_per_ls;
1224 totltime = double(totlsec*90);
1227 if (totltime==0.) totltime=1.;
1236 for (std::vector<std::string>::const_iterator ic = folderNames_.begin(); ic != folderNames_.end(); ic++) {
1238 DirName = dirName+*ic;
1239 makeRatePlot(DirName,totltime);
1240 for ( std::vector<GenericTriggerEventFlag *>::const_iterator it = triggerFolderEventFlag_.begin(); it!= triggerFolderEventFlag_.end(); it++) {
1241 int pos = it - triggerFolderEventFlag_.begin();
1243 makeRatePlot(DirName+
"/"+triggerFolderLabels_[pos],totltime);
1262 if ( meMET && mMETRate){
1267 tMETRate = (TH1F*) tMET->Clone(
"METRateHist");
1268 for (
int i = tMETRate->GetNbinsX()-1;
i>=0;
i--){
1269 mMETRate->
setBinContent(
i+1,tMETRate->GetBinContent(
i+2)+tMET->GetBinContent(
i+1));
1271 for (
int i = 0;
i<tMETRate->GetNbinsX();
i++){
1272 mMETRate->
setBinContent(
i+1,tMETRate->GetBinContent(
i+1)/double(totltime));
1284 int myLuminosityBlock;
1286 if(fill_met_high_level_histo){
1287 lumisecME=map_dijet_MEs[
"JetMET/lumisec"];
if(lumisecME && lumisecME->getRootObject()) lumisecME->Fill(myLuminosityBlock);
1290 if (myLuminosityBlock<LSBegin_)
return;
1291 if (myLuminosityBlock>LSEnd_ && LSEnd_>0)
return;
1293 if (verbose_)
std::cout <<
"METAnalyzer analyze" << std::endl;
1312 std::vector<int> triggerFolderDecisions;
1313 triggerFolderDecisions_ = std::vector<int> (triggerFolderEventFlag_.size(), 0);
1318 if( triggerResults.isValid()) {
1321 int ntrigs = (*triggerResults).size();
1322 if (verbose_)
std::cout <<
"ntrigs=" << ntrigs << std::endl;
1324 for (std::vector<GenericTriggerEventFlag *>::const_iterator it = triggerFolderEventFlag_.begin(); it!=triggerFolderEventFlag_.end();it++) {
1325 unsigned int pos = it - triggerFolderEventFlag_.begin();
1326 bool fd = (*it)->accept(iEvent, iSetup);
1327 triggerFolderDecisions_[pos] =
fd;
1329 allTriggerDecisions_.clear();
1330 for (
unsigned int i=0;
i<allTriggerNames_.size();++
i) {
1331 allTriggerDecisions_.push_back((*triggerResults).accept(
i));
1350 iEvent.
getByToken(caloMetToken_, calometcoll);
1351 if(!calometcoll.
isValid())
return;
1355 if(!pfmetcoll.
isValid())
return;
1359 if(!patmetcoll.
isValid())
return;
1370 met=&(pfmetcoll->front());
1371 pfmet=&(pfmetcoll->front());
1374 met=&(calometcoll->front());
1375 calomet=&(calometcoll->front());
1378 met=&(patmetcoll->front());
1379 patmet=&(patmetcoll->front());
1382 LogTrace(
"METAnalyzer")<<
"[METAnalyzer] Call to the MET analyzer";
1404 bool bJetID =
false;
1405 bool bDiJetID =
false;
1419 LogDebug(
"") <<
"METAnalyzer: Could not find calojet product" << std::endl;
1420 if (verbose_)
std::cout <<
"METAnalyzer: Could not find calojet product" << std::endl;
1422 collsize=caloJets->size();
1437 iEvent.
getByToken(jetID_ValueMapToken_,jetID_ValueMap_Handle);
1441 if (isMiniAODMet_){ iEvent.
getByToken(patJetsToken_, patJets);
1443 LogDebug(
"") <<
"METAnalyzer: Could not find patjet product" << std::endl;
1444 if (verbose_)
std::cout <<
"METAnalyzer: Could not find patjet product" << std::endl;
1446 collsize=patJets->size();
1449 if (isPFMet_){ iEvent.
getByToken(pfJetsToken_, pfJets);
1451 LogDebug(
"") <<
"METAnalyzer: Could not find pfjet product" << std::endl;
1452 if (verbose_)
std::cout <<
"METAnalyzer: Could not find pfjet product" << std::endl;
1454 collsize=pfJets->size();
1457 unsigned int ind1=-1;
1459 bool pass_jetID1=
false;
1460 unsigned int ind2=-1;
1462 bool pass_jetID2=
false;
1465 bool pass_correction_flag=
false;
1467 iEvent.
getByToken(jetCorrectorToken_, jetCorr);
1469 pass_correction_flag=
true;
1472 pass_correction_flag=
true;
1475 for (
int ijet=0; ijet<collsize; ijet++) {
1478 bool iscleaned=
false;
1479 if (pass_correction_flag) {
1481 scale = jetCorr->correction((*caloJets)[ijet]);
1487 scale = jetCorr->correction((*pfJets)[ijet]);
1491 pt_jet=scale*(*caloJets)[ijet].pt();
1492 if(pt_jet> ptThreshold_){
1495 reco::JetID jetID = (*jetID_ValueMap_Handle)[calojetref];
1496 iscleaned = jetIDFunctorLoose((*caloJets)[ijet], jetID);
1518 pt_jet=scale*(*pfJets)[ijet].pt();
1519 if(pt_jet> ptThreshold_){
1520 iscleaned = pfjetIDFunctorLoose((*pfJets)[ijet]);
1524 pt_jet=(*patJets)[ijet].pt();
1525 if(pt_jet> ptThreshold_){
1526 pat::strbitset stringbitset=pfjetIDFunctorLoose.getBitTemplate();
1527 iscleaned = pfjetIDFunctorLoose((*patJets)[ijet],stringbitset);
1536 pass_jetID2=pass_jetID1;
1539 pass_jetID1=iscleaned;
1540 }
else if (pt_jet>pt2){
1543 pass_jetID2=iscleaned;
1546 if(pass_jetID1 && pass_jetID2){
1549 dphi=fabs((*caloJets)[ind1].
phi()-(*caloJets)[ind2].
phi());
1555 dphi=fabs((*pfJets)[ind1].
phi()-(*pfJets)[ind2].
phi());
1558 dphi=fabs((*patJets)[0].
phi()-(*patJets)[1].
phi());
1561 dphi=2*acos(-1.)-dphi;
1571 iEvent.
getByToken(vertexToken_, vertexHandle);
1573 if (!vertexHandle.
isValid()) {
1574 LogDebug(
"") <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1575 if (verbose_)
std::cout <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1578 if ( vertexHandle.
isValid() ){
1580 numPV_ = vertexCollection.size();
1582 bool bPrimaryVertex = (bypassAllPVChecks_ || (numPV_>0));
1584 bool bZJets =
false;
1596 if(Muons.
isValid() && Muons->size()>1){
1597 for (
unsigned int i=0;
i<Muons->size();
i++){
1598 bool pass_muon_id=
false;
1599 bool pass_muon_iso=
false;
1600 double dxy=fabs((*Muons)[
i].muonBestTrack()->dxy());
1601 double dz=fabs((*Muons)[
i].muonBestTrack()->dz());
1603 dxy=fabs((*Muons)[
i].muonBestTrack()->dxy((*vertexHandle)[0].
position()));
1604 dz=fabs((*Muons)[
i].muonBestTrack()->dz((*vertexHandle)[0].
position()));
1606 if((*Muons)[
i].
pt()>20 && fabs((*Muons)[
i].
eta())<2.3){
1607 if((*Muons)[
i].isGlobalMuon() && (*Muons)[
i].isPFMuon() && (*Muons)[
i].globalTrack()->hitPattern().numberOfValidMuonHits() > 0 && (*Muons)[
i].numberOfMatchedStations() > 1 && dxy < 0.2 && (*Muons)[
i].numberOfMatchedStations() > 1 && dz<0.5 && (*Muons)[
i].
innerTrack()->hitPattern().numberOfValidPixelHits() > 0 && (*Muons)[
i].
innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5){
1611 float muonIsoPFdb = ((*Muons)[
i].pfIsolationR04().sumChargedHadronPt
1612 +
std::max(0., (*Muons)[
i].pfIsolationR04().sumNeutralHadronEt + (*Muons)[
i].pfIsolationR04().sumPhotonEt - 0.5*(*Muons)[
i].pfIsolationR04().sumPUPt))/(*Muons)[
i].pt();
1613 if(muonIsoPFdb<0.12){
1617 if(pass_muon_id && pass_muon_iso){
1618 if((*Muons)[
i].
pt()>pt_muon0){
1619 mu_index1=mu_index0;
1622 pt_muon0=(*Muons)[
i].pt();
1623 }
else if ((*Muons)[
i].
pt()>pt_muon1){
1625 pt_muon1=(*Muons)[
i].pt();
1630 if(mu_index0>=0 && mu_index1>=0){
1631 if((*Muons)[mu_index0].
charge()*(*Muons)[mu_index1].
charge()<0){
1632 zCand=(*Muons)[mu_index0].polarP4()+(*Muons)[mu_index1].polarP4();
1633 if(fabs(zCand.M()-91.)<20){
1643 iEvent.
getByToken( gtToken_, gtReadoutRecord);
1645 bool techTriggerResultBxM2 =
false;
1646 bool techTriggerResultBxM1 =
false;
1647 bool techTriggerResultBx0 =
false;
1649 if (!gtReadoutRecord.
isValid()) {
1650 LogDebug(
"") <<
"METAnalyzer: Could not find GT readout record" << std::endl;
1651 if (verbose_)
std::cout <<
"METAnalyzer: Could not find GT readout record product" << std::endl;
1654 const TechnicalTriggerWord& technicalTriggerWordBeforeMaskBxM2 = gtReadoutRecord->technicalTriggerWord(-2);
1655 const TechnicalTriggerWord& technicalTriggerWordBeforeMaskBxM1 = gtReadoutRecord->technicalTriggerWord(-1);
1656 const TechnicalTriggerWord& technicalTriggerWordBeforeMaskBx0 = gtReadoutRecord->technicalTriggerWord();
1659 if (m_bitAlgTechTrig_ > -1 && technicalTriggerWordBeforeMaskBx0.size() > 0) {
1660 techTriggerResultBx0 = technicalTriggerWordBeforeMaskBx0.at(m_bitAlgTechTrig_);
1661 if(techTriggerResultBx0!=0){
1662 techTriggerResultBxM2 = technicalTriggerWordBeforeMaskBxM2.at(m_bitAlgTechTrig_);
1663 techTriggerResultBxM1 = technicalTriggerWordBeforeMaskBxM1.at(m_bitAlgTechTrig_);
1668 std::vector<bool>trigger_flag(4,
false);
1669 if(techTriggerResultBx0 && techTriggerResultBxM2 && techTriggerResultBxM1){
1670 trigger_flag[0]=
true;
1672 if(techTriggerResultBx0 && techTriggerResultBxM1){
1673 trigger_flag[1]=
true; }
1674 if(techTriggerResultBx0 && !techTriggerResultBxM1){
1675 trigger_flag[2]=
true;
1677 if(techTriggerResultBx0 && !techTriggerResultBxM2 && !techTriggerResultBxM1){
1678 trigger_flag[3]=
true;
1680 std::vector<bool>filter_decisions(8,
false);
1681 if(!isMiniAODMet_ && !runcosmics_){
1683 iEvent.
getByToken(hbheNoiseFilterResultToken_, HBHENoiseFilterResultHandle);
1684 if (!HBHENoiseFilterResultHandle.
isValid()) {
1685 LogDebug(
"") <<
"METAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
1686 if (verbose_)
std::cout <<
"METAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
1688 filter_decisions[0]= *HBHENoiseFilterResultHandle;
1690 iEvent.
getByToken(CSCHaloResultToken_, CSCTightHaloFilterResultHandle);
1691 if (!CSCTightHaloFilterResultHandle.
isValid()) {
1692 LogDebug(
"") <<
"METAnalyzer: Could not find CSCTightHaloFilterResultHandle" << std::endl;
1693 if (verbose_)
std::cout <<
"METAnalyzer: CSCTightHaloFilterResultHandle" << std::endl;
1695 filter_decisions[1]= *CSCTightHaloFilterResultHandle;
1697 iEvent.
getByToken(eeBadScFilterToken_, eeBadScFilterResultHandle);
1698 if (!eeBadScFilterResultHandle.
isValid()) {
1699 LogDebug(
"") <<
"METAnalyzer: Could not find eeBadScFilterResultHandle" << std::endl;
1700 if (verbose_)
std::cout <<
"METAnalyzer: eeBadScFilterResultHandle" << std::endl;
1702 filter_decisions[2]= *eeBadScFilterResultHandle;
1704 iEvent.
getByToken(hbheIsoNoiseFilterResultToken_, HBHENoiseIsoFilterResultHandle);
1705 if (!HBHENoiseIsoFilterResultHandle.
isValid()) {
1706 LogDebug(
"") <<
"METAnalyzer: Could not find HBHENoiseIsoFilterResult" << std::endl;
1707 if (verbose_)
std::cout <<
"METAnalyzer: Could not find HBHENoiseIsoFilterResult" << std::endl;
1709 filter_decisions[3]= *HBHENoiseIsoFilterResultHandle;
1711 iEvent.
getByToken(CSCHalo2015ResultToken_, CSCTightHalo2015FilterResultHandle);
1712 if (!CSCTightHalo2015FilterResultHandle.
isValid()) {
1713 LogDebug(
"") <<
"METAnalyzer: Could not find CSCTightHalo2015FilterResultHandle" << std::endl;
1714 if (verbose_)
std::cout <<
"METAnalyzer: CSCTightHalo2015FilterResultHandle" << std::endl;
1716 filter_decisions[4]= *CSCTightHalo2015FilterResultHandle;
1718 iEvent.
getByToken(EcalDeadCellTriggerToken_, EcalDeadCellTriggerFilterResultHandle);
1719 if (!EcalDeadCellTriggerFilterResultHandle.
isValid()) {
1720 LogDebug(
"") <<
"METAnalyzer: Could not find EcalDeadCellTriggerFilterResultHandle" << std::endl;
1721 if (verbose_)
std::cout <<
"METAnalyzer: EcalDeadCellTriggerFilterResultHandle" << std::endl;
1723 filter_decisions[5]= *EcalDeadCellTriggerFilterResultHandle;
1725 iEvent.
getByToken(EcalDeadCellBoundaryToken_, EcalDeadCellBoundaryHandle);
1726 if (!EcalDeadCellBoundaryHandle.
isValid()) {
1727 LogDebug(
"") <<
"METAnalyzer: Could not find EcalDeadCellBoundaryHandle" << std::endl;
1728 if (verbose_)
std::cout <<
"METAnalyzer: EcalDeadCellBoundaryHandle" << std::endl;
1730 filter_decisions[6]= *EcalDeadCellBoundaryHandle;
1732 iEvent.
getByToken(HcalStripHaloToken_, HcalStripHaloFilterHandle);
1733 if (!HcalStripHaloFilterHandle.
isValid()) {
1734 LogDebug(
"") <<
"METAnalyzer: Could not find CSCTightHalo2015FilterResultHandle" << std::endl;
1735 if (verbose_)
std::cout <<
"METAnalyzer: CSCTightHalo2015FilterResultHandle" << std::endl;
1737 filter_decisions[7]= *HcalStripHaloFilterHandle;
1740 iEvent.
getByToken(METFilterMiniAODToken_, metFilterResults);
1741 if(metFilterResults.
isValid()){
1742 if(miniaodFilterIndex_[0]!=-1){
1743 filter_decisions[0] = metFilterResults->accept(miniaodFilterIndex_[0]);
1745 if(miniaodFilterIndex_[1]!=-1){
1746 filter_decisions[1] = metFilterResults->accept(miniaodFilterIndex_[1]);
1748 if(miniaodFilterIndex_[2]!=-1){
1749 filter_decisions[2] = metFilterResults->accept(miniaodFilterIndex_[2]);
1751 if(miniaodFilterIndex_[3]!=-1){
1752 filter_decisions[3] = metFilterResults->accept(miniaodFilterIndex_[3]);
1754 if(miniaodFilterIndex_[4]!=-1){
1755 filter_decisions[4] = metFilterResults->accept(miniaodFilterIndex_[4]);
1757 if(miniaodFilterIndex_[5]!=-1){
1758 filter_decisions[5] = metFilterResults->accept(miniaodFilterIndex_[5]);
1760 if(miniaodFilterIndex_[6]!=-1){
1761 filter_decisions[6] = metFilterResults->accept(miniaodFilterIndex_[6]);
1763 if(miniaodFilterIndex_[7]!=-1){
1764 filter_decisions[7] = metFilterResults->accept(miniaodFilterIndex_[7]);
1767 iEvent.
getByToken(METFilterMiniAODToken2_, metFilterResults);
1768 if(metFilterResults.
isValid()){
1769 if(miniaodFilterIndex_[0]!=-1){
1770 filter_decisions[0] = metFilterResults->accept(miniaodFilterIndex_[0]);
1772 if(miniaodFilterIndex_[1]!=-1){
1773 filter_decisions[1] = metFilterResults->accept(miniaodFilterIndex_[1]);
1775 if(miniaodFilterIndex_[2]!=-1){
1776 filter_decisions[2] = metFilterResults->accept(miniaodFilterIndex_[2]);
1778 if(miniaodFilterIndex_[3]!=-1){
1779 filter_decisions[3] = metFilterResults->accept(miniaodFilterIndex_[3]);
1781 if(miniaodFilterIndex_[4]!=-1){
1782 filter_decisions[4] = metFilterResults->accept(miniaodFilterIndex_[4]);
1784 if(miniaodFilterIndex_[5]!=-1){
1785 filter_decisions[5] = metFilterResults->accept(miniaodFilterIndex_[5]);
1787 if(miniaodFilterIndex_[6]!=-1){
1788 filter_decisions[6] = metFilterResults->accept(miniaodFilterIndex_[6]);
1790 if(miniaodFilterIndex_[7]!=-1){
1791 filter_decisions[7] = metFilterResults->accept(miniaodFilterIndex_[7]);
1796 bool HBHENoiseFilterResultFlag=filter_decisions[0];
1799 bool bHBHENoiseFilter = HBHENoiseFilterResultFlag;
1802 bool bDCSFilter = (bypassAllDCSChecks_ || DCSFilter_->filter(iEvent, iSetup));
1806 for (std::vector<std::string>::const_iterator ic = folderNames_.begin();
1807 ic != folderNames_.end(); ic++){
1808 bool pass_selection =
false;
1809 if ((*ic==
"Uncleaned") &&(isCaloMet_ || bPrimaryVertex)){
1810 fillMESet(iEvent, DirName_old+
"/"+*ic, *met,*patmet, *pfmet,*calomet,zCand, map_dijet_MEs,trigger_flag,filter_decisions);
1811 pass_selection =
true;
1814 if ((*ic==
"Cleaned") && bDCSFilter && bHBHENoiseFilter && bPrimaryVertex && bJetID){
1815 fillMESet(iEvent, DirName_old+
"/"+*ic, *met,*patmet,*pfmet,*calomet,zCand,map_dijet_MEs,trigger_flag,filter_decisions);
1816 pass_selection=
true;
1818 if ((*ic==
"DiJet" ) &&bDCSFilter&&bHBHENoiseFilter&& bPrimaryVertex&& bDiJetID){
1819 fillMESet(iEvent, DirName_old+
"/"+*ic, *met,*patmet,*pfmet,*calomet,zCand,map_dijet_MEs,trigger_flag,filter_decisions);
1820 pass_selection=
true;
1822 if ((*ic==
"ZJets" ) &&bDCSFilter&&bHBHENoiseFilter&& bPrimaryVertex&& bZJets){
1823 fillMESet(iEvent, DirName_old+
"/"+*ic, *met,*patmet,*pfmet,*calomet,zCand,map_dijet_MEs,trigger_flag,filter_decisions);
1824 pass_selection=
true;
1826 if(pass_selection && isPFMet_){
1827 DirName=DirName_old+
"/"+*ic;
1838 bool bLumiSecPlot=fill_met_high_level_histo;
1839 bool fillPFCandidatePlots=
false;
1840 if (DirName.find(
"Cleaned")!=std::string::npos) {
1841 fillPFCandidatePlots=
true;
1842 fillMonitorElement(iEvent, DirName,
std::string(
""), met, patmet, pfmet, calomet, zCand, map_of_MEs,bLumiSecPlot,fillPFCandidatePlots,techTriggerCase,METFilterDecision);
1843 for (
unsigned int i = 0;
i<triggerFolderLabels_.size();
i++) {
1844 fillPFCandidatePlots=
false;
1845 if (triggerFolderDecisions_[
i]){
1846 fillMonitorElement(iEvent, DirName, triggerFolderLabels_[i], met, patmet, pfmet, calomet, zCand, map_of_MEs,bLumiSecPlot,fillPFCandidatePlots,techTriggerCase,METFilterDecision);
1849 }
else if (DirName.find(
"DiJet")!=std::string::npos) {
1850 fillMonitorElement(iEvent, DirName,
std::string(
""), met, patmet, pfmet, calomet, zCand, map_of_MEs,bLumiSecPlot,fillPFCandidatePlots,techTriggerCase,METFilterDecision);
1851 for (
unsigned int i = 0;
i<triggerFolderLabels_.size();
i++) {
1852 if (triggerFolderDecisions_[
i]) fillMonitorElement(iEvent, DirName, triggerFolderLabels_[i], met, patmet, pfmet, calomet, zCand, map_of_MEs,bLumiSecPlot,fillPFCandidatePlots,techTriggerCase,METFilterDecision);
1854 }
else if (DirName.find(
"ZJets")!=std::string::npos) {
1855 fillMonitorElement(iEvent, DirName,
std::string(
""), met, patmet, pfmet, calomet, zCand, map_of_MEs,bLumiSecPlot,fillPFCandidatePlots,techTriggerCase,METFilterDecision);
1857 fillMonitorElement(iEvent, DirName,
std::string(
""), met, patmet, pfmet, calomet, zCand, map_of_MEs,bLumiSecPlot,fillPFCandidatePlots,techTriggerCase,METFilterDecision);
1866 const reco::MET&
met,
const pat::MET & patmet,
const reco::PFMET & pfmet,
const reco::CaloMET &calomet,
const reco::Candidate::PolarLorentzVector& zCand,std::map<std::string,MonitorElement*>& map_of_MEs,
bool bLumiSecPlot,
bool fillPFCandidatePlots,std::vector<bool> techTriggerCase,std::vector<bool> METFilterDecision)
1868 bool do_only_Z_histograms=
false;
1869 if(DirName.find(
"ZJets")!=std::string::npos) {
1870 do_only_Z_histograms=
true;
1873 double u_x=-met.
px()-zCand.Px();
1874 double u_y=-met.
py()-zCand.Py();
1878 double u_perp=
sqrt(u_x*u_x + u_y*u_y);
1882 e_Z_x = zCand.Px()/zCand.Pt();
1883 e_Z_y = zCand.Py()/zCand.Pt();
1885 u_par=u_x*e_Z_x+u_y*e_Z_y;
1886 u_perp= -e_Z_y * u_x + e_Z_x*u_y;
1888 meZJets_u_par = map_of_MEs[DirName+
"/"+
"u_parallel_Z_inc"];
if ( meZJets_u_par && meZJets_u_par->getRootObject()) meZJets_u_par ->
Fill(u_par);
1890 meZJets_u_par_ZPt_0_15 = map_of_MEs[DirName+
"/"+
"u_parallel_ZPt_0_15"];
if ( meZJets_u_par_ZPt_0_15 && meZJets_u_par_ZPt_0_15->getRootObject()) meZJets_u_par_ZPt_0_15 ->
Fill(u_par);
1891 }
else if(zCand.Pt()<30){
1892 meZJets_u_par_ZPt_15_30 = map_of_MEs[DirName+
"/"+
"u_parallel_ZPt_15_30"];
if ( meZJets_u_par_ZPt_15_30 && meZJets_u_par_ZPt_15_30->getRootObject()) meZJets_u_par_ZPt_15_30 ->
Fill(u_par);
1893 }
else if(zCand.Pt()<55){
1894 meZJets_u_par_ZPt_30_55 = map_of_MEs[DirName+
"/"+
"u_parallel_ZPt_30_55"];
if ( meZJets_u_par_ZPt_30_55 && meZJets_u_par_ZPt_30_55->getRootObject()) meZJets_u_par_ZPt_30_55 ->
Fill(u_par);
1895 }
else if(zCand.Pt()<75){
1896 meZJets_u_par_ZPt_55_75 = map_of_MEs[DirName+
"/"+
"u_parallel_ZPt_55_75"];
if ( meZJets_u_par_ZPt_55_75 && meZJets_u_par_ZPt_55_75->getRootObject()) meZJets_u_par_ZPt_55_75 ->
Fill(u_par);
1897 }
else if(zCand.Pt()<150){
1898 meZJets_u_par_ZPt_75_150 = map_of_MEs[DirName+
"/"+
"u_parallel_ZPt_75_150"];
if ( meZJets_u_par_ZPt_75_150 && meZJets_u_par_ZPt_75_150->getRootObject()) meZJets_u_par_ZPt_75_150 ->
Fill(u_par);
1899 }
else if(zCand.Pt()<290){
1900 meZJets_u_par_ZPt_150_290 = map_of_MEs[DirName+
"/"+
"u_parallel_ZPt_150_290"];
if ( meZJets_u_par_ZPt_150_290 && meZJets_u_par_ZPt_150_290->getRootObject()) meZJets_u_par_ZPt_150_290 ->
Fill(u_par);
1902 meZJets_u_par_ZPt_290 = map_of_MEs[DirName+
"/"+
"u_parallel_ZPt_290"];
if ( meZJets_u_par_ZPt_290 && meZJets_u_par_ZPt_290->getRootObject()) meZJets_u_par_ZPt_290 ->
Fill(u_par);
1905 meZJets_u_perp = map_of_MEs[DirName+
"/"+
"u_perp_Z_inc"];
if ( meZJets_u_perp && meZJets_u_perp->getRootObject()) meZJets_u_perp ->
Fill(u_perp);
1907 meZJets_u_perp_ZPt_0_15 = map_of_MEs[DirName+
"/"+
"u_perp_ZPt_0_15"];
if ( meZJets_u_perp_ZPt_0_15 && meZJets_u_perp_ZPt_0_15->getRootObject()) meZJets_u_perp_ZPt_0_15 ->
Fill(u_perp);
1908 }
else if(zCand.Pt()<30){
1909 meZJets_u_perp_ZPt_15_30 = map_of_MEs[DirName+
"/"+
"u_perp_ZPt_15_30"];
if ( meZJets_u_perp_ZPt_15_30 && meZJets_u_perp_ZPt_15_30->getRootObject()) meZJets_u_perp_ZPt_15_30 ->
Fill(u_perp);
1910 }
else if(zCand.Pt()<55){
1911 meZJets_u_perp_ZPt_30_55 = map_of_MEs[DirName+
"/"+
"u_perp_ZPt_30_55"];
if ( meZJets_u_perp_ZPt_30_55 && meZJets_u_perp_ZPt_30_55->getRootObject()) meZJets_u_perp_ZPt_30_55 ->
Fill(u_perp);
1912 }
else if(zCand.Pt()<75){
1913 meZJets_u_perp_ZPt_55_75 = map_of_MEs[DirName+
"/"+
"u_perp_ZPt_55_75"];
if ( meZJets_u_perp_ZPt_55_75 && meZJets_u_perp_ZPt_55_75->getRootObject()) meZJets_u_perp_ZPt_55_75 ->
Fill(u_perp);
1914 }
else if(zCand.Pt()<150){
1915 meZJets_u_perp_ZPt_75_150 = map_of_MEs[DirName+
"/"+
"u_perp_ZPt_75_150"];
if ( meZJets_u_perp_ZPt_75_150 && meZJets_u_perp_ZPt_75_150->getRootObject()) meZJets_u_perp_ZPt_75_150 ->
Fill(u_perp);
1916 }
else if(zCand.Pt()<290){
1917 meZJets_u_perp_ZPt_150_290 = map_of_MEs[DirName+
"/"+
"u_perp_ZPt_150_290"];
if ( meZJets_u_perp_ZPt_150_290 && meZJets_u_perp_ZPt_150_290->getRootObject()) meZJets_u_perp_ZPt_150_290 ->
Fill(u_perp);
1919 meZJets_u_perp_ZPt_290 = map_of_MEs[DirName+
"/"+
"u_perp_ZPt_290"];
if ( meZJets_u_perp_ZPt_290 && meZJets_u_perp_ZPt_290->getRootObject()) meZJets_u_perp_ZPt_290 ->
Fill(u_perp);
1924 if(!do_only_Z_histograms){
1926 double SumET = met.
sumEt();
1927 double METSig = met.
mEtSig();
1929 double MET = met.
pt();
1930 double MEx = met.
px();
1931 double MEy = met.
py();
1932 double METPhi = met.
phi();
1934 int myLuminosityBlock;
1938 if (subFolderName!=
""){
1939 DirName = DirName +
"/"+subFolderName;
1942 hTrigger = map_of_MEs[DirName+
"/triggerResults"];
1943 if (hTrigger && hTrigger->getRootObject()) {
1944 for (
unsigned int i = 0;
i<allTriggerDecisions_.size();
i++){
1945 if(
i<(
unsigned int)hTrigger->getNbinsX()){
1946 hTrigger->Fill(
i + .5, allTriggerDecisions_[
i]);
1952 hMEx = map_of_MEs[DirName+
"/"+
"MEx"];
if (hMEx && hMEx->getRootObject()) hMEx ->
Fill(MEx);
1953 hMEy = map_of_MEs[DirName+
"/"+
"MEy"];
if (hMEy && hMEy->getRootObject()) hMEy ->
Fill(MEy);
1954 hMET = map_of_MEs[DirName+
"/"+
"MET"];
if (hMET && hMET->getRootObject()) hMET ->
Fill(MET);
1955 hMET_2 = map_of_MEs[DirName+
"/"+
"MET_2"];
if (hMET_2 && hMET_2->getRootObject()) hMET_2 ->
Fill(MET);
1967 bool HBHENoiseFilterResult=
false;
1968 bool CSCTightHaloFilterResult=
false;
1969 bool eeBadScFilterResult=
false;
1970 bool HBHEIsoNoiseFilterResult=
false;
1971 bool CSCTightHalo2015FilterResult=
false;
1972 bool EcalDeadCellTriggerFilterResult=
false;
1973 bool EcalDeadCellBoundaryFilterResult=
false;
1974 bool HcalStripHaloFilterResult=
false;
1975 HBHENoiseFilterResult = METFilterDecision[0];
1976 if(HBHENoiseFilterResult){
1977 hMET_HBHENoiseFilter = map_of_MEs[DirName+
"/"+
"MET_HBHENoiseFilter"];
if (hMET_HBHENoiseFilter && hMET_HBHENoiseFilter->getRootObject()) hMET_HBHENoiseFilter ->
Fill(MET);
1979 CSCTightHaloFilterResult = METFilterDecision[1];
1980 if(CSCTightHaloFilterResult){
1981 hMET_CSCTightHaloFilter = map_of_MEs[DirName+
"/"+
"MET_CSCTightHaloFilter"];
if (hMET_CSCTightHaloFilter && hMET_CSCTightHaloFilter->getRootObject()) hMET_CSCTightHaloFilter ->
Fill(MET);
1983 eeBadScFilterResult = METFilterDecision[2];
1984 if(eeBadScFilterResult){
1985 hMET_eeBadScFilter = map_of_MEs[DirName+
"/"+
"MET_eeBadScFilter"];
if (hMET_eeBadScFilter && hMET_eeBadScFilter->getRootObject()) hMET_eeBadScFilter ->
Fill(MET);
1987 HBHEIsoNoiseFilterResult = METFilterDecision[3];
1988 if(HBHEIsoNoiseFilterResult){
1989 hMET_HBHEIsoNoiseFilter = map_of_MEs[DirName+
"/"+
"MET_HBHEIsoNoiseFilter"];
if (hMET_HBHEIsoNoiseFilter && hMET_HBHEIsoNoiseFilter->getRootObject()) hMET_HBHEIsoNoiseFilter ->
Fill(MET);
1991 CSCTightHalo2015FilterResult = METFilterDecision[4];
1992 if(CSCTightHalo2015FilterResult){
1993 hMET_CSCTightHalo2015Filter = map_of_MEs[DirName+
"/"+
"MET_CSCTightHalo2015Filter"];
if (hMET_CSCTightHalo2015Filter && hMET_CSCTightHalo2015Filter->getRootObject()) hMET_CSCTightHalo2015Filter ->
Fill(MET);
1995 EcalDeadCellTriggerFilterResult = METFilterDecision[5];
1996 if(EcalDeadCellTriggerFilterResult){
1997 hMET_EcalDeadCellTriggerFilter = map_of_MEs[DirName+
"/"+
"MET_EcalDeadCellTriggerFilter"];
if (hMET_EcalDeadCellTriggerFilter && hMET_EcalDeadCellTriggerFilter->getRootObject()) hMET_EcalDeadCellTriggerFilter ->
Fill(MET);
1999 EcalDeadCellBoundaryFilterResult = METFilterDecision[6];
2000 if(EcalDeadCellBoundaryFilterResult){
2001 hMET_EcalDeadCellBoundaryFilter = map_of_MEs[DirName+
"/"+
"MET_EcalDeadCellBoundaryFilter"];
if (hMET_EcalDeadCellBoundaryFilter && hMET_EcalDeadCellBoundaryFilter->getRootObject()) hMET_EcalDeadCellBoundaryFilter ->
Fill(MET);
2003 HcalStripHaloFilterResult = METFilterDecision[7];
2004 if(HcalStripHaloFilterResult){
2005 hMET_HcalStripHaloFilter = map_of_MEs[DirName+
"/"+
"MET_HcalStripHaloFilter"];
if (hMET_HcalStripHaloFilter && hMET_HcalStripHaloFilter->getRootObject()) hMET_HcalStripHaloFilter ->
Fill(MET);
2007 hMETPhi = map_of_MEs[DirName+
"/"+
"METPhi"];
if (hMETPhi && hMETPhi->getRootObject()) hMETPhi ->
Fill(METPhi);
2008 hSumET = map_of_MEs[DirName+
"/"+
"SumET"];
if (hSumET && hSumET->getRootObject()) hSumET ->
Fill(SumET);
2009 hMETSig = map_of_MEs[DirName+
"/"+
"METSig"];
if (hMETSig && hMETSig->getRootObject()) hMETSig ->
Fill(METSig);
2010 hMET_logx = map_of_MEs[DirName+
"/"+
"MET_logx"];
if (hMET_logx && hMET_logx->getRootObject()) hMET_logx->Fill(log10(MET));
2011 hSumET_logx = map_of_MEs[DirName+
"/"+
"SumET_logx"];
if (hSumET_logx && hSumET_logx->getRootObject()) hSumET_logx->Fill(log10(SumET));
2015 meMEx_profile = map_of_MEs[DirName +
"/MEx_profile"];
2016 meMEy_profile = map_of_MEs[DirName +
"/MEy_profile"];
2017 meMET_profile = map_of_MEs[DirName +
"/MET_profile"];
2018 meSumET_profile = map_of_MEs[DirName +
"/SumET_profile"];
2020 if (meMEx_profile && meMEx_profile ->getRootObject()) meMEx_profile ->Fill(numPV_, MEx);
2021 if (meMEy_profile && meMEy_profile ->getRootObject()) meMEy_profile ->Fill(numPV_, MEy);
2022 if (meMET_profile && meMET_profile ->getRootObject()) meMET_profile ->Fill(numPV_, MET);
2023 if (meSumET_profile && meSumET_profile->getRootObject()) meSumET_profile->Fill(numPV_, SumET);
2033 double caloHadEtInHB = calomet.
hadEtInHB();
2034 double caloHadEtInHO = calomet.
hadEtInHO();
2035 double caloHadEtInHE = calomet.
hadEtInHE();
2036 double caloHadEtInHF = calomet.
hadEtInHF();
2037 double caloEmEtInEB = calomet.
emEtInEB();
2038 double caloEmEtInEE = calomet.
emEtInEE();
2039 double caloEmEtInHF = calomet.
emEtInHF();
2041 hCaloHadEtInHB = map_of_MEs[DirName+
"/"+
"CaloHadEtInHB"];
if (hCaloHadEtInHB && hCaloHadEtInHB->getRootObject()) hCaloHadEtInHB->Fill(caloHadEtInHB);
2042 hCaloHadEtInHO = map_of_MEs[DirName+
"/"+
"CaloHadEtInHO"];
if (hCaloHadEtInHO && hCaloHadEtInHO->getRootObject()) hCaloHadEtInHO->Fill(caloHadEtInHO);
2043 hCaloHadEtInHE = map_of_MEs[DirName+
"/"+
"CaloHadEtInHE"];
if (hCaloHadEtInHE && hCaloHadEtInHE->getRootObject()) hCaloHadEtInHE->Fill(caloHadEtInHE);
2044 hCaloHadEtInHF = map_of_MEs[DirName+
"/"+
"CaloHadEtInHF"];
if (hCaloHadEtInHF && hCaloHadEtInHF->getRootObject()) hCaloHadEtInHF->Fill(caloHadEtInHF);
2045 hCaloEmEtInEB = map_of_MEs[DirName+
"/"+
"CaloEmEtInEB"];
if (hCaloEmEtInEB && hCaloEmEtInEB->getRootObject()) hCaloEmEtInEB->Fill(caloEmEtInEB);
2046 hCaloEmEtInEE = map_of_MEs[DirName+
"/"+
"CaloEmEtInEE"];
if (hCaloEmEtInEE && hCaloEmEtInEE->getRootObject()) hCaloEmEtInEE->Fill(caloEmEtInEE);
2047 hCaloEmEtInHF = map_of_MEs[DirName+
"/"+
"CaloEmEtInHF"];
if (hCaloEmEtInHF && hCaloEmEtInHF->getRootObject()) hCaloEmEtInHF->Fill(caloEmEtInHF);
2049 hCaloMETPhi020 = map_of_MEs[DirName+
"/"+
"CaloMETPhi020"];
if (MET> 20. && hCaloMETPhi020 && hCaloMETPhi020->getRootObject()) { hCaloMETPhi020->Fill(METPhi);}
2052 hCaloEtFractionHadronic = map_of_MEs[DirName+
"/"+
"CaloEtFractionHadronic"];
if (hCaloEtFractionHadronic && hCaloEtFractionHadronic->getRootObject()) hCaloEtFractionHadronic->Fill(caloEtFractionHadronic);
2053 hCaloEmEtFraction = map_of_MEs[DirName+
"/"+
"CaloEmEtFraction"];
if (hCaloEmEtFraction && hCaloEmEtFraction->getRootObject()) hCaloEmEtFraction->Fill(caloEmEtFraction);
2054 hCaloEmEtFraction020 = map_of_MEs[DirName+
"/"+
"CaloEmEtFraction020"];
if (MET> 20. && hCaloEmEtFraction020 && hCaloEmEtFraction020->getRootObject()) hCaloEmEtFraction020->Fill(caloEmEtFraction);
2080 if(fillPFCandidatePlots && fillCandidateMap_histos){
2083 iEvent.
getByToken(pflowPackedToken_, packedParticleFlow);
2085 for (
unsigned int i = 0;
i < packedParticleFlow->size();
i++) {
2087 for (
unsigned int j=0;
j<typePFCand_.size();
j++) {
2088 if (fabs(c.
pdgId())==typePFCand_[
j]) {
2090 if ( ((c.
eta()>etaMinPFCand_[
j]) && (c.
eta()<etaMaxPFCand_[
j])) || ((c.
eta()> (-etaMaxPFCand_[
j])) && (c.
eta()< (-etaMinPFCand_[
j]))) ){
2091 ptPFCand_[
j] = map_of_MEs[DirName +
"/"+ptPFCand_name_[
j]];
2092 if ( ptPFCand_[j] && ptPFCand_[j]->getRootObject()) ptPFCand_[j]->
Fill(c.
eta(), c.
phi(), c.
pt()*c.
puppiWeight());
2093 occupancyPFCand_[
j] = map_of_MEs[DirName +
"/"+occupancyPFCand_name_[
j]];
2094 if ( occupancyPFCand_[j] && occupancyPFCand_[j]->getRootObject()) occupancyPFCand_[j]->
Fill(c.
eta(), c.
phi()),c.
puppiWeight();
2095 ptPFCand_puppiNolepWeight_[
j] = map_of_MEs[DirName +
"/"+ptPFCand_name_puppiNolepWeight_[
j]];
2096 if ( ptPFCand_puppiNolepWeight_[j] && ptPFCand_puppiNolepWeight_[j]->getRootObject()) ptPFCand_puppiNolepWeight_[
j]->Fill(c.
eta(), c.
phi(), c.
pt()*c.
puppiWeightNoLep());
2097 occupancyPFCand_puppiNolepWeight_[
j] = map_of_MEs[DirName +
"/"+occupancyPFCand_name_puppiNolepWeight_[
j]];
2098 if ( occupancyPFCand_puppiNolepWeight_[j] && occupancyPFCand_puppiNolepWeight_[j]->getRootObject()) occupancyPFCand_puppiNolepWeight_[j]->
Fill(c.
eta(), c.
phi()),c.
puppiWeightNoLep();
2107 if(fillPFCandidatePlots && fillCandidateMap_histos){
2108 for (
unsigned int i=0;
i<countsPFCand_.size();
i++) {
2116 iEvent.
getByToken(pflowToken_, particleFlow);
2118 float pt_sum_CHF_Barrel=0;
2119 float pt_sum_CHF_Endcap_plus=0;
2120 float pt_sum_CHF_Endcap_minus=0;
2121 float pt_sum_NHF_Barrel=0;
2122 float pt_sum_NHF_Endcap_plus=0;
2123 float pt_sum_NHF_Endcap_minus=0;
2124 float pt_sum_PhF_Barrel=0;
2125 float pt_sum_PhF_Endcap_plus=0;
2126 float pt_sum_PhF_Endcap_minus=0;
2127 float pt_sum_HFH_plus=0;
2128 float pt_sum_HFH_minus=0;
2129 float pt_sum_HFE_plus=0;
2130 float pt_sum_HFE_minus=0;
2132 float px_chargedHadronsBarrel=0;
2133 float py_chargedHadronsBarrel=0;
2134 float px_chargedHadronsEndcapPlus=0;
2135 float py_chargedHadronsEndcapPlus=0;
2136 float px_chargedHadronsEndcapMinus=0;
2137 float py_chargedHadronsEndcapMinus=0;
2138 float px_neutralHadronsBarrel=0;
2139 float py_neutralHadronsBarrel=0;
2140 float px_neutralHadronsEndcapPlus=0;
2141 float py_neutralHadronsEndcapPlus=0;
2142 float px_neutralHadronsEndcapMinus=0;
2143 float py_neutralHadronsEndcapMinus=0;
2144 float px_PhotonsBarrel=0;
2145 float py_PhotonsBarrel=0;
2146 float px_PhotonsEndcapPlus=0;
2147 float py_PhotonsEndcapPlus=0;
2148 float px_PhotonsEndcapMinus=0;
2149 float py_PhotonsEndcapMinus=0;
2150 float px_HFHadronsPlus=0;
2151 float py_HFHadronsPlus=0;
2152 float px_HFHadronsMinus=0;
2153 float py_HFHadronsMinus=0;
2154 float px_HFEGammasPlus=0;
2155 float py_HFEGammasPlus=0;
2156 float px_HFEGammasMinus=0;
2157 float py_HFEGammasMinus=0;
2158 for (
unsigned int i = 0;
i < particleFlow->size();
i++) {
2162 if(c.
eta()>(-3.0)&& c.
eta()<(-1.392)){
2163 px_chargedHadronsEndcapMinus-=c.
px();
2164 py_chargedHadronsEndcapMinus-=c.
py();
2165 pt_sum_CHF_Endcap_minus+=c.
et();
2166 }
else if(c.
eta()>=(-1.392)&& c.
eta()<=1.392){
2167 px_chargedHadronsBarrel-=c.
px();
2168 py_chargedHadronsBarrel-=c.
py();
2169 pt_sum_CHF_Barrel+=c.
et();
2170 }
else if(c.
eta()>1.392 && c.
eta()<3.0){
2171 px_chargedHadronsEndcapPlus-=c.
px();
2172 py_chargedHadronsEndcapPlus-=c.
py();
2173 pt_sum_CHF_Endcap_plus+=c.
et();
2178 if(c.
eta()>(-3.0)&& c.
eta()<(-1.392)){
2179 px_neutralHadronsEndcapMinus-=c.
px();
2180 py_neutralHadronsEndcapMinus-=c.
py();
2181 pt_sum_NHF_Endcap_minus+=c.
et();
2182 }
else if(c.
eta()>=(-1.392)&& c.
eta()<=1.392){
2183 px_neutralHadronsBarrel-=c.
px();
2184 py_neutralHadronsBarrel-=c.
py();
2185 pt_sum_NHF_Barrel+=c.
et();
2186 }
else if(c.
eta()>1.392 && c.
eta()<3.0){
2187 px_neutralHadronsEndcapPlus-=c.
px();
2188 py_neutralHadronsEndcapPlus-=c.
py();
2189 pt_sum_NHF_Endcap_plus+=c.
et();
2194 if(c.
eta()>(-3.0)&& c.
eta()<(-1.479)){
2195 px_PhotonsEndcapMinus-=c.
px();
2196 py_PhotonsEndcapMinus-=c.
py();
2197 pt_sum_PhF_Endcap_minus+=c.
et();
2198 }
else if(c.
eta()>=(-1.479)&& c.
eta()<=1.479){
2199 px_PhotonsBarrel-=c.
px();
2200 py_PhotonsBarrel-=c.
py();
2201 pt_sum_PhF_Barrel+=c.
et();
2202 }
else if(c.
eta()>1.479 && c.
eta()<3.0){
2203 px_PhotonsEndcapPlus-=c.
px();
2204 py_PhotonsEndcapPlus-=c.
py();
2205 pt_sum_PhF_Endcap_plus+=c.
et();
2210 if(c.
eta()>(-5.20)&& c.
eta()<-2.901376){
2211 pt_sum_HFH_minus+=c.
et();
2212 px_HFHadronsMinus-=c.
px();
2213 py_HFHadronsMinus-=c.
py();
2214 }
else if(c.
eta()>2.901376 && c.
eta()<5.20){
2215 px_HFHadronsPlus-=c.
px();
2216 py_HFHadronsPlus-=c.
py();
2217 pt_sum_HFH_plus+=c.
et();
2222 if(c.
eta()>(-5.20)&& c.
eta()<-2.901376){
2223 pt_sum_HFE_minus+=c.
et();
2224 px_HFEGammasMinus-=c.
px();
2225 py_HFEGammasMinus-=c.
py();
2226 }
else if(c.
eta()>2.901376 && c.
eta()<5.20){
2227 px_HFEGammasPlus-=c.
px();
2228 py_HFEGammasPlus-=c.
py();
2229 pt_sum_HFE_plus+=c.
et();
2232 for (
unsigned int j=0;
j<typePFCand_.size();
j++) {
2235 if ( ((c.
eta()>etaMinPFCand_[
j]) && (c.
eta()<etaMaxPFCand_[
j])) || ((c.
eta()> (-etaMaxPFCand_[
j])) && (c.
eta()< (-etaMinPFCand_[
j]))) ){
2236 countsPFCand_[
j]+=1;
2237 MExPFCand_[
j]-=c.
px();
2238 MEyPFCand_[
j]-=c.
py();
2239 ptPFCand_[
j] = map_of_MEs[DirName +
"/"+ptPFCand_name_[
j]];
2240 if ( ptPFCand_[
j] && ptPFCand_[
j]->getRootObject()) ptPFCand_[
j]->Fill(c.
eta(), c.
phi(), c.
pt());
2241 occupancyPFCand_[
j] = map_of_MEs[DirName +
"/"+occupancyPFCand_name_[
j]];
2242 if ( occupancyPFCand_[
j] && occupancyPFCand_[
j]->getRootObject()) occupancyPFCand_[
j]->
Fill(c.
eta(), c.
phi());
2252 if ( (ecalRaw + hcalRaw) > hcalMin_ ){
2254 if( theElements.empty() )
continue;
2255 unsigned int iTrack=-999;
2256 std::vector<unsigned int> iECAL;
2257 std::vector<unsigned int> iHCAL;
2262 for(
unsigned int iEle=0; iEle<elements.
size(); iEle++) {
2266 case PFBlockElement::TRACK:
2271 iECAL.push_back( iEle );
2274 iHCAL.push_back( iEle );
2280 if ( nTracks == 1 ){
2283 mProfileIsoPFChHad_TrackOccupancy=map_of_MEs[DirName+
"/"+
"IsoPfChHad_Track_profile"];
2284 if (mProfileIsoPFChHad_TrackOccupancy && mProfileIsoPFChHad_TrackOccupancy->getRootObject()) mProfileIsoPFChHad_TrackOccupancy->Fill(et.
trackRef()->eta(),et.
trackRef()->phi());
2285 mProfileIsoPFChHad_TrackPt=map_of_MEs[DirName+
"/"+
"IsoPfChHad_TrackPt"];
2286 if (mProfileIsoPFChHad_TrackPt && mProfileIsoPFChHad_TrackPt->getRootObject()) mProfileIsoPFChHad_TrackPt->Fill(et.
trackRef()->eta(),et.
trackRef()->phi(),et.
trackRef()->pt());
2288 for(
unsigned int ii=0;
ii<iECAL.size();
ii++) {
2291 mProfileIsoPFChHad_EcalOccupancyCentral=map_of_MEs[DirName+
"/"+
"IsoPfChHad_ECAL_profile_central"];
2292 if (mProfileIsoPFChHad_EcalOccupancyCentral && mProfileIsoPFChHad_EcalOccupancyCentral->getRootObject()) mProfileIsoPFChHad_EcalOccupancyCentral->Fill(eecal.
clusterRef()->eta(),eecal.
clusterRef()->phi());
2293 mProfileIsoPFChHad_EMPtCentral=map_of_MEs[DirName+
"/"+
"IsoPfChHad_EMPt_central"];
2294 if (mProfileIsoPFChHad_EMPtCentral && mProfileIsoPFChHad_EMPtCentral->getRootObject()) mProfileIsoPFChHad_EMPtCentral->Fill(eecal.
clusterRef()->eta(),eecal.
clusterRef()->phi(),eecal.
clusterRef()->pt());
2296 mProfileIsoPFChHad_EcalOccupancyEndcap=map_of_MEs[DirName+
"/"+
"IsoPfChHad_ECAL_profile_endcap"];
2297 if (mProfileIsoPFChHad_EcalOccupancyEndcap && mProfileIsoPFChHad_EcalOccupancyEndcap->getRootObject()) mProfileIsoPFChHad_EcalOccupancyEndcap->Fill(eecal.
clusterRef()->eta(),eecal.
clusterRef()->phi());
2298 mProfileIsoPFChHad_EMPtEndcap=map_of_MEs[DirName+
"/"+
"IsoPfChHad_EMPt_endcap"];
2299 if (mProfileIsoPFChHad_EMPtEndcap && mProfileIsoPFChHad_EMPtEndcap->getRootObject()) mProfileIsoPFChHad_EMPtEndcap->Fill(eecal.
clusterRef()->eta(),eecal.
clusterRef()->phi(),eecal.
clusterRef()->pt());
2304 for(
unsigned int ii=0;
ii<iHCAL.size();
ii++) {
2307 mProfileIsoPFChHad_HcalOccupancyCentral=map_of_MEs[DirName+
"/"+
"IsoPfChHad_HCAL_profile_central"];
2308 if (mProfileIsoPFChHad_HcalOccupancyCentral && mProfileIsoPFChHad_HcalOccupancyCentral->getRootObject()) mProfileIsoPFChHad_HcalOccupancyCentral->Fill(ehcal.
clusterRef()->eta(),ehcal.
clusterRef()->phi());
2309 mProfileIsoPFChHad_HadPtCentral=map_of_MEs[DirName+
"/"+
"IsoPfChHad_HadPt_central"];
2310 if (mProfileIsoPFChHad_HadPtCentral && mProfileIsoPFChHad_HadPtCentral->getRootObject()) mProfileIsoPFChHad_HadPtCentral->Fill(ehcal.
clusterRef()->eta(),ehcal.
clusterRef()->phi(),ehcal.
clusterRef()->pt());
2312 mProfileIsoPFChHad_HcalOccupancyEndcap=map_of_MEs[DirName+
"/"+
"IsoPfChHad_HCAL_profile_endcap"];
2313 if (mProfileIsoPFChHad_HcalOccupancyEndcap && mProfileIsoPFChHad_HcalOccupancyEndcap->getRootObject()) mProfileIsoPFChHad_HcalOccupancyEndcap->Fill(ehcal.
clusterRef()->eta(),ehcal.
clusterRef()->phi());
2314 mProfileIsoPFChHad_HadPtEndcap=map_of_MEs[DirName+
"/"+
"IsoPfChHad_HadPt_endcap"];
2315 if (mProfileIsoPFChHad_HadPtEndcap && mProfileIsoPFChHad_HadPtEndcap->getRootObject()) mProfileIsoPFChHad_HadPtEndcap->Fill(ehcal.
clusterRef()->eta(),ehcal.
clusterRef()->phi(),ehcal.
clusterRef()->pt());
2323 for (
unsigned int j=0;
j<countsPFCand_.size();
j++) {
2324 profilePFCand_x_[
j] = map_of_MEs[DirName +
"/"+profilePFCand_x_name_[
j]];
2325 if(profilePFCand_x_[
j] && profilePFCand_x_[
j]->getRootObject()) profilePFCand_x_[
j]->
Fill(countsPFCand_[
j], MExPFCand_[j]);
2326 profilePFCand_y_[
j] = map_of_MEs[DirName +
"/"+profilePFCand_y_name_[
j]];
2327 if(profilePFCand_y_[j] && profilePFCand_y_[j]->getRootObject()) profilePFCand_y_[j]->
Fill(countsPFCand_[j], MEyPFCand_[j]);
2328 multiplicityPFCand_[
j] = map_of_MEs[DirName +
"/"+multiplicityPFCand_name_[
j]];
2329 if(multiplicityPFCand_[j] && multiplicityPFCand_[j]->getRootObject()) multiplicityPFCand_[j]->
Fill(countsPFCand_[j]);
2331 meCHF_Barrel=map_of_MEs[DirName+
"/"+
"PfChargedHadronEtFractionBarrel"];
if(meCHF_Barrel && meCHF_Barrel->getRootObject()) meCHF_Barrel->Fill(pt_sum_CHF_Barrel/pfmet.
sumEt());
2332 meCHF_EndcapPlus=map_of_MEs[DirName+
"/"+
"PfChargedHadronEtFractionEndcapPlus"];
if(meCHF_EndcapPlus && meCHF_EndcapPlus->getRootObject()) meCHF_EndcapPlus->Fill(pt_sum_CHF_Endcap_plus/pfmet.
sumEt());
2333 meCHF_EndcapMinus=map_of_MEs[DirName+
"/"+
"PfChargedHadronEtFractionEndcapMinus"];
if(meCHF_EndcapMinus && meCHF_EndcapMinus->getRootObject()) meCHF_EndcapMinus->Fill(pt_sum_CHF_Endcap_minus/pfmet.
sumEt());
2334 meNHF_Barrel=map_of_MEs[DirName+
"/"+
"PfNeutralHadronEtFractionBarrel"];
if(meNHF_Barrel && meNHF_Barrel->getRootObject()) meNHF_Barrel->Fill(pt_sum_NHF_Barrel/pfmet.
sumEt());
2335 meNHF_EndcapPlus=map_of_MEs[DirName+
"/"+
"PfNeutralHadronEtFractionEndcapPlus"];
if(meNHF_EndcapPlus && meNHF_EndcapPlus->getRootObject()) meNHF_EndcapPlus->Fill(pt_sum_NHF_Endcap_plus/pfmet.
sumEt());
2336 meNHF_EndcapMinus=map_of_MEs[DirName+
"/"+
"PfNeutralHadronEtFractionEndcapMinus"];
if(meNHF_EndcapMinus && meNHF_EndcapMinus->getRootObject()) meNHF_EndcapMinus->Fill(pt_sum_NHF_Endcap_minus/pfmet.
sumEt());
2337 mePhF_Barrel=map_of_MEs[DirName+
"/"+
"PfPhotonEtFractionBarrel"];
if(mePhF_Barrel && mePhF_Barrel->getRootObject()) mePhF_Barrel->Fill(pt_sum_PhF_Barrel/pfmet.
sumEt());
2338 mePhF_EndcapPlus=map_of_MEs[DirName+
"/"+
"PfPhotonEtFractionEndcapPlus"];
if(mePhF_EndcapPlus && mePhF_EndcapPlus->getRootObject()) mePhF_EndcapPlus->Fill(pt_sum_PhF_Endcap_plus/pfmet.
sumEt());
2339 mePhF_EndcapMinus=map_of_MEs[DirName+
"/"+
"PfPhotonEtFractionEndcapMinus"];
if(mePhF_EndcapMinus && mePhF_EndcapMinus->getRootObject()) mePhF_EndcapMinus->Fill(pt_sum_PhF_Endcap_minus/pfmet.
sumEt());
2340 meHFHadF_Plus=map_of_MEs[DirName+
"/"+
"PfHFHadronEtFractionPlus"];
if(meHFHadF_Plus && meHFHadF_Plus->getRootObject()) meHFHadF_Plus->Fill(pt_sum_HFH_plus/pfmet.
sumEt());
2341 meHFHadF_Minus=map_of_MEs[DirName+
"/"+
"PfHFHadronEtFractionMinus"];
if(meHFHadF_Minus && meHFHadF_Minus->getRootObject()) meHFHadF_Minus->Fill(pt_sum_HFH_minus/pfmet.
sumEt());
2342 meHFEMF_Plus=map_of_MEs[DirName+
"/"+
"PfHFEMEtFractionPlus"];
if(meHFEMF_Plus && meHFEMF_Plus->getRootObject()) meHFEMF_Plus->Fill(pt_sum_HFE_plus/pfmet.
sumEt());
2343 meHFEMF_Minus=map_of_MEs[DirName+
"/"+
"PfHFEMEtFractionMinus"];
if(meHFEMF_Minus && meHFEMF_Minus->getRootObject()) meHFEMF_Minus->Fill(pt_sum_HFE_minus/pfmet.
sumEt());
2347 if(pt_sum_CHF_Barrel){
2348 meMETPhiChargedHadronsBarrel = map_of_MEs[DirName+
"/"+
"METPhiChargedHadronsBarrel"];
if(meMETPhiChargedHadronsBarrel && meMETPhiChargedHadronsBarrel ->getRootObject())meMETPhiChargedHadronsBarrel->Fill(atan2(py_chargedHadronsBarrel,px_chargedHadronsBarrel));
2350 if(pt_sum_CHF_Endcap_plus){
2351 meMETPhiChargedHadronsEndcapPlus = map_of_MEs[DirName+
"/"+
"METPhiChargedHadronsEndcapPlus"];
if(meMETPhiChargedHadronsEndcapPlus && meMETPhiChargedHadronsEndcapPlus ->getRootObject())meMETPhiChargedHadronsEndcapPlus->Fill(atan2(py_chargedHadronsEndcapPlus,px_chargedHadronsEndcapPlus));
2353 if(pt_sum_CHF_Endcap_minus){
2354 meMETPhiChargedHadronsEndcapMinus = map_of_MEs[DirName+
"/"+
"METPhiChargedHadronsEndcapMinus"];
if(meMETPhiChargedHadronsEndcapMinus && meMETPhiChargedHadronsEndcapMinus ->getRootObject())meMETPhiChargedHadronsEndcapMinus->Fill(atan2(py_chargedHadronsEndcapMinus,px_chargedHadronsEndcapMinus));
2356 if(pt_sum_NHF_Barrel){
2357 meMETPhiNeutralHadronsBarrel = map_of_MEs[DirName+
"/"+
"METPhiNeutralHadronsBarrel"];
if(meMETPhiNeutralHadronsBarrel && meMETPhiNeutralHadronsBarrel ->getRootObject())meMETPhiNeutralHadronsBarrel->Fill(atan2(py_neutralHadronsBarrel,px_neutralHadronsBarrel));
2359 if(pt_sum_NHF_Endcap_plus){
2360 meMETPhiNeutralHadronsEndcapPlus = map_of_MEs[DirName+
"/"+
"METPhiNeutralHadronsEndcapPlus"];
if(meMETPhiNeutralHadronsEndcapPlus && meMETPhiNeutralHadronsEndcapPlus ->getRootObject())meMETPhiNeutralHadronsEndcapPlus->Fill(atan2(py_neutralHadronsEndcapPlus,px_neutralHadronsEndcapPlus));
2362 if(pt_sum_NHF_Endcap_minus){
2363 meMETPhiNeutralHadronsEndcapMinus = map_of_MEs[DirName+
"/"+
"METPhiNeutralHadronsEndcapMinus"];
if(meMETPhiNeutralHadronsEndcapMinus && meMETPhiNeutralHadronsEndcapMinus ->getRootObject())meMETPhiNeutralHadronsEndcapMinus->Fill(atan2(py_neutralHadronsEndcapMinus,px_neutralHadronsEndcapMinus));
2365 if(pt_sum_PhF_Barrel){
2366 meMETPhiPhotonsBarrel = map_of_MEs[DirName+
"/"+
"METPhiPhotonsBarrel"];
if(meMETPhiPhotonsBarrel && meMETPhiPhotonsBarrel ->getRootObject())meMETPhiPhotonsBarrel->Fill(atan2(py_PhotonsBarrel,px_PhotonsBarrel));
2368 if(pt_sum_PhF_Endcap_plus){
2369 meMETPhiPhotonsEndcapPlus = map_of_MEs[DirName+
"/"+
"METPhiPhotonsEndcapPlus"];
if(meMETPhiPhotonsEndcapPlus && meMETPhiPhotonsEndcapPlus ->getRootObject())meMETPhiPhotonsEndcapPlus->Fill(atan2(py_PhotonsEndcapPlus,px_PhotonsEndcapPlus));
2371 if(pt_sum_PhF_Endcap_minus){
2372 meMETPhiPhotonsEndcapMinus = map_of_MEs[DirName+
"/"+
"METPhiPhotonsEndcapMinus"];
if(meMETPhiPhotonsEndcapMinus && meMETPhiPhotonsEndcapMinus ->getRootObject())meMETPhiPhotonsEndcapMinus->Fill(atan2(py_PhotonsEndcapMinus,px_PhotonsEndcapMinus));
2374 if(pt_sum_HFH_plus){
2375 meMETPhiHFHadronsPlus = map_of_MEs[DirName+
"/"+
"METPhiHFHadronsPlus"];
if(meMETPhiHFHadronsPlus && meMETPhiHFHadronsPlus ->getRootObject())meMETPhiHFHadronsPlus->Fill(atan2(py_HFHadronsPlus,px_HFHadronsPlus));
2377 if(pt_sum_HFH_minus){
2378 meMETPhiHFHadronsMinus = map_of_MEs[DirName+
"/"+
"METPhiHFHadronsMinus"];
if(meMETPhiHFHadronsMinus && meMETPhiHFHadronsMinus ->getRootObject())meMETPhiHFHadronsMinus->Fill(atan2(py_HFHadronsMinus,px_HFHadronsMinus));
2380 if(pt_sum_HFE_plus){
2381 meMETPhiHFEGammasPlus = map_of_MEs[DirName+
"/"+
"METPhiHFEGammasPlus"];
if(meMETPhiHFEGammasPlus && meMETPhiHFEGammasPlus ->getRootObject())meMETPhiHFEGammasPlus->Fill(atan2(py_HFEGammasPlus,px_HFEGammasPlus));
2383 if(pt_sum_HFE_minus){
2384 meMETPhiHFEGammasMinus = map_of_MEs[DirName+
"/"+
"METPhiHFEGammasMinus"];
if(meMETPhiHFEGammasMinus && meMETPhiHFEGammasMinus ->getRootObject())meMETPhiHFEGammasMinus->Fill(atan2(py_HFEGammasMinus,px_HFEGammasMinus));
2446 if(techTriggerCase[1]){
2447 meCHF_Barrel_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfChargedHadronEtFractionBarrel_BXm1Filled"];
if(meCHF_Barrel_BXm1Filled && meCHF_Barrel_BXm1Filled->getRootObject()) meCHF_Barrel_BXm1Filled->Fill(pt_sum_CHF_Barrel/pfmet.
sumEt());
2448 meCHF_EndcapPlus_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfChargedHadronEtFractionEndcapPlus_BXm1Filled"];
if(meCHF_EndcapPlus_BXm1Filled && meCHF_EndcapPlus_BXm1Filled->getRootObject()) meCHF_EndcapPlus_BXm1Filled->Fill(pt_sum_CHF_Endcap_plus/pfmet.
sumEt());
2449 meCHF_EndcapMinus_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfChargedHadronEtFractionEndcapMinus_BXm1Filled"];
if(meCHF_EndcapMinus_BXm1Filled && meCHF_EndcapMinus_BXm1Filled->getRootObject()) meCHF_EndcapMinus_BXm1Filled->Fill(pt_sum_CHF_Endcap_minus/pfmet.
sumEt());
2450 meNHF_Barrel_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfNeutralHadronEtFractionBarrel_BXm1Filled"];
if(meNHF_Barrel_BXm1Filled && meNHF_Barrel_BXm1Filled->getRootObject()) meNHF_Barrel_BXm1Filled->Fill(pt_sum_NHF_Barrel/pfmet.
sumEt());
2451 meNHF_EndcapPlus_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfNeutralHadronEtFractionEndcapPlus_BXm1Filled"];
if(meNHF_EndcapPlus_BXm1Filled && meNHF_EndcapPlus_BXm1Filled->getRootObject()) meNHF_EndcapPlus_BXm1Filled->Fill(pt_sum_NHF_Endcap_plus/pfmet.
sumEt());
2452 meNHF_EndcapMinus_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfNeutralHadronEtFractionEndcapMinus_BXm1Filled"];
if(meNHF_EndcapMinus_BXm1Filled && meNHF_EndcapMinus_BXm1Filled->getRootObject()) meNHF_EndcapMinus_BXm1Filled->Fill(pt_sum_NHF_Endcap_minus/pfmet.
sumEt());
2453 mePhF_Barrel_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfPhotonEtFractionBarrel_BXm1Filled"];
if(mePhF_Barrel_BXm1Filled && mePhF_Barrel_BXm1Filled->getRootObject()) mePhF_Barrel_BXm1Filled->Fill(pt_sum_PhF_Barrel/pfmet.
sumEt());
2454 mePhF_EndcapPlus_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfPhotonEtFractionEndcapPlus_BXm1Filled"];
if(mePhF_EndcapPlus_BXm1Filled && mePhF_EndcapPlus_BXm1Filled->getRootObject()) mePhF_EndcapPlus_BXm1Filled->Fill(pt_sum_PhF_Endcap_plus/pfmet.
sumEt());
2455 mePhF_EndcapMinus_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfPhotonEtFractionEndcapMinus_BXm1Filled"];
if(mePhF_EndcapMinus_BXm1Filled && mePhF_EndcapMinus_BXm1Filled->getRootObject()) mePhF_EndcapMinus_BXm1Filled->Fill(pt_sum_PhF_Endcap_minus/pfmet.
sumEt());
2456 meHFHadF_Plus_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfHFHadronEtFractionPlus_BXm1Filled"];
if(meHFHadF_Plus_BXm1Filled && meHFHadF_Plus_BXm1Filled->getRootObject()) meHFHadF_Plus_BXm1Filled->Fill(pt_sum_HFH_plus/pfmet.
sumEt());
2457 meHFHadF_Minus_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfHFHadronEtFractionMinus_BXm1Filled"];
if(meHFHadF_Minus_BXm1Filled && meHFHadF_Minus_BXm1Filled->getRootObject()) meHFHadF_Minus_BXm1Filled->Fill(pt_sum_HFH_minus/pfmet.
sumEt());
2458 meHFEMF_Plus_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfHFEMEtFractionPlus_BXm1Filled"];
if(meHFEMF_Plus_BXm1Filled && meHFEMF_Plus_BXm1Filled->getRootObject()) meHFEMF_Plus_BXm1Filled->Fill(pt_sum_HFE_plus/pfmet.
sumEt());
2459 meHFEMF_Minus_BXm1Filled=map_of_MEs[DirName+
"/"+
"PfHFEMEtFractionMinus_BXm1Filled"];
if(meHFEMF_Minus_BXm1Filled && meHFEMF_Minus_BXm1Filled->getRootObject()) meHFEMF_Minus_BXm1Filled->Fill(pt_sum_HFE_minus/pfmet.
sumEt());
2460 mePhotonEtFraction_BXm1Filled = map_of_MEs[DirName+
"/"+
"PfPhotonEtFraction_BXm1Filled"];
if ( mePhotonEtFraction_BXm1Filled && mePhotonEtFraction_BXm1Filled ->getRootObject()) mePhotonEtFraction_BXm1Filled ->
Fill(pfmet.
photonEtFraction());
2461 meNeutralHadronEtFraction_BXm1Filled = map_of_MEs[DirName+
"/"+
"PfNeutralHadronEtFraction_BXm1Filled"];
if ( meNeutralHadronEtFraction_BXm1Filled && meNeutralHadronEtFraction_BXm1Filled ->getRootObject()) meNeutralHadronEtFraction_BXm1Filled ->
Fill(pfmet.
neutralHadronEtFraction());
2462 meChargedHadronEtFraction_BXm1Filled = map_of_MEs[DirName+
"/"+
"PfChargedHadronEtFraction_BXm1Filled"];
if ( meChargedHadronEtFraction_BXm1Filled && meChargedHadronEtFraction_BXm1Filled ->getRootObject()) meChargedHadronEtFraction_BXm1Filled ->
Fill(pfmet.
chargedHadronEtFraction());
2463 meMET_BXm1Filled = map_of_MEs[DirName+
"/"+
"MET_BXm1Filled"];
if ( meMET_BXm1Filled && meMET_BXm1Filled ->getRootObject()) meMET_BXm1Filled ->
Fill(pfmet.
pt());
2464 meSumET_BXm1Filled = map_of_MEs[DirName+
"/"+
"SumET_BXm1Filled"];
if ( meSumET_BXm1Filled && meSumET_BXm1Filled ->getRootObject()) meSumET_BXm1Filled ->
Fill(pfmet.
sumEt());
2465 if(pt_sum_CHF_Barrel){
2466 meMETPhiChargedHadronsBarrel_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiChargedHadronsBarrel_BXm1Filled"];
if(meMETPhiChargedHadronsBarrel_BXm1Filled && meMETPhiChargedHadronsBarrel_BXm1Filled ->getRootObject())meMETPhiChargedHadronsBarrel_BXm1Filled->Fill(atan2(py_chargedHadronsBarrel,px_chargedHadronsBarrel));
2468 if(pt_sum_CHF_Endcap_plus){
2469 meMETPhiChargedHadronsEndcapPlus_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiChargedHadronsEndcapPlus_BXm1Filled"];
if(meMETPhiChargedHadronsEndcapPlus_BXm1Filled && meMETPhiChargedHadronsEndcapPlus_BXm1Filled ->getRootObject())meMETPhiChargedHadronsEndcapPlus_BXm1Filled->Fill(atan2(py_chargedHadronsEndcapPlus,px_chargedHadronsEndcapPlus));
2471 if(pt_sum_CHF_Endcap_minus){
2472 meMETPhiChargedHadronsEndcapMinus_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiChargedHadronsEndcapMinus_BXm1Filled"];
if(meMETPhiChargedHadronsEndcapMinus_BXm1Filled && meMETPhiChargedHadronsEndcapMinus_BXm1Filled ->getRootObject())meMETPhiChargedHadronsEndcapMinus_BXm1Filled->Fill(atan2(py_chargedHadronsEndcapMinus,px_chargedHadronsEndcapMinus));
2474 if(pt_sum_NHF_Barrel){
2475 meMETPhiNeutralHadronsBarrel_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiNeutralHadronsBarrel_BXm1Filled"];
if(meMETPhiNeutralHadronsBarrel_BXm1Filled && meMETPhiNeutralHadronsBarrel_BXm1Filled ->getRootObject())meMETPhiNeutralHadronsBarrel_BXm1Filled->Fill(atan2(py_neutralHadronsBarrel,px_neutralHadronsBarrel));
2477 if(pt_sum_NHF_Endcap_plus){
2478 meMETPhiNeutralHadronsEndcapPlus_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiNeutralHadronsEndcapPlus_BXm1Filled"];
if(meMETPhiNeutralHadronsEndcapPlus_BXm1Filled && meMETPhiNeutralHadronsEndcapPlus_BXm1Filled ->getRootObject())meMETPhiNeutralHadronsEndcapPlus_BXm1Filled->Fill(atan2(py_neutralHadronsEndcapPlus,px_neutralHadronsEndcapPlus));
2480 if(pt_sum_NHF_Endcap_minus){
2481 meMETPhiNeutralHadronsEndcapMinus_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiNeutralHadronsEndcapMinus_BXm1Filled"];
if(meMETPhiNeutralHadronsEndcapMinus_BXm1Filled && meMETPhiNeutralHadronsEndcapMinus_BXm1Filled ->getRootObject())meMETPhiNeutralHadronsEndcapMinus_BXm1Filled->Fill(atan2(py_neutralHadronsEndcapMinus,px_neutralHadronsEndcapMinus));
2483 if(pt_sum_PhF_Barrel){
2484 meMETPhiPhotonsBarrel_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiPhotonsBarrel_BXm1Filled"];
if(meMETPhiPhotonsBarrel_BXm1Filled && meMETPhiPhotonsBarrel_BXm1Filled ->getRootObject())meMETPhiPhotonsBarrel_BXm1Filled->Fill(atan2(py_PhotonsBarrel,px_PhotonsBarrel));
2486 if(pt_sum_PhF_Endcap_plus){
2487 meMETPhiPhotonsEndcapPlus_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiPhotonsEndcapPlus_BXm1Filled"];
if(meMETPhiPhotonsEndcapPlus_BXm1Filled && meMETPhiPhotonsEndcapPlus_BXm1Filled ->getRootObject())meMETPhiPhotonsEndcapPlus_BXm1Filled->Fill(atan2(py_PhotonsEndcapPlus,px_PhotonsEndcapPlus));
2489 if(pt_sum_PhF_Endcap_minus){
2490 meMETPhiPhotonsEndcapMinus_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiPhotonsEndcapMinus_BXm1Filled"];
if(meMETPhiPhotonsEndcapMinus_BXm1Filled && meMETPhiPhotonsEndcapMinus_BXm1Filled ->getRootObject())meMETPhiPhotonsEndcapMinus_BXm1Filled->Fill(atan2(py_PhotonsEndcapMinus,px_PhotonsEndcapMinus));
2492 if(pt_sum_HFH_plus){
2493 meMETPhiHFHadronsPlus_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiHFHadronsPlus_BXm1Filled"];
if(meMETPhiHFHadronsPlus_BXm1Filled && meMETPhiHFHadronsPlus_BXm1Filled ->getRootObject())meMETPhiHFHadronsPlus_BXm1Filled->Fill(atan2(py_HFHadronsPlus,px_HFHadronsPlus));
2495 if(pt_sum_HFH_minus){
2496 meMETPhiHFHadronsMinus_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiHFHadronsMinus_BXm1Filled"];
if(meMETPhiHFHadronsMinus_BXm1Filled && meMETPhiHFHadronsMinus_BXm1Filled ->getRootObject())meMETPhiHFHadronsMinus_BXm1Filled->Fill(atan2(py_HFHadronsMinus,px_HFHadronsMinus));
2498 if(pt_sum_HFE_plus){
2499 meMETPhiHFEGammasPlus_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiHFEGammasPlus_BXm1Filled"];
if(meMETPhiHFEGammasPlus_BXm1Filled && meMETPhiHFEGammasPlus_BXm1Filled ->getRootObject())meMETPhiHFEGammasPlus_BXm1Filled->Fill(atan2(py_HFEGammasPlus,px_HFEGammasPlus));
2501 if(pt_sum_HFE_minus){
2502 meMETPhiHFEGammasMinus_BXm1Filled = map_of_MEs[DirName+
"/"+
"METPhiHFEGammasMinus_BXm1Filled"];
if(meMETPhiHFEGammasMinus_BXm1Filled && meMETPhiHFEGammasMinus_BXm1Filled ->getRootObject())meMETPhiHFEGammasMinus_BXm1Filled->Fill(atan2(py_HFEGammasMinus,px_HFEGammasMinus));
2564 if(techTriggerCase[2]){
2565 meCHF_Barrel_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfChargedHadronEtFractionBarrel_BXm1Empty"];
if(meCHF_Barrel_BXm1Empty && meCHF_Barrel_BXm1Empty->getRootObject()) meCHF_Barrel_BXm1Empty->Fill(pt_sum_CHF_Barrel/pfmet.
sumEt());
2566 meCHF_EndcapPlus_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfChargedHadronEtFractionEndcapPlus_BXm1Empty"];
if(meCHF_EndcapPlus_BXm1Empty && meCHF_EndcapPlus_BXm1Empty->getRootObject()) meCHF_EndcapPlus_BXm1Empty->Fill(pt_sum_CHF_Endcap_plus/pfmet.
sumEt());
2567 meCHF_EndcapMinus_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfChargedHadronEtFractionEndcapMinus_BXm1Empty"];
if(meCHF_EndcapMinus_BXm1Empty && meCHF_EndcapMinus_BXm1Empty->getRootObject()) meCHF_EndcapMinus_BXm1Empty->Fill(pt_sum_CHF_Endcap_minus/pfmet.
sumEt());
2568 meNHF_Barrel_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfNeutralHadronEtFractionBarrel_BXm1Empty"];
if(meNHF_Barrel_BXm1Empty && meNHF_Barrel_BXm1Empty->getRootObject()) meNHF_Barrel_BXm1Empty->Fill(pt_sum_NHF_Barrel/pfmet.
sumEt());
2569 meNHF_EndcapPlus_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfNeutralHadronEtFractionEndcapPlus_BXm1Empty"];
if(meNHF_EndcapPlus_BXm1Empty && meNHF_EndcapPlus_BXm1Empty->getRootObject()) meNHF_EndcapPlus_BXm1Empty->Fill(pt_sum_NHF_Endcap_plus/pfmet.
sumEt());
2570 meNHF_EndcapMinus_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfNeutralHadronEtFractionEndcapMinus_BXm1Empty"];
if(meNHF_EndcapMinus_BXm1Empty && meNHF_EndcapMinus_BXm1Empty->getRootObject()) meNHF_EndcapMinus_BXm1Empty->Fill(pt_sum_NHF_Endcap_minus/pfmet.
sumEt());
2571 mePhF_Barrel_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfPhotonEtFractionBarrel_BXm1Empty"];
if(mePhF_Barrel_BXm1Empty && mePhF_Barrel_BXm1Empty->getRootObject()) mePhF_Barrel_BXm1Empty->Fill(pt_sum_PhF_Barrel/pfmet.
sumEt());
2572 mePhF_EndcapPlus_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfPhotonEtFractionEndcapPlus_BXm1Empty"];
if(mePhF_EndcapPlus_BXm1Empty && mePhF_EndcapPlus_BXm1Empty->getRootObject()) mePhF_EndcapPlus_BXm1Empty->Fill(pt_sum_PhF_Endcap_plus/pfmet.
sumEt());
2573 mePhF_EndcapMinus_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfPhotonEtFractionEndcapMinus_BXm1Empty"];
if(mePhF_EndcapMinus_BXm1Empty && mePhF_EndcapMinus_BXm1Empty->getRootObject()) mePhF_EndcapMinus_BXm1Empty->Fill(pt_sum_PhF_Endcap_minus/pfmet.
sumEt());
2574 meHFHadF_Plus_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfHFHadronEtFractionPlus_BXm1Empty"];
if(meHFHadF_Plus_BXm1Empty && meHFHadF_Plus_BXm1Empty->getRootObject()) meHFHadF_Plus_BXm1Empty->Fill(pt_sum_HFH_plus/pfmet.
sumEt());
2575 meHFHadF_Minus_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfHFHadronEtFractionMinus_BXm1Empty"];
if(meHFHadF_Minus_BXm1Empty && meHFHadF_Minus_BXm1Empty->getRootObject()) meHFHadF_Minus_BXm1Empty->Fill(pt_sum_HFH_minus/pfmet.
sumEt());
2576 meHFEMF_Plus_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfHFEMEtFractionPlus_BXm1Empty"];
if(meHFEMF_Plus_BXm1Empty && meHFEMF_Plus_BXm1Empty->getRootObject()) meHFEMF_Plus_BXm1Empty->Fill(pt_sum_HFE_plus/pfmet.
sumEt());
2577 meHFEMF_Minus_BXm1Empty=map_of_MEs[DirName+
"/"+
"PfHFEMEtFractionMinus_BXm1Empty"];
if(meHFEMF_Minus_BXm1Empty && meHFEMF_Minus_BXm1Empty->getRootObject()) meHFEMF_Minus_BXm1Empty->Fill(pt_sum_HFE_minus/pfmet.
sumEt());
2578 mePhotonEtFraction_BXm1Empty = map_of_MEs[DirName+
"/"+
"PfPhotonEtFraction_BXm1Empty"];
if ( mePhotonEtFraction_BXm1Empty && mePhotonEtFraction_BXm1Empty ->getRootObject()) mePhotonEtFraction_BXm1Empty ->
Fill(pfmet.
photonEtFraction());
2579 meNeutralHadronEtFraction_BXm1Empty = map_of_MEs[DirName+
"/"+
"PfNeutralHadronEtFraction_BXm1Empty"];
if ( meNeutralHadronEtFraction_BXm1Empty && meNeutralHadronEtFraction_BXm1Empty ->getRootObject()) meNeutralHadronEtFraction_BXm1Empty ->
Fill(pfmet.
neutralHadronEtFraction());
2580 meChargedHadronEtFraction_BXm1Empty = map_of_MEs[DirName+
"/"+
"PfChargedHadronEtFraction_BXm1Empty"];
if ( meChargedHadronEtFraction_BXm1Empty && meChargedHadronEtFraction_BXm1Empty ->getRootObject()) meChargedHadronEtFraction_BXm1Empty ->
Fill(pfmet.
chargedHadronEtFraction());
2581 meMET_BXm1Empty = map_of_MEs[DirName+
"/"+
"MET_BXm1Empty"];
if ( meMET_BXm1Empty && meMET_BXm1Empty ->getRootObject()) meMET_BXm1Empty ->
Fill(pfmet.
pt());
2582 meSumET_BXm1Empty = map_of_MEs[DirName+
"/"+
"SumET_BXm1Empty"];
if ( meSumET_BXm1Empty && meSumET_BXm1Empty ->getRootObject()) meSumET_BXm1Empty ->
Fill(pfmet.
sumEt());
2583 if(pt_sum_CHF_Barrel){
2584 meMETPhiChargedHadronsBarrel_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiChargedHadronsBarrel_BXm1Empty"];
if(meMETPhiChargedHadronsBarrel_BXm1Empty && meMETPhiChargedHadronsBarrel_BXm1Empty ->getRootObject())meMETPhiChargedHadronsBarrel_BXm1Empty->Fill(atan2(py_chargedHadronsBarrel,px_chargedHadronsBarrel));
2586 if(pt_sum_CHF_Endcap_plus){
2587 meMETPhiChargedHadronsEndcapPlus_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiChargedHadronsEndcapPlus_BXm1Empty"];
if(meMETPhiChargedHadronsEndcapPlus_BXm1Empty && meMETPhiChargedHadronsEndcapPlus_BXm1Empty ->getRootObject())meMETPhiChargedHadronsEndcapPlus_BXm1Empty->Fill(atan2(py_chargedHadronsEndcapPlus,px_chargedHadronsEndcapPlus));
2589 if(pt_sum_CHF_Endcap_minus){
2590 meMETPhiChargedHadronsEndcapMinus_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiChargedHadronsEndcapMinus_BXm1Empty"];
if(meMETPhiChargedHadronsEndcapMinus_BXm1Empty && meMETPhiChargedHadronsEndcapMinus_BXm1Empty ->getRootObject())meMETPhiChargedHadronsEndcapMinus_BXm1Empty->Fill(atan2(py_chargedHadronsEndcapMinus,px_chargedHadronsEndcapMinus));
2592 if(pt_sum_NHF_Barrel){
2593 meMETPhiNeutralHadronsBarrel_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiNeutralHadronsBarrel_BXm1Empty"];
if(meMETPhiNeutralHadronsBarrel_BXm1Empty && meMETPhiNeutralHadronsBarrel_BXm1Empty ->getRootObject())meMETPhiNeutralHadronsBarrel_BXm1Empty->Fill(atan2(py_neutralHadronsBarrel,px_neutralHadronsBarrel));
2595 if(pt_sum_NHF_Endcap_plus){
2596 meMETPhiNeutralHadronsEndcapPlus_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiNeutralHadronsEndcapPlus_BXm1Empty"];
if(meMETPhiNeutralHadronsEndcapPlus_BXm1Empty && meMETPhiNeutralHadronsEndcapPlus_BXm1Empty ->getRootObject())meMETPhiNeutralHadronsEndcapPlus_BXm1Empty->Fill(atan2(py_neutralHadronsEndcapPlus,px_neutralHadronsEndcapPlus));
2598 if(pt_sum_NHF_Endcap_minus){
2599 meMETPhiNeutralHadronsEndcapMinus_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiNeutralHadronsEndcapMinus_BXm1Empty"];
if(meMETPhiNeutralHadronsEndcapMinus_BXm1Empty && meMETPhiNeutralHadronsEndcapMinus_BXm1Empty ->getRootObject())meMETPhiNeutralHadronsEndcapMinus_BXm1Empty->Fill(atan2(py_neutralHadronsEndcapMinus,px_neutralHadronsEndcapMinus));
2601 if(pt_sum_PhF_Barrel){
2602 meMETPhiPhotonsBarrel_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiPhotonsBarrel_BXm1Empty"];
if(meMETPhiPhotonsBarrel_BXm1Empty && meMETPhiPhotonsBarrel_BXm1Empty ->getRootObject())meMETPhiPhotonsBarrel_BXm1Empty->Fill(atan2(py_PhotonsBarrel,px_PhotonsBarrel));
2604 if(pt_sum_PhF_Endcap_plus){
2605 meMETPhiPhotonsEndcapPlus_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiPhotonsEndcapPlus_BXm1Empty"];
if(meMETPhiPhotonsEndcapPlus_BXm1Empty && meMETPhiPhotonsEndcapPlus_BXm1Empty ->getRootObject())meMETPhiPhotonsEndcapPlus_BXm1Empty->Fill(atan2(py_PhotonsEndcapPlus,px_PhotonsEndcapPlus));
2607 if(pt_sum_PhF_Endcap_minus){
2608 meMETPhiPhotonsEndcapMinus_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiPhotonsEndcapMinus_BXm1Empty"];
if(meMETPhiPhotonsEndcapMinus_BXm1Empty && meMETPhiPhotonsEndcapMinus_BXm1Empty ->getRootObject())meMETPhiPhotonsEndcapMinus_BXm1Empty->Fill(atan2(py_PhotonsEndcapMinus,px_PhotonsEndcapMinus));
2610 if(pt_sum_HFH_plus){
2611 meMETPhiHFHadronsPlus_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiHFHadronsPlus_BXm1Empty"];
if(meMETPhiHFHadronsPlus_BXm1Empty && meMETPhiHFHadronsPlus_BXm1Empty ->getRootObject())meMETPhiHFHadronsPlus_BXm1Empty->Fill(atan2(py_HFHadronsPlus,px_HFHadronsPlus));
2613 if(pt_sum_HFH_minus){
2614 meMETPhiHFHadronsMinus_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiHFHadronsMinus_BXm1Empty"];
if(meMETPhiHFHadronsMinus_BXm1Empty && meMETPhiHFHadronsMinus_BXm1Empty ->getRootObject())meMETPhiHFHadronsMinus_BXm1Empty->Fill(atan2(py_HFHadronsMinus,px_HFHadronsMinus));
2616 if(pt_sum_HFE_plus){
2617 meMETPhiHFEGammasPlus_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiHFEGammasPlus_BXm1Empty"];
if(meMETPhiHFEGammasPlus_BXm1Empty && meMETPhiHFEGammasPlus_BXm1Empty ->getRootObject())meMETPhiHFEGammasPlus_BXm1Empty->Fill(atan2(py_HFEGammasPlus,px_HFEGammasPlus)); \
2619 if(pt_sum_HFE_minus){
2620 meMETPhiHFEGammasMinus_BXm1Empty = map_of_MEs[DirName+
"/"+
"METPhiHFEGammasMinus_BXm1Empty"];
if(meMETPhiHFEGammasMinus_BXm1Empty && meMETPhiHFEGammasMinus_BXm1Empty ->getRootObject())meMETPhiHFEGammasMinus_BXm1Empty->Fill(atan2(py_HFEGammasMinus,px_HFEGammasMinus));
2629 double pfPhotonEt = pfmet.
photonEt();
2635 double pfMuonEt = pfmet.
muonEt();
2639 double pfHFEMEt = pfmet.
HFEMEt();
2640 mePhotonEtFraction = map_of_MEs[DirName +
"/PfPhotonEtFraction"];
2641 mePhotonEt = map_of_MEs[DirName +
"/PfPhotonEt"];
2642 meNeutralHadronEtFraction = map_of_MEs[DirName +
"/PfNeutralHadronEtFraction"];
2643 meNeutralHadronEt = map_of_MEs[DirName +
"/PfNeutralHadronEt"];
2644 meElectronEt = map_of_MEs[DirName +
"/PfElectronEt"];
2645 meChargedHadronEtFraction = map_of_MEs[DirName +
"/PfChargedHadronEtFraction"];
2646 meChargedHadronEt = map_of_MEs[DirName +
"/PfChargedHadronEt"];
2647 meMuonEt = map_of_MEs[DirName +
"/PfMuonEt"];
2648 meHFHadronEtFraction = map_of_MEs[DirName +
"/PfHFHadronEtFraction"];
2649 meHFHadronEt = map_of_MEs[DirName +
"/PfHFHadronEt"];
2650 meHFEMEtFraction = map_of_MEs[DirName +
"/PfHFEMEtFraction"];
2651 meHFEMEt = map_of_MEs[DirName +
"/PfHFEMEt"];
2653 if (mePhotonEtFraction && mePhotonEtFraction ->getRootObject()) mePhotonEtFraction ->Fill(pfPhotonEtFraction);
2654 if (mePhotonEt && mePhotonEt ->getRootObject()) mePhotonEt ->Fill(pfPhotonEt);
2655 if (meNeutralHadronEtFraction && meNeutralHadronEtFraction->getRootObject()) meNeutralHadronEtFraction->Fill(pfNeutralHadronEtFraction);
2656 if (meNeutralHadronEt && meNeutralHadronEt ->getRootObject()) meNeutralHadronEt ->Fill(pfNeutralHadronEt);
2657 if (meElectronEt && meElectronEt ->getRootObject()) meElectronEt ->Fill(pfElectronEt);
2658 if (meChargedHadronEtFraction && meChargedHadronEtFraction->getRootObject()) meChargedHadronEtFraction->Fill(pfChargedHadronEtFraction);
2659 if (meChargedHadronEt && meChargedHadronEt ->getRootObject()) meChargedHadronEt ->Fill(pfChargedHadronEt);
2660 if (meMuonEt && meMuonEt ->getRootObject()) meMuonEt ->Fill(pfMuonEt);
2661 if (meHFHadronEtFraction && meHFHadronEtFraction ->getRootObject()) meHFHadronEtFraction ->Fill(pfHFHadronEtFraction);
2662 if (meHFHadronEt && meHFHadronEt ->getRootObject()) meHFHadronEt ->Fill(pfHFHadronEt);
2663 if (meHFEMEtFraction && meHFEMEtFraction ->getRootObject()) meHFEMEtFraction ->Fill(pfHFEMEtFraction);
2664 if (meHFEMEt && meHFEMEt ->getRootObject()) meHFEMEt ->Fill(pfHFEMEt);
2668 mePhotonEtFraction_profile = map_of_MEs[DirName +
"/PfPhotonEtFraction_profile"];
2669 mePhotonEt_profile = map_of_MEs[DirName +
"/PfPhotonEt_profile"];
2670 meNeutralHadronEtFraction_profile = map_of_MEs[DirName +
"/PfNeutralHadronEtFraction_profile"];
2671 meNeutralHadronEt_profile = map_of_MEs[DirName +
"/PfNeutralHadronEt_profile"];
2672 meChargedHadronEtFraction_profile = map_of_MEs[DirName +
"/PfChargedHadronEtFraction_profile"];
2673 meChargedHadronEt_profile = map_of_MEs[DirName +
"/PfChargedHadronEt_profile"];
2674 meHFHadronEtFraction_profile = map_of_MEs[DirName +
"/PfHFHadronEtFraction_profile"];
2675 meHFHadronEt_profile = map_of_MEs[DirName +
"/PfHFHadronEt_profile"];
2676 meHFEMEtFraction_profile = map_of_MEs[DirName +
"/PfHFEMEtFraction_profile"];
2677 meHFEMEt_profile = map_of_MEs[DirName +
"/PfHFEMEt_profile"];
2679 if (mePhotonEtFraction_profile && mePhotonEtFraction_profile ->getRootObject()) mePhotonEtFraction_profile ->Fill(numPV_, pfPhotonEtFraction);
2680 if (mePhotonEt_profile && mePhotonEt_profile ->getRootObject()) mePhotonEt_profile ->Fill(numPV_, pfPhotonEt);
2681 if (meNeutralHadronEtFraction_profile && meNeutralHadronEtFraction_profile->getRootObject()) meNeutralHadronEtFraction_profile->Fill(numPV_, pfNeutralHadronEtFraction);
2682 if (meNeutralHadronEt_profile && meNeutralHadronEt_profile ->getRootObject()) meNeutralHadronEt_profile ->Fill(numPV_, pfNeutralHadronEt);
2683 if (meChargedHadronEtFraction_profile && meChargedHadronEtFraction_profile->getRootObject()) meChargedHadronEtFraction_profile->Fill(numPV_, pfChargedHadronEtFraction);
2684 if (meChargedHadronEt_profile && meChargedHadronEt_profile ->getRootObject()) meChargedHadronEt_profile ->Fill(numPV_, pfChargedHadronEt);
2685 if (meHFHadronEtFraction_profile && meHFHadronEtFraction_profile ->getRootObject()) meHFHadronEtFraction_profile ->Fill(numPV_, pfHFHadronEtFraction);
2686 if (meHFHadronEt_profile && meHFHadronEt_profile ->getRootObject()) meHFHadronEt_profile ->Fill(numPV_, pfHFHadronEt);
2687 if (meHFEMEtFraction_profile && meHFEMEtFraction_profile ->getRootObject()) meHFEMEtFraction_profile ->Fill(numPV_, pfHFEMEtFraction);
2688 if (meHFEMEt_profile && meHFEMEt_profile ->getRootObject()) meHFEMEt_profile ->Fill(numPV_, pfHFEMEt);
2692 mePhotonEtFraction = map_of_MEs[DirName +
"/PfPhotonEtFraction"];
2693 meNeutralHadronEtFraction = map_of_MEs[DirName +
"/PfNeutralHadronEtFraction"];
2694 meChargedHadronEtFraction = map_of_MEs[DirName +
"/PfChargedHadronEtFraction"];
2695 meHFHadronEtFraction = map_of_MEs[DirName +
"/PfHFHadronEtFraction"];
2696 meHFEMEtFraction = map_of_MEs[DirName +
"/PfHFEMEtFraction"];
2698 if (mePhotonEtFraction && mePhotonEtFraction ->getRootObject()) mePhotonEtFraction ->Fill(patmet.
NeutralEMFraction());
2699 if (meNeutralHadronEtFraction && meNeutralHadronEtFraction->getRootObject()) meNeutralHadronEtFraction->Fill(patmet.
NeutralHadEtFraction());
2700 if (meChargedHadronEtFraction && meChargedHadronEtFraction->getRootObject()) meChargedHadronEtFraction->Fill(patmet.
ChargedHadEtFraction());
2701 if (meHFHadronEtFraction && meHFHadronEtFraction ->getRootObject()) meHFHadronEtFraction ->
Fill(patmet.
Type6EtFraction());
2702 if (meHFEMEtFraction && meHFEMEtFraction ->getRootObject()) meHFEMEtFraction ->
Fill(patmet.
Type7EtFraction());
2705 mePhotonEtFraction_profile = map_of_MEs[DirName +
"/PfPhotonEtFraction_profile"];
2706 meNeutralHadronEtFraction_profile = map_of_MEs[DirName +
"/PfNeutralHadronEtFraction_profile"];
2707 meChargedHadronEtFraction_profile = map_of_MEs[DirName +
"/PfChargedHadronEtFraction_profile"];
2708 meHFHadronEtFraction_profile = map_of_MEs[DirName +
"/PfHFHadronEtFraction_profile"];
2709 meHFEMEtFraction_profile = map_of_MEs[DirName +
"/PfHFEMEtFraction_profile"];
2711 if (mePhotonEtFraction_profile && mePhotonEtFraction_profile ->getRootObject()) mePhotonEtFraction_profile ->
Fill(numPV_, patmet.
NeutralEMFraction());
2712 if (meNeutralHadronEtFraction_profile && meNeutralHadronEtFraction_profile->getRootObject()) meNeutralHadronEtFraction_profile->Fill(numPV_, patmet.
NeutralHadEtFraction());
2713 if (meChargedHadronEtFraction_profile && meChargedHadronEtFraction_profile->getRootObject()) meChargedHadronEtFraction_profile->Fill(numPV_, patmet.
ChargedHadEtFraction());
2714 if (meHFHadronEtFraction_profile && meHFHadronEtFraction_profile ->getRootObject()) meHFHadronEtFraction_profile ->
Fill(numPV_, patmet.
Type6EtFraction());
2715 if (meHFEMEtFraction_profile && meHFEMEtFraction_profile ->getRootObject()) meHFEMEtFraction_profile ->
Fill(numPV_, patmet.
Type7EtFraction());
2717 mePhotonEt = map_of_MEs[DirName +
"/PfPhotonEt"];
2718 meNeutralHadronEt = map_of_MEs[DirName +
"/PfNeutralHadronEt"];
2719 meChargedHadronEt = map_of_MEs[DirName +
"/PfChargedHadronEt"];
2720 meHFHadronEt = map_of_MEs[DirName +
"/PfHFHadronEt"];
2721 meHFEMEt = map_of_MEs[DirName +
"/PfHFEMEt"];
2722 meMuonEt = map_of_MEs[DirName +
"/PfMuonEt"];
2723 meElectronEt = map_of_MEs[DirName +
"/PfElectronEt"];
2726 if (meNeutralHadronEt && meNeutralHadronEt->getRootObject()) meNeutralHadronEt->Fill(patmet.
NeutralHadEtFraction()*patmet.
sumEt());
2727 if (meChargedHadronEt && meChargedHadronEt->getRootObject()) meChargedHadronEt->Fill(patmet.
ChargedHadEtFraction()*patmet.
sumEt());
2734 mePhotonEt_profile = map_of_MEs[DirName +
"/PfPhotonEt_profile"];
2735 meNeutralHadronEt_profile = map_of_MEs[DirName +
"/PfNeutralHadronEt_profile"];
2736 meChargedHadronEt_profile = map_of_MEs[DirName +
"/PfChargedHadronEt_profile"];
2737 meHFHadronEt_profile = map_of_MEs[DirName +
"/PfHFHadronEt_profile"];
2738 meHFEMEt_profile = map_of_MEs[DirName +
"/PfHFEMEt_profile"];
2740 if (mePhotonEt_profile && mePhotonEt_profile ->getRootObject()) mePhotonEt_profile ->
Fill(numPV_, patmet.
NeutralEMFraction()*patmet.
sumEt());
2741 if (meNeutralHadronEt_profile && meNeutralHadronEt_profile->getRootObject()) meNeutralHadronEt_profile->Fill(numPV_, patmet.
NeutralHadEtFraction()*patmet.
sumEt());
2742 if (meChargedHadronEt_profile && meChargedHadronEt_profile->getRootObject()) meChargedHadronEt_profile->Fill(numPV_, patmet.
ChargedHadEtFraction()*patmet.
sumEt());
2743 if (meHFHadronEt_profile && meHFHadronEt_profile ->getRootObject()) meHFHadronEt_profile ->
Fill(numPV_, patmet.
Type6EtFraction()*patmet.
sumEt());
2744 if (meHFEMEt_profile && meHFEMEt_profile ->getRootObject()) meHFEMEt_profile ->
Fill(numPV_, patmet.
Type7EtFraction()*patmet.
sumEt());
2749 if (fill_met_high_level_histo){
2750 hMExLS = map_of_MEs[DirName+
"/"+
"MExLS"];
if (hMExLS && hMExLS->getRootObject()) hMExLS->Fill(MEx,myLuminosityBlock);
2751 hMEyLS = map_of_MEs[DirName+
"/"+
"MEyLS"];
if (hMEyLS && hMEyLS->getRootObject()) hMEyLS->Fill(MEy,myLuminosityBlock);
float puppiWeight() const
Set both weights at once (with option for only full PUPPI)
Analysis-level MET class.
T getParameter(std::string const &) const
const reco::TrackRef & trackRef() const
const unsigned int nTracks(const reco::Vertex &sv)
void setBinContent(int binx, double content)
set content of bin (1-D)
double HFEMEtFraction() const
double rawEcalEnergy() const
return corrected Ecal energy
double ChargedHadEtFraction() const
float puppiWeightNoLep() const
Weight from full PUPPI.
double MuonEtFraction() const
MonitorElement * bookProfile(Args &&...args)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
std::vector< T >::const_iterator search(const cond::Time_t &val, const std::vector< T > &container)
std::vector< ParameterSet > VParameterSet
double neutralHadronEtFraction() const
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)
edm::LuminosityBlockNumber_t luminosityBlock() const
virtual double phi() const final
momentum azimuthal angle
std::vector< Vertex > VertexCollection
collection of Vertex objects
const PFClusterRef & clusterRef() const
double NeutralEMFraction() const
Provides a code based selection for trigger and DCS information in order to have no failing filters i...
double ChargedEMEtFraction() const
virtual double pt() const
transverse momentum
double Type6EtFraction() const
double Type7EtFraction() const
virtual double eta() const
momentum pseudorapidity
double NeutralHadEtFraction() const
std::vector< ElementInBlock > ElementsInBlocks
triggerResultsToken_(consumes< edm::TriggerResults >(edm::InputTag("TriggerResults")))
virtual double phi() const
momentum azimuthal angle
double HFHadronEtFraction() const
double HFHadronEt() const
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
double photonEtFraction() const
MonitorElement * book1D(Args &&...args)
std::vector< bool > TechnicalTriggerWord
technical trigger bits (64 bits)
virtual double py() const final
y coordinate of momentum vector
PF Jet selector for pat::Jets.
double electronEt() const
Jet selector for pat::Jets and for CaloJets.
TObject * getRootObject(void) const
void setCurrentFolder(const std::string &fullpath)
double etFractionHadronic() const
T const * product() const
MonitorElement * book2D(Args &&...args)
static std::string const triggerResults("TriggerResults")
T const * product() const
TH1F * getTH1F(void) const
Geom::Phi< T > phi() const
double chargedHadronEtFraction() const
Particle reconstructed by the particle flow algorithm.
static int position[264][3]
virtual double px() const final
x coordinate of momentum vector
virtual double et() const final
transverse energy
double neutralHadronEt() const
virtual double eta() const final
momentum pseudorapidity
virtual ParticleType particleId() const
virtual int pdgId() const
PDG identifier.
const ElementsInBlocks & elementsInBlocks() const
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
double emEtFraction() const
double chargedHadronEt() const
double rawHcalEnergy() const
return raw Hcal energy
virtual double pt() const final
transverse momentum
math::PtEtaPhiMLorentzVector PolarLorentzVector
Lorentz vector.
1.8.5
