|
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
|
#include "DQM/GEM/plugins/GEMEfficiencyHarvester.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/isFinite.h"
GEMEfficiencyHarvester::GEMEfficiencyHarvester(const edm::ParameterSet& ps)
: GEMDQMEfficiencyClientBase(ps), kFolders_(ps.getUntrackedParameter<std::vector<std::string> >("folders")) {}
GEMEfficiencyHarvester::~GEMEfficiencyHarvester() {}
void GEMEfficiencyHarvester::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
edm::ParameterSetDescription desc;
// GEMDQMEfficiencyClientBase
desc.addUntracked<double>("confidenceLevel", 0.683); // 1-sigma
desc.addUntracked<std::string>("logCategory", "GEMEfficiencyHarvester");
// GEMEfficiencyHarvester
desc.addUntracked<std::vector<std::string> >("folders", {"GEM/Efficiency/muonSTA"});
descriptions.add("gemEfficiencyHarvester", desc);
}
// boook MEs for
// - efficiency vs camber id
// - efficiency vs ieta
// by projecting MEs used for 2d chamber-ieta efficiency.
void GEMEfficiencyHarvester::bookDetector1DEfficiency(DQMStore::IBooker& ibooker,
DQMStore::IGetter& igetter,
const std::string& folder) {
const std::map<std::string, MEPair> me_pairs = makeEfficiencySourcePair(ibooker, igetter, folder, "chamber_ieta_");
for (const auto& [key, value] : me_pairs) {
// numerator and denominator
const auto& [me_num, me_den] = value;
if (me_num->kind() != MonitorElement::Kind::TH2F) {
edm::LogError(kLogCategory_) << key << "expected TH2F but got "; // TODO
continue;
}
const TH2F* h_num = me_num->getTH2F();
if (h_num == nullptr) {
edm::LogError(kLogCategory_) << "numerator: failed to get TH2F from MonitorElement" << key;
continue;
}
const TH2F* h_den = me_den->getTH2F();
if (h_den == nullptr) {
edm::LogError(kLogCategory_) << "denominator: failed to get TH2F from MonitorElement" << key;
continue;
}
const auto [parsing_success, var_name, gem_name, is_matched] = parseEfficiencySourceName(me_den->getName());
if (not parsing_success) {
edm::LogError(kLogCategory_) << "failed to parse " << me_den->getName();
continue;
}
// TODO sanity-check
const TH1F* h_chamber_num = projectHistogram(h_num, TH1::kXaxis);
const TH1F* h_chamber_den = projectHistogram(h_den, TH1::kXaxis);
const char* eff_chamber_name = Form("eff_chamber_%s", gem_name.c_str());
if (TH1F* eff = makeEfficiency(h_chamber_num, h_chamber_den, eff_chamber_name)) {
ibooker.book1D(eff_chamber_name, eff);
}
const TH1F* h_ieta_num = projectHistogram(h_num, TH1::kYaxis);
const TH1F* h_ieta_den = projectHistogram(h_den, TH1::kYaxis);
const char* eff_ieta_name = Form("eff_ieta_%s", gem_name.c_str());
if (TH1F* eff = makeEfficiency(h_ieta_num, h_ieta_den, eff_ieta_name)) {
ibooker.book1D(eff_ieta_name, eff);
}
} // pairs
}
void GEMEfficiencyHarvester::dqmEndJob(DQMStore::IBooker& ibooker, DQMStore::IGetter& igetter) {
for (const std::string& folder : kFolders_) {
bookEfficiencyAuto(ibooker, igetter, folder);
bookDetector1DEfficiency(ibooker, igetter, folder);
}
}
|