CMS 3D CMS Logo
DigiTask.cc
Go to the documentation of this file.
1 #include "DQM/HcalTasks/interface/DigiTask.h"
2 
3 using namespace hcaldqm;
4 using namespace hcaldqm::constants;
5 using namespace hcaldqm::filter;
6 
7 DigiTask::DigiTask(edm::ParameterSet const& ps) : DQTask(ps) {
8  _tagQIE11 = ps.getUntrackedParameter<edm::InputTag>("tagHE",
9  edm::InputTag("hcalDigis"));
10  _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO",
11  edm::InputTag("hcalDigis"));
12  _tagQIE10 = ps.getUntrackedParameter<edm::InputTag>("tagHF",
13  edm::InputTag("hcalDigis"));
14 
15  _tokQIE11 = consumes<QIE11DigiCollection>(_tagQIE11);
16  _tokHO = consumes<HODigiCollection>(_tagHO);
17  _tokQIE10 = consumes<QIE10DigiCollection>(_tagQIE10);
18 
19  _cutSumQ_HBHE = ps.getUntrackedParameter<double>("cutSumQ_HBHE", 20);
20  _cutSumQ_HO = ps.getUntrackedParameter<double>("cutSumQ_HO", 20);
21  _cutSumQ_HF = ps.getUntrackedParameter<double>("cutSumQ_HF", 20);
22  _thresh_unihf = ps.getUntrackedParameter<double>("thresh_unihf", 0.2);
23  _thresh_led = ps.getUntrackedParameter<double>("thresh_led", 20);
24 
25  _vflags.resize(nDigiFlag);
26  _vflags[fUni] = hcaldqm::flag::Flag("UniSlotHF");
27  _vflags[fDigiSize] = hcaldqm::flag::Flag("DigiSize");
28  _vflags[fNChsHF] = hcaldqm::flag::Flag("NChsHF");
29  _vflags[fUnknownIds] = hcaldqm::flag::Flag("UnknownIds");
30  _vflags[fLED] = hcaldqm::flag::Flag("LEDMisfire");
31  _vflags[fCapId] = hcaldqm::flag::Flag("BadCapId");
32 
33  _qie10InConditions = ps.getUntrackedParameter<bool>("qie10InConditions", true);
34 
35  // Get reference digi sizes. Convert from unsigned to signed int, because <digi>::size()/samples() return ints for some reason.
36  std::vector<uint32_t> vrefDigiSize = ps.getUntrackedParameter<std::vector<uint32_t>>("refDigiSize");
37  _refDigiSize[HcalBarrel] = (int)vrefDigiSize[0];
38  _refDigiSize[HcalEndcap] = (int)vrefDigiSize[1];
39  _refDigiSize[HcalOuter] = (int)vrefDigiSize[2];
40  _refDigiSize[HcalForward] = (int)vrefDigiSize[3];
41 
42  // (capid - BX) % 4 to 1
43  _capidmbx[HcalBarrel] = 1;
44  _capidmbx[HcalEndcap] = 1;
45  _capidmbx[HcalOuter] = 1;
46  _capidmbx[HcalForward] = 1;
47 
48  // LED calibration channels
49  std::vector<edm::ParameterSet> vLedCalibChannels =
50  ps.getParameter<std::vector<edm::ParameterSet>>("ledCalibrationChannels");
51  for (int i = 0; i <= 3; ++i) {
52  HcalSubdetector this_subdet = HcalEmpty;
53  switch (i) {
54  case 0:
55  this_subdet = HcalBarrel;
56  break;
57  case 1:
58  this_subdet = HcalEndcap;
59  break;
60  case 2:
61  this_subdet = HcalOuter;
62  break;
63  case 3:
64  this_subdet = HcalForward;
65  break;
66  default:
67  this_subdet = HcalEmpty;
68  break;
69  }
70  std::vector<int32_t> subdet_calib_ietas = vLedCalibChannels[i].getUntrackedParameter<std::vector<int32_t>>("ieta");
71  std::vector<int32_t> subdet_calib_iphis = vLedCalibChannels[i].getUntrackedParameter<std::vector<int32_t>>("iphi");
72  std::vector<int32_t> subdet_calib_depths =
73  vLedCalibChannels[i].getUntrackedParameter<std::vector<int32_t>>("depth");
74  for (unsigned int ichannel = 0; ichannel < subdet_calib_ietas.size(); ++ichannel) {
75  _ledCalibrationChannels[this_subdet].push_back(HcalDetId(
76  HcalOther, subdet_calib_ietas[ichannel], subdet_calib_iphis[ichannel], subdet_calib_depths[ichannel]));
77  }
78  }
79 }
80 
81 /* virtual */ void DigiTask::bookHistograms(DQMStore::IBooker& ib, edm::Run const& r, edm::EventSetup const& es) {
82  DQTask::bookHistograms(ib, r, es);
83 
84  // GET WHAT YOU NEED
85  edm::ESHandle<HcalDbService> dbs;
86  es.get<HcalDbRecord>().get(dbs);
87  _emap = dbs->getHcalMapping();
88  std::vector<uint32_t> vVME;
89  std::vector<uint32_t> vuTCA;
90  vVME.push_back(
91  HcalElectronicsId(constants::FIBERCH_MIN, constants::FIBER_VME_MIN, SPIGOT_MIN, CRATE_VME_MIN).rawId());
92  vuTCA.push_back(HcalElectronicsId(CRATE_uTCA_MIN, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
93  _filter_VME.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics, vVME);
94  _filter_uTCA.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics, vuTCA);
95 
96  // Filters for QIE8 vs QIE10/11
97  std::vector<uint32_t> vhashQIE8;
98  vhashQIE8.push_back(hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdet](HcalDetId(HcalOuter, 1,1,4)));
99  _filter_QIE8.initialize(filter::fPreserver, hcaldqm::hashfunctions::fSubdet,
100  vhashQIE8);
101 
102  std::vector<uint32_t> vhashQIE1011;
103  vhashQIE1011.push_back(hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdet](HcalDetId(HcalBarrel, 1,1,1)));
104  vhashQIE1011.push_back(hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdet](HcalDetId(HcalEndcap, 20,1,1)));
105  vhashQIE1011.push_back(hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdet](HcalDetId(HcalForward, 29,1,1)));
106  _filter_QIE1011.initialize(filter::fPreserver, hcaldqm::hashfunctions::fSubdet,
107  vhashQIE1011);
108 
109  std::vector<uint32_t> vhash_TDC2bit;
110  vhash_TDC2bit.push_back(
111  hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdetPM](HcalDetId(HcalBarrel, 1, 1, 1)));
112  vhash_TDC2bit.push_back(
113  hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdetPM](HcalDetId(HcalBarrel, -11, 1, 1)));
114  _filter_TDC2bit.initialize(filter::fPreserver, hcaldqm::hashfunctions::fSubdetPM, vhash_TDC2bit);
115 
116  std::vector<uint32_t> vhash_TDC6bit;
117  vhash_TDC6bit.push_back(
118  hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdetPM](HcalDetId(HcalEndcap, 20, 1, 1)));
119  vhash_TDC6bit.push_back(
120  hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdetPM](HcalDetId(HcalEndcap, -20, 1, 1)));
121  vhash_TDC6bit.push_back(
122  hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdetPM](HcalDetId(HcalForward, 29, 1, 1)));
123  vhash_TDC6bit.push_back(
124  hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdetPM](HcalDetId(HcalForward, -29, 1, 1)));
125  _filter_TDC6bit.initialize(filter::fPreserver, hcaldqm::hashfunctions::fSubdetPM, vhash_TDC6bit);
126 
127  // INITIALIZE FIRST
128  _cADC_SubdetPM.initialize(_name,
129  "ADC",
130  hcaldqm::hashfunctions::fSubdetPM,
131  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_128),
132  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
133  0);
134  _cfC_SubdetPM.initialize(_name,
135  "fC",
136  hcaldqm::hashfunctions::fSubdetPM,
137  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),
138  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
139  0);
140  _cSumQ_SubdetPM.initialize(_name,
141  "SumQ",
142  hcaldqm::hashfunctions::fSubdetPM,
143  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),
144  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
145  0);
146  _cSumQ_depth.initialize(_name,
147  "SumQ",
148  hcaldqm::hashfunctions::fdepth,
149  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
150  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
151  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_400000, true),
152  0);
153  _cSumQvsLS_SubdetPM.initialize(_name,
154  "SumQvsLS",
155  hcaldqm::hashfunctions::fSubdetPM,
156  new hcaldqm::quantity::LumiSection(_maxLS),
157  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),
158  0);
159 
160  _cADC_SubdetPM_QIE1011.initialize(_name,
161  "ADC",
162  hcaldqm::hashfunctions::fSubdetPM,
163  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
164  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
165  0);
166  _cfC_SubdetPM_QIE1011.initialize(_name,
167  "fC",
168  hcaldqm::hashfunctions::fSubdetPM,
169  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000),
170  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
171  0);
172  _cSumQ_SubdetPM_QIE1011.initialize(_name,
173  "SumQ",
174  hcaldqm::hashfunctions::fSubdetPM,
175  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000),
176  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
177  0);
178  _cSumQvsLS_SubdetPM_QIE1011.initialize(_name,
179  "SumQvsLS",
180  hcaldqm::hashfunctions::fSubdetPM,
181  new hcaldqm::quantity::LumiSection(_maxLS),
182  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000),
183  0);
184 
185  _cTimingCut_SubdetPM.initialize(_name,
186  "TimingCut",
187  hcaldqm::hashfunctions::fSubdetPM,
188  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
189  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
190  0);
191  _cTimingCut_depth.initialize(_name,
192  "TimingCut",
193  hcaldqm::hashfunctions::fdepth,
194  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
195  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
196  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
197  0);
198  _cTimingCutvsLS_SubdetPM.initialize(_name,
199  "TimingvsLS",
200  hcaldqm::hashfunctions::fSubdetPM,
201  new hcaldqm::quantity::LumiSection(_maxLS),
202  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
203  0);
204 
205  // Occupancy w/o a cut
206  _cOccupancyvsLS_Subdet.initialize(_name,
207  "OccupancyvsLS",
208  hcaldqm::hashfunctions::fSubdet,
209  new hcaldqm::quantity::LumiSection(_maxLS),
210  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),
211  0);
212  _cOccupancy_depth.initialize(_name,
213  "Occupancy",
214  hcaldqm::hashfunctions::fdepth,
215  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
216  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
217  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
218  0);
219 
220  // Occupancy w/ a cut
221  _cOccupancyCutvsLS_Subdet.initialize(_name,
222  "OccupancyCutvsLS",
223  hcaldqm::hashfunctions::fSubdet,
224  new hcaldqm::quantity::LumiSection(_maxLS),
225  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),
226  0);
227  _cOccupancyCut_depth.initialize(_name,
228  "OccupancyCut",
229  hcaldqm::hashfunctions::fdepth,
230  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
231  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
232  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
233  0);
234 
235  // Digi size
236  _cDigiSize_Crate.initialize(_name,
237  "DigiSize",
238  hcaldqm::hashfunctions::fCrate,
239  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDigiSize),
240  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
241  0);
242  _cADCvsTS_SubdetPM.initialize(_name,
243  "ADCvsTS",
244  hcaldqm::hashfunctions::fSubdetPM,
245  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
246  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_128),
247  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
248  0);
249  _cADCvsTS_SubdetPM_QIE1011.initialize(_name,
250  "ADCvsTS",
251  hcaldqm::hashfunctions::fSubdetPM,
252  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
253  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
254  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
255  0);
256 
257  _cLETDCTimevsADC_SubdetPM.initialize(_name,
258  "LETDCTimevsADC",
259  hcaldqm::hashfunctions::fSubdetPM,
260  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
261  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250_coarse),
262  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
263  _cLETDCvsADC_2bit_SubdetPM.initialize(_name,
264  "LETDCvsADC",
265  hcaldqm::hashfunctions::fSubdetPM,
266  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
267  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_4),
268  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
269  _cLETDCvsADC_6bit_SubdetPM.initialize(_name,
270  "LETDCvsADC",
271  hcaldqm::hashfunctions::fSubdetPM,
272  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
273  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
274  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
275  _cLETDCvsTS_2bit_SubdetPM.initialize(_name,
276  "LETDCvsTS",
277  hcaldqm::hashfunctions::fSubdetPM,
278  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
279  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_4),
280  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
281  _cLETDCvsTS_6bit_SubdetPM.initialize(_name,
282  "LETDCvsTS",
283  hcaldqm::hashfunctions::fSubdetPM,
284  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
285  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
286  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
287 
288  _cLETDCTime_SubdetPM.initialize(_name,
289  "LETDCTime",
290  hcaldqm::hashfunctions::fSubdetPM,
291  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250_coarse),
292  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
293  _cLETDCTime_depth.initialize(_name,
294  "LETDCTime",
295  hcaldqm::hashfunctions::fdepth,
296  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
297  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
298  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
299  0);
300 
301  _cBadTDCValues_SubdetPM.initialize(_name,
302  "BadTDCValues",
303  hcaldqm::hashfunctions::fSubdetPM,
304  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBadTDC),
305  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
306  _cBadTDCvsBX_SubdetPM.initialize(_name,
307  "BadTDCvsBX",
308  hcaldqm::hashfunctions::fSubdetPM,
309  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
310  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
311  _cBadTDCvsLS_SubdetPM.initialize(_name,
312  "BadTDCvsLS",
313  hcaldqm::hashfunctions::fSubdetPM,
314  new hcaldqm::quantity::LumiSection(_maxLS),
315  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
316  _cBadTDCCount_depth.initialize(_name,
317  "BadTDCCount",
318  hcaldqm::hashfunctions::fdepth,
319  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
320  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
321  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
322  0);
323 
324  if (_ptype == fOnline || _ptype == fLocal) {
325  _cOccupancy_Crate.initialize(_name,
326  "Occupancy",
327  hashfunctions::fCrate,
328  new quantity::ElectronicsQuantity(quantity::fSlotuTCA),
329  new quantity::ElectronicsQuantity(quantity::fFiberuTCAFiberCh),
330  new quantity::ValueQuantity(quantity::fN),
331  0);
332  _cOccupancy_CrateSlot.initialize(_name,
333  "Occupancy",
334  hashfunctions::fCrateSlot,
335  new quantity::ElectronicsQuantity(quantity::fFiberuTCA),
336  new quantity::ElectronicsQuantity(quantity::fFiberCh),
337  new quantity::ValueQuantity(quantity::fN),
338  0);
339  }
340 
341  // INITIALIZE HISTOGRAMS that are only for Online
342  if (_ptype == fOnline) {
343  // Charge sharing
344  _cQ2Q12CutvsLS_FEDHF.initialize(_name,
345  "Q2Q12vsLS",
346  hcaldqm::hashfunctions::fFED,
347  new hcaldqm::quantity::LumiSection(_maxLS),
348  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fRatio_0to2),
349  0);
350  _cSumQvsBX_SubdetPM.initialize(_name,
351  "SumQvsBX",
352  hcaldqm::hashfunctions::fSubdetPM,
353  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
354  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),
355  0);
356  _cSumQvsBX_SubdetPM_QIE1011.initialize(_name,
357  "SumQvsBX",
358  hcaldqm::hashfunctions::fSubdetPM,
359  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
360  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_10000),
361  0);
362  _cDigiSizevsLS_FED.initialize(_name,
363  "DigiSizevsLS",
364  hcaldqm::hashfunctions::fFED,
365  new hcaldqm::quantity::LumiSection(_maxLS),
366  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDigiSize),
367  0);
368  _cTimingCutvsiphi_SubdetPM.initialize(_name,
369  "TimingCutvsiphi",
370  hcaldqm::hashfunctions::fSubdetPM,
371  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
372  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
373  0);
374  _cTimingCutvsieta_Subdet.initialize(_name,
375  "TimingCutvsieta",
376  hcaldqm::hashfunctions::fSubdet,
377  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
378  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
379  0);
380  _cOccupancyvsiphi_SubdetPM.initialize(_name,
381  "Occupancyvsiphi",
382  hcaldqm::hashfunctions::fSubdetPM,
383  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
384  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
385  0);
386  _cOccupancyvsieta_Subdet.initialize(_name,
387  "Occupancyvsieta",
388  hcaldqm::hashfunctions::fSubdet,
389  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
390  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
391  0);
392  _cOccupancyCutvsiphi_SubdetPM.initialize(_name,
393  "OccupancyCutvsiphi",
394  hcaldqm::hashfunctions::fSubdetPM,
395  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
396  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
397  0);
398  _cOccupancyCutvsieta_Subdet.initialize(_name,
399  "OccupancyCutvsieta",
400  hcaldqm::hashfunctions::fSubdet,
401  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
402  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
403  0);
404  _cOccupancyCutvsLS_Subdet.initialize(_name,
405  "OccupancyCutvsLS",
406  hcaldqm::hashfunctions::fSubdet,
407  new hcaldqm::quantity::LumiSection(_maxLS),
408  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),
409  0);
410  _cOccupancyCutvsBX_Subdet.initialize(_name,
411  "OccupancyCutvsBX",
412  hcaldqm::hashfunctions::fSubdet,
413  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
414  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),
415  0);
416  // _cOccupancyCutvsSlotvsLS_HFPM.initialize(_name,
417  // "OccupancyCutvsSlotvsLS", hcaldqm::hashfunctions::fSubdetPM,
418  // new hcaldqm::quantity::LumiSection(_maxLS),
419  // new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
420  // new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
421  _cOccupancyCutvsiphivsLS_SubdetPM.initialize(_name,
422  "OccupancyCutvsiphivsLS",
423  hcaldqm::hashfunctions::fSubdetPM,
424  new hcaldqm::quantity::LumiSection(_maxLS),
425  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
426  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
427  0);
428  }
429  _cCapidMinusBXmod4_SubdetPM.initialize(_name,
430  "CapID",
431  hcaldqm::hashfunctions::fSubdetPM,
432  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fCapidMinusBXmod4),
433  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
434 
435  for (int i = 0; i < 4; ++i) {
436  _cCapidMinusBXmod4_CrateSlotuTCA[i].initialize(_name,
437  "CapID",
438  new quantity::ElectronicsQuantity(quantity::fCrateuTCA),
439  new quantity::ElectronicsQuantity(quantity::fSlotuTCA),
440  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
441  0);
442  _cCapidMinusBXmod4_CrateSlotVME[i].initialize(_name,
443  "CapID",
444  new quantity::ElectronicsQuantity(quantity::fCrateVME),
445  new quantity::ElectronicsQuantity(quantity::fSlotVME),
446  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
447  0);
448  }
449 
450  if (_ptype != fOffline) { // hidefed2crate
451  std::vector<int> vFEDs = hcaldqm::utilities::getFEDList(_emap);
452  std::vector<int> vFEDsVME = hcaldqm::utilities::getFEDVMEList(_emap);
453  std::vector<int> vFEDsuTCA = hcaldqm::utilities::getFEDuTCAList(_emap);
454 
455  if (_ptype == fOnline) {
456  _cCapid_BadvsFEDvsLS.initialize(_name,
457  "CapID",
458  new hcaldqm::quantity::LumiSectionCoarse(_maxLS, 10),
459  new hcaldqm::quantity::FEDQuantity(vFEDs),
460  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
461  0);
462 
463  _cCapid_BadvsFEDvsLSmod60.initialize(_name,
464  "CapID",
465  new hcaldqm::quantity::LumiSection(60),
466  new hcaldqm::quantity::FEDQuantity(vFEDs),
467  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
468  0);
469  }
470 
471  std::vector<uint32_t> vFEDHF;
472  vFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
473  vFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN + 6, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
474  vFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
475  vFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN + 6, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
476  vFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
477  vFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN + 6, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
478 
479  // initialize filters
480  _filter_FEDHF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fFED, vFEDHF);
481 
482  // push the rawIds of each fed into the vector...
483  for (std::vector<int>::const_iterator it = vFEDsVME.begin(); it != vFEDsVME.end(); ++it)
484  _vhashFEDs.push_back(
485  HcalElectronicsId(constants::FIBERCH_MIN, FIBER_VME_MIN, SPIGOT_MIN, (*it) - FED_VME_MIN).rawId());
486  for (std::vector<int>::const_iterator it = vFEDsuTCA.begin(); it != vFEDsuTCA.end(); ++it) {
487  std::pair<uint16_t, uint16_t> cspair = utilities::fed2crate(*it);
488  _vhashFEDs.push_back(HcalElectronicsId(cspair.first, cspair.second, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
489  }
490 
491  _cShapeCut_FED.initialize(_name,
492  "ShapeCut",
493  hcaldqm::hashfunctions::fFED,
494  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
495  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),
496  0);
497 
498  _cTimingCut_FEDVME.initialize(_name,
499  "TimingCut",
500  hcaldqm::hashfunctions::fFED,
501  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
502  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
503  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
504  0);
505  _cTimingCut_FEDuTCA.initialize(_name,
506  "TimingCut",
507  hcaldqm::hashfunctions::fFED,
508  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
509  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
510  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
511  0);
512  _cTimingCut_ElectronicsVME.initialize(_name,
513  "TimingCut",
514  hcaldqm::hashfunctions::fElectronics,
515  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
516  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
517  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
518  0);
519  _cTimingCut_ElectronicsuTCA.initialize(_name,
520  "TimingCut",
521  hcaldqm::hashfunctions::fElectronics,
522  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
523  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
524  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
525  0);
526  _cTimingCutvsLS_FED.initialize(_name,
527  "TimingvsLS",
528  hcaldqm::hashfunctions::fFED,
529  new hcaldqm::quantity::LumiSection(_maxLS),
530  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
531  0);
532 
533  _cOccupancy_FEDVME.initialize(_name,
534  "Occupancy",
535  hcaldqm::hashfunctions::fFED,
536  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
537  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
538  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
539  0);
540  _cOccupancy_FEDuTCA.initialize(_name,
541  "Occupancy",
542  hcaldqm::hashfunctions::fFED,
543  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
544  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
545  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
546  0);
547  _cOccupancy_ElectronicsVME.initialize(_name,
548  "Occupancy",
549  hcaldqm::hashfunctions::fElectronics,
550  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
551  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
552  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
553  0);
554  _cOccupancy_ElectronicsuTCA.initialize(_name,
555  "Occupancy",
556  hcaldqm::hashfunctions::fElectronics,
557  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
558  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
559  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
560  0);
561 
562  _cOccupancyCut_FEDVME.initialize(_name,
563  "OccupancyCut",
564  hcaldqm::hashfunctions::fFED,
565  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
566  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
567  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
568  0);
569  _cOccupancyCut_FEDuTCA.initialize(_name,
570  "OccupancyCut",
571  hcaldqm::hashfunctions::fFED,
572  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
573  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
574  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
575  0);
576  _cOccupancyCut_ElectronicsVME.initialize(_name,
577  "OccupancyCut",
578  hcaldqm::hashfunctions::fElectronics,
579  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
580  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
581  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
582  0);
583  _cOccupancyCut_ElectronicsuTCA.initialize(_name,
584  "OccupancyCut",
585  hcaldqm::hashfunctions::fElectronics,
586  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
587  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
588  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
589  0);
590 
591  _cDigiSize_FED.initialize(_name,
592  "DigiSize",
593  hcaldqm::hashfunctions::fFED,
594  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDigiSize),
595  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
596  0);
597 
598  if (_ptype == fOnline) {
599  _cSummaryvsLS_FED.initialize(_name,
600  "SummaryvsLS",
601  hcaldqm::hashfunctions::fFED,
602  new hcaldqm::quantity::LumiSection(_maxLS),
603  new hcaldqm::quantity::FlagQuantity(_vflags),
604  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),
605  0);
606  _cSummaryvsLS.initialize(_name,
607  "SummaryvsLS",
608  new hcaldqm::quantity::LumiSection(_maxLS),
609  new hcaldqm::quantity::FEDQuantity(vFEDs),
610  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),
611  0);
612 
613  _xUniHF.initialize(hcaldqm::hashfunctions::fFEDSlot);
614  _xUni.initialize(hcaldqm::hashfunctions::fFED);
615  _xDigiSize.initialize(hcaldqm::hashfunctions::fFED);
616  _xNChs.initialize(hcaldqm::hashfunctions::fFED);
617  _xNChsNominal.initialize(hcaldqm::hashfunctions::fFED);
618  _xBadCapid.initialize(hcaldqm::hashfunctions::fFED);
619  }
620  }
621  if (_ptype != fLocal) {
622  _LED_ADCvsBX_Subdet.initialize(_name,
623  "LED_ADCvsBX",
624  hcaldqm::hashfunctions::fSubdet,
625  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX_36),
626  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_256_4),
627  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
628  0);
629 
630  _LED_CUCountvsLS_Subdet.initialize(_name,
631  "LED_CUCountvsLS",
632  hcaldqm::hashfunctions::fSubdet,
633  new hcaldqm::quantity::LumiSection(_maxLS),
634  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
635  0);
636  if (_ptype == fOnline) {
637  _LED_CUCountvsLSmod60_Subdet.initialize(_name,
638  "LED_CUCountvsLSmod60",
639  hcaldqm::hashfunctions::fSubdet,
640  new hcaldqm::quantity::LumiSection(60),
641  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
642  0);
643  }
644  }
645 
646  // BOOK HISTOGRAMS
647  char cutstr[200];
648  sprintf(cutstr, "_SumQHBHE%dHO%dHF%d", int(_cutSumQ_HBHE), int(_cutSumQ_HO), int(_cutSumQ_HF));
649  char cutstr2[200];
650  sprintf(cutstr2, "_SumQHF%d", int(_cutSumQ_HF));
651 
652  _cADC_SubdetPM.book(ib, _emap, _filter_QIE8, _subsystem);
653  _cADC_SubdetPM_QIE1011.book(ib, _emap, _filter_QIE1011, _subsystem);
654  _cfC_SubdetPM.book(ib, _emap, _filter_QIE8, _subsystem);
655  _cfC_SubdetPM_QIE1011.book(ib, _emap, _filter_QIE1011, _subsystem);
656  _cSumQ_SubdetPM.book(ib, _emap, _filter_QIE8, _subsystem);
657  _cSumQ_SubdetPM_QIE1011.book(ib, _emap, _filter_QIE1011, _subsystem);
658  _cSumQ_depth.book(ib, _emap, _subsystem);
659  _cSumQvsLS_SubdetPM.book(ib, _emap, _filter_QIE8, _subsystem);
660  _cSumQvsLS_SubdetPM_QIE1011.book(ib, _emap, _filter_QIE1011, _subsystem);
661  _cDigiSize_Crate.book(ib, _emap, _subsystem);
662  _cADCvsTS_SubdetPM.book(ib, _emap, _filter_QIE8, _subsystem);
663  _cADCvsTS_SubdetPM_QIE1011.book(ib, _emap, _filter_QIE1011, _subsystem);
664 
665  if (_ptype != fOffline) { // hidefed2crate
666  _cShapeCut_FED.book(ib, _emap, _subsystem);
667  _cTimingCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
668  _cTimingCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
669  _cTimingCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
670  _cTimingCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
671  _cTimingCutvsLS_FED.book(ib, _emap, _subsystem);
672  _cOccupancy_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
673  _cOccupancy_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
674  _cOccupancy_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
675  _cOccupancy_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
676  _cOccupancyCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
677  _cOccupancyCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
678  _cOccupancyCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
679  _cOccupancyCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
680  _cDigiSize_FED.book(ib, _emap, _subsystem);
681  }
682 
683  _cTimingCut_SubdetPM.book(ib, _emap, _subsystem);
684  _cTimingCut_depth.book(ib, _emap, _subsystem);
685  _cTimingCutvsLS_SubdetPM.book(ib, _emap, _subsystem);
686 
687  _cOccupancyvsLS_Subdet.book(ib, _emap, _subsystem);
688  _cOccupancy_depth.book(ib, _emap, _subsystem);
689  _cOccupancyCut_depth.book(ib, _emap, _subsystem);
690 
691  _cLETDCTimevsADC_SubdetPM.book(ib, _emap, _subsystem);
692  _cLETDCvsADC_2bit_SubdetPM.book(ib, _emap, _filter_TDC2bit, _subsystem);
693  _cLETDCvsADC_6bit_SubdetPM.book(ib, _emap, _filter_TDC6bit, _subsystem);
694  _cLETDCvsTS_2bit_SubdetPM.book(ib, _emap, _filter_TDC2bit, _subsystem);
695  _cLETDCvsTS_6bit_SubdetPM.book(ib, _emap, _filter_TDC6bit, _subsystem);
696  _cLETDCTime_SubdetPM.book(ib, _emap, _subsystem);
697  _cLETDCTime_depth.book(ib, _emap, _subsystem);
698 
699  _cBadTDCValues_SubdetPM.book(ib, _emap, _subsystem);
700  _cBadTDCvsBX_SubdetPM.book(ib, _emap, _subsystem);
701  _cBadTDCvsLS_SubdetPM.book(ib, _emap, _subsystem);
702  _cBadTDCCount_depth.book(ib, _emap, _subsystem);
703 
704  _cCapidMinusBXmod4_SubdetPM.book(ib, _emap, _subsystem);
705  if (_ptype == fOnline) {
706  _cCapid_BadvsFEDvsLS.book(ib, _subsystem, "BadvsLS");
707  _cCapid_BadvsFEDvsLSmod60.book(ib, _subsystem, "BadvsLSmod60");
708  }
709  for (int i = 0; i < 4; ++i) {
710  constexpr unsigned int kSize = 16;
711  char aux[kSize];
712  snprintf(aux, kSize, "%d_uTCA", i);
713  _cCapidMinusBXmod4_CrateSlotuTCA[i].book(ib, _subsystem, aux);
714 
715  snprintf(aux, kSize, "%d_VME", i);
716  _cCapidMinusBXmod4_CrateSlotVME[i].book(ib, _subsystem, aux);
717  }
718 
719  if (_ptype != fLocal) {
720  _LED_ADCvsBX_Subdet.book(ib, _emap, _subsystem);
721  _LED_CUCountvsLS_Subdet.book(ib, _emap, _subsystem);
722  if (_ptype == fOnline) {
723  _LED_CUCountvsLSmod60_Subdet.book(ib, _emap, _subsystem);
724  }
725  }
726 
727  // BOOK HISTOGRAMS that are only for Online
728  _ehashmap.initialize(_emap, electronicsmap::fD2EHashMap);
729  _dhashmap.initialize(_emap, electronicsmap::fE2DHashMap);
730 
731  if (_ptype == fOnline || _ptype == fLocal) {
732  _cOccupancy_Crate.book(ib, _emap, _subsystem);
733  _cOccupancy_CrateSlot.book(ib, _emap, _subsystem);
734  }
735 
736  if (_ptype == fOnline) {
737  _cQ2Q12CutvsLS_FEDHF.book(ib, _emap, _filter_FEDHF, _subsystem);
738  _cSumQvsBX_SubdetPM.book(ib, _emap, _filter_QIE8, _subsystem);
739  _cSumQvsBX_SubdetPM_QIE1011.book(ib, _emap, _filter_QIE1011, _subsystem);
740  _cDigiSizevsLS_FED.book(ib, _emap, _subsystem);
741  _cTimingCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
742  _cTimingCutvsieta_Subdet.book(ib, _emap, _subsystem);
743  _cOccupancyCutvsLS_Subdet.book(ib, _emap, _subsystem);
744  _cOccupancyCutvsBX_Subdet.book(ib, _emap, _subsystem);
745  _cOccupancyvsiphi_SubdetPM.book(ib, _emap, _subsystem);
746  _cOccupancyvsieta_Subdet.book(ib, _emap, _subsystem);
747  _cOccupancyCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
748  _cOccupancyCutvsieta_Subdet.book(ib, _emap, _subsystem);
749  // _cOccupancyCutvsSlotvsLS_HFPM.book(ib, _emap, _filter_QIE1011, _subsystem);
750  _cOccupancyCutvsiphivsLS_SubdetPM.book(ib, _emap, _subsystem);
751  _cSummaryvsLS_FED.book(ib, _emap, _subsystem);
752  _cSummaryvsLS.book(ib, _subsystem);
753 
754  _xUniHF.book(_emap, _filter_FEDHF);
755  _xNChs.book(_emap);
756  _xNChsNominal.book(_emap);
757  _xUni.book(_emap);
758  _xDigiSize.book(_emap);
759  _xBadCapid.book(_emap);
760 
761  // just PER HF FED RECORD THE #CHANNELS
762  // ONLY WAY TO DO THAT AUTOMATICALLY AND W/O HARDCODING 1728
763  // or ANY OTHER VALUES LIKE 2592, 2192
764  std::vector<HcalGenericDetId> gids = _emap->allPrecisionId();
765  for (std::vector<HcalGenericDetId>::const_iterator it = gids.begin(); it != gids.end(); ++it) {
766  if (!it->isHcalDetId())
767  continue;
768  HcalDetId did(it->rawId());
769  if (_xQuality.exists(did)) {
770  HcalChannelStatus cs(it->rawId(), _xQuality.get(HcalDetId(*it)));
771  if (cs.isBitSet(HcalChannelStatus::HcalCellMask) || cs.isBitSet(HcalChannelStatus::HcalCellDead))
772  continue;
773  }
774  HcalElectronicsId eid = HcalElectronicsId(_ehashmap.lookup(did));
775  _xNChsNominal.get(eid)++; // he will know the nominal #channels per FED
776  }
777  }
778 
779  // MARK THESE HISTOGRAMS AS LUMI BASED FOR OFFLINE PROCESSING
780  if (_ptype == fOffline) {
781  _cDigiSize_Crate.setLumiFlag();
782  //_cDigiSize_FED.setLumiFlag();
783  _cOccupancy_depth.setLumiFlag();
784  }
785 
786  // book Number of Events vs LS histogram
787  ib.setCurrentFolder(_subsystem + "/RunInfo");
788  meNumEvents1LS = ib.book1D("NumberOfEvents", "NumberOfEvents", 1, 0, 1);
789  meNumEvents1LS->setLumiFlag();
790 
791  // book the flag for unknown ids and the online guy as well
792  ib.setCurrentFolder(_subsystem + "/" + _name);
793  meUnknownIds1LS = ib.book1D("UnknownIds", "UnknownIds", 1, 0, 1);
794  _unknownIdsPresent = false;
795  meUnknownIds1LS->setLumiFlag();
796 }
797 
798 /* virtual */ void DigiTask::_resetMonitors(hcaldqm::UpdateFreq uf) {
799  DQTask::_resetMonitors(uf);
800 
801  switch (uf) {
802  case hcaldqm::f1LS:
803  _unknownIdsPresent = false;
804  break;
805  case hcaldqm::f50LS:
806  // ^^^ONLINE ONLY!
807  if (_ptype == fOnline)
808  _cOccupancyvsiphi_SubdetPM.reset();
809  // ^^^
810  break;
811  default:
812  break;
813  }
814 }
815 
816 /* virtual */ void DigiTask::_process(edm::Event const& e, edm::EventSetup const&) {
817  edm::Handle<QIE11DigiCollection> c_QIE11;
818  edm::Handle<HODigiCollection> c_ho;
819  edm::Handle<QIE10DigiCollection> c_QIE10;
820 
821  if (!e.getByToken(_tokQIE11, c_QIE11))
822  _logger.dqmthrow("Collection QIE11DigiCollection isn't available" + _tagQIE11.label() + " " + _tagQIE11.instance());
823  if (!e.getByToken(_tokHO, c_ho))
824  _logger.dqmthrow("Collection HODigiCollection isn't available" + _tagHO.label() + " " + _tagHO.instance());
825  if (!e.getByToken(_tokQIE10, c_QIE10))
826  _logger.dqmthrow("Collection QIE10DigiCollection isn't available" + _tagQIE10.label() + " " + _tagQIE10.instance());
827 
828  // extract some info per event
829  int bx = e.bunchCrossing();
830  meNumEvents1LS->Fill(0.5); // just increment
831 
832  // To fill histograms outside of the loop, you need to determine if there were
833  // any valid det ids first
834  uint32_t rawidValid = 0;
835  uint32_t rawidHBValid = 0;
836  uint32_t rawidHEValid = 0;
837 
838  // HB collection
839  int numChs = 0;
840  int numChsCut = 0;
841  int numChsHE = 0;
842  int numChsCutHE = 0;
843 
844  // HB+HE QIE11 collection
845  for (QIE11DigiCollection::const_iterator it=c_QIE11->begin(); it!=c_QIE11->end(); ++it) {
846  const QIE11DataFrame digi = static_cast<const QIE11DataFrame>(*it);
847 
848  // Explicit check on the DetIds present in the Collection
849  HcalDetId const& did = digi.detid();
850  if ((did.subdet() != HcalBarrel) && (did.subdet() != HcalEndcap)) {
851  // LED monitoring from calibration channels
852  if (_ptype != fLocal) {
853  if (did.subdet() == HcalOther) {
854  HcalOtherDetId hodid(digi.detid());
855  if (hodid.subdet() == HcalCalibration) {
856  // New method: use configurable list of channels
857  if (std::find(_ledCalibrationChannels[HcalEndcap].begin(),
858  _ledCalibrationChannels[HcalEndcap].end(),
859  did) != _ledCalibrationChannels[HcalEndcap].end()) {
860  bool channelLEDSignalPresent = false;
861  for (int i = 0; i < digi.samples(); i++) {
862  _LED_ADCvsBX_Subdet.fill(HcalDetId(HcalEndcap, 16, 1, 1), bx, digi[i].adc());
863 
864  if (digi[i].adc() > _thresh_led) {
865  channelLEDSignalPresent = true;
866  }
867  }
868  if (channelLEDSignalPresent) {
869  _LED_CUCountvsLS_Subdet.fill(HcalDetId(HcalEndcap, 16, 1, 1), _currentLS);
870  if (_ptype == fOnline) {
871  _LED_CUCountvsLSmod60_Subdet.fill(HcalDetId(HcalEndcap, 16, 1, 1), _currentLS % 60);
872  }
873  }
874  }
875  }
876  }
877  }
878  continue;
879  }
880 
881  uint32_t rawid = _ehashmap.lookup(did);
882  if (rawid == 0) {
883  meUnknownIds1LS->Fill(1);
884  _unknownIdsPresent = true;
885  continue;
886  } else {
887  if (did.subdet() == HcalBarrel) {
888  rawidHBValid = did.rawId();
889  } else if (did.subdet() == HcalEndcap) {
890  rawidHEValid = did.rawId();
891  }
892  }
893  HcalElectronicsId const& eid(rawid);
894 
895  // filter out channels that are masked out
896  if (_xQuality.exists(did)) {
897  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
898  if (cs.isBitSet(HcalChannelStatus::HcalCellMask) || cs.isBitSet(HcalChannelStatus::HcalCellDead))
899  continue;
900  }
901 
902  // (capid - BX) % 4
903  if (_ptype == fOnline) {
904  short soi = -1;
905  for (int i = 0; i < digi.samples(); i++) {
906  if (digi[i].soi()) {
907  soi = i;
908  break;
909  }
910  }
911  short this_capidmbx = (digi[soi].capid() - bx) % 4;
912  if (this_capidmbx < 0) {
913  this_capidmbx += 4;
914  }
915  _cCapidMinusBXmod4_SubdetPM.fill(did, this_capidmbx);
916  bool good_capidmbx = (_capidmbx[did.subdet()] == this_capidmbx);
917  if (!good_capidmbx) {
918  _xBadCapid.get(eid)++;
919  _cCapid_BadvsFEDvsLS.fill(eid, _currentLS);
920  _cCapid_BadvsFEDvsLSmod60.fill(eid, _currentLS % 60);
921  }
922  if (eid.isVMEid()) {
923  _cCapidMinusBXmod4_CrateSlotVME[this_capidmbx].fill(eid);
924 
925  } else {
926  _cCapidMinusBXmod4_CrateSlotuTCA[this_capidmbx].fill(eid);
927  }
928  }
929 
930  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE11DataFrame>(_dbService, did, digi);
931  double sumQ = hcaldqm::utilities::sumQDB<QIE11DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples() - 1);
932 
933  _cSumQ_SubdetPM_QIE1011.fill(did, sumQ);
934  _cOccupancy_depth.fill(did);
935  if (_ptype == fOnline || _ptype == fLocal) {
936  _cOccupancy_Crate.fill(eid);
937  _cOccupancy_CrateSlot.fill(eid);
938  }
939  if (_ptype == fOnline) {
940  _cDigiSizevsLS_FED.fill(eid, _currentLS, digi.samples());
941  digi.samples() != _refDigiSize[did.subdet()] ? _xDigiSize.get(eid)++ : _xDigiSize.get(eid) += 0;
942  _cOccupancyvsiphi_SubdetPM.fill(did);
943  _cOccupancyvsieta_Subdet.fill(did);
944  }
945  _cDigiSize_Crate.fill(eid, digi.samples());
946  if (_ptype != fOffline) { // hidefed2crate
947  _cDigiSize_FED.fill(eid, digi.samples());
948  if (eid.isVMEid()) {
949  _cOccupancy_FEDVME.fill(eid);
950  _cOccupancy_ElectronicsVME.fill(eid);
951  } else {
952  _cOccupancy_FEDuTCA.fill(eid);
953  _cOccupancy_ElectronicsuTCA.fill(eid);
954  /*
955  if (!digi.validate(0, digi.size()))
956  {
957  _cCapIdRots_depth.fill(did);
958  _cCapIdRots_FEDuTCA.fill(eid, 1);
959  }*/
960  }
961  }
962  for (int i = 0; i < digi.samples(); i++) {
963  double q = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE11DataFrame>(_dbService, digi_fC, did, digi, i);
964  _cADC_SubdetPM_QIE1011.fill(did, digi[i].adc());
965  _cfC_SubdetPM_QIE1011.fill(did, q);
966 
967  if (did.subdet() == HcalBarrel) {
968  _cLETDCvsADC_2bit_SubdetPM.fill(did, digi[i].adc(), digi[i].tdc());
969  _cLETDCvsTS_2bit_SubdetPM.fill(did, (int)i, digi[i].tdc());
970 
971  if (digi[i].tdc() < 2) {
972  double time = i * 25. + (digi[i].tdc() * 12.5);
973  _cLETDCTime_SubdetPM.fill(did, time);
974  _cLETDCTime_depth.fill(did, time);
975  _cLETDCTimevsADC_SubdetPM.fill(did, digi[i].adc(), time);
976  }
977  } else if (did.subdet() == HcalEndcap) {
978  _cLETDCvsADC_6bit_SubdetPM.fill(did, digi[i].adc(), digi[i].tdc());
979  _cLETDCvsTS_6bit_SubdetPM.fill(did, (int)i, digi[i].tdc());
980  if (digi[i].tdc() < 50) {
981  double time = i * 25. + (digi[i].tdc() / 2.);
982  _cLETDCTime_SubdetPM.fill(did, time);
983  _cLETDCTime_depth.fill(did, time);
984  _cLETDCTimevsADC_SubdetPM.fill(did, digi[i].adc(), time);
985  }
986  // Bad TDC values: 50-61 should never happen in QIE10 or QIE11, but we saw some in 2017 data.
987  if ((50 <= digi[i].tdc()) && (digi[i].tdc() <= 61)) {
988  _cBadTDCValues_SubdetPM.fill(did, digi[i].tdc());
989  _cBadTDCvsBX_SubdetPM.fill(did, bx);
990  _cBadTDCvsLS_SubdetPM.fill(did, _currentLS);
991  _cBadTDCCount_depth.fill(did);
992  }
993  }
994  if (_ptype != fOffline) { // hidefed2crate
995  _cADCvsTS_SubdetPM_QIE1011.fill(did, i, digi[i].adc());
996  if (sumQ > _cutSumQ_HBHE) {
997  _cShapeCut_FED.fill(eid, i, q);
998  }
999  }
1000  }
1001 
1002  if (sumQ > _cutSumQ_HBHE) {
1003  //double timing = hcaldqm::utilities::aveTS_v10<QIE11DataFrame>(digi, 2.5, 0,digi.samples()-1);
1004  double timing =
1005  hcaldqm::utilities::aveTSDB<QIE11DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples() - 1);
1006  _cTimingCut_SubdetPM.fill(did, timing);
1007  _cTimingCut_depth.fill(did, timing);
1008  _cOccupancyCut_depth.fill(did);
1009  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
1010  if (_ptype != fOffline) { // hidefed2crate
1011  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
1012  }
1013  _cSumQ_depth.fill(did, sumQ);
1014  _cSumQvsLS_SubdetPM_QIE1011.fill(did, _currentLS, sumQ);
1015  if (_ptype == fOnline) {
1016  _cSumQvsBX_SubdetPM_QIE1011.fill(did, bx, sumQ);
1017  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
1018  _cTimingCutvsieta_Subdet.fill(did, timing);
1019  _cOccupancyCutvsiphi_SubdetPM.fill(did);
1020  _cOccupancyCutvsieta_Subdet.fill(did);
1021  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
1022  }
1023  if (_ptype != fOffline) { // hidefed2crate
1024  if (eid.isVMEid()) {
1025  _cTimingCut_FEDVME.fill(eid, timing);
1026  _cTimingCut_ElectronicsVME.fill(eid, timing);
1027  _cOccupancyCut_FEDVME.fill(eid);
1028  _cOccupancyCut_ElectronicsVME.fill(eid);
1029  } else {
1030  _cTimingCut_FEDuTCA.fill(eid, timing);
1031  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
1032  _cOccupancyCut_FEDuTCA.fill(eid);
1033  _cOccupancyCut_ElectronicsuTCA.fill(eid);
1034  }
1035  }
1036  did.subdet() == HcalBarrel ? numChsCut++ : numChsCutHE++;
1037  }
1038  did.subdet() == HcalBarrel ? numChs++ : numChsHE++;
1039  }
1040 
1041  if (rawidHBValid != 0 && rawidHEValid != 0) {
1042  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHBValid), _currentLS, numChs);
1043  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHEValid), _currentLS, numChsHE);
1044  // ONLINE ONLY!
1045  if (_ptype == fOnline) {
1046  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHBValid), _currentLS, numChsCut);
1047  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHBValid), bx, numChsCut);
1048  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHEValid), _currentLS, numChsCutHE);
1049  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHEValid), bx, numChsCutHE);
1050  }
1051  // ^^^ONLINE ONLY!
1052  }
1053  numChs = 0;
1054  numChsCut = 0;
1055 
1056  // reset
1057  rawidValid = 0;
1058 
1059  // HO collection
1060  for (HODigiCollection::const_iterator it = c_ho->begin(); it != c_ho->end(); ++it) {
1061  // Explicit check on the DetIds present in the Collection
1062  HcalDetId const& did = it->id();
1063  if (did.subdet() != HcalOuter) {
1064  continue;
1065  }
1066  uint32_t rawid = _ehashmap.lookup(did);
1067  if (rawid == 0) {
1068  meUnknownIds1LS->Fill(1);
1069  _unknownIdsPresent = true;
1070  continue;
1071  } else {
1072  rawidValid = did.rawId();
1073  }
1074  HcalElectronicsId const& eid(rawid);
1075 
1076  // filter out channels that are masked out
1077  if (_xQuality.exists(did)) {
1078  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
1079  if (cs.isBitSet(HcalChannelStatus::HcalCellMask) || cs.isBitSet(HcalChannelStatus::HcalCellDead))
1080  continue;
1081  }
1082 
1083  if (_ptype == fOnline) {
1084  short this_capidmbx = (it->sample(it->presamples()).capid() - bx) % 4;
1085  if (this_capidmbx < 0) {
1086  this_capidmbx += 4;
1087  }
1088  _cCapidMinusBXmod4_SubdetPM.fill(did, this_capidmbx);
1089  bool good_capidmbx = (_capidmbx[did.subdet()] == this_capidmbx);
1090  if (!good_capidmbx) {
1091  _xBadCapid.get(eid)++;
1092  _cCapid_BadvsFEDvsLS.fill(eid, _currentLS);
1093  _cCapid_BadvsFEDvsLSmod60.fill(eid, _currentLS % 60);
1094  }
1095  if (eid.isVMEid()) {
1096  _cCapidMinusBXmod4_CrateSlotVME[this_capidmbx].fill(eid);
1097 
1098  } else {
1099  _cCapidMinusBXmod4_CrateSlotuTCA[this_capidmbx].fill(eid);
1100  }
1101  }
1102 
1103  //double sumQ = hcaldqm::utilities::sumQ<HODataFrame>(*it, 8.5, 0, it->size()-1);
1104  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<HODataFrame>(_dbService, did, *it);
1105  double sumQ = hcaldqm::utilities::sumQDB<HODataFrame>(_dbService, digi_fC, did, *it, 0, it->size() - 1);
1106 
1107  _cSumQ_SubdetPM.fill(did, sumQ);
1108  _cOccupancy_depth.fill(did);
1109  if (_ptype == fOnline) {
1110  _cDigiSizevsLS_FED.fill(eid, _currentLS, it->size());
1111  it->size() != _refDigiSize[did.subdet()] ? _xDigiSize.get(eid)++ : _xDigiSize.get(eid) += 0;
1112  _cOccupancyvsiphi_SubdetPM.fill(did);
1113  _cOccupancyvsieta_Subdet.fill(did);
1114  }
1115  _cDigiSize_Crate.fill(eid, it->size());
1116  if (_ptype != fOffline) { // hidefed2crate
1117  _cDigiSize_FED.fill(eid, it->size());
1118  if (eid.isVMEid()) {
1119  _cOccupancy_FEDVME.fill(eid);
1120  _cOccupancy_ElectronicsVME.fill(eid);
1121  /*
1122  if (!it->validate(0, it->size()))
1123  _cCapIdRots_FEDVME.fill(eid, 1);
1124  */
1125  } else {
1126  _cOccupancy_FEDuTCA.fill(eid);
1127  _cOccupancy_ElectronicsuTCA.fill(eid);
1128  /*
1129  if (!it->validate(0, it->size()))
1130  _cCapIdRots_FEDuTCA.fill(eid, 1);*/
1131  }
1132  }
1133 
1134  for (int i = 0; i < it->size(); i++) {
1135  _cADC_SubdetPM.fill(did, it->sample(i).adc());
1136  _cfC_SubdetPM.fill(did, it->sample(i).nominal_fC());
1137  if (_ptype != fOffline) { // hidefed2crate
1138  _cADCvsTS_SubdetPM.fill(did, i, it->sample(i).adc());
1139  if (sumQ > _cutSumQ_HO)
1140  _cShapeCut_FED.fill(eid, i, it->sample(i).nominal_fC());
1141  }
1142  }
1143 
1144  if (sumQ > _cutSumQ_HO) {
1145  //double timing = hcaldqm::utilities::aveTS<HODataFrame>(*it, 8.5, 0,it->size()-1);
1146  double timing = hcaldqm::utilities::aveTSDB<HODataFrame>(_dbService, digi_fC, did, *it, 0, it->size() - 1);
1147  _cSumQ_depth.fill(did, sumQ);
1148  _cSumQvsLS_SubdetPM.fill(did, _currentLS, sumQ);
1149  _cOccupancyCut_depth.fill(did);
1150  _cTimingCut_SubdetPM.fill(did, timing);
1151  _cTimingCut_depth.fill(did, timing);
1152  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
1153  if (_ptype != fOffline) { // hidefed2crate
1154  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
1155  }
1156  if (_ptype == fOnline) {
1157  _cSumQvsBX_SubdetPM.fill(did, bx, sumQ);
1158  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
1159  _cTimingCutvsieta_Subdet.fill(did, timing);
1160  _cOccupancyCutvsiphi_SubdetPM.fill(did);
1161  _cOccupancyCutvsieta_Subdet.fill(did);
1162  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
1163  }
1164  if (_ptype != fOffline) { // hidefed2crate
1165  if (eid.isVMEid()) {
1166  _cTimingCut_FEDVME.fill(eid, timing);
1167  _cTimingCut_ElectronicsVME.fill(eid, timing);
1168  _cOccupancyCut_FEDVME.fill(eid);
1169  _cOccupancyCut_ElectronicsVME.fill(eid);
1170  } else {
1171  _cTimingCut_FEDuTCA.fill(eid, timing);
1172  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
1173  _cOccupancyCut_FEDuTCA.fill(eid);
1174  _cOccupancyCut_ElectronicsuTCA.fill(eid);
1175  }
1176  }
1177  numChsCut++;
1178  }
1179  numChs++;
1180  }
1181 
1182  if (rawidValid != 0) {
1183  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS, numChs);
1184 
1185  if (_ptype == fOnline) {
1186  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS, numChsCut);
1187  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid), bx, numChsCut);
1188  }
1189  }
1190  numChs = 0;
1191  numChsCut = 0;
1192 
1193  // reset
1194  rawidValid = 0;
1195 
1196  // HF collection
1197  if (_qie10InConditions) {
1198  for (QIE10DigiCollection::const_iterator it = c_QIE10->begin(); it != c_QIE10->end(); ++it) {
1199  const QIE10DataFrame digi = static_cast<const QIE10DataFrame>(*it);
1200 
1201  // Explicit check on the DetIds present in the Collection
1202  HcalDetId const& did = digi.detid();
1203  if (did.subdet() != HcalForward) {
1204  // LED monitoring from calibration channels
1205  if (_ptype != fLocal) {
1206  if (did.subdet() == HcalOther) {
1207  HcalOtherDetId hodid(digi.detid());
1208  if (hodid.subdet() == HcalCalibration) {
1209  // New method: use configurable list of channels
1210  if (std::find(_ledCalibrationChannels[HcalForward].begin(),
1211  _ledCalibrationChannels[HcalForward].end(),
1212  did) != _ledCalibrationChannels[HcalForward].end()) {
1213  bool channelLEDSignalPresent = false;
1214  for (int i = 0; i < digi.samples(); i++) {
1215  _LED_ADCvsBX_Subdet.fill(HcalDetId(HcalForward, 16, 1, 1), bx, digi[i].adc());
1216 
1217  if (digi[i].adc() > _thresh_led) {
1218  channelLEDSignalPresent = true;
1219  }
1220  }
1221  if (channelLEDSignalPresent) {
1222  _LED_CUCountvsLS_Subdet.fill(HcalDetId(HcalForward, 16, 1, 1), _currentLS);
1223  if (_ptype == fOnline) {
1224  _LED_CUCountvsLSmod60_Subdet.fill(HcalDetId(HcalForward, 16, 1, 1), _currentLS % 60);
1225  }
1226  }
1227  }
1228  }
1229  }
1230  }
1231  continue;
1232  }
1233 
1234  uint32_t rawid = _ehashmap.lookup(did);
1235  if (rawid == 0) {
1236  meUnknownIds1LS->Fill(1);
1237  _unknownIdsPresent = true;
1238  continue;
1239  } else {
1240  rawidValid = did.rawId();
1241  }
1242  HcalElectronicsId const& eid(rawid);
1243 
1244  // filter out channels that are masked out
1245  if (_xQuality.exists(did)) {
1246  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
1247  if (cs.isBitSet(HcalChannelStatus::HcalCellMask) || cs.isBitSet(HcalChannelStatus::HcalCellDead))
1248  continue;
1249  }
1250 
1251  // (capid - BX) % 4
1252  if (_ptype == fOnline) {
1253  short soi = -1;
1254  for (int i = 0; i < digi.samples(); i++) {
1255  if (digi[i].soi()) {
1256  soi = i;
1257  break;
1258  }
1259  }
1260  short this_capidmbx = (digi[soi].capid() - bx) % 4;
1261  if (this_capidmbx < 0) {
1262  this_capidmbx += 4;
1263  }
1264  _cCapidMinusBXmod4_SubdetPM.fill(did, this_capidmbx);
1265  bool good_capidmbx = (_capidmbx[did.subdet()] == this_capidmbx);
1266  if (!good_capidmbx) {
1267  _xBadCapid.get(eid)++;
1268  _cCapid_BadvsFEDvsLS.fill(eid, _currentLS);
1269  _cCapid_BadvsFEDvsLSmod60.fill(eid, _currentLS % 60);
1270  }
1271  if (eid.isVMEid()) {
1272  _cCapidMinusBXmod4_CrateSlotVME[this_capidmbx].fill(eid);
1273 
1274  } else {
1275  _cCapidMinusBXmod4_CrateSlotuTCA[this_capidmbx].fill(eid);
1276  }
1277  }
1278 
1279  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE10DataFrame>(_dbService, did, digi);
1280  double sumQ = hcaldqm::utilities::sumQDB<QIE10DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples() - 1);
1281  //double sumQ = hcaldqm::utilities::sumQ_v10<QIE10DataFrame>(digi, 2.5, 0, digi.samples()-1);
1282 
1283  //if (!_filter_QIE1011.filter(did)) {
1284  _cSumQ_SubdetPM_QIE1011.fill(did, sumQ);
1285  //}
1286 
1287  _cOccupancy_depth.fill(did);
1288  if (_ptype == fOnline) {
1289  _xNChs.get(eid)++;
1290  _cDigiSizevsLS_FED.fill(eid, _currentLS, digi.samples());
1291  digi.samples() != _refDigiSize[did.subdet()] ? _xDigiSize.get(eid)++ : _xDigiSize.get(eid) += 0;
1292  _cOccupancyvsiphi_SubdetPM.fill(did);
1293  _cOccupancyvsieta_Subdet.fill(did);
1294  }
1295  _cDigiSize_Crate.fill(eid, digi.samples());
1296  if (_ptype != fOffline) { // hidefed2crate
1297  _cDigiSize_FED.fill(eid, digi.samples());
1298  if (eid.isVMEid()) {
1299  _cOccupancy_FEDVME.fill(eid);
1300  _cOccupancy_ElectronicsVME.fill(eid);
1301  /*
1302  if (!it->validate(0, it->size()))
1303  _cCapIdRots_FEDVME.fill(eid, 1);*/
1304  } else {
1305  _cOccupancy_FEDuTCA.fill(eid);
1306  _cOccupancy_ElectronicsuTCA.fill(eid);
1307  /*
1308  if (!it->validate(0, it->size()))
1309  _cCapIdRots_FEDuTCA.fill(eid, 1);*/
1310  }
1311  }
1312 
1313  for (int i = 0; i < digi.samples(); i++) {
1314  double q = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, i);
1315  //if (!_filter_QIE1011.filter(did)) {
1316  _cADC_SubdetPM_QIE1011.fill(did, digi[i].adc());
1317  _cfC_SubdetPM_QIE1011.fill(did, q);
1318  _cLETDCvsADC_6bit_SubdetPM.fill(did, digi[i].adc(), digi[i].le_tdc());
1319  _cLETDCvsTS_6bit_SubdetPM.fill(did, (int)i, digi[i].le_tdc());
1320  if (digi[i].le_tdc() < 50) {
1321  double time = i * 25. + (digi[i].le_tdc() / 2.);
1322  _cLETDCTime_SubdetPM.fill(did, time);
1323  _cLETDCTime_depth.fill(did, time);
1324  _cLETDCTimevsADC_SubdetPM.fill(did, digi[i].adc(), time);
1325  }
1326 
1327  // Bad TDC values: 50-61 should never happen in QIE10 or QIE11, but we are seeing some in 2017 data.
1328  if ((50 <= digi[i].le_tdc()) && (digi[i].le_tdc() <= 61)) {
1329  _cBadTDCValues_SubdetPM.fill(did, digi[i].le_tdc());
1330  _cBadTDCvsBX_SubdetPM.fill(did, bx);
1331  _cBadTDCvsLS_SubdetPM.fill(did, _currentLS);
1332  _cBadTDCCount_depth.fill(did);
1333  }
1334  if (_ptype != fOffline) { // hidefed2crate
1335  _cADCvsTS_SubdetPM_QIE1011.fill(did, (int)i, digi[i].adc());
1336  if (sumQ > _cutSumQ_HF)
1337  _cShapeCut_FED.fill(eid, (int)i, q);
1338  }
1339  //}
1340  }
1341 
1342  if (sumQ > _cutSumQ_HF) {
1343  double timing = hcaldqm::utilities::aveTS_v10<QIE10DataFrame>(digi, 2.5, 0, digi.samples() - 1);
1344  double q1 = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, 1);
1345  double q2 = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, 2);
1346  double q2q12 = q2 / (q1 + q2);
1347  _cSumQ_depth.fill(did, sumQ);
1348  //if (!_filter_QIE1011.filter(did)) {
1349  _cSumQvsLS_SubdetPM_QIE1011.fill(did, _currentLS, sumQ);
1350  //}
1351  _cTimingCut_SubdetPM.fill(did, timing);
1352  _cTimingCut_depth.fill(did, timing);
1353  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
1354  if (_ptype == fOnline) {
1355  //if (!_filter_QIE1011.filter(did)) {
1356  _cSumQvsBX_SubdetPM_QIE1011.fill(did, bx, sumQ);
1357  //}
1358  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
1359  _cTimingCutvsieta_Subdet.fill(did, timing);
1360  _cOccupancyCutvsiphi_SubdetPM.fill(did);
1361  _cOccupancyCutvsieta_Subdet.fill(did);
1362  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
1363  // _cOccupancyCutvsSlotvsLS_HFPM.fill(did, _currentLS);
1364  _xUniHF.get(eid)++;
1365  }
1366  if (_ptype != fOffline) { // hidefed2crate
1367  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
1368  }
1369  _cOccupancyCut_depth.fill(did);
1370  if (!eid.isVMEid())
1371  if (_ptype == fOnline)
1372  _cQ2Q12CutvsLS_FEDHF.fill(eid, _currentLS, q2q12);
1373  if (_ptype != fOffline) { // hidefed2crate
1374  if (eid.isVMEid()) {
1375  _cTimingCut_FEDVME.fill(eid, timing);
1376  _cTimingCut_ElectronicsVME.fill(eid, timing);
1377  _cOccupancyCut_FEDVME.fill(eid);
1378  _cOccupancyCut_ElectronicsVME.fill(eid);
1379  } else {
1380  _cTimingCut_FEDuTCA.fill(eid, timing);
1381  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
1382  _cOccupancyCut_FEDuTCA.fill(eid);
1383  _cOccupancyCut_ElectronicsuTCA.fill(eid);
1384  }
1385  }
1386  numChsCut++;
1387  }
1388  numChs++;
1389  }
1390  }
1391 
1392  if (rawidValid != 0) {
1393  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS, numChs);
1394 
1395  if (_ptype == fOnline) {
1396  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS, numChsCut);
1397  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid), bx, numChsCut);
1398  }
1399  }
1400 }
1401 
1402 /* virtual */ void DigiTask::beginLuminosityBlock(edm::LuminosityBlock const& lb, edm::EventSetup const& es) {
1403  DQTask::beginLuminosityBlock(lb, es);
1404  if (_ptype == fOnline) {
1405  // Reset the bin for _cCapid_BadvsFEDvsLSmod60
1406  for (std::vector<uint32_t>::const_iterator it = _vhashFEDs.begin(); it != _vhashFEDs.end(); ++it) {
1407  HcalElectronicsId eid = HcalElectronicsId(*it);
1408  _cCapid_BadvsFEDvsLSmod60.setBinContent(eid, _currentLS % 50, 0);
1409  }
1410  }
1411 }
1412 
1413 /* virtual */ void DigiTask::endLuminosityBlock(edm::LuminosityBlock const& lb, edm::EventSetup const& es) {
1414  if (_ptype != fOnline)
1415  return;
1416 
1417  for (uintCompactMap::const_iterator it = _xUniHF.begin(); it != _xUniHF.end(); ++it) {
1418  uint32_t hash1 = it->first;
1419  HcalElectronicsId eid1(hash1);
1420  double x1 = it->second;
1421 
1422  for (uintCompactMap::const_iterator jt = _xUniHF.begin(); jt != _xUniHF.end(); ++jt) {
1423  if (jt == it)
1424  continue;
1425  double x2 = jt->second;
1426  if (x2 == 0)
1427  continue;
1428  if (x1 / x2 < _thresh_unihf)
1429  _xUni.get(eid1)++;
1430  }
1431  }
1432 
1433  if (_ptype != fOffline) { // hidefed2crate
1434  for (std::vector<uint32_t>::const_iterator it = _vhashFEDs.begin(); it != _vhashFEDs.end(); ++it) {
1435  hcaldqm::flag::Flag fSum("DIGI");
1436  HcalElectronicsId eid = HcalElectronicsId(*it);
1437 
1438  std::vector<uint32_t>::const_iterator cit = std::find(_vcdaqEids.begin(), _vcdaqEids.end(), *it);
1439  if (cit == _vcdaqEids.end()) {
1440  // not @cDAQ
1441  for (uint32_t iflag = 0; iflag < _vflags.size(); iflag++)
1442  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag), int(hcaldqm::flag::fNCDAQ));
1443  _cSummaryvsLS.setBinContent(eid, _currentLS, int(hcaldqm::flag::fNCDAQ));
1444  continue;
1445  }
1446 
1447  // FED is @cDAQ
1448  if (hcaldqm::utilities::isFEDHBHE(eid) || hcaldqm::utilities::isFEDHF(eid) || hcaldqm::utilities::isFEDHO(eid)) {
1449  if (_xDigiSize.get(eid) > 0)
1450  _vflags[fDigiSize]._state = hcaldqm::flag::fBAD;
1451  else
1452  _vflags[fDigiSize]._state = hcaldqm::flag::fGOOD;
1453 
1454  if (_xBadCapid.get(eid) > 0) {
1455  _vflags[fCapId]._state = hcaldqm::flag::fBAD;
1456  } else {
1457  _vflags[fCapId]._state = hcaldqm::flag::fGOOD;
1458  }
1459 
1460  if (hcaldqm::utilities::isFEDHF(eid)) {
1461  double fr = double(_xNChs.get(eid)) / double(_xNChsNominal.get(eid) * _evsPerLS);
1462  if (_runkeyVal == 0 || _runkeyVal == 4) {
1463  // only for pp or hi
1464  if (_xUni.get(eid) > 0)
1465  _vflags[fUni]._state = hcaldqm::flag::fPROBLEMATIC;
1466  else
1467  _vflags[fUni]._state = hcaldqm::flag::fGOOD;
1468  }
1469  if (fr < 0.95)
1470  _vflags[fNChsHF]._state = hcaldqm::flag::fBAD;
1471  else if (fr < 1.0)
1472  _vflags[fNChsHF]._state = hcaldqm::flag::fPROBLEMATIC;
1473  else
1474  _vflags[fNChsHF]._state = hcaldqm::flag::fGOOD;
1475  }
1476  }
1477  if (_unknownIdsPresent)
1478  _vflags[fUnknownIds]._state = hcaldqm::flag::fBAD;
1479  else
1480  _vflags[fUnknownIds]._state = hcaldqm::flag::fGOOD;
1481 
1482  // LED misfires
1483  if (_ptype != fLocal) {
1484  if (hcaldqm::utilities::isFEDHBHE(eid)) {
1485  HcalDetId did_hb(hcaldqm::hashfunctions::hash_Subdet(HcalDetId(HcalBarrel, 1, 1, 1)));
1486  HcalDetId did_he(hcaldqm::hashfunctions::hash_Subdet(HcalDetId(HcalEndcap, 16, 1, 1)));
1487 
1488  if (_LED_CUCountvsLS_Subdet.getBinContent(did_hb, _currentLS) > 0 ||
1489  _LED_CUCountvsLS_Subdet.getBinContent(did_he, _currentLS) > 0) {
1490  _vflags[fLED]._state = hcaldqm::flag::fBAD;
1491  } else {
1492  _vflags[fLED]._state = hcaldqm::flag::fGOOD;
1493  }
1494  } else if (hcaldqm::utilities::isFEDHF(eid)) {
1495  HcalDetId did_hf(hcaldqm::hashfunctions::hash_Subdet(HcalDetId(HcalForward, 29, 1, 1)));
1496  if (_LED_CUCountvsLS_Subdet.getBinContent(did_hf, _currentLS) > 0) {
1497  _vflags[fLED]._state = hcaldqm::flag::fBAD;
1498  } else {
1499  _vflags[fLED]._state = hcaldqm::flag::fGOOD;
1500  }
1501  } else if (hcaldqm::utilities::isFEDHO(eid)) {
1502  HcalDetId did_ho(hcaldqm::hashfunctions::hash_Subdet(HcalDetId(HcalOuter, 1, 1, 1)));
1503  if (_LED_CUCountvsLS_Subdet.getBinContent(did_ho, _currentLS) > 0) {
1504  _vflags[fLED]._state = hcaldqm::flag::fBAD;
1505  } else {
1506  _vflags[fLED]._state = hcaldqm::flag::fGOOD;
1507  }
1508  }
1509  }
1510 
1511  int iflag = 0;
1512  for (std::vector<hcaldqm::flag::Flag>::iterator ft = _vflags.begin(); ft != _vflags.end(); ++ft) {
1513  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, iflag, int(ft->_state));
1514  fSum += (*ft);
1515  iflag++;
1516 
1517  // reset!
1518  ft->reset();
1519  }
1520  _cSummaryvsLS.setBinContent(eid, _currentLS, fSum._state);
1521  }
1522  }
1523 
1524  _xDigiSize.reset();
1525  _xUniHF.reset();
1526  _xUni.reset();
1527  _xNChs.reset();
1528  _xBadCapid.reset();
1529 
1530  // in the end always do the DQTask::endLumi
1531  DQTask::endLuminosityBlock(lb, es);
1532 }
1533 
1534 DEFINE_FWK_MODULE(DigiTask);
QIE10DataFrame::samples
int samples() const
total number of samples in the digi
Definition: QIE10DataFrame.h:88
hcaldqm::ContainerProf1D::initialize
void initialize(std::string const &folder, hashfunctions::HashType, quantity::Quantity *, quantity::Quantity *, int debug=0) override
Definition: ContainerProf1D.cc:22
QIE11DataFrame::detid
DetId detid() const
Get the detector id.
Definition: QIE11DataFrame.h:60
DigiTask::_cTimingCutvsiphi_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsiphi_SubdetPM
Definition: DigiTask.h:105
HcalElectronicsId::rawId
constexpr uint32_t rawId() const
Definition: HcalElectronicsId.h:65
DigiTask::_process
void _process(edm::Event const &, edm::EventSetup const &) override
Definition: DigiTask.cc:816
DigiTask::_xBadCapid
hcaldqm::ContainerXXX< uint32_t > _xBadCapid
Definition: DigiTask.h:149
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
DigiTask::_qie10InConditions
bool _qie10InConditions
Definition: DigiTask.h:186
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
DigiTask::_thresh_led
double _thresh_led
Definition: DigiTask.h:191
DigiTask::meNumEvents1LS
MonitorElement * meNumEvents1LS
Definition: DigiTask.h:179
DigiTask::_LED_CUCountvsLSmod60_Subdet
hcaldqm::Container1D _LED_CUCountvsLSmod60_Subdet
Definition: DigiTask.h:195
DigiTask::_tagQIE10
edm::InputTag _tagQIE10
Definition: DigiTask.h:44
DigiTask::_cADC_SubdetPM
hcaldqm::Container1D _cADC_SubdetPM
Definition: DigiTask.h:76
DigiTask::_cADC_SubdetPM_QIE1011
hcaldqm::Container1D _cADC_SubdetPM_QIE1011
Definition: DigiTask.h:84
DigiTask::_LED_CUCountvsLS_Subdet
hcaldqm::Container1D _LED_CUCountvsLS_Subdet
Definition: DigiTask.h:194
DigiTask::endLuminosityBlock
void endLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &) override
Definition: DigiTask.cc:1413
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
DigiTask::_cOccupancyCutvsieta_Subdet
hcaldqm::Container1D _cOccupancyCutvsieta_Subdet
Definition: DigiTask.h:132
DigiTask::_cLETDCvsADC_2bit_SubdetPM
hcaldqm::Container2D _cLETDCvsADC_2bit_SubdetPM
Definition: DigiTask.h:152
DigiTask::_cSumQvsBX_SubdetPM
hcaldqm::ContainerProf1D _cSumQvsBX_SubdetPM
Definition: DigiTask.h:81
hcaldqm::ContainerXXX::exists
virtual bool exists(HcalDetId const &)
Definition: ContainerXXX.h:212
mps_fire.i
i
Definition: mps_fire.py:338
hcaldqm::constants::CRATE_VME_MIN
int const CRATE_VME_MIN
Definition: Constants.h:91
edm::DataFrameContainer::const_iterator
boost::transform_iterator< IterHelp, boost::counting_iterator< int > > const_iterator
Definition: DataFrameContainer.h:58
HcalDetId::subdet
HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:146
hcaldqm::ContainerProf2D::fill
void fill(HcalDetId const &) override
Definition: ContainerProf2D.cc:228
DigiTask::_tagHO
edm::InputTag _tagHO
Definition: DigiTask.h:43
DigiTask::_refDigiSize
std::map< HcalSubdetector, int > _refDigiSize
Definition: DigiTask.h:59
DigiTask::_cADCvsTS_SubdetPM
hcaldqm::Container2D _cADCvsTS_SubdetPM
Definition: DigiTask.h:92
hcaldqm::utilities::fed2crate
std::pair< uint16_t, uint16_t > fed2crate(int fed)
Definition: Utilities.cc:10
DigiTask::beginLuminosityBlock
void beginLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &) override
Definition: DigiTask.cc:1402
DigiTask::_cQ2Q12CutvsLS_FEDHF
hcaldqm::ContainerProf1D _cQ2Q12CutvsLS_FEDHF
Definition: DigiTask.h:109
DigiTask::meUnknownIds1LS
MonitorElement * meUnknownIds1LS
Definition: DigiTask.h:180
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
hcaldqm::electronicsmap::ElectronicsMap::initialize
void initialize(HcalElectronicsMap const *, ElectronicsMapType etype=fHcalElectronicsMap)
Definition: ElectronicsMap.cc:6
fwrapper::cs
unique_ptr< ClusterSequence > cs
Definition: fastjetfortran_madfks.cc:45
DigiTask::_cLETDCTimevsADC_SubdetPM
hcaldqm::Container2D _cLETDCTimevsADC_SubdetPM
Definition: DigiTask.h:158
DigiTask::_cutSumQ_HO
double _cutSumQ_HO
Definition: DigiTask.h:49
hcaldqm::ContainerXXX::begin
virtual CompactMap::const_iterator begin()
Definition: ContainerXXX.h:69
DigiTask::_thresh_unihf
double _thresh_unihf
Definition: DigiTask.h:50
hcaldqm::UpdateFreq
UpdateFreq
Definition: DQTask.h:16
hcaldqm::flag::Flag::_state
State _state
Definition: Flag.h:62
DigiTask::_capidmbx
std::map< HcalSubdetector, short > _capidmbx
Definition: DigiTask.h:188
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:50
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:82
DigiTask::_cLETDCvsTS_2bit_SubdetPM
hcaldqm::Container2D _cLETDCvsTS_2bit_SubdetPM
Definition: DigiTask.h:154
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:64
DigiTask::_cCapidMinusBXmod4_CrateSlotuTCA
hcaldqm::ContainerSingle2D _cCapidMinusBXmod4_CrateSlotuTCA[4]
Definition: DigiTask.h:172
DigiTask::_tokQIE10
edm::EDGetTokenT< QIE10DigiCollection > _tokQIE10
Definition: DigiTask.h:47
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:116
hcaldqm::ContainerProf1D::book
void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="") override
Definition: ContainerProf1D.cc:40
DigiTask::_cSumQ_depth
hcaldqm::ContainerProf2D _cSumQ_depth
Definition: DigiTask.h:79
DigiTask::_tokHO
edm::EDGetTokenT< HODigiCollection > _tokHO
Definition: DigiTask.h:46
DigiTask::_cOccupancyCutvsiphi_SubdetPM
hcaldqm::Container1D _cOccupancyCutvsiphi_SubdetPM
Definition: DigiTask.h:131
edm::DataFrameContainer::begin
const_iterator begin() const
Definition: DataFrameContainer.h:159
hcaldqm::ContainerXXX::reset
virtual void reset()
Definition: ContainerXXX.h:292
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
DigiTask::_cSumQvsLS_SubdetPM_QIE1011
hcaldqm::ContainerProf1D _cSumQvsLS_SubdetPM_QIE1011
Definition: DigiTask.h:87
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:20
DigiTask::_cSummaryvsLS
hcaldqm::ContainerSingle2D _cSummaryvsLS
Definition: DigiTask.h:184
DigiTask::_cTimingCut_FEDVME
hcaldqm::ContainerProf2D _cTimingCut_FEDVME
Definition: DigiTask.h:98
DigiTask::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: DigiTask.cc:81
hcaldqm::constants::FIBER_VME_MIN
int const FIBER_VME_MIN
Definition: Constants.h:121
DigiTask::_cOccupancyCutvsiphivsLS_SubdetPM
hcaldqm::Container2D _cOccupancyCutvsiphivsLS_SubdetPM
Definition: DigiTask.h:134
DigiTask::_xDigiSize
hcaldqm::ContainerXXX< uint32_t > _xDigiSize
Definition: DigiTask.h:145
DigiTask::_cTimingCut_FEDuTCA
hcaldqm::ContainerProf2D _cTimingCut_FEDuTCA
Definition: DigiTask.h:99
q2
double q2[4]
Definition: TauolaWrapper.h:88
DigiTask::_cBadTDCvsLS_SubdetPM
hcaldqm::Container1D _cBadTDCvsLS_SubdetPM
Definition: DigiTask.h:163
hcaldqm::constants::FIBERCH_MIN
int const FIBERCH_MIN
Definition: Constants.h:130
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:166
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
DigiTask::_filter_VME
hcaldqm::filter::HashFilter _filter_VME
Definition: DigiTask.h:66
hcaldqm::ContainerXXX::end
virtual CompactMap::const_iterator end()
Definition: ContainerXXX.h:70
hcaldqm::DQModule::_ptype
ProcessingType _ptype
Definition: DQModule.h:43
hcaldqm::Container1D::fill
virtual void fill(uint32_t)
Definition: Container1D.cc:73
hcaldqm::Container1D::reset
virtual void reset()
Definition: Container1D.cc:60
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
DigiTask::_cLETDCvsADC_6bit_SubdetPM
hcaldqm::Container2D _cLETDCvsADC_6bit_SubdetPM
Definition: DigiTask.h:153
hcaldqm::flag::Flag::reset
void reset()
Definition: Flag.h:59
MonitorElement::setLumiFlag
void setLumiFlag()
this ME is meant to be stored for each luminosity section
Definition: MonitorElement.h:150
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:268
DigiTask::_cOccupancy_CrateSlot
hcaldqm::Container2D _cOccupancy_CrateSlot
Definition: DigiTask.h:118
hcaldqm::ContainerSingle2D::setBinContent
virtual void setBinContent(int, int, int)
Definition: ContainerSingle2D.cc:154
HcalElectronicsId::isVMEid
constexpr bool isVMEid() const
Definition: HcalElectronicsId.h:67
Logger::dqmthrow
void dqmthrow(std::string const &msg) const
Definition: Logger.h:12
hcaldqm::DQTask::_xQuality
ContainerXXX< uint32_t > _xQuality
Definition: DQTask.h:48
DigiTask::_cCapid_BadvsFEDvsLS
hcaldqm::ContainerSingle2D _cCapid_BadvsFEDvsLS
Definition: DigiTask.h:174
hcaldqm::constants::FED_VME_MIN
int const FED_VME_MIN
Definition: Constants.h:78
DigiTask::_cADCvsTS_SubdetPM_QIE1011
hcaldqm::Container2D _cADCvsTS_SubdetPM_QIE1011
Definition: DigiTask.h:93
DigiTask::_cOccupancyCut_ElectronicsVME
hcaldqm::Container2D _cOccupancyCut_ElectronicsVME
Definition: DigiTask.h:128
DigiTask::_cDigiSize_Crate
hcaldqm::Container1D _cDigiSize_Crate
Definition: DigiTask.h:142
DigiTask::_cShapeCut_FED
hcaldqm::Container1D _cShapeCut_FED
Definition: DigiTask.h:91
DigiTask::_cOccupancyCut_depth
hcaldqm::Container2D _cOccupancyCut_depth
Definition: DigiTask.h:130
DigiTask::_cOccupancyCutvsBX_Subdet
hcaldqm::ContainerProf1D _cOccupancyCutvsBX_Subdet
Definition: DigiTask.h:139
DigiTask::_cutSumQ_HF
double _cutSumQ_HF
Definition: DigiTask.h:49
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:106
HcalSubdetector
HcalSubdetector
Definition: HcalAssistant.h:31
DigiTask::_tokQIE11
edm::EDGetTokenT< QIE11DigiCollection > _tokQIE11
Definition: DigiTask.h:45
hcaldqm::constants::FIBER_uTCA_MIN1
int const FIBER_uTCA_MIN1
Definition: Constants.h:124
DigiTask::_cfC_SubdetPM_QIE1011
hcaldqm::Container1D _cfC_SubdetPM_QIE1011
Definition: DigiTask.h:85
DigiTask::_cutSumQ_HBHE
double _cutSumQ_HBHE
Definition: DigiTask.h:49
DigiTask::_ledCalibrationChannels
std::map< HcalSubdetector, std::vector< HcalDetId > > _ledCalibrationChannels
Definition: DigiTask.h:192
hcaldqm::hashfunctions::hash_did
hash_function_did const hash_did[nHashType_did]
Definition: HashFunctions.h:163
DigiTask::_cTimingCut_depth
hcaldqm::ContainerProf2D _cTimingCut_depth
Definition: DigiTask.h:104
QIE10DataFrame::detid
DetId detid() const
Get the detector id.
Definition: QIE10DataFrame.h:78
HcalElectronicsMap::allPrecisionId
std::vector< HcalGenericDetId > allPrecisionId() const
Definition: HcalElectronicsMap.cc:139
hcaldqm::DQModule::_logger
Logger _logger
Definition: DQModule.h:54
DigiTask::_cBadTDCValues_SubdetPM
hcaldqm::Container1D _cBadTDCValues_SubdetPM
Definition: DigiTask.h:161
end
#define end
Definition: vmac.h:39
DigiTask::_cOccupancy_FEDVME
hcaldqm::Container2D _cOccupancy_FEDVME
Definition: DigiTask.h:113
ecalMGPA::adc
constexpr int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:11
DigiTask::_cTimingCutvsLS_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsLS_SubdetPM
Definition: DigiTask.h:103
DigiTask::_cOccupancyCut_FEDVME
hcaldqm::Container2D _cOccupancyCut_FEDVME
Definition: DigiTask.h:126
DigiTask::_filter_TDC2bit
hcaldqm::filter::HashFilter _filter_TDC2bit
Definition: DigiTask.h:71
DigiTask::_cSumQvsBX_SubdetPM_QIE1011
hcaldqm::ContainerProf1D _cSumQvsBX_SubdetPM_QIE1011
Definition: DigiTask.h:88
DigiTask::_cCapidMinusBXmod4_CrateSlotVME
hcaldqm::ContainerSingle2D _cCapidMinusBXmod4_CrateSlotVME[4]
Definition: DigiTask.h:173
hcaldqm::hashfunctions::hash_Subdet
uint32_t hash_Subdet(HcalDetId const &)
Definition: HashFunctions.cc:10
DigiTask::_cOccupancy_depth
hcaldqm::Container2D _cOccupancy_depth
Definition: DigiTask.h:119
hcaldqm::DQModule::_name
std::string _name
Definition: DQModule.h:41
hcaldqm::DQTask::_emap
HcalElectronicsMap const * _emap
Definition: DQTask.h:59
globals_cff.x1
x1
Definition: globals_cff.py:29
DigiTask::_cSummaryvsLS_FED
hcaldqm::Container2D _cSummaryvsLS_FED
Definition: DigiTask.h:183
DigiTask::_filter_FEDHF
hcaldqm::filter::HashFilter _filter_FEDHF
Definition: DigiTask.h:68
hcaldqm::constants::CRATE_uTCA_MIN
int const CRATE_uTCA_MIN
Definition: Constants.h:96
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:310
protons_cff.time
time
Definition: protons_cff.py:35
DigiTask::DigiTask
DigiTask(edm::ParameterSet const &)
Definition: DigiTask.cc:7
DigiTask::_cOccupancyvsLS_Subdet
hcaldqm::ContainerProf1D _cOccupancyvsLS_Subdet
Definition: DigiTask.h:137
DigiTask::_resetMonitors
void _resetMonitors(hcaldqm::UpdateFreq) override
Definition: DigiTask.cc:798
DigiTask::_cSumQvsLS_SubdetPM
hcaldqm::ContainerProf1D _cSumQvsLS_SubdetPM
Definition: DigiTask.h:80
DigiTask::_xUni
hcaldqm::ContainerXXX< uint32_t > _xUni
Definition: DigiTask.h:146
DigiTask::_cTimingCutvsLS_FED
hcaldqm::ContainerProf1D _cTimingCutvsLS_FED
Definition: DigiTask.h:102
globals_cff.x2
x2
Definition: globals_cff.py:30
DigiTask::_cOccupancyCut_FEDuTCA
hcaldqm::Container2D _cOccupancyCut_FEDuTCA
Definition: DigiTask.h:127
DigiTask::_cOccupancy_ElectronicsVME
hcaldqm::Container2D _cOccupancy_ElectronicsVME
Definition: DigiTask.h:115
hcaldqm::utilities::isFEDHBHE
bool isFEDHBHE(HcalElectronicsId const &)
Definition: Utilities.cc:142
hcaldqm::ContainerXXX::book
virtual void book(HcalElectronicsMap const *)
Definition: ContainerXXX.h:87
CaloSamples::size
int size() const
get the size
Definition: CaloSamples.h:24
DigiTask::_cDigiSize_FED
hcaldqm::Container1D _cDigiSize_FED
Definition: DigiTask.h:143
DigiTask::_filter_QIE8
hcaldqm::filter::HashFilter _filter_QIE8
Definition: DigiTask.h:70
q1
double q1[4]
Definition: TauolaWrapper.h:87
hcaldqm::Container1D::book
virtual void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="")
Definition: Container1D.cc:592
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
DigiTask::_LED_ADCvsBX_Subdet
hcaldqm::Container2D _LED_ADCvsBX_Subdet
Definition: DigiTask.h:196
DigiTask::_xNChsNominal
hcaldqm::ContainerXXX< uint32_t > _xNChsNominal
Definition: DigiTask.h:148
DigiTask::_cSumQ_SubdetPM
hcaldqm::Container1D _cSumQ_SubdetPM
Definition: DigiTask.h:78
DigiTask::_cOccupancyvsiphi_SubdetPM
hcaldqm::Container1D _cOccupancyvsiphi_SubdetPM
Definition: DigiTask.h:120
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
DigiTask::_cTimingCutvsieta_Subdet
hcaldqm::ContainerProf1D _cTimingCutvsieta_Subdet
Definition: DigiTask.h:106
DigiTask::_cCapid_BadvsFEDvsLSmod60
hcaldqm::ContainerSingle2D _cCapid_BadvsFEDvsLSmod60
Definition: DigiTask.h:176
DigiTask::_cLETDCTime_SubdetPM
hcaldqm::Container1D _cLETDCTime_SubdetPM
Definition: DigiTask.h:156
DigiTask::_dhashmap
hcaldqm::electronicsmap::ElectronicsMap _dhashmap
Definition: DigiTask.h:63
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
DigiTask::_vflags
std::vector< hcaldqm::flag::Flag > _vflags
Definition: DigiTask.h:53
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
edm::DataFrameContainer::end
const_iterator end() const
Definition: DataFrameContainer.h:163
DigiTask::_cOccupancyCut_ElectronicsuTCA
hcaldqm::Container2D _cOccupancyCut_ElectronicsuTCA
Definition: DigiTask.h:129
DigiTask::_cTimingCut_ElectronicsVME
hcaldqm::ContainerProf2D _cTimingCut_ElectronicsVME
Definition: DigiTask.h:100
hcaldqm::ContainerXXX::get
virtual STDTYPE & get(HcalDetId const &)
Definition: ContainerXXX.h:197
begin
#define begin
Definition: vmac.h:32
hcaldqm::DQTask::_dbService
edm::ESHandle< HcalDbService > _dbService
Definition: DQTask.h:58
DigiTask::_cBadTDCvsBX_SubdetPM
hcaldqm::Container1D _cBadTDCvsBX_SubdetPM
Definition: DigiTask.h:162
hcaldqm::Container1D::setLumiFlag
virtual void setLumiFlag()
Definition: Container1D.cc:769
hcaldqm::Container1D::getBinContent
virtual double getBinContent(HcalDetId const &)
Definition: Container1D.cc:177
hcaldqm::utilities::getFEDList
std::vector< int > getFEDList(HcalElectronicsMap const *)
Definition: Utilities.cc:83
DigiTask::_cOccupancy_Crate
hcaldqm::Container2D _cOccupancy_Crate
Definition: DigiTask.h:117
edm::EventSetup::get
T get() const
Definition: EventSetup.h:71
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:309
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
DigiTask::_filter_uTCA
hcaldqm::filter::HashFilter _filter_uTCA
Definition: DigiTask.h:67
hcaldqm::utilities::isFEDHO
bool isFEDHO(HcalElectronicsId const &)
Definition: Utilities.cc:180
hcaldqm::DQTask::_vcdaqEids
std::vector< uint32_t > _vcdaqEids
Definition: DQTask.h:51
DigiTask::_cLETDCvsTS_6bit_SubdetPM
hcaldqm::Container2D _cLETDCvsTS_6bit_SubdetPM
Definition: DigiTask.h:155
DigiTask::_cDigiSizevsLS_FED
hcaldqm::ContainerProf1D _cDigiSizevsLS_FED
Definition: DigiTask.h:144
DigiTask::_xNChs
hcaldqm::ContainerXXX< uint32_t > _xNChs
Definition: DigiTask.h:147
DigiTask::_filter_TDC6bit
hcaldqm::filter::HashFilter _filter_TDC6bit
Definition: DigiTask.h:72
hcaldqm::ContainerSingle2D::fill
virtual void fill(int, int)
Definition: ContainerSingle2D.cc:105
DigiTask::_ehashmap
hcaldqm::electronicsmap::ElectronicsMap _ehashmap
Definition: DigiTask.h:62
DigiTask::_cOccupancyCutvsLS_Subdet
hcaldqm::ContainerProf1D _cOccupancyCutvsLS_Subdet
Definition: DigiTask.h:138
hcaldqm::ContainerXXX::initialize
virtual void initialize(hashfunctions::HashType, int debug=0)
Definition: ContainerXXX.h:81
DigiTask::_filter_QIE1011
hcaldqm::filter::HashFilter _filter_QIE1011
Definition: DigiTask.h:69
hcaldqm::Container2D::book
void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="") override
Definition: Container2D.cc:558
DigiTask::_cBadTDCCount_depth
hcaldqm::Container2D _cBadTDCCount_depth
Definition: DigiTask.h:164
DigiTask::_cLETDCTime_depth
hcaldqm::ContainerProf2D _cLETDCTime_depth
Definition: DigiTask.h:157
DigiTask::_vhashFEDs
std::vector< uint32_t > _vhashFEDs
Definition: DigiTask.h:57
QIE11DataFrame::samples
int samples() const
total number of samples in the digi
Definition: QIE11DataFrame.h:70
DigiTask::_tagQIE11
edm::InputTag _tagQIE11
Definition: DigiTask.h:42
DigiTask::_cOccupancyvsieta_Subdet
hcaldqm::Container1D _cOccupancyvsieta_Subdet
Definition: DigiTask.h:121
HcalElectronicsId
Readout chain identification for Hcal.
Definition: HcalElectronicsId.h:32
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:37
nanoDQM_cff.Flag
Flag
Definition: nanoDQM_cff.py:98
hcaldqm::utilities::isFEDHF
bool isFEDHF(HcalElectronicsId const &)
Definition: Utilities.cc:156
DigiTask::_cTimingCut_ElectronicsuTCA
hcaldqm::ContainerProf2D _cTimingCut_ElectronicsuTCA
Definition: DigiTask.h:101
constexpr
#define constexpr
Definition: GCC11Compatibility.h:35
hcaldqm::constants::SLOT_uTCA_MIN
int const SLOT_uTCA_MIN
Definition: Constants.h:103
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:303
DigiTask::_cSumQ_SubdetPM_QIE1011
hcaldqm::Container1D _cSumQ_SubdetPM_QIE1011
Definition: DigiTask.h:86
DigiTask::_cTimingCut_SubdetPM
hcaldqm::Container1D _cTimingCut_SubdetPM
Definition: DigiTask.h:97
edm::Run
Definition: Run.h:45
DigiTask::_unknownIdsPresent
bool _unknownIdsPresent
Definition: DigiTask.h:181
DigiTask::_cOccupancy_FEDuTCA
hcaldqm::Container2D _cOccupancy_FEDuTCA
Definition: DigiTask.h:114
DigiTask::_cCapidMinusBXmod4_SubdetPM
hcaldqm::Container1D _cCapidMinusBXmod4_SubdetPM
Definition: DigiTask.h:171
cuy.ib
ib
Definition: cuy.py:662
DigiTask::_xUniHF
hcaldqm::ContainerXXX< uint32_t > _xUniHF
Definition: DigiTask.h:146
hcaldqm::utilities::sumQ
double sumQ(DIGI const &digi, double ped, int i=0, int j=3)
Definition: Utilities.h:126
DigiTask::_cfC_SubdetPM
hcaldqm::Container1D _cfC_SubdetPM
Definition: DigiTask.h:77
DigiTask::_cOccupancy_ElectronicsuTCA
hcaldqm::Container2D _cOccupancy_ElectronicsuTCA
Definition: DigiTask.h:116
Generated for CMSSW Reference Manual by doxygen 1.8.11