|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
|
#include "DQMOffline/Trigger/plugins/HMesonGammaDQM.h"
#include "DataFormats/Math/interface/deltaR.h"
HMesonGammaDQM::HMesonGammaDQM() = default;
HMesonGammaDQM::~HMesonGammaDQM() = default;
void HMesonGammaDQM::initialise(const edm::ParameterSet& iConfig) {
gammapt_variable_binning_ =
iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("gammaptBinning");
mesonpt_variable_binning_ =
iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("mesonptBinning");
eta_binning_ =
getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("hmgetaPSet"));
ls_binning_ =
getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("hmglsPSet"));
}
void HMesonGammaDQM::bookHistograms(DQMStore::IBooker& ibooker) {
std::string histname, histtitle;
histname = "gammapt";
histtitle = "Gamma pT";
bookME(ibooker, gammaptME_, histname, histtitle, gammapt_variable_binning_);
setMETitle(gammaptME_, "Gamma p_{T} [GeV]", "events / [GeV]");
histname = "mesonpt";
histtitle = "Meson pT";
bookME(ibooker, mesonptME_, histname, histtitle, mesonpt_variable_binning_);
setMETitle(mesonptME_, "Meson p_{T} [GeV]", "events / [GeV]");
histname = "gammaeta";
histtitle = "Gamma eta";
bookME(ibooker, gammaetaME_, histname, histtitle, eta_binning_.nbins, eta_binning_.xmin, eta_binning_.xmax);
setMETitle(gammaetaME_, "Gamma #eta", "events");
histname = "mesoneta";
histtitle = "Meson eta";
bookME(ibooker, mesonetaME_, histname, histtitle, eta_binning_.nbins, eta_binning_.xmin, eta_binning_.xmax);
setMETitle(mesonetaME_, "Meson #eta", "events");
histname = "gammaetaVsLS";
histtitle = "Gamma eta vs LS";
bookME(ibooker,
gammaetaVsLS_,
histname,
histtitle,
ls_binning_.nbins,
ls_binning_.xmin,
ls_binning_.xmax,
eta_binning_.xmin,
eta_binning_.xmax);
setMETitle(gammaetaVsLS_, "LS", "Gamma #eta");
}
void HMesonGammaDQM::fillHistograms(const reco::PhotonCollection& photons,
const std::vector<TLorentzVector>& mesons,
const int ls,
const bool passCond) {
// filling histograms (denominator)
if (!photons.empty()) {
double eta1 = photons[0].eta();
gammaptME_.denominator->Fill(photons[0].pt());
gammaetaME_.denominator->Fill(eta1);
gammaetaVsLS_.denominator->Fill(ls, eta1);
}
if (!mesons.empty()) {
double eta2 = mesons[0].Eta();
mesonptME_.denominator->Fill(mesons[0].Pt());
mesonetaME_.denominator->Fill(eta2);
}
// applying selection for numerator
if (passCond) {
if (!photons.empty()) {
double eta1 = photons[0].eta();
gammaptME_.numerator->Fill(photons[0].pt());
gammaetaME_.numerator->Fill(eta1);
gammaetaVsLS_.numerator->Fill(ls, eta1);
}
if (!mesons.empty()) {
double eta2 = mesons[0].Eta();
mesonptME_.numerator->Fill(mesons[0].Pt());
mesonetaME_.numerator->Fill(eta2);
}
}
}
void HMesonGammaDQM::fillHmgDescription(edm::ParameterSetDescription& histoPSet) {
edm::ParameterSetDescription hmgetaPSet;
fillHistoPSetDescription(hmgetaPSet);
histoPSet.add<edm::ParameterSetDescription>("hmgetaPSet", hmgetaPSet);
std::vector<double> pt1bins = {0., 20., 40., 60., 80., 90., 100., 110., 120., 130., 140., 150., 160.,
180., 210., 240., 270., 300., 330., 360., 400., 450., 500., 750., 1000., 1500.};
histoPSet.add<std::vector<double> >("gammaptBinning", pt1bins);
std::vector<double> pt2bins = {0., 20., 40., 45., 50., 55., 60., 65., 70., 80., 90., 100.,
110., 120., 150., 180., 210., 240., 270., 300., 350., 400., 1000.};
histoPSet.add<std::vector<double> >("mesonptBinning", pt2bins);
edm::ParameterSetDescription lsPSet;
fillHistoLSPSetDescription(lsPSet);
histoPSet.add<edm::ParameterSetDescription>("hmglsPSet", lsPSet);
}
|