Inizialize parameters for histo binning.
156 int etaNBinsPFCand =
v->getParameter<
int>(
"etaNBins");
157 double etaMinPFCand =
v->getParameter<
double>(
"etaMin");
158 double etaMaxPFCand =
v->getParameter<
double>(
"etaMax");
159 int phiNBinsPFCand =
v->getParameter<
int>(
"phiNBins");
160 double phiMinPFCand =
v->getParameter<
double>(
"phiMin");
161 double phiMaxPFCand =
v->getParameter<
double>(
"phiMax");
162 int nMinPFCand =
v->getParameter<
int>(
"nMin");
163 int nMaxPFCand =
v->getParameter<
int>(
"nMax");
164 int nbinsPFCand =
v->getParameter<
double>(
"nbins");
177 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
188 if (etaMinPFCand * etaMaxPFCand < 0) {
207 const int nbins_eta_endcap = 2 * (etaNBinsPFCand + 1);
208 double eta_limits_endcap[nbins_eta_endcap];
209 for (
int i = 0;
i < nbins_eta_endcap;
i++) {
210 if (
i < (etaNBinsPFCand + 1)) {
211 eta_limits_endcap[
i] = etaMinPFCand +
i * (etaMaxPFCand - etaMinPFCand) / (
double)etaNBinsPFCand;
213 eta_limits_endcap[
i] =
214 -etaMaxPFCand + (
i - (etaNBinsPFCand + 1)) * (etaMaxPFCand - etaMinPFCand) / (
double)etaNBinsPFCand;
219 nbins_eta_endcap - 1,
228 nbins_eta_endcap - 1,
236 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
239 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
244 "IsoPfChHad_Track_profile",
"Isolated PFChHadron Tracker_occupancy", 108, -2.7, 2.7, 160, -
M_PI,
M_PI);
245 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"IsoPfChHad_Track_profile",
248 ibooker.book2D(
"IsoPfChHad_TrackPt",
"Isolated PFChHadron TrackPt", 108, -2.7, 2.7, 160, -
M_PI,
M_PI);
253 "IsolatedPFChHa ECAL occupancy (Barrel)",
260 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"IsoPfChHad_ECAL_profile_central",
263 "IsoPfChHad_EMPt_central",
"Isolated PFChHadron HadPt_central", 180, -1.479, 1.479, 360, -
M_PI,
M_PI);
264 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"IsoPfChHad_EMPt_central",
268 "IsoPfChHad_ECAL_profile_endcap",
"IsolatedPFChHa ECAL occupancy (Endcap)", 110, -2.75, 2.75, 125, -
M_PI,
M_PI);
269 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"IsoPfChHad_ECAL_profile_endcap",
272 ibooker.book2D(
"IsoPfChHad_EMPt_endcap",
"Isolated PFChHadron EMPt_endcap", 110, -2.75, 2.75, 125, -
M_PI,
M_PI);
273 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"IsoPfChHad_EMPt_endcap",
295 TH2F* hist_temp_HCAL =
new TH2F(
"IsoPfChHad_HCAL_profile_endcap",
296 "IsolatedPFChHa HCAL occupancy (outer endcap)",
302 TH2F* hist_tempPt_HCAL = (TH2F*)hist_temp_HCAL->Clone(
"Isolated PFCHHadron HadPt (outer endcap)");
305 "IsolatedPFChHa HCAL occupancy (Central Part)",
312 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"IsoPfChHad_HCAL_profile_central",
315 "IsoPfChHad_HadPt_central",
"Isolated PFChHadron HadPt_central", 40, -1.740, 1.740, 72, -
M_PI,
M_PI);
316 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"IsoPfChHad_HadPt_central",
320 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"IsoPfChHad_HCAL_profile_endcap",
323 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"IsoPfChHad_HadPt_endcap",
335 const int nbins_eta_hcal_total = 28;
336 double eta_limits_hcal_total[nbins_eta_hcal_total] = {
337 -2.650, -2.322, -2.043, -1.830, -1.653, -1.479, -1.305, -1.131, -0.957, -0.783, -0.609, -0.435, -0.261, -0.087,
338 0.087, 0.261, 0.435, 0.609, 0.783, 0.957, 1.131, 1.305, 1.479, 1.653, 1.830, 2.043, 2.322, 2.650};
339 float eta_limits_hcal_total_f[nbins_eta_hcal_total];
340 float log_bin_spacing =
log(200.) / 40.;
341 const int nbins_pt_total_hcal = 41;
342 double pt_limits_hcal[nbins_pt_total_hcal];
343 float pt_limits_hcal_f[nbins_pt_total_hcal];
344 for (
int i = 0;
i < nbins_pt_total_hcal;
i++) {
345 pt_limits_hcal[
i] =
exp(
i * log_bin_spacing);
346 pt_limits_hcal_f[
i] =
exp(
i * log_bin_spacing);
348 for (
int i = 0;
i < nbins_eta_hcal_total;
i++) {
349 eta_limits_hcal_total_f[
i] = (
float)eta_limits_hcal_total[
i];
352 "HOverTrackP_trackPtVsEta",
353 nbins_pt_total_hcal - 1,
355 nbins_eta_hcal_total - 1,
356 eta_limits_hcal_total_f);
358 "HOverTrackPVsTrackP_Barrel",
"HOverTrackPVsTrackP_Barrel", nbins_pt_total_hcal - 1, pt_limits_hcal, 0, 4,
" ");
360 "HOverTrackPVsTrackP_EndCap",
"HOverTrackPVsTrackP_EndCap", nbins_pt_total_hcal - 1, pt_limits_hcal, 0, 4,
" ");
362 "HOverTrackPVsTrackPt_Barrel",
363 nbins_pt_total_hcal - 1,
369 "HOverTrackPVsTrackPt_EndCap",
370 nbins_pt_total_hcal - 1,
377 "HOverTrackPVsEta, 1<hPt<10 GeV",
378 nbins_eta_hcal_total - 1,
379 eta_limits_hcal_total,
384 "HOverTrackPVsEta, 10<hPt<20 GeV",
385 nbins_eta_hcal_total - 1,
386 eta_limits_hcal_total,
391 "HOverTrackPVsEta, 20<hPt<50 GeV",
392 nbins_eta_hcal_total - 1,
393 eta_limits_hcal_total,
398 "HOverTrackPVsEta, hPt>50 GeV",
399 nbins_eta_hcal_total - 1,
400 eta_limits_hcal_total,
406 ibooker.book1D(
"HOverTrackP_Barrel_hPt_1_10",
"HOverTrackP_B, 1<hPt<10 GeV", 50, 0, 4);
408 ibooker.book1D(
"HOverTrackP_Barrel_hPt_10_20",
"HOverTrackP_B, 10<hPt<20 GeV", 50, 0, 4);
410 ibooker.book1D(
"HOverTrackP_Barrel_hPt_20_50",
"HOverTrackP_B, 20<hPt<50 GeV", 50, 0, 4);
414 ibooker.book1D(
"HOverTrackP_EndCap_hPt_1_10",
"HOverTrackP_E, 1<hPt<10 GeV", 50, 0, 4);
416 ibooker.book1D(
"HOverTrackP_EndCap_hPt_10_20",
"HOverTrackP_E, 10<hPt<20 GeV", 50, 0, 4);
418 ibooker.book1D(
"HOverTrackP_EndCap_hPt_20_50",
"HOverTrackP_E, 20<hPt<50 GeV", 50, 0, 4);
421 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackP_trackPtVsEta",
423 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackPVsTrackP_Barrel",
425 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackPVsTrackP_EndCap",
427 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackPVsTrackPt_Barrel",
429 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackPVsTrackPt_EndCap",
431 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackPVsEta_hPt_1_10",
433 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackPVsEta_hPt_10_20",
435 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackPVsEta_hPt_20_50",
439 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackP_Barrel_hPt_1_10",
441 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackP_Barrel_hPt_10_20",
443 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackP_Barrel_hPt_20_50",
445 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackP_Barrel_hPt_50",
447 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackP_EndCap_hPt_1_10",
449 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackP_EndCap_hPt_10_20",
451 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackP_EndCap_hPt_20_50",
453 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
DirName +
"/" +
"HOverTrackP_EndCap_hPt_50",
487 int etaNBinsPFCand =
v->getParameter<
int>(
"etaNBins");
488 double etaMinPFCand =
v->getParameter<
double>(
"etaMin");
489 double etaMaxPFCand =
v->getParameter<
double>(
"etaMax");
490 int phiNBinsPFCand =
v->getParameter<
int>(
"phiNBins");
491 double phiMinPFCand =
v->getParameter<
double>(
"phiMin");
492 double phiMaxPFCand =
v->getParameter<
double>(
"phiMax");
493 int nMinPFCand =
v->getParameter<
int>(
"nMin");
494 int nMaxPFCand =
v->getParameter<
int>(
"nMax");
495 int nbinsPFCand =
v->getParameter<
double>(
"nbins");
510 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
525 if (etaMinPFCand * etaMaxPFCand < 0) {
563 const int nbins_eta_endcap = 2 * (etaNBinsPFCand + 1);
564 double eta_limits_endcap[nbins_eta_endcap];
565 for (
int i = 0;
i < nbins_eta_endcap;
i++) {
566 if (
i < (etaNBinsPFCand + 1)) {
567 eta_limits_endcap[
i] = etaMinPFCand +
i * (etaMaxPFCand - etaMinPFCand) / (
double)etaNBinsPFCand;
569 eta_limits_endcap[
i] =
570 -etaMaxPFCand + (
i - (etaNBinsPFCand + 1)) * (etaMaxPFCand - etaMinPFCand) / (
double)etaNBinsPFCand;
575 nbins_eta_endcap - 1,
584 nbins_eta_endcap - 1,
590 TH2F* hist_temp_occup_puppiNolepWeight =
new TH2F(
593 nbins_eta_endcap - 1,
600 hist_temp_occup_puppiNolepWeight));
601 TH2F* hist_temp_pt_puppiNolepWeight =
604 nbins_eta_endcap - 1,
611 hist_temp_pt_puppiNolepWeight));
613 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
616 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
619 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
622 map_of_MEs.insert(std::pair<std::string, MonitorElement*>(