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