test
CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
LaserTask.cc
Go to the documentation of this file.
1 
2 #include "DQM/HcalTasks/interface/LaserTask.h"
3 
4 using namespace hcaldqm;
5 using namespace hcaldqm::constants;
6 using namespace hcaldqm::filter;
7 LaserTask::LaserTask(edm::ParameterSet const& ps):
8  DQTask(ps)
9 {
10  _nevents = ps.getUntrackedParameter<int>("nevents", 2000);
11 
12  // tags
13  _tagHBHE = ps.getUntrackedParameter<edm::InputTag>("tagHBHE",
14  edm::InputTag("hcalDigis"));
15  _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO",
16  edm::InputTag("hcalDigis"));
17  _tagHF = ps.getUntrackedParameter<edm::InputTag>("tagHF",
18  edm::InputTag("hcalDigis"));
19  _taguMN = ps.getUntrackedParameter<edm::InputTag>("taguMN",
20  edm::InputTag("hcalDigis"));
21  _tokHBHE = consumes<HBHEDigiCollection>(_tagHBHE);
22  _tokHO = consumes<HODigiCollection>(_tagHO);
23  _tokHF = consumes<HFDigiCollection>(_tagHF);
24  _tokuMN = consumes<HcalUMNioDigi>(_taguMN);
25 
26  // constants
27  _lowHBHE = ps.getUntrackedParameter<double>("lowHBHE",
28  20);
29  _lowHO = ps.getUntrackedParameter<double>("lowHO",
30  20);
31  _lowHF = ps.getUntrackedParameter<double>("lowHF",
32  20);
33  _laserType = (uint32_t)ps.getUntrackedParameter<uint32_t>("laserType");
34 }
35 
36 /* virtual */ void LaserTask::bookHistograms(DQMStore::IBooker &ib,
37  edm::Run const& r, edm::EventSetup const& es)
38 {
39  if (_ptype==fLocal)
40  if (r.runAuxiliary().run()==1)
41  return;
42 
43  DQTask::bookHistograms(ib, r, es);
44 
45  edm::ESHandle<HcalDbService> dbService;
46  es.get<HcalDbRecord>().get(dbService);
47  _emap = dbService->getHcalMapping();
48 
49  std::vector<uint32_t> vhashVME;
50  std::vector<uint32_t> vhashuTCA;
51  std::vector<uint32_t> vhashC36;
52  vhashVME.push_back(HcalElectronicsId(constants::FIBERCH_MIN,
53  constants::FIBER_VME_MIN, SPIGOT_MIN, CRATE_VME_MIN).rawId());
54  vhashuTCA.push_back(HcalElectronicsId(CRATE_uTCA_MIN, SLOT_uTCA_MIN,
55  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
56  _filter_VME.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics,
57  vhashVME);
58  _filter_uTCA.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics,
59  vhashuTCA);
60 
61  // INITIALIZE
62  _cSignalMean_Subdet.initialize(_name, "SignalMean",
63  hcaldqm::hashfunctions::fSubdet,
64  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),
65  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
66  _cSignalRMS_Subdet.initialize(_name, "SignalRMS",
67  hcaldqm::hashfunctions::fSubdet,
68  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_1000),
69  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
70  _cTimingMean_Subdet.initialize(_name, "TimingMean",
71  hcaldqm::hashfunctions::fSubdet,
72  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
73  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
74  _cTimingRMS_Subdet.initialize(_name, "TimingRMS",
75  hcaldqm::hashfunctions::fSubdet,
76  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
77  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
78 
79  _cADC_SubdetPM.initialize(_name, "ADC",
80  hcaldqm::hashfunctions::fSubdetPM,
81  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_128),
82  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
83 
84  _cSignalMean_FEDVME.initialize(_name, "SignalMean",
85  hcaldqm::hashfunctions::fFED,
86  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
87  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
88  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000));
89  _cSignalMean_FEDuTCA.initialize(_name, "SignalMean",
90  hcaldqm::hashfunctions::fFED,
91  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
92  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
93  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000));
94  _cSignalRMS_FEDVME.initialize(_name, "SignalRMS",
95  hcaldqm::hashfunctions::fFED,
96  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
97  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
98  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000));
99  _cSignalRMS_FEDuTCA.initialize(_name, "SignalRMS",
100  hcaldqm::hashfunctions::fFED,
101  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
102  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
103  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000));
104  _cTimingMean_FEDVME.initialize(_name, "TimingMean",
105  hcaldqm::hashfunctions::fFED,
106  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
107  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
108  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
109  _cTimingMean_FEDuTCA.initialize(_name, "TimingMean",
110  hcaldqm::hashfunctions::fFED,
111  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
112  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
113  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
114  _cTimingRMS_FEDVME.initialize(_name, "TimingRMS",
115  hcaldqm::hashfunctions::fFED,
116  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
117  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
118  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
119  _cTimingRMS_FEDuTCA.initialize(_name, "TimingRMS",
120  hcaldqm::hashfunctions::fFED,
121  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
122  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
123  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
124 
125  _cShapeCut_FEDSlot.initialize(_name, "Shape",
126  hcaldqm::hashfunctions::fFEDSlot,
127  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
128  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000));
129  _cTimingvsEvent_SubdetPM.initialize(_name, "TimingvsEvent",
130  hcaldqm::hashfunctions::fSubdetPM,
131  new hcaldqm::quantity::EventNumber(_nevents),
132  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
133  _cSignalvsEvent_SubdetPM.initialize(_name, "SignalvsEvent",
134  hcaldqm::hashfunctions::fSubdetPM,
135  new hcaldqm::quantity::EventNumber(_nevents),
136  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000));
137  _cTimingvsLS_SubdetPM.initialize(_name, "TimingvsLS",
138  hcaldqm::hashfunctions::fSubdetPM,
139  new hcaldqm::quantity::LumiSection(_maxLS),
140  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
141  _cSignalvsLS_SubdetPM.initialize(_name, "SignalvsLS",
142  hcaldqm::hashfunctions::fSubdetPM,
143  new hcaldqm::quantity::LumiSection(_maxLS),
144  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000));
145  _cTimingvsBX_SubdetPM.initialize(_name, "TimingvsBX",
146  hcaldqm::hashfunctions::fSubdetPM,
147  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
148  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
149  _cSignalvsBX_SubdetPM.initialize(_name, "SignalvsBX",
150  hcaldqm::hashfunctions::fSubdetPM,
151  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
152  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000));
153 
154  _cSignalMean_depth.initialize(_name, "SignalMean",
155  hcaldqm::hashfunctions::fdepth,
156  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
157  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
158  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000));
159  _cSignalRMS_depth.initialize(_name, "SignalRMS",
160  hcaldqm::hashfunctions::fdepth,
161  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
162  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
163  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_1000));
164  _cTimingMean_depth.initialize(_name, "TimingMean",
165  hcaldqm::hashfunctions::fdepth,
166  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
167  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
168  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
169  _cTimingRMS_depth.initialize(_name, "TimingRMS",
170  hcaldqm::hashfunctions::fdepth,
171  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
172  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
173  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
174 
175  _cMissing_depth.initialize(_name, "Missing",
176  hcaldqm::hashfunctions::fdepth,
177  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
178  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
179  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
180  _cMissing_FEDVME.initialize(_name, "Missing",
181  hcaldqm::hashfunctions::fFED,
182  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
183  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
184  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
185  _cMissing_FEDuTCA.initialize(_name, "Missing",
186  hcaldqm::hashfunctions::fFED,
187  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
188  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
189  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
190 
191 
192  // initialize compact containers
193  _xSignalSum.initialize(hcaldqm::hashfunctions::fDChannel);
194  _xSignalSum2.initialize(hcaldqm::hashfunctions::fDChannel);
195  _xTimingSum.initialize(hcaldqm::hashfunctions::fDChannel);
196  _xTimingSum2.initialize(hcaldqm::hashfunctions::fDChannel);
197  _xEntries.initialize(hcaldqm::hashfunctions::fDChannel);
198 
199  // BOOK
200  _cSignalMean_Subdet.book(ib, _emap, _subsystem);
201  _cSignalRMS_Subdet.book(ib, _emap, _subsystem);
202  _cTimingMean_Subdet.book(ib, _emap, _subsystem);
203  _cTimingRMS_Subdet.book(ib, _emap, _subsystem);
204 
205  _cSignalMean_depth.book(ib, _emap, _subsystem);
206  _cSignalRMS_depth.book(ib, _emap, _subsystem);
207  _cTimingMean_depth.book(ib, _emap, _subsystem);
208  _cTimingRMS_depth.book(ib, _emap, _subsystem);
209 
210  if (_ptype==fLocal)
211  {
212  _cTimingvsEvent_SubdetPM.book(ib, _emap, _subsystem);
213  _cSignalvsEvent_SubdetPM.book(ib, _emap, _subsystem);
214  }
215  else
216  {
217  _cTimingvsLS_SubdetPM.book(ib, _emap, _subsystem);
218  _cSignalvsLS_SubdetPM.book(ib, _emap, _subsystem);
219  _cTimingvsBX_SubdetPM.book(ib, _emap, _subsystem);
220  _cSignalvsBX_SubdetPM.book(ib, _emap, _subsystem);
221  }
222 
223  _cSignalMean_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
224  _cSignalMean_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
225  _cSignalRMS_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
226  _cSignalRMS_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
227  _cTimingMean_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
228  _cTimingMean_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
229  _cTimingRMS_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
230  _cTimingRMS_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
231  _cADC_SubdetPM.book(ib, _emap, _subsystem);
232 
233  _cShapeCut_FEDSlot.book(ib, _emap, _subsystem);
234  _cMissing_depth.book(ib, _emap,_subsystem);
235  _cMissing_FEDVME.book(ib, _emap, _filter_uTCA,_subsystem);
236  _cMissing_FEDuTCA.book(ib, _emap, _filter_VME,_subsystem);
237 
238  _xSignalSum.book(_emap);
239  _xSignalSum2.book(_emap);
240  _xEntries.book(_emap);
241  _xTimingSum.book(_emap);
242  _xTimingSum2.book(_emap);
243 
244  _ehashmap.initialize(_emap, electronicsmap::fD2EHashMap);
245 }
246 
247 /* virtual */ void LaserTask::_resetMonitors(hcaldqm::UpdateFreq uf)
248 {
249  DQTask::_resetMonitors(uf);
250 }
251 
252 /* virtual */ void LaserTask::_dump()
253 {
254  _cSignalMean_Subdet.reset();
255  _cSignalRMS_Subdet.reset();
256  _cTimingMean_Subdet.reset();
257  _cTimingRMS_Subdet.reset();
258  _cSignalMean_depth.reset();
259  _cSignalRMS_depth.reset();
260  _cTimingMean_depth.reset();
261  _cTimingRMS_depth.reset();
262 
263  _cSignalMean_FEDVME.reset();
264  _cSignalMean_FEDuTCA.reset();
265  _cSignalRMS_FEDVME.reset();
266  _cSignalRMS_FEDuTCA.reset();
267  _cTimingMean_FEDVME.reset();
268  _cTimingMean_FEDuTCA.reset();
269  _cTimingRMS_FEDVME.reset();
270  _cTimingRMS_FEDuTCA.reset();
271 
272  std::vector<HcalGenericDetId> dids = _emap->allPrecisionId();
273  for (std::vector<HcalGenericDetId>::const_iterator it=dids.begin();
274  it!=dids.end(); ++it)
275  {
276  if (!it->isHcalDetId())
277  continue;
278  HcalDetId did = HcalDetId(it->rawId());
279  HcalElectronicsId eid(_ehashmap.lookup(*it));
280  int n = _xEntries.get(did);
281  double msig = _xSignalSum.get(did)/n;
282  double mtim = _xTimingSum.get(did)/n;
283  double rsig = sqrt(_xSignalSum2.get(did)/n-msig*msig);
284  double rtim = sqrt(_xTimingSum2.get(did)/n-mtim*mtim);
285 
286  // channels missing or low signal
287  if (n==0)
288  {
289  _cMissing_depth.fill(did);
290  if (eid.isVMEid())
291  _cMissing_FEDVME.fill(eid);
292  else
293  _cMissing_FEDuTCA.fill(eid);
294  continue;
295  }
296  _cSignalMean_Subdet.fill(did, msig);
297  _cSignalMean_depth.fill(did, msig);
298  _cSignalRMS_Subdet.fill(did, rsig);
299  _cSignalRMS_depth.fill(did, rsig);
300  _cTimingMean_Subdet.fill(did, mtim);
301  _cTimingMean_depth.fill(did, mtim);
302  _cTimingRMS_Subdet.fill(did, rtim);
303  _cTimingRMS_depth.fill(did, rtim);
304  if (eid.isVMEid())
305  {
306  _cSignalMean_FEDVME.fill(eid, msig);
307  _cSignalRMS_FEDVME.fill(eid, rsig);
308  _cTimingMean_FEDVME.fill(eid, mtim);
309  _cTimingRMS_FEDVME.fill(eid, rtim);
310  }
311  else
312  {
313  _cSignalMean_FEDuTCA.fill(eid, msig);
314  _cSignalRMS_FEDuTCA.fill(eid, rsig);
315  _cTimingMean_FEDuTCA.fill(eid, mtim);
316  _cTimingRMS_FEDuTCA.fill(eid, rtim);
317  }
318  }
319 }
320 
321 /* virtual */ void LaserTask::_process(edm::Event const& e,
322  edm::EventSetup const& es)
323 {
324  edm::Handle<HBHEDigiCollection> chbhe;
325  edm::Handle<HODigiCollection> cho;
326  edm::Handle<HFDigiCollection> chf;
327 
328  if (!e.getByToken(_tokHBHE, chbhe))
329  _logger.dqmthrow("Collection HBHEDigiCollection isn't available "
330  + _tagHBHE.label() + " " + _tagHBHE.instance());
331  if (!e.getByToken(_tokHO, cho))
332  _logger.dqmthrow("Collection HODigiCollection isn't available "
333  + _tagHO.label() + " " + _tagHO.instance());
334  if (!e.getByToken(_tokHF, chf))
335  _logger.dqmthrow("Collection HFDigiCollection isn't available "
336  + _tagHF.label() + " " + _tagHF.instance());
337 
338 // int currentEvent = e.eventAuxiliary().id().event();
339  int bx = e.bunchCrossing();
340 
341  for (HBHEDigiCollection::const_iterator it=chbhe->begin();
342  it!=chbhe->end(); ++it)
343  {
344  const HBHEDataFrame digi = (const HBHEDataFrame)(*it);
345  double sumQ = hcaldqm::utilities::sumQ<HBHEDataFrame>(digi, 2.5, 0,
346  digi.size()-1);
347  if (sumQ<_lowHBHE)
348  continue;
349  HcalDetId did = digi.id();
350  HcalElectronicsId eid = digi.elecId();
351 
352  double aveTS = hcaldqm::utilities::aveTS<HBHEDataFrame>(digi, 2.5, 0,
353  digi.size()-1);
354  _xSignalSum.get(did)+=sumQ;
355  _xSignalSum2.get(did)+=sumQ*sumQ;
356  _xTimingSum.get(did)+=aveTS;
357  _xTimingSum2.get(did)+=aveTS*aveTS;
358  _xEntries.get(did)++;
359 
360  for (int i=0; i<digi.size(); i++)
361  {
362  _cShapeCut_FEDSlot.fill(eid, i,
363  digi.sample(i).nominal_fC()-2.5);
364  _cADC_SubdetPM.fill(did, digi.sample(i).adc());
365  }
366 
367  // select based on local global
368  if (_ptype==fLocal)
369  {
370  int currentEvent = e.eventAuxiliary().id().event();
371  _cTimingvsEvent_SubdetPM.fill(did, currentEvent, aveTS);
372  _cSignalvsEvent_SubdetPM.fill(did, currentEvent, sumQ);
373  }
374  else
375  {
376  _cTimingvsLS_SubdetPM.fill(did, _currentLS, aveTS);
377  _cSignalvsLS_SubdetPM.fill(did, _currentLS, sumQ);
378  _cTimingvsBX_SubdetPM.fill(did, bx, aveTS);
379  _cSignalvsBX_SubdetPM.fill(did, bx, sumQ);
380  }
381  }
382  for (HODigiCollection::const_iterator it=cho->begin();
383  it!=cho->end(); ++it)
384  {
385  const HODataFrame digi = (const HODataFrame)(*it);
386  double sumQ = hcaldqm::utilities::sumQ<HODataFrame>(digi, 8.5, 0,
387  digi.size()-1);
388  if (sumQ<_lowHO)
389  continue;
390  HcalDetId did = digi.id();
391  HcalElectronicsId eid = digi.elecId();
392 
393  double aveTS = hcaldqm::utilities::aveTS<HODataFrame>(digi, 8.5, 0,
394  digi.size()-1);
395  _xSignalSum.get(did)+=sumQ;
396  _xSignalSum2.get(did)+=sumQ*sumQ;
397  _xTimingSum.get(did)+=aveTS;
398  _xTimingSum2.get(did)+=aveTS*aveTS;
399  _xEntries.get(did)++;
400 
401  for (int i=0; i<digi.size(); i++)
402  {
403  _cShapeCut_FEDSlot.fill(eid, i,
404  digi.sample(i).nominal_fC()-8.5);
405  _cADC_SubdetPM.fill(did, digi.sample(i).adc());
406  }
407 
408  // select based on local global
409  if (_ptype==fLocal)
410  {
411  int currentEvent = e.eventAuxiliary().id().event();
412  _cTimingvsEvent_SubdetPM.fill(did, currentEvent, aveTS);
413  _cSignalvsEvent_SubdetPM.fill(did, currentEvent, sumQ);
414  }
415  else
416  {
417  _cTimingvsLS_SubdetPM.fill(did, _currentLS, aveTS);
418  _cSignalvsLS_SubdetPM.fill(did, _currentLS, sumQ);
419  _cTimingvsBX_SubdetPM.fill(did, bx, aveTS);
420  _cSignalvsBX_SubdetPM.fill(did, bx, sumQ);
421  }
422  }
423  for (HFDigiCollection::const_iterator it=chf->begin();
424  it!=chf->end(); ++it)
425  {
426  const HFDataFrame digi = (const HFDataFrame)(*it);
427  double sumQ = hcaldqm::utilities::sumQ<HFDataFrame>(digi, 2.5, 0,
428  digi.size()-1);
429  if (sumQ<_lowHF)
430  continue;
431  HcalDetId did = digi.id();
432  HcalElectronicsId eid = digi.elecId();
433 
434  double aveTS = hcaldqm::utilities::aveTS<HFDataFrame>(digi, 2.5, 0,
435  digi.size()-1);
436  _xSignalSum.get(did)+=sumQ;
437  _xSignalSum2.get(did)+=sumQ*sumQ;
438  _xTimingSum.get(did)+=aveTS;
439  _xTimingSum2.get(did)+=aveTS*aveTS;
440  _xEntries.get(did)++;
441 
442  for (int i=0; i<digi.size(); i++)
443  {
444  _cShapeCut_FEDSlot.fill(eid, i,
445  digi.sample(i).nominal_fC()-2.5);
446  _cADC_SubdetPM.fill(did, digi.sample(i).adc());
447  }
448 
449  // select based on local global
450  if (_ptype==fLocal)
451  {
452  int currentEvent = e.eventAuxiliary().id().event();
453  _cTimingvsEvent_SubdetPM.fill(did, currentEvent, aveTS);
454  _cSignalvsEvent_SubdetPM.fill(did, currentEvent, sumQ);
455  }
456  else
457  {
458  _cTimingvsLS_SubdetPM.fill(did, _currentLS, aveTS);
459  _cSignalvsLS_SubdetPM.fill(did, _currentLS, sumQ);
460  _cTimingvsBX_SubdetPM.fill(did, bx, aveTS);
461  _cSignalvsBX_SubdetPM.fill(did, bx, sumQ);
462  }
463  }
464 }
465 
466 /* virtual */ void LaserTask::endLuminosityBlock(edm::LuminosityBlock const& lb,
467  edm::EventSetup const& es)
468 {
469  if (_ptype==fLocal)
470  return;
471  this->_dump();
472 
473  DQTask::endLuminosityBlock(lb, es);
474 }
475 
476 /* virtual */ bool LaserTask::_isApplicable(edm::Event const& e)
477 {
478  if (_ptype!=fOnline)
479  return true;
480  else
481  {
482  // fOnline mode
483  edm::Handle<HcalUMNioDigi> cumn;
484  if (!e.getByToken(_tokuMN, cumn))
485  return false;
486 
487  // event type check first
488  uint8_t eventType = cumn->eventType();
489  if (eventType!=constants::EVENTTYPE_LASER)
490  return false;
491 
492  // check if this analysis task is of the right laser type
493  uint32_t laserType = cumn->valueUserWord(0);
494  if (laserType==_laserType) return true;
495  }
496 
497  return false;
498 }
499 
500 DEFINE_FWK_MODULE(LaserTask);
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:41
LaserTask::_cShapeCut_FEDSlot
hcaldqm::ContainerProf1D _cShapeCut_FEDSlot
Definition: LaserTask.h:89
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
hcaldqm::utilities::aveTS
double aveTS(DIGI const &digi, double ped=0, int i=0, int j=3)
Definition: Utilities.h:65
LaserTask::_cSignalvsEvent_SubdetPM
hcaldqm::ContainerProf1D _cSignalvsEvent_SubdetPM
Definition: LaserTask.h:91
hcaldqm::Container1D::initialize
virtual void initialize(std::string const &folder, hashfunctions::HashType, Quantity *, Quantity *qy=new ValueQuantity(quantity::fN), int debug=0)
Definition: Container1D.cc:40
LaserTask::_lowHF
double _lowHF
Definition: LaserTask.h:70
hcaldqm::constants::CRATE_VME_MIN
int const CRATE_VME_MIN
Definition: Constants.h:60
cuy.ib
int ib
Definition: cuy.py:660
LaserTask::_cTimingMean_FEDuTCA
hcaldqm::ContainerProf2D _cTimingMean_FEDuTCA
Definition: LaserTask.h:106
LaserTask::_xSignalSum2
hcaldqm::ContainerXXX< double > _xSignalSum2
Definition: LaserTask.h:75
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:462
hcaldqm::electronicsmap::ElectronicsMap::initialize
void initialize(HcalElectronicsMap const *, ElectronicsMapType etype=fHcalElectronicsMap)
Definition: ElectronicsMap.cc:8
LaserTask::_xEntries
hcaldqm::ContainerXXX< int > _xEntries
Definition: LaserTask.h:76
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
LaserTask::_cTimingRMS_Subdet
hcaldqm::Container1D _cTimingRMS_Subdet
Definition: LaserTask.h:84
HBHEDataFrame::size
int size() const
total number of samples in the digi
Definition: HBHEDataFrame.h:26
LaserTask::_lowHBHE
double _lowHBHE
Definition: LaserTask.h:68
LaserTask::_laserType
uint32_t _laserType
Definition: LaserTask.h:71
hcaldqm::UpdateFreq
UpdateFreq
Definition: DQTask.h:17
LaserTask::_dump
virtual void _dump()
Definition: LaserTask.cc:252
edm::SortedCollection< HBHEDataFrame >::const_iterator
std::vector< HBHEDataFrame >::const_iterator const_iterator
Definition: SortedCollection.h:82
LaserTask::_cTimingvsEvent_SubdetPM
hcaldqm::ContainerProf1D _cTimingvsEvent_SubdetPM
Definition: LaserTask.h:90
LaserTask::_cSignalRMS_FEDVME
hcaldqm::ContainerProf2D _cSignalRMS_FEDVME
Definition: LaserTask.h:107
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:65
stage2BMTFBufferRaw_cfi.eventType
tuple eventType
Definition: stage2BMTFBufferRaw_cfi.py:61
hcaldqm::constants::SPIGOT_MIN
int const SPIGOT_MIN
Definition: Constants.h:85
alcazmumu_cfi.filter
tuple filter
Definition: alcazmumu_cfi.py:34
hcaldqm::constants::FIBER_VME_MIN
int const FIBER_VME_MIN
Definition: Constants.h:90
LaserTask::_xSignalSum
hcaldqm::ContainerXXX< double > _xSignalSum
Definition: LaserTask.h:74
LaserTask::bookHistograms
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
Definition: LaserTask.cc:36
LaserTask::_xTimingSum
hcaldqm::ContainerXXX< double > _xTimingSum
Definition: LaserTask.h:77
LaserTask::_cMissing_depth
hcaldqm::Container2D _cMissing_depth
Definition: LaserTask.h:113
LaserTask::_cMissing_FEDVME
hcaldqm::Container2D _cMissing_FEDVME
Definition: LaserTask.h:114
bookHistograms
void bookHistograms(fwlite::EventContainer &eventCont)
Definition: secvtxMassTemplates.cc:293
LaserTask::_cSignalMean_depth
hcaldqm::ContainerProf2D _cSignalMean_depth
Definition: LaserTask.h:98
hcaldqm::constants::FIBERCH_MIN
int const FIBERCH_MIN
Definition: Constants.h:100
LaserTask::_tagHO
edm::InputTag _tagHO
Definition: LaserTask.h:52
LaserTask::_cADC_SubdetPM
hcaldqm::Container1D _cADC_SubdetPM
Definition: LaserTask.h:86
LaserTask::_emap
HcalElectronicsMap const * _emap
Definition: LaserTask.h:61
HLT_25ns10e33_v2_cff.InputTag
tuple InputTag
Definition: HLT_25ns10e33_v2_cff.py:58467
LaserTask::_cTimingRMS_FEDVME
hcaldqm::ContainerProf2D _cTimingRMS_FEDVME
Definition: LaserTask.h:109
hcaldqm::Container2D::fill
virtual void fill(HcalDetId const &)
Definition: Container2D.cc:59
LaserTask::_cSignalMean_FEDVME
hcaldqm::ContainerProf2D _cSignalMean_FEDVME
Definition: LaserTask.h:103
LaserTask::_tagHF
edm::InputTag _tagHF
Definition: LaserTask.h:53
LaserTask::_cTimingMean_Subdet
hcaldqm::Container1D _cTimingMean_Subdet
Definition: LaserTask.h:83
hcaldqm::DQModule::_ptype
ProcessingType _ptype
Definition: DQModule.h:59
hcaldqm::Container1D::fill
virtual void fill(uint32_t)
Definition: Container1D.cc:82
hcaldqm::Container1D::reset
virtual void reset()
Definition: Container1D.cc:65
hcaldqm::ContainerProf1D::initialize
virtual void initialize(std::string const &folder, hashfunctions::HashType, Quantity *, Quantity *, int debug=0)
Definition: ContainerProf1D.cc:26
LaserTask::_cTimingRMS_FEDuTCA
hcaldqm::ContainerProf2D _cTimingRMS_FEDuTCA
Definition: LaserTask.h:110
LaserTask::_cSignalRMS_FEDuTCA
hcaldqm::ContainerProf2D _cSignalRMS_FEDuTCA
Definition: LaserTask.h:108
Logger::dqmthrow
void dqmthrow(std::string const &msg) const
Definition: Logger.h:15
LaserTask::_filter_uTCA
hcaldqm::filter::HashFilter _filter_uTCA
Definition: LaserTask.h:63
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
hcaldqm::ContainerProf2D::initialize
virtual void initialize(std::string const &folder, hashfunctions::HashType, Quantity *, Quantity *, Quantity *qz=new ValueQuantity(quantity::fEnergy), int debug=0)
Definition: ContainerProf2D.cc:27
LaserTask::_cTimingMean_FEDVME
hcaldqm::ContainerProf2D _cTimingMean_FEDVME
Definition: LaserTask.h:105
LaserTask::_nevents
int _nevents
Definition: LaserTask.h:67
HODataFrame::size
int size() const
total number of samples in the digi
Definition: HODataFrame.h:27
hcaldqm::Container2D::initialize
virtual void initialize(std::string const &folder, hashfunctions::HashType, Quantity *, Quantity *, Quantity *qz=new ValueQuantity(quantity::fN), int debug=0)
Definition: Container2D.cc:32
hcaldqm::constants::FIBER_uTCA_MIN1
int const FIBER_uTCA_MIN1
Definition: Constants.h:93
LaserTask::_tokHO
edm::EDGetTokenT< HODigiCollection > _tokHO
Definition: LaserTask.h:56
HcalElectronicsMap::allPrecisionId
std::vector< HcalGenericDetId > allPrecisionId() const
Definition: HcalElectronicsMap.cc:168
hcaldqm::DQModule::_logger
Logger _logger
Definition: DQModule.h:70
LaserTask::_tokHF
edm::EDGetTokenT< HFDigiCollection > _tokHF
Definition: LaserTask.h:57
hcaldqm::DQModule::_name
std::string _name
Definition: DQModule.h:57
hcaldqm::ContainerProf1D::book
virtual void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="")
Definition: ContainerProf1D.cc:47
LaserTask::_cSignalvsBX_SubdetPM
hcaldqm::ContainerProf1D _cSignalvsBX_SubdetPM
Definition: LaserTask.h:95
LaserTask::_lowHO
double _lowHO
Definition: LaserTask.h:69
LaserTask::_isApplicable
virtual bool _isApplicable(edm::Event const &)
Definition: LaserTask.cc:476
hcaldqm::constants::CRATE_uTCA_MIN
int const CRATE_uTCA_MIN
Definition: Constants.h:65
LaserTask::_cMissing_FEDuTCA
hcaldqm::Container2D _cMissing_FEDuTCA
Definition: LaserTask.h:115
LaserTask::_cSignalMean_Subdet
hcaldqm::Container1D _cSignalMean_Subdet
Definition: LaserTask.h:81
LaserTask::_tokHBHE
edm::EDGetTokenT< HBHEDigiCollection > _tokHBHE
Definition: LaserTask.h:55
edm::Event::eventAuxiliary
EventAuxiliary const & eventAuxiliary() const
Definition: Event.h:77
LaserTask::_process
virtual void _process(edm::Event const &, edm::EventSetup const &)
Definition: LaserTask.cc:321
LaserTask::_cSignalvsLS_SubdetPM
hcaldqm::ContainerProf1D _cSignalvsLS_SubdetPM
Definition: LaserTask.h:93
LaserTask::_cTimingvsLS_SubdetPM
hcaldqm::ContainerProf1D _cTimingvsLS_SubdetPM
Definition: LaserTask.h:92
HFDataFrame::size
int size() const
total number of samples in the digi
Definition: HFDataFrame.h:26
LaserTask::_tagHBHE
edm::InputTag _tagHBHE
Definition: LaserTask.h:51
hcaldqm::Container2D::book
virtual void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="")
Definition: Container2D.cc:895
hcaldqm::ContainerXXX::book
virtual void book(HcalElectronicsMap const *)
Definition: ContainerXXX.h:99
LaserTask::_cSignalMean_FEDuTCA
hcaldqm::ContainerProf2D _cSignalMean_FEDuTCA
Definition: LaserTask.h:104
hcaldqm::Container1D::book
virtual void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="")
Definition: Container1D.cc:957
hcaldqm::constants::EVENTTYPE_LASER
uint8_t const EVENTTYPE_LASER
Definition: Constants.h:274
LaserTask::_cTimingMean_depth
hcaldqm::ContainerProf2D _cTimingMean_depth
Definition: LaserTask.h:100
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
LaserTask::_ehashmap
hcaldqm::electronicsmap::ElectronicsMap _ehashmap
Definition: LaserTask.h:62
LaserTask::_taguMN
edm::InputTag _taguMN
Definition: LaserTask.h:54
edm::Run::runAuxiliary
RunAuxiliary const & runAuxiliary() const
Definition: Run.h:60
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
edm::EventAuxiliary::id
EventID const & id() const
Definition: EventAuxiliary.h:59
hcaldqm::ContainerXXX::get
virtual STDTYPE & get(HcalDetId const &)
Definition: ContainerXXX.h:249
LaserTask::_tokuMN
edm::EDGetTokenT< HcalUMNioDigi > _tokuMN
Definition: LaserTask.h:58
LaserTask::_cSignalRMS_Subdet
hcaldqm::Container1D _cSignalRMS_Subdet
Definition: LaserTask.h:82
LaserTask::endLuminosityBlock
virtual void endLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &)
Definition: LaserTask.cc:466
hcaldqm::DQModule::_subsystem
std::string _subsystem
Definition: DQModule.h:64
LaserTask::_xTimingSum2
hcaldqm::ContainerXXX< double > _xTimingSum2
Definition: LaserTask.h:78
hcaldqm::ContainerXXX::initialize
virtual void initialize(hashfunctions::HashType, int debug=0)
Definition: ContainerXXX.h:91
hcaldqm::filter::HashFilter::initialize
virtual void initialize(FilterType ftype, HashType htype, std::vector< uint32_t > const &)
Definition: HashFilter.cc:26
hcaldqm::ContainerProf2D::book
virtual void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="")
Definition: ContainerProf2D.cc:51
LaserTask::_cSignalRMS_depth
hcaldqm::ContainerProf2D _cSignalRMS_depth
Definition: LaserTask.h:99
LaserTask::_resetMonitors
virtual void _resetMonitors(hcaldqm::UpdateFreq)
Definition: LaserTask.cc:247
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:94
edm::RunAuxiliary::run
RunNumber_t run() const
Definition: RunAuxiliary.h:41
LaserTask::LaserTask
LaserTask(edm::ParameterSet const &)
Definition: LaserTask.cc:7
hcaldqm::constants::SLOT_uTCA_MIN
int const SLOT_uTCA_MIN
Definition: Constants.h:72
edm::Run
Definition: Run.h:43
LaserTask::_cTimingRMS_depth
hcaldqm::ContainerProf2D _cTimingRMS_depth
Definition: LaserTask.h:101
hcaldqm::utilities::sumQ
double sumQ(DIGI const &digi, double ped, int i=0, int j=3)
Definition: Utilities.h:79
Generated for CMSSW Reference Manual by doxygen 1.8.5