|
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
|
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DQMServices/Core/interface/DQMEDHarvester.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
class DQMHcalIsoTrackPostProcessor : public DQMEDHarvester {
public:
DQMHcalIsoTrackPostProcessor(const edm::ParameterSet &pset);
~DQMHcalIsoTrackPostProcessor() override {}
// void analyze(const edm::Event& event, const edm::EventSetup& eventSetup)
// override {};
void dqmEndJob(DQMStore::IBooker &,
DQMStore::IGetter &) override; // performed in the endJob
private:
std::string subDir_;
};
DQMHcalIsoTrackPostProcessor::DQMHcalIsoTrackPostProcessor(const edm::ParameterSet &pset) {
subDir_ = pset.getUntrackedParameter<std::string>("subDir");
}
void DQMHcalIsoTrackPostProcessor::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) {
if (igetter.dirExists(subDir_)) {
igetter.cd(subDir_);
} else {
edm::LogWarning("DQMHcalIsoTrackPostProcessor") << "cannot find directory: " << subDir_ << " , skipping";
return;
}
MonitorElement *hSumEta[4], *hSumPhi[4], *hPurityEta[3], *hPurityPhi[3];
std::string types[4] = {"L2", "L2x", "L3", "Off"};
char name[100], title[200];
for (int i = 0; i < 4; ++i) {
sprintf(name, "/heta%s", types[i].c_str());
std::string hname1 = ibooker.pwd() + std::string(name);
int nbinEta = igetter.get(hname1)->getTH1F()->GetNbinsX();
double xminEta = igetter.get(hname1)->getTH1F()->GetBinLowEdge(1);
double xmaxEta =
igetter.get(hname1)->getTH1F()->GetBinLowEdge(nbinEta) + igetter.get(hname1)->getTH1F()->GetBinWidth(nbinEta);
sprintf(name, "/hphi%s", types[i].c_str());
std::string hname2 = ibooker.pwd() + std::string(name);
int nbinPhi = igetter.get(hname2)->getTH1F()->GetNbinsX();
double xminPhi = igetter.get(hname2)->getTH1F()->GetBinLowEdge(1);
double xmaxPhi =
igetter.get(hname2)->getTH1F()->GetBinLowEdge(nbinEta) + igetter.get(hname2)->getTH1F()->GetBinWidth(nbinEta);
sprintf(name, "hSum%sEta", types[i].c_str());
hSumEta[i] = ibooker.book1D(name, name, nbinEta, xminEta, xmaxEta);
sprintf(name, "hSum%sPhi", types[i].c_str());
hSumPhi[i] = ibooker.book1D(name, name, nbinPhi, xminPhi, xmaxPhi);
if (i < 3) {
sprintf(name, "hPurity%sEta", types[i].c_str());
sprintf(title, "Purity of %s sample vs #eta", types[i].c_str());
hPurityEta[i] = ibooker.book1D(name, title, nbinEta, xminEta, xmaxEta);
sprintf(name, "hPurity%sPhi", types[i].c_str());
sprintf(title, "Purity of %s sample vs #phi", types[i].c_str());
hPurityPhi[i] = ibooker.book1D(name, title, nbinPhi, xminPhi, xmaxPhi);
}
}
for (int i = 0; i < 4; ++i) {
sprintf(name, "/heta%s", types[i].c_str());
std::string hname1 = ibooker.pwd() + std::string(name);
edm::LogVerbatim("DQMHcal") << "PostProcesor " << hname1 << " " << igetter.get(hname1);
hSumEta[i]->getTH1F()->Add(igetter.get(hname1)->getTH1F(), 1);
sprintf(name, "/hphi%s", types[i].c_str());
std::string hname2 = ibooker.pwd() + std::string(name);
edm::LogVerbatim("DQMHcal") << "PostProcesor " << hname2 << " " << igetter.get(hname2);
hSumPhi[i]->getTH1F()->Add(igetter.get(hname2)->getTH1F(), 1);
}
for (int i = 0; i < 3; ++i) {
hPurityEta[i]->getTH1F()->Divide(hSumEta[i + 1]->getTH1F(), hSumEta[i]->getTH1F(), 1, 1);
hPurityPhi[i]->getTH1F()->Divide(hSumPhi[i + 1]->getTH1F(), hSumPhi[i]->getTH1F(), 1, 1);
}
}
DEFINE_FWK_MODULE(DQMHcalIsoTrackPostProcessor);
|