CMS 3D CMS Logo
MaterialBudgetMtdHistos.cc
Go to the documentation of this file.
1 #include "Validation/Geometry/interface/MaterialBudgetMtdHistos.h"
2 #include "Validation/Geometry/interface/MaterialBudgetData.h"
3 
4 #include "CLHEP/Units/GlobalSystemOfUnits.h"
5 
6 MaterialBudgetMtdHistos::MaterialBudgetMtdHistos(std::shared_ptr<MaterialBudgetData> data,
7  std::shared_ptr<TestHistoMgr> mgr,
8  const std::string& fileName)
9  : MaterialBudgetFormat(data), hmgr(mgr) {
10  theFileName = fileName;
11  book();
12 }
13 
14 void MaterialBudgetMtdHistos::book() {
15  edm::LogInfo("MaterialBudget") << "MaterialBudgetMtdHistos: Booking user histos";
16 
17  static constexpr double minEta = -5.;
18  static constexpr double maxEta = 5.;
19  static constexpr double minPhi = -3.1416;
20  static constexpr double maxPhi = 3.1416;
21  static constexpr int nbinEta = 250;
22  static constexpr int nbinPhi = 180;
23 
24  static constexpr double minEtaBTLZoom = 0.;
25  static constexpr double maxEtaBTLZoom = 0.087;
26  static constexpr double minPhiBTLZoom = 0.;
27  static constexpr double maxPhiBTLZoom = 0.35;
28  static constexpr int nbinEtaBTLZoom = 64;
29  static constexpr int nbinPhiBTLZoom = 20;
30 
31  static constexpr double minMB = 0.;
32  static constexpr double maxMBetl = 0.025;
33  static constexpr double maxMBbtl = 0.5;
34  static constexpr int nbinMB = 25;
35 
36  // Material budget: radiation length
37  // total X0
38  hmgr->addHistoProf1(new TProfile("10", "MB prof Eta [Total];#eta;x/X_{0} ", nbinEta, minEta, maxEta));
39  hmgr->addHisto1(new TH1F("11", "Eta ", nbinEta, minEta, maxEta));
40  hmgr->addHistoProf1(new TProfile("20", "MB prof Phi [Total];#varphi [rad];x/X_{0} ", nbinPhi, minPhi, maxPhi));
41  hmgr->addHisto1(new TH1F("21", "Phi ", nbinPhi, minPhi, maxPhi));
42  hmgr->addHistoProf2(new TProfile2D(
43  "30", "MB prof Eta Phi [Total];#eta;#varphi;x/X_{0} ", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
44  hmgr->addHisto2(new TH2F("31", "Eta vs Phi ", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
45 
46  // Support
47  hmgr->addHistoProf1(new TProfile("110", "MB prof Eta [Support];#eta;x/X_{0}", nbinEta, minEta, maxEta));
48  hmgr->addHisto1(new TH1F("111", "Eta [Support]", nbinEta, minEta, maxEta));
49  hmgr->addHistoProf1(new TProfile("120", "MB prof Phi [Support];#varphi [rad];x/X_{0}", nbinPhi, minPhi, maxPhi));
50  hmgr->addHisto1(new TH1F("121", "Phi [Support]", nbinPhi, minPhi, maxPhi));
51  hmgr->addHistoProf2(new TProfile2D(
52  "130", "MB prof Eta Phi [Support];#eta;#varphi;x/X_{0}", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
53  hmgr->addHisto2(new TH2F("131", "Eta vs Phi [Support]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
54 
55  // Sensitive
56  hmgr->addHistoProf1(new TProfile("210", "MB prof Eta [Sensitive];#eta;x/X_{0}", nbinEta, minEta, maxEta));
57  hmgr->addHisto1(new TH1F("211", "Eta [Sensitive]", nbinEta, minEta, maxEta));
58  hmgr->addHistoProf1(new TProfile("220", "MB prof Phi [Sensitive];#varphi [rad];x/X_{0}", nbinPhi, minPhi, maxPhi));
59  hmgr->addHisto1(new TH1F("221", "Phi [Sensitive]", nbinPhi, minPhi, maxPhi));
60  hmgr->addHistoProf2(new TProfile2D(
61  "230", "MB prof Eta Phi [Sensitive];#eta;#varphi;x/X_{0}", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
62 
63  hmgr->addHistoProf2(new TProfile2D("10230",
64  "MB prof Eta Phi [Sensitive];#eta;#varphi;x/X_{0}",
65  nbinEtaBTLZoom,
66  minEtaBTLZoom,
67  maxEtaBTLZoom,
68  nbinPhiBTLZoom,
69  minPhiBTLZoom,
70  maxPhiBTLZoom));
71  hmgr->addHisto2(
72  new TH2F("10234", "MB vs Eta [Sensitive];#eta;x/X_{0}", nbinEta, minEta, maxEta, nbinMB, minMB, maxMBetl));
73  hmgr->addHisto2(new TH2F(
74  "20234", "MB along z vs Eta [Sensitive];#eta;x/X_{0}", nbinEta, minEta, maxEta, nbinMB, minMB, maxMBetl));
75  hmgr->addHisto2(
76  new TH2F("10235", "MB vs Eta [Sensitive];#eta;x/X_{0}", nbinEta, minEta, maxEta, nbinMB, minMB, maxMBbtl));
77 
78  hmgr->addHisto2(new TH2F("231", "Eta vs Phi [Sensitive]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
79 
80  // Cables
81  hmgr->addHistoProf1(new TProfile("310", "MB prof Eta [Cables];#eta;x/X_{0}", nbinEta, minEta, maxEta));
82  hmgr->addHisto1(new TH1F("311", "Eta [Cables]", nbinEta, minEta, maxEta));
83  hmgr->addHistoProf1(new TProfile("320", "MB prof Phi [Cables];#varphi [rad];x/X_{0}", nbinPhi, minPhi, maxPhi));
84  hmgr->addHisto1(new TH1F("321", "Phi [Cables]", nbinPhi, minPhi, maxPhi));
85  hmgr->addHistoProf2(new TProfile2D(
86  "330", "MB prof Eta Phi [Cables];#eta;#varphi;x/X_{0}", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
87  hmgr->addHisto2(new TH2F("331", "Eta vs Phi [Cables]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
88 
89  // Cooling
90  hmgr->addHistoProf1(new TProfile("410", "MB prof Eta [Cooling];#eta;x/X_{0}", nbinEta, minEta, maxEta));
91  hmgr->addHisto1(new TH1F("411", "Eta [Cooling]", nbinEta, minEta, maxEta));
92  hmgr->addHistoProf1(new TProfile("420", "MB prof Phi [Cooling];#varphi [rad];x/X_{0}", nbinPhi, minPhi, maxPhi));
93  hmgr->addHisto1(new TH1F("421", "Phi [Cooling]", nbinPhi, minPhi, maxPhi));
94  hmgr->addHistoProf2(new TProfile2D(
95  "430", "MB prof Eta Phi [Cooling];#eta;#varphi;x/X_{0}", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
96  hmgr->addHisto2(new TH2F("431", "Eta vs Phi [Cooling]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
97 
98  // Electronics
99  hmgr->addHistoProf1(new TProfile("510", "MB prof Eta [Electronics];#eta;x/X_{0}", nbinEta, minEta, maxEta));
100  hmgr->addHisto1(new TH1F("511", "Eta [Electronics]", nbinEta, minEta, maxEta));
101  hmgr->addHistoProf1(new TProfile("520", "MB prof Phi [Electronics];#varphi [rad];x/X_{0}", nbinPhi, minPhi, maxPhi));
102  hmgr->addHisto1(new TH1F("521", "Phi [Electronics]", nbinPhi, minPhi, maxPhi));
103  hmgr->addHistoProf2(new TProfile2D(
104  "530", "MB prof Eta Phi [Electronics];#eta;#varphi;x/X_{0}", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
105  hmgr->addHisto2(new TH2F("531", "Eta vs Phi [Electronics]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
106 
107  // Other
108  hmgr->addHistoProf1(new TProfile("610", "MB prof Eta [Other];#eta;x/X_{0}", nbinEta, minEta, maxEta));
109  hmgr->addHisto1(new TH1F("611", "Eta [Other]", nbinEta, minEta, maxEta));
110  hmgr->addHistoProf1(new TProfile("620", "MB prof Phi [Other];#varphi [rad];x/X_{0}", nbinPhi, minPhi, maxPhi));
111  hmgr->addHisto1(new TH1F("621", "Phi [Other]", nbinPhi, minPhi, maxPhi));
112  hmgr->addHistoProf2(new TProfile2D(
113  "630", "MB prof Eta Phi [Other];#eta;#varphi;x/X_{0}", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
114  hmgr->addHisto2(new TH2F("631", "Eta vs Phi [Other]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
115 
116  // Material budget: interaction length
117  // total Lambda0
118  hmgr->addHistoProf1(new TProfile("1010", "MB prof Eta [Total];#eta;#lambda/#lambda_{0} ", nbinEta, minEta, maxEta));
119  hmgr->addHisto1(new TH1F("1011", "Eta ", nbinEta, minEta, maxEta));
120  hmgr->addHistoProf1(
121  new TProfile("1020", "MB prof Phi [Total];#varphi [rad];#lambda/#lambda_{0} ", nbinPhi, minPhi, maxPhi));
122  hmgr->addHisto1(new TH1F("1021", "Phi ", nbinPhi, minPhi, maxPhi));
123  hmgr->addHistoProf2(new TProfile2D("1030",
124  "MB prof Eta Phi [Total];#eta;#varphi;#lambda/#lambda_{0} ",
125  nbinEta,
126  minEta,
127  maxEta,
128  nbinPhi,
129  minPhi,
130  maxPhi));
131  hmgr->addHisto2(new TH2F("1031", "Eta vs Phi ", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
132 
133  // Support
134  hmgr->addHistoProf1(new TProfile("1110", "MB prof Eta [Support];#eta;#lambda/#lambda_{0}", nbinEta, minEta, maxEta));
135  hmgr->addHisto1(new TH1F("1111", "Eta [Support]", nbinEta, minEta, maxEta));
136  hmgr->addHistoProf1(
137  new TProfile("1120", "MB prof Phi [Support];#varphi [rad];#lambda/#lambda_{0}", nbinPhi, minPhi, maxPhi));
138  hmgr->addHisto1(new TH1F("1121", "Phi [Support]", nbinPhi, minPhi, maxPhi));
139  hmgr->addHistoProf2(new TProfile2D("1130",
140  "MB prof Eta Phi [Support];#eta;#varphi;#lambda/#lambda_{0}",
141  nbinEta,
142  minEta,
143  maxEta,
144  nbinPhi,
145  minPhi,
146  maxPhi));
147  hmgr->addHisto2(new TH2F("1131", "Eta vs Phi [Support]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
148 
149  // Sensitive
150  hmgr->addHistoProf1(
151  new TProfile("1210", "MB prof Eta [Sensitive];#eta;#lambda/#lambda_{0}", nbinEta, minEta, maxEta));
152  hmgr->addHisto1(new TH1F("1211", "Eta [Sensitive]", nbinEta, minEta, maxEta));
153  hmgr->addHistoProf1(
154  new TProfile("1220", "MB prof Phi [Sensitive];#varphi [rad];#lambda/#lambda_{0}", nbinPhi, minPhi, maxPhi));
155  hmgr->addHisto1(new TH1F("1221", "Phi [Sensitive]", nbinPhi, minPhi, maxPhi));
156  hmgr->addHistoProf2(new TProfile2D("1230",
157  "MB prof Eta Phi [Sensitive];#eta;#varphi;#lambda/#lambda_{0}",
158  nbinEta,
159  minEta,
160  maxEta,
161  nbinPhi,
162  minPhi,
163  maxPhi));
164  hmgr->addHisto2(new TH2F("1231", "Eta vs Phi [Sensitive]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
165 
166  // Cables
167  hmgr->addHistoProf1(new TProfile("1310", "MB prof Eta [Cables];#eta;#lambda/#lambda_{0}", nbinEta, minEta, maxEta));
168  hmgr->addHisto1(new TH1F("1311", "Eta [Cables]", nbinEta, minEta, maxEta));
169  hmgr->addHistoProf1(
170  new TProfile("1320", "MB prof Phi [Cables];#varphi [rad];#lambda/#lambda_{0}", nbinPhi, minPhi, maxPhi));
171  hmgr->addHisto1(new TH1F("1321", "Phi [Cables]", nbinPhi, minPhi, maxPhi));
172  hmgr->addHistoProf2(new TProfile2D("1330",
173  "MB prof Eta Phi [Cables];#eta;#varphi;#lambda/#lambda_{0}",
174  nbinEta,
175  minEta,
176  maxEta,
177  nbinPhi,
178  minPhi,
179  maxPhi));
180  hmgr->addHisto2(new TH2F("1331", "Eta vs Phi [Cables]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
181 
182  // Cooling
183  hmgr->addHistoProf1(new TProfile("1410", "MB prof Eta [Cooling];#eta;#lambda/#lambda_{0}", nbinEta, minEta, maxEta));
184  hmgr->addHisto1(new TH1F("1411", "Eta [Cooling]", nbinEta, minEta, maxEta));
185  hmgr->addHistoProf1(
186  new TProfile("1420", "MB prof Phi [Cooling];#varphi [rad];#lambda/#lambda_{0}", nbinPhi, minPhi, maxPhi));
187  hmgr->addHisto1(new TH1F("1421", "Phi [Cooling]", nbinPhi, minPhi, maxPhi));
188  hmgr->addHistoProf2(new TProfile2D("1430",
189  "MB prof Eta Phi [Cooling];#eta;#varphi;#lambda/#lambda_{0}",
190  nbinEta,
191  minEta,
192  maxEta,
193  nbinPhi,
194  minPhi,
195  maxPhi));
196  hmgr->addHisto2(new TH2F("1431", "Eta vs Phi [Cooling]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
197 
198  // Electronics
199  hmgr->addHistoProf1(
200  new TProfile("1510", "MB prof Eta [Electronics];#eta;#lambda/#lambda_{0}", nbinEta, minEta, maxEta));
201  hmgr->addHisto1(new TH1F("1511", "Eta [Electronics]", nbinEta, minEta, maxEta));
202  hmgr->addHistoProf1(
203  new TProfile("1520", "MB prof Phi [Electronics];#varphi [rad];#lambda/#lambda_{0}", nbinPhi, minPhi, maxPhi));
204  hmgr->addHisto1(new TH1F("1521", "Phi [Electronics]", nbinPhi, minPhi, maxPhi));
205  hmgr->addHistoProf2(new TProfile2D("1530",
206  "MB prof Eta Phi [Electronics];#eta;#varphi;#lambda/#lambda_{0}",
207  nbinEta,
208  minEta,
209  maxEta,
210  nbinPhi,
211  minPhi,
212  maxPhi));
213  hmgr->addHisto2(new TH2F("1531", "Eta vs Phi [Electronics]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
214 
215  // Other
216  hmgr->addHistoProf1(new TProfile("1610", "MB prof Eta [Other];#eta;#lambda/#lambda_{0}", nbinEta, minEta, maxEta));
217  hmgr->addHisto1(new TH1F("1611", "Eta [Other]", nbinEta, minEta, maxEta));
218  hmgr->addHistoProf1(
219  new TProfile("1620", "MB prof Phi [Other];#varphi [rad];#lambda/#lambda_{0}", nbinPhi, minPhi, maxPhi));
220  hmgr->addHisto1(new TH1F("1621", "Phi [Other]", nbinPhi, minPhi, maxPhi));
221  hmgr->addHistoProf2(new TProfile2D("1630",
222  "MB prof Eta Phi [Other];#eta;#varphi;#lambda/#lambda_{0}",
223  nbinEta,
224  minEta,
225  maxEta,
226  nbinPhi,
227  minPhi,
228  maxPhi));
229  hmgr->addHisto2(new TH2F("1631", "Eta vs Phi [Other]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
230 
231  edm::LogInfo("MaterialBudget") << "MaterialBudgetMtdHistos: booking user histos done";
232 }
233 
234 void MaterialBudgetMtdHistos::fillStartTrack() {}
235 
236 void MaterialBudgetMtdHistos::fillPerStep() {}
237 
238 void MaterialBudgetMtdHistos::fillEndTrack() {
239  // Total X0
240  hmgr->getHisto1(11)->Fill(theData->getEta());
241  hmgr->getHisto1(21)->Fill(theData->getPhi());
242  hmgr->getHisto2(31)->Fill(theData->getEta(), theData->getPhi());
243 
244  hmgr->getHistoProf1(10)->Fill(theData->getEta(), theData->getTotalMB());
245  hmgr->getHistoProf1(20)->Fill(theData->getPhi(), theData->getTotalMB());
246  hmgr->getHistoProf2(30)->Fill(theData->getEta(), theData->getPhi(), theData->getTotalMB());
247 
248  // Support
249  hmgr->getHisto1(111)->Fill(theData->getEta());
250  hmgr->getHisto1(121)->Fill(theData->getPhi());
251  hmgr->getHisto2(131)->Fill(theData->getEta(), theData->getPhi());
252 
253  hmgr->getHistoProf1(110)->Fill(theData->getEta(), theData->getSupportMB());
254  hmgr->getHistoProf1(120)->Fill(theData->getPhi(), theData->getSupportMB());
255  hmgr->getHistoProf2(130)->Fill(theData->getEta(), theData->getPhi(), theData->getSupportMB());
256 
257  // Sensitive
258  hmgr->getHisto1(211)->Fill(theData->getEta());
259  hmgr->getHisto1(221)->Fill(theData->getPhi());
260  hmgr->getHisto2(231)->Fill(theData->getEta(), theData->getPhi());
261 
262  hmgr->getHistoProf1(210)->Fill(theData->getEta(), theData->getSensitiveMB());
263  hmgr->getHistoProf1(220)->Fill(theData->getPhi(), theData->getSensitiveMB());
264  hmgr->getHistoProf2(230)->Fill(theData->getEta(), theData->getPhi(), theData->getSensitiveMB());
265 
266  static constexpr double bfTransitionEta = 1.55;
267 
268  if (std::abs(theData->getEta()) > bfTransitionEta) {
269  hmgr->getHisto2(10234)->Fill(theData->getEta(), theData->getSensitiveMB());
270  double norma = std::tanh(std::abs(theData->getEta()));
271  hmgr->getHisto2(20234)->Fill(theData->getEta(), theData->getSensitiveMB() * norma);
272  } else {
273  hmgr->getHistoProf2(10230)->Fill(theData->getEta(), theData->getPhi(), theData->getSensitiveMB());
274  hmgr->getHisto2(10235)->Fill(theData->getEta(), theData->getSensitiveMB());
275  }
276 
277  // Cables
278  hmgr->getHisto1(311)->Fill(theData->getEta());
279  hmgr->getHisto1(321)->Fill(theData->getPhi());
280  hmgr->getHisto2(331)->Fill(theData->getEta(), theData->getPhi());
281 
282  hmgr->getHistoProf1(310)->Fill(theData->getEta(), theData->getCablesMB());
283  hmgr->getHistoProf1(320)->Fill(theData->getPhi(), theData->getCablesMB());
284  hmgr->getHistoProf2(330)->Fill(theData->getEta(), theData->getPhi(), theData->getCablesMB());
285 
286  // Cooling
287  hmgr->getHisto1(411)->Fill(theData->getEta());
288  hmgr->getHisto1(421)->Fill(theData->getPhi());
289  hmgr->getHisto2(431)->Fill(theData->getEta(), theData->getPhi());
290 
291  hmgr->getHistoProf1(410)->Fill(theData->getEta(), theData->getCoolingMB());
292  hmgr->getHistoProf1(420)->Fill(theData->getPhi(), theData->getCoolingMB());
293  hmgr->getHistoProf2(430)->Fill(theData->getEta(), theData->getPhi(), theData->getCoolingMB());
294 
295  // Electronics
296  hmgr->getHisto1(511)->Fill(theData->getEta());
297  hmgr->getHisto1(521)->Fill(theData->getPhi());
298  hmgr->getHisto2(531)->Fill(theData->getEta(), theData->getPhi());
299 
300  hmgr->getHistoProf1(510)->Fill(theData->getEta(), theData->getElectronicsMB());
301  hmgr->getHistoProf1(520)->Fill(theData->getPhi(), theData->getElectronicsMB());
302  hmgr->getHistoProf2(530)->Fill(theData->getEta(), theData->getPhi(), theData->getElectronicsMB());
303 
304  // Other
305  hmgr->getHisto1(611)->Fill(theData->getEta());
306  hmgr->getHisto1(621)->Fill(theData->getPhi());
307  hmgr->getHisto2(631)->Fill(theData->getEta(), theData->getPhi());
308 
309  hmgr->getHistoProf1(610)->Fill(theData->getEta(), theData->getOtherMB());
310  hmgr->getHistoProf1(620)->Fill(theData->getPhi(), theData->getOtherMB());
311  hmgr->getHistoProf2(630)->Fill(theData->getEta(), theData->getPhi(), theData->getOtherMB());
312 
313  // Total Lambda0
314  hmgr->getHisto1(1011)->Fill(theData->getEta());
315  hmgr->getHisto1(1021)->Fill(theData->getPhi());
316  hmgr->getHisto2(1031)->Fill(theData->getEta(), theData->getPhi());
317 
318  hmgr->getHistoProf1(1010)->Fill(theData->getEta(), theData->getTotalIL());
319  hmgr->getHistoProf1(1020)->Fill(theData->getPhi(), theData->getTotalIL());
320  hmgr->getHistoProf2(1030)->Fill(theData->getEta(), theData->getPhi(), theData->getTotalIL());
321 
322  // Support
323  hmgr->getHisto1(1111)->Fill(theData->getEta());
324  hmgr->getHisto1(1121)->Fill(theData->getPhi());
325  hmgr->getHisto2(1131)->Fill(theData->getEta(), theData->getPhi());
326 
327  hmgr->getHistoProf1(1110)->Fill(theData->getEta(), theData->getSupportIL());
328  hmgr->getHistoProf1(1120)->Fill(theData->getPhi(), theData->getSupportIL());
329  hmgr->getHistoProf2(1130)->Fill(theData->getEta(), theData->getPhi(), theData->getSupportIL());
330 
331  // Sensitive
332  hmgr->getHisto1(1211)->Fill(theData->getEta());
333  hmgr->getHisto1(1221)->Fill(theData->getPhi());
334  hmgr->getHisto2(1231)->Fill(theData->getEta(), theData->getPhi());
335 
336  hmgr->getHistoProf1(1210)->Fill(theData->getEta(), theData->getSensitiveIL());
337  hmgr->getHistoProf1(1220)->Fill(theData->getPhi(), theData->getSensitiveIL());
338  hmgr->getHistoProf2(1230)->Fill(theData->getEta(), theData->getPhi(), theData->getSensitiveIL());
339 
340  // Cables
341  hmgr->getHisto1(1311)->Fill(theData->getEta());
342  hmgr->getHisto1(1321)->Fill(theData->getPhi());
343  hmgr->getHisto2(1331)->Fill(theData->getEta(), theData->getPhi());
344 
345  hmgr->getHistoProf1(1310)->Fill(theData->getEta(), theData->getCablesIL());
346  hmgr->getHistoProf1(1320)->Fill(theData->getPhi(), theData->getCablesIL());
347  hmgr->getHistoProf2(1330)->Fill(theData->getEta(), theData->getPhi(), theData->getCablesIL());
348 
349  // Cooling
350  hmgr->getHisto1(1411)->Fill(theData->getEta());
351  hmgr->getHisto1(1421)->Fill(theData->getPhi());
352  hmgr->getHisto2(1431)->Fill(theData->getEta(), theData->getPhi());
353 
354  hmgr->getHistoProf1(1410)->Fill(theData->getEta(), theData->getCoolingIL());
355  hmgr->getHistoProf1(1420)->Fill(theData->getPhi(), theData->getCoolingIL());
356  hmgr->getHistoProf2(1430)->Fill(theData->getEta(), theData->getPhi(), theData->getCoolingIL());
357 
358  // Electronics
359  hmgr->getHisto1(1511)->Fill(theData->getEta());
360  hmgr->getHisto1(1521)->Fill(theData->getPhi());
361  hmgr->getHisto2(1531)->Fill(theData->getEta(), theData->getPhi());
362 
363  hmgr->getHistoProf1(1510)->Fill(theData->getEta(), theData->getElectronicsIL());
364  hmgr->getHistoProf1(1520)->Fill(theData->getPhi(), theData->getElectronicsIL());
365  hmgr->getHistoProf2(1530)->Fill(theData->getEta(), theData->getPhi(), theData->getElectronicsIL());
366 
367  // Other
368  hmgr->getHisto1(1611)->Fill(theData->getEta());
369  hmgr->getHisto1(1621)->Fill(theData->getPhi());
370  hmgr->getHisto2(1631)->Fill(theData->getEta(), theData->getPhi());
371 
372  hmgr->getHistoProf1(1610)->Fill(theData->getEta(), theData->getOtherIL());
373  hmgr->getHistoProf1(1620)->Fill(theData->getPhi(), theData->getOtherIL());
374  hmgr->getHistoProf2(1630)->Fill(theData->getEta(), theData->getPhi(), theData->getOtherIL());
375 }
376 
377 void MaterialBudgetMtdHistos::endOfRun() {
378  edm::LogInfo("MaterialBudget") << "MaterialBudgetMtdHistos: Writing histos ROOT file to:" << theFileName;
379  hmgr->save(theFileName);
380 }
HLT_2024v12_cff.minPhi
minPhi
Definition: HLT_2024v12_cff.py:40350
HLT_2024v12_cff.maxPhi
maxPhi
Definition: HLT_2024v12_cff.py:40351
nnet::tanh
void tanh(data_T data[CONFIG_T::n_in], res_T res[CONFIG_T::n_in])
Definition: nnet_activation.h:227
ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr
if constexpr(n > 3)
Definition: BrokenLine.h:164
MaterialBudgetMtdHistos::hmgr
std::shared_ptr< TestHistoMgr > hmgr
Definition: MaterialBudgetMtdHistos.h:27
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
razorScouting_cff.maxEta
maxEta
Definition: razorScouting_cff.py:6
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
EgHLTOffEleSelection_cfi.minEta
minEta
Definition: EgHLTOffEleSelection_cfi.py:11
MillePedeFileConverter_cfg.fileName
fileName
Definition: MillePedeFileConverter_cfg.py:32
data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:80
MaterialBudgetMtdHistos::MaterialBudgetMtdHistos
MaterialBudgetMtdHistos(std::shared_ptr< MaterialBudgetData > data, std::shared_ptr< TestHistoMgr > mgr, const std::string &fileName)
Definition: MaterialBudgetMtdHistos.cc:6
MaterialBudgetFormat::theData
std::shared_ptr< MaterialBudgetData > theData
Definition: MaterialBudgetFormat.h:20
Generated for CMSSW Reference Manual by doxygen 1.8.14