Belle II Software development
TOPPhotonYieldsAlgorithm.cc
1/**************************************************************************
2 * basf2 (Belle II Analysis Software Framework) *
3 * Author: The Belle II Collaboration *
4 * *
5 * See git log for contributors and copyright holders. *
6 * This file is licensed under LGPL-3.0, see LICENSE.md. *
7 **************************************************************************/
8
9#include <top/calibration/TOPPhotonYieldsAlgorithm.h>
10#include <top/geometry/TOPGeometryPar.h>
11#include <top/dbobjects/TOPCalPhotonYields.h>
12#include <top/dbobjects/TOPCalChannelRQE.h>
13#include <TROOT.h>
14#include <TFile.h>
15#include <TProfile.h>
16#include <string>
17
18using namespace std;
19
20namespace Belle2 {
25 namespace TOP {
26
27 TOPPhotonYieldsAlgorithm::TOPPhotonYieldsAlgorithm():
28 CalibrationAlgorithm("TOPPhotonYieldsCollector")
29 {
30 setDescription("Calibration algorithm for photon pixel yields aimed for PMT ageing studies, for RQE calibration "
31 "and for finding optically decoupled PMT's");
32 }
33
34
35 CalibrationAlgorithm::EResult TOPPhotonYieldsAlgorithm::calibrate()
36 {
37 gROOT->SetBatch();
38
39 // IoV for getting geometry payloads - they differ only between Run1 and Run2
40
41 const auto& expRun = getRunList();
42 updateDBObjPtrs(1, expRun[0].second, expRun[0].first);
43
44 // geometry parameters
45
46 auto* topgp = TOPGeometryPar::Instance();
47 if (not topgp->isValid()) topgp->Initialize();
48 if (not topgp->isValid()) {
49 B2ERROR("TOPPhotonYieldsAlgorithm: cannot initialize geometry parameters");
50 return c_Failure;
51 }
52
53 // construct file name and open output root file for storing merged histograms
54
55 string expNo = to_string(expRun[0].first);
56 while (expNo.length() < 4) expNo.insert(0, "0");
57 string runNo = to_string(expRun[0].second);
58 while (runNo.length() < 5) runNo.insert(0, "0");
59 string outputFileName = "photonYields-e" + expNo + "-r" + runNo + ".root";
60 auto* file = TFile::Open(outputFileName.c_str(), "recreate");
61
62 // create some additional histograms for monitoring of RQE calibration
63
64 auto* h_rqe = new TH1F("RQE", "RQE distribution; RQE; entries/bin", 100, 0, m_maxRQE);
65 auto* h_err = new TH1F("RQEerr", "RQE uncertainty distribution; RQE uncertainty; entries/bin", 100, 0, m_maxErrorRQE);
66
67 // get basic histograms
68
69 auto numTracks = getObjectPtr<TH1F>("numTracks");
70 if (not numTracks) {
71 B2ERROR("TOPPhotonYieldsAlgorithm: histogram 'numTracks' not found");
72 file->Close();
73 return c_Failure;
74 }
75 numTracks->Write();
76
77 auto timeStamp = getObjectPtr<TProfile>("timeStamp");
78 if (not timeStamp) {
79 B2ERROR("TOPPhotonYieldsAlgorithm: histogram 'timeStamp' not found");
80 file->Close();
81 return c_Failure;
82 }
83 timeStamp->Write();
84
85 // construct DB objects for storing the results and set the time stamp
86
87 auto* dbPhotonYields = new TOPCalPhotonYields();
88 dbPhotonYields->setTimeStamp(timeStamp->GetBinContent(1), timeStamp->GetBinError(1));
89 auto* dbRQE = new TOPCalChannelRQE();
90
91 // for each slot determine photon pixel yields, equalized alpha ratio and RQE, and store them in DB objects
92
93 int numModules = numTracks->GetNbinsX();
94 for (int slot = 1; slot <= numModules; slot++) {
95 string slotName = (slot < 10) ? "_0" + to_string(slot) : "_" + to_string(slot);
96
97 auto signalHits = getObjectPtr<TH1F>("signalHits" + slotName);
98 if (not signalHits) {
99 B2ERROR("TOPPhotonYieldsAlgorithm: histogram 'signalHits' for slot " << slot << " not found");
100 file->Close();
101 delete dbPhotonYields;
102 delete dbRQE;
103 return c_Failure;
104 }
105 signalHits->Write();
106
107 auto bkgHits = getObjectPtr<TH1F>("bkgHits" + slotName);
108 if (not bkgHits) {
109 B2ERROR("TOPPhotonYieldsAlgorithm: histogram 'bkgHits' for slot " << slot << " not found");
110 file->Close();
111 delete dbPhotonYields;
112 delete dbRQE;
113 return c_Failure;
114 }
115 bkgHits->Write();
116
117 auto effectiveSignalHits = getObjectPtr<TH1F>("effectiveSignalHits" + slotName);
118 if (not effectiveSignalHits) {
119 B2ERROR("TOPPhotonYieldsAlgorithm: histogram 'effectiveSignalHits' for slot " << slot << " not found");
120 file->Close();
121 delete dbPhotonYields;
122 delete dbRQE;
123 return c_Failure;
124 }
125 effectiveSignalHits->Write();
126
127 auto activePixels = getObjectPtr<TH1F>("activePixels" + slotName);
128 if (not activePixels) {
129 B2ERROR("TOPPhotonYieldsAlgorithm: histogram 'activePixels' for slot " << slot << " not found");
130 file->Close();
131 delete dbPhotonYields;
132 delete dbRQE;
133 return c_Failure;
134 }
135 activePixels->Write();
136
137 auto alphaLow = getObjectPtr<TH1F>("alphaLow" + slotName);
138 if (not alphaLow) {
139 B2ERROR("TOPPhotonYieldsAlgorithm: histogram 'alphaLow' for slot " << slot << " not found");
140 file->Close();
141 delete dbPhotonYields;
142 delete dbRQE;
143 return c_Failure;
144 }
145 alphaLow->Sumw2();
146 alphaLow->Write();
147
148 auto alphaHigh = getObjectPtr<TH1F>("alphaHigh" + slotName);
149 if (not alphaHigh) {
150 B2ERROR("TOPPhotonYieldsAlgorithm: histogram 'alphaHigh' for slot " << slot << " not found");
151 file->Close();
152 delete dbPhotonYields;
153 delete dbRQE;
154 return c_Failure;
155 }
156 alphaHigh->Sumw2();
157 alphaHigh->Write();
158
159 auto pulseHeights = getObjectPtr<TH2F>("pulseHeights" + slotName);
160 if (not pulseHeights) {
161 B2ERROR("TOPPhotonYieldsAlgorithm: histogram 'pulseHeights' for slot " << slot << " not found");
162 file->Close();
163 delete dbPhotonYields;
164 delete dbRQE;
165 return c_Failure;
166 }
167 pulseHeights->Write();
168
169 auto muonZ = getObjectPtr<TH1F>("muonZ" + slotName);
170 if (not muonZ) {
171 B2ERROR("TOPPhotonYieldsAlgorithm: histogram 'muonZ' for slot " << slot << " not found");
172 file->Close();
173 delete dbPhotonYields;
174 delete dbRQE;
175 return c_Failure;
176 }
177 muonZ->Write();
178
179 auto* photonYields = (TH1F*) signalHits->Clone("tmp1");
180 photonYields->Add(signalHits.get(), bkgHits.get(), 1, -1); // subtract background
181 for (int bin = 1; bin <= activePixels->GetNbinsX(); bin++) activePixels->SetBinError(bin, 0);
182 photonYields->Divide(activePixels.get()); // normalize
183
184 auto* bkgYields = (TH1F*) bkgHits->Clone("tmp2");
185 bkgYields->Divide(activePixels.get()); // normalize
186
187 auto* alphaRatio = (TH1F*) alphaHigh->Clone("tmp2");
188 alphaRatio->Divide(alphaHigh.get(), alphaLow.get());
189 equalize(alphaRatio);
190
191 dbPhotonYields->set(slot, photonYields, bkgYields, alphaRatio, activePixels.get(), pulseHeights.get(), muonZ.get());
192
193 effectiveSignalHits->Sumw2();
194 effectiveSignalHits->Divide(activePixels.get()); // normalize
195 const auto& channelMapper = topgp->getChannelMapper();
196 for (int row = 0; row < 8; row++) {
197 double muonCorr = getMuonCorrection(muonZ.get(), row);
198 for (int col = 0; col < 64; col++) {
199 int pixelID = row * 64 + col + 1;
200 double mcYield = getNominalYield(slot, row, col) * topgp->getRelativePDEonMC(slot, pixelID) * muonCorr;
201 if (mcYield == 0) continue; // this should not happen
202 double rqe = effectiveSignalHits->GetBinContent(pixelID) / mcYield;
203 double err = effectiveSignalHits->GetBinError(pixelID) / mcYield;
204 if (rqe == 0) continue; // pixel dead or masked-out
205 h_rqe->Fill(rqe);
206 h_err->Fill(err);
207 if (rqe > m_maxRQE or err > m_maxErrorRQE) continue; // bad calibration, do not save it in payload
208 unsigned channel = channelMapper.getChannel(pixelID);
209 dbRQE->setRQE(slot, channel, rqe);
210 }
211 }
212
213 delete photonYields;
214 delete bkgYields;
215 delete alphaRatio;
216 }
217
218 h_rqe->Write();
219 h_err->Write();
220 file->Close();
221
222 saveCalibration(dbPhotonYields);
223 saveCalibration(dbRQE);
224
225 return c_OK;
226 }
227
228
229 double TOPPhotonYieldsAlgorithm::getEqualizingValue(int bin) const
230 {
231 // these constants are determined from fits to MC
232 const double a[] = {1.50301, 1.16907, 1.00912, 0.943571, 0.889658, 0.93418, 1.01846, 0.953715};
233 const double b[] = {-1.02531, -0.111269, 0.0731595, -0.121119, -0.79083, -1.00261, -1.41012, -1.23368};
234 const double c[] = {8.47995, 3.18639, 1.9513, 2.58279, 4.79672, 5.5659, 7.558, 6.71336};
235 const double d[] = {1.49743, 1.17224, 1.02211, 0.928502, 0.829334, 0.849355, 0.921871, 0.850427};
236
237 int row = (bin - 1) / 64;
238 if (row < 0 or row > 7) return 1;
239
240 int col = (bin - 1) % 64 + 1;
241 if (col == 1 or col == 64) return d[row];
242
243 double x = (col - 32.5) / 64.;
244 return a[row] + b[row] * pow(x, 2) + c[row] * pow(x, 4);
245 }
246
247
248 void TOPPhotonYieldsAlgorithm::equalize(TH1F* h) const
249 {
250 for (int bin = 1; bin <= h->GetNbinsX(); bin++) {
251 double s = getEqualizingValue(bin);
252 h->SetBinContent(bin, h->GetBinContent(bin) / s);
253 h->SetBinError(bin, h->GetBinError(bin) / s);
254 }
255 }
256
257
258 double TOPPhotonYieldsAlgorithm::getNominalYield(int slot, int row, int col) const
259 {
260 // these constants are determined from fits to MC
261 const double par[8][4] = {{0.0404913, -2.10514e-07, -1.01752e-08, 1.00656e-11},
262 {0.0406668, -3.045e-06, -3.77908e-09, 6.4395e-12},
263 {0.041261, -1.18e-06, -7.12563e-09, 8.68068e-12},
264 {0.0439403, -1.38318e-06, -4.55235e-09, 6.7621e-12},
265 {0.04622, 2.92406e-06, -1.4338e-08, 1.43519e-11},
266 {0.0469347, 1.30621e-06, -1.10208e-08, 1.23847e-11},
267 {0.0476307, 7.59058e-07, -1.1142e-08, 1.28852e-11},
268 {0.0410036, 1.2299e-06, -1.26732e-08, 1.32523e-11}
269 };
270 const double slot_sf[16] = {1.0296, 1.00401, 0.980722, 0.997997, 1.01303, 1.00586, 0.993214, 0.978076,
271 0.981156, 0.992147, 0.999084, 0.964364, 1.00546, 1.02583, 1.01906, 1.01039
272 };
273
274 if (slot < 1 or slot > 16) return 0;
275 if (row < 0 or row > 7) return 0;
276
277 double x = col - 31.5;
278 x *= x;
279 const double* p = par[row];
280
281 return (p[0] + p[1] * x + p[2] * x * x + p[3] * x * x * x) * slot_sf[slot - 1];
282 }
283
284
285 double TOPPhotonYieldsAlgorithm::getMuonCorrection(const TH1F* h_mu, int row) const
286 {
287 // these constants are determined from MC
288 const double numFot[8][100] = {{
289 2.35, 3.41, 5.22, 5.57, 4.26, 4.47, 6.13, 5.01, 4.93, 6.35, 5.07, 5.71, 4.88, 4.98, 4.56, 4.06,
290 3.86, 3.38, 2.69, 2.11, 1.8, 1.44, 0.99, 0.52, 0.32, 0.26, 0.27, 0.28, 0.29, 0.36, 0.63, 1.0,
291 1.29, 1.53, 1.89, 2.38, 2.76, 3.05, 3.34, 3.63, 3.79, 4.04, 4.2, 4.33, 4.3, 4.13, 4.0, 3.86,
292 3.72, 3.62, 3.39, 3.14, 2.89, 2.68, 2.54, 2.44, 2.37, 2.31, 2.2, 2.09, 1.89, 1.88, 1.82, 1.76,
293 1.7, 1.71, 1.71, 1.74, 1.7, 1.73, 1.67, 1.62, 1.49, 1.45, 1.44, 1.45, 1.36, 1.38, 1.35, 1.4,
294 1.42, 1.46, 1.45, 1.47, 1.36, 1.24, 1.24, 1.34, 1.35, 1.33, 1.27, 1.29, 1.38, 1.53, 1.69, 1.75,
295 1.77, 1.78, 1.75, 0.66
296 },
297 {
298 2.41, 3.37, 5.34, 4.54, 3.28, 4.59, 4.64, 3.62, 4.89, 4.22, 4.52, 4.75, 4.94, 4.58, 4.28, 4.21,
299 3.54, 3.09, 2.82, 2.63, 2.46, 2.17, 1.92, 1.54, 1.07, 0.52, 0.35, 0.41, 0.66, 1.05, 1.35, 1.6,
300 1.82, 1.99, 2.12, 2.32, 2.57, 2.98, 3.25, 3.48, 3.59, 3.65, 3.55, 3.39, 3.3, 3.23, 3.2, 3.17,
301 3.18, 3.19, 3.17, 3.19, 3.17, 3.06, 2.84, 2.68, 2.5, 2.39, 2.25, 2.15, 1.97, 1.93, 1.89, 1.77,
302 1.71, 1.73, 1.73, 1.77, 1.74, 1.72, 1.68, 1.63, 1.54, 1.51, 1.51, 1.46, 1.38, 1.38, 1.4, 1.41,
303 1.45, 1.49, 1.49, 1.55, 1.46, 1.33, 1.36, 1.41, 1.49, 1.51, 1.53, 1.72, 2.04, 2.42, 2.69, 3.01,
304 3.13, 3.13, 3.02, 1.37
305 },
306 {
307 2.61, 3.9, 4.86, 3.35, 3.08, 4.39, 3.66, 3.21, 4.28, 3.04, 3.95, 3.33, 3.71, 3.47, 3.49, 3.4,
308 3.1, 2.95, 2.95, 2.88, 2.69, 2.54, 2.39, 2.22, 1.97, 1.75, 1.49, 1.38, 1.47, 1.64, 1.79, 1.91,
309 2.03, 2.13, 2.23, 2.31, 2.39, 2.47, 2.61, 2.74, 2.7, 2.64, 2.67, 2.65, 2.66, 2.74, 2.74, 2.77,
310 2.81, 2.81, 2.88, 2.93, 2.97, 3.01, 3.11, 3.14, 3.1, 2.97, 2.8, 2.54, 2.25, 2.23, 2.12, 2.01,
311 1.97, 1.93, 1.91, 1.91, 1.88, 1.87, 1.86, 1.83, 1.7, 1.63, 1.62, 1.63, 1.54, 1.56, 1.55, 1.6,
312 1.64, 1.75, 1.8, 1.96, 2.08, 2.26, 2.38, 2.59, 2.79, 2.91, 2.99, 3.14, 3.52, 3.73, 3.95, 4.1,
313 4.13, 4.01, 3.86, 1.78
314 },
315 {
316 3.16, 4.53, 3.66, 2.69, 3.28, 3.9, 3.01, 3.16, 3.17, 2.91, 3.2, 2.18, 2.72, 2.11, 2.32, 2.58,
317 2.89, 3.11, 3.31, 3.25, 3.16, 3.02, 2.89, 2.87, 3.07, 3.27, 3.27, 3.18, 2.75, 2.37, 2.22, 2.26,
318 2.35, 2.42, 2.5, 2.53, 2.36, 2.12, 1.8, 1.63, 1.65, 1.8, 1.99, 2.06, 2.21, 2.24, 2.36, 2.45,
319 2.51, 2.55, 2.64, 2.7, 2.77, 2.83, 2.88, 2.97, 3.08, 3.16, 3.24, 3.22, 3.15, 2.94, 2.73, 2.51,
320 2.38, 2.34, 2.27, 2.2, 2.22, 2.16, 2.15, 2.07, 2.02, 1.94, 1.89, 1.94, 1.95, 2.0, 2.2, 2.45,
321 2.76, 3.01, 3.28, 3.55, 3.54, 3.71, 3.68, 3.88, 3.96, 4.05, 4.08, 4.26, 4.47, 4.7, 4.74, 4.88,
322 4.86, 4.81, 4.52, 2.24
323 },
324 {
325 4.53, 2.88, 1.1, 2.47, 3.57, 2.03, 1.38, 2.58, 1.56, 1.4, 1.42, 1.52, 1.54, 1.64, 1.76, 1.85,
326 1.96, 2.12, 2.29, 2.57, 2.98, 3.48, 3.95, 4.43, 4.82, 4.95, 5.01, 4.85, 4.75, 4.21, 3.53, 2.93,
327 2.45, 2.02, 1.82, 1.6, 1.45, 1.29, 1.19, 1.07, 0.92, 0.84, 0.77, 0.73, 0.81, 1.04, 1.33, 1.55,
328 1.76, 1.88, 2.05, 2.13, 2.2, 2.35, 2.37, 2.49, 2.55, 2.67, 2.76, 2.9, 2.96, 3.05, 3.1, 3.16,
329 3.34, 3.34, 3.46, 3.49, 3.69, 3.69, 3.83, 3.91, 3.99, 4.14, 4.14, 4.22, 4.37, 4.42, 4.55, 4.55,
330 4.74, 4.9, 5.06, 5.11, 5.27, 5.28, 5.34, 5.34, 5.37, 5.56, 5.59, 5.68, 5.78, 5.87, 5.96, 6.22,
331 6.29, 6.11, 5.57, 2.84
332 },
333 {
334 3.6, 1.42, 0.89, 2.24, 2.72, 1.27, 1.14, 1.17, 1.22, 1.23, 1.34, 1.32, 1.44, 1.49, 1.59, 1.71,
335 1.79, 1.93, 2.11, 2.32, 2.58, 2.9, 3.46, 4.3, 5.16, 5.49, 5.58, 5.52, 5.22, 4.48, 3.44, 2.69,
336 2.18, 1.86, 1.6, 1.39, 1.2, 1.03, 0.92, 0.82, 0.73, 0.65, 0.62, 0.57, 0.54, 0.52, 0.54, 0.73,
337 1.04, 1.34, 1.55, 1.7, 1.84, 1.98, 2.04, 2.2, 2.31, 2.39, 2.51, 2.61, 2.71, 2.79, 2.93, 3.12,
338 3.42, 3.78, 4.14, 4.5, 4.77, 5.02, 5.23, 5.49, 5.58, 5.61, 5.56, 5.55, 5.44, 5.44, 5.42, 5.49,
339 5.54, 5.63, 5.76, 5.94, 5.99, 6.04, 5.96, 6.05, 6.15, 6.15, 6.23, 6.32, 6.42, 6.44, 6.6, 6.83,
340 6.93, 6.9, 6.15, 3.12
341 },
342 {
343 2.05, 0.84, 0.9, 1.14, 1.09, 0.99, 1.08, 1.12, 1.15, 1.21, 1.23, 1.3, 1.36, 1.42, 1.53, 1.58,
344 1.68, 1.84, 2.07, 2.3, 2.59, 2.87, 3.33, 3.96, 4.66, 5.27, 5.7, 5.58, 5.04, 4.32, 3.37, 2.71,
345 2.29, 1.9, 1.54, 1.25, 1.0, 0.86, 0.75, 0.68, 0.61, 0.55, 0.53, 0.5, 0.47, 0.45, 0.44, 0.44,
346 0.47, 0.57, 0.76, 1.04, 1.29, 1.55, 1.7, 1.82, 1.94, 2.1, 2.23, 2.49, 2.74, 3.14, 3.43, 3.86,
347 4.12, 4.48, 4.74, 4.91, 5.16, 5.35, 5.57, 5.79, 5.99, 6.16, 6.33, 6.49, 6.61, 6.69, 6.75, 6.72,
348 6.7, 6.69, 6.71, 6.64, 6.73, 6.79, 6.89, 6.75, 6.85, 6.86, 6.99, 6.98, 7.06, 7.03, 7.08, 7.35,
349 7.44, 7.31, 6.48, 3.32
350 },
351 {
352 0.54, 0.54, 0.66, 0.66, 0.68, 0.77, 0.79, 0.82, 0.92, 0.9, 0.93, 1.0, 1.01, 1.07, 1.13, 1.19,
353 1.25, 1.42, 1.54, 1.78, 2.0, 2.31, 2.7, 3.22, 3.6, 3.68, 3.64, 3.68, 3.78, 3.6, 2.83, 2.27,
354 1.83, 1.43, 1.1, 0.89, 0.72, 0.58, 0.49, 0.45, 0.43, 0.39, 0.37, 0.33, 0.34, 0.32, 0.32, 0.31,
355 0.34, 0.32, 0.33, 0.36, 0.43, 0.62, 0.83, 1.07, 1.37, 1.74, 2.14, 2.55, 2.88, 3.24, 3.47, 3.74,
356 4.01, 4.29, 4.52, 4.74, 4.91, 5.1, 5.31, 5.46, 5.62, 5.84, 5.97, 6.14, 6.26, 6.41, 6.52, 6.62,
357 6.78, 6.82, 6.98, 7.03, 7.0, 6.99, 7.0, 6.84, 6.85, 6.65, 6.59, 6.49, 6.44, 6.31, 6.12, 6.05,
358 6.05, 6.01, 5.23, 2.58
359 }
360 };
361 const double nfotMC[8] = {2.5503, 2.5504, 2.6065, 2.7892, 2.9729, 3.0088, 3.045, 2.6138};
362
363 if (row < 0 or row > 7) return 0;
364
365 if (h_mu->GetNbinsX() != 100 or h_mu->GetXaxis()->GetXmin() != -130 or h_mu->GetXaxis()->GetXmax() != 130) {
366 B2ERROR("Histogram " << h_mu->GetTitle() << ": wrong number of bins or wrong range, no yield corrections possible");
367 return 0;
368 }
369
370 const double* nFot = numFot[row];
371 double nfot = 0;
372 double suma = 0;
373 for (int i = 0; i < 100; i++) {
374 double binContent = h_mu->GetBinContent(i + 1);
375 nfot += nFot[i] * binContent;
376 suma += binContent;
377 }
378 if (suma != 0) nfot /= suma;
379
380 return nfot / nfotMC[row];
381 }
382
383
384
385 } // end namespace TOP
387} // end namespace Belle2
Belle2::CalibrationAlgorithm::saveCalibration
void saveCalibration(TClonesArray *data, const std::string &name)
Store DBArray payload with given name with default IOV.
Definition CalibrationAlgorithm.cc:297
Belle2::CalibrationAlgorithm::updateDBObjPtrs
void updateDBObjPtrs(const unsigned int event, const int run, const int experiment)
Updates any DBObjPtrs by calling update(event) for DBStore.
Definition CalibrationAlgorithm.cc:404
Belle2::CalibrationAlgorithm::setDescription
void setDescription(const std::string &description)
Set algorithm description (in constructor)
Definition CalibrationAlgorithm.h:321
Belle2::CalibrationAlgorithm::getRunList
const std::vector< Calibration::ExpRun > & getRunList() const
Get the list of runs for which calibration is called.
Definition CalibrationAlgorithm.h:266
Belle2::CalibrationAlgorithm::EResult
EResult
The result of calibration.
Definition CalibrationAlgorithm.h:40
Belle2::CalibrationAlgorithm::c_OK
@ c_OK
Finished successfully =0 in Python.
Definition CalibrationAlgorithm.h:41
Belle2::CalibrationAlgorithm::CalibrationAlgorithm
CalibrationAlgorithm(const std::string &collectorModuleName)
Constructor - sets the prefix for collected objects (won't be accesses until execute(....
Definition CalibrationAlgorithm.h:140
Belle2::TOPCalChannelRQE
Class to store relative quantum efficiency of channels w.r.t initial one measured in PMT QA QE is exp...
Definition TOPCalChannelRQE.h:24
Belle2::TOPCalPhotonYields
Class to store photon pixel yields for PMT ageing studies, and equalized alpha ratios for finding opt...
Definition TOPCalPhotonYields.h:26
Belle2::TOP::TOPGeometryPar::Instance
static TOPGeometryPar * Instance()
Static method to obtain the pointer to its instance.
Definition TOPGeometryPar.cc:43
Belle2::TOP::TOPPhotonYieldsAlgorithm::m_maxErrorRQE
double m_maxErrorRQE
maximal error on RQE required for good calibration
Definition TOPPhotonYieldsAlgorithm.h:83
Belle2::TOP::TOPPhotonYieldsAlgorithm::calibrate
virtual EResult calibrate() final
algorithm implementation
Definition TOPPhotonYieldsAlgorithm.cc:35
Belle2::TOP::TOPPhotonYieldsAlgorithm::getEqualizingValue
double getEqualizingValue(int bin) const
Returns equalizing value for alpha ratio.
Definition TOPPhotonYieldsAlgorithm.cc:229
Belle2::TOP::TOPPhotonYieldsAlgorithm::equalize
void equalize(TH1F *h) const
Equalize alpha ratio histogram.
Definition TOPPhotonYieldsAlgorithm.cc:248
Belle2::TOP::TOPPhotonYieldsAlgorithm::getMuonCorrection
double getMuonCorrection(const TH1F *h_mu, int row) const
Returns photon yield correction factor according to muon local-z distribution.
Definition TOPPhotonYieldsAlgorithm.cc:285
Belle2::TOP::TOPPhotonYieldsAlgorithm::getNominalYield
double getNominalYield(int slot, int row, int col) const
Returns photon yield of nominal QE for given slot and pixel.
Definition TOPPhotonYieldsAlgorithm.cc:258
Belle2::TOP::TOPPhotonYieldsAlgorithm::m_maxRQE
double m_maxRQE
maximal RQE required for good calibration
Definition TOPPhotonYieldsAlgorithm.h:82
Belle2::CalibrationAlgorithm::getObjectPtr
std::shared_ptr< T > getObjectPtr(const std::string &name, const std::vector< Calibration::ExpRun > &requestedRuns)
Get calibration data object by name and list of runs, the Merge function will be called to generate t...
Definition CalibrationAlgorithm.h:414
Belle2
Abstract base class for different kinds of events.
Definition MillepedeAlgorithm.h:17
std
STL namespace.