|
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
107
108
109
110
|
#include "CondFormats/PPSObjects/interface/PPSDirectSimulationData.h"
#include "DataFormats/CTPPSDetId/interface/CTPPSDetId.h"
#include "DataFormats/CTPPSDetId/interface/TotemRPDetId.h"
#include "DataFormats/CTPPSDetId/interface/CTPPSPixelDetId.h"
#include "DataFormats/CTPPSDetId/interface/CTPPSDiamondDetId.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "FWCore/Utilities/interface/FileInPath.h"
#include "TFile.h"
#include <iostream>
#include <string>
PPSDirectSimulationData::PPSDirectSimulationData()
: empiricalAperture45_(""),
empiricalAperture56_(""),
timeResolutionDiamonds45_(""),
timeResolutionDiamonds56_("") {}
PPSDirectSimulationData::~PPSDirectSimulationData() {}
// Getters
const std::string &PPSDirectSimulationData::getEmpiricalAperture45() const { return empiricalAperture45_; }
const std::string &PPSDirectSimulationData::getEmpiricalAperture56() const { return empiricalAperture56_; }
const std::string &PPSDirectSimulationData::getTimeResolutionDiamonds45() const { return timeResolutionDiamonds45_; }
const std::string &PPSDirectSimulationData::getTimeResolutionDiamonds56() const { return timeResolutionDiamonds56_; }
std::map<unsigned int, PPSDirectSimulationData::FileObject> &PPSDirectSimulationData::getEfficienciesPerRP() {
return efficienciesPerRP_;
}
std::map<unsigned int, PPSDirectSimulationData::FileObject> &PPSDirectSimulationData::getEfficienciesPerPlane() {
return efficienciesPerPlane_;
};
// Setters
void PPSDirectSimulationData::setEmpiricalAperture45(std::string s) { empiricalAperture45_ = s; }
void PPSDirectSimulationData::setEmpiricalAperture56(std::string s) { empiricalAperture56_ = s; }
void PPSDirectSimulationData::setTimeResolutionDiamonds45(std::string s) { timeResolutionDiamonds45_ = s; }
void PPSDirectSimulationData::setTimeResolutionDiamonds56(std::string s) { timeResolutionDiamonds56_ = s; }
std::map<unsigned int, std::unique_ptr<TH2F>> PPSDirectSimulationData::loadEffeciencyHistogramsPerRP() const {
std::map<unsigned int, std::unique_ptr<TH2F>> result;
for (const auto &it : efficienciesPerRP_)
result[it.first] = loadObject(it.second.first, it.second.second);
return result;
}
std::map<unsigned int, std::unique_ptr<TH2F>> PPSDirectSimulationData::loadEffeciencyHistogramsPerPlane() const {
std::map<unsigned int, std::unique_ptr<TH2F>> result;
for (const auto &it : efficienciesPerPlane_) {
CTPPSDetId rpId(it.first);
if (rpId.subdetId() == CTPPSDetId::sdTrackingStrip) {
for (unsigned int pl = 0; pl < 10; ++pl) {
TotemRPDetId plId(rpId.arm(), rpId.station(), rpId.rp(), pl);
result[plId] = loadObject(it.second.first, replace(it.second.second, "<detid>", std::to_string(pl)));
}
}
if (rpId.subdetId() == CTPPSDetId::sdTrackingPixel) {
for (unsigned int pl = 0; pl < 6; ++pl) {
CTPPSPixelDetId plId(rpId.arm(), rpId.station(), rpId.rp(), pl);
result[plId] = loadObject(it.second.first, replace(it.second.second, "<detid>", std::to_string(pl)));
}
}
if (rpId.subdetId() == CTPPSDetId::sdTimingDiamond) {
for (unsigned int pl = 0; pl < 4; ++pl) {
CTPPSDiamondDetId plId(rpId.arm(), rpId.station(), rpId.rp(), pl);
result[plId] = loadObject(it.second.first, replace(it.second.second, "<detid>", std::to_string(pl)));
}
}
}
return result;
}
std::unique_ptr<TH2F> PPSDirectSimulationData::loadObject(const std::string &file, const std::string &object) {
edm::FileInPath fip(file.c_str());
TFile *f_in = TFile::Open(fip.fullPath().c_str());
if (!f_in)
throw cms::Exception("PPS") << "Cannot open file '" << fip.fullPath() << "'.";
TH2F *o_in = (TH2F *)f_in->Get(object.c_str());
if (!o_in)
throw cms::Exception("PPS") << "Cannot load object '" << object << "' from file '" << fip.fullPath() << "'.";
// disassociate histogram from the file
o_in->SetDirectory(nullptr);
delete f_in;
return std::unique_ptr<TH2F>(o_in);
}
std::string PPSDirectSimulationData::replace(std::string input, const std::string &from, const std::string &to) {
size_t start_pos = 0;
while ((start_pos = input.find(from, start_pos)) != std::string::npos) {
input.replace(start_pos, from.length(), to);
start_pos += to.length();
}
return input;
}
|