CMS 3D CMS Logo
LaserTask.cc
Go to the documentation of this file.
1 #include "DQM/HcalTasks/interface/LaserTask.h"
2 
3 using namespace hcaldqm;
4 using namespace hcaldqm::constants;
5 using namespace hcaldqm::filter;
6 LaserTask::LaserTask(edm::ParameterSet const& ps) : DQTask(ps) {
7  _nevents = ps.getUntrackedParameter<int>("nevents", 2000);
8 
9  // tags
10  _tagQIE11 = ps.getUntrackedParameter<edm::InputTag>("tagHE", edm::InputTag("hcalDigis"));
11  _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO", edm::InputTag("hcalDigis"));
12  _tagQIE10 = ps.getUntrackedParameter<edm::InputTag>("tagHF", edm::InputTag("hcalDigis"));
13  _taguMN = ps.getUntrackedParameter<edm::InputTag>("taguMN", edm::InputTag("hcalDigis"));
14  _tagLaserMon = ps.getUntrackedParameter<edm::InputTag>("tagLaserMon", edm::InputTag("LASERMON"));
15  _tokQIE11 = consumes<QIE11DigiCollection>(_tagQIE11);
16  _tokHO = consumes<HODigiCollection>(_tagHO);
17  _tokQIE10 = consumes<QIE10DigiCollection>(_tagQIE10);
18  _tokuMN = consumes<HcalUMNioDigi>(_taguMN);
19  _tokLaserMon = consumes<QIE10DigiCollection>(_tagLaserMon);
20 
21  _vflags.resize(nLaserFlag);
22  _vflags[fBadTiming] = hcaldqm::flag::Flag("BadTiming");
23  _vflags[fMissingLaserMon] = hcaldqm::flag::Flag("MissingLaserMon");
24 
25  // constants
26  _lowHBHE = ps.getUntrackedParameter<double>("lowHBHE", 50);
27  _lowHO = ps.getUntrackedParameter<double>("lowHO", 100);
28  _lowHF = ps.getUntrackedParameter<double>("lowHF", 50);
29  _laserType = (uint32_t)ps.getUntrackedParameter<uint32_t>("laserType");
30 
31  // Laser mon digi ordering list
32  _vLaserMonIPhi = ps.getUntrackedParameter<std::vector<int>>("vLaserMonIPhi");
33  _laserMonIEta = ps.getUntrackedParameter<int>("laserMonIEta");
34  _laserMonCBox = ps.getUntrackedParameter<int>("laserMonCBox");
35  _laserMonDigiOverlap = ps.getUntrackedParameter<int>("laserMonDigiOverlap");
36  _laserMonTS0 = ps.getUntrackedParameter<int>("laserMonTS0");
37  _laserMonThreshold = ps.getUntrackedParameter<double>("laserMonThreshold", 1.e5);
38  _thresh_frac_timingreflm = ps.getUntrackedParameter<double>("_thresh_frac_timingreflm", 0.01);
39  _thresh_min_lmsumq = ps.getUntrackedParameter<double>("thresh_min_lmsumq", 50000.);
40 
41  std::vector<double> vTimingRangeHB =
42  ps.getUntrackedParameter<std::vector<double>>("thresh_timingreflm_HB", std::vector<double>({-70, -10.}));
43  std::vector<double> vTimingRangeHE =
44  ps.getUntrackedParameter<std::vector<double>>("thresh_timingreflm_HE", std::vector<double>({-60., 0.}));
45  std::vector<double> vTimingRangeHO =
46  ps.getUntrackedParameter<std::vector<double>>("thresh_timingreflm_HO", std::vector<double>({-50., 20.}));
47  std::vector<double> vTimingRangeHF =
48  ps.getUntrackedParameter<std::vector<double>>("thresh_timingreflm_HF", std::vector<double>({-50., 20.}));
49  _thresh_timingreflm[HcalBarrel] = std::make_pair(vTimingRangeHB[0], vTimingRangeHB[1]);
50  _thresh_timingreflm[HcalEndcap] = std::make_pair(vTimingRangeHE[0], vTimingRangeHE[1]);
51  _thresh_timingreflm[HcalOuter] = std::make_pair(vTimingRangeHO[0], vTimingRangeHO[1]);
52  _thresh_timingreflm[HcalForward] = std::make_pair(vTimingRangeHF[0], vTimingRangeHF[1]);
53 }
54 
55 /* virtual */ void LaserTask::bookHistograms(DQMStore::IBooker& ib, edm::Run const& r, edm::EventSetup const& es) {
56  if (_ptype == fLocal)
57  if (r.runAuxiliary().run() == 1)
58  return;
59 
60  DQTask::bookHistograms(ib, r, es);
61 
62  edm::ESHandle<HcalDbService> dbService;
63  es.get<HcalDbRecord>().get(dbService);
64  _emap = dbService->getHcalMapping();
65 
66  std::vector<uint32_t> vhashVME;
67  std::vector<uint32_t> vhashuTCA;
68  std::vector<uint32_t> vhashC36;
69  vhashVME.push_back(
70  HcalElectronicsId(constants::FIBERCH_MIN, constants::FIBER_VME_MIN, SPIGOT_MIN, CRATE_VME_MIN).rawId());
71  vhashuTCA.push_back(HcalElectronicsId(CRATE_uTCA_MIN, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
72  _filter_VME.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics, vhashVME);
73  _filter_uTCA.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics, vhashuTCA);
74 
75  // INITIALIZE
76  _cSignalMean_Subdet.initialize(_name,
77  "SignalMean",
78  hcaldqm::hashfunctions::fSubdet,
79  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_100000Coarse),
80  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
81  0);
82  _cSignalRMS_Subdet.initialize(_name,
83  "SignalRMS",
84  hcaldqm::hashfunctions::fSubdet,
85  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_1000),
86  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
87  0);
88  _cTimingMean_Subdet.initialize(_name,
89  "TimingMean",
90  hcaldqm::hashfunctions::fSubdet,
91  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
92  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
93  0);
94  _cTimingRMS_Subdet.initialize(_name,
95  "TimingRMS",
96  hcaldqm::hashfunctions::fSubdet,
97  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
98  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
99  0);
100 
101  _cADC_SubdetPM.initialize(_name,
102  "ADC",
103  hcaldqm::hashfunctions::fSubdetPM,
104  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_128),
105  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
106  0);
107 
108  if (_ptype == fLocal) { // hidefed2crate
109  _cSignalMean_FEDVME.initialize(_name,
110  "SignalMean",
111  hcaldqm::hashfunctions::fFED,
112  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
113  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
114  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_100000Coarse),
115  0);
116  _cSignalMean_FEDuTCA.initialize(_name,
117  "SignalMean",
118  hcaldqm::hashfunctions::fFED,
119  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
120  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
121  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_100000Coarse),
122  0);
123  _cSignalRMS_FEDVME.initialize(_name,
124  "SignalRMS",
125  hcaldqm::hashfunctions::fFED,
126  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
127  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
128  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),
129  0);
130  _cSignalRMS_FEDuTCA.initialize(_name,
131  "SignalRMS",
132  hcaldqm::hashfunctions::fFED,
133  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
134  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
135  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),
136  0);
137  _cTimingMean_FEDVME.initialize(_name,
138  "TimingMean",
139  hcaldqm::hashfunctions::fFED,
140  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
141  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
142  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
143  0);
144  _cTimingMean_FEDuTCA.initialize(_name,
145  "TimingMean",
146  hcaldqm::hashfunctions::fFED,
147  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
148  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
149  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
150  0);
151  _cTimingRMS_FEDVME.initialize(_name,
152  "TimingRMS",
153  hcaldqm::hashfunctions::fFED,
154  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
155  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
156  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
157  0);
158  _cTimingRMS_FEDuTCA.initialize(_name,
159  "TimingRMS",
160  hcaldqm::hashfunctions::fFED,
161  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
162  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
163  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
164  0);
165  _cMissing_FEDVME.initialize(_name,
166  "Missing",
167  hcaldqm::hashfunctions::fFED,
168  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
169  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
170  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
171  0);
172  _cMissing_FEDuTCA.initialize(_name,
173  "Missing",
174  hcaldqm::hashfunctions::fFED,
175  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
176  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
177  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
178  0);
179 
180  _cShapeCut_FEDSlot.initialize(_name,
181  "Shape",
182  hcaldqm::hashfunctions::fFEDSlot,
183  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
184  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),
185  0);
186  }
187  _cTimingvsEvent_SubdetPM.initialize(_name,
188  "TimingvsEvent",
189  hcaldqm::hashfunctions::fSubdetPM,
190  new hcaldqm::quantity::EventNumber(_nevents),
191  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
192  0);
193  _cSignalvsEvent_SubdetPM.initialize(_name,
194  "SignalvsEvent",
195  hcaldqm::hashfunctions::fSubdetPM,
196  new hcaldqm::quantity::EventNumber(_nevents),
197  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),
198  0);
199  _cTimingvsLS_SubdetPM.initialize(_name,
200  "TimingvsLS",
201  hcaldqm::hashfunctions::fSubdetPM,
202  new hcaldqm::quantity::LumiSection(_maxLS),
203  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
204  0);
205  _cSignalvsLS_SubdetPM.initialize(_name,
206  "SignalvsLS",
207  hcaldqm::hashfunctions::fSubdetPM,
208  new hcaldqm::quantity::LumiSection(_maxLS),
209  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),
210  0);
211  _cTimingvsBX_SubdetPM.initialize(_name,
212  "TimingvsBX",
213  hcaldqm::hashfunctions::fSubdetPM,
214  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
215  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
216  0);
217  _cSignalvsBX_SubdetPM.initialize(_name,
218  "SignalvsBX",
219  hcaldqm::hashfunctions::fSubdetPM,
220  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
221  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),
222  0);
223 
224  _cSignalMean_depth.initialize(_name,
225  "SignalMean",
226  hcaldqm::hashfunctions::fdepth,
227  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
228  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
229  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_100000Coarse),
230  0);
231  _cSignalRMS_depth.initialize(_name,
232  "SignalRMS",
233  hcaldqm::hashfunctions::fdepth,
234  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
235  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
236  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_1000),
237  0);
238  _cTimingMean_depth.initialize(_name,
239  "TimingMean",
240  hcaldqm::hashfunctions::fdepth,
241  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
242  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
243  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
244  0);
245  _cTimingRMS_depth.initialize(_name,
246  "TimingRMS",
247  hcaldqm::hashfunctions::fdepth,
248  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
249  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
250  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
251  0);
252 
253  _cMissing_depth.initialize(_name,
254  "Missing",
255  hcaldqm::hashfunctions::fdepth,
256  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
257  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
258  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
259  0);
260 
261  // initialize compact containers
262  _xSignalSum.initialize(hcaldqm::hashfunctions::fDChannel);
263  _xSignalSum2.initialize(hcaldqm::hashfunctions::fDChannel);
264  _xTimingSum.initialize(hcaldqm::hashfunctions::fDChannel);
265  _xTimingSum2.initialize(hcaldqm::hashfunctions::fDChannel);
266  _xEntries.initialize(hcaldqm::hashfunctions::fDChannel);
267 
268  // LaserMon containers
269  if (_ptype == fOnline || _ptype == fLocal) {
270  _cLaserMonSumQ.initialize(_name,
271  "LaserMonSumQ",
272  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_1000000),
273  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
274  0);
275  _cLaserMonSumQ.showOverflowX(true);
276 
277  _cLaserMonTiming.initialize(_name,
278  "LaserMonTiming",
279  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250),
280  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
281  0);
282  _cLaserMonTiming.showOverflowX(true);
283 
284  if (_ptype == fOnline) {
285  _cLaserMonSumQ_LS.initialize(_name,
286  "LaserMonSumQ_LS",
287  new hcaldqm::quantity::LumiSectionCoarse(_maxLS, 10),
288  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_1000000));
289  _cLaserMonTiming_LS.initialize(_name,
290  "LaserMonTiming_LS",
291  new hcaldqm::quantity::LumiSectionCoarse(_maxLS, 10),
292  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_100TS));
293  _cTimingDiffLS_SubdetPM.initialize(_name,
294  "TimingDiff_DigiMinusLaserMon",
295  hcaldqm::hashfunctions::fSubdetPM,
296  new hcaldqm::quantity::LumiSectionCoarse(_maxLS, 10),
297  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fRBX),
298  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTimingDiff_ns),
299  0);
300  _cTimingDiffLS_SubdetPM.showOverflowY(true);
301  } else if (_ptype == fLocal) {
302  _cLaserMonSumQ_Event.initialize(_name,
303  "LaserMonSumQ_Event",
304  new hcaldqm::quantity::EventNumber(_nevents),
305  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_1000000));
306  _cLaserMonTiming_Event.initialize(_name,
307  "LaserMonTiming_Event",
308  new hcaldqm::quantity::EventNumber(_nevents),
309  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_100TS));
310  _cTimingDiffEvent_SubdetPM.initialize(_name,
311  "TimingDiff_DigiMinusLaserMon",
312  hcaldqm::hashfunctions::fSubdetPM,
313  new hcaldqm::quantity::EventNumber(_nevents),
314  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fRBX),
315  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTimingDiff_ns),
316  0);
317  _cTimingDiffEvent_SubdetPM.showOverflowY(true);
318  }
319  _cTiming_DigivsLaserMon_SubdetPM.initialize(
320  _name,
321  "Timing_DigivsLaserMon",
322  hcaldqm::hashfunctions::fSubdetPM,
323  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250_coarse),
324  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250_coarse),
325  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
326  0);
327  _cTiming_DigivsLaserMon_SubdetPM.showOverflowX(true);
328  _cTiming_DigivsLaserMon_SubdetPM.showOverflowY(true);
329 
330  _xTimingRefLMSum.initialize(hcaldqm::hashfunctions::fDChannel);
331  _xTimingRefLMSum2.initialize(hcaldqm::hashfunctions::fDChannel);
332  _xNBadTimingRefLM.initialize(hcaldqm::hashfunctions::fFED);
333  _xNChs.initialize(hcaldqm::hashfunctions::fFED);
334  _xMissingLaserMon = 0;
335 
336  std::vector<int> vFEDs = hcaldqm::utilities::getFEDList(_emap);
337  std::vector<int> vFEDsVME = hcaldqm::utilities::getFEDVMEList(_emap);
338  std::vector<int> vFEDsuTCA = hcaldqm::utilities::getFEDuTCAList(_emap);
339  for (std::vector<int>::const_iterator it = vFEDsVME.begin(); it != vFEDsVME.end(); ++it)
340  _vhashFEDs.push_back(
341  HcalElectronicsId(constants::FIBERCH_MIN, FIBER_VME_MIN, SPIGOT_MIN, (*it) - FED_VME_MIN).rawId());
342  for (std::vector<int>::const_iterator it = vFEDsuTCA.begin(); it != vFEDsuTCA.end(); ++it) {
343  std::pair<uint16_t, uint16_t> cspair = utilities::fed2crate(*it);
344  _vhashFEDs.push_back(HcalElectronicsId(cspair.first, cspair.second, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
345  }
346 
347  _cSummaryvsLS_FED.initialize(_name,
348  "SummaryvsLS",
349  hcaldqm::hashfunctions::fFED,
350  new hcaldqm::quantity::LumiSection(_maxLS),
351  new hcaldqm::quantity::FlagQuantity(_vflags),
352  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),
353  0);
354  _cSummaryvsLS.initialize(_name,
355  "SummaryvsLS",
356  new hcaldqm::quantity::LumiSection(_maxLS),
357  new hcaldqm::quantity::FEDQuantity(vFEDs),
358  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),
359  0);
360  } // End if (_ptype == fOnline || _ptype == fLocal) {
361 
362  // BOOK
363  _cSignalMean_Subdet.book(ib, _emap, _subsystem);
364  _cSignalRMS_Subdet.book(ib, _emap, _subsystem);
365  _cTimingMean_Subdet.book(ib, _emap, _subsystem);
366  _cTimingRMS_Subdet.book(ib, _emap, _subsystem);
367 
368  _cSignalMean_depth.book(ib, _emap, _subsystem);
369  _cSignalRMS_depth.book(ib, _emap, _subsystem);
370  _cTimingMean_depth.book(ib, _emap, _subsystem);
371  _cTimingRMS_depth.book(ib, _emap, _subsystem);
372 
373  if (_ptype == fLocal) {
374  _cTimingvsEvent_SubdetPM.book(ib, _emap, _subsystem);
375  _cSignalvsEvent_SubdetPM.book(ib, _emap, _subsystem);
376  } else {
377  _cTimingvsLS_SubdetPM.book(ib, _emap, _subsystem);
378  _cSignalvsLS_SubdetPM.book(ib, _emap, _subsystem);
379  _cTimingvsBX_SubdetPM.book(ib, _emap, _subsystem);
380  _cSignalvsBX_SubdetPM.book(ib, _emap, _subsystem);
381  }
382 
383  if (_ptype == fLocal) { // hidefed2crate
384  _cSignalMean_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
385  _cSignalMean_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
386  _cSignalRMS_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
387  _cSignalRMS_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
388  _cTimingMean_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
389  _cTimingMean_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
390  _cTimingRMS_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
391  _cTimingRMS_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
392  _cMissing_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
393  _cMissing_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
394  _cShapeCut_FEDSlot.book(ib, _emap, _subsystem);
395  }
396  _cADC_SubdetPM.book(ib, _emap, _subsystem);
397 
398  _cMissing_depth.book(ib, _emap, _subsystem);
399 
400  _xSignalSum.book(_emap);
401  _xSignalSum2.book(_emap);
402  _xEntries.book(_emap);
403  _xTimingSum.book(_emap);
404  _xTimingSum2.book(_emap);
405 
406  if (_ptype == fOnline || _ptype == fLocal) {
407  _cLaserMonSumQ.book(ib, _subsystem);
408  _cLaserMonTiming.book(ib, _subsystem);
409  if (_ptype == fOnline) {
410  _cLaserMonSumQ_LS.book(ib, _subsystem);
411  _cLaserMonTiming_LS.book(ib, _subsystem);
412  _cTimingDiffLS_SubdetPM.book(ib, _emap, _subsystem);
413  } else if (_ptype == fLocal) {
414  _cLaserMonSumQ_Event.book(ib, _subsystem);
415  _cLaserMonTiming_Event.book(ib, _subsystem);
416  _cTimingDiffEvent_SubdetPM.book(ib, _emap, _subsystem);
417  }
418  _cTiming_DigivsLaserMon_SubdetPM.book(ib, _emap, _subsystem);
419  _xTimingRefLMSum.book(_emap);
420  _xTimingRefLMSum2.book(_emap);
421  _xNBadTimingRefLM.book(_emap);
422  _xNChs.book(_emap);
423  _cSummaryvsLS_FED.book(ib, _emap, _subsystem);
424  _cSummaryvsLS.book(ib, _subsystem);
425  }
426 
427  _ehashmap.initialize(_emap, electronicsmap::fD2EHashMap);
428 }
429 
430 /* virtual */ void LaserTask::_resetMonitors(hcaldqm::UpdateFreq uf) { DQTask::_resetMonitors(uf); }
431 
432 /* virtual */ void LaserTask::_dump() {
433  _cSignalMean_Subdet.reset();
434  _cSignalRMS_Subdet.reset();
435  _cTimingMean_Subdet.reset();
436  _cTimingRMS_Subdet.reset();
437  _cSignalMean_depth.reset();
438  _cSignalRMS_depth.reset();
439  _cTimingMean_depth.reset();
440  _cTimingRMS_depth.reset();
441 
442  if (_ptype == fLocal) { // hidefed2crate
443  _cSignalMean_FEDVME.reset();
444  _cSignalMean_FEDuTCA.reset();
445  _cSignalRMS_FEDVME.reset();
446  _cSignalRMS_FEDuTCA.reset();
447  _cTimingMean_FEDVME.reset();
448  _cTimingMean_FEDuTCA.reset();
449  _cTimingRMS_FEDVME.reset();
450  _cTimingRMS_FEDuTCA.reset();
451  }
452  if (_ptype != fOffline) {
453  _xNChs.reset();
454  }
455 
456  std::vector<HcalGenericDetId> dids = _emap->allPrecisionId();
457  for (std::vector<HcalGenericDetId>::const_iterator it = dids.begin(); it != dids.end(); ++it) {
458  if (!it->isHcalDetId())
459  continue;
460  HcalDetId did = HcalDetId(it->rawId());
461  HcalElectronicsId eid(_ehashmap.lookup(*it));
462  int n = _xEntries.get(did);
463  // channels missing or low signal
464  if (n == 0) {
465  _cMissing_depth.fill(did);
466  if (_ptype == fLocal) { // hidefed2crate
467  if (eid.isVMEid())
468  _cMissing_FEDVME.fill(eid);
469  else
470  _cMissing_FEDuTCA.fill(eid);
471  }
472  continue;
473  }
474 
475  _xNChs.get(eid)++;
476 
477  double msig = _xSignalSum.get(did) / n;
478  double mtim = _xTimingSum.get(did) / n;
479  double rsig = sqrt(_xSignalSum2.get(did) / n - msig * msig);
480  double rtim = sqrt(_xTimingSum2.get(did) / n - mtim * mtim);
481 
482  _cSignalMean_Subdet.fill(did, msig);
483  _cSignalMean_depth.fill(did, msig);
484  _cSignalRMS_Subdet.fill(did, rsig);
485  _cSignalRMS_depth.fill(did, rsig);
486  _cTimingMean_Subdet.fill(did, mtim);
487  _cTimingMean_depth.fill(did, mtim);
488  _cTimingRMS_Subdet.fill(did, rtim);
489  _cTimingRMS_depth.fill(did, rtim);
490  if (_ptype == fLocal) { // hidefed2crate
491  if (eid.isVMEid()) {
492  _cSignalMean_FEDVME.fill(eid, msig);
493  _cSignalRMS_FEDVME.fill(eid, rsig);
494  _cTimingMean_FEDVME.fill(eid, mtim);
495  _cTimingRMS_FEDVME.fill(eid, rtim);
496  } else {
497  _cSignalMean_FEDuTCA.fill(eid, msig);
498  _cSignalRMS_FEDuTCA.fill(eid, rsig);
499  _cTimingMean_FEDuTCA.fill(eid, mtim);
500  _cTimingRMS_FEDuTCA.fill(eid, rtim);
501  }
502  }
503 
504  // Bad timing
505 
506  double timingreflm_mean = _xTimingRefLMSum.get(did) / n;
507  //double timingreflm_rms = sqrt(_xTimingRefLMSum2.get(did) / n - timingreflm_mean * timingreflm_mean);
508 
509  if ((timingreflm_mean < _thresh_timingreflm[did.subdet()].first) ||
510  (timingreflm_mean > _thresh_timingreflm[did.subdet()].second)) {
511  _xNBadTimingRefLM.get(eid)++;
512  }
513  }
514  if (_ptype != fOffline) { // hidefed2crate
515  for (std::vector<uint32_t>::const_iterator it = _vhashFEDs.begin(); it != _vhashFEDs.end(); ++it) {
516  hcaldqm::flag::Flag fSum("LASER");
517  HcalElectronicsId eid = HcalElectronicsId(*it);
518  std::vector<uint32_t>::const_iterator jt = std::find(_vcdaqEids.begin(), _vcdaqEids.end(), (*it));
519  if (jt == _vcdaqEids.end()) {
520  // not @cDAQ
521  for (uint32_t iflag = 0; iflag < _vflags.size(); iflag++)
522  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag), int(hcaldqm::flag::fNCDAQ));
523  _cSummaryvsLS.setBinContent(eid, _currentLS, int(hcaldqm::flag::fNCDAQ));
524  continue;
525  }
526  // @cDAQ
527  if (hcaldqm::utilities::isFEDHBHE(eid) || hcaldqm::utilities::isFEDHO(eid) || hcaldqm::utilities::isFEDHF(eid)) {
528  int min_nchs = 10;
529  if (_xNChs.get(eid) < min_nchs) {
530  _vflags[fBadTiming]._state = hcaldqm::flag::fNA;
531  } else {
532  double frbadtimingreflm = double(_xNBadTimingRefLM.get(eid)) / double(_xNChs.get(eid));
533  if (frbadtimingreflm > _thresh_frac_timingreflm) {
534  _vflags[fBadTiming]._state = hcaldqm::flag::fBAD;
535  } else {
536  _vflags[fBadTiming]._state = hcaldqm::flag::fGOOD;
537  }
538  }
539  if (_xMissingLaserMon) {
540  _vflags[fMissingLaserMon]._state = hcaldqm::flag::fBAD;
541  } else {
542  _vflags[fMissingLaserMon]._state = hcaldqm::flag::fGOOD;
543  }
544  }
545 
546  // Set SummaryVsLS bins
547  int iflag = 0;
548  for (std::vector<hcaldqm::flag::Flag>::iterator ft = _vflags.begin(); ft != _vflags.end(); ++ft) {
549  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, iflag, int(ft->_state));
550  fSum += (*ft);
551  iflag++;
552  ft->reset();
553  }
554  _cSummaryvsLS.setBinContent(eid, _currentLS, fSum._state);
555  } // End loop over FEDs
556 
557  // Reset per-LS containers
558  _xNBadTimingRefLM.reset();
559  _xMissingLaserMon = 0;
560  } // End if _ptype != fOffline
561 }
562 
563 /* virtual */ void LaserTask::_process(edm::Event const& e, edm::EventSetup const& es) {
564  edm::Handle<QIE11DigiCollection> c_QIE11;
565  edm::Handle<HODigiCollection> c_ho;
566  edm::Handle<QIE10DigiCollection> c_QIE10;
567 
568  if (!e.getByToken(_tokQIE11, c_QIE11))
569  _logger.dqmthrow("Collection QIE11DigiCollection isn't available " + _tagQIE11.label() + " " + _tagQIE11.instance());
570  if (!e.getByToken(_tokHO, c_ho))
571  _logger.dqmthrow("Collection HODigiCollection isn't available " + _tagHO.label() + " " + _tagHO.instance());
572  if (!e.getByToken(_tokQIE10, c_QIE10))
573  _logger.dqmthrow("Collection QIE10DigiCollection isn't available " + _tagQIE10.label() + " " + _tagQIE10.instance());
574 
575  // int currentEvent = e.eventAuxiliary().id().event();
576  int bx = e.bunchCrossing();
577 
578  // LASERMON
579  edm::Handle<QIE10DigiCollection> cLaserMon;
580  if (!e.getByToken(_tokLaserMon, cLaserMon)) {
581  _logger.dqmthrow("QIE10DigiCollection for laserMonDigis isn't available.");
582  }
583  std::vector<int> laserMonADC;
584  processLaserMon(cLaserMon, laserMonADC);
585 
586  // SumQ = peak +/- 3 TSes
587  // Timing = fC-weighted average (TS-TS0) * 25 ns, also in peak +/- 3 TSes
588  int peakTS = -1;
589  double peakLaserMonADC = -1;
590  for (unsigned int iTS = 0; iTS < laserMonADC.size(); ++iTS) {
591  if (laserMonADC[iTS] > peakLaserMonADC) {
592  peakLaserMonADC = laserMonADC[iTS];
593  peakTS = iTS;
594  }
595  }
596 
597  double laserMonSumQ = 0;
598  double laserMonTiming = 0.;
599 
600  if (peakTS >= 0) {
601  int minTS = std::max(0, peakTS - 3);
602  int maxTS = std::min(int(laserMonADC.size() - 1), peakTS + 3);
603  for (int iTS = minTS; iTS <= maxTS; ++iTS) {
604  double this_fC = hcaldqm::constants::adc2fC[laserMonADC[iTS]];
605  laserMonSumQ += this_fC;
606  laserMonTiming += 25. * (iTS - _laserMonTS0) * this_fC;
607  }
608  }
609  if (laserMonSumQ > 0.) {
610  laserMonTiming = laserMonTiming / laserMonSumQ;
611  } else {
612  ++_xMissingLaserMon;
613  }
614 
615  if (laserMonSumQ > _laserMonThreshold) {
616  _cLaserMonSumQ.fill(laserMonSumQ);
617  _cLaserMonTiming.fill(laserMonTiming);
618  if (_ptype == fOnline) {
619  _cLaserMonSumQ_LS.fill(_currentLS, laserMonSumQ);
620  _cLaserMonTiming_LS.fill(_currentLS, laserMonTiming);
621  } else if (_ptype == fLocal) {
622  int currentEvent = e.eventAuxiliary().id().event();
623  _cLaserMonSumQ_Event.fill(currentEvent, laserMonSumQ);
624  _cLaserMonTiming_Event.fill(currentEvent, laserMonTiming);
625  }
626  }
627 
628  for (QIE11DigiCollection::const_iterator it = c_QIE11->begin(); it != c_QIE11->end(); ++it) {
629  const QIE11DataFrame digi = static_cast<const QIE11DataFrame>(*it);
630  HcalDetId const& did = digi.detid();
631  if ((did.subdet() != HcalBarrel) && (did.subdet() != HcalEndcap)) {
632  continue;
633  }
634  uint32_t rawid = _ehashmap.lookup(did);
635  HcalElectronicsId const& eid(rawid);
636 
637  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE11DataFrame>(_dbService, did, digi);
638  //double sumQ = hcaldqm::utilities::sumQ_v10<QIE11DataFrame>(digi, 2.5, 0, digi.samples()-1);
639  double sumQ = hcaldqm::utilities::sumQDB<QIE11DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples() - 1);
640  if (sumQ < _lowHBHE)
641  continue;
642 
643  //double aveTS = hcaldqm::utilities::aveTS_v10<QIE11DataFrame>(digi, 2.5, 0,digi.samples()-1);
644  double aveTS = hcaldqm::utilities::aveTSDB<QIE11DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples() - 1);
645  _xSignalSum.get(did) += sumQ;
646  _xSignalSum2.get(did) += sumQ * sumQ;
647  _xTimingSum.get(did) += aveTS;
648  _xTimingSum2.get(did) += aveTS * aveTS;
649  _xEntries.get(did)++;
650 
651  for (int i = 0; i < digi.samples(); i++) {
652  if (_ptype == fLocal) {
653  _cShapeCut_FEDSlot.fill(
654  eid, i, hcaldqm::utilities::adc2fCDBMinusPedestal<QIE11DataFrame>(_dbService, digi_fC, did, digi, i));
655  }
656  _cADC_SubdetPM.fill(did, digi[i].adc());
657  }
658 
659  // select based on local global
660  double digiTimingSOI = (aveTS - digi.presamples()) * 25.;
661  double deltaTiming = digiTimingSOI - laserMonTiming;
662  _cTiming_DigivsLaserMon_SubdetPM.fill(did, laserMonTiming, digiTimingSOI);
663  _xTimingRefLMSum.get(did) += deltaTiming;
664  _xTimingRefLMSum2.get(did) += deltaTiming * deltaTiming;
665  if (_ptype == fLocal) {
666  int currentEvent = e.eventAuxiliary().id().event();
667  _cTimingvsEvent_SubdetPM.fill(did, currentEvent, aveTS);
668  _cSignalvsEvent_SubdetPM.fill(did, currentEvent, sumQ);
669  _cTimingDiffEvent_SubdetPM.fill(did, currentEvent, hcaldqm::utilities::getRBX(did.iphi()), deltaTiming);
670  } else {
671  _cTimingvsLS_SubdetPM.fill(did, _currentLS, aveTS);
672  _cSignalvsLS_SubdetPM.fill(did, _currentLS, sumQ);
673  _cTimingvsBX_SubdetPM.fill(did, bx, aveTS);
674  _cSignalvsBX_SubdetPM.fill(did, bx, sumQ);
675  _cTimingDiffLS_SubdetPM.fill(did, _currentLS, hcaldqm::utilities::getRBX(did.iphi()), deltaTiming);
676  }
677  }
678  for (HODigiCollection::const_iterator it = c_ho->begin(); it != c_ho->end(); ++it) {
679  const HODataFrame digi = (const HODataFrame)(*it);
680  double sumQ = hcaldqm::utilities::sumQ<HODataFrame>(digi, 8.5, 0, digi.size() - 1);
681  if (sumQ < _lowHO)
682  continue;
683  HcalDetId did = digi.id();
684  HcalElectronicsId eid = digi.elecId();
685 
686  double aveTS = hcaldqm::utilities::aveTS<HODataFrame>(digi, 8.5, 0, digi.size() - 1);
687  _xSignalSum.get(did) += sumQ;
688  _xSignalSum2.get(did) += sumQ * sumQ;
689  _xTimingSum.get(did) += aveTS;
690  _xTimingSum2.get(did) += aveTS * aveTS;
691  _xEntries.get(did)++;
692 
693  for (int i = 0; i < digi.size(); i++) {
694  if (_ptype == fLocal) { // hidefed2crate
695  _cShapeCut_FEDSlot.fill(eid, i, digi.sample(i).nominal_fC() - 8.5);
696  }
697  _cADC_SubdetPM.fill(did, digi.sample(i).adc());
698  }
699 
700  // select based on local global
701  double digiTimingSOI = (aveTS - digi.presamples()) * 25.;
702  double deltaTiming = digiTimingSOI - laserMonTiming;
703  _cTiming_DigivsLaserMon_SubdetPM.fill(did, laserMonTiming, digiTimingSOI);
704  _xTimingRefLMSum.get(did) += deltaTiming;
705  _xTimingRefLMSum2.get(did) += deltaTiming * deltaTiming;
706  if (_ptype == fLocal) {
707  int currentEvent = e.eventAuxiliary().id().event();
708  _cTimingvsEvent_SubdetPM.fill(did, currentEvent, aveTS);
709  _cSignalvsEvent_SubdetPM.fill(did, currentEvent, sumQ);
710  _cTimingDiffEvent_SubdetPM.fill(did, currentEvent, hcaldqm::utilities::getRBX(did.iphi()), deltaTiming);
711  } else {
712  _cTimingvsLS_SubdetPM.fill(did, _currentLS, aveTS);
713  _cSignalvsLS_SubdetPM.fill(did, _currentLS, sumQ);
714  _cTimingvsBX_SubdetPM.fill(did, bx, aveTS);
715  _cSignalvsBX_SubdetPM.fill(did, bx, sumQ);
716  _cTimingDiffLS_SubdetPM.fill(did, _currentLS, hcaldqm::utilities::getRBX(did.iphi()), deltaTiming);
717  }
718  }
719  for (QIE10DigiCollection::const_iterator it = c_QIE10->begin(); it != c_QIE10->end(); ++it) {
720  const QIE10DataFrame digi = (const QIE10DataFrame)(*it);
721  HcalDetId did = digi.detid();
722  if (did.subdet() != HcalForward) {
723  continue;
724  }
725  HcalElectronicsId eid = HcalElectronicsId(_ehashmap.lookup(did));
726 
727  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE10DataFrame>(_dbService, did, digi);
728  double sumQ = hcaldqm::utilities::sumQDB<QIE10DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples() - 1);
729  //double sumQ = hcaldqm::utilities::sumQ_v10<QIE10DataFrame>(digi, 2.5, 0, digi.samples()-1);
730  if (sumQ < _lowHF)
731  continue;
732 
733  //double aveTS = hcaldqm::utilities::aveTS_v10<QIE10DataFrame>(digi, 2.5, 0, digi.samples()-1);
734  double aveTS = hcaldqm::utilities::aveTSDB<QIE10DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples() - 1);
735 
736  _xSignalSum.get(did) += sumQ;
737  _xSignalSum2.get(did) += sumQ * sumQ;
738  _xTimingSum.get(did) += aveTS;
739  _xTimingSum2.get(did) += aveTS * aveTS;
740  _xEntries.get(did)++;
741 
742  for (int i = 0; i < digi.samples(); i++) {
743  if (_ptype == fLocal) { // hidefed2crate
744  _cShapeCut_FEDSlot.fill(
745  eid, (int)i, hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, i));
746  }
747  _cADC_SubdetPM.fill(did, digi[i].adc());
748  }
749 
750  // select based on local global
751  double digiTimingSOI = (aveTS - digi.presamples()) * 25.;
752  double deltaTiming = digiTimingSOI - laserMonTiming;
753  _cTiming_DigivsLaserMon_SubdetPM.fill(did, laserMonTiming, digiTimingSOI);
754  _xTimingRefLMSum.get(did) += deltaTiming;
755  _xTimingRefLMSum2.get(did) += deltaTiming * deltaTiming;
756  if (_ptype == fLocal) {
757  int currentEvent = e.eventAuxiliary().id().event();
758  _cTimingvsEvent_SubdetPM.fill(did, currentEvent, aveTS);
759  _cSignalvsEvent_SubdetPM.fill(did, currentEvent, sumQ);
760  _cTimingDiffEvent_SubdetPM.fill(did, currentEvent, hcaldqm::utilities::getRBX(did.iphi()), deltaTiming);
761  } else {
762  _cTimingvsLS_SubdetPM.fill(did, _currentLS, aveTS);
763  _cSignalvsLS_SubdetPM.fill(did, _currentLS, sumQ);
764  _cTimingvsBX_SubdetPM.fill(did, bx, aveTS);
765  _cSignalvsBX_SubdetPM.fill(did, bx, sumQ);
766  _cTimingDiffLS_SubdetPM.fill(did, _currentLS, hcaldqm::utilities::getRBX(did.iphi()), deltaTiming);
767  }
768  }
769 }
770 
771 void LaserTask::processLaserMon(edm::Handle<QIE10DigiCollection>& col, std::vector<int>& iLaserMonADC) {
772  for (QIE10DigiCollection::const_iterator it = col->begin(); it != col->end(); ++it) {
773  const QIE10DataFrame digi = (const QIE10DataFrame)(*it);
774  HcalCalibDetId hcdid(digi.id());
775 
776  if ((hcdid.ieta() != _laserMonIEta) || (hcdid.cboxChannel() != _laserMonCBox)) {
777  continue;
778  }
779 
780  unsigned int digiIndex =
781  std::find(_vLaserMonIPhi.begin(), _vLaserMonIPhi.end(), hcdid.iphi()) - _vLaserMonIPhi.begin();
782  if (digiIndex == _vLaserMonIPhi.size()) {
783  continue;
784  }
785 
786  // First digi: initialize the vectors to -1 (need to know the length of the digi)
787  if (iLaserMonADC.empty()) {
788  int totalNSamples = (digi.samples() - _laserMonDigiOverlap) * _vLaserMonIPhi.size();
789  for (int i = 0; i < totalNSamples; ++i) {
790  iLaserMonADC.push_back(-1);
791  }
792  }
793 
794  for (int subindex = 0; subindex < digi.samples() - _laserMonDigiOverlap; ++subindex) {
795  int totalIndex = (digi.samples() - _laserMonDigiOverlap) * digiIndex + subindex;
796  iLaserMonADC[totalIndex] = (digi[subindex].ok() ? digi[subindex].adc() : -1);
797  }
798  }
799 }
800 
801 /* virtual */ void LaserTask::endLuminosityBlock(edm::LuminosityBlock const& lb, edm::EventSetup const& es) {
802  if (_ptype == fLocal)
803  return;
804  this->_dump();
805 
806  DQTask::endLuminosityBlock(lb, es);
807 }
808 
809 /* virtual */ bool LaserTask::_isApplicable(edm::Event const& e) {
810  if (_ptype != fOnline)
811  return true;
812  else {
813  // fOnline mode
814  edm::Handle<HcalUMNioDigi> cumn;
815  if (!e.getByToken(_tokuMN, cumn))
816  return false;
817 
818  // event type check first
819  uint8_t eventType = cumn->eventType();
820  if (eventType != constants::EVENTTYPE_LASER)
821  return false;
822 
823  // check if this analysis task is of the right laser type
824  uint32_t laserType = cumn->valueUserWord(0);
825  if (laserType == _laserType)
826  return true;
827  }
828 
829  return false;
830 }
831 
832 DEFINE_FWK_MODULE(LaserTask);
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
HcalElectronicsId::rawId
constexpr uint32_t rawId() const
Definition: HcalElectronicsId.h:65
LaserTask::_tagQIE11
edm::InputTag _tagQIE11
Definition: LaserTask.h:49
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:41
LaserTask::_cShapeCut_FEDSlot
hcaldqm::ContainerProf1D _cShapeCut_FEDSlot
Definition: LaserTask.h:97
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
hcaldqm::utilities::aveTS
double aveTS(DIGI const &digi, double ped=0, int i=0, int j=3)
Definition: Utilities.h:114
LaserTask::_cSignalvsEvent_SubdetPM
hcaldqm::ContainerProf1D _cSignalvsEvent_SubdetPM
Definition: LaserTask.h:99
LaserTask::_isApplicable
bool _isApplicable(edm::Event const &) override
Definition: LaserTask.cc:809
LaserTask::_cLaserMonTiming
hcaldqm::ContainerSingle1D _cLaserMonTiming
Definition: LaserTask.h:145
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
edm::Event::eventAuxiliary
EventAuxiliary const & eventAuxiliary() const override
Definition: Event.h:92
LaserTask::_vflags
std::vector< hcaldqm::flag::Flag > _vflags
Definition: LaserTask.h:59
mps_fire.i
i
Definition: mps_fire.py:338
LaserTask::_lowHF
double _lowHF
Definition: LaserTask.h:72
hcaldqm::constants::CRATE_VME_MIN
int const CRATE_VME_MIN
Definition: Constants.h:91
hcaldqm::utilities::maxTS
int maxTS(DIGI const &digi, double ped=0)
Definition: Utilities.h:102
LaserTask::_cTimingMean_FEDuTCA
hcaldqm::ContainerProf2D _cTimingMean_FEDuTCA
Definition: LaserTask.h:118
edm::DataFrameContainer::const_iterator
boost::transform_iterator< IterHelp, boost::counting_iterator< int > > const_iterator
Definition: DataFrameContainer.h:58
LaserTask::_xNBadTimingRefLM
hcaldqm::ContainerXXX< int > _xNBadTimingRefLM
Definition: LaserTask.h:83
HcalDetId::subdet
HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:146
hcaldqm::ContainerSingle1D::book
virtual void book(DQMStore::IBooker &, std::string subsystem="Hcal", std::string aux="")
Definition: ContainerSingle1D.cc:57
LaserTask::_xSignalSum2
hcaldqm::ContainerXXX< double > _xSignalSum2
Definition: LaserTask.h:77
hcaldqm::ContainerProf2D::fill
void fill(HcalDetId const &) override
Definition: ContainerProf2D.cc:228
hcaldqm::utilities::fed2crate
std::pair< uint16_t, uint16_t > fed2crate(int fed)
Definition: Utilities.cc:10
LaserTask::_vhashFEDs
std::vector< uint32_t > _vhashFEDs
Definition: LaserTask.h:65
LaserTask::_xTimingRefLMSum
hcaldqm::ContainerXXX< double > _xTimingRefLMSum
Definition: LaserTask.h:81
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
LaserTask::_xEntries
hcaldqm::ContainerXXX< int > _xEntries
Definition: LaserTask.h:78
LaserTask::_cTimingRMS_Subdet
hcaldqm::Container1D _cTimingRMS_Subdet
Definition: LaserTask.h:92
LaserTask::_lowHBHE
double _lowHBHE
Definition: LaserTask.h:69
LaserTask::_laserType
uint32_t _laserType
Definition: LaserTask.h:73
hcaldqm::UpdateFreq
UpdateFreq
Definition: DQTask.h:16
LaserTask::_laserMonThreshold
double _laserMonThreshold
Definition: LaserTask.h:138
LaserTask::_dump
virtual void _dump()
Definition: LaserTask.cc:432
hcaldqm::flag::Flag::_state
State _state
Definition: Flag.h:62
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:82
LaserTask::_cTimingvsEvent_SubdetPM
hcaldqm::ContainerProf1D _cTimingvsEvent_SubdetPM
Definition: LaserTask.h:98
LaserTask::_cSignalRMS_FEDVME
hcaldqm::ContainerProf2D _cSignalRMS_FEDVME
Definition: LaserTask.h:119
LaserTask::_thresh_frac_timingreflm
double _thresh_frac_timingreflm
Definition: LaserTask.h:140
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:64
QIE10DataFrame::id
edm::DataFrame::id_type id() const
Definition: QIE10DataFrame.h:79
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
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
LaserTask::_resetMonitors
void _resetMonitors(hcaldqm::UpdateFreq) override
Definition: LaserTask.cc:430
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:20
hcaldqm::constants::FIBER_VME_MIN
int const FIBER_VME_MIN
Definition: Constants.h:121
LaserTask::_xSignalSum
hcaldqm::ContainerXXX< double > _xSignalSum
Definition: LaserTask.h:76
LaserTask::_xTimingSum
hcaldqm::ContainerXXX< double > _xTimingSum
Definition: LaserTask.h:79
LaserTask::_cMissing_depth
hcaldqm::Container2D _cMissing_depth
Definition: LaserTask.h:125
hcaldqm::ContainerSingleProf1D::initialize
void initialize(std::string const &folder, quantity::Quantity *, quantity::Quantity *qy=new quantity::ValueQuantity(quantity::fN), int debug=0) override
Definition: ContainerSingleProf1D.cc:17
LaserTask::_cMissing_FEDVME
hcaldqm::Container2D _cMissing_FEDVME
Definition: LaserTask.h:126
LaserTask::_cSignalMean_depth
hcaldqm::ContainerProf2D _cSignalMean_depth
Definition: LaserTask.h:108
LaserTask::_cLaserMonTiming_LS
hcaldqm::ContainerSingleProf1D _cLaserMonTiming_LS
Definition: LaserTask.h:147
hcaldqm::constants::FIBERCH_MIN
int const FIBERCH_MIN
Definition: Constants.h:130
LaserTask::_laserMonTS0
int _laserMonTS0
Definition: LaserTask.h:137
LaserTask::_tagHO
edm::InputTag _tagHO
Definition: LaserTask.h:50
LaserTask::_cADC_SubdetPM
hcaldqm::Container1D _cADC_SubdetPM
Definition: LaserTask.h:94
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
LaserTask::_cSummaryvsLS_FED
hcaldqm::Container2D _cSummaryvsLS_FED
Definition: LaserTask.h:155
LaserTask::_cTimingRMS_FEDVME
hcaldqm::ContainerProf2D _cTimingRMS_FEDVME
Definition: LaserTask.h:121
LaserTask::_cTimingDiffLS_SubdetPM
hcaldqm::ContainerProf2D _cTimingDiffLS_SubdetPM
Definition: LaserTask.h:151
LaserTask::_cSignalMean_FEDVME
hcaldqm::ContainerProf2D _cSignalMean_FEDVME
Definition: LaserTask.h:115
LaserTask::_cTimingMean_Subdet
hcaldqm::Container1D _cTimingMean_Subdet
Definition: LaserTask.h:91
LaserTask::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: LaserTask.cc:55
LaserTask::_laserMonDigiOverlap
int _laserMonDigiOverlap
Definition: LaserTask.h:136
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
LaserTask::_cTimingRMS_FEDuTCA
hcaldqm::ContainerProf2D _cTimingRMS_FEDuTCA
Definition: LaserTask.h:122
LaserTask::_cSignalRMS_FEDuTCA
hcaldqm::ContainerProf2D _cSignalRMS_FEDuTCA
Definition: LaserTask.h:120
LaserTask::_tokLaserMon
edm::EDGetTokenT< QIE10DigiCollection > _tokLaserMon
Definition: LaserTask.h:131
LaserTask::_thresh_timingreflm
std::map< HcalSubdetector, std::pair< double, double > > _thresh_timingreflm
Definition: LaserTask.h:139
hcaldqm::ContainerSingle1D::initialize
virtual void initialize(std::string const &folder, quantity::Quantity *, quantity::Quantity *qy=new quantity::ValueQuantity(quantity::fN), int debug=0)
Definition: ContainerSingle1D.cc:37
hcaldqm::flag::Flag::reset
void reset()
Definition: Flag.h:59
LaserTask::processLaserMon
void processLaserMon(edm::Handle< QIE10DigiCollection > &col, std::vector< int > &iLaserMonADC)
Definition: LaserTask.cc:771
hcaldqm::ContainerSingle2D::setBinContent
virtual void setBinContent(int, int, int)
Definition: ContainerSingle2D.cc:154
LaserTask::_cTimingDiffEvent_SubdetPM
hcaldqm::ContainerProf2D _cTimingDiffEvent_SubdetPM
Definition: LaserTask.h:152
LaserTask::_cLaserMonSumQ_LS
hcaldqm::ContainerSingleProf1D _cLaserMonSumQ_LS
Definition: LaserTask.h:146
Logger::dqmthrow
void dqmthrow(std::string const &msg) const
Definition: Logger.h:12
HcalUMNioDigi::valueUserWord
uint32_t valueUserWord(int iword) const
Definition: HcalUMNioDigi.cc:54
LaserTask::_filter_uTCA
hcaldqm::filter::HashFilter _filter_uTCA
Definition: LaserTask.h:63
LaserTask::_xNChs
hcaldqm::ContainerXXX< int > _xNChs
Definition: LaserTask.h:84
hcaldqm::constants::FED_VME_MIN
int const FED_VME_MIN
Definition: Constants.h:78
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
LaserTask::_tokQIE11
edm::EDGetTokenT< QIE11DigiCollection > _tokQIE11
Definition: LaserTask.h:53
LaserTask::_cLaserMonTiming_Event
hcaldqm::ContainerSingleProf1D _cLaserMonTiming_Event
Definition: LaserTask.h:149
LaserTask::_cTimingMean_FEDVME
hcaldqm::ContainerProf2D _cTimingMean_FEDVME
Definition: LaserTask.h:117
LaserTask::_nevents
int _nevents
Definition: LaserTask.h:68
hcaldqm::Container1D::showOverflowX
virtual void showOverflowX(bool showOverflow)
Definition: Container1D.cc:775
HODataFrame::size
int size() const
total number of samples in the digi
Definition: HODataFrame.h:31
LaserTask::_cLaserMonSumQ
hcaldqm::ContainerSingle1D _cLaserMonSumQ
Definition: LaserTask.h:144
HcalUMNioDigi::eventType
uint8_t eventType() const
Definition: HcalUMNioDigi.cc:32
hcaldqm::constants::FIBER_uTCA_MIN1
int const FIBER_uTCA_MIN1
Definition: Constants.h:124
LaserTask::_tokHO
edm::EDGetTokenT< HODigiCollection > _tokHO
Definition: LaserTask.h:54
QIE10DataFrame::detid
DetId detid() const
Get the detector id.
Definition: QIE10DataFrame.h:78
HcalElectronicsMap::allPrecisionId
std::vector< HcalGenericDetId > allPrecisionId() const
Definition: HcalElectronicsMap.cc:139
LaserTask::endLuminosityBlock
void endLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &) override
Definition: LaserTask.cc:801
hcaldqm::DQModule::_logger
Logger _logger
Definition: DQModule.h:54
LaserTask::_thresh_min_lmsumq
double _thresh_min_lmsumq
Definition: LaserTask.h:141
hcaldqm::Container1D::showOverflowY
virtual void showOverflowY(bool showOverflow)
Definition: Container1D.cc:777
min
T min(T a, T b)
Definition: MathUtil.h:58
ecalMGPA::adc
constexpr int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:11
edm::Run::runAuxiliary
RunAuxiliary const & runAuxiliary() const override
Definition: Run.h:61
LaserTask::_xTimingRefLMSum2
hcaldqm::ContainerXXX< double > _xTimingRefLMSum2
Definition: LaserTask.h:82
hcaldqm::ContainerSingle1D::showOverflowX
virtual void showOverflowX(bool showOverflow)
Definition: ContainerSingle1D.cc:203
hcaldqm::DQModule::_name
std::string _name
Definition: DQModule.h:41
LaserTask::_laserMonCBox
int _laserMonCBox
Definition: LaserTask.h:135
hcaldqm::DQTask::_emap
HcalElectronicsMap const * _emap
Definition: DQTask.h:59
LaserTask::_cSignalvsBX_SubdetPM
hcaldqm::ContainerProf1D _cSignalvsBX_SubdetPM
Definition: LaserTask.h:104
LaserTask::_cSummaryvsLS
hcaldqm::ContainerSingle2D _cSummaryvsLS
Definition: LaserTask.h:156
LaserTask::_tokQIE10
edm::EDGetTokenT< QIE10DigiCollection > _tokQIE10
Definition: LaserTask.h:55
LaserTask::_lowHO
double _lowHO
Definition: LaserTask.h:71
hcaldqm::constants::CRATE_uTCA_MIN
int const CRATE_uTCA_MIN
Definition: Constants.h:96
LaserTask::_cMissing_FEDuTCA
hcaldqm::Container2D _cMissing_FEDuTCA
Definition: LaserTask.h:127
LaserTask::_process
void _process(edm::Event const &, edm::EventSetup const &) override
Definition: LaserTask.cc:563
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:310
LaserTask::_cSignalMean_Subdet
hcaldqm::Container1D _cSignalMean_Subdet
Definition: LaserTask.h:87
hcaldqm::constants::adc2fC
double const adc2fC[256]
Definition: Constants.h:248
HcalDetId::iphi
int iphi() const
get the cell iphi
Definition: HcalDetId.h:161
LaserTask::_cSignalvsLS_SubdetPM
hcaldqm::ContainerProf1D _cSignalvsLS_SubdetPM
Definition: LaserTask.h:101
LaserTask::_cTimingvsLS_SubdetPM
hcaldqm::ContainerProf1D _cTimingvsLS_SubdetPM
Definition: LaserTask.h:100
LaserTask::_laserMonIEta
int _laserMonIEta
Definition: LaserTask.h:134
hcaldqm::utilities::isFEDHBHE
bool isFEDHBHE(HcalElectronicsId const &)
Definition: Utilities.cc:142
hcaldqm::ContainerXXX::book
virtual void book(HcalElectronicsMap const *)
Definition: ContainerXXX.h:87
LaserTask::_cSignalMean_FEDuTCA
hcaldqm::ContainerProf2D _cSignalMean_FEDuTCA
Definition: LaserTask.h:116
LaserTask::_cTiming_DigivsLaserMon_SubdetPM
hcaldqm::Container2D _cTiming_DigivsLaserMon_SubdetPM
Definition: LaserTask.h:150
hcaldqm::Container1D::book
virtual void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="")
Definition: Container1D.cc:592
hcaldqm::constants::EVENTTYPE_LASER
uint8_t const EVENTTYPE_LASER
Definition: Constants.h:288
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
LaserTask::_cTimingMean_depth
hcaldqm::ContainerProf2D _cTimingMean_depth
Definition: LaserTask.h:112
hcaldqm::ContainerSingleProf1D::book
void book(DQMStore::IBooker &, std::string subsystem="Hcal", std::string aux="") override
Definition: ContainerSingleProf1D.cc:33
LaserTask::_ehashmap
hcaldqm::electronicsmap::ElectronicsMap _ehashmap
Definition: LaserTask.h:62
LaserTask::_tagLaserMon
edm::InputTag _tagLaserMon
Definition: LaserTask.h:130
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
LaserTask::_taguMN
edm::InputTag _taguMN
Definition: LaserTask.h:52
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
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
LaserTask::_vLaserMonIPhi
std::vector< int > _vLaserMonIPhi
Definition: LaserTask.h:133
edm::EventAuxiliary::id
EventID const & id() const
Definition: EventAuxiliary.h:59
LaserTask::_xMissingLaserMon
int _xMissingLaserMon
Definition: LaserTask.h:142
hcaldqm::ContainerXXX::get
virtual STDTYPE & get(HcalDetId const &)
Definition: ContainerXXX.h:197
hcaldqm::DQTask::_dbService
edm::ESHandle< HcalDbService > _dbService
Definition: DQTask.h:58
hcaldqm::utilities::getFEDList
std::vector< int > getFEDList(HcalElectronicsMap const *)
Definition: Utilities.cc:83
edm::EventSetup::get
T get() const
Definition: EventSetup.h:71
LaserTask::_tokuMN
edm::EDGetTokenT< HcalUMNioDigi > _tokuMN
Definition: LaserTask.h:56
hcaldqm::utilities::getFEDuTCAList
std::vector< int > getFEDuTCAList(HcalElectronicsMap const *)
Definition: Utilities.cc:121
LaserTask::_cSignalRMS_Subdet
hcaldqm::Container1D _cSignalRMS_Subdet
Definition: LaserTask.h:88
cuy.col
col
Definition: cuy.py:1010
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
hcaldqm::utilities::getRBX
int getRBX(uint32_t iphi)
Definition: Utilities.cc:233
hcaldqm::utilities::isFEDHO
bool isFEDHO(HcalElectronicsId const &)
Definition: Utilities.cc:180
hcaldqm::DQTask::_vcdaqEids
std::vector< uint32_t > _vcdaqEids
Definition: DQTask.h:51
LaserTask::_cLaserMonSumQ_Event
hcaldqm::ContainerSingleProf1D _cLaserMonSumQ_Event
Definition: LaserTask.h:148
LaserTask::_xTimingSum2
hcaldqm::ContainerXXX< double > _xTimingSum2
Definition: LaserTask.h:80
hcaldqm::ContainerXXX::initialize
virtual void initialize(hashfunctions::HashType, int debug=0)
Definition: ContainerXXX.h:81
hcaldqm::Container2D::book
void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="") override
Definition: Container2D.cc:558
QIE11DataFrame::samples
int samples() const
total number of samples in the digi
Definition: QIE11DataFrame.h:70
LaserTask::_cSignalRMS_depth
hcaldqm::ContainerProf2D _cSignalRMS_depth
Definition: LaserTask.h:109
LaserTask::_filter_VME
hcaldqm::filter::HashFilter _filter_VME
Definition: LaserTask.h:64
HcalElectronicsId
Readout chain identification for Hcal.
Definition: HcalElectronicsId.h:32
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:37
LaserTask::_cTimingvsBX_SubdetPM
hcaldqm::ContainerProf1D _cTimingvsBX_SubdetPM
Definition: LaserTask.h:103
nanoDQM_cff.Flag
Flag
Definition: nanoDQM_cff.py:98
hcaldqm::utilities::isFEDHF
bool isFEDHF(HcalElectronicsId const &)
Definition: Utilities.cc:156
edm::RunAuxiliary::run
RunNumber_t run() const
Definition: RunAuxiliary.h:41
LaserTask::LaserTask
LaserTask(edm::ParameterSet const &)
Definition: LaserTask.cc:6
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
edm::Run
Definition: Run.h:45
cuy.ib
ib
Definition: cuy.py:662
LaserTask::_cTimingRMS_depth
hcaldqm::ContainerProf2D _cTimingRMS_depth
Definition: LaserTask.h:113
LaserTask::_tagQIE10
edm::InputTag _tagQIE10
Definition: LaserTask.h:51
hcaldqm::utilities::sumQ
double sumQ(DIGI const &digi, double ped, int i=0, int j=3)
Definition: Utilities.h:126
Generated for CMSSW Reference Manual by doxygen 1.8.11