CMS 3D CMS Logo
DigiTask.cc
Go to the documentation of this file.
1 #include "DQM/HcalTasks/interface/DigiTask.h"
2 
3 using namespace hcaldqm;
4 using namespace hcaldqm::constants;
5 using namespace hcaldqm::filter;
6 
7 DigiTask::DigiTask(edm::ParameterSet const& ps):
8  DQTask(ps)
9 {
10  _tagHBHE = ps.getUntrackedParameter<edm::InputTag>("tagHBHE",
11  edm::InputTag("hcalDigis"));
12  _tagHEP17 = ps.getUntrackedParameter<edm::InputTag>("tagHEP17",
13  edm::InputTag("hcalDigis"));
14  _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO",
15  edm::InputTag("hcalDigis"));
16  _tagHF = ps.getUntrackedParameter<edm::InputTag>("tagHF",
17  edm::InputTag("hcalDigis"));
18 
19  _tokHBHE = consumes<HBHEDigiCollection>(_tagHBHE);
20  _tokHEP17 = consumes<QIE11DigiCollection>(_tagHEP17);
21  _tokHO = consumes<HODigiCollection>(_tagHO);
22  _tokHF = consumes<QIE10DigiCollection>(_tagHF);
23 
24  _cutSumQ_HBHE = ps.getUntrackedParameter<double>("cutSumQ_HBHE", 20);
25  _cutSumQ_HEP17 = ps.getUntrackedParameter<double>("cutSumQ_HEP17", 20);
26  _cutSumQ_HO = ps.getUntrackedParameter<double>("cutSumQ_HO", 20);
27  _cutSumQ_HF = ps.getUntrackedParameter<double>("cutSumQ_HF", 20);
28  _thresh_unihf = ps.getUntrackedParameter<double>("thresh_unihf", 0.2);
29 
30  _vflags.resize(nDigiFlag);
31  _vflags[fUni]=hcaldqm::flag::Flag("UniSlotHF");
32  _vflags[fDigiSize]=hcaldqm::flag::Flag("DigiSize");
33  _vflags[fNChsHF]=hcaldqm::flag::Flag("NChsHF");
34  _vflags[fUnknownIds]=hcaldqm::flag::Flag("UnknownIds");
35 
36  _qie10InConditions = ps.getUntrackedParameter<bool>("qie10InConditions", true);
37 }
38 
39 /* virtual */ void DigiTask::bookHistograms(DQMStore::IBooker& ib,
40  edm::Run const& r, edm::EventSetup const& es)
41 {
42  DQTask::bookHistograms(ib,r,es);
43 
44  // GET WHAT YOU NEED
45  edm::ESHandle<HcalDbService> dbs;
46  es.get<HcalDbRecord>().get(dbs);
47  _emap = dbs->getHcalMapping();
48  std::vector<uint32_t> vVME;
49  std::vector<uint32_t> vuTCA;
50  vVME.push_back(HcalElectronicsId(constants::FIBERCH_MIN,
51  constants::FIBER_VME_MIN, SPIGOT_MIN, CRATE_VME_MIN).rawId());
52  vuTCA.push_back(HcalElectronicsId(CRATE_uTCA_MIN, SLOT_uTCA_MIN,
53  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
54  _filter_VME.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics,
55  vVME);
56  _filter_uTCA.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics,
57  vuTCA);
58 
59  // Filters for HEP17 and HF, aka QIE10/11
60  std::vector<uint32_t> vhashHF;
61  vhashHF.push_back(hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdet](HcalDetId(HcalForward, 29,1,1)));
62  _filter_HF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fSubdet,
63  vhashHF);
64  _filter_notHF.initialize(filter::fFilter, hcaldqm::hashfunctions::fSubdet,
65  vhashHF);
66 
67  std::vector<uint32_t> vhashHEP17;
68  vhashHEP17.push_back(HcalDetId(HcalEndcap, 1, 63, 1));
69  vhashHEP17.push_back(HcalDetId(HcalEndcap, 1, 64, 1));
70  vhashHEP17.push_back(HcalDetId(HcalEndcap, 1, 65, 1));
71  vhashHEP17.push_back(HcalDetId(HcalEndcap, 1, 66, 1));
72  _filter_HEP17.initialize(filter::fPreserver, hcaldqm::hashfunctions::fSubdetPMiphi, vhashHEP17);
73 
74  // INITIALIZE FIRST
75  _cADC_SubdetPM.initialize(_name, "ADC", hcaldqm::hashfunctions::fSubdetPM,
76  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_128),
77  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
78  _cfC_SubdetPM.initialize(_name, "fC", hcaldqm::hashfunctions::fSubdetPM,
79  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),
80  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
81  _cSumQ_SubdetPM.initialize(_name, "SumQ", hcaldqm::hashfunctions::fSubdetPM,
82  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),
83  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
84  _cSumQ_depth.initialize(_name, "SumQ", hcaldqm::hashfunctions::fdepth,
85  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
86  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
87  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_400000, true),0);
88  _cSumQvsLS_SubdetPM.initialize(_name, "SumQvsLS",
89  hcaldqm::hashfunctions::fSubdetPM,
90  new hcaldqm::quantity::LumiSection(_maxLS),
91  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),0);
92 
93  _cADC_SubdetPM_HF.initialize(_name, "ADC", hcaldqm::hashfunctions::fSubdetPM,
94  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
95  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
96  _cfC_SubdetPM_HF.initialize(_name, "fC", hcaldqm::hashfunctions::fSubdetPM,
97  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_10000),
98  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
99  _cSumQ_SubdetPM_HF.initialize(_name, "SumQ", hcaldqm::hashfunctions::fSubdetPM,
100  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_10000),
101  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
102  _cSumQvsLS_SubdetPM_HF.initialize(_name, "SumQvsLS",
103  hcaldqm::hashfunctions::fSubdetPM,
104  new hcaldqm::quantity::LumiSection(_maxLS),
105  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_10000),0);
106 
107  _cTimingCut_SubdetPM.initialize(_name, "TimingCut",
108  hcaldqm::hashfunctions::fSubdetPM,
109  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
110  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
111  _cTimingCut_depth.initialize(_name, "TimingCut",
112  hcaldqm::hashfunctions::fdepth,
113  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
114  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
115  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
116  _cTimingCutvsLS_SubdetPM.initialize(_name, "TimingvsLS",
117  hcaldqm::hashfunctions::fSubdetPM,
118  new hcaldqm::quantity::LumiSection(_maxLS),
119  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
120 
121  // Occupancy w/o a cut
122  _cOccupancyvsLS_Subdet.initialize(_name, "OccupancyvsLS",
123  hcaldqm::hashfunctions::fSubdet,
124  new hcaldqm::quantity::LumiSection(_maxLS),
125  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),0);
126  _cOccupancy_depth.initialize(_name, "Occupancy",
127  hcaldqm::hashfunctions::fdepth,
128  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
129  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
130  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
131 
132  // Occupancy w/ a cut
133  _cOccupancyCutvsLS_Subdet.initialize(_name, "OccupancyCutvsLS",
134  hcaldqm::hashfunctions::fSubdet,
135  new hcaldqm::quantity::LumiSection(_maxLS),
136  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),0);
137  _cOccupancyCut_depth.initialize(_name, "OccupancyCut",
138  hcaldqm::hashfunctions::fdepth,
139  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
140  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
141  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
142 
143  // Digi size
144  _cDigiSize_Crate.initialize(_name, "DigiSize",
145  hcaldqm::hashfunctions::fCrate,
146  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDigiSize),
147  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
148  _cADCvsTS_SubdetPM.initialize(_name, "ADCvsTS",
149  hcaldqm::hashfunctions::fSubdetPM,
150  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
151  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_128),
152  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
153  _cADCvsTS_SubdetPM_HF.initialize(_name, "ADCvsTS",
154  hcaldqm::hashfunctions::fSubdetPM,
155  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
156  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
157  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
158  _cLETDCTimevsADC_SubdetPM.initialize(_name, "LETDCTimevsADC",
159  hcaldqm::hashfunctions::fSubdetPM,
160  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
161  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250),
162  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
163  _cLETDCvsADC_SubdetPM.initialize(_name, "LETDCvsADC",
164  hcaldqm::hashfunctions::fSubdetPM,
165  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
166  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
167  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
168  _cLETDCvsTS_SubdetPM.initialize(_name, "LETDCvsTS",
169  hcaldqm::hashfunctions::fSubdetPM,
170  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
171  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
172  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
173  _cLETDCTime_SubdetPM.initialize(_name, "LETDCTime",
174  hcaldqm::hashfunctions::fSubdetPM,
175  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250),
176  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
177 
178  if (_ptype == fOnline || _ptype == fLocal) {
179  _cOccupancy_Crate.initialize(_name,
180  "Occupancy", hashfunctions::fCrate,
181  new quantity::ElectronicsQuantity(quantity::fSlotuTCA),
182  new quantity::ElectronicsQuantity(quantity::fFiberuTCAFiberCh),
183  new quantity::ValueQuantity(quantity::fN),0);
184  _cOccupancy_CrateSlot.initialize(_name,
185  "Occupancy", hashfunctions::fCrateSlot,
186  new quantity::ElectronicsQuantity(quantity::fFiberuTCA),
187  new quantity::ElectronicsQuantity(quantity::fFiberCh),
188  new quantity::ValueQuantity(quantity::fN),0);
189  }
190 
191  // INITIALIZE HISTOGRAMS that are only for Online
192  if (_ptype==fOnline)
193  {
194  // Charge sharing
195  _cQ2Q12CutvsLS_FEDHF.initialize(_name, "Q2Q12vsLS",
196  hcaldqm::hashfunctions::fFED,
197  new hcaldqm::quantity::LumiSection(_maxLS),
198  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fRatio_0to2),0);
199  _cSumQvsBX_SubdetPM.initialize(_name, "SumQvsBX",
200  hcaldqm::hashfunctions::fSubdetPM,
201  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
202  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),0);
203  _cSumQvsBX_SubdetPM_HF.initialize(_name, "SumQvsBX",
204  hcaldqm::hashfunctions::fSubdetPM,
205  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
206  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_10000),0);
207  _cDigiSizevsLS_FED.initialize(_name, "DigiSizevsLS",
208  hcaldqm::hashfunctions::fFED,
209  new hcaldqm::quantity::LumiSection(_maxLS),
210  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDigiSize),0);
211  _cTimingCutvsiphi_SubdetPM.initialize(_name, "TimingCutvsiphi",
212  hcaldqm::hashfunctions::fSubdetPM,
213  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
214  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
215  _cTimingCutvsieta_Subdet.initialize(_name, "TimingCutvsieta",
216  hcaldqm::hashfunctions::fSubdet,
217  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
218  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
219  _cOccupancyvsiphi_SubdetPM.initialize(_name, "Occupancyvsiphi",
220  hcaldqm::hashfunctions::fSubdetPM,
221  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
222  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
223  _cOccupancyvsieta_Subdet.initialize(_name, "Occupancyvsieta",
224  hcaldqm::hashfunctions::fSubdet,
225  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
226  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
227  _cOccupancyCutvsiphi_SubdetPM.initialize(_name, "OccupancyCutvsiphi",
228  hcaldqm::hashfunctions::fSubdetPM,
229  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
230  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
231  _cOccupancyCutvsieta_Subdet.initialize(_name, "OccupancyCutvsieta",
232  hcaldqm::hashfunctions::fSubdet,
233  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
234  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
235  _cOccupancyCutvsLS_Subdet.initialize(_name, "OccupancyCutvsLS",
236  hcaldqm::hashfunctions::fSubdet,
237  new hcaldqm::quantity::LumiSection(_maxLS),
238  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),0);
239  _cOccupancyCutvsBX_Subdet.initialize(_name, "OccupancyCutvsBX",
240  hcaldqm::hashfunctions::fSubdet,
241  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
242  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),0);
243 // _cOccupancyCutvsSlotvsLS_HFPM.initialize(_name,
244 // "OccupancyCutvsSlotvsLS", hcaldqm::hashfunctions::fSubdetPM,
245 // new hcaldqm::quantity::LumiSection(_maxLS),
246 // new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
247 // new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
248  _cOccupancyCutvsiphivsLS_SubdetPM.initialize(_name,
249  "OccupancyCutvsiphivsLS", hcaldqm::hashfunctions::fSubdetPM,
250  new hcaldqm::quantity::LumiSection(_maxLS),
251  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
252  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
253  }
254  if (_ptype != fOffline) { // hidefed2crate
255  std::vector<int> vFEDs = hcaldqm::utilities::getFEDList(_emap);
256  std::vector<int> vFEDsVME = hcaldqm::utilities::getFEDVMEList(_emap);
257  std::vector<int> vFEDsuTCA = hcaldqm::utilities::getFEDuTCAList(_emap);
258 
259  std::vector<uint32_t> vFEDHF;
260  vFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN,
261  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
262  vFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN+6,
263  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
264  vFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN,
265  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
266  vFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN+6,
267  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
268  vFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN,
269  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
270  vFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN+6,
271  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
272 
273  // initialize filters
274  _filter_FEDHF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fFED,
275  vFEDHF);
276 
277  // push the rawIds of each fed into the vector...
278  for (std::vector<int>::const_iterator it=vFEDsVME.begin();
279  it!=vFEDsVME.end(); ++it)
280  _vhashFEDs.push_back(HcalElectronicsId(
281  constants::FIBERCH_MIN, FIBER_VME_MIN, SPIGOT_MIN,
282  (*it)-FED_VME_MIN).rawId());
283  for (std::vector<int>::const_iterator it=vFEDsuTCA.begin();
284  it!=vFEDsuTCA.end(); ++it)
285  {
286  std::pair<uint16_t, uint16_t> cspair = utilities::fed2crate(*it);
287  _vhashFEDs.push_back(HcalElectronicsId(
288  cspair.first, cspair.second, FIBER_uTCA_MIN1,
289  FIBERCH_MIN, false).rawId());
290  }
291 
292  _cShapeCut_FED.initialize(_name, "ShapeCut",
293  hcaldqm::hashfunctions::fFED,
294  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
295  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),0);
296 
297  _cTimingCut_FEDVME.initialize(_name, "TimingCut",
298  hcaldqm::hashfunctions::fFED,
299  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
300  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
301  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
302  _cTimingCut_FEDuTCA.initialize(_name, "TimingCut",
303  hcaldqm::hashfunctions::fFED,
304  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
305  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
306  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
307  _cTimingCut_ElectronicsVME.initialize(_name, "TimingCut",
308  hcaldqm::hashfunctions::fElectronics,
309  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
310  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
311  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
312  _cTimingCut_ElectronicsuTCA.initialize(_name, "TimingCut",
313  hcaldqm::hashfunctions::fElectronics,
314  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
315  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
316  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
317  _cTimingCutvsLS_FED.initialize(_name, "TimingvsLS",
318  hcaldqm::hashfunctions::fFED,
319  new hcaldqm::quantity::LumiSection(_maxLS),
320  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
321 
322  _cOccupancy_FEDVME.initialize(_name, "Occupancy",
323  hcaldqm::hashfunctions::fFED,
324  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
325  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
326  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
327  _cOccupancy_FEDuTCA.initialize(_name, "Occupancy",
328  hcaldqm::hashfunctions::fFED,
329  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
330  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
331  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
332  _cOccupancy_ElectronicsVME.initialize(_name, "Occupancy",
333  hcaldqm::hashfunctions::fElectronics,
334  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
335  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
336  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
337  _cOccupancy_ElectronicsuTCA.initialize(_name, "Occupancy",
338  hcaldqm::hashfunctions::fElectronics,
339  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
340  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
341  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
342 
343  _cOccupancyCut_FEDVME.initialize(_name, "OccupancyCut",
344  hcaldqm::hashfunctions::fFED,
345  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
346  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
347  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
348  _cOccupancyCut_FEDuTCA.initialize(_name, "OccupancyCut",
349  hcaldqm::hashfunctions::fFED,
350  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
351  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
352  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
353  _cOccupancyCut_ElectronicsVME.initialize(_name, "OccupancyCut",
354  hcaldqm::hashfunctions::fElectronics,
355  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
356  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
357  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
358  _cOccupancyCut_ElectronicsuTCA.initialize(_name, "OccupancyCut",
359  hcaldqm::hashfunctions::fElectronics,
360  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
361  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
362  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
363 
364  _cDigiSize_FED.initialize(_name, "DigiSize",
365  hcaldqm::hashfunctions::fFED,
366  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDigiSize),
367  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
368 
369  if (_ptype == fOnline) {
370  _cSummaryvsLS_FED.initialize(_name, "SummaryvsLS",
371  hcaldqm::hashfunctions::fFED,
372  new hcaldqm::quantity::LumiSection(_maxLS),
373  new hcaldqm::quantity::FlagQuantity(_vflags),
374  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),0);
375  _cSummaryvsLS.initialize(_name, "SummaryvsLS",
376  new hcaldqm::quantity::LumiSection(_maxLS),
377  new hcaldqm::quantity::FEDQuantity(vFEDs),
378  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),0);
379 
380  _xUniHF.initialize(hcaldqm::hashfunctions::fFEDSlot);
381  _xUni.initialize(hcaldqm::hashfunctions::fFED);
382  _xDigiSize.initialize(hcaldqm::hashfunctions::fFED);
383  _xNChs.initialize(hcaldqm::hashfunctions::fFED);
384  _xNChsNominal.initialize(hcaldqm::hashfunctions::fFED);
385  }
386  }
387 
388  // BOOK HISTOGRAMS
389  char cutstr[200];
390  sprintf(cutstr, "_SumQHBHE%dHO%dHF%d", int(_cutSumQ_HBHE),
391  int(_cutSumQ_HO), int(_cutSumQ_HF));
392  char cutstr2[200];
393  sprintf(cutstr2, "_SumQHF%d", int(_cutSumQ_HF));
394 
395  _cADC_SubdetPM.book(ib, _emap, _filter_notHF, _subsystem);
396  _cADC_SubdetPM_HF.book(ib, _emap, _filter_HF, _subsystem);
397  _cfC_SubdetPM.book(ib, _emap, _filter_notHF, _subsystem);
398  _cfC_SubdetPM_HF.book(ib, _emap, _filter_HF, _subsystem);
399  _cSumQ_SubdetPM.book(ib, _emap, _filter_notHF, _subsystem);
400  _cSumQ_SubdetPM_HF.book(ib, _emap, _filter_HF, _subsystem);
401  _cSumQ_depth.book(ib, _emap, _subsystem);
402  _cSumQvsLS_SubdetPM.book(ib, _emap, _filter_notHF, _subsystem);
403  _cSumQvsLS_SubdetPM_HF.book(ib, _emap, _filter_HF, _subsystem);
404  _cDigiSize_Crate.book(ib, _emap, _subsystem);
405  _cADCvsTS_SubdetPM.book(ib, _emap, _filter_notHF, _subsystem);
406  _cADCvsTS_SubdetPM_HF.book(ib, _emap, _filter_HF, _subsystem);
407 
408  if (_ptype != fOffline) { // hidefed2crate
409  _cShapeCut_FED.book(ib, _emap, _subsystem);
410  _cTimingCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
411  _cTimingCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
412  _cTimingCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
413  _cTimingCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
414  _cTimingCutvsLS_FED.book(ib, _emap, _subsystem);
415  _cOccupancy_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
416  _cOccupancy_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
417  _cOccupancy_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
418  _cOccupancy_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
419  _cOccupancyCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
420  _cOccupancyCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
421  _cOccupancyCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
422  _cOccupancyCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
423  _cDigiSize_FED.book(ib, _emap, _subsystem);
424  }
425 
426  _cTimingCut_SubdetPM.book(ib, _emap, _subsystem);
427  _cTimingCut_depth.book(ib, _emap, _subsystem);
428  _cTimingCutvsLS_SubdetPM.book(ib, _emap, _subsystem);
429 
430  _cOccupancyvsLS_Subdet.book(ib, _emap, _subsystem);
431  _cOccupancy_depth.book(ib, _emap, _subsystem);
432  _cOccupancyCut_depth.book(ib, _emap, _subsystem);
433 
434  _cLETDCTimevsADC_SubdetPM.book(ib, _emap, _filter_HF, _subsystem);
435  _cLETDCvsADC_SubdetPM.book(ib, _emap, _filter_HF, _subsystem);
436  _cLETDCvsTS_SubdetPM.book(ib, _emap, _filter_HF, _subsystem);
437  _cLETDCTime_SubdetPM.book(ib, _emap, _filter_HF, _subsystem);
438 
439  // BOOK HISTOGRAMS that are only for Online
440  _ehashmap.initialize(_emap, electronicsmap::fD2EHashMap);
441  _dhashmap.initialize(_emap, electronicsmap::fE2DHashMap);
442 
443  if (_ptype == fOnline || _ptype == fLocal) {
444  _cOccupancy_Crate.book(ib, _emap, _subsystem);
445  _cOccupancy_CrateSlot.book(ib, _emap, _subsystem);
446  }
447 
448  if (_ptype==fOnline)
449  {
450  _cQ2Q12CutvsLS_FEDHF.book(ib, _emap, _filter_FEDHF, _subsystem);
451  _cSumQvsBX_SubdetPM.book(ib, _emap, _filter_notHF, _subsystem);
452  _cSumQvsBX_SubdetPM_HF.book(ib, _emap, _filter_HF, _subsystem);
453  _cDigiSizevsLS_FED.book(ib, _emap, _subsystem);
454  _cTimingCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
455  _cTimingCutvsieta_Subdet.book(ib, _emap, _subsystem);
456  _cOccupancyCutvsLS_Subdet.book(ib, _emap, _subsystem);
457  _cOccupancyCutvsBX_Subdet.book(ib, _emap, _subsystem);
458  _cOccupancyvsiphi_SubdetPM.book(ib, _emap, _subsystem);
459  _cOccupancyvsieta_Subdet.book(ib, _emap, _subsystem);
460  _cOccupancyCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
461  _cOccupancyCutvsieta_Subdet.book(ib, _emap, _subsystem);
462 // _cOccupancyCutvsSlotvsLS_HFPM.book(ib, _emap, _filter_HF, _subsystem);
463  _cOccupancyCutvsiphivsLS_SubdetPM.book(ib, _emap, _subsystem);
464  _cSummaryvsLS_FED.book(ib, _emap, _subsystem);
465  _cSummaryvsLS.book(ib, _subsystem);
466 
467  _xUniHF.book(_emap, _filter_FEDHF);
468  _xNChs.book(_emap);
469  _xNChsNominal.book(_emap);
470  _xUni.book(_emap);
471  _xDigiSize.book(_emap);
472 
473  // just PER HF FED RECORD THE #CHANNELS
474  // ONLY WAY TO DO THAT AUTOMATICALLY AND W/O HARDCODING 1728
475  // or ANY OTHER VALUES LIKE 2592, 2192
476  std::vector<HcalGenericDetId> gids = _emap->allPrecisionId();
477  for (std::vector<HcalGenericDetId>::const_iterator it=gids.begin();
478  it!=gids.end(); ++it)
479  {
480  if (!it->isHcalDetId())
481  continue;
482  HcalDetId did(it->rawId());
483  if (_xQuality.exists(did))
484  {
485  HcalChannelStatus cs(it->rawId(), _xQuality.get(
486  HcalDetId(*it)));
487  if (
488  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
489  cs.isBitSet(HcalChannelStatus::HcalCellDead))
490  continue;
491  }
492  HcalElectronicsId eid = HcalElectronicsId(_ehashmap.lookup(did));
493  _xNChsNominal.get(eid)++; // he will know the nominal #channels per FED
494  }
495  }
496 
497  // MARK THESE HISTOGRAMS AS LUMI BASED FOR OFFLINE PROCESSING
498  if (_ptype==fOffline)
499  {
500  _cDigiSize_Crate.setLumiFlag();
501  //_cDigiSize_FED.setLumiFlag();
502  _cOccupancy_depth.setLumiFlag();
503  }
504 
505  // book Number of Events vs LS histogram
506  ib.setCurrentFolder(_subsystem+"/RunInfo");
507  meNumEvents1LS = ib.book1D("NumberOfEvents", "NumberOfEvents",
508  1, 0, 1);
509  meNumEvents1LS->setLumiFlag();
510 
511  // book the flag for unknown ids and the online guy as well
512  ib.setCurrentFolder(_subsystem+"/"+_name);
513  meUnknownIds1LS = ib.book1D("UnknownIds", "UnknownIds",
514  1, 0, 1);
515  _unknownIdsPresent = false;
516  meUnknownIds1LS->setLumiFlag();
517 }
518 
519 /* virtual */ void DigiTask::_resetMonitors(hcaldqm::UpdateFreq uf)
520 {
521  DQTask::_resetMonitors(uf);
522 
523  switch(uf)
524  {
525  case hcaldqm::f1LS:
526  _unknownIdsPresent = false;
527  break;
528  case hcaldqm::f50LS:
529  // ^^^ONLINE ONLY!
530  if (_ptype==fOnline)
531  _cOccupancyvsiphi_SubdetPM.reset();
532  // ^^^
533  break;
534  default:
535  break;
536  }
537 }
538 
539 /* virtual */ void DigiTask::_process(edm::Event const& e,
540  edm::EventSetup const&)
541 {
542  edm::Handle<HBHEDigiCollection> chbhe;
543  edm::Handle<QIE11DigiCollection> chep17;
544  edm::Handle<HODigiCollection> cho;
545  edm::Handle<QIE10DigiCollection> chf;
546 
547  if (!e.getByToken(_tokHBHE, chbhe))
548  _logger.dqmthrow("Collection HBHEDigiCollection isn't available"
549  + _tagHBHE.label() + " " + _tagHBHE.instance());
550  if (!e.getByToken(_tokHEP17, chep17))
551  _logger.dqmthrow("Collection HEP17DigiCollection isn't available"
552  + _tagHEP17.label() + " " + _tagHEP17.instance());
553  if (!e.getByToken(_tokHO, cho))
554  _logger.dqmthrow("Collection HODigiCollection isn't available"
555  + _tagHO.label() + " " + _tagHO.instance());
556  if (!e.getByToken(_tokHF, chf))
557  _logger.dqmthrow("Collection QIE10DigiCollection isn't available"
558  + _tagHF.label() + " " + _tagHF.instance());
559 
560  // extract some info per event
561  int bx = e.bunchCrossing();
562  meNumEvents1LS->Fill(0.5); // just increment
563 
564  // To fill histograms outside of the loop, you need to determine if there were
565  // any valid det ids first
566  uint32_t rawidValid = 0;
567  uint32_t rawidHBValid = 0;
568  uint32_t rawidHEValid = 0;
569 
570  // HB collection
571  int numChs = 0;
572  int numChsCut = 0;
573  int numChsHE = 0;
574  int numChsCutHE = 0;
575  for (HBHEDigiCollection::const_iterator it=chbhe->begin(); it!=chbhe->end();
576  ++it)
577  {
578  // Explicit check on the DetIds present in the Collection
579  HcalDetId const& did = it->id();
580  uint32_t rawid = _ehashmap.lookup(did);
581  if (rawid==0)
582  {meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
583  HcalElectronicsId const& eid(rawid);
584  if (did.subdet()==HcalBarrel)
585  rawidHBValid = did.rawId();
586  else if (did.subdet()==HcalEndcap)
587  rawidHEValid = did.rawId();
588 
589  //double sumQ = hcaldqm::utilities::sumQ<HBHEDataFrame>(*it, 2.5, 0, it->size()-1);
590  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<HBHEDataFrame>(_dbService, did, *it);
591  double sumQ = hcaldqm::utilities::sumQDB<HBHEDataFrame>(_dbService, digi_fC, did, *it, 0, it->size()-1);
592 
593  // filter out channels that are masked out
594  if (_xQuality.exists(did))
595  {
596  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
597  if (
598  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
599  cs.isBitSet(HcalChannelStatus::HcalCellDead))
600  continue;
601  }
602 
603  _cSumQ_SubdetPM.fill(did, sumQ);
604  _cOccupancy_depth.fill(did);
605  if (_ptype == fOnline || _ptype == fLocal) {
606  _cOccupancy_Crate.fill(eid);
607  _cOccupancy_CrateSlot.fill(eid);
608  }
609  if (_ptype==fOnline)
610  {
611  _cDigiSizevsLS_FED.fill(eid, _currentLS, it->size());
612  it->size()!=constants::DIGISIZE[did.subdet()-1]?
613  _xDigiSize.get(eid)++:_xDigiSize.get(eid)+=0;
614  _cOccupancyvsiphi_SubdetPM.fill(did);
615  _cOccupancyvsieta_Subdet.fill(did);
616  }
617  _cDigiSize_Crate.fill(eid, it->size());
618  if (_ptype != fOffline) { // hidefed2crate
619  _cDigiSize_FED.fill(eid, it->size());
620  if (eid.isVMEid())
621  {
622  _cOccupancy_FEDVME.fill(eid);
623  _cOccupancy_ElectronicsVME.fill(eid);
624  }
625  else
626  {
627  _cOccupancy_FEDuTCA.fill(eid);
628  _cOccupancy_ElectronicsuTCA.fill(eid);
629  /*
630  if (!it->validate(0, it->size()))
631  {
632  _cCapIdRots_depth.fill(did);
633  _cCapIdRots_FEDuTCA.fill(eid, 1);
634  }*/
635  }
636  }
637 
638  for (int i=0; i<it->size(); i++)
639  {
640  _cADC_SubdetPM.fill(did, it->sample(i).adc());
641  _cfC_SubdetPM.fill(did, it->sample(i).nominal_fC());
642  if (_ptype != fOffline) { // hidefed2crate
643  _cADCvsTS_SubdetPM.fill(did, i, it->sample(i).nominal_fC());
644  if (sumQ>_cutSumQ_HBHE) {
645  _cShapeCut_FED.fill(eid, i, it->sample(i).nominal_fC());
646  }
647  }
648  }
649 
650  if (sumQ>_cutSumQ_HBHE)
651  {
652  //double timing = hcaldqm::utilities::aveTS<HBHEDataFrame>(*it, 2.5, 0, it->size()-1);
653  double timing = hcaldqm::utilities::aveTSDB<HBHEDataFrame>(_dbService, digi_fC, did, *it, 0, it->size()-1);
654  _cTimingCut_SubdetPM.fill(did, timing);
655  _cTimingCut_depth.fill(did, timing);
656  _cOccupancyCut_depth.fill(did);
657  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
658  if (_ptype != fOffline) { // hidefed2crate
659  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
660  }
661  _cSumQ_depth.fill(did, sumQ);
662  _cSumQvsLS_SubdetPM.fill(did, _currentLS, sumQ);
663  if (_ptype==fOnline)
664  {
665  _cSumQvsBX_SubdetPM.fill(did, bx, sumQ);
666  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
667  _cTimingCutvsieta_Subdet.fill(did, timing);
668  _cOccupancyCutvsiphi_SubdetPM.fill(did);
669  _cOccupancyCutvsieta_Subdet.fill(did);
670  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
671  }
672  if (_ptype != fOffline) { // hidefed2crate
673  if (eid.isVMEid())
674  {
675  _cTimingCut_FEDVME.fill(eid, timing);
676  _cTimingCut_ElectronicsVME.fill(eid, timing);
677  _cOccupancyCut_FEDVME.fill(eid);
678  _cOccupancyCut_ElectronicsVME.fill(eid);
679  }
680  else
681  {
682  _cTimingCut_FEDuTCA.fill(eid, timing);
683  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
684  _cOccupancyCut_FEDuTCA.fill(eid);
685  _cOccupancyCut_ElectronicsuTCA.fill(eid);
686  }
687  }
688  did.subdet()==HcalBarrel?numChsCut++:numChsCutHE++;
689  }
690  did.subdet()==HcalBarrel?numChs++:numChsHE++;
691  }
692 
693  // HEP17 collection
694  // The following are filled w.r.t. HBHE digis
695  // - All eta-phi maps
696  // - Occupancy in electronics coordinates
697  // - Digi size
698  //
699  // The following are not filled:
700  // - ADC, fC, sumQ, timing. These are different for QIE11 vs. QIE8. Find them in QIE11Task instead.
701  for (QIE11DigiCollection::const_iterator it=chep17->begin(); it!=chep17->end();
702  ++it)
703  {
704  const QIE11DataFrame digi = static_cast<const QIE11DataFrame>(*it);
705 
706  // Explicit check on the DetIds present in the Collection
707  HcalDetId const& did = digi.detid();
708  uint32_t rawid = _ehashmap.lookup(did);
709  if (rawid==0)
710  {meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
711  HcalElectronicsId const& eid(rawid);
712  if (did.subdet()==HcalBarrel) // Note: since this is HEP17, we obviously expect did.subdet() always to be HcalEndcap, but QIE11DigiCollection may someday expand.
713  rawidHBValid = did.rawId();
714  else if (did.subdet()==HcalEndcap)
715  rawidHEValid = did.rawId();
716 
717  //double sumQ = hcaldqm::utilities::sumQ_v10<QIE11DataFrame>(digi, 2.5, 0, digi.samples()-1);
718  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE11DataFrame>(_dbService, did, digi);
719  double sumQ = hcaldqm::utilities::sumQDB<QIE11DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples()-1);
720 
721  // filter out channels that are masked out
722  if (_xQuality.exists(did))
723  {
724  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
725  if (
726  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
727  cs.isBitSet(HcalChannelStatus::HcalCellDead))
728  continue;
729  }
730 
731  _cOccupancy_depth.fill(did);
732  if (_ptype==fOnline)
733  {
734  _cDigiSizevsLS_FED.fill(eid, _currentLS, digi.samples());
735  digi.samples()!=constants::DIGISIZE[did.subdet()-1]?
736  _xDigiSize.get(eid)++:_xDigiSize.get(eid)+=0;
737  _cOccupancyvsiphi_SubdetPM.fill(did);
738  _cOccupancyvsieta_Subdet.fill(did);
739  }
740  _cDigiSize_Crate.fill(eid, digi.samples());
741  if (_ptype != fOffline) { // hidefed2crate
742  _cDigiSize_FED.fill(eid, digi.samples());
743  if (eid.isVMEid())
744  {
745  _cOccupancy_FEDVME.fill(eid);
746  _cOccupancy_ElectronicsVME.fill(eid);
747  }
748  else
749  {
750  _cOccupancy_FEDuTCA.fill(eid);
751  _cOccupancy_ElectronicsuTCA.fill(eid);
752  /*
753  if (!digi.validate(0, digi.size()))
754  {
755  _cCapIdRots_depth.fill(did);
756  _cCapIdRots_FEDuTCA.fill(eid, 1);
757  }*/
758  }
759  }
760 
761  if (sumQ>_cutSumQ_HEP17)
762  {
763  //double timing = hcaldqm::utilities::aveTS_v10<QIE11DataFrame>(digi, 2.5, 0,digi.samples()-1);
764  double timing = hcaldqm::utilities::aveTSDB<QIE11DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples()-1);
765 
766  _cOccupancyCut_depth.fill(did);
767  _cTimingCut_depth.fill(did, timing);
768  _cSumQ_depth.fill(did, sumQ);
769  if (_ptype==fOnline)
770  {
771  _cOccupancyCutvsiphi_SubdetPM.fill(did);
772  _cOccupancyCutvsieta_Subdet.fill(did);
773  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
774  }
775  if (_ptype != fOffline) { // hidefed2crate
776  if (eid.isVMEid())
777  {
778  _cOccupancyCut_FEDVME.fill(eid);
779  _cOccupancyCut_ElectronicsVME.fill(eid);
780  }
781  else
782  {
783  _cOccupancyCut_FEDuTCA.fill(eid);
784  _cOccupancyCut_ElectronicsuTCA.fill(eid);
785  }
786  }
787  did.subdet()==HcalBarrel?numChsCut++:numChsCutHE++;
788  }
789  did.subdet()==HcalBarrel?numChs++:numChsHE++;
790  }
791 
792  if (rawidHBValid!=0 && rawidHEValid!=0)
793  {
794  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHBValid), _currentLS,
795  numChs);
796  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHEValid), _currentLS,
797  numChsHE);
798  // ONLINE ONLY!
799  if (_ptype==fOnline)
800  {
801  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHBValid),
802  _currentLS, numChsCut);
803  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHBValid), bx,
804  numChsCut);
805  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHEValid),
806  _currentLS, numChsCutHE);
807  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHEValid), bx,
808  numChsCutHE);
809  }
810  // ^^^ONLINE ONLY!
811  }
812  numChs=0;
813  numChsCut = 0;
814 
815  // reset
816  rawidValid = 0;
817 
818 
819  // HO collection
820  for (HODigiCollection::const_iterator it=cho->begin(); it!=cho->end();
821  ++it)
822  {
823  // Explicit check on the DetIds present in the Collection
824  HcalDetId const& did = it->id();
825  uint32_t rawid = _ehashmap.lookup(did);
826  if (rawid==0)
827  {meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
828  HcalElectronicsId const& eid(rawid);
829  if (did.subdet()==HcalOuter)
830  rawidValid = did.rawId();
831 
832  //double sumQ = hcaldqm::utilities::sumQ<HODataFrame>(*it, 8.5, 0, it->size()-1);
833  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<HODataFrame>(_dbService, did, *it);
834  double sumQ = hcaldqm::utilities::sumQDB<HODataFrame>(_dbService, digi_fC, did, *it, 0, it->size()-1);
835 
836  // filter out channels that are masked out
837  if (_xQuality.exists(did))
838  {
839  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
840  if (
841  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
842  cs.isBitSet(HcalChannelStatus::HcalCellDead))
843  continue;
844  }
845 
846  _cSumQ_SubdetPM.fill(did, sumQ);
847  _cOccupancy_depth.fill(did);
848  if (_ptype==fOnline)
849  {
850  _cDigiSizevsLS_FED.fill(eid, _currentLS, it->size());
851  it->size()!=constants::DIGISIZE[did.subdet()-1]?
852  _xDigiSize.get(eid)++:_xDigiSize.get(eid)+=0;
853  _cOccupancyvsiphi_SubdetPM.fill(did);
854  _cOccupancyvsieta_Subdet.fill(did);
855  }
856  _cDigiSize_Crate.fill(eid, it->size());
857  if (_ptype != fOffline) { // hidefed2crate
858  _cDigiSize_FED.fill(eid, it->size());
859  if (eid.isVMEid())
860  {
861  _cOccupancy_FEDVME.fill(eid);
862  _cOccupancy_ElectronicsVME.fill(eid);
863  /*
864  if (!it->validate(0, it->size()))
865  _cCapIdRots_FEDVME.fill(eid, 1);
866  */
867  }
868  else
869  {
870  _cOccupancy_FEDuTCA.fill(eid);
871  _cOccupancy_ElectronicsuTCA.fill(eid);
872  /*
873  if (!it->validate(0, it->size()))
874  _cCapIdRots_FEDuTCA.fill(eid, 1);*/
875  }
876  }
877 
878  for (int i=0; i<it->size(); i++)
879  {
880  _cADC_SubdetPM.fill(did, it->sample(i).adc());
881  _cfC_SubdetPM.fill(did, it->sample(i).nominal_fC());
882  if (_ptype != fOffline) { // hidefed2crate
883  _cADCvsTS_SubdetPM.fill(did, i, it->sample(i).nominal_fC());
884  if (sumQ>_cutSumQ_HO)
885  _cShapeCut_FED.fill(eid, i, it->sample(i).nominal_fC());
886  }
887  }
888 
889  if (sumQ>_cutSumQ_HO)
890  {
891  //double timing = hcaldqm::utilities::aveTS<HODataFrame>(*it, 8.5, 0,it->size()-1);
892  double timing = hcaldqm::utilities::aveTSDB<HODataFrame>(_dbService, digi_fC, did, *it, 0, it->size()-1);
893  _cSumQ_depth.fill(did, sumQ);
894  _cSumQvsLS_SubdetPM.fill(did, _currentLS, sumQ);
895  _cOccupancyCut_depth.fill(did);
896  _cTimingCut_SubdetPM.fill(did, timing);
897  _cTimingCut_depth.fill(did, timing);
898  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
899  if (_ptype != fOffline) { // hidefed2crate
900  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
901  }
902  if (_ptype==fOnline)
903  {
904  _cSumQvsBX_SubdetPM.fill(did, bx, sumQ);
905  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
906  _cTimingCutvsieta_Subdet.fill(did, timing);
907  _cOccupancyCutvsiphi_SubdetPM.fill(did);
908  _cOccupancyCutvsieta_Subdet.fill(did);
909  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
910  }
911  if (_ptype != fOffline) { // hidefed2crate
912  if (eid.isVMEid())
913  {
914  _cTimingCut_FEDVME.fill(eid, timing);
915  _cTimingCut_ElectronicsVME.fill(eid, timing);
916  _cOccupancyCut_FEDVME.fill(eid);
917  _cOccupancyCut_ElectronicsVME.fill(eid);
918  }
919  else
920  {
921  _cTimingCut_FEDuTCA.fill(eid, timing);
922  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
923  _cOccupancyCut_FEDuTCA.fill(eid);
924  _cOccupancyCut_ElectronicsuTCA.fill(eid);
925  }
926  }
927  numChsCut++;
928  }
929  numChs++;
930  }
931 
932  if (rawidValid!=0)
933  {
934  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS,
935  numChs);
936 
937  if (_ptype==fOnline)
938  {
939  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid),
940  _currentLS, numChsCut);
941  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid), bx,
942  numChsCut);
943  }
944  }
945  numChs=0; numChsCut=0;
946 
947  // reset
948  rawidValid = 0;
949 
950  // HF collection
951  if (_qie10InConditions) {
952  for (QIE10DigiCollection::const_iterator it=chf->begin(); it!=chf->end(); ++it) {
953  const QIE10DataFrame digi = static_cast<const QIE10DataFrame>(*it);
954 
955  // Explicit check on the DetIds present in the Collection
956  HcalDetId const& did = digi.detid();
957 
958  // Require subdet == HF. In 2017, calibration channels are included in QIE10DigiCollection :(
959  if (did.subdet() != HcalForward) {
960  continue;
961  }
962 
963  uint32_t rawid = _ehashmap.lookup(did);
964  if (rawid==0) {
965  meUnknownIds1LS->Fill(1);
966  _unknownIdsPresent=true;
967  continue;
968  }
969  HcalElectronicsId const& eid(rawid);
970  if (did.subdet()==HcalForward)
971  rawidValid = did.rawId();
972 
973  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE10DataFrame>(_dbService, did, digi);
974  double sumQ = hcaldqm::utilities::sumQDB<QIE10DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples()-1);
975  //double sumQ = hcaldqm::utilities::sumQ_v10<QIE10DataFrame>(digi, 2.5, 0, digi.samples()-1);
976 
977 
978  // filter out channels that are masked out
979  if (_xQuality.exists(did))
980  {
981  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
982  if (
983  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
984  cs.isBitSet(HcalChannelStatus::HcalCellDead))
985  continue;
986  }
987  if (!_filter_HF.filter(did)) {
988  _cSumQ_SubdetPM_HF.fill(did, sumQ);
989  }
990  _cOccupancy_depth.fill(did);
991  if (_ptype==fOnline)
992  {
993  _xNChs.get(eid)++;
994  _cDigiSizevsLS_FED.fill(eid, _currentLS, digi.samples());
995  digi.samples()!=constants::DIGISIZE[did.subdet()-1]?
996  _xDigiSize.get(eid)++:_xDigiSize.get(eid)+=0;
997  _cOccupancyvsiphi_SubdetPM.fill(did);
998  _cOccupancyvsieta_Subdet.fill(did);
999  }
1000  _cDigiSize_Crate.fill(eid, digi.samples());
1001  if (_ptype != fOffline) { // hidefed2crate
1002  _cDigiSize_FED.fill(eid, digi.samples());
1003  if (eid.isVMEid())
1004  {
1005  _cOccupancy_FEDVME.fill(eid);
1006  _cOccupancy_ElectronicsVME.fill(eid);
1007  /*
1008  if (!it->validate(0, it->size()))
1009  _cCapIdRots_FEDVME.fill(eid, 1);*/
1010  }
1011  else
1012  {
1013  _cOccupancy_FEDuTCA.fill(eid);
1014  _cOccupancy_ElectronicsuTCA.fill(eid);
1015  /*
1016  if (!it->validate(0, it->size()))
1017  _cCapIdRots_FEDuTCA.fill(eid, 1);*/
1018  }
1019  }
1020 
1021  for (int i=0; i<digi.samples(); i++)
1022  {
1023  double q = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, i);
1024  if (!_filter_HF.filter(did)) {
1025  _cADC_SubdetPM_HF.fill(did, digi[i].adc());
1026  _cfC_SubdetPM_HF.fill(did, q);
1027  _cLETDCvsADC_SubdetPM.fill(did, digi[i].adc(), digi[i].le_tdc());
1028  _cLETDCvsTS_SubdetPM.fill(did, (int)i, digi[i].le_tdc());
1029  if (digi[i].le_tdc() <50) {
1030  double time = i*25. + (digi[i].le_tdc() / 2.);
1031  _cLETDCTime_SubdetPM.fill(did, time);
1032  _cLETDCTimevsADC_SubdetPM.fill(did, digi[i].adc(), time);
1033  }
1034  }
1035 
1036  if (_ptype != fOffline) { // hidefed2crate
1037  _cADCvsTS_SubdetPM_HF.fill(did, (int)i, q);
1038  if (sumQ>_cutSumQ_HF)
1039  _cShapeCut_FED.fill(eid, (int)i, q);
1040  }
1041  }
1042 
1043  if (sumQ>_cutSumQ_HF)
1044  {
1045  double timing = hcaldqm::utilities::aveTS_v10<QIE10DataFrame>(digi, 2.5, 0,
1046  digi.samples()-1);
1047  double q1 = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, 1);
1048  double q2 = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, 2);
1049  double q2q12 = q2/(q1+q2);
1050  _cSumQ_depth.fill(did, sumQ);
1051  if (!_filter_HF.filter(did)) {
1052  _cSumQvsLS_SubdetPM_HF.fill(did, _currentLS, sumQ);
1053  }
1054  if (_ptype==fOnline)
1055  {
1056  if (!_filter_HF.filter(did)) {
1057  _cSumQvsBX_SubdetPM_HF.fill(did, bx, sumQ);
1058  }
1059  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
1060  _cTimingCutvsieta_Subdet.fill(did, timing);
1061  _cOccupancyCutvsiphi_SubdetPM.fill(did);
1062  _cOccupancyCutvsieta_Subdet.fill(did);
1063  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
1064  // _cOccupancyCutvsSlotvsLS_HFPM.fill(did, _currentLS);
1065  _xUniHF.get(eid)++;
1066  }
1067  _cTimingCut_SubdetPM.fill(did, timing);
1068  _cTimingCut_depth.fill(did, timing);
1069  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
1070  if (_ptype != fOffline) { // hidefed2crate
1071  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
1072  }
1073  _cOccupancyCut_depth.fill(did);
1074  if (!eid.isVMEid())
1075  if (_ptype==fOnline)
1076  _cQ2Q12CutvsLS_FEDHF.fill(eid, _currentLS, q2q12);
1077  if (_ptype != fOffline) { // hidefed2crate
1078  if (eid.isVMEid())
1079  {
1080  _cTimingCut_FEDVME.fill(eid, timing);
1081  _cTimingCut_ElectronicsVME.fill(eid, timing);
1082  _cOccupancyCut_FEDVME.fill(eid);
1083  _cOccupancyCut_ElectronicsVME.fill(eid);
1084  }
1085  else
1086  {
1087  _cTimingCut_FEDuTCA.fill(eid, timing);
1088  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
1089  _cOccupancyCut_FEDuTCA.fill(eid);
1090  _cOccupancyCut_ElectronicsuTCA.fill(eid);
1091  }
1092  }
1093  numChsCut++;
1094  }
1095  numChs++;
1096  }
1097  }
1098 
1099  if (rawidValid!=0)
1100  {
1101  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS,
1102  numChs);
1103 
1104  if (_ptype==fOnline)
1105  {
1106  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid),
1107  _currentLS, numChsCut);
1108  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid), bx,
1109  numChsCut);
1110  }
1111  }
1112 }
1113 
1114 /* virtual */ void DigiTask::beginLuminosityBlock(
1115  edm::LuminosityBlock const& lb, edm::EventSetup const& es)
1116 {
1117  DQTask::beginLuminosityBlock(lb, es);
1118 }
1119 
1120 /* virtual */ void DigiTask::endLuminosityBlock(edm::LuminosityBlock const& lb,
1121  edm::EventSetup const& es)
1122 {
1123  if (_ptype!=fOnline)
1124  return;
1125 
1126  for (uintCompactMap::const_iterator it=_xUniHF.begin();
1127  it!=_xUniHF.end(); ++it)
1128  {
1129  uint32_t hash1 = it->first;
1130  HcalElectronicsId eid1(hash1);
1131  double x1 = it->second;
1132 
1133  for (uintCompactMap::const_iterator jt=_xUniHF.begin();
1134  jt!=_xUniHF.end(); ++jt)
1135  {
1136  if (jt==it)
1137  continue;
1138  double x2 = jt->second;
1139  if (x2==0)
1140  continue;
1141  if (x1/x2<_thresh_unihf)
1142  _xUni.get(eid1)++;
1143  }
1144  }
1145 
1146  if (_ptype != fOffline) { // hidefed2crate
1147  for (std::vector<uint32_t>::const_iterator it=_vhashFEDs.begin();
1148  it!=_vhashFEDs.end(); ++it)
1149  {
1150  hcaldqm::flag::Flag fSum("DIGI");
1151  HcalElectronicsId eid = HcalElectronicsId(*it);
1152 
1153  std::vector<uint32_t>::const_iterator cit=std::find(
1154  _vcdaqEids.begin(), _vcdaqEids.end(), *it);
1155  if (cit==_vcdaqEids.end())
1156  {
1157  // not @cDAQ
1158  for (uint32_t iflag=0; iflag<_vflags.size(); iflag++)
1159  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag),
1160  int(hcaldqm::flag::fNCDAQ));
1161  _cSummaryvsLS.setBinContent(eid, _currentLS, int(hcaldqm::flag::fNCDAQ));
1162  continue;
1163  }
1164 
1165  // FED is @cDAQ
1166  if (hcaldqm::utilities::isFEDHBHE(eid) || hcaldqm::utilities::isFEDHF(eid) ||
1167  hcaldqm::utilities::isFEDHO(eid))
1168  {
1169  if (_xDigiSize.get(eid)>0)
1170  _vflags[fDigiSize]._state = hcaldqm::flag::fBAD;
1171  else
1172  _vflags[fDigiSize]._state = hcaldqm::flag::fGOOD;
1173  if (hcaldqm::utilities::isFEDHF(eid))
1174  {
1175  double fr = double(_xNChs.get(eid))/double(
1176  _xNChsNominal.get(eid)*_evsPerLS);
1177  if (_runkeyVal==0 || _runkeyVal==4)
1178  {
1179  // only for pp or hi
1180  if (_xUni.get(eid)>0)
1181  _vflags[fUni]._state = hcaldqm::flag::fPROBLEMATIC;
1182  else
1183  _vflags[fUni]._state = hcaldqm::flag::fGOOD;
1184  }
1185  if (fr<0.95)
1186  _vflags[fNChsHF]._state = hcaldqm::flag::fBAD;
1187  else if (fr<1.0)
1188  _vflags[fNChsHF]._state = hcaldqm::flag::fPROBLEMATIC;
1189  else
1190  _vflags[fNChsHF]._state = hcaldqm::flag::fGOOD;
1191  }
1192  }
1193  if (_unknownIdsPresent)
1194  _vflags[fUnknownIds]._state = hcaldqm::flag::fBAD;
1195  else
1196  _vflags[fUnknownIds]._state = hcaldqm::flag::fGOOD;
1197 
1198  int iflag=0;
1199  for (std::vector<hcaldqm::flag::Flag>::iterator ft=_vflags.begin();
1200  ft!=_vflags.end(); ++ft)
1201  {
1202  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, iflag,
1203  int(ft->_state));
1204  fSum+=(*ft);
1205  iflag++;
1206 
1207  // reset!
1208  ft->reset();
1209  }
1210  _cSummaryvsLS.setBinContent(eid, _currentLS, fSum._state);
1211  }
1212  }
1213 
1214  _xDigiSize.reset(); _xUniHF.reset(); _xUni.reset();
1215  _xNChs.reset();
1216 
1217  // in the end always do the DQTask::endLumi
1218  DQTask::endLuminosityBlock(lb, es);
1219 }
1220 
1221 DEFINE_FWK_MODULE(DigiTask);
1222 
ecalMGPA::adc
int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:11
QIE10DataFrame::samples
int samples() const
total number of samples in the digi
Definition: QIE10DataFrame.h:66
QIE11DataFrame::detid
DetId detid() const
Get the detector id.
Definition: QIE11DataFrame.h:42
DigiTask::_cTimingCutvsiphi_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsiphi_SubdetPM
Definition: DigiTask.h:116
DigiTask::_cADC_SubdetPM_HF
hcaldqm::Container1D _cADC_SubdetPM_HF
Definition: DigiTask.h:94
DigiTask::_qie10InConditions
bool _qie10InConditions
Definition: DigiTask.h:175
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
DigiTask::meNumEvents1LS
MonitorElement * meNumEvents1LS
Definition: DigiTask.h:168
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:86
DigiTask::_cOccupancyCutvsieta_Subdet
hcaldqm::Container1D _cOccupancyCutvsieta_Subdet
Definition: DigiTask.h:143
fwrapper::cs
auto_ptr< ClusterSequence > cs
Definition: fastjetfortran_madfks.cc:45
DigiTask::_cSumQvsBX_SubdetPM
hcaldqm::ContainerProf1D _cSumQvsBX_SubdetPM
Definition: DigiTask.h:91
hcaldqm::ContainerXXX::exists
virtual bool exists(HcalDetId const &)
Definition: ContainerXXX.h:267
mps_fire.i
i
Definition: mps_fire.py:156
hcaldqm::constants::CRATE_VME_MIN
int const CRATE_VME_MIN
Definition: Constants.h:111
edm::DataFrameContainer::const_iterator
boost::transform_iterator< IterHelp, boost::counting_iterator< int > > const_iterator
Definition: DataFrameContainer.h:58
HcalDetId::subdet
HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:49
DigiTask::_tagHO
edm::InputTag _tagHO
Definition: DigiTask.h:48
DigiTask::_cADCvsTS_SubdetPM
hcaldqm::Container2D _cADCvsTS_SubdetPM
Definition: DigiTask.h:103
hcaldqm::utilities::fed2crate
std::pair< uint16_t, uint16_t > fed2crate(int fed)
Definition: Utilities.cc:12
HcalElectronicsId::rawId
uint32_t rawId() const
Definition: HcalElectronicsId.h:47
DigiTask::_cfC_SubdetPM_HF
hcaldqm::Container1D _cfC_SubdetPM_HF
Definition: DigiTask.h:95
DigiTask::_cQ2Q12CutvsLS_FEDHF
hcaldqm::ContainerProf1D _cQ2Q12CutvsLS_FEDHF
Definition: DigiTask.h:120
DigiTask::meUnknownIds1LS
MonitorElement * meUnknownIds1LS
Definition: DigiTask.h:169
DigiTask::_tokHBHE
edm::EDGetTokenT< HBHEDigiCollection > _tokHBHE
Definition: DigiTask.h:50
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:460
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::_cLETDCTimevsADC_SubdetPM
hcaldqm::Container2D _cLETDCTimevsADC_SubdetPM
Definition: DigiTask.h:162
DigiTask::_cutSumQ_HO
double _cutSumQ_HO
Definition: DigiTask.h:55
hcaldqm::ContainerXXX::begin
virtual CompactMap::const_iterator begin()
Definition: ContainerXXX.h:75
DigiTask::_thresh_unihf
double _thresh_unihf
Definition: DigiTask.h:56
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:66
hcaldqm::constants::SPIGOT_MIN
int const SPIGOT_MIN
Definition: Constants.h:136
DigiTask::_cSumQ_depth
hcaldqm::ContainerProf2D _cSumQ_depth
Definition: DigiTask.h:89
DigiTask::_tokHO
edm::EDGetTokenT< HODigiCollection > _tokHO
Definition: DigiTask.h:52
DigiTask::_cOccupancyCutvsiphi_SubdetPM
hcaldqm::Container1D _cOccupancyCutvsiphi_SubdetPM
Definition: DigiTask.h:142
edm::DataFrameContainer::begin
const_iterator begin() const
Definition: DataFrameContainer.h:159
hcaldqm::ContainerXXX::reset
virtual void reset()
Definition: ContainerXXX.h:366
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:20
DigiTask::_cSummaryvsLS
hcaldqm::ContainerSingle2D _cSummaryvsLS
Definition: DigiTask.h:173
DigiTask::_cTimingCut_FEDVME
hcaldqm::ContainerProf2D _cTimingCut_FEDVME
Definition: DigiTask.h:109
hcaldqm::constants::FIBER_VME_MIN
int const FIBER_VME_MIN
Definition: Constants.h:141
DigiTask::_cOccupancyCutvsiphivsLS_SubdetPM
hcaldqm::Container2D _cOccupancyCutvsiphivsLS_SubdetPM
Definition: DigiTask.h:145
DigiTask::_xDigiSize
hcaldqm::ContainerXXX< uint32_t > _xDigiSize
Definition: DigiTask.h:156
DigiTask::_cTimingCut_FEDuTCA
hcaldqm::ContainerProf2D _cTimingCut_FEDuTCA
Definition: DigiTask.h:110
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:151
DigiTask::beginLuminosityBlock
virtual void beginLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &)
Definition: DigiTask.cc:1114
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:119
DigiTask::_filter_VME
hcaldqm::filter::HashFilter _filter_VME
Definition: DigiTask.h:77
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
hcaldqm::ContainerXXX::end
virtual CompactMap::const_iterator end()
Definition: ContainerXXX.h:77
DigiTask::_cLETDCvsTS_SubdetPM
hcaldqm::Container2D _cLETDCvsTS_SubdetPM
Definition: DigiTask.h:164
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
DigiTask::_cOccupancy_CrateSlot
hcaldqm::Container2D _cOccupancy_CrateSlot
Definition: DigiTask.h:129
hcaldqm::ContainerSingle2D::setBinContent
virtual void setBinContent(int, int, int)
Definition: ContainerSingle2D.cc:209
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:98
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:139
DigiTask::_cDigiSize_Crate
hcaldqm::Container1D _cDigiSize_Crate
Definition: DigiTask.h:153
DigiTask::_cShapeCut_FED
hcaldqm::Container1D _cShapeCut_FED
Definition: DigiTask.h:102
DigiTask::_tagHEP17
edm::InputTag _tagHEP17
Definition: DigiTask.h:47
DigiTask::_cOccupancyCut_depth
hcaldqm::Container2D _cOccupancyCut_depth
Definition: DigiTask.h:141
DigiTask::_cOccupancyCutvsBX_Subdet
hcaldqm::ContainerProf1D _cOccupancyCutvsBX_Subdet
Definition: DigiTask.h:150
DigiTask::_cutSumQ_HF
double _cutSumQ_HF
Definition: DigiTask.h:55
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:144
DigiTask::_resetMonitors
virtual void _resetMonitors(hcaldqm::UpdateFreq)
Definition: DigiTask.cc:519
DigiTask::_cutSumQ_HBHE
double _cutSumQ_HBHE
Definition: DigiTask.h:55
hcaldqm::hashfunctions::hash_did
hash_function_did const hash_did[nHashType_did]
Definition: HashFunctions.h:167
DigiTask::_cTimingCut_depth
hcaldqm::ContainerProf2D _cTimingCut_depth
Definition: DigiTask.h:115
QIE10DataFrame::detid
DetId detid() const
Get the detector id.
Definition: QIE10DataFrame.h:56
HcalElectronicsMap::allPrecisionId
std::vector< HcalGenericDetId > allPrecisionId() const
Definition: HcalElectronicsMap.cc:139
hcaldqm::DQModule::_logger
Logger _logger
Definition: DQModule.h:70
DigiTask::_cOccupancy_FEDVME
hcaldqm::Container2D _cOccupancy_FEDVME
Definition: DigiTask.h:124
DigiTask::_cTimingCutvsLS_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsLS_SubdetPM
Definition: DigiTask.h:114
DigiTask::_cOccupancyCut_FEDVME
hcaldqm::Container2D _cOccupancyCut_FEDVME
Definition: DigiTask.h:137
DigiTask::_cOccupancy_depth
hcaldqm::Container2D _cOccupancy_depth
Definition: DigiTask.h:130
hcaldqm::DQModule::_name
std::string _name
Definition: DQModule.h:57
hcaldqm::DQTask::_emap
HcalElectronicsMap const * _emap
Definition: DQTask.h:73
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:172
DigiTask::_tokHEP17
edm::EDGetTokenT< QIE11DigiCollection > _tokHEP17
Definition: DigiTask.h:51
DigiTask::_filter_FEDHF
hcaldqm::filter::HashFilter _filter_FEDHF
Definition: DigiTask.h:79
hcaldqm::constants::CRATE_uTCA_MIN
int const CRATE_uTCA_MIN
Definition: Constants.h:116
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:310
DigiTask::_tokHF
edm::EDGetTokenT< QIE10DigiCollection > _tokHF
Definition: DigiTask.h:53
DigiTask::DigiTask
DigiTask(edm::ParameterSet const &)
Definition: DigiTask.cc:7
DigiTask::_cOccupancyvsLS_Subdet
hcaldqm::ContainerProf1D _cOccupancyvsLS_Subdet
Definition: DigiTask.h:148
DigiTask::_cSumQvsLS_SubdetPM
hcaldqm::ContainerProf1D _cSumQvsLS_SubdetPM
Definition: DigiTask.h:90
DigiTask::_xUni
hcaldqm::ContainerXXX< uint32_t > _xUni
Definition: DigiTask.h:157
DigiTask::_cTimingCutvsLS_FED
hcaldqm::ContainerProf1D _cTimingCutvsLS_FED
Definition: DigiTask.h:113
HcalElectronicsId::isVMEid
bool isVMEid() const
Definition: HcalElectronicsId.h:49
DigiTask::_cOccupancyCut_FEDuTCA
hcaldqm::Container2D _cOccupancyCut_FEDuTCA
Definition: DigiTask.h:138
DigiTask::_cOccupancy_ElectronicsVME
hcaldqm::Container2D _cOccupancy_ElectronicsVME
Definition: DigiTask.h:126
DigiTask::_tagHF
edm::InputTag _tagHF
Definition: DigiTask.h:49
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:277
hcaldqm::utilities::isFEDHBHE
bool isFEDHBHE(HcalElectronicsId const &)
Definition: Utilities.cc:172
hcaldqm::ContainerXXX::book
virtual void book(HcalElectronicsMap const *)
Definition: ContainerXXX.h:99
CaloSamples::size
int size() const
get the size
Definition: CaloSamples.h:24
DigiTask::_cDigiSize_FED
hcaldqm::Container1D _cDigiSize_FED
Definition: DigiTask.h:154
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:159
DigiTask::_cSumQ_SubdetPM
hcaldqm::Container1D _cSumQ_SubdetPM
Definition: DigiTask.h:88
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
DigiTask::_cOccupancyvsiphi_SubdetPM
hcaldqm::Container1D _cOccupancyvsiphi_SubdetPM
Definition: DigiTask.h:131
hcaldqm::constants::DIGISIZE
int const DIGISIZE[SUBDET_NUM]
Definition: Constants.h:183
DigiTask::_cTimingCutvsieta_Subdet
hcaldqm::ContainerProf1D _cTimingCutvsieta_Subdet
Definition: DigiTask.h:117
DigiTask::_filter_HEP17
hcaldqm::filter::HashFilter _filter_HEP17
Definition: DigiTask.h:82
DigiTask::_cLETDCTime_SubdetPM
hcaldqm::Container1D _cLETDCTime_SubdetPM
Definition: DigiTask.h:165
DigiTask::_dhashmap
hcaldqm::electronicsmap::ElectronicsMap _dhashmap
Definition: DigiTask.h:74
DigiTask::_vflags
std::vector< hcaldqm::flag::Flag > _vflags
Definition: DigiTask.h:59
hcaldqm::Container2D::setBinContent
virtual void setBinContent(HcalDetId const &, int) override
Definition: Container2D.cc:235
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
edm::DataFrameContainer::end
const_iterator end() const
Definition: DataFrameContainer.h:163
DigiTask::_cOccupancyCut_ElectronicsuTCA
hcaldqm::Container2D _cOccupancyCut_ElectronicsuTCA
Definition: DigiTask.h:140
DigiTask::_cTimingCut_ElectronicsVME
hcaldqm::ContainerProf2D _cTimingCut_ElectronicsVME
Definition: DigiTask.h:111
DigiTask::_cLETDCvsADC_SubdetPM
hcaldqm::Container2D _cLETDCvsADC_SubdetPM
Definition: DigiTask.h:163
hcaldqm::ContainerXXX::get
virtual STDTYPE & get(HcalDetId const &)
Definition: ContainerXXX.h:249
hcaldqm::DQTask::_dbService
edm::ESHandle< HcalDbService > _dbService
Definition: DQTask.h:72
hcaldqm::Container1D::setLumiFlag
virtual void setLumiFlag()
Definition: Container1D.cc:1360
hcaldqm::utilities::getFEDList
std::vector< int > getFEDList(HcalElectronicsMap const *)
Definition: Utilities.cc:95
DigiTask::_cOccupancy_Crate
hcaldqm::Container2D _cOccupancy_Crate
Definition: DigiTask.h:128
hcaldqm::utilities::getFEDuTCAList
std::vector< int > getFEDuTCAList(HcalElectronicsMap const *)
Definition: Utilities.cc:145
hcaldqm::Container2D::fill
virtual void fill(HcalDetId const &) override
Definition: Container2D.cc:59
HcalDbService::getHcalMapping
const HcalElectronicsMap * getHcalMapping() const
Definition: HcalDbService.cc:267
hcaldqm::filter::HashFilter::filter
virtual bool filter(HcalDetId const &) const
Definition: HashFilter.cc:36
hcaldqm::DQModule::_subsystem
std::string _subsystem
Definition: DQModule.h:64
DigiTask::_filter_uTCA
hcaldqm::filter::HashFilter _filter_uTCA
Definition: DigiTask.h:78
DigiTask::endLuminosityBlock
virtual void endLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &)
Definition: DigiTask.cc:1120
DigiTask::_cSumQvsLS_SubdetPM_HF
hcaldqm::ContainerProf1D _cSumQvsLS_SubdetPM_HF
Definition: DigiTask.h:97
hcaldqm::utilities::isFEDHO
bool isFEDHO(HcalElectronicsId const &)
Definition: Utilities.cc:212
hcaldqm::DQTask::_vcdaqEids
std::vector< uint32_t > _vcdaqEids
Definition: DQTask.h:65
DigiTask::_cDigiSizevsLS_FED
hcaldqm::ContainerProf1D _cDigiSizevsLS_FED
Definition: DigiTask.h:155
DigiTask::_xNChs
hcaldqm::ContainerXXX< uint32_t > _xNChs
Definition: DigiTask.h:158
DigiTask::_ehashmap
hcaldqm::electronicsmap::ElectronicsMap _ehashmap
Definition: DigiTask.h:73
DigiTask::_cOccupancyCutvsLS_Subdet
hcaldqm::ContainerProf1D _cOccupancyCutvsLS_Subdet
Definition: DigiTask.h:149
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::Container2D::book
virtual void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="") override
Definition: Container2D.cc:895
hcaldqm::ContainerProf2D::book
virtual void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="")
Definition: ContainerProf2D.cc:51
DigiTask::_cSumQ_SubdetPM_HF
hcaldqm::Container1D _cSumQ_SubdetPM_HF
Definition: DigiTask.h:96
DigiTask::_vhashFEDs
std::vector< uint32_t > _vhashFEDs
Definition: DigiTask.h:70
DigiTask::_filter_HF
hcaldqm::filter::HashFilter _filter_HF
Definition: DigiTask.h:80
QIE11DataFrame::samples
int samples() const
total number of samples in the digi
Definition: QIE11DataFrame.h:52
DigiTask::bookHistograms
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
Definition: DigiTask.cc:39
DigiTask::_cSumQvsBX_SubdetPM_HF
hcaldqm::ContainerProf1D _cSumQvsBX_SubdetPM_HF
Definition: DigiTask.h:98
DigiTask::_process
virtual void _process(edm::Event const &, edm::EventSetup const &)
Definition: DigiTask.cc:539
DigiTask::_cOccupancyvsieta_Subdet
hcaldqm::Container1D _cOccupancyvsieta_Subdet
Definition: DigiTask.h:132
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:187
DigiTask::_cTimingCut_ElectronicsuTCA
hcaldqm::ContainerProf2D _cTimingCut_ElectronicsuTCA
Definition: DigiTask.h:112
DigiTask::_filter_notHF
hcaldqm::filter::HashFilter _filter_notHF
Definition: DigiTask.h:81
hcaldqm::constants::SLOT_uTCA_MIN
int const SLOT_uTCA_MIN
Definition: Constants.h:123
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:303
DigiTask::_cTimingCut_SubdetPM
hcaldqm::Container1D _cTimingCut_SubdetPM
Definition: DigiTask.h:108
edm::Run
Definition: Run.h:42
DigiTask::_unknownIdsPresent
bool _unknownIdsPresent
Definition: DigiTask.h:170
DigiTask::_cOccupancy_FEDuTCA
hcaldqm::Container2D _cOccupancy_FEDuTCA
Definition: DigiTask.h:125
cuy.ib
ib
Definition: cuy.py:660
DigiTask::_xUniHF
hcaldqm::ContainerXXX< uint32_t > _xUniHF
Definition: DigiTask.h:157
DigiTask::_cADCvsTS_SubdetPM_HF
hcaldqm::Container2D _cADCvsTS_SubdetPM_HF
Definition: DigiTask.h:104
DigiTask::_cutSumQ_HEP17
double _cutSumQ_HEP17
Definition: DigiTask.h:55
hcaldqm::utilities::sumQ
double sumQ(DIGI const &digi, double ped, int i=0, int j=3)
Definition: Utilities.h:122
DigiTask::_cfC_SubdetPM
hcaldqm::Container1D _cfC_SubdetPM
Definition: DigiTask.h:87
DigiTask::_cOccupancy_ElectronicsuTCA
hcaldqm::Container2D _cOccupancy_ElectronicsuTCA
Definition: DigiTask.h:127
Generated for CMSSW Reference Manual by doxygen 1.8.11