CMS 3D CMS Logo
RecHitTask.cc
Go to the documentation of this file.
1 #include "DQM/HcalTasks/interface/RecHitTask.h"
2 #include <cmath>
3 
4 using namespace hcaldqm;
5 using namespace hcaldqm::constants;
6 using namespace hcaldqm::filter;
7 
8 RecHitTask::RecHitTask(edm::ParameterSet const& ps) : DQTask(ps) {
9  _tagHBHE = ps.getUntrackedParameter<edm::InputTag>("tagHBHE", edm::InputTag("hbhereco"));
10  _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO", edm::InputTag("horeco"));
11  _tagHF = ps.getUntrackedParameter<edm::InputTag>("tagHF", edm::InputTag("hfreco"));
12  _tagPreHF = ps.getUntrackedParameter<edm::InputTag>("tagPreHF", edm::InputTag(""));
13  _hfPreRecHitsAvailable = ps.getUntrackedParameter<bool>("hfPreRecHitsAvailable", false);
14 
15  _tokHBHE = consumes<HBHERecHitCollection>(_tagHBHE);
16  _tokHO = consumes<HORecHitCollection>(_tagHO);
17  _tokHF = consumes<HFRecHitCollection>(_tagHF);
18  _tokPreHF = consumes<HFPreRecHitCollection>(_tagPreHF);
19 
20  _cutE_HBHE = ps.getUntrackedParameter<double>("cutE_HBHE", 5);
21  _cutE_HO = ps.getUntrackedParameter<double>("cutE_HO", 5);
22  _cutE_HF = ps.getUntrackedParameter<double>("cutE_HF", 5);
23  _thresh_unihf = ps.getUntrackedParameter<double>("thresh_unihf", 0.2);
24 
25  // order must be the same as in RecoFlag enum
26  _vflags.resize(nRecoFlag);
27  _vflags[fUni] = flag::Flag("UniSlotHF");
28  _vflags[fTCDS] = flag::Flag("TCDS");
29  _vflags[fUnknownIds] = flag::Flag("UnknownIds");
30 }
31 
32 /* virtual */ void RecHitTask::bookHistograms(DQMStore::IBooker& ib, edm::Run const& r, edm::EventSetup const& es) {
33  DQTask::bookHistograms(ib, r, es);
34 
35  // GET WHAT YOU NEED
36  edm::ESHandle<HcalDbService> dbs;
37  es.get<HcalDbRecord>().get(dbs);
38  _emap = dbs->getHcalMapping();
39 
40  std::vector<uint32_t> vVME;
41  std::vector<uint32_t> vuTCA;
42  vVME.push_back(
43  HcalElectronicsId(constants::FIBERCH_MIN, constants::FIBER_VME_MIN, SPIGOT_MIN, CRATE_VME_MIN).rawId());
44  vuTCA.push_back(HcalElectronicsId(CRATE_uTCA_MIN, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
45  _filter_VME.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics, vVME);
46  _filter_uTCA.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics, vuTCA);
47  std::vector<uint32_t> vhashHF;
48  vhashHF.push_back(
49  hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdet](HcalDetId(HcalForward, 29, 1, 1)));
50 
51  _filter_HF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fSubdet, vhashHF);
52 
53  // INITIALIZE FIRST
54  // Energy
55  _cEnergy_Subdet.initialize(_name,
56  "Energy",
57  hcaldqm::hashfunctions::fSubdet,
58  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy),
59  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
60  0);
61  _cEnergy_depth.initialize(_name,
62  "Energy",
63  hcaldqm::hashfunctions::fdepth,
64  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
65  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
66  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy, true),
67  0);
68 
69  // Timing
70  _cTimingCut_SubdetPM.initialize(_name,
71  "TimingCut",
72  hcaldqm::hashfunctions::fSubdetPM,
73  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
74  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
75  0);
76  _cTimingvsEnergy_SubdetPM.initialize(_name,
77  "TimingvsEnergy",
78  hcaldqm::hashfunctions::fSubdetPM,
79  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy, true),
80  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
81  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
82  0);
83 
84  _cTimingCut_HBHEPartition.initialize(_name,
85  "TimingCut",
86  hcaldqm::hashfunctions::fHBHEPartition,
87  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
88  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
89  0);
90  _cTimingCut_depth.initialize(_name,
91  "TimingCut",
92  hcaldqm::hashfunctions::fdepth,
93  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
94  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
95  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
96  0);
97  _cTimingCutvsLS_SubdetPM.initialize(_name,
98  "TimingCutvsLS",
99  hcaldqm::hashfunctions::fSubdetPM,
100  new hcaldqm::quantity::LumiSection(_maxLS),
101  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
102  0);
103 
104  // Occupancy
105  _cOccupancy_depth.initialize(_name,
106  "Occupancy",
107  hcaldqm::hashfunctions::fdepth,
108  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
109  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
110  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
111  0);
112 
113  _cOccupancyvsLS_Subdet.initialize(_name,
114  "OccupancyvsLS",
115  hcaldqm::hashfunctions::fSubdet,
116  new hcaldqm::quantity::LumiSection(_maxLS),
117  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),
118  0);
119 
120  _cOccupancyCut_depth.initialize(_name,
121  "OccupancyCut",
122  hcaldqm::hashfunctions::fdepth,
123  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
124  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
125  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
126  0);
127 
128  if (_hfPreRecHitsAvailable) {
129  _cDAAsymmetryVsCharge_SubdetPM.initialize(
130  _name,
131  "ChargeVsAsymmetry",
132  hcaldqm::hashfunctions::fSubdetPM,
133  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDualAnodeAsymmetry),
134  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000),
135  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
136  0);
137  _cDAAsymmetryMean_cut_depth.initialize(_name,
138  "AsymmetryMean",
139  hcaldqm::hashfunctions::fdepth,
140  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
141  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
142  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDualAnodeAsymmetry),
143  0);
144  _cDAAsymmetry_cut_SubdetPM.initialize(_name,
145  "Asymmetry",
146  hcaldqm::hashfunctions::fSubdetPM,
147  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDualAnodeAsymmetry),
148  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
149  0);
150  }
151 
152  // INITIALIZE HISTOGRAMS to be used only in Online
153  if (_ptype == fOnline) {
154  _cEnergyvsieta_Subdet.initialize(_name,
155  "Energyvsieta",
156  hcaldqm::hashfunctions::fSubdet,
157  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
158  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),
159  0);
160  _cEnergyvsiphi_SubdetPM.initialize(_name,
161  "Energyvsiphi",
162  hcaldqm::hashfunctions::fSubdetPM,
163  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
164  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),
165  0);
166  _cEnergyvsLS_SubdetPM.initialize(_name,
167  "EnergyvsLS",
168  hcaldqm::hashfunctions::fSubdetPM,
169  new hcaldqm::quantity::LumiSection(_maxLS),
170  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),
171  0);
172  _cEnergyvsBX_SubdetPM.initialize(_name,
173  "EnergyvsBX",
174  hcaldqm::hashfunctions::fSubdetPM,
175  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
176  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),
177  0);
178  _cTimingCutvsieta_Subdet.initialize(_name,
179  "TimingCutvsieta",
180  hcaldqm::hashfunctions::fSubdet,
181  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
182  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
183  0);
184  _cTimingCutvsiphi_SubdetPM.initialize(_name,
185  "TimingCutvsiphi",
186  hcaldqm::hashfunctions::fSubdetPM,
187  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
188  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
189  0);
190  _cTimingCutvsBX_SubdetPM.initialize(_name,
191  "TimingCutvsBX",
192  hcaldqm::hashfunctions::fSubdetPM,
193  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
194  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
195  0);
196  _cOccupancyvsiphi_SubdetPM.initialize(_name,
197  "Occupancyvsiphi",
198  hcaldqm::hashfunctions::fSubdetPM,
199  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
200  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
201  0);
202  _cOccupancyvsieta_Subdet.initialize(_name,
203  "Occupancyvsieta",
204  hcaldqm::hashfunctions::fSubdet,
205  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
206  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
207  0);
208  _cOccupancyCutvsiphi_SubdetPM.initialize(_name,
209  "OccupancyCutvsiphi",
210  hcaldqm::hashfunctions::fSubdetPM,
211  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
212  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
213  0);
214  _cOccupancyCutvsieta_Subdet.initialize(_name,
215  "OccupancyCutvsieta",
216  hcaldqm::hashfunctions::fSubdet,
217  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
218  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
219  0);
220  _cOccupancyCutvsBX_Subdet.initialize(_name,
221  "OccupancyCutvsBX",
222  hcaldqm::hashfunctions::fSubdet,
223  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
224  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
225  0);
226  _cOccupancyCutvsiphivsLS_SubdetPM.initialize(_name,
227  "OccupancyCutvsiphivsLS",
228  hcaldqm::hashfunctions::fSubdetPM,
229  new hcaldqm::quantity::LumiSection(_maxLS),
230  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
231  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
232  0);
233  _cOccupancyCutvsLS_Subdet.initialize(_name,
234  "OccupancyCutvsLS",
235  hcaldqm::hashfunctions::fSubdet,
236  new hcaldqm::quantity::LumiSection(_maxLS),
237  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),
238  0);
239  }
240 
241  // FED-based plots
242  if (_ptype != fOffline) { // hidefed2crate
243  std::vector<int> vFEDs = hcaldqm::utilities::getFEDList(_emap);
244  std::vector<int> vFEDsVME = hcaldqm::utilities::getFEDVMEList(_emap);
245  std::vector<int> vFEDsuTCA = hcaldqm::utilities::getFEDuTCAList(_emap);
246 
247  // push the rawIds of each fed into the vector
248  for (std::vector<int>::const_iterator it = vFEDsVME.begin(); it != vFEDsVME.end(); ++it)
249  _vhashFEDs.push_back(HcalElectronicsId(FIBERCH_MIN, FIBER_VME_MIN, SPIGOT_MIN, (*it) - FED_VME_MIN).rawId());
250  for (std::vector<int>::const_iterator it = vFEDsuTCA.begin(); it != vFEDsuTCA.end(); ++it) {
251  std::pair<uint16_t, uint16_t> cspair = utilities::fed2crate(*it);
252  _vhashFEDs.push_back(HcalElectronicsId(cspair.first, cspair.second, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
253  }
254 
255  _cTimingCut_FEDVME.initialize(_name,
256  "TimingCut",
257  hcaldqm::hashfunctions::fFED,
258  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
259  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
260  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
261  0);
262  _cTimingCut_FEDuTCA.initialize(_name,
263  "TimingCut",
264  hcaldqm::hashfunctions::fFED,
265  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
266  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
267  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
268  0);
269  _cTimingCutvsLS_FED.initialize(_name,
270  "TimingCutvsLS",
271  hcaldqm::hashfunctions::fFED,
272  new hcaldqm::quantity::LumiSection(_maxLS),
273  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
274  0);
275  _cTimingCut_ElectronicsVME.initialize(_name,
276  "TimingCut",
277  hcaldqm::hashfunctions::fElectronics,
278  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
279  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
280  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
281  0);
282  _cTimingCut_ElectronicsuTCA.initialize(_name,
283  "TimingCut",
284  hcaldqm::hashfunctions::fElectronics,
285  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
286  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
287  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
288  0);
289 
290  _cOccupancy_FEDVME.initialize(_name,
291  "Occupancy",
292  hcaldqm::hashfunctions::fFED,
293  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
294  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
295  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
296  0);
297  _cOccupancy_FEDuTCA.initialize(_name,
298  "Occupancy",
299  hcaldqm::hashfunctions::fFED,
300  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
301  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
302  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
303  0);
304  _cOccupancy_ElectronicsVME.initialize(_name,
305  "Occupancy",
306  hcaldqm::hashfunctions::fElectronics,
307  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
308  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
309  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
310  0);
311  _cOccupancy_ElectronicsuTCA.initialize(_name,
312  "Occupancy",
313  hcaldqm::hashfunctions::fElectronics,
314  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
315  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
316  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
317  0);
318 
319  _cOccupancyCut_FEDVME.initialize(_name,
320  "OccupancyCut",
321  hcaldqm::hashfunctions::fFED,
322  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
323  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
324  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
325  0);
326  _cOccupancyCut_FEDuTCA.initialize(_name,
327  "OccupancyCut",
328  hcaldqm::hashfunctions::fFED,
329  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
330  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
331  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
332  0);
333  _cOccupancyCut_ElectronicsVME.initialize(_name,
334  "OccupancyCut",
335  hcaldqm::hashfunctions::fElectronics,
336  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
337  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
338  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
339  0);
340  _cOccupancyCut_ElectronicsuTCA.initialize(_name,
341  "OccupancyCut",
342  hcaldqm::hashfunctions::fElectronics,
343  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
344  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
345  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
346  0);
347  if (_ptype == fOnline) {
348  _cSummaryvsLS_FED.initialize(_name,
349  "SummaryvsLS",
350  hcaldqm::hashfunctions::fFED,
351  new hcaldqm::quantity::LumiSection(_maxLS),
352  new hcaldqm::quantity::FlagQuantity(_vflags),
353  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),
354  0);
355  _cSummaryvsLS.initialize(_name,
356  "SummaryvsLS",
357  new hcaldqm::quantity::LumiSection(_maxLS),
358  new hcaldqm::quantity::FEDQuantity(vFEDs),
359  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),
360  0);
361 
362  _xUniHF.initialize(hcaldqm::hashfunctions::fFEDSlot);
363  _xUni.initialize(hcaldqm::hashfunctions::fFED);
364  }
365  }
366 
367  // BOOK HISTOGRAMS
368  char cutstr[200];
369  sprintf(cutstr, "_EHBHE%dHO%dHF%d", int(_cutE_HBHE), int(_cutE_HO), int(_cutE_HF));
370  char cutstr2[200];
371  sprintf(cutstr2, "_EHF%d", int(_cutE_HF));
372 
373  // Energy
374  _cEnergy_Subdet.book(ib, _emap, _subsystem);
375  _cEnergy_depth.book(ib, _emap, _subsystem);
376 
377  // Timing
378  _cTimingCut_SubdetPM.book(ib, _emap, _subsystem);
379  _cTimingvsEnergy_SubdetPM.book(ib, _emap, _subsystem);
380  if (_ptype != fOffline) { // hidefed2crate
381  _cTimingCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
382  _cTimingCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
383  _cTimingCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
384  _cTimingCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
385  }
386  _cTimingCut_HBHEPartition.book(ib, _emap, _subsystem);
387  _cTimingCut_depth.book(ib, _emap, _subsystem);
388  _cTimingCutvsLS_SubdetPM.book(ib, _emap, _subsystem);
389  if (_ptype != fOffline) { // hidefed2crate
390  _cTimingCutvsLS_FED.book(ib, _emap, _subsystem);
391  }
392 
393  // Occupancy
394  _cOccupancy_depth.book(ib, _emap, _subsystem);
395  if (_ptype != fOffline) { // hidefed2crate
396  _cOccupancy_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
397  _cOccupancy_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
398  _cOccupancy_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
399  _cOccupancy_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
400  }
401  _cOccupancyvsLS_Subdet.book(ib, _emap, _subsystem);
402  _cOccupancyCut_depth.book(ib, _emap, _subsystem);
403  if (_ptype != fOffline) { // hidefed2crate
404  _cOccupancyCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
405  _cOccupancyCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
406  _cOccupancyCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
407  _cOccupancyCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
408  }
409 
410  if (_hfPreRecHitsAvailable) {
411  _cDAAsymmetryVsCharge_SubdetPM.book(ib, _emap, _filter_HF, _subsystem);
412  _cDAAsymmetryMean_cut_depth.book(ib, _emap, _filter_HF, _subsystem);
413  _cDAAsymmetry_cut_SubdetPM.book(ib, _emap, _filter_HF, _subsystem);
414  }
415 
416  // BOOK HISTOGRAMS to be used only in Online
417  if (_ptype == fOnline) {
418  _cEnergyvsLS_SubdetPM.book(ib, _emap, _subsystem);
419  _cEnergyvsieta_Subdet.book(ib, _emap, _subsystem);
420  _cEnergyvsiphi_SubdetPM.book(ib, _emap, _subsystem);
421  _cEnergyvsBX_SubdetPM.book(ib, _emap, _subsystem);
422  _cTimingCutvsieta_Subdet.book(ib, _emap, _subsystem);
423  _cTimingCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
424  _cTimingCutvsBX_SubdetPM.book(ib, _emap, _subsystem);
425  _cOccupancyvsiphi_SubdetPM.book(ib, _emap, _subsystem);
426  _cOccupancyvsieta_Subdet.book(ib, _emap, _subsystem);
427  _cOccupancyCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
428  _cOccupancyCutvsieta_Subdet.book(ib, _emap, _subsystem);
429  _cOccupancyCutvsBX_Subdet.book(ib, _emap, _subsystem);
430  _cOccupancyCutvsiphivsLS_SubdetPM.book(ib, _emap, _subsystem);
431  _cOccupancyCutvsLS_Subdet.book(ib, _emap, _subsystem);
432  _cSummaryvsLS_FED.book(ib, _emap, _subsystem);
433  _cSummaryvsLS.book(ib, _subsystem);
434 
435  std::vector<uint32_t> vhashFEDHF;
436  vhashFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
437  vhashFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
438  vhashFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
439  vhashFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN + 6, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
440  vhashFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN + 6, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
441  vhashFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN + 6, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
442  HashFilter filter_FEDHF;
443  filter_FEDHF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fFED, vhashFEDHF);
444 
445  _gids = _emap->allPrecisionId();
446  if (_ptype != fOffline) { // hidefed2crate
447  _xUniHF.book(_emap, filter_FEDHF);
448  _xUni.book(_emap);
449  }
450  }
451 
452  // initialize hash map
453  _ehashmap.initialize(_emap, hcaldqm::electronicsmap::fD2EHashMap);
454 
455  // book some mes...
456  ib.setCurrentFolder(_subsystem + "/" + _name);
457  meUnknownIds1LS = ib.book1D("UnknownIds", "UnknownIds", 1, 0, 1);
458  _unknownIdsPresent = false;
459  meUnknownIds1LS->setLumiFlag();
460 }
461 
462 /* virtual */ void RecHitTask::_resetMonitors(hcaldqm::UpdateFreq uf) {
463  switch (uf) {
464  case hcaldqm::f1LS:
465  _unknownIdsPresent = false;
466  break;
467  default:
468  break;
469  }
470 
471  DQTask::_resetMonitors(uf);
472 }
473 
474 /* virtual */ void RecHitTask::_process(edm::Event const& e, edm::EventSetup const&) {
475  edm::Handle<HBHERecHitCollection> chbhe;
476  edm::Handle<HORecHitCollection> cho;
477  edm::Handle<HFRecHitCollection> chf;
478 
479  if (!(e.getByToken(_tokHBHE, chbhe)))
480  _logger.dqmthrow("Collection HBHERecHitCollection not available " + _tagHBHE.label() + " " + _tagHBHE.instance());
481  if (!(e.getByToken(_tokHO, cho)))
482  _logger.dqmthrow("Collection HORecHitCollection not available " + _tagHO.label() + " " + _tagHO.instance());
483  if (!(e.getByToken(_tokHF, chf)))
484  _logger.dqmthrow("Collection HFRecHitCollection not available " + _tagHF.label() + " " + _tagHF.instance());
485 
486  edm::Handle<HFPreRecHitCollection> cprehf;
487  if (_hfPreRecHitsAvailable) {
488  if (!(e.getByToken(_tokPreHF, cprehf)))
489  _logger.dqmthrow("Collection HFPreRecHitCollection not available " + _tagPreHF.label() + " " +
490  _tagPreHF.instance());
491  }
492 
493  // extract some info per event
494  int bx = e.bunchCrossing();
495 
496  // To fill histograms outside of the loop, you need to determine if there were
497  // any valid det ids first
498  uint32_t rawidValid = 0;
499  uint32_t rawidHBValid = 0;
500  uint32_t rawidHEValid = 0;
501 
502  double ehbm = 0;
503  double ehbp = 0;
504  double ehem = 0;
505  double ehep = 0;
506  int nChsHB = 0;
507  int nChsHE = 0;
508  int nChsHBCut = 0;
509  int nChsHECut = 0;
510  for (HBHERecHitCollection::const_iterator it = chbhe->begin(); it != chbhe->end(); ++it) {
511  double energy = it->energy();
512  double timing = it->time();
513 
514  // Explicit check on the DetIds present in the Collection
515  HcalDetId did = it->id();
516  uint32_t rawid = _ehashmap.lookup(did);
517  /*
518  * Needs to be removed as DetIds that belong to the HEP17 after combination
519  * are not present in the emap
520  * Removed until further notice!
521  *
522  */
523  //if (rawid==0)
524  //{meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
525 
526  HcalElectronicsId const& eid(rawid);
527  rawidValid = did.rawId();
528  if (did.subdet() == HcalBarrel)
529  rawidHBValid = did.rawId();
530  else if (did.subdet() == HcalEndcap)
531  rawidHEValid = did.rawId();
532 
533  _cEnergy_Subdet.fill(did, energy);
534  _cTimingvsEnergy_SubdetPM.fill(did, energy, timing);
535  _cOccupancy_depth.fill(did);
536  did.subdet() == HcalBarrel ? did.ieta() > 0 ? ehbp += energy : ehbm += energy
537  : did.ieta() > 0 ? ehep += energy : ehem += energy;
538 
539  // ONLINE ONLY!
540  if (_ptype == fOnline) {
541  _cOccupancyvsiphi_SubdetPM.fill(did);
542  _cOccupancyvsieta_Subdet.fill(did);
543  }
544  // ^^^ONLINE ONLY!
545  //
546  if (_ptype != fOffline) { // hidefed2crate
547  // Also, for these electronics plots, require that the channel was found in the emap.
548  if (rawid != 0) {
549  if (eid.isVMEid()) {
550  _cOccupancy_FEDVME.fill(eid);
551  _cOccupancy_ElectronicsVME.fill(eid);
552  } else if (eid.isUTCAid()) {
553  _cOccupancy_FEDuTCA.fill(eid);
554  _cOccupancy_ElectronicsuTCA.fill(eid);
555  }
556  }
557  }
558 
559  if (energy > _cutE_HBHE) {
560  // ONLINE ONLY!
561  if (_ptype == fOnline) {
562  _cEnergyvsLS_SubdetPM.fill(did, _currentLS, energy);
563  _cEnergyvsBX_SubdetPM.fill(did, bx, energy);
564  _cEnergyvsieta_Subdet.fill(did, energy);
565  _cEnergyvsiphi_SubdetPM.fill(did, energy);
566  _cTimingCutvsieta_Subdet.fill(did, timing);
567  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
568  _cTimingCutvsBX_SubdetPM.fill(did, bx, timing);
569  _cOccupancyCutvsiphi_SubdetPM.fill(did);
570  _cOccupancyCutvsieta_Subdet.fill(did);
571  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
572  }
573  // ^^^ONLINE ONLY!
574  _cEnergy_depth.fill(did, energy);
575  _cTimingCut_SubdetPM.fill(did, timing);
576  _cTimingCut_HBHEPartition.fill(did, timing);
577  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
578 
579  // ONLINE
580  if (_ptype == fOnline) {
581  if (rawid != 0) {
582  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
583  }
584  _cTimingCut_depth.fill(did, timing);
585  // ^^^ONLINE
586  } else {
587  if (_ptype != fOffline) { // hidefed2crate
588  if (rawid != 0) {
589  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
590  }
591  }
592  _cTimingCut_depth.fill(did, timing);
593  }
594  _cOccupancyCut_depth.fill(did);
595  if (_ptype != fOffline) { // hidefed2crate
596  if (rawid != 0) {
597  if (eid.isVMEid()) {
598  // ONLINE
599  if (_ptype == fOnline) {
600  _cTimingCut_FEDVME.fill(eid, timing);
601  _cTimingCut_ElectronicsVME.fill(eid, timing);
602  } // ^^^ ONLINE
603  else {
604  _cTimingCut_FEDVME.fill(eid, timing);
605  _cTimingCut_ElectronicsVME.fill(eid, timing);
606  }
607  // ^^^ONLINE
608 
609  _cOccupancyCut_FEDVME.fill(eid);
610  _cOccupancyCut_ElectronicsVME.fill(eid);
611  } else if (eid.isUTCAid()) {
612  if (_ptype == fOnline) {
613  // time constraints are explicit!
614  _cTimingCut_FEDuTCA.fill(eid, timing);
615  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
616  } else {
617  _cTimingCut_FEDuTCA.fill(eid, timing);
618  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
619  }
620  _cOccupancyCut_FEDuTCA.fill(eid);
621  _cOccupancyCut_ElectronicsuTCA.fill(eid);
622  }
623  }
624  }
625  did.subdet() == HcalBarrel ? nChsHBCut++ : nChsHECut++;
626  }
627  did.subdet() == HcalBarrel ? nChsHB++ : nChsHE++;
628  }
629 
630  if (rawidHBValid != 0 && rawidHEValid != 0) {
631  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHBValid), _currentLS, nChsHB);
632  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHEValid), _currentLS, nChsHE);
633 
634  // ONLINE ONLY!
635  if (_ptype == fOnline) {
636  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHBValid), bx, nChsHBCut);
637  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHEValid), bx, nChsHECut);
638  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHBValid), _currentLS, nChsHBCut);
639  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHEValid), _currentLS, nChsHECut);
640  }
641  // ^^^ONLINE ONLY!
642  }
643 
644  // reset
645  rawidValid = 0;
646 
647  int nChsHO = 0;
648  int nChsHOCut = 0;
649  double ehop = 0;
650  double ehom = 0;
651  for (HORecHitCollection::const_iterator it = cho->begin(); it != cho->end(); ++it) {
652  double energy = it->energy();
653  double timing = it->time();
654 
655  // Explicit check on the DetIds present in the Collection
656  HcalDetId did = it->id();
657  uint32_t rawid = _ehashmap.lookup(did);
658  if (rawid == 0) {
659  meUnknownIds1LS->Fill(1);
660  _unknownIdsPresent = true;
661  continue;
662  }
663  HcalElectronicsId const& eid(rawid);
664  if (did.subdet() == HcalOuter)
665  rawidValid = did.rawId();
666 
667  _cEnergy_Subdet.fill(did, energy);
668  _cTimingvsEnergy_SubdetPM.fill(did, energy, timing);
669  _cOccupancy_depth.fill(did);
670  did.ieta() > 0 ? ehop += energy : ehom += energy;
671 
672  // IMPORTANT: ONLINE ONLY!
673  if (_ptype == fOnline) {
674  _cOccupancyvsiphi_SubdetPM.fill(did);
675  _cOccupancyvsieta_Subdet.fill(did);
676  }
677  // ONLINE ONLY!
678 
679  if (_ptype != fOffline) { // hidefed2crate
680  if (eid.isVMEid()) {
681  _cOccupancy_FEDVME.fill(eid);
682  _cOccupancy_ElectronicsVME.fill(eid);
683  } else {
684  _cOccupancy_FEDuTCA.fill(eid);
685  _cOccupancy_ElectronicsuTCA.fill(eid);
686  }
687  }
688 
689  if (energy > _cutE_HO) {
690  // ONLINE ONLY!
691  if (_ptype == fOnline) {
692  _cEnergyvsLS_SubdetPM.fill(did, _currentLS, energy);
693  _cEnergyvsBX_SubdetPM.fill(did, bx, energy);
694  _cEnergyvsieta_Subdet.fill(did, energy);
695  _cEnergyvsiphi_SubdetPM.fill(did, energy);
696  _cTimingCutvsieta_Subdet.fill(did, timing);
697  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
698  _cTimingCutvsBX_SubdetPM.fill(did, bx, timing);
699  _cOccupancyCutvsiphi_SubdetPM.fill(did);
700  _cOccupancyCutvsieta_Subdet.fill(did);
701  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
702  }
703  // ^^^ONLINE ONLY!
704 
705  _cEnergy_depth.fill(did, energy);
706  _cTimingCut_SubdetPM.fill(did, timing);
707  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
708  if (_ptype != fOffline) { // hidefed2crate
709  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
710  }
711  _cOccupancyCut_depth.fill(did);
712  _cTimingCut_depth.fill(did, timing);
713  if (_ptype != fOffline) { // hidefed2crate
714  if (eid.isVMEid()) {
715  _cTimingCut_FEDVME.fill(eid, timing);
716  _cTimingCut_ElectronicsVME.fill(eid, timing);
717  _cOccupancyCut_FEDVME.fill(eid);
718  _cOccupancyCut_ElectronicsVME.fill(eid);
719  } else {
720  _cTimingCut_FEDuTCA.fill(eid, timing);
721  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
722  _cOccupancyCut_FEDuTCA.fill(eid);
723  _cOccupancyCut_ElectronicsuTCA.fill(eid);
724  }
725  }
726  nChsHOCut++;
727  }
728  nChsHO++;
729  }
730 
731  if (rawidValid != 0) {
732  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS, nChsHO);
733  // ONLINE ONLY!
734  if (_ptype == fOnline) {
735  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid), bx, nChsHOCut);
736  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS, nChsHOCut);
737  }
738  // ^^^ONLINE ONLY!
739  }
740 
741  //reset
742  rawidValid = 0;
743 
744  int nChsHF = 0;
745  int nChsHFCut = 0;
746  double ehfp = 0;
747  double ehfm = 0;
748  for (HFRecHitCollection::const_iterator it = chf->begin(); it != chf->end(); ++it) {
749  double energy = it->energy();
750  double timing = it->time();
751 
752  // Explicit check on the DetIds present in the Collection
753  HcalDetId did = it->id();
754  uint32_t rawid = _ehashmap.lookup(did);
755  if (rawid == 0) {
756  meUnknownIds1LS->Fill(1);
757  _unknownIdsPresent = true;
758  continue;
759  }
760  HcalElectronicsId const& eid(rawid);
761  if (did.subdet() == HcalForward)
762  rawidValid = did.rawId();
763 
764  _cEnergy_Subdet.fill(did, energy);
765  _cTimingvsEnergy_SubdetPM.fill(did, energy, timing);
766  _cOccupancy_depth.fill(did);
767  did.ieta() > 0 ? ehfp += energy : ehfm += energy;
768 
769  // IMPORTANT:
770  // only for Online Processing
771  //
772  if (_ptype == fOnline) {
773  _cOccupancyvsiphi_SubdetPM.fill(did);
774  _cOccupancyvsieta_Subdet.fill(did);
775  }
776  // ONLINE ONLY!
777 
778  if (_ptype != fOffline) { // hidefed2crate
779  if (eid.isVMEid()) {
780  _cOccupancy_FEDVME.fill(eid);
781  _cOccupancy_ElectronicsVME.fill(eid);
782  } else {
783  _cOccupancy_FEDuTCA.fill(eid);
784  _cOccupancy_ElectronicsuTCA.fill(eid);
785  }
786  }
787 
788  if (energy > _cutE_HF) {
789  // ONLINE ONLY!
790  if (_ptype == fOnline) {
791  _cEnergyvsLS_SubdetPM.fill(did, _currentLS, energy);
792  _cEnergyvsBX_SubdetPM.fill(did, bx, energy);
793  _cEnergyvsieta_Subdet.fill(did, energy);
794  _cEnergyvsiphi_SubdetPM.fill(did, energy);
795  _cTimingCutvsieta_Subdet.fill(did, timing);
796  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
797  _cTimingCutvsBX_SubdetPM.fill(did, bx, timing);
798  _cOccupancyCutvsiphi_SubdetPM.fill(did);
799  _cOccupancyCutvsieta_Subdet.fill(did);
800  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
801  if (_ptype != fOffline) { // hidefed2crate
802  _xUniHF.get(eid)++;
803  }
804  }
805  // ^^^ONLINE ONLY!
806  _cEnergy_depth.fill(did, energy);
807  _cTimingCut_SubdetPM.fill(did, timing);
808  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
809  if (_ptype != fOffline) { // hidefed2crate
810  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
811  }
812  _cOccupancyCut_depth.fill(did);
813  _cTimingCut_depth.fill(did, timing);
814  if (_ptype != fOffline) { // hidefed2crate
815  if (eid.isVMEid()) {
816  _cTimingCut_FEDVME.fill(eid, timing);
817  _cTimingCut_ElectronicsVME.fill(eid, timing);
818  _cOccupancyCut_FEDVME.fill(eid);
819  _cOccupancyCut_ElectronicsVME.fill(eid);
820  } else {
821  _cTimingCut_FEDuTCA.fill(eid, timing);
822  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
823  _cOccupancyCut_FEDuTCA.fill(eid);
824  _cOccupancyCut_ElectronicsuTCA.fill(eid);
825  }
826  }
827  nChsHFCut++;
828  }
829  nChsHF++;
830  }
831 
832  if (rawidValid != 0) {
833  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS, nChsHF);
834  // ONLINE ONLY!
835  if (_ptype == fOnline) {
836  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid), bx, nChsHFCut);
837  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS, nChsHFCut);
838  }
839  // ^^^ONLINE ONLY!
840  }
841 
842  // Loop over HFPreRecHits to get charge and charge asymmetry
843  if (_hfPreRecHitsAvailable) {
844  for (HFPreRecHitCollection::const_iterator it = cprehf->begin(); it != cprehf->end(); ++it) {
845  HcalDetId did = it->id();
846  if (_filter_HF.filter(did)) {
847  continue;
848  }
849  std::pair<float, bool> chargeAsymmetry = it->chargeAsymmetry(0.);
850  std::pair<float, bool> chargeAsymmetryCut = it->chargeAsymmetry(20.);
851 
852  if (chargeAsymmetry.second) {
853  _cDAAsymmetryVsCharge_SubdetPM.fill(did, chargeAsymmetry.first, it->charge());
854  }
855  if (chargeAsymmetryCut.second) {
856  _cDAAsymmetryMean_cut_depth.fill(did, chargeAsymmetryCut.first);
857  _cDAAsymmetry_cut_SubdetPM.fill(did, chargeAsymmetryCut.first);
858  }
859  }
860  }
861 }
862 
863 /* virtual */ void RecHitTask::dqmBeginLuminosityBlock(edm::LuminosityBlock const& lb, edm::EventSetup const& es) {
864  DQTask::dqmBeginLuminosityBlock(lb, es);
865 }
866 
867 /* virtual */ void RecHitTask::dqmEndLuminosityBlock(edm::LuminosityBlock const& lb, edm::EventSetup const& es) {
868  if (_ptype != fOnline)
869  return;
870 
871  //
872  // GENERATE STATUS ONLY FOR ONLINE
873  //
874  // for (std::vector<HcalGenericDetId>::const_iterator it=gids.begin();
875  // it!=gids.end(); ++it)
876  // {}
877 
878  for (uintCompactMap::const_iterator it = _xUniHF.begin(); it != _xUniHF.end(); ++it) {
879  uint32_t hash1 = it->first;
880  HcalElectronicsId eid1(hash1);
881  double x1 = it->second;
882 
883  for (uintCompactMap::const_iterator jt = _xUniHF.begin(); jt != _xUniHF.end(); ++jt) {
884  if (jt == it)
885  continue;
886  double x2 = jt->second;
887  if (x2 == 0)
888  continue;
889  if (x1 / x2 < _thresh_unihf)
890  _xUni.get(eid1)++;
891  }
892  }
893 
894  if (_ptype != fOffline) { // hidefed2crate
895  for (std::vector<uint32_t>::const_iterator it = _vhashFEDs.begin(); it != _vhashFEDs.end(); ++it) {
896  flag::Flag fSum("RECO");
897  HcalElectronicsId eid = HcalElectronicsId(*it);
898 
899  std::vector<uint32_t>::const_iterator cit = std::find(_vcdaqEids.begin(), _vcdaqEids.end(), *it);
900  if (cit == _vcdaqEids.end()) {
901  // not @cDAQ
902  for (uint32_t iflag = 0; iflag < _vflags.size(); iflag++)
903  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag), int(flag::fNCDAQ));
904  _cSummaryvsLS.setBinContent(eid, _currentLS, int(flag::fNCDAQ));
905  continue;
906  }
907 
908  // FED is @cDAQ
909  if (hcaldqm::utilities::isFEDHF(eid) && (_runkeyVal == 0 || _runkeyVal == 4)) {
910  if (_xUni.get(eid) > 0)
911  _vflags[fUni]._state = flag::fPROBLEMATIC;
912  else
913  _vflags[fUni]._state = flag::fGOOD;
914  }
915 
916  if (_unknownIdsPresent)
917  _vflags[fUnknownIds]._state = flag::fBAD;
918  else
919  _vflags[fUnknownIds]._state = flag::fGOOD;
920 
921  int iflag = 0;
922  for (std::vector<flag::Flag>::iterator ft = _vflags.begin(); ft != _vflags.end(); ++ft) {
923  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag), int(ft->_state));
924  fSum += (*ft);
925  iflag++;
926 
927  // reset after using
928  ft->reset();
929  }
930  _cSummaryvsLS.setBinContent(eid, _currentLS, fSum._state);
931  }
932  _xUniHF.reset();
933  _xUni.reset();
934  }
935 
936  // in the end always do the DQTask::endLumi
937  DQTask::dqmEndLuminosityBlock(lb, es);
938 }
939 
940 DEFINE_FWK_MODULE(RecHitTask);
hcaldqm::ContainerProf1D::initialize
void initialize(std::string const &folder, hashfunctions::HashType, quantity::Quantity *, quantity::Quantity *, int debug=0) override
Definition: ContainerProf1D.cc:22
RecHitTask::_xUniHF
hcaldqm::ContainerXXX< uint32_t > _xUniHF
Definition: RecHitTask.h:115
HcalElectronicsId::rawId
constexpr uint32_t rawId() const
Definition: HcalElectronicsId.h:59
dqm::dqmstoreimpl::DQMStore::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX)
Definition: DQMStore.cc:239
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
RecHitTask::_cOccupancyCut_depth
hcaldqm::Container2D _cOccupancyCut_depth
Definition: RecHitTask.h:110
RecHitTask::_cSummaryvsLS_FED
hcaldqm::Container2D _cSummaryvsLS_FED
Definition: RecHitTask.h:127
RecHitTask::_vhashFEDs
std::vector< uint32_t > _vhashFEDs
Definition: RecHitTask.h:53
RecHitTask::_cOccupancyvsLS_Subdet
hcaldqm::ContainerProf1D _cOccupancyvsLS_Subdet
Definition: RecHitTask.h:100
HLT_2018_cff.InputTag
InputTag
Definition: HLT_2018_cff.py:78913
RecHitTask::_tokHBHE
edm::EDGetTokenT< HBHERecHitCollection > _tokHBHE
Definition: RecHitTask.h:44
RecHitTask::_cEnergy_depth
hcaldqm::ContainerProf2D _cEnergy_depth
Definition: RecHitTask.h:73
hcaldqm::Container1D::initialize
virtual void initialize(std::string const &folder, hashfunctions::HashType, quantity::Quantity *, quantity::Quantity *qy=new quantity::ValueQuantity(quantity::fN), int debug=0)
Definition: Container1D.cc:33
RecHitTask::_cTimingvsEnergy_SubdetPM
hcaldqm::Container2D _cTimingvsEnergy_SubdetPM
Definition: RecHitTask.h:78
hcaldqm::constants::CRATE_VME_MIN
int const CRATE_VME_MIN
Definition: Constants.h:83
RecHitTask::_cOccupancyCut_FEDVME
hcaldqm::Container2D _cOccupancyCut_FEDVME
Definition: RecHitTask.h:105
RecHitTask::_cOccupancy_FEDVME
hcaldqm::Container2D _cOccupancy_FEDVME
Definition: RecHitTask.h:96
HcalDetId::subdet
HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:138
RecHitTask::_cutE_HO
double _cutE_HO
Definition: RecHitTask.h:49
hcaldqm::ContainerProf2D::fill
void fill(HcalDetId const &) override
Definition: ContainerProf2D.cc:228
hcaldqm::utilities::fed2crate
std::pair< uint16_t, uint16_t > fed2crate(int fed)
Definition: Utilities.cc:10
RecHitTask::_tokHF
edm::EDGetTokenT< HFRecHitCollection > _tokHF
Definition: RecHitTask.h:46
RecHitTask::_xUni
hcaldqm::ContainerXXX< uint32_t > _xUni
Definition: RecHitTask.h:115
l1GtPatternGenerator_cfi.bx
bx
Definition: l1GtPatternGenerator_cfi.py:18
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:525
hcaldqm::electronicsmap::ElectronicsMap::initialize
void initialize(HcalElectronicsMap const *, ElectronicsMapType etype=fHcalElectronicsMap)
Definition: ElectronicsMap.cc:6
RecHitTask::_resetMonitors
void _resetMonitors(hcaldqm::UpdateFreq) override
Definition: RecHitTask.cc:462
RecHitTask::_cOccupancyCutvsiphi_SubdetPM
hcaldqm::Container1D _cOccupancyCutvsiphi_SubdetPM
Definition: RecHitTask.h:111
RecHitTask::_cEnergyvsieta_Subdet
hcaldqm::ContainerProf1D _cEnergyvsieta_Subdet
Definition: RecHitTask.h:71
RecHitTask::_cEnergyvsiphi_SubdetPM
hcaldqm::ContainerProf1D _cEnergyvsiphi_SubdetPM
Definition: RecHitTask.h:72
RecHitTask::_cOccupancyCutvsiphivsLS_SubdetPM
hcaldqm::Container2D _cOccupancyCutvsiphivsLS_SubdetPM
Definition: RecHitTask.h:114
dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:418
RecHitTask::_gids
std::vector< HcalGenericDetId > _gids
Definition: RecHitTask.h:126
hcaldqm::ContainerXXX::begin
virtual CompactMap::const_iterator begin()
Definition: ContainerXXX.h:69
hcaldqm::UpdateFreq
UpdateFreq
Definition: DQTask.h:16
hcaldqm::flag::Flag::_state
State _state
Definition: Flag.h:62
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:82
dqm::impl::MonitorElement::setLumiFlag
void setLumiFlag()
this ME is meant to be stored for each luminosity section
Definition: MonitorElement.h:225
RecHitTask::_cTimingCut_SubdetPM
hcaldqm::Container1D _cTimingCut_SubdetPM
Definition: RecHitTask.h:81
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:64
RecHitTask::RecHitTask
RecHitTask(edm::ParameterSet const &)
Definition: RecHitTask.cc:8
hcaldqm::ContainerSingle2D::initialize
virtual void initialize(std::string const &folder, quantity::Quantity *, quantity::Quantity *, quantity::Quantity *qz=new quantity::ValueQuantity(quantity::fN), int debug=0)
Definition: ContainerSingle2D.cc:46
hcaldqm::constants::SPIGOT_MIN
int const SPIGOT_MIN
Definition: Constants.h:108
hcaldqm::ContainerProf1D::book
void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="") override
Definition: ContainerProf1D.cc:40
hcaldqm::ContainerXXX::reset
virtual void reset()
Definition: ContainerXXX.h:292
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
RecHitTask::_cTimingCutvsieta_Subdet
hcaldqm::ContainerProf1D _cTimingCutvsieta_Subdet
Definition: RecHitTask.h:90
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
hcaldqm::constants::FIBER_VME_MIN
int const FIBER_VME_MIN
Definition: Constants.h:113
RecHitTask::_cEnergyvsBX_SubdetPM
hcaldqm::ContainerProf1D _cEnergyvsBX_SubdetPM
Definition: RecHitTask.h:75
RecHitTask::_cOccupancyCut_ElectronicsVME
hcaldqm::Container2D _cOccupancyCut_ElectronicsVME
Definition: RecHitTask.h:107
hcaldqm::constants::FIBERCH_MIN
int const FIBERCH_MIN
Definition: Constants.h:122
hcaldqm::ContainerSingle2D::book
virtual void book(DQMStore::IBooker &, std::string subsystem="Hcal", std::string aux="")
Definition: ContainerSingle2D.cc:68
hcaldqm::utilities::getFEDVMEList
std::vector< int > getFEDVMEList(HcalElectronicsMap const *)
Definition: Utilities.cc:101
RecHitTask::_cTimingCut_ElectronicsuTCA
hcaldqm::ContainerProf2D _cTimingCut_ElectronicsuTCA
Definition: RecHitTask.h:86
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:242
hcaldqm::ContainerXXX::end
virtual CompactMap::const_iterator end()
Definition: ContainerXXX.h:70
HcalElectronicsId::isUTCAid
constexpr bool isUTCAid() const
Definition: HcalElectronicsId.h:62
RecHitTask::_vflags
std::vector< hcaldqm::flag::Flag > _vflags
Definition: RecHitTask.h:56
RecHitTask::_cutE_HBHE
double _cutE_HBHE
Definition: RecHitTask.h:49
hcaldqm::DQModule::_ptype
ProcessingType _ptype
Definition: DQModule.h:43
hcaldqm::Container1D::fill
virtual void fill(uint32_t)
Definition: Container1D.cc:73
RecHitTask::_cEnergy_Subdet
hcaldqm::Container1D _cEnergy_Subdet
Definition: RecHitTask.h:70
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
RecHitTask::_cTimingCutvsLS_FED
hcaldqm::ContainerProf1D _cTimingCutvsLS_FED
Definition: RecHitTask.h:88
RecHitTask::_cTimingCut_FEDVME
hcaldqm::ContainerProf2D _cTimingCut_FEDVME
Definition: RecHitTask.h:83
RecHitTask::_cTimingCutvsLS_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsLS_SubdetPM
Definition: RecHitTask.h:89
hcaldqm::flag::Flag::reset
void reset()
Definition: Flag.h:59
RecHitTask::_cOccupancyCut_ElectronicsuTCA
hcaldqm::Container2D _cOccupancyCut_ElectronicsuTCA
Definition: RecHitTask.h:108
hcaldqm::ContainerSingle2D::setBinContent
virtual void setBinContent(int, int, int)
Definition: ContainerSingle2D.cc:154
RecHitTask::_cOccupancyCutvsBX_Subdet
hcaldqm::ContainerProf1D _cOccupancyCutvsBX_Subdet
Definition: RecHitTask.h:113
HcalElectronicsId::isVMEid
constexpr bool isVMEid() const
Definition: HcalElectronicsId.h:61
Logger::dqmthrow
void dqmthrow(std::string const &msg) const
Definition: Logger.h:12
testProducerWithPsetDescEmpty_cfi.x2
x2
Definition: testProducerWithPsetDescEmpty_cfi.py:28
hcaldqm::constants::FED_VME_MIN
int const FED_VME_MIN
Definition: Constants.h:70
RecHitTask::_tagPreHF
edm::InputTag _tagPreHF
Definition: RecHitTask.h:42
HcalDetId::ieta
int ieta() const
get the cell ieta
Definition: HcalDetId.h:155
RecHitTask::_cOccupancy_ElectronicsVME
hcaldqm::Container2D _cOccupancy_ElectronicsVME
Definition: RecHitTask.h:98
hcaldqm::constants::FIBER_uTCA_MIN1
int const FIBER_uTCA_MIN1
Definition: Constants.h:116
RecHitTask::_cutE_HF
double _cutE_HF
Definition: RecHitTask.h:49
ecalMatacq_cfi.timing
timing
Definition: ecalMatacq_cfi.py:26
RecHitTask::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: RecHitTask.cc:32
hcaldqm::hashfunctions::hash_did
hash_function_did const hash_did[nHashType_did]
Definition: HashFunctions.h:163
RecHitTask::dqmEndLuminosityBlock
void dqmEndLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &) override
Definition: RecHitTask.cc:867
HcalElectronicsMap::allPrecisionId
std::vector< HcalGenericDetId > allPrecisionId() const
Definition: HcalElectronicsMap.cc:139
hcaldqm::DQModule::_logger
Logger _logger
Definition: DQModule.h:54
RecHitTask::_cOccupancyCut_FEDuTCA
hcaldqm::Container2D _cOccupancyCut_FEDuTCA
Definition: RecHitTask.h:106
RecHitTask::_ehashmap
hcaldqm::electronicsmap::ElectronicsMap _ehashmap
Definition: RecHitTask.h:60
RecHitTask::meUnknownIds1LS
MonitorElement * meUnknownIds1LS
Definition: RecHitTask.h:123
RecHitTask::_tokPreHF
edm::EDGetTokenT< HFPreRecHitCollection > _tokPreHF
Definition: RecHitTask.h:47
hcaldqm::DQModule::_name
std::string _name
Definition: DQModule.h:41
hcaldqm::DQTask::_emap
HcalElectronicsMap const * _emap
Definition: DQTask.h:59
RecHitTask::_cSummaryvsLS
hcaldqm::ContainerSingle2D _cSummaryvsLS
Definition: RecHitTask.h:128
RecHitTask::_thresh_unihf
double _thresh_unihf
Definition: RecHitTask.h:50
hcaldqm::constants::CRATE_uTCA_MIN
int const CRATE_uTCA_MIN
Definition: Constants.h:88
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:269
HCALHighEnergyHPDFilter_cfi.energy
energy
Definition: HCALHighEnergyHPDFilter_cfi.py:5
RecHitTask::_cOccupancy_ElectronicsuTCA
hcaldqm::Container2D _cOccupancy_ElectronicsuTCA
Definition: RecHitTask.h:99
hcaldqm::ContainerXXX::book
virtual void book(HcalElectronicsMap const *)
Definition: ContainerXXX.h:87
RecHitTask::_tagHO
edm::InputTag _tagHO
Definition: RecHitTask.h:40
hcaldqm::Container1D::book
virtual void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="")
Definition: Container1D.cc:592
RecHitTask::_filter_HF
hcaldqm::filter::HashFilter _filter_HF
Definition: RecHitTask.h:67
hcaldqm::Container2D::initialize
virtual void initialize(std::string const &folder, hashfunctions::HashType, quantity::Quantity *, quantity::Quantity *, quantity::Quantity *qz=new quantity::ValueQuantity(quantity::fN), int debug=0)
Definition: Container2D.cc:25
RecHitTask::dqmBeginLuminosityBlock
void dqmBeginLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &) override
Definition: RecHitTask.cc:863
RecHitTask::_cOccupancy_FEDuTCA
hcaldqm::Container2D _cOccupancy_FEDuTCA
Definition: RecHitTask.h:97
hcaldqm::ContainerProf2D::book
void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="") override
Definition: ContainerProf2D.cc:47
bookHistograms
example_stream void bookHistograms(DQMStore::IBooker &,@example_stream edm::Run const &,@example_stream edm::EventSetup const &) override
RecHitTask::_tagHF
edm::InputTag _tagHF
Definition: RecHitTask.h:41
hcaldqm::ContainerProf2D::initialize
void initialize(std::string const &folder, hashfunctions::HashType, quantity::Quantity *, quantity::Quantity *, quantity::Quantity *qz=new quantity::ValueQuantity(quantity::fEnergy), int debug=0) override
Definition: ContainerProf2D.cc:22
RecHitTask::_cDAAsymmetryVsCharge_SubdetPM
hcaldqm::Container2D _cDAAsymmetryVsCharge_SubdetPM
Definition: RecHitTask.h:118
hcaldqm::Container2D::setBinContent
void setBinContent(HcalDetId const &, int) override
Definition: Container2D.cc:163
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
RecHitTask::_cTimingCut_HBHEPartition
hcaldqm::Container1D _cTimingCut_HBHEPartition
Definition: RecHitTask.h:82
RecHitTask::_filter_uTCA
hcaldqm::filter::HashFilter _filter_uTCA
Definition: RecHitTask.h:64
RecHitTask::_cOccupancyCutvsLS_Subdet
hcaldqm::ContainerProf1D _cOccupancyCutvsLS_Subdet
Definition: RecHitTask.h:109
hcaldqm::ContainerXXX::get
virtual STDTYPE & get(HcalDetId const &)
Definition: ContainerXXX.h:197
RecHitTask::_cOccupancyvsieta_Subdet
hcaldqm::Container1D _cOccupancyvsieta_Subdet
Definition: RecHitTask.h:102
hcaldqm::utilities::getFEDList
std::vector< int > getFEDList(HcalElectronicsMap const *)
Definition: Utilities.cc:83
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
hcaldqm::utilities::getFEDuTCAList
std::vector< int > getFEDuTCAList(HcalElectronicsMap const *)
Definition: Utilities.cc:121
hcaldqm::Container2D::fill
void fill(HcalDetId const &) override
Definition: Container2D.cc:52
HcalDbService::getHcalMapping
const HcalElectronicsMap * getHcalMapping() const
Definition: HcalDbService.cc:343
hcaldqm::filter::HashFilter::filter
virtual bool filter(HcalDetId const &) const
Definition: HashFilter.cc:24
RecHitTask::_hfPreRecHitsAvailable
bool _hfPreRecHitsAvailable
Definition: RecHitTask.h:43
hcaldqm::filter::HashFilter::initialize
virtual void initialize(FilterType ftype, hashfunctions::HashType htype, std::vector< uint32_t > const &)
Definition: HashFilter.cc:17
hcaldqm::DQModule::_subsystem
std::string _subsystem
Definition: DQModule.h:48
RecHitTask::_cTimingCutvsBX_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsBX_SubdetPM
Definition: RecHitTask.h:92
RecHitTask::_tokHO
edm::EDGetTokenT< HORecHitCollection > _tokHO
Definition: RecHitTask.h:45
hcaldqm::DQTask::_vcdaqEids
std::vector< uint32_t > _vcdaqEids
Definition: DQTask.h:51
RecHitTask::_process
void _process(edm::Event const &, edm::EventSetup const &) override
Definition: RecHitTask.cc:474
RecHitTask::_tagHBHE
edm::InputTag _tagHBHE
Definition: RecHitTask.h:39
RecHitTask::_cTimingCut_FEDuTCA
hcaldqm::ContainerProf2D _cTimingCut_FEDuTCA
Definition: RecHitTask.h:84
hcaldqm::ContainerXXX::initialize
virtual void initialize(hashfunctions::HashType, int debug=0)
Definition: ContainerXXX.h:81
hcaldqm::Container2D::book
void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="") override
Definition: Container2D.cc:558
RecHitTask::_cOccupancy_depth
hcaldqm::Container2D _cOccupancy_depth
Definition: RecHitTask.h:95
RecHitTask::_cDAAsymmetryMean_cut_depth
hcaldqm::ContainerProf2D _cDAAsymmetryMean_cut_depth
Definition: RecHitTask.h:119
RecHitTask::_cTimingCutvsiphi_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsiphi_SubdetPM
Definition: RecHitTask.h:91
testProducerWithPsetDescEmpty_cfi.x1
x1
Definition: testProducerWithPsetDescEmpty_cfi.py:33
RecHitTask::_cTimingCut_depth
hcaldqm::ContainerProf2D _cTimingCut_depth
Definition: RecHitTask.h:87
HcalElectronicsId
Readout chain identification for Hcal.
Definition: HcalElectronicsId.h:32
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:37
RecHitTask::_cTimingCut_ElectronicsVME
hcaldqm::ContainerProf2D _cTimingCut_ElectronicsVME
Definition: RecHitTask.h:85
RecHitTask::_filter_VME
hcaldqm::filter::HashFilter _filter_VME
Definition: RecHitTask.h:63
nanoDQM_cff.Flag
Flag
Definition: nanoDQM_cff.py:65
RecHitTask::_cOccupancyCutvsieta_Subdet
hcaldqm::Container1D _cOccupancyCutvsieta_Subdet
Definition: RecHitTask.h:112
hcaldqm::utilities::isFEDHF
bool isFEDHF(HcalElectronicsId const &)
Definition: Utilities.cc:156
RecHitTask::_cOccupancyvsiphi_SubdetPM
hcaldqm::Container1D _cOccupancyvsiphi_SubdetPM
Definition: RecHitTask.h:101
hcaldqm::constants::SLOT_uTCA_MIN
int const SLOT_uTCA_MIN
Definition: Constants.h:95
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:264
edm::Run
Definition: Run.h:45
RecHitTask::_unknownIdsPresent
bool _unknownIdsPresent
Definition: RecHitTask.h:124
cuy.ib
ib
Definition: cuy.py:662
RecHitTask::_cDAAsymmetry_cut_SubdetPM
hcaldqm::Container1D _cDAAsymmetry_cut_SubdetPM
Definition: RecHitTask.h:120
RecHitTask::_cEnergyvsLS_SubdetPM
hcaldqm::ContainerProf1D _cEnergyvsLS_SubdetPM
Definition: RecHitTask.h:74
Generated for CMSSW Reference Manual by doxygen 1.8.11