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