test
CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
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):
8  DQTask(ps)
9 {
10  _tagHBHE = ps.getUntrackedParameter<edm::InputTag>("tagHBHE",
11  edm::InputTag("hcalDigis"));
12  _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO",
13  edm::InputTag("hcalDigis"));
14  _tagHF = ps.getUntrackedParameter<edm::InputTag>("tagHF",
15  edm::InputTag("hcalDigis"));
16 
17  _tokHBHE = consumes<HBHEDigiCollection>(_tagHBHE);
18  _tokHO = consumes<HODigiCollection>(_tagHO);
19  _tokHF = consumes<HFDigiCollection>(_tagHF);
20 
21  _cutSumQ_HBHE = ps.getUntrackedParameter<double>("cutSumQ_HBHE", 20);
22  _cutSumQ_HO = ps.getUntrackedParameter<double>("cutSumQ_HO", 20);
23  _cutSumQ_HF = ps.getUntrackedParameter<double>("cutSumQ_HF", 20);
24  _thresh_unihf = ps.getUntrackedParameter<double>("thresh_unihf", 0.2);
25 
26  _vflags.resize(nDigiFlag);
27  _vflags[fUni]=hcaldqm::flag::Flag("UniSlotHF");
28  _vflags[fDigiSize]=hcaldqm::flag::Flag("DigiSize");
29  _vflags[fNChsHF]=hcaldqm::flag::Flag("NChsHF");
30  _vflags[fUnknownIds]=hcaldqm::flag::Flag("UnknownIds");
31 }
32 
33 /* virtual */ void DigiTask::bookHistograms(DQMStore::IBooker& ib,
34  edm::Run const& r, edm::EventSetup const& es)
35 {
36  DQTask::bookHistograms(ib,r,es);
37 
38  // GET WHAT YOU NEED
39  edm::ESHandle<HcalDbService> dbs;
40  es.get<HcalDbRecord>().get(dbs);
41  _emap = dbs->getHcalMapping();
42  std::vector<int> vFEDs = hcaldqm::utilities::getFEDList(_emap);
43  std::vector<int> vFEDsVME = hcaldqm::utilities::getFEDVMEList(_emap);
44  std::vector<int> vFEDsuTCA = hcaldqm::utilities::getFEDuTCAList(_emap);
45  std::vector<uint32_t> vVME;
46  std::vector<uint32_t> vuTCA;
47  std::vector<uint32_t> vFEDHF;
48  vVME.push_back(HcalElectronicsId(constants::FIBERCH_MIN,
49  constants::FIBER_VME_MIN, SPIGOT_MIN, CRATE_VME_MIN).rawId());
50  vuTCA.push_back(HcalElectronicsId(CRATE_uTCA_MIN, SLOT_uTCA_MIN,
51  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
52  _filter_VME.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics,
53  vVME);
54  _filter_uTCA.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics,
55  vuTCA);
56  vFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN,
57  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
58  vFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN,
59  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
60  vFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN,
61  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
62 
63  // initialize filters
64  _filter_FEDHF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fFED,
65  vFEDHF);
66 
67  // push the rawIds of each fed into the vector...
68  for (std::vector<int>::const_iterator it=vFEDsVME.begin();
69  it!=vFEDsVME.end(); ++it)
70  _vhashFEDs.push_back(HcalElectronicsId(
71  constants::FIBERCH_MIN, FIBER_VME_MIN, SPIGOT_MIN,
72  (*it)-FED_VME_MIN).rawId());
73  for (std::vector<int>::const_iterator it=vFEDsuTCA.begin();
74  it!=vFEDsuTCA.end(); ++it)
75  _vhashFEDs.push_back(HcalElectronicsId(
76  hcaldqm::utilities::fed2crate(*it), SLOT_uTCA_MIN, FIBER_uTCA_MIN1,
77  FIBERCH_MIN, false).rawId());
78 
79  // INITIALIZE FIRST
80  _cADC_SubdetPM.initialize(_name, "ADC", hcaldqm::hashfunctions::fSubdetPM,
81  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_128),
82  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
83  _cfC_SubdetPM.initialize(_name, "fC", hcaldqm::hashfunctions::fSubdetPM,
84  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),
85  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
86  _cSumQ_SubdetPM.initialize(_name, "SumQ", hcaldqm::hashfunctions::fSubdetPM,
87  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),
88  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
89  _cSumQ_depth.initialize(_name, "SumQ", hcaldqm::hashfunctions::fdepth,
90  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
91  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
92  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000));
93  _cSumQvsLS_SubdetPM.initialize(_name, "SumQvsLS",
94  hcaldqm::hashfunctions::fSubdetPM,
95  new hcaldqm::quantity::LumiSection(_maxLS),
96  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000));
97  _cShapeCut_FED.initialize(_name, "Shape",
98  hcaldqm::hashfunctions::fFED,
99  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
100  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000));
101  _cTimingCut_SubdetPM.initialize(_name, "TimingCut",
102  hcaldqm::hashfunctions::fSubdetPM,
103  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
104  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
105  _cTimingCut_FEDVME.initialize(_name, "TimingCut",
106  hcaldqm::hashfunctions::fFED,
107  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
108  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
109  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
110  _cTimingCut_FEDuTCA.initialize(_name, "TimingCut",
111  hcaldqm::hashfunctions::fFED,
112  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
113  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
114  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
115  _cTimingCut_ElectronicsVME.initialize(_name, "TimingCut",
116  hcaldqm::hashfunctions::fElectronics,
117  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
118  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
119  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
120  _cTimingCut_ElectronicsuTCA.initialize(_name, "TimingCut",
121  hcaldqm::hashfunctions::fElectronics,
122  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
123  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
124  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
125  _cTimingCutvsLS_FED.initialize(_name, "TimingvsLS",
126  hcaldqm::hashfunctions::fFED,
127  new hcaldqm::quantity::LumiSection(_maxLS),
128  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
129  _cTimingCut_depth.initialize(_name, "TimingCut",
130  hcaldqm::hashfunctions::fdepth,
131  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
132  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
133  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
134 
135  // Occupancy w/o a cut
136  _cOccupancy_FEDVME.initialize(_name, "Occupancy",
137  hcaldqm::hashfunctions::fFED,
138  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
139  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
140  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
141  _cOccupancy_FEDuTCA.initialize(_name, "Occupancy",
142  hcaldqm::hashfunctions::fFED,
143  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
144  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
145  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
146  _cOccupancy_ElectronicsVME.initialize(_name, "Occupancy",
147  hcaldqm::hashfunctions::fElectronics,
148  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
149  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
150  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
151  _cOccupancy_ElectronicsuTCA.initialize(_name, "Occupancy",
152  hcaldqm::hashfunctions::fElectronics,
153  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
154  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
155  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
156  _cOccupancyvsLS_Subdet.initialize(_name, "OccupancyvsLS",
157  hcaldqm::hashfunctions::fSubdet,
158  new hcaldqm::quantity::LumiSection(_maxLS),
159  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to3000));
160  _cOccupancy_depth.initialize(_name, "Occupancy",
161  hcaldqm::hashfunctions::fdepth,
162  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
163  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
164  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
165 
166  // Occupancy w/ a cut
167  _cOccupancyCut_FEDVME.initialize(_name, "OccupancyCut",
168  hcaldqm::hashfunctions::fFED,
169  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
170  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
171  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
172  _cOccupancyCut_FEDuTCA.initialize(_name, "OccupancyCut",
173  hcaldqm::hashfunctions::fFED,
174  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
175  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
176  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
177  _cOccupancyCut_ElectronicsVME.initialize(_name, "OccupancyCut",
178  hcaldqm::hashfunctions::fElectronics,
179  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
180  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
181  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
182  _cOccupancyCut_ElectronicsuTCA.initialize(_name, "OccupancyCut",
183  hcaldqm::hashfunctions::fElectronics,
184  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
185  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
186  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
187  _cOccupancyCutvsLS_Subdet.initialize(_name, "OccupancyCutvsLS",
188  hcaldqm::hashfunctions::fSubdet,
189  new hcaldqm::quantity::LumiSection(_maxLS),
190  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to3000));
191  _cOccupancyCut_depth.initialize(_name, "OccupancyCut",
192  hcaldqm::hashfunctions::fdepth,
193  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
194  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
195  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
196 
197  _cDigiSize_FED.initialize(_name, "DigiSize",
198  hcaldqm::hashfunctions::fFED,
199  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDigiSize),
200  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
201 
202  // INITIALIZE HISTOGRAMS that are only for Online
203  if (_ptype==fOnline)
204  {
205  std::vector<uint32_t> vhashHF;
206  vhashHF.push_back(HcalDetId(HcalForward, 31,1,1).rawId());
207  _filter_HF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fSubdet,
208  vhashHF);
209 
210  // Charge sharing
211  _cQ2Q12CutvsLS_FEDHF.initialize(_name, "Q2Q12vsLS",
212  hcaldqm::hashfunctions::fFED,
213  new hcaldqm::quantity::LumiSection(_maxLS),
214  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fRatio_0to2));
215  _cSumQvsBX_SubdetPM.initialize(_name, "SumQvsBX",
216  hcaldqm::hashfunctions::fSubdetPM,
217  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
218  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000));
219  _cDigiSizevsLS_FED.initialize(_name, "DigiSizevsLS",
220  hcaldqm::hashfunctions::fFED,
221  new hcaldqm::quantity::LumiSection(_maxLS),
222  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDigiSize));
223  _cTimingCutvsiphi_SubdetPM.initialize(_name, "TimingCutvsiphi",
224  hcaldqm::hashfunctions::fSubdetPM,
225  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
226  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
227  _cTimingCutvsieta_Subdet.initialize(_name, "TimingCutvsieta",
228  hcaldqm::hashfunctions::fSubdet,
229  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
230  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200));
231  _cOccupancyvsiphi_SubdetPM.initialize(_name, "Occupancyvsiphi",
232  hcaldqm::hashfunctions::fSubdetPM,
233  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
234  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
235  _cOccupancyvsieta_Subdet.initialize(_name, "Occupancyvsieta",
236  hcaldqm::hashfunctions::fSubdet,
237  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
238  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
239  _cOccupancyCutvsiphi_SubdetPM.initialize(_name, "OccupancyCutvsiphi",
240  hcaldqm::hashfunctions::fSubdetPM,
241  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
242  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
243  _cOccupancyCutvsieta_Subdet.initialize(_name, "OccupancyCutvsieta",
244  hcaldqm::hashfunctions::fSubdet,
245  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
246  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
247  _cOccupancyCutvsLS_Subdet.initialize(_name, "OccupancyCutvsLS",
248  hcaldqm::hashfunctions::fSubdet,
249  new hcaldqm::quantity::LumiSection(_maxLS),
250  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to3000));
251  _cOccupancyCutvsBX_Subdet.initialize(_name, "OccupancyCutvsBX",
252  hcaldqm::hashfunctions::fSubdet,
253  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
254  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to3000));
255 // _cOccupancyCutvsSlotvsLS_HFPM.initialize(_name,
256 // "OccupancyCutvsSlotvsLS", hcaldqm::hashfunctions::fSubdetPM,
257 // new hcaldqm::quantity::LumiSection(_maxLS),
258 // new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
259 // new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
260  _cOccupancyCutvsiphivsLS_SubdetPM.initialize(_name,
261  "OccupancyCutvsiphivsLS", hcaldqm::hashfunctions::fSubdetPM,
262  new hcaldqm::quantity::LumiSection(_maxLS),
263  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
264  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
265  _cSummaryvsLS_FED.initialize(_name, "SummaryvsLS",
266  hcaldqm::hashfunctions::fFED,
267  new hcaldqm::quantity::LumiSection(_maxLS),
268  new hcaldqm::quantity::FlagQuantity(_vflags),
269  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState));
270  _cSummaryvsLS.initialize(_name, "SummaryvsLS",
271  new hcaldqm::quantity::LumiSection(_maxLS),
272  new hcaldqm::quantity::FEDQuantity(vFEDs),
273  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState));
274 
275  _xUniHF.initialize(hcaldqm::hashfunctions::fFEDSlot);
276  _xUni.initialize(hcaldqm::hashfunctions::fFED);
277  _xDigiSize.initialize(hcaldqm::hashfunctions::fFED);
278  _xNChs.initialize(hcaldqm::hashfunctions::fFED);
279  _xNChsNominal.initialize(hcaldqm::hashfunctions::fFED);
280  }
281 
282  // BOOK HISTOGRAMS
283  char cutstr[200];
284  sprintf(cutstr, "_SumQHBHE%dHO%dHF%d", int(_cutSumQ_HBHE),
285  int(_cutSumQ_HO), int(_cutSumQ_HF));
286  char cutstr2[200];
287  sprintf(cutstr2, "_SumQHF%d", int(_cutSumQ_HF));
288 
289  _cADC_SubdetPM.book(ib, _emap, _subsystem);
290 
291  _cfC_SubdetPM.book(ib, _emap, _subsystem);
292  _cSumQ_SubdetPM.book(ib, _emap, _subsystem);
293  _cSumQ_depth.book(ib, _emap, _subsystem);
294  _cSumQvsLS_SubdetPM.book(ib, _emap, _subsystem);
295 
296  _cShapeCut_FED.book(ib, _emap, _subsystem);
297 
298  _cTimingCut_SubdetPM.book(ib, _emap, _subsystem);
299  _cTimingCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
300  _cTimingCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
301  _cTimingCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
302  _cTimingCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
303  _cTimingCutvsLS_FED.book(ib, _emap, _subsystem);
304  _cTimingCut_depth.book(ib, _emap, _subsystem);
305 
306  _cOccupancy_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
307  _cOccupancy_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
308  _cOccupancy_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
309  _cOccupancy_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
310  _cOccupancyvsLS_Subdet.book(ib, _emap, _subsystem);
311  _cOccupancy_depth.book(ib, _emap, _subsystem);
312  _cOccupancyCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
313  _cOccupancyCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
314  _cOccupancyCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
315  _cOccupancyCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
316  _cOccupancyCut_depth.book(ib, _emap, _subsystem);
317 
318  _cDigiSize_FED.book(ib, _emap, _subsystem);
319 
320  // BOOK HISTOGRAMS that are only for Online
321  _ehashmap.initialize(_emap, electronicsmap::fD2EHashMap);
322  _dhashmap.initialize(_emap, electronicsmap::fE2DHashMap);
323  if (_ptype==fOnline)
324  {
325  _cQ2Q12CutvsLS_FEDHF.book(ib, _emap, _filter_FEDHF, _subsystem);
326  _cSumQvsBX_SubdetPM.book(ib, _emap, _subsystem);
327  _cDigiSizevsLS_FED.book(ib, _emap, _subsystem);
328  _cTimingCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
329  _cTimingCutvsieta_Subdet.book(ib, _emap, _subsystem);
330  _cOccupancyCutvsLS_Subdet.book(ib, _emap, _subsystem);
331  _cOccupancyCutvsBX_Subdet.book(ib, _emap, _subsystem);
332  _cOccupancyvsiphi_SubdetPM.book(ib, _emap, _subsystem);
333  _cOccupancyvsieta_Subdet.book(ib, _emap, _subsystem);
334  _cOccupancyCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
335  _cOccupancyCutvsieta_Subdet.book(ib, _emap, _subsystem);
336 // _cOccupancyCutvsSlotvsLS_HFPM.book(ib, _emap, _filter_HF, _subsystem);
337  _cOccupancyCutvsiphivsLS_SubdetPM.book(ib, _emap, _subsystem);
338  _cSummaryvsLS_FED.book(ib, _emap, _subsystem);
339  _cSummaryvsLS.book(ib, _subsystem);
340 
341  _xUniHF.book(_emap, _filter_FEDHF);
342  _xNChs.book(_emap);
343  _xNChsNominal.book(_emap);
344  _xUni.book(_emap);
345  _xDigiSize.book(_emap);
346 
347  // just PER HF FED RECORD THE #CHANNELS
348  // ONLY WAY TO DO THAT AUTOMATICALLY AND W/O HARDCODING 1728
349  // or ANY OTHER VALUES LIKE 2592, 2192
350  std::vector<HcalGenericDetId> gids = _emap->allPrecisionId();
351  for (std::vector<HcalGenericDetId>::const_iterator it=gids.begin();
352  it!=gids.end(); ++it)
353  {
354  if (!it->isHcalDetId())
355  continue;
356  HcalDetId did(it->rawId());
357  if (_xQuality.exists(did))
358  {
359  HcalChannelStatus cs(it->rawId(), _xQuality.get(
360  HcalDetId(*it)));
361  if (
362  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
363  cs.isBitSet(HcalChannelStatus::HcalCellDead))
364  continue;
365  }
366  HcalElectronicsId eid = HcalElectronicsId(_ehashmap.lookup(did));
367  _xNChsNominal.get(eid)++; // he will know the nominal #channels per FED
368  }
369  }
370 
371  // MARK THESE HISTOGRAMS AS LUMI BASED FOR OFFLINE PROCESSING
372  if (_ptype==fOffline)
373  {
374  _cDigiSize_FED.setLumiFlag();
375  _cOccupancy_depth.setLumiFlag();
376  }
377 
378  // book Number of Events vs LS histogram
379  ib.setCurrentFolder(_subsystem+"/RunInfo");
380  meNumEvents1LS = ib.book1D("NumberOfEvents", "NumberOfEvents",
381  1, 0, 1);
382  meNumEvents1LS->setLumiFlag();
383 
384  // book the flag for unknown ids and the online guy as well
385  ib.setCurrentFolder(_subsystem+"/"+_name);
386  meUnknownIds1LS = ib.book1D("UnknownIds", "UnknownIds",
387  1, 0, 1);
388  _unknownIdsPresent = false;
389  meUnknownIds1LS->setLumiFlag();
390 }
391 
392 /* virtual */ void DigiTask::_resetMonitors(hcaldqm::UpdateFreq uf)
393 {
394  DQTask::_resetMonitors(uf);
395 
396  switch(uf)
397  {
398  case hcaldqm::f1LS:
399  _unknownIdsPresent = false;
400  break;
401  case hcaldqm::f50LS:
402  // ^^^ONLINE ONLY!
403  if (_ptype==fOnline)
404  _cOccupancyvsiphi_SubdetPM.reset();
405  // ^^^
406  break;
407  default:
408  break;
409  }
410 }
411 
412 /* virtual */ void DigiTask::_process(edm::Event const& e,
413  edm::EventSetup const&)
414 {
415  edm::Handle<HBHEDigiCollection> chbhe;
416  edm::Handle<HODigiCollection> cho;
417  edm::Handle<HFDigiCollection> chf;
418 
419  if (!e.getByToken(_tokHBHE, chbhe))
420  _logger.dqmthrow("Collection HBHEDigiCollection isn't available"
421  + _tagHBHE.label() + " " + _tagHBHE.instance());
422  if (!e.getByToken(_tokHO, cho))
423  _logger.dqmthrow("Collection HODigiCollection isn't available"
424  + _tagHO.label() + " " + _tagHO.instance());
425  if (!e.getByToken(_tokHF, chf))
426  _logger.dqmthrow("Collection HFDigiCollection isn't available"
427  + _tagHF.label() + " " + _tagHF.instance());
428 
429  // extract some info per event
430  int bx = e.bunchCrossing();
431  meNumEvents1LS->Fill(0.5); // just increment
432 
433  // To fill histograms outside of the loop, you need to determine if there were
434  // any valid det ids first
435  uint32_t rawidValid = 0;
436  uint32_t rawidHBValid = 0;
437  uint32_t rawidHEValid = 0;
438 
439  // HB collection
440  int numChs = 0;
441  int numChsCut = 0;
442  int numChsHE = 0;
443  int numChsCutHE = 0;
444  for (HBHEDigiCollection::const_iterator it=chbhe->begin(); it!=chbhe->end();
445  ++it)
446  {
447  double sumQ = hcaldqm::utilities::sumQ<HBHEDataFrame>(*it, 2.5, 0, it->size()-1);
448 
449  // Explicit check on the DetIds present in the Collection
450  HcalDetId const& did = it->id();
451  uint32_t rawid = _ehashmap.lookup(did);
452  if (rawid==0)
453  {meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
454  HcalElectronicsId const& eid(rawid);
455  if (did.subdet()==HcalBarrel)
456  rawidHBValid = did.rawId();
457  else if (did.subdet()==HcalEndcap)
458  rawidHEValid = did.rawId();
459 
460  // filter out channels that are masked out
461  if (_xQuality.exists(did))
462  {
463  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
464  if (
465  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
466  cs.isBitSet(HcalChannelStatus::HcalCellDead))
467  continue;
468  }
469 
470  _cSumQ_SubdetPM.fill(did, sumQ);
471  _cOccupancy_depth.fill(did);
472  if (_ptype==fOnline)
473  {
474  _cDigiSizevsLS_FED.fill(eid, _currentLS, it->size());
475  it->size()!=constants::DIGISIZE[did.subdet()-1]?
476  _xDigiSize.get(eid)++:_xDigiSize.get(eid)+=0;
477  _cOccupancyvsiphi_SubdetPM.fill(did);
478  _cOccupancyvsieta_Subdet.fill(did);
479  }
480  _cDigiSize_FED.fill(eid, it->size());
481  if (eid.isVMEid())
482  {
483  _cOccupancy_FEDVME.fill(eid);
484  _cOccupancy_ElectronicsVME.fill(eid);
485  }
486  else
487  {
488  _cOccupancy_FEDuTCA.fill(eid);
489  _cOccupancy_ElectronicsuTCA.fill(eid);
490  /*
491  if (!it->validate(0, it->size()))
492  {
493  _cCapIdRots_depth.fill(did);
494  _cCapIdRots_FEDuTCA.fill(eid, 1);
495  }*/
496  }
497 
498  for (int i=0; i<it->size(); i++)
499  {
500  _cADC_SubdetPM.fill(did, it->sample(i).adc());
501  _cfC_SubdetPM.fill(did, it->sample(i).nominal_fC());
502  if (sumQ>_cutSumQ_HBHE)
503  _cShapeCut_FED.fill(eid, i, it->sample(i).nominal_fC());
504  }
505 
506  if (sumQ>_cutSumQ_HBHE)
507  {
508  double timing = hcaldqm::utilities::aveTS<HBHEDataFrame>(*it, 2.5, 0,
509  it->size()-1);
510  _cTimingCut_SubdetPM.fill(did, timing);
511  _cTimingCut_depth.fill(did, timing);
512  _cOccupancyCut_depth.fill(did);
513  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
514  _cSumQ_depth.fill(did, sumQ);
515  _cSumQvsLS_SubdetPM.fill(did, _currentLS, sumQ);
516  if (_ptype==fOnline)
517  {
518  _cSumQvsBX_SubdetPM.fill(did, bx, sumQ);
519  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
520  _cTimingCutvsieta_Subdet.fill(did, timing);
521  _cOccupancyCutvsiphi_SubdetPM.fill(did);
522  _cOccupancyCutvsieta_Subdet.fill(did);
523  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
524  }
525  if (eid.isVMEid())
526  {
527  _cTimingCut_FEDVME.fill(eid, timing);
528  _cTimingCut_ElectronicsVME.fill(eid, timing);
529  _cOccupancyCut_FEDVME.fill(eid);
530  _cOccupancyCut_ElectronicsVME.fill(eid);
531  }
532  else
533  {
534  _cTimingCut_FEDuTCA.fill(eid, timing);
535  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
536  _cOccupancyCut_FEDuTCA.fill(eid);
537  _cOccupancyCut_ElectronicsuTCA.fill(eid);
538  }
539  did.subdet()==HcalBarrel?numChsCut++:numChsCutHE++;
540  }
541  did.subdet()==HcalBarrel?numChs++:numChsHE++;
542  }
543 
544  if (rawidHBValid!=0 && rawidHEValid!=0)
545  {
546  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHBValid), _currentLS,
547  numChs);
548  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHEValid), _currentLS,
549  numChsHE);
550  // ONLINE ONLY!
551  if (_ptype==fOnline)
552  {
553  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHBValid),
554  _currentLS, numChsCut);
555  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHBValid), bx,
556  numChsCut);
557  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHEValid),
558  _currentLS, numChsCutHE);
559  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHEValid), bx,
560  numChsCutHE);
561  }
562  // ^^^ONLINE ONLY!
563  }
564  numChs=0;
565  numChsCut = 0;
566 
567  // reset
568  rawidValid = 0;
569 
570  // HO collection
571  for (HODigiCollection::const_iterator it=cho->begin(); it!=cho->end();
572  ++it)
573  {
574  double sumQ = hcaldqm::utilities::sumQ<HODataFrame>(*it, 8.5, 0, it->size()-1);
575 
576  // Explicit check on the DetIds present in the Collection
577  HcalDetId const& did = it->id();
578  uint32_t rawid = _ehashmap.lookup(did);
579  if (rawid==0)
580  {meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
581  HcalElectronicsId const& eid(rawid);
582  if (did.subdet()==HcalOuter)
583  rawidValid = did.rawId();
584 
585  // filter out channels that are masked out
586  if (_xQuality.exists(did))
587  {
588  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
589  if (
590  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
591  cs.isBitSet(HcalChannelStatus::HcalCellDead))
592  continue;
593  }
594 
595  _cSumQ_SubdetPM.fill(did, sumQ);
596  _cOccupancy_depth.fill(did);
597  if (_ptype==fOnline)
598  {
599  _cDigiSizevsLS_FED.fill(eid, _currentLS, it->size());
600  it->size()!=constants::DIGISIZE[did.subdet()-1]?
601  _xDigiSize.get(eid)++:_xDigiSize.get(eid)+=0;
602  _cOccupancyvsiphi_SubdetPM.fill(did);
603  _cOccupancyvsieta_Subdet.fill(did);
604  }
605  _cDigiSize_FED.fill(eid, it->size());
606  if (eid.isVMEid())
607  {
608  _cOccupancy_FEDVME.fill(eid);
609  _cOccupancy_ElectronicsVME.fill(eid);
610  /*
611  if (!it->validate(0, it->size()))
612  _cCapIdRots_FEDVME.fill(eid, 1);
613  */
614  }
615  else
616  {
617  _cOccupancy_FEDuTCA.fill(eid);
618  _cOccupancy_ElectronicsuTCA.fill(eid);
619  /*
620  if (!it->validate(0, it->size()))
621  _cCapIdRots_FEDuTCA.fill(eid, 1);*/
622  }
623 
624  for (int i=0; i<it->size(); i++)
625  {
626  _cADC_SubdetPM.fill(did, it->sample(i).adc());
627  _cfC_SubdetPM.fill(did, it->sample(i).nominal_fC());
628  if (sumQ>_cutSumQ_HO)
629  _cShapeCut_FED.fill(eid, i, it->sample(i).nominal_fC());
630  }
631 
632  if (sumQ>_cutSumQ_HO)
633  {
634  double timing = hcaldqm::utilities::aveTS<HODataFrame>(*it, 8.5, 0,
635  it->size()-1);
636  _cSumQ_depth.fill(did, sumQ);
637  _cSumQvsLS_SubdetPM.fill(did, _currentLS, sumQ);
638  _cOccupancyCut_depth.fill(did);
639  _cTimingCut_SubdetPM.fill(did, timing);
640  _cTimingCut_depth.fill(did, timing);
641  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
642  if (_ptype==fOnline)
643  {
644  _cSumQvsBX_SubdetPM.fill(did, bx, sumQ);
645  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
646  _cTimingCutvsieta_Subdet.fill(did, timing);
647  _cOccupancyCutvsiphi_SubdetPM.fill(did);
648  _cOccupancyCutvsieta_Subdet.fill(did);
649  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
650  }
651  if (eid.isVMEid())
652  {
653  _cTimingCut_FEDVME.fill(eid, timing);
654  _cTimingCut_ElectronicsVME.fill(eid, timing);
655  _cOccupancyCut_FEDVME.fill(eid);
656  _cOccupancyCut_ElectronicsVME.fill(eid);
657  }
658  else
659  {
660  _cTimingCut_FEDuTCA.fill(eid, timing);
661  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
662  _cOccupancyCut_FEDuTCA.fill(eid);
663  _cOccupancyCut_ElectronicsuTCA.fill(eid);
664  }
665  numChsCut++;
666  }
667  numChs++;
668  }
669 
670  if (rawidValid!=0)
671  {
672  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS,
673  numChs);
674 
675  if (_ptype==fOnline)
676  {
677  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid),
678  _currentLS, numChsCut);
679  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid), bx,
680  numChsCut);
681  }
682  }
683  numChs=0; numChsCut=0;
684 
685  // reset
686  rawidValid = 0;
687 
688  // HF collection
689  for (HFDigiCollection::const_iterator it=chf->begin(); it!=chf->end();
690  ++it)
691  {
692  double sumQ = hcaldqm::utilities::sumQ<HFDataFrame>(*it, 2.5, 0, it->size()-1);
693 
694  // Explicit check on the DetIds present in the Collection
695  HcalDetId const& did = it->id();
696  uint32_t rawid = _ehashmap.lookup(did);
697  if (rawid==0)
698  {meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
699  HcalElectronicsId const& eid(rawid);
700  if (did.subdet()==HcalForward)
701  rawidValid = did.rawId();
702 
703  // filter out channels that are masked out
704  if (_xQuality.exists(did))
705  {
706  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
707  if (
708  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
709  cs.isBitSet(HcalChannelStatus::HcalCellDead))
710  continue;
711  }
712 
713  _cSumQ_SubdetPM.fill(did, sumQ);
714  _cOccupancy_depth.fill(did);
715  if (_ptype==fOnline)
716  {
717  _xNChs.get(eid)++;
718  _cDigiSizevsLS_FED.fill(eid, _currentLS, it->size());
719  it->size()!=constants::DIGISIZE[did.subdet()-1]?
720  _xDigiSize.get(eid)++:_xDigiSize.get(eid)+=0;
721  _cOccupancyvsiphi_SubdetPM.fill(did);
722  _cOccupancyvsieta_Subdet.fill(did);
723  }
724  _cDigiSize_FED.fill(eid, it->size());
725  if (eid.isVMEid())
726  {
727  _cOccupancy_FEDVME.fill(eid);
728  _cOccupancy_ElectronicsVME.fill(eid);
729  /*
730  if (!it->validate(0, it->size()))
731  _cCapIdRots_FEDVME.fill(eid, 1);*/
732  }
733  else
734  {
735  _cOccupancy_FEDuTCA.fill(eid);
736  _cOccupancy_ElectronicsuTCA.fill(eid);
737  /*
738  if (!it->validate(0, it->size()))
739  _cCapIdRots_FEDuTCA.fill(eid, 1);*/
740  }
741 
742  for (int i=0; i<it->size(); i++)
743  {
744  _cADC_SubdetPM.fill(did, it->sample(i).adc());
745  _cfC_SubdetPM.fill(did, it->sample(i).nominal_fC());
746  if (sumQ>_cutSumQ_HF)
747  _cShapeCut_FED.fill(eid, i, it->sample(i).nominal_fC());
748  }
749 
750  if (sumQ>_cutSumQ_HF)
751  {
752  double timing = hcaldqm::utilities::aveTS<HFDataFrame>(*it, 2.5, 0,
753  it->size()-1);
754  double q1 = it->sample(1).nominal_fC()-2.5;
755  double q2 = it->sample(2).nominal_fC()-2.5;
756  double q2q12 = q2/(q1+q2);
757  _cSumQ_depth.fill(did, sumQ);
758  _cSumQvsLS_SubdetPM.fill(did, _currentLS, sumQ);
759  if (_ptype==fOnline)
760  {
761  _cSumQvsBX_SubdetPM.fill(did, bx, sumQ);
762  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
763  _cTimingCutvsieta_Subdet.fill(did, timing);
764  _cOccupancyCutvsiphi_SubdetPM.fill(did);
765  _cOccupancyCutvsieta_Subdet.fill(did);
766  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
767 // _cOccupancyCutvsSlotvsLS_HFPM.fill(did, _currentLS);
768  _xUniHF.get(eid)++;
769  }
770  _cTimingCut_SubdetPM.fill(did, timing);
771  _cTimingCut_depth.fill(did, timing);
772  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
773  _cOccupancyCut_depth.fill(did);
774  if (!eid.isVMEid())
775  if (_ptype==fOnline)
776  _cQ2Q12CutvsLS_FEDHF.fill(eid, _currentLS, q2q12);
777  if (eid.isVMEid())
778  {
779  _cTimingCut_FEDVME.fill(eid, timing);
780  _cTimingCut_ElectronicsVME.fill(eid, timing);
781  _cOccupancyCut_FEDVME.fill(eid);
782  _cOccupancyCut_ElectronicsVME.fill(eid);
783  }
784  else
785  {
786  _cTimingCut_FEDuTCA.fill(eid, timing);
787  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
788  _cOccupancyCut_FEDuTCA.fill(eid);
789  _cOccupancyCut_ElectronicsuTCA.fill(eid);
790  }
791  numChsCut++;
792  }
793  numChs++;
794  }
795 
796  if (rawidValid!=0)
797  {
798  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS,
799  numChs);
800 
801  if (_ptype==fOnline)
802  {
803  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid),
804  _currentLS, numChsCut);
805  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid), bx,
806  numChsCut);
807  }
808  }
809 }
810 
811 /* virtual */ void DigiTask::beginLuminosityBlock(
812  edm::LuminosityBlock const& lb, edm::EventSetup const& es)
813 {
814  DQTask::beginLuminosityBlock(lb, es);
815 }
816 
817 /* virtual */ void DigiTask::endLuminosityBlock(edm::LuminosityBlock const& lb,
818  edm::EventSetup const& es)
819 {
820  if (_ptype!=fOnline)
821  return;
822 
823  for (uintCompactMap::const_iterator it=_xUniHF.begin();
824  it!=_xUniHF.end(); ++it)
825  {
826  uint32_t hash1 = it->first;
827  HcalElectronicsId eid1(hash1);
828  double x1 = it->second;
829 
830  for (uintCompactMap::const_iterator jt=_xUniHF.begin();
831  jt!=_xUniHF.end(); ++jt)
832  {
833  if (jt==it)
834  continue;
835  double x2 = jt->second;
836  if (x2==0)
837  continue;
838  if (x1/x2<_thresh_unihf)
839  _xUni.get(eid1)++;
840  }
841  }
842 
843  for (std::vector<uint32_t>::const_iterator it=_vhashFEDs.begin();
844  it!=_vhashFEDs.end(); ++it)
845  {
846  hcaldqm::flag::Flag fSum("DIGI");
847  HcalElectronicsId eid = HcalElectronicsId(*it);
848 
849  std::vector<uint32_t>::const_iterator cit=std::find(
850  _vcdaqEids.begin(), _vcdaqEids.end(), *it);
851  if (cit==_vcdaqEids.end())
852  {
853  // not @cDAQ
854  for (uint32_t iflag=0; iflag<_vflags.size(); iflag++)
855  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag),
856  int(hcaldqm::flag::fNCDAQ));
857  _cSummaryvsLS.setBinContent(eid, _currentLS, int(hcaldqm::flag::fNCDAQ));
858  continue;
859  }
860 
861  // FED is @cDAQ
862  if (hcaldqm::utilities::isFEDHBHE(eid) || hcaldqm::utilities::isFEDHF(eid) ||
863  hcaldqm::utilities::isFEDHO(eid))
864  {
865  if (_xDigiSize.get(eid)>0)
866  _vflags[fDigiSize]._state = hcaldqm::flag::fBAD;
867  else
868  _vflags[fDigiSize]._state = hcaldqm::flag::fGOOD;
869  if (hcaldqm::utilities::isFEDHF(eid))
870  {
871  double fr = double(_xNChs.get(eid))/double(
872  _xNChsNominal.get(eid)*_evsPerLS);
873  if (_runkeyVal==0 || _runkeyVal==4)
874  {
875  // only for pp or hi
876  if (_xUni.get(eid)>0)
877  _vflags[fUni]._state = hcaldqm::flag::fPROBLEMATIC;
878  else
879  _vflags[fUni]._state = hcaldqm::flag::fGOOD;
880  }
881  if (fr<0.95)
882  _vflags[fNChsHF]._state = hcaldqm::flag::fBAD;
883  else if (fr<1.0)
884  _vflags[fNChsHF]._state = hcaldqm::flag::fPROBLEMATIC;
885  else
886  _vflags[fNChsHF]._state = hcaldqm::flag::fGOOD;
887  }
888  }
889  if (_unknownIdsPresent)
890  _vflags[fUnknownIds]._state = hcaldqm::flag::fBAD;
891  else
892  _vflags[fUnknownIds]._state = hcaldqm::flag::fGOOD;
893 
894  int iflag=0;
895  for (std::vector<hcaldqm::flag::Flag>::iterator ft=_vflags.begin();
896  ft!=_vflags.end(); ++ft)
897  {
898  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, iflag,
899  int(ft->_state));
900  fSum+=(*ft);
901  iflag++;
902 
903  // reset!
904  ft->reset();
905  }
906  _cSummaryvsLS.setBinContent(eid, _currentLS, fSum._state);
907  }
908 
909  _xDigiSize.reset(); _xUniHF.reset(); _xUni.reset();
910  _xNChs.reset();
911 
912  // in the end always do the DQTask::endLumi
913  DQTask::endLuminosityBlock(lb, es);
914 }
915 
916 DEFINE_FWK_MODULE(DigiTask);
917 
DigiTask::_cTimingCutvsiphi_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsiphi_SubdetPM
Definition: DigiTask.h:102
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
DigiTask::meNumEvents1LS
MonitorElement * meNumEvents1LS
Definition: DigiTask.h:145
i
int i
Definition: DBlmapReader.cc:9
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
DigiTask::_cADC_SubdetPM
hcaldqm::Container1D _cADC_SubdetPM
Definition: DigiTask.h:83
DigiTask::_cOccupancyCutvsieta_Subdet
hcaldqm::Container1D _cOccupancyCutvsieta_Subdet
Definition: DigiTask.h:127
fwrapper::cs
auto_ptr< ClusterSequence > cs
Definition: fastjetfortran_madfks.cc:45
DigiTask::_cSumQvsBX_SubdetPM
hcaldqm::ContainerProf1D _cSumQvsBX_SubdetPM
Definition: DigiTask.h:88
hcaldqm::ContainerXXX::exists
virtual bool exists(HcalDetId const &)
Definition: ContainerXXX.h:267
hcaldqm::Container2D::setBinContent
virtual void setBinContent(HcalDetId const &, int)
Definition: Container2D.cc:235
hcaldqm::constants::CRATE_VME_MIN
int const CRATE_VME_MIN
Definition: Constants.h:60
cuy.ib
int ib
Definition: cuy.py:660
hcaldqm::utilities::fed2crate
uint16_t fed2crate(int fed)
Definition: Utilities.cc:12
DigiTask::_tagHO
edm::InputTag _tagHO
Definition: DigiTask.h:47
DigiTask::_cQ2Q12CutvsLS_FEDHF
hcaldqm::ContainerProf1D _cQ2Q12CutvsLS_FEDHF
Definition: DigiTask.h:106
DigiTask::meUnknownIds1LS
MonitorElement * meUnknownIds1LS
Definition: DigiTask.h:146
DigiTask::_tokHBHE
edm::EDGetTokenT< HBHEDigiCollection > _tokHBHE
Definition: DigiTask.h:49
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
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
DigiTask::_tokHF
edm::EDGetTokenT< HFDigiCollection > _tokHF
Definition: DigiTask.h:51
DigiTask::_cutSumQ_HO
double _cutSumQ_HO
Definition: DigiTask.h:53
hcaldqm::ContainerXXX::begin
virtual CompactMap::const_iterator begin()
Definition: ContainerXXX.h:75
DigiTask::_thresh_unihf
double _thresh_unihf
Definition: DigiTask.h:54
hcaldqm::UpdateFreq
UpdateFreq
Definition: DQTask.h:17
hcaldqm::flag::Flag::_state
State _state
Definition: Flag.h:79
edm::SortedCollection< HBHEDataFrame >::const_iterator
std::vector< HBHEDataFrame >::const_iterator const_iterator
Definition: SortedCollection.h:82
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:65
hcaldqm::constants::SPIGOT_MIN
int const SPIGOT_MIN
Definition: Constants.h:85
DigiTask::_cSumQ_depth
hcaldqm::ContainerProf2D _cSumQ_depth
Definition: DigiTask.h:86
DigiTask::_tokHO
edm::EDGetTokenT< HODigiCollection > _tokHO
Definition: DigiTask.h:50
DigiTask::_cOccupancyCutvsiphi_SubdetPM
hcaldqm::Container1D _cOccupancyCutvsiphi_SubdetPM
Definition: DigiTask.h:126
hcaldqm::ContainerXXX::reset
virtual void reset()
Definition: ContainerXXX.h:366
alcazmumu_cfi.filter
tuple filter
Definition: alcazmumu_cfi.py:34
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
DigiTask::_cSummaryvsLS
hcaldqm::ContainerSingle2D _cSummaryvsLS
Definition: DigiTask.h:150
DigiTask::_cTimingCut_FEDVME
hcaldqm::ContainerProf2D _cTimingCut_FEDVME
Definition: DigiTask.h:96
hcaldqm::constants::FIBER_VME_MIN
int const FIBER_VME_MIN
Definition: Constants.h:90
DigiTask::_cOccupancyCutvsiphivsLS_SubdetPM
hcaldqm::Container2D _cOccupancyCutvsiphivsLS_SubdetPM
Definition: DigiTask.h:129
DigiTask::_xDigiSize
hcaldqm::ContainerXXX< uint32_t > _xDigiSize
Definition: DigiTask.h:139
DigiTask::_cTimingCut_FEDuTCA
hcaldqm::ContainerProf2D _cTimingCut_FEDuTCA
Definition: DigiTask.h:97
bookHistograms
void bookHistograms(fwlite::EventContainer &eventCont)
Definition: secvtxMassTemplates.cc:293
q2
double q2[4]
Definition: TauolaWrapper.h:88
hcaldqm::constants::FIBERCH_MIN
int const FIBERCH_MIN
Definition: Constants.h:100
DigiTask::beginLuminosityBlock
virtual void beginLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &)
Definition: DigiTask.cc:811
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
hcaldqm::ContainerSingle2D::book
virtual void book(DQMStore::IBooker &, std::string subsystem="Hcal", std::string aux="")
Definition: ContainerSingle2D.cc:81
hcaldqm::utilities::getFEDVMEList
std::vector< int > getFEDVMEList(HcalElectronicsMap const *)
Definition: Utilities.cc:71
DigiTask::_filter_VME
hcaldqm::filter::HashFilter _filter_VME
Definition: DigiTask.h:76
hcaldqm::ContainerXXX::end
virtual CompactMap::const_iterator end()
Definition: ContainerXXX.h:77
HLT_25ns10e33_v2_cff.InputTag
tuple InputTag
Definition: HLT_25ns10e33_v2_cff.py:58467
hcaldqm::Container2D::fill
virtual void fill(HcalDetId const &)
Definition: Container2D.cc:59
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
DigiTask::_tagHBHE
edm::InputTag _tagHBHE
Definition: DigiTask.h:46
hcaldqm::flag::Flag::reset
void reset()
Definition: Flag.h:76
hcaldqm::ContainerSingle2D::setBinContent
virtual void setBinContent(int, int, int)
Definition: ContainerSingle2D.cc:207
Logger::dqmthrow
void dqmthrow(std::string const &msg) const
Definition: Logger.h:15
hcaldqm::DQTask::_xQuality
ContainerXXX< uint32_t > _xQuality
Definition: DQTask.h:62
hcaldqm::constants::FED_VME_MIN
int const FED_VME_MIN
Definition: Constants.h:47
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
DigiTask::_cOccupancyCut_ElectronicsVME
hcaldqm::Container2D _cOccupancyCut_ElectronicsVME
Definition: DigiTask.h:123
DigiTask::_cShapeCut_FED
hcaldqm::Container1D _cShapeCut_FED
Definition: DigiTask.h:91
DigiTask::_cOccupancyCut_depth
hcaldqm::Container2D _cOccupancyCut_depth
Definition: DigiTask.h:125
DigiTask::_cOccupancyCutvsBX_Subdet
hcaldqm::ContainerProf1D _cOccupancyCutvsBX_Subdet
Definition: DigiTask.h:134
DigiTask::_cutSumQ_HF
double _cutSumQ_HF
Definition: DigiTask.h:53
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
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
DigiTask::_resetMonitors
virtual void _resetMonitors(hcaldqm::UpdateFreq)
Definition: DigiTask.cc:392
DigiTask::_cutSumQ_HBHE
double _cutSumQ_HBHE
Definition: DigiTask.h:53
DigiTask::_cTimingCut_depth
hcaldqm::ContainerProf2D _cTimingCut_depth
Definition: DigiTask.h:101
HcalElectronicsMap::allPrecisionId
std::vector< HcalGenericDetId > allPrecisionId() const
Definition: HcalElectronicsMap.cc:168
hcaldqm::DQModule::_logger
Logger _logger
Definition: DQModule.h:70
DigiTask::_cOccupancy_FEDVME
hcaldqm::Container2D _cOccupancy_FEDVME
Definition: DigiTask.h:110
DigiTask::_cOccupancyCut_FEDVME
hcaldqm::Container2D _cOccupancyCut_FEDVME
Definition: DigiTask.h:121
DigiTask::_cOccupancy_depth
hcaldqm::Container2D _cOccupancy_depth
Definition: DigiTask.h:114
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
DigiTask::_cSummaryvsLS_FED
hcaldqm::Container2D _cSummaryvsLS_FED
Definition: DigiTask.h:149
DigiTask::_filter_FEDHF
hcaldqm::filter::HashFilter _filter_FEDHF
Definition: DigiTask.h:78
hcaldqm::constants::CRATE_uTCA_MIN
int const CRATE_uTCA_MIN
Definition: Constants.h:65
DigiTask::DigiTask
DigiTask(edm::ParameterSet const &)
Definition: DigiTask.cc:7
DigiTask::_cOccupancyvsLS_Subdet
hcaldqm::ContainerProf1D _cOccupancyvsLS_Subdet
Definition: DigiTask.h:132
DigiTask::_cSumQvsLS_SubdetPM
hcaldqm::ContainerProf1D _cSumQvsLS_SubdetPM
Definition: DigiTask.h:87
DigiTask::_xUni
hcaldqm::ContainerXXX< uint32_t > _xUni
Definition: DigiTask.h:140
DigiTask::_cTimingCutvsLS_FED
hcaldqm::ContainerProf1D _cTimingCutvsLS_FED
Definition: DigiTask.h:100
HcalElectronicsId::isVMEid
bool isVMEid() const
Definition: HcalElectronicsId.h:49
DigiTask::_cOccupancyCut_FEDuTCA
hcaldqm::Container2D _cOccupancyCut_FEDuTCA
Definition: DigiTask.h:122
DigiTask::_cOccupancy_ElectronicsVME
hcaldqm::Container2D _cOccupancy_ElectronicsVME
Definition: DigiTask.h:112
DigiTask::_tagHF
edm::InputTag _tagHF
Definition: DigiTask.h:48
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:276
hcaldqm::utilities::isFEDHBHE
bool isFEDHBHE(HcalElectronicsId const &)
Definition: Utilities.cc:124
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
DigiTask::_cDigiSize_FED
hcaldqm::Container1D _cDigiSize_FED
Definition: DigiTask.h:137
hcaldqm::ContainerSingle2D::initialize
virtual void initialize(std::string const &folder, Quantity *, Quantity *, Quantity *qz=new ValueQuantity(quantity::fN), int debug=0)
Definition: ContainerSingle2D.cc:55
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:957
DigiTask::_xNChsNominal
hcaldqm::ContainerXXX< uint32_t > _xNChsNominal
Definition: DigiTask.h:142
DigiTask::_cSumQ_SubdetPM
hcaldqm::Container1D _cSumQ_SubdetPM
Definition: DigiTask.h:85
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
DigiTask::_cOccupancyvsiphi_SubdetPM
hcaldqm::Container1D _cOccupancyvsiphi_SubdetPM
Definition: DigiTask.h:115
hcaldqm::constants::DIGISIZE
int const DIGISIZE[SUBDET_NUM]
Definition: Constants.h:132
DigiTask::_cTimingCutvsieta_Subdet
hcaldqm::ContainerProf1D _cTimingCutvsieta_Subdet
Definition: DigiTask.h:103
DigiTask::_dhashmap
hcaldqm::electronicsmap::ElectronicsMap _dhashmap
Definition: DigiTask.h:73
DigiTask::_vflags
std::vector< hcaldqm::flag::Flag > _vflags
Definition: DigiTask.h:57
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
DigiTask::_cOccupancyCut_ElectronicsuTCA
hcaldqm::Container2D _cOccupancyCut_ElectronicsuTCA
Definition: DigiTask.h:124
DigiTask::_cTimingCut_ElectronicsVME
hcaldqm::ContainerProf2D _cTimingCut_ElectronicsVME
Definition: DigiTask.h:98
hcaldqm::ContainerXXX::get
virtual STDTYPE & get(HcalDetId const &)
Definition: ContainerXXX.h:249
hcaldqm::Container1D::setLumiFlag
virtual void setLumiFlag()
Definition: Container1D.cc:1359
hcaldqm::utilities::getFEDList
std::vector< int > getFEDList(HcalElectronicsMap const *)
Definition: Utilities.cc:47
hcaldqm::utilities::getFEDuTCAList
std::vector< int > getFEDuTCAList(HcalElectronicsMap const *)
Definition: Utilities.cc:97
hcaldqm::DQModule::_subsystem
std::string _subsystem
Definition: DQModule.h:64
DigiTask::_filter_uTCA
hcaldqm::filter::HashFilter _filter_uTCA
Definition: DigiTask.h:77
DigiTask::endLuminosityBlock
virtual void endLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &)
Definition: DigiTask.cc:817
hcaldqm::utilities::isFEDHO
bool isFEDHO(HcalElectronicsId const &)
Definition: Utilities.cc:171
hcaldqm::DQTask::_vcdaqEids
std::vector< uint32_t > _vcdaqEids
Definition: DQTask.h:65
DigiTask::_cDigiSizevsLS_FED
hcaldqm::ContainerProf1D _cDigiSizevsLS_FED
Definition: DigiTask.h:138
DigiTask::_xNChs
hcaldqm::ContainerXXX< uint32_t > _xNChs
Definition: DigiTask.h:141
DigiTask::_ehashmap
hcaldqm::electronicsmap::ElectronicsMap _ehashmap
Definition: DigiTask.h:72
DigiTask::_cOccupancyCutvsLS_Subdet
hcaldqm::ContainerProf1D _cOccupancyCutvsLS_Subdet
Definition: DigiTask.h:133
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
MonitorElement::setLumiFlag
void setLumiFlag(void)
this ME is meant to be stored for each luminosity section
Definition: MonitorElement.h:154
hcaldqm::ContainerProf2D::book
virtual void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="")
Definition: ContainerProf2D.cc:51
DigiTask::_emap
HcalElectronicsMap const * _emap
Definition: DigiTask.h:71
DigiTask::_vhashFEDs
std::vector< uint32_t > _vhashFEDs
Definition: DigiTask.h:68
DigiTask::_filter_HF
hcaldqm::filter::HashFilter _filter_HF
Definition: DigiTask.h:79
DigiTask::bookHistograms
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
Definition: DigiTask.cc:33
DigiTask::_process
virtual void _process(edm::Event const &, edm::EventSetup const &)
Definition: DigiTask.cc:412
DigiTask::_cOccupancyvsieta_Subdet
hcaldqm::Container1D _cOccupancyvsieta_Subdet
Definition: DigiTask.h:116
HcalElectronicsId
Readout chain identification for Hcal.
Definition: HcalElectronicsId.h:32
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:37
hcaldqm::utilities::isFEDHF
bool isFEDHF(HcalElectronicsId const &)
Definition: Utilities.cc:146
DigiTask::_cTimingCut_ElectronicsuTCA
hcaldqm::ContainerProf2D _cTimingCut_ElectronicsuTCA
Definition: DigiTask.h:99
hcaldqm::constants::SLOT_uTCA_MIN
int const SLOT_uTCA_MIN
Definition: Constants.h:72
DigiTask::_cTimingCut_SubdetPM
hcaldqm::Container1D _cTimingCut_SubdetPM
Definition: DigiTask.h:95
edm::Run
Definition: Run.h:43
DigiTask::_unknownIdsPresent
bool _unknownIdsPresent
Definition: DigiTask.h:147
DigiTask::_cOccupancy_FEDuTCA
hcaldqm::Container2D _cOccupancy_FEDuTCA
Definition: DigiTask.h:111
DigiTask::_xUniHF
hcaldqm::ContainerXXX< uint32_t > _xUniHF
Definition: DigiTask.h:140
hcaldqm::utilities::sumQ
double sumQ(DIGI const &digi, double ped, int i=0, int j=3)
Definition: Utilities.h:79
DigiTask::_cfC_SubdetPM
hcaldqm::Container1D _cfC_SubdetPM
Definition: DigiTask.h:84
DigiTask::_cOccupancy_ElectronicsuTCA
hcaldqm::Container2D _cOccupancy_ElectronicsuTCA
Definition: DigiTask.h:113
Generated for CMSSW Reference Manual by doxygen 1.8.5