|
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
106
|
// Migrated to use DQMEDHarvester by: Jyothsna Rani Komaragiri, Oct 2014
#include "DQMOffline/Trigger/interface/JetMETHLTOfflineClient.h"
#include "FWCore/Framework/interface/Run.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/DQMStore.h"
JetMETHLTOfflineClient::JetMETHLTOfflineClient(const edm::ParameterSet& iConfig) : conf_(iConfig) {
debug_ = false;
verbose_ = false;
processname_ = iConfig.getParameter<std::string>("processname");
hltTag_ = iConfig.getParameter<std::string>("hltTag");
if (debug_)
std::cout << hltTag_ << std::endl;
dirName_ = iConfig.getParameter<std::string>("DQMDirName");
}
JetMETHLTOfflineClient::~JetMETHLTOfflineClient() = default;
void JetMETHLTOfflineClient::dqmEndJob(DQMStore::IBooker& ibooker, DQMStore::IGetter& igetter) {
ibooker.setCurrentFolder(dirName_);
LogDebug("JetMETHLTOfflineClient") << "dqmEndJob" << std::endl;
if (debug_)
std::cout << "dqmEndJob" << std::endl;
std::vector<MonitorElement*> hltMEs;
// Look at all folders, go to the subfolder which includes the string "Eff"
std::vector<std::string> fullPathHLTFolders = igetter.getSubdirs();
for (auto& fullPathHLTFolder : fullPathHLTFolders) {
// Move on only if the folder name contains "Eff" Or "Trigger Summary"
if (debug_)
std::cout << fullPathHLTFolder << std::endl;
if ((fullPathHLTFolder.find("Eff") != std::string::npos)) {
ibooker.setCurrentFolder(fullPathHLTFolder);
} else {
continue;
}
// Look at all subfolders, go to the subfolder which includes the string "Eff"
std::vector<std::string> fullSubPathHLTFolders = igetter.getSubdirs();
for (auto& fullSubPathHLTFolder : fullSubPathHLTFolders) {
if (debug_)
std::cout << fullSubPathHLTFolder << std::endl;
ibooker.setCurrentFolder(fullSubPathHLTFolder);
// Look at all MonitorElements in this folder
hltMEs = igetter.getContents(fullSubPathHLTFolder);
LogDebug("JetMETHLTOfflineClient") << "Number of MEs for this HLT path = " << hltMEs.size() << std::endl;
for (unsigned int k = 0; k < hltMEs.size(); k++) {
if (debug_)
std::cout << hltMEs[k]->getName() << std::endl;
//-----
if ((hltMEs[k]->getName().find("ME_Numerator") != std::string::npos) &&
hltMEs[k]->getName().find("ME_Numerator") == 0) {
std::string name = hltMEs[k]->getName();
name.erase(0, 12); // Removed "ME_Numerator"
if (debug_)
std::cout << "==name==" << name << std::endl;
// if( name.find("EtaPhi") !=std::string::npos ) continue; // do not consider EtaPhi 2D plots
for (unsigned int l = 0; l < hltMEs.size(); l++) {
if (hltMEs[l]->getName() == "ME_Denominator" + name) {
// found denominator too
if (name.find("EtaPhi") != std::string::npos) {
TH2F* tNumerator = hltMEs[k]->getTH2F();
TH2F* tDenominator = hltMEs[l]->getTH2F();
std::string title = "Eff_" + hltMEs[k]->getTitle();
auto* teff = (TH2F*)tNumerator->Clone(title.c_str());
teff->Divide(tNumerator, tDenominator, 1, 1);
ibooker.book2D("ME_Eff_" + name, teff);
delete teff;
} else {
TH1F* tNumerator = hltMEs[k]->getTH1F();
TH1F* tDenominator = hltMEs[l]->getTH1F();
std::string title = "Eff_" + hltMEs[k]->getTitle();
auto* teff = (TH1F*)tNumerator->Clone(title.c_str());
teff->Divide(tNumerator, tDenominator, 1, 1);
ibooker.book1D("ME_Eff_" + name, teff);
delete teff;
}
} // Denominator
} // Loop-l
} // Numerator
} // Loop-k
} // fullSubPath
} // fullPath
}
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(JetMETHLTOfflineClient);
|