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  _tagHE = ps.getUntrackedParameter<edm::InputTag>("tagHE",
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  _tokHE = consumes<QIE11DigiCollection>(_tagHE);
21  _tokHO = consumes<HODigiCollection>(_tagHO);
22  _tokHF = consumes<QIE10DigiCollection>(_tagHF);
23 
24  _cutSumQ_HBHE = ps.getUntrackedParameter<double>("cutSumQ_HBHE", 20);
25  _cutSumQ_HE = ps.getUntrackedParameter<double>("cutSumQ_HE", 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  _thresh_led = ps.getUntrackedParameter<double>("thresh_led", 20);
30 
31  _vflags.resize(nDigiFlag);
32  _vflags[fUni]=hcaldqm::flag::Flag("UniSlotHF");
33  _vflags[fDigiSize]=hcaldqm::flag::Flag("DigiSize");
34  _vflags[fNChsHF]=hcaldqm::flag::Flag("NChsHF");
35  _vflags[fUnknownIds]=hcaldqm::flag::Flag("UnknownIds");
36  _vflags[fLED]=hcaldqm::flag::Flag("LEDMisfire");
37  _vflags[fCapId]=hcaldqm::flag::Flag("BadCapId");
38 
39  _qie10InConditions = ps.getUntrackedParameter<bool>("qie10InConditions", true);
40 
41  // Get reference digi sizes. Convert from unsigned to signed int, because <digi>::size()/samples() return ints for some reason.
42  std::vector<uint32_t> vrefDigiSize = ps.getUntrackedParameter<std::vector<uint32_t>>("refDigiSize");
43  _refDigiSize[HcalBarrel] = (int)vrefDigiSize[0];
44  _refDigiSize[HcalEndcap] = (int)vrefDigiSize[1];
45  _refDigiSize[HcalOuter] = (int)vrefDigiSize[2];
46  _refDigiSize[HcalForward] = (int)vrefDigiSize[3];
47 
48  // (capid - BX) % 4 to 1
49  _capidmbx[HcalBarrel] = 1;
50  _capidmbx[HcalEndcap] = 1;
51  _capidmbx[HcalOuter] = 1;
52  _capidmbx[HcalForward] = 1;
53 }
54 
55 /* virtual */ void DigiTask::bookHistograms(DQMStore::IBooker& ib,
56  edm::Run const& r, edm::EventSetup const& es)
57 {
58  DQTask::bookHistograms(ib,r,es);
59 
60  // GET WHAT YOU NEED
61  edm::ESHandle<HcalDbService> dbs;
62  es.get<HcalDbRecord>().get(dbs);
63  _emap = dbs->getHcalMapping();
64  std::vector<uint32_t> vVME;
65  std::vector<uint32_t> vuTCA;
66  vVME.push_back(HcalElectronicsId(constants::FIBERCH_MIN,
67  constants::FIBER_VME_MIN, SPIGOT_MIN, CRATE_VME_MIN).rawId());
68  vuTCA.push_back(HcalElectronicsId(CRATE_uTCA_MIN, SLOT_uTCA_MIN,
69  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
70  _filter_VME.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics,
71  vVME);
72  _filter_uTCA.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics,
73  vuTCA);
74 
75  // Filters for QIE8 vs QIE10/11
76  std::vector<uint32_t> vhashQIE1011;
77  vhashQIE1011.push_back(hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdet](HcalDetId(HcalEndcap, 20,1,1)));
78  vhashQIE1011.push_back(hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdet](HcalDetId(HcalForward, 29,1,1)));
79  _filter_QIE1011.initialize(filter::fPreserver, hcaldqm::hashfunctions::fSubdet,
80  vhashQIE1011);
81  _filter_QIE8.initialize(filter::fFilter, hcaldqm::hashfunctions::fSubdet,
82  vhashQIE1011);
83 
84  // INITIALIZE FIRST
85  _cADC_SubdetPM.initialize(_name, "ADC", hcaldqm::hashfunctions::fSubdetPM,
86  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_128),
87  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
88  _cfC_SubdetPM.initialize(_name, "fC", hcaldqm::hashfunctions::fSubdetPM,
89  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),
90  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
91  _cSumQ_SubdetPM.initialize(_name, "SumQ", hcaldqm::hashfunctions::fSubdetPM,
92  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),
93  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
94  _cSumQ_depth.initialize(_name, "SumQ", hcaldqm::hashfunctions::fdepth,
95  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
96  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
97  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_400000, true),0);
98  _cSumQvsLS_SubdetPM.initialize(_name, "SumQvsLS",
99  hcaldqm::hashfunctions::fSubdetPM,
100  new hcaldqm::quantity::LumiSection(_maxLS),
101  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),0);
102 
103  _cADC_SubdetPM_QIE1011.initialize(_name, "ADC", hcaldqm::hashfunctions::fSubdetPM,
104  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
105  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
106  _cfC_SubdetPM_QIE1011.initialize(_name, "fC", hcaldqm::hashfunctions::fSubdetPM,
107  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000),
108  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
109  _cSumQ_SubdetPM_QIE1011.initialize(_name, "SumQ", hcaldqm::hashfunctions::fSubdetPM,
110  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000),
111  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
112  _cSumQvsLS_SubdetPM_QIE1011.initialize(_name, "SumQvsLS",
113  hcaldqm::hashfunctions::fSubdetPM,
114  new hcaldqm::quantity::LumiSection(_maxLS),
115  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000),0);
116 
117  _cTimingCut_SubdetPM.initialize(_name, "TimingCut",
118  hcaldqm::hashfunctions::fSubdetPM,
119  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
120  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
121  _cTimingCut_depth.initialize(_name, "TimingCut",
122  hcaldqm::hashfunctions::fdepth,
123  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
124  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
125  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
126  _cTimingCutvsLS_SubdetPM.initialize(_name, "TimingvsLS",
127  hcaldqm::hashfunctions::fSubdetPM,
128  new hcaldqm::quantity::LumiSection(_maxLS),
129  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
130 
131  // Occupancy w/o a cut
132  _cOccupancyvsLS_Subdet.initialize(_name, "OccupancyvsLS",
133  hcaldqm::hashfunctions::fSubdet,
134  new hcaldqm::quantity::LumiSection(_maxLS),
135  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),0);
136  _cOccupancy_depth.initialize(_name, "Occupancy",
137  hcaldqm::hashfunctions::fdepth,
138  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
139  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
140  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
141 
142  // Occupancy w/ a cut
143  _cOccupancyCutvsLS_Subdet.initialize(_name, "OccupancyCutvsLS",
144  hcaldqm::hashfunctions::fSubdet,
145  new hcaldqm::quantity::LumiSection(_maxLS),
146  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),0);
147  _cOccupancyCut_depth.initialize(_name, "OccupancyCut",
148  hcaldqm::hashfunctions::fdepth,
149  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
150  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
151  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
152 
153  // Digi size
154  _cDigiSize_Crate.initialize(_name, "DigiSize",
155  hcaldqm::hashfunctions::fCrate,
156  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDigiSize),
157  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
158  _cADCvsTS_SubdetPM.initialize(_name, "ADCvsTS",
159  hcaldqm::hashfunctions::fSubdetPM,
160  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
161  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_128),
162  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
163  _cADCvsTS_SubdetPM_QIE1011.initialize(_name, "ADCvsTS",
164  hcaldqm::hashfunctions::fSubdetPM,
165  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
166  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
167  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
168 
169  _cLETDCTimevsADC_SubdetPM.initialize(_name, "LETDCTimevsADC",
170  hcaldqm::hashfunctions::fSubdetPM,
171  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
172  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250_coarse),
173  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
174  _cLETDCvsADC_SubdetPM.initialize(_name, "LETDCvsADC",
175  hcaldqm::hashfunctions::fSubdetPM,
176  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
177  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
178  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
179  _cLETDCvsTS_SubdetPM.initialize(_name, "LETDCvsTS",
180  hcaldqm::hashfunctions::fSubdetPM,
181  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
182  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
183  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
184  _cLETDCTime_SubdetPM.initialize(_name, "LETDCTime",
185  hcaldqm::hashfunctions::fSubdetPM,
186  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250),
187  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
188  _cLETDCTime_depth.initialize(_name, "LETDCTime",
189  hcaldqm::hashfunctions::fdepth,
190  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
191  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
192  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
193 
194  _cBadTDCValues_SubdetPM.initialize(_name, "BadTDCValues",
195  hcaldqm::hashfunctions::fSubdetPM,
196  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBadTDC),
197  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
198  _cBadTDCvsBX_SubdetPM.initialize(_name, "BadTDCvsBX",
199  hcaldqm::hashfunctions::fSubdetPM,
200  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
201  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
202  _cBadTDCvsLS_SubdetPM.initialize(_name, "BadTDCvsLS",
203  hcaldqm::hashfunctions::fSubdetPM,
204  new hcaldqm::quantity::LumiSection(_maxLS),
205  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
206  _cBadTDCCount_depth.initialize(_name, "BadTDCCount",
207  hcaldqm::hashfunctions::fdepth,
208  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
209  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
210  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
211 
212  if (_ptype == fOnline || _ptype == fLocal) {
213  _cOccupancy_Crate.initialize(_name,
214  "Occupancy", hashfunctions::fCrate,
215  new quantity::ElectronicsQuantity(quantity::fSlotuTCA),
216  new quantity::ElectronicsQuantity(quantity::fFiberuTCAFiberCh),
217  new quantity::ValueQuantity(quantity::fN),0);
218  _cOccupancy_CrateSlot.initialize(_name,
219  "Occupancy", hashfunctions::fCrateSlot,
220  new quantity::ElectronicsQuantity(quantity::fFiberuTCA),
221  new quantity::ElectronicsQuantity(quantity::fFiberCh),
222  new quantity::ValueQuantity(quantity::fN),0);
223  }
224 
225  // INITIALIZE HISTOGRAMS that are only for Online
226  if (_ptype==fOnline)
227  {
228  // Charge sharing
229  _cQ2Q12CutvsLS_FEDHF.initialize(_name, "Q2Q12vsLS",
230  hcaldqm::hashfunctions::fFED,
231  new hcaldqm::quantity::LumiSection(_maxLS),
232  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fRatio_0to2),0);
233  _cSumQvsBX_SubdetPM.initialize(_name, "SumQvsBX",
234  hcaldqm::hashfunctions::fSubdetPM,
235  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
236  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),0);
237  _cSumQvsBX_SubdetPM_QIE1011.initialize(_name, "SumQvsBX",
238  hcaldqm::hashfunctions::fSubdetPM,
239  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
240  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_10000),0);
241  _cDigiSizevsLS_FED.initialize(_name, "DigiSizevsLS",
242  hcaldqm::hashfunctions::fFED,
243  new hcaldqm::quantity::LumiSection(_maxLS),
244  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDigiSize),0);
245  _cTimingCutvsiphi_SubdetPM.initialize(_name, "TimingCutvsiphi",
246  hcaldqm::hashfunctions::fSubdetPM,
247  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
248  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
249  _cTimingCutvsieta_Subdet.initialize(_name, "TimingCutvsieta",
250  hcaldqm::hashfunctions::fSubdet,
251  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
252  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
253  _cOccupancyvsiphi_SubdetPM.initialize(_name, "Occupancyvsiphi",
254  hcaldqm::hashfunctions::fSubdetPM,
255  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
256  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
257  _cOccupancyvsieta_Subdet.initialize(_name, "Occupancyvsieta",
258  hcaldqm::hashfunctions::fSubdet,
259  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
260  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
261  _cOccupancyCutvsiphi_SubdetPM.initialize(_name, "OccupancyCutvsiphi",
262  hcaldqm::hashfunctions::fSubdetPM,
263  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
264  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
265  _cOccupancyCutvsieta_Subdet.initialize(_name, "OccupancyCutvsieta",
266  hcaldqm::hashfunctions::fSubdet,
267  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
268  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
269  _cOccupancyCutvsLS_Subdet.initialize(_name, "OccupancyCutvsLS",
270  hcaldqm::hashfunctions::fSubdet,
271  new hcaldqm::quantity::LumiSection(_maxLS),
272  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),0);
273  _cOccupancyCutvsBX_Subdet.initialize(_name, "OccupancyCutvsBX",
274  hcaldqm::hashfunctions::fSubdet,
275  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
276  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),0);
277 // _cOccupancyCutvsSlotvsLS_HFPM.initialize(_name,
278 // "OccupancyCutvsSlotvsLS", hcaldqm::hashfunctions::fSubdetPM,
279 // new hcaldqm::quantity::LumiSection(_maxLS),
280 // new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
281 // new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
282  _cOccupancyCutvsiphivsLS_SubdetPM.initialize(_name,
283  "OccupancyCutvsiphivsLS", hcaldqm::hashfunctions::fSubdetPM,
284  new hcaldqm::quantity::LumiSection(_maxLS),
285  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
286  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
287  }
288  _cCapidMinusBXmod4_SubdetPM.initialize(_name,
289  "CapID", hcaldqm::hashfunctions::fSubdetPM,
290  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fCapidMinusBXmod4),
291  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
292 
293  for (int i = 0; i < 4; ++i) {
294  _cCapidMinusBXmod4_CrateSlotuTCA[i].initialize(_name, "CapID",
295  new quantity::ElectronicsQuantity(quantity::fCrateuTCA),
296  new quantity::ElectronicsQuantity(quantity::fSlotuTCA),
297  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
298  _cCapidMinusBXmod4_CrateSlotVME[i].initialize(_name, "CapID",
299  new quantity::ElectronicsQuantity(quantity::fCrateVME),
300  new quantity::ElectronicsQuantity(quantity::fSlotVME),
301  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
302  }
303 
304  if (_ptype != fOffline) { // hidefed2crate
305  std::vector<int> vFEDs = hcaldqm::utilities::getFEDList(_emap);
306  std::vector<int> vFEDsVME = hcaldqm::utilities::getFEDVMEList(_emap);
307  std::vector<int> vFEDsuTCA = hcaldqm::utilities::getFEDuTCAList(_emap);
308 
309  std::vector<uint32_t> vFEDHF;
310  vFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN,
311  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
312  vFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN+6,
313  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
314  vFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN,
315  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
316  vFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN+6,
317  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
318  vFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN,
319  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
320  vFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN+6,
321  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
322 
323  // initialize filters
324  _filter_FEDHF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fFED,
325  vFEDHF);
326 
327  // push the rawIds of each fed into the vector...
328  for (std::vector<int>::const_iterator it=vFEDsVME.begin();
329  it!=vFEDsVME.end(); ++it)
330  _vhashFEDs.push_back(HcalElectronicsId(
331  constants::FIBERCH_MIN, FIBER_VME_MIN, SPIGOT_MIN,
332  (*it)-FED_VME_MIN).rawId());
333  for (std::vector<int>::const_iterator it=vFEDsuTCA.begin();
334  it!=vFEDsuTCA.end(); ++it)
335  {
336  std::pair<uint16_t, uint16_t> cspair = utilities::fed2crate(*it);
337  _vhashFEDs.push_back(HcalElectronicsId(
338  cspair.first, cspair.second, FIBER_uTCA_MIN1,
339  FIBERCH_MIN, false).rawId());
340  }
341 
342  _cShapeCut_FED.initialize(_name, "ShapeCut",
343  hcaldqm::hashfunctions::fFED,
344  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
345  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_10000),0);
346 
347  _cTimingCut_FEDVME.initialize(_name, "TimingCut",
348  hcaldqm::hashfunctions::fFED,
349  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
350  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
351  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
352  _cTimingCut_FEDuTCA.initialize(_name, "TimingCut",
353  hcaldqm::hashfunctions::fFED,
354  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
355  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
356  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
357  _cTimingCut_ElectronicsVME.initialize(_name, "TimingCut",
358  hcaldqm::hashfunctions::fElectronics,
359  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
360  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
361  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
362  _cTimingCut_ElectronicsuTCA.initialize(_name, "TimingCut",
363  hcaldqm::hashfunctions::fElectronics,
364  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
365  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
366  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
367  _cTimingCutvsLS_FED.initialize(_name, "TimingvsLS",
368  hcaldqm::hashfunctions::fFED,
369  new hcaldqm::quantity::LumiSection(_maxLS),
370  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
371 
372  _cOccupancy_FEDVME.initialize(_name, "Occupancy",
373  hcaldqm::hashfunctions::fFED,
374  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
375  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
376  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
377  _cOccupancy_FEDuTCA.initialize(_name, "Occupancy",
378  hcaldqm::hashfunctions::fFED,
379  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
380  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
381  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
382  _cOccupancy_ElectronicsVME.initialize(_name, "Occupancy",
383  hcaldqm::hashfunctions::fElectronics,
384  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
385  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
386  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
387  _cOccupancy_ElectronicsuTCA.initialize(_name, "Occupancy",
388  hcaldqm::hashfunctions::fElectronics,
389  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
390  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
391  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
392 
393  _cOccupancyCut_FEDVME.initialize(_name, "OccupancyCut",
394  hcaldqm::hashfunctions::fFED,
395  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
396  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
397  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
398  _cOccupancyCut_FEDuTCA.initialize(_name, "OccupancyCut",
399  hcaldqm::hashfunctions::fFED,
400  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
401  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
402  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
403  _cOccupancyCut_ElectronicsVME.initialize(_name, "OccupancyCut",
404  hcaldqm::hashfunctions::fElectronics,
405  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
406  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
407  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
408  _cOccupancyCut_ElectronicsuTCA.initialize(_name, "OccupancyCut",
409  hcaldqm::hashfunctions::fElectronics,
410  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
411  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
412  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
413 
414  _cDigiSize_FED.initialize(_name, "DigiSize",
415  hcaldqm::hashfunctions::fFED,
416  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDigiSize),
417  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
418 
419  if (_ptype == fOnline) {
420  _cSummaryvsLS_FED.initialize(_name, "SummaryvsLS",
421  hcaldqm::hashfunctions::fFED,
422  new hcaldqm::quantity::LumiSection(_maxLS),
423  new hcaldqm::quantity::FlagQuantity(_vflags),
424  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),0);
425  _cSummaryvsLS.initialize(_name, "SummaryvsLS",
426  new hcaldqm::quantity::LumiSection(_maxLS),
427  new hcaldqm::quantity::FEDQuantity(vFEDs),
428  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),0);
429 
430  _xUniHF.initialize(hcaldqm::hashfunctions::fFEDSlot);
431  _xUni.initialize(hcaldqm::hashfunctions::fFED);
432  _xDigiSize.initialize(hcaldqm::hashfunctions::fFED);
433  _xNChs.initialize(hcaldqm::hashfunctions::fFED);
434  _xNChsNominal.initialize(hcaldqm::hashfunctions::fFED);
435  _xBadCapid.initialize(hcaldqm::hashfunctions::fFED);
436  }
437  }
438 
439  // BOOK HISTOGRAMS
440  char cutstr[200];
441  sprintf(cutstr, "_SumQHBHE%dHO%dHF%d", int(_cutSumQ_HBHE),
442  int(_cutSumQ_HO), int(_cutSumQ_HF));
443  char cutstr2[200];
444  sprintf(cutstr2, "_SumQHF%d", int(_cutSumQ_HF));
445 
446  _cADC_SubdetPM.book(ib, _emap, _filter_QIE8, _subsystem);
447  _cADC_SubdetPM_QIE1011.book(ib, _emap, _filter_QIE1011, _subsystem);
448  _cfC_SubdetPM.book(ib, _emap, _filter_QIE8, _subsystem);
449  _cfC_SubdetPM_QIE1011.book(ib, _emap, _filter_QIE1011, _subsystem);
450  _cSumQ_SubdetPM.book(ib, _emap, _filter_QIE8, _subsystem);
451  _cSumQ_SubdetPM_QIE1011.book(ib, _emap, _filter_QIE1011, _subsystem);
452  _cSumQ_depth.book(ib, _emap, _subsystem);
453  _cSumQvsLS_SubdetPM.book(ib, _emap, _filter_QIE8, _subsystem);
454  _cSumQvsLS_SubdetPM_QIE1011.book(ib, _emap, _filter_QIE1011, _subsystem);
455  _cDigiSize_Crate.book(ib, _emap, _subsystem);
456  _cADCvsTS_SubdetPM.book(ib, _emap, _filter_QIE8, _subsystem);
457  _cADCvsTS_SubdetPM_QIE1011.book(ib, _emap, _filter_QIE1011, _subsystem);
458 
459  if (_ptype != fOffline) { // hidefed2crate
460  _cShapeCut_FED.book(ib, _emap, _subsystem);
461  _cTimingCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
462  _cTimingCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
463  _cTimingCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
464  _cTimingCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
465  _cTimingCutvsLS_FED.book(ib, _emap, _subsystem);
466  _cOccupancy_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
467  _cOccupancy_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
468  _cOccupancy_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
469  _cOccupancy_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
470  _cOccupancyCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
471  _cOccupancyCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
472  _cOccupancyCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
473  _cOccupancyCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
474  _cDigiSize_FED.book(ib, _emap, _subsystem);
475  }
476 
477  _cTimingCut_SubdetPM.book(ib, _emap, _subsystem);
478  _cTimingCut_depth.book(ib, _emap, _subsystem);
479  _cTimingCutvsLS_SubdetPM.book(ib, _emap, _subsystem);
480 
481  _cOccupancyvsLS_Subdet.book(ib, _emap, _subsystem);
482  _cOccupancy_depth.book(ib, _emap, _subsystem);
483  _cOccupancyCut_depth.book(ib, _emap, _subsystem);
484 
485  _cLETDCTimevsADC_SubdetPM.book(ib, _emap, _subsystem);
486  _cLETDCvsADC_SubdetPM.book(ib, _emap, _subsystem);
487  _cLETDCvsTS_SubdetPM.book(ib, _emap, _subsystem);
488  _cLETDCTime_SubdetPM.book(ib, _emap, _subsystem);
489  _cLETDCTime_depth.book(ib, _emap, _subsystem);
490 
491  _cBadTDCValues_SubdetPM.book(ib, _emap, _subsystem);
492  _cBadTDCvsBX_SubdetPM.book(ib, _emap, _subsystem);
493  _cBadTDCvsLS_SubdetPM.book(ib, _emap, _subsystem);
494  _cBadTDCCount_depth.book(ib, _emap, _subsystem);
495 
496  _cCapidMinusBXmod4_SubdetPM.book(ib, _emap, _subsystem);
497  for (int i = 0; i < 4; ++i) {
498  char aux[10];
499  sprintf(aux, "%d_uTCA", i);
500  _cCapidMinusBXmod4_CrateSlotuTCA[i].book(ib, _subsystem, aux);
501 
502  sprintf(aux, "%d_VME", i);
503  _cCapidMinusBXmod4_CrateSlotVME[i].book(ib, _subsystem, aux);
504  }
505 
506  // BOOK HISTOGRAMS that are only for Online
507  _ehashmap.initialize(_emap, electronicsmap::fD2EHashMap);
508  _dhashmap.initialize(_emap, electronicsmap::fE2DHashMap);
509 
510  if (_ptype == fOnline || _ptype == fLocal) {
511  _cOccupancy_Crate.book(ib, _emap, _subsystem);
512  _cOccupancy_CrateSlot.book(ib, _emap, _subsystem);
513  }
514 
515  if (_ptype==fOnline)
516  {
517  _cQ2Q12CutvsLS_FEDHF.book(ib, _emap, _filter_FEDHF, _subsystem);
518  _cSumQvsBX_SubdetPM.book(ib, _emap, _filter_QIE8, _subsystem);
519  _cSumQvsBX_SubdetPM_QIE1011.book(ib, _emap, _filter_QIE1011, _subsystem);
520  _cDigiSizevsLS_FED.book(ib, _emap, _subsystem);
521  _cTimingCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
522  _cTimingCutvsieta_Subdet.book(ib, _emap, _subsystem);
523  _cOccupancyCutvsLS_Subdet.book(ib, _emap, _subsystem);
524  _cOccupancyCutvsBX_Subdet.book(ib, _emap, _subsystem);
525  _cOccupancyvsiphi_SubdetPM.book(ib, _emap, _subsystem);
526  _cOccupancyvsieta_Subdet.book(ib, _emap, _subsystem);
527  _cOccupancyCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
528  _cOccupancyCutvsieta_Subdet.book(ib, _emap, _subsystem);
529 // _cOccupancyCutvsSlotvsLS_HFPM.book(ib, _emap, _filter_QIE1011, _subsystem);
530  _cOccupancyCutvsiphivsLS_SubdetPM.book(ib, _emap, _subsystem);
531  _cSummaryvsLS_FED.book(ib, _emap, _subsystem);
532  _cSummaryvsLS.book(ib, _subsystem);
533 
534  _xUniHF.book(_emap, _filter_FEDHF);
535  _xNChs.book(_emap);
536  _xNChsNominal.book(_emap);
537  _xUni.book(_emap);
538  _xDigiSize.book(_emap);
539  _xBadCapid.book(_emap);
540 
541  // Manually book LED monitoring histogram, to get custom axis
542  ib.setCurrentFolder(_subsystem+"/"+_name+"/LED");
543  _meLEDMon = ib.book2D("LED_ADCvsBX", "Pin diode ADC vs BX", 99, -0.5, 3564-0.5, 64, -0.5, 255.5);
544  _meLEDMon->setAxisTitle("BX", 1);
545  _meLEDMon->setAxisTitle("ADC", 2);
546 
547  // just PER HF FED RECORD THE #CHANNELS
548  // ONLY WAY TO DO THAT AUTOMATICALLY AND W/O HARDCODING 1728
549  // or ANY OTHER VALUES LIKE 2592, 2192
550  std::vector<HcalGenericDetId> gids = _emap->allPrecisionId();
551  for (std::vector<HcalGenericDetId>::const_iterator it=gids.begin();
552  it!=gids.end(); ++it)
553  {
554  if (!it->isHcalDetId())
555  continue;
556  HcalDetId did(it->rawId());
557  if (_xQuality.exists(did))
558  {
559  HcalChannelStatus cs(it->rawId(), _xQuality.get(
560  HcalDetId(*it)));
561  if (
562  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
563  cs.isBitSet(HcalChannelStatus::HcalCellDead))
564  continue;
565  }
566  HcalElectronicsId eid = HcalElectronicsId(_ehashmap.lookup(did));
567  _xNChsNominal.get(eid)++; // he will know the nominal #channels per FED
568  }
569  }
570 
571  // MARK THESE HISTOGRAMS AS LUMI BASED FOR OFFLINE PROCESSING
572  if (_ptype==fOffline)
573  {
574  _cDigiSize_Crate.setLumiFlag();
575  //_cDigiSize_FED.setLumiFlag();
576  _cOccupancy_depth.setLumiFlag();
577  }
578 
579  // book Number of Events vs LS histogram
580  ib.setCurrentFolder(_subsystem+"/RunInfo");
581  meNumEvents1LS = ib.book1D("NumberOfEvents", "NumberOfEvents",
582  1, 0, 1);
583  meNumEvents1LS->setLumiFlag();
584 
585  // book the flag for unknown ids and the online guy as well
586  ib.setCurrentFolder(_subsystem+"/"+_name);
587  meUnknownIds1LS = ib.book1D("UnknownIds", "UnknownIds",
588  1, 0, 1);
589  _unknownIdsPresent = false;
590  meUnknownIds1LS->setLumiFlag();
591 
592  _ledSignalPresent = false;
593  ib.setCurrentFolder(_subsystem+"/"+_name+"/LED");
594  _meLEDEventCount = ib.book1D("LEDEventCount", "LEDEventCount", 1, 0, 2);
595  _meLEDEventCount->setLumiFlag();
596 }
597 
598 /* virtual */ void DigiTask::_resetMonitors(hcaldqm::UpdateFreq uf)
599 {
600  DQTask::_resetMonitors(uf);
601 
602  switch(uf)
603  {
604  case hcaldqm::f1LS:
605  _unknownIdsPresent = false;
606  break;
607  case hcaldqm::f50LS:
608  // ^^^ONLINE ONLY!
609  if (_ptype==fOnline)
610  _cOccupancyvsiphi_SubdetPM.reset();
611  // ^^^
612  break;
613  default:
614  break;
615  }
616 }
617 
618 /* virtual */ void DigiTask::_process(edm::Event const& e,
619  edm::EventSetup const&)
620 {
621  edm::Handle<HBHEDigiCollection> chbhe;
622  edm::Handle<QIE11DigiCollection> che_qie11;
623  edm::Handle<HODigiCollection> cho;
624  edm::Handle<QIE10DigiCollection> chf;
625 
626  if (!e.getByToken(_tokHBHE, chbhe))
627  _logger.dqmthrow("Collection HBHEDigiCollection isn't available"
628  + _tagHBHE.label() + " " + _tagHBHE.instance());
629  if (!e.getByToken(_tokHE, che_qie11))
630  _logger.dqmthrow("Collection QIE11DigiCollection isn't available"
631  + _tagHE.label() + " " + _tagHE.instance());
632  if (!e.getByToken(_tokHO, cho))
633  _logger.dqmthrow("Collection HODigiCollection isn't available"
634  + _tagHO.label() + " " + _tagHO.instance());
635  if (!e.getByToken(_tokHF, chf))
636  _logger.dqmthrow("Collection QIE10DigiCollection isn't available"
637  + _tagHF.label() + " " + _tagHF.instance());
638 
639  // extract some info per event
640  int bx = e.bunchCrossing();
641  meNumEvents1LS->Fill(0.5); // just increment
642 
643  // To fill histograms outside of the loop, you need to determine if there were
644  // any valid det ids first
645  uint32_t rawidValid = 0;
646  uint32_t rawidHBValid = 0;
647  uint32_t rawidHEValid = 0;
648 
649  // HB collection
650  int numChs = 0;
651  int numChsCut = 0;
652  int numChsHE = 0;
653  int numChsCutHE = 0;
654  for (HBHEDigiCollection::const_iterator it=chbhe->begin(); it!=chbhe->end();
655  ++it)
656  {
657  // Explicit check on the DetIds present in the Collection
658  HcalDetId const& did = it->id();
659  if (did.subdet() != HcalBarrel) {
660  continue;
661  }
662  uint32_t rawid = _ehashmap.lookup(did);
663  if (rawid == 0) {
664  meUnknownIds1LS->Fill(1);
665  _unknownIdsPresent=true;
666  continue;
667  } else {
668  if (did.subdet()==HcalBarrel) {
669  rawidHBValid = did.rawId();
670  } else if (did.subdet()==HcalEndcap) {
671  rawidHEValid = did.rawId();
672  }
673  }
674  HcalElectronicsId const& eid(rawid);
675 
676  // filter out channels that are masked out
677  if (_xQuality.exists(did))
678  {
679  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
680  if (
681  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
682  cs.isBitSet(HcalChannelStatus::HcalCellDead))
683  continue;
684  }
685 
686  if (_ptype == fOnline) {
687  short this_capidmbx = (it->sample(it->presamples()).capid() - bx) % 4;
688  if (this_capidmbx < 0) {
689  this_capidmbx += 4;
690  }
691  _cCapidMinusBXmod4_SubdetPM.fill(did, this_capidmbx);
692  bool good_capidmbx = (_capidmbx[did.subdet()] == this_capidmbx);
693  if (!good_capidmbx) {
694  _xBadCapid.get(eid)++;
695  }
696  if (eid.isVMEid()) {
697  _cCapidMinusBXmod4_CrateSlotVME[this_capidmbx].fill(eid);
698 
699  } else {
700  _cCapidMinusBXmod4_CrateSlotuTCA[this_capidmbx].fill(eid);
701  }
702  }
703 
704  //double sumQ = hcaldqm::utilities::sumQ<HBHEDataFrame>(*it, 2.5, 0, it->size()-1);
705  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<HBHEDataFrame>(_dbService, did, *it);
706  double sumQ = hcaldqm::utilities::sumQDB<HBHEDataFrame>(_dbService, digi_fC, did, *it, 0, it->size()-1);
707 
708  _cSumQ_SubdetPM.fill(did, sumQ);
709  _cOccupancy_depth.fill(did);
710  if (_ptype == fOnline || _ptype == fLocal) {
711  _cOccupancy_Crate.fill(eid);
712  _cOccupancy_CrateSlot.fill(eid);
713  }
714  if (_ptype==fOnline)
715  {
716  _cDigiSizevsLS_FED.fill(eid, _currentLS, it->size());
717  it->size()!=_refDigiSize[did.subdet()]?
718  _xDigiSize.get(eid)++:_xDigiSize.get(eid)+=0;
719  _cOccupancyvsiphi_SubdetPM.fill(did);
720  _cOccupancyvsieta_Subdet.fill(did);
721  }
722  _cDigiSize_Crate.fill(eid, it->size());
723  if (_ptype != fOffline) { // hidefed2crate
724  _cDigiSize_FED.fill(eid, it->size());
725  if (eid.isVMEid())
726  {
727  _cOccupancy_FEDVME.fill(eid);
728  _cOccupancy_ElectronicsVME.fill(eid);
729  }
730  else
731  {
732  _cOccupancy_FEDuTCA.fill(eid);
733  _cOccupancy_ElectronicsuTCA.fill(eid);
734  /*
735  if (!it->validate(0, it->size()))
736  {
737  _cCapIdRots_depth.fill(did);
738  _cCapIdRots_FEDuTCA.fill(eid, 1);
739  }*/
740  }
741  }
742 
743  for (int i=0; i<it->size(); i++)
744  {
745  _cADC_SubdetPM.fill(did, it->sample(i).adc());
746  _cfC_SubdetPM.fill(did, it->sample(i).nominal_fC());
747  if (_ptype != fOffline) { // hidefed2crate
748  _cADCvsTS_SubdetPM.fill(did, i, it->sample(i).adc());
749  if (sumQ>_cutSumQ_HBHE) {
750  _cShapeCut_FED.fill(eid, i, it->sample(i).nominal_fC());
751  }
752  }
753  }
754 
755  if (sumQ>_cutSumQ_HBHE)
756  {
757  //double timing = hcaldqm::utilities::aveTS<HBHEDataFrame>(*it, 2.5, 0, it->size()-1);
758  double timing = hcaldqm::utilities::aveTSDB<HBHEDataFrame>(_dbService, digi_fC, did, *it, 0, it->size()-1);
759  _cTimingCut_SubdetPM.fill(did, timing);
760  _cTimingCut_depth.fill(did, timing);
761  _cOccupancyCut_depth.fill(did);
762  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
763  if (_ptype != fOffline) { // hidefed2crate
764  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
765  }
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  }
777  if (_ptype != fOffline) { // hidefed2crate
778  if (eid.isVMEid())
779  {
780  _cTimingCut_FEDVME.fill(eid, timing);
781  _cTimingCut_ElectronicsVME.fill(eid, timing);
782  _cOccupancyCut_FEDVME.fill(eid);
783  _cOccupancyCut_ElectronicsVME.fill(eid);
784  }
785  else
786  {
787  _cTimingCut_FEDuTCA.fill(eid, timing);
788  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
789  _cOccupancyCut_FEDuTCA.fill(eid);
790  _cOccupancyCut_ElectronicsuTCA.fill(eid);
791  }
792  }
793  did.subdet()==HcalBarrel?numChsCut++:numChsCutHE++;
794  }
795  did.subdet()==HcalBarrel?numChs++:numChsHE++;
796  }
797 
798  // HE QIE11 collection
799  for (QIE11DigiCollection::const_iterator it=che_qie11->begin(); it!=che_qie11->end();
800  ++it)
801  {
802  const QIE11DataFrame digi = static_cast<const QIE11DataFrame>(*it);
803 
804  // Explicit check on the DetIds present in the Collection
805  HcalDetId const& did = digi.detid();
806  if (did.subdet() != HcalEndcap) {
807  // LED monitoring from calibration channels
808  if (did.subdet() == HcalOther) {
809  HcalOtherDetId hodid(digi.detid());
810  if (hodid.subdet() == HcalCalibration) {
811  if (did.depth() == 10) {
812  _ledSignalPresent = false;
813  for (int i=0; i<digi.samples(); i++) {
814  if (_ptype == fOnline) {
815  _meLEDMon->Fill(bx, digi[i].adc());
816  }
817  if (digi[i].adc() > _thresh_led) {
818  _ledSignalPresent = true;
819  }
820  }
821  if (_ledSignalPresent) {
822  _meLEDEventCount->Fill(1);
823  }
824  }
825  }
826  }
827 
828  continue;
829  }
830  uint32_t rawid = _ehashmap.lookup(did);
831  if (rawid == 0) {
832  meUnknownIds1LS->Fill(1);
833  _unknownIdsPresent=true;
834  continue;
835  } else {
836  if (did.subdet()==HcalBarrel) { // Note: since this is HE, we obviously expect did.subdet() always to be HcalEndcap, but QIE11DigiCollection will have HB for Run 3.
837  rawidHBValid = did.rawId();
838  } else if (did.subdet()==HcalEndcap) {
839  rawidHEValid = did.rawId();
840  }
841  }
842  HcalElectronicsId const& eid(rawid);
843 
844  // filter out channels that are masked out
845  if (_xQuality.exists(did))
846  {
847  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
848  if (
849  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
850  cs.isBitSet(HcalChannelStatus::HcalCellDead))
851  continue;
852  }
853 
854  // (capid - BX) % 4
855  if (_ptype == fOnline) {
856  short soi = -1;
857  for (int i=0; i<digi.samples(); i++) {
858  if (digi[i].soi()) {
859  soi = i;
860  break;
861  }
862  }
863  short this_capidmbx = (digi[soi].capid() - bx) % 4;
864  if (this_capidmbx < 0) {
865  this_capidmbx += 4;
866  }
867  _cCapidMinusBXmod4_SubdetPM.fill(did, this_capidmbx);
868  bool good_capidmbx = (_capidmbx[did.subdet()] == this_capidmbx);
869  if (!good_capidmbx) {
870  _xBadCapid.get(eid)++;
871  }
872  if (eid.isVMEid()) {
873  _cCapidMinusBXmod4_CrateSlotVME[this_capidmbx].fill(eid);
874 
875  } else {
876  _cCapidMinusBXmod4_CrateSlotuTCA[this_capidmbx].fill(eid);
877  }
878  }
879 
880  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE11DataFrame>(_dbService, did, digi);
881  double sumQ = hcaldqm::utilities::sumQDB<QIE11DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples()-1);
882 
883  _cSumQ_SubdetPM_QIE1011.fill(did, sumQ);
884  _cOccupancy_depth.fill(did);
885  if (_ptype == fOnline || _ptype == fLocal) {
886  _cOccupancy_Crate.fill(eid);
887  _cOccupancy_CrateSlot.fill(eid);
888  }
889  if (_ptype==fOnline)
890  {
891  _cDigiSizevsLS_FED.fill(eid, _currentLS, digi.samples());
892  digi.samples()!=_refDigiSize[did.subdet()]?
893  _xDigiSize.get(eid)++:_xDigiSize.get(eid)+=0;
894  _cOccupancyvsiphi_SubdetPM.fill(did);
895  _cOccupancyvsieta_Subdet.fill(did);
896  }
897  _cDigiSize_Crate.fill(eid, digi.samples());
898  if (_ptype != fOffline) { // hidefed2crate
899  _cDigiSize_FED.fill(eid, digi.samples());
900  if (eid.isVMEid())
901  {
902  _cOccupancy_FEDVME.fill(eid);
903  _cOccupancy_ElectronicsVME.fill(eid);
904  }
905  else
906  {
907  _cOccupancy_FEDuTCA.fill(eid);
908  _cOccupancy_ElectronicsuTCA.fill(eid);
909  /*
910  if (!digi.validate(0, digi.size()))
911  {
912  _cCapIdRots_depth.fill(did);
913  _cCapIdRots_FEDuTCA.fill(eid, 1);
914  }*/
915  }
916  }
917  for (int i=0; i<digi.samples(); i++) {
918  double q = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE11DataFrame>(_dbService, digi_fC, did, digi, i);
919  _cADC_SubdetPM_QIE1011.fill(did, digi[i].adc());
920  _cfC_SubdetPM_QIE1011.fill(did, q);
921  _cLETDCvsADC_SubdetPM.fill(did, digi[i].adc(), digi[i].tdc());
922  _cLETDCvsTS_SubdetPM.fill(did, (int)i, digi[i].tdc());
923  if (digi[i].tdc() <50) {
924  double time = i*25. + (digi[i].tdc() / 2.);
925  _cLETDCTime_SubdetPM.fill(did, time);
926  _cLETDCTime_depth.fill(did, time);
927  _cLETDCTimevsADC_SubdetPM.fill(did, digi[i].adc(), time);
928  }
929  // Bad TDC values: 50-61 should never happen in QIE10 or QIE11, but we saw some in 2017 data.
930  if ((50 <= digi[i].tdc()) && (digi[i].tdc() <= 61)) {
931  _cBadTDCValues_SubdetPM.fill(did, digi[i].tdc());
932  _cBadTDCvsBX_SubdetPM.fill(did, bx);
933  _cBadTDCvsLS_SubdetPM.fill(did, _currentLS);
934  _cBadTDCCount_depth.fill(did);
935  }
936  if (_ptype != fOffline) { // hidefed2crate
937  _cADCvsTS_SubdetPM_QIE1011.fill(did, i, digi[i].adc());
938  if (sumQ>_cutSumQ_HE) {
939  _cShapeCut_FED.fill(eid, i, q);
940  }
941  }
942  }
943 
944  if (sumQ>_cutSumQ_HE)
945  {
946  //double timing = hcaldqm::utilities::aveTS_v10<QIE11DataFrame>(digi, 2.5, 0,digi.samples()-1);
947  double timing = hcaldqm::utilities::aveTSDB<QIE11DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples()-1);
948  _cTimingCut_SubdetPM.fill(did, timing);
949  _cTimingCut_depth.fill(did, timing);
950  _cOccupancyCut_depth.fill(did);
951  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
952  if (_ptype != fOffline) { // hidefed2crate
953  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
954  }
955  _cSumQ_depth.fill(did, sumQ);
956  _cSumQvsLS_SubdetPM_QIE1011.fill(did, _currentLS, sumQ);
957  if (_ptype==fOnline)
958  {
959  _cSumQvsBX_SubdetPM_QIE1011.fill(did, bx, sumQ);
960  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
961  _cTimingCutvsieta_Subdet.fill(did, timing);
962  _cOccupancyCutvsiphi_SubdetPM.fill(did);
963  _cOccupancyCutvsieta_Subdet.fill(did);
964  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
965  }
966  if (_ptype != fOffline) { // hidefed2crate
967  if (eid.isVMEid())
968  {
969  _cTimingCut_FEDVME.fill(eid, timing);
970  _cTimingCut_ElectronicsVME.fill(eid, timing);
971  _cOccupancyCut_FEDVME.fill(eid);
972  _cOccupancyCut_ElectronicsVME.fill(eid);
973  }
974  else
975  {
976  _cTimingCut_FEDuTCA.fill(eid, timing);
977  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
978  _cOccupancyCut_FEDuTCA.fill(eid);
979  _cOccupancyCut_ElectronicsuTCA.fill(eid);
980  }
981  }
982  did.subdet()==HcalBarrel?numChsCut++:numChsCutHE++;
983  }
984  did.subdet()==HcalBarrel?numChs++:numChsHE++;
985  }
986 
987  if (rawidHBValid!=0 && rawidHEValid!=0)
988  {
989  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHBValid), _currentLS,
990  numChs);
991  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHEValid), _currentLS,
992  numChsHE);
993  // ONLINE ONLY!
994  if (_ptype==fOnline)
995  {
996  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHBValid),
997  _currentLS, numChsCut);
998  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHBValid), bx,
999  numChsCut);
1000  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHEValid),
1001  _currentLS, numChsCutHE);
1002  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHEValid), bx,
1003  numChsCutHE);
1004  }
1005  // ^^^ONLINE ONLY!
1006  }
1007  numChs=0;
1008  numChsCut = 0;
1009 
1010  // reset
1011  rawidValid = 0;
1012 
1013 
1014  // HO collection
1015  for (HODigiCollection::const_iterator it=cho->begin(); it!=cho->end();
1016  ++it)
1017  {
1018  // Explicit check on the DetIds present in the Collection
1019  HcalDetId const& did = it->id();
1020  if (did.subdet() != HcalOuter) {
1021  continue;
1022  }
1023  uint32_t rawid = _ehashmap.lookup(did);
1024  if (rawid == 0) {
1025  meUnknownIds1LS->Fill(1);
1026  _unknownIdsPresent = true;
1027  continue;
1028  } else {
1029  rawidValid = did.rawId();
1030  }
1031  HcalElectronicsId const& eid(rawid);
1032 
1033  // filter out channels that are masked out
1034  if (_xQuality.exists(did))
1035  {
1036  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
1037  if (
1038  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
1039  cs.isBitSet(HcalChannelStatus::HcalCellDead))
1040  continue;
1041  }
1042 
1043  if (_ptype == fOnline) {
1044  short this_capidmbx = (it->sample(it->presamples()).capid() - bx) % 4;
1045  if (this_capidmbx < 0) {
1046  this_capidmbx += 4;
1047  }
1048  _cCapidMinusBXmod4_SubdetPM.fill(did, this_capidmbx);
1049  bool good_capidmbx = (_capidmbx[did.subdet()] == this_capidmbx);
1050  if (!good_capidmbx) {
1051  _xBadCapid.get(eid)++;
1052  }
1053  if (eid.isVMEid()) {
1054  _cCapidMinusBXmod4_CrateSlotVME[this_capidmbx].fill(eid);
1055 
1056  } else {
1057  _cCapidMinusBXmod4_CrateSlotuTCA[this_capidmbx].fill(eid);
1058  }
1059  }
1060 
1061  //double sumQ = hcaldqm::utilities::sumQ<HODataFrame>(*it, 8.5, 0, it->size()-1);
1062  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<HODataFrame>(_dbService, did, *it);
1063  double sumQ = hcaldqm::utilities::sumQDB<HODataFrame>(_dbService, digi_fC, did, *it, 0, it->size()-1);
1064 
1065  _cSumQ_SubdetPM.fill(did, sumQ);
1066  _cOccupancy_depth.fill(did);
1067  if (_ptype==fOnline)
1068  {
1069  _cDigiSizevsLS_FED.fill(eid, _currentLS, it->size());
1070  it->size()!=_refDigiSize[did.subdet()]?
1071  _xDigiSize.get(eid)++:_xDigiSize.get(eid)+=0;
1072  _cOccupancyvsiphi_SubdetPM.fill(did);
1073  _cOccupancyvsieta_Subdet.fill(did);
1074  }
1075  _cDigiSize_Crate.fill(eid, it->size());
1076  if (_ptype != fOffline) { // hidefed2crate
1077  _cDigiSize_FED.fill(eid, it->size());
1078  if (eid.isVMEid())
1079  {
1080  _cOccupancy_FEDVME.fill(eid);
1081  _cOccupancy_ElectronicsVME.fill(eid);
1082  /*
1083  if (!it->validate(0, it->size()))
1084  _cCapIdRots_FEDVME.fill(eid, 1);
1085  */
1086  }
1087  else
1088  {
1089  _cOccupancy_FEDuTCA.fill(eid);
1090  _cOccupancy_ElectronicsuTCA.fill(eid);
1091  /*
1092  if (!it->validate(0, it->size()))
1093  _cCapIdRots_FEDuTCA.fill(eid, 1);*/
1094  }
1095  }
1096 
1097  for (int i=0; i<it->size(); i++)
1098  {
1099  _cADC_SubdetPM.fill(did, it->sample(i).adc());
1100  _cfC_SubdetPM.fill(did, it->sample(i).nominal_fC());
1101  if (_ptype != fOffline) { // hidefed2crate
1102  _cADCvsTS_SubdetPM.fill(did, i, it->sample(i).adc());
1103  if (sumQ>_cutSumQ_HO)
1104  _cShapeCut_FED.fill(eid, i, it->sample(i).nominal_fC());
1105  }
1106  }
1107 
1108  if (sumQ>_cutSumQ_HO)
1109  {
1110  //double timing = hcaldqm::utilities::aveTS<HODataFrame>(*it, 8.5, 0,it->size()-1);
1111  double timing = hcaldqm::utilities::aveTSDB<HODataFrame>(_dbService, digi_fC, did, *it, 0, it->size()-1);
1112  _cSumQ_depth.fill(did, sumQ);
1113  _cSumQvsLS_SubdetPM.fill(did, _currentLS, sumQ);
1114  _cOccupancyCut_depth.fill(did);
1115  _cTimingCut_SubdetPM.fill(did, timing);
1116  _cTimingCut_depth.fill(did, timing);
1117  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
1118  if (_ptype != fOffline) { // hidefed2crate
1119  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
1120  }
1121  if (_ptype==fOnline)
1122  {
1123  _cSumQvsBX_SubdetPM.fill(did, bx, sumQ);
1124  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
1125  _cTimingCutvsieta_Subdet.fill(did, timing);
1126  _cOccupancyCutvsiphi_SubdetPM.fill(did);
1127  _cOccupancyCutvsieta_Subdet.fill(did);
1128  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
1129  }
1130  if (_ptype != fOffline) { // hidefed2crate
1131  if (eid.isVMEid())
1132  {
1133  _cTimingCut_FEDVME.fill(eid, timing);
1134  _cTimingCut_ElectronicsVME.fill(eid, timing);
1135  _cOccupancyCut_FEDVME.fill(eid);
1136  _cOccupancyCut_ElectronicsVME.fill(eid);
1137  }
1138  else
1139  {
1140  _cTimingCut_FEDuTCA.fill(eid, timing);
1141  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
1142  _cOccupancyCut_FEDuTCA.fill(eid);
1143  _cOccupancyCut_ElectronicsuTCA.fill(eid);
1144  }
1145  }
1146  numChsCut++;
1147  }
1148  numChs++;
1149  }
1150 
1151  if (rawidValid!=0)
1152  {
1153  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS,
1154  numChs);
1155 
1156  if (_ptype==fOnline)
1157  {
1158  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid),
1159  _currentLS, numChsCut);
1160  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid), bx,
1161  numChsCut);
1162  }
1163  }
1164  numChs=0; numChsCut=0;
1165 
1166  // reset
1167  rawidValid = 0;
1168 
1169  // HF collection
1170  if (_qie10InConditions) {
1171  for (QIE10DigiCollection::const_iterator it=chf->begin(); it!=chf->end(); ++it) {
1172  const QIE10DataFrame digi = static_cast<const QIE10DataFrame>(*it);
1173 
1174  // Explicit check on the DetIds present in the Collection
1175  HcalDetId const& did = digi.detid();
1176  if (did.subdet() != HcalForward) {
1177  continue;
1178  }
1179 
1180  uint32_t rawid = _ehashmap.lookup(did);
1181  if (rawid == 0) {
1182  meUnknownIds1LS->Fill(1);
1183  _unknownIdsPresent=true;
1184  continue;
1185  } else {
1186  rawidValid = did.rawId();
1187  }
1188  HcalElectronicsId const& eid(rawid);
1189 
1190  // filter out channels that are masked out
1191  if (_xQuality.exists(did))
1192  {
1193  HcalChannelStatus cs(did.rawId(), _xQuality.get(did));
1194  if (
1195  cs.isBitSet(HcalChannelStatus::HcalCellMask) ||
1196  cs.isBitSet(HcalChannelStatus::HcalCellDead))
1197  continue;
1198  }
1199 
1200  // (capid - BX) % 4
1201  if (_ptype == fOnline) {
1202  short soi = -1;
1203  for (int i=0; i<digi.samples(); i++) {
1204  if (digi[i].soi()) {
1205  soi = i;
1206  break;
1207  }
1208  }
1209  short this_capidmbx = (digi[soi].capid() - bx) % 4;
1210  if (this_capidmbx < 0) {
1211  this_capidmbx += 4;
1212  }
1213  _cCapidMinusBXmod4_SubdetPM.fill(did, this_capidmbx);
1214  bool good_capidmbx = (_capidmbx[did.subdet()] == this_capidmbx);
1215  if (!good_capidmbx) {
1216  _xBadCapid.get(eid)++;
1217  }
1218  if (eid.isVMEid()) {
1219  _cCapidMinusBXmod4_CrateSlotVME[this_capidmbx].fill(eid);
1220 
1221  } else {
1222  _cCapidMinusBXmod4_CrateSlotuTCA[this_capidmbx].fill(eid);
1223  }
1224  }
1225 
1226  CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE10DataFrame>(_dbService, did, digi);
1227  double sumQ = hcaldqm::utilities::sumQDB<QIE10DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples()-1);
1228  //double sumQ = hcaldqm::utilities::sumQ_v10<QIE10DataFrame>(digi, 2.5, 0, digi.samples()-1);
1229 
1230  //if (!_filter_QIE1011.filter(did)) {
1231  _cSumQ_SubdetPM_QIE1011.fill(did, sumQ);
1232  //}
1233  _cOccupancy_depth.fill(did);
1234  if (_ptype==fOnline)
1235  {
1236  _xNChs.get(eid)++;
1237  _cDigiSizevsLS_FED.fill(eid, _currentLS, digi.samples());
1238  digi.samples()!=_refDigiSize[did.subdet()]?
1239  _xDigiSize.get(eid)++:_xDigiSize.get(eid)+=0;
1240  _cOccupancyvsiphi_SubdetPM.fill(did);
1241  _cOccupancyvsieta_Subdet.fill(did);
1242  }
1243  _cDigiSize_Crate.fill(eid, digi.samples());
1244  if (_ptype != fOffline) { // hidefed2crate
1245  _cDigiSize_FED.fill(eid, digi.samples());
1246  if (eid.isVMEid())
1247  {
1248  _cOccupancy_FEDVME.fill(eid);
1249  _cOccupancy_ElectronicsVME.fill(eid);
1250  /*
1251  if (!it->validate(0, it->size()))
1252  _cCapIdRots_FEDVME.fill(eid, 1);*/
1253  }
1254  else
1255  {
1256  _cOccupancy_FEDuTCA.fill(eid);
1257  _cOccupancy_ElectronicsuTCA.fill(eid);
1258  /*
1259  if (!it->validate(0, it->size()))
1260  _cCapIdRots_FEDuTCA.fill(eid, 1);*/
1261  }
1262  }
1263 
1264  for (int i=0; i<digi.samples(); i++)
1265  {
1266  double q = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, i);
1267  //if (!_filter_QIE1011.filter(did)) {
1268  _cADC_SubdetPM_QIE1011.fill(did, digi[i].adc());
1269  _cfC_SubdetPM_QIE1011.fill(did, q);
1270  _cLETDCvsADC_SubdetPM.fill(did, digi[i].adc(), digi[i].le_tdc());
1271  _cLETDCvsTS_SubdetPM.fill(did, (int)i, digi[i].le_tdc());
1272  if (digi[i].le_tdc() <50) {
1273  double time = i*25. + (digi[i].le_tdc() / 2.);
1274  _cLETDCTime_SubdetPM.fill(did, time);
1275  _cLETDCTime_depth.fill(did, time);
1276  _cLETDCTimevsADC_SubdetPM.fill(did, digi[i].adc(), time);
1277  }
1278 
1279  // Bad TDC values: 50-61 should never happen in QIE10 or QIE11, but we are seeing some in 2017 data.
1280  if ((50 <= digi[i].le_tdc()) && (digi[i].le_tdc() <= 61)) {
1281  _cBadTDCValues_SubdetPM.fill(did, digi[i].le_tdc());
1282  _cBadTDCvsBX_SubdetPM.fill(did, bx);
1283  _cBadTDCvsLS_SubdetPM.fill(did, _currentLS);
1284  _cBadTDCCount_depth.fill(did);
1285  }
1286  if (_ptype != fOffline) { // hidefed2crate
1287  _cADCvsTS_SubdetPM_QIE1011.fill(did, (int)i, digi[i].adc());
1288  if (sumQ>_cutSumQ_HF)
1289  _cShapeCut_FED.fill(eid, (int)i, q);
1290  }
1291  //}
1292  }
1293 
1294  if (sumQ>_cutSumQ_HF)
1295  {
1296  double timing = hcaldqm::utilities::aveTS_v10<QIE10DataFrame>(digi, 2.5, 0,
1297  digi.samples()-1);
1298  double q1 = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, 1);
1299  double q2 = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, 2);
1300  double q2q12 = q2/(q1+q2);
1301  _cSumQ_depth.fill(did, sumQ);
1302  //if (!_filter_QIE1011.filter(did)) {
1303  _cSumQvsLS_SubdetPM_QIE1011.fill(did, _currentLS, sumQ);
1304  //}
1305  _cTimingCut_SubdetPM.fill(did, timing);
1306  _cTimingCut_depth.fill(did, timing);
1307  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
1308  if (_ptype==fOnline)
1309  {
1310  //if (!_filter_QIE1011.filter(did)) {
1311  _cSumQvsBX_SubdetPM_QIE1011.fill(did, bx, sumQ);
1312  //}
1313  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
1314  _cTimingCutvsieta_Subdet.fill(did, timing);
1315  _cOccupancyCutvsiphi_SubdetPM.fill(did);
1316  _cOccupancyCutvsieta_Subdet.fill(did);
1317  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
1318  // _cOccupancyCutvsSlotvsLS_HFPM.fill(did, _currentLS);
1319  _xUniHF.get(eid)++;
1320  }
1321  if (_ptype != fOffline) { // hidefed2crate
1322  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
1323  }
1324  _cOccupancyCut_depth.fill(did);
1325  if (!eid.isVMEid())
1326  if (_ptype==fOnline)
1327  _cQ2Q12CutvsLS_FEDHF.fill(eid, _currentLS, q2q12);
1328  if (_ptype != fOffline) { // hidefed2crate
1329  if (eid.isVMEid())
1330  {
1331  _cTimingCut_FEDVME.fill(eid, timing);
1332  _cTimingCut_ElectronicsVME.fill(eid, timing);
1333  _cOccupancyCut_FEDVME.fill(eid);
1334  _cOccupancyCut_ElectronicsVME.fill(eid);
1335  }
1336  else
1337  {
1338  _cTimingCut_FEDuTCA.fill(eid, timing);
1339  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
1340  _cOccupancyCut_FEDuTCA.fill(eid);
1341  _cOccupancyCut_ElectronicsuTCA.fill(eid);
1342  }
1343  }
1344  numChsCut++;
1345  }
1346  numChs++;
1347  }
1348  }
1349 
1350  if (rawidValid!=0)
1351  {
1352  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS,
1353  numChs);
1354 
1355  if (_ptype==fOnline)
1356  {
1357  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid),
1358  _currentLS, numChsCut);
1359  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid), bx,
1360  numChsCut);
1361  }
1362  }
1363 }
1364 
1365 /* virtual */ void DigiTask::beginLuminosityBlock(
1366  edm::LuminosityBlock const& lb, edm::EventSetup const& es)
1367 {
1368  DQTask::beginLuminosityBlock(lb, es);
1369 }
1370 
1371 /* virtual */ void DigiTask::endLuminosityBlock(edm::LuminosityBlock const& lb,
1372  edm::EventSetup const& es)
1373 {
1374  if (_ptype!=fOnline)
1375  return;
1376 
1377  for (uintCompactMap::const_iterator it=_xUniHF.begin();
1378  it!=_xUniHF.end(); ++it)
1379  {
1380  uint32_t hash1 = it->first;
1381  HcalElectronicsId eid1(hash1);
1382  double x1 = it->second;
1383 
1384  for (uintCompactMap::const_iterator jt=_xUniHF.begin();
1385  jt!=_xUniHF.end(); ++jt)
1386  {
1387  if (jt==it)
1388  continue;
1389  double x2 = jt->second;
1390  if (x2==0)
1391  continue;
1392  if (x1/x2<_thresh_unihf)
1393  _xUni.get(eid1)++;
1394  }
1395  }
1396 
1397  if (_ptype != fOffline) { // hidefed2crate
1398  for (std::vector<uint32_t>::const_iterator it=_vhashFEDs.begin();
1399  it!=_vhashFEDs.end(); ++it)
1400  {
1401  hcaldqm::flag::Flag fSum("DIGI");
1402  HcalElectronicsId eid = HcalElectronicsId(*it);
1403 
1404  std::vector<uint32_t>::const_iterator cit=std::find(
1405  _vcdaqEids.begin(), _vcdaqEids.end(), *it);
1406  if (cit==_vcdaqEids.end())
1407  {
1408  // not @cDAQ
1409  for (uint32_t iflag=0; iflag<_vflags.size(); iflag++)
1410  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag),
1411  int(hcaldqm::flag::fNCDAQ));
1412  _cSummaryvsLS.setBinContent(eid, _currentLS, int(hcaldqm::flag::fNCDAQ));
1413  continue;
1414  }
1415 
1416  // FED is @cDAQ
1417  if (hcaldqm::utilities::isFEDHBHE(eid) || hcaldqm::utilities::isFEDHF(eid) ||
1418  hcaldqm::utilities::isFEDHO(eid))
1419  {
1420  if (_xDigiSize.get(eid)>0)
1421  _vflags[fDigiSize]._state = hcaldqm::flag::fBAD;
1422  else
1423  _vflags[fDigiSize]._state = hcaldqm::flag::fGOOD;
1424 
1425  if (_xBadCapid.get(eid) > 0) {
1426  _vflags[fCapId]._state = hcaldqm::flag::fBAD;
1427  } else {
1428  _vflags[fCapId]._state = hcaldqm::flag::fGOOD;
1429  }
1430 
1431  if (hcaldqm::utilities::isFEDHF(eid))
1432  {
1433  double fr = double(_xNChs.get(eid))/double(
1434  _xNChsNominal.get(eid)*_evsPerLS);
1435  if (_runkeyVal==0 || _runkeyVal==4)
1436  {
1437  // only for pp or hi
1438  if (_xUni.get(eid)>0)
1439  _vflags[fUni]._state = hcaldqm::flag::fPROBLEMATIC;
1440  else
1441  _vflags[fUni]._state = hcaldqm::flag::fGOOD;
1442  }
1443  if (fr<0.95)
1444  _vflags[fNChsHF]._state = hcaldqm::flag::fBAD;
1445  else if (fr<1.0)
1446  _vflags[fNChsHF]._state = hcaldqm::flag::fPROBLEMATIC;
1447  else
1448  _vflags[fNChsHF]._state = hcaldqm::flag::fGOOD;
1449  }
1450  }
1451  if (_unknownIdsPresent)
1452  _vflags[fUnknownIds]._state = hcaldqm::flag::fBAD;
1453  else
1454  _vflags[fUnknownIds]._state = hcaldqm::flag::fGOOD;
1455 
1456  if (_ledSignalPresent)
1457  _vflags[fLED]._state = hcaldqm::flag::fBAD;
1458  else
1459  _vflags[fLED]._state = hcaldqm::flag::fGOOD;
1460 
1461  int iflag=0;
1462  for (std::vector<hcaldqm::flag::Flag>::iterator ft=_vflags.begin();
1463  ft!=_vflags.end(); ++ft)
1464  {
1465  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, iflag,
1466  int(ft->_state));
1467  fSum+=(*ft);
1468  iflag++;
1469 
1470  // reset!
1471  ft->reset();
1472  }
1473  _cSummaryvsLS.setBinContent(eid, _currentLS, fSum._state);
1474  }
1475  }
1476 
1477  _xDigiSize.reset(); _xUniHF.reset(); _xUni.reset();
1478  _xNChs.reset();
1479  _xBadCapid.reset();
1480 
1481  // in the end always do the DQTask::endLumi
1482  DQTask::endLuminosityBlock(lb, es);
1483 }
1484 
1485 DEFINE_FWK_MODULE(DigiTask);
1486 
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:88
hcaldqm::ContainerProf1D::initialize
void initialize(std::string const &folder, hashfunctions::HashType, quantity::Quantity *, quantity::Quantity *, int debug=0) override
Definition: ContainerProf1D.cc:26
QIE11DataFrame::detid
DetId detid() const
Get the detector id.
Definition: QIE11DataFrame.h:47
DigiTask::_cTimingCutvsiphi_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsiphi_SubdetPM
Definition: DigiTask.h:120
HcalElectronicsId::rawId
constexpr uint32_t rawId() const
Definition: HcalElectronicsId.h:65
DigiTask::_process
void _process(edm::Event const &, edm::EventSetup const &) override
Definition: DigiTask.cc:618
DigiTask::_xBadCapid
hcaldqm::ContainerXXX< uint32_t > _xBadCapid
Definition: DigiTask.h:164
DigiTask::_qie10InConditions
bool _qie10InConditions
Definition: DigiTask.h:202
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
DigiTask::_thresh_led
double _thresh_led
Definition: DigiTask.h:196
DigiTask::meNumEvents1LS
MonitorElement * meNumEvents1LS
Definition: DigiTask.h:192
DigiTask::_cADC_SubdetPM
hcaldqm::Container1D _cADC_SubdetPM
Definition: DigiTask.h:90
DigiTask::_cADC_SubdetPM_QIE1011
hcaldqm::Container1D _cADC_SubdetPM_QIE1011
Definition: DigiTask.h:98
DigiTask::endLuminosityBlock
void endLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &) override
Definition: DigiTask.cc:1371
hcaldqm::Container1D::initialize
virtual void initialize(std::string const &folder, hashfunctions::HashType, quantity::Quantity *, quantity::Quantity *qy=new quantity::ValueQuantity(quantity::fN), int debug=0)
Definition: Container1D.cc:41
DigiTask::_cOccupancyCutvsieta_Subdet
hcaldqm::Container1D _cOccupancyCutvsieta_Subdet
Definition: DigiTask.h:147
fwrapper::cs
auto_ptr< ClusterSequence > cs
Definition: fastjetfortran_madfks.cc:45
DigiTask::_cSumQvsBX_SubdetPM
hcaldqm::ContainerProf1D _cSumQvsBX_SubdetPM
Definition: DigiTask.h:95
hcaldqm::ContainerXXX::exists
virtual bool exists(HcalDetId const &)
Definition: ContainerXXX.h:265
mps_fire.i
i
Definition: mps_fire.py:269
hcaldqm::constants::CRATE_VME_MIN
int const CRATE_VME_MIN
Definition: Constants.h:111
DigiTask::_meLEDMon
MonitorElement * _meLEDMon
Definition: DigiTask.h:189
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:142
hcaldqm::ContainerProf2D::fill
void fill(HcalDetId const &) override
Definition: ContainerProf2D.cc:403
DigiTask::_tagHO
edm::InputTag _tagHO
Definition: DigiTask.h:48
DigiTask::_refDigiSize
std::map< HcalSubdetector, int > _refDigiSize
Definition: DigiTask.h:74
DigiTask::_cADCvsTS_SubdetPM
hcaldqm::Container2D _cADCvsTS_SubdetPM
Definition: DigiTask.h:107
hcaldqm::utilities::fed2crate
std::pair< uint16_t, uint16_t > fed2crate(int fed)
Definition: Utilities.cc:12
DigiTask::beginLuminosityBlock
void beginLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &) override
Definition: DigiTask.cc:1365
DigiTask::_tokHE
edm::EDGetTokenT< QIE11DigiCollection > _tokHE
Definition: DigiTask.h:51
DigiTask::_cutSumQ_HE
double _cutSumQ_HE
Definition: DigiTask.h:55
DigiTask::_cQ2Q12CutvsLS_FEDHF
hcaldqm::ContainerProf1D _cQ2Q12CutvsLS_FEDHF
Definition: DigiTask.h:124
DigiTask::meUnknownIds1LS
MonitorElement * meUnknownIds1LS
Definition: DigiTask.h:193
DigiTask::_tokHBHE
edm::EDGetTokenT< HBHEDigiCollection > _tokHBHE
Definition: DigiTask.h:50
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:579
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::_meLEDEventCount
MonitorElement * _meLEDEventCount
Definition: DigiTask.h:195
DigiTask::_cLETDCTimevsADC_SubdetPM
hcaldqm::Container2D _cLETDCTimevsADC_SubdetPM
Definition: DigiTask.h:167
DigiTask::_cutSumQ_HO
double _cutSumQ_HO
Definition: DigiTask.h:55
hcaldqm::ContainerXXX::begin
virtual CompactMap::const_iterator begin()
Definition: ContainerXXX.h:73
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
DigiTask::_capidmbx
std::map< HcalSubdetector, short > _capidmbx
Definition: DigiTask.h:204
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:47
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
DigiTask::_cCapidMinusBXmod4_CrateSlotuTCA
hcaldqm::ContainerSingle2D _cCapidMinusBXmod4_CrateSlotuTCA[4]
Definition: DigiTask.h:185
hcaldqm::ContainerSingle2D::initialize
virtual void initialize(std::string const &folder, quantity::Quantity *, quantity::Quantity *, quantity::Quantity *qz=new quantity::ValueQuantity(quantity::fN), int debug=0)
Definition: ContainerSingle2D.cc:57
hcaldqm::constants::SPIGOT_MIN
int const SPIGOT_MIN
Definition: Constants.h:136
hcaldqm::ContainerProf1D::book
void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="") override
Definition: ContainerProf1D.cc:47
DigiTask::_cSumQ_depth
hcaldqm::ContainerProf2D _cSumQ_depth
Definition: DigiTask.h:93
DigiTask::_tokHO
edm::EDGetTokenT< HODigiCollection > _tokHO
Definition: DigiTask.h:52
DigiTask::_cOccupancyCutvsiphi_SubdetPM
hcaldqm::Container1D _cOccupancyCutvsiphi_SubdetPM
Definition: DigiTask.h:146
edm::DataFrameContainer::begin
const_iterator begin() const
Definition: DataFrameContainer.h:159
hcaldqm::ContainerXXX::reset
virtual void reset()
Definition: ContainerXXX.h:364
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
DigiTask::_cSumQvsLS_SubdetPM_QIE1011
hcaldqm::ContainerProf1D _cSumQvsLS_SubdetPM_QIE1011
Definition: DigiTask.h:101
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:200
DigiTask::_cTimingCut_FEDVME
hcaldqm::ContainerProf2D _cTimingCut_FEDVME
Definition: DigiTask.h:113
DigiTask::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: DigiTask.cc:55
hcaldqm::constants::FIBER_VME_MIN
int const FIBER_VME_MIN
Definition: Constants.h:141
DigiTask::_cOccupancyCutvsiphivsLS_SubdetPM
hcaldqm::Container2D _cOccupancyCutvsiphivsLS_SubdetPM
Definition: DigiTask.h:149
DigiTask::_xDigiSize
hcaldqm::ContainerXXX< uint32_t > _xDigiSize
Definition: DigiTask.h:160
DigiTask::_cTimingCut_FEDuTCA
hcaldqm::ContainerProf2D _cTimingCut_FEDuTCA
Definition: DigiTask.h:114
q2
double q2[4]
Definition: TauolaWrapper.h:88
DigiTask::_cBadTDCvsLS_SubdetPM
hcaldqm::Container1D _cBadTDCvsLS_SubdetPM
Definition: DigiTask.h:176
hcaldqm::constants::FIBERCH_MIN
int const FIBERCH_MIN
Definition: Constants.h:151
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:166
hcaldqm::ContainerSingle2D::book
virtual void book(DQMStore::IBooker &, std::string subsystem="Hcal", std::string aux="")
Definition: ContainerSingle2D.cc:83
hcaldqm::utilities::getFEDVMEList
std::vector< int > getFEDVMEList(HcalElectronicsMap const *)
Definition: Utilities.cc:119
DigiTask::_filter_VME
hcaldqm::filter::HashFilter _filter_VME
Definition: DigiTask.h:81
hcaldqm::ContainerXXX::end
virtual CompactMap::const_iterator end()
Definition: ContainerXXX.h:75
DigiTask::_cLETDCvsTS_SubdetPM
hcaldqm::Container2D _cLETDCvsTS_SubdetPM
Definition: DigiTask.h:169
HcalDetId::depth
int depth() const
get the tower depth
Definition: HcalDetId.h:162
hcaldqm::DQModule::_ptype
ProcessingType _ptype
Definition: DQModule.h:59
hcaldqm::Container1D::fill
virtual void fill(uint32_t)
Definition: Container1D.cc:83
hcaldqm::Container1D::reset
virtual void reset()
Definition: Container1D.cc:66
DigiTask::_ledSignalPresent
bool _ledSignalPresent
Definition: DigiTask.h:197
DigiTask::_tagHBHE
edm::InputTag _tagHBHE
Definition: DigiTask.h:46
hcaldqm::flag::Flag::reset
void reset()
Definition: Flag.h:76
MonitorElement::setLumiFlag
void setLumiFlag()
this ME is meant to be stored for each luminosity section
Definition: MonitorElement.h:150
DigiTask::_cOccupancy_CrateSlot
hcaldqm::Container2D _cOccupancy_CrateSlot
Definition: DigiTask.h:133
hcaldqm::ContainerSingle2D::setBinContent
virtual void setBinContent(int, int, int)
Definition: ContainerSingle2D.cc:211
HcalElectronicsId::isVMEid
constexpr bool isVMEid() const
Definition: HcalElectronicsId.h:67
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
DigiTask::_cADCvsTS_SubdetPM_QIE1011
hcaldqm::Container2D _cADCvsTS_SubdetPM_QIE1011
Definition: DigiTask.h:108
DigiTask::_cOccupancyCut_ElectronicsVME
hcaldqm::Container2D _cOccupancyCut_ElectronicsVME
Definition: DigiTask.h:143
DigiTask::_cDigiSize_Crate
hcaldqm::Container1D _cDigiSize_Crate
Definition: DigiTask.h:157
DigiTask::_cShapeCut_FED
hcaldqm::Container1D _cShapeCut_FED
Definition: DigiTask.h:106
DigiTask::_cOccupancyCut_depth
hcaldqm::Container2D _cOccupancyCut_depth
Definition: DigiTask.h:145
DigiTask::_cOccupancyCutvsBX_Subdet
hcaldqm::ContainerProf1D _cOccupancyCutvsBX_Subdet
Definition: DigiTask.h:154
DigiTask::_cutSumQ_HF
double _cutSumQ_HF
Definition: DigiTask.h:55
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:118
bookHistograms
virtual example_stream void bookHistograms(DQMStore::IBooker &,@example_stream edm::Run const &,@example_stream edm::EventSetup const &) override
hcaldqm::constants::FIBER_uTCA_MIN1
int const FIBER_uTCA_MIN1
Definition: Constants.h:144
DigiTask::_cfC_SubdetPM_QIE1011
hcaldqm::Container1D _cfC_SubdetPM_QIE1011
Definition: DigiTask.h:99
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:166
DigiTask::_cTimingCut_depth
hcaldqm::ContainerProf2D _cTimingCut_depth
Definition: DigiTask.h:119
QIE10DataFrame::detid
DetId detid() const
Get the detector id.
Definition: QIE10DataFrame.h:78
HcalElectronicsMap::allPrecisionId
std::vector< HcalGenericDetId > allPrecisionId() const
Definition: HcalElectronicsMap.cc:139
hcaldqm::DQModule::_logger
Logger _logger
Definition: DQModule.h:70
DigiTask::_cBadTDCValues_SubdetPM
hcaldqm::Container1D _cBadTDCValues_SubdetPM
Definition: DigiTask.h:174
DigiTask::_cOccupancy_FEDVME
hcaldqm::Container2D _cOccupancy_FEDVME
Definition: DigiTask.h:128
DigiTask::_cTimingCutvsLS_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsLS_SubdetPM
Definition: DigiTask.h:118
DigiTask::_cOccupancyCut_FEDVME
hcaldqm::Container2D _cOccupancyCut_FEDVME
Definition: DigiTask.h:141
DigiTask::_cSumQvsBX_SubdetPM_QIE1011
hcaldqm::ContainerProf1D _cSumQvsBX_SubdetPM_QIE1011
Definition: DigiTask.h:102
DigiTask::_cCapidMinusBXmod4_CrateSlotVME
hcaldqm::ContainerSingle2D _cCapidMinusBXmod4_CrateSlotVME[4]
Definition: DigiTask.h:186
DigiTask::_cOccupancy_depth
hcaldqm::Container2D _cOccupancy_depth
Definition: DigiTask.h:134
hcaldqm::DQModule::_name
std::string _name
Definition: DQModule.h:57
hcaldqm::DQTask::_emap
HcalElectronicsMap const * _emap
Definition: DQTask.h:73
globals_cff.x1
x1
Definition: globals_cff.py:28
DigiTask::_cSummaryvsLS_FED
hcaldqm::Container2D _cSummaryvsLS_FED
Definition: DigiTask.h:199
DigiTask::_filter_FEDHF
hcaldqm::filter::HashFilter _filter_FEDHF
Definition: DigiTask.h:83
hcaldqm::constants::CRATE_uTCA_MIN
int const CRATE_uTCA_MIN
Definition: Constants.h:116
DigiTask::_tagHE
edm::InputTag _tagHE
Definition: DigiTask.h:47
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:152
DigiTask::_resetMonitors
void _resetMonitors(hcaldqm::UpdateFreq) override
Definition: DigiTask.cc:598
DigiTask::_cSumQvsLS_SubdetPM
hcaldqm::ContainerProf1D _cSumQvsLS_SubdetPM
Definition: DigiTask.h:94
DigiTask::_xUni
hcaldqm::ContainerXXX< uint32_t > _xUni
Definition: DigiTask.h:161
DigiTask::_cTimingCutvsLS_FED
hcaldqm::ContainerProf1D _cTimingCutvsLS_FED
Definition: DigiTask.h:117
globals_cff.x2
x2
Definition: globals_cff.py:29
DigiTask::_cOccupancyCut_FEDuTCA
hcaldqm::Container2D _cOccupancyCut_FEDuTCA
Definition: DigiTask.h:142
DigiTask::_cOccupancy_ElectronicsVME
hcaldqm::Container2D _cOccupancy_ElectronicsVME
Definition: DigiTask.h:130
DigiTask::_tagHF
edm::InputTag _tagHF
Definition: DigiTask.h:49
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:274
hcaldqm::utilities::isFEDHBHE
bool isFEDHBHE(HcalElectronicsId const &)
Definition: Utilities.cc:172
hcaldqm::ContainerXXX::book
virtual void book(HcalElectronicsMap const *)
Definition: ContainerXXX.h:97
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:136
CaloSamples::size
int size() const
get the size
Definition: CaloSamples.h:24
DigiTask::_cDigiSize_FED
hcaldqm::Container1D _cDigiSize_FED
Definition: DigiTask.h:158
DigiTask::_filter_QIE8
hcaldqm::filter::HashFilter _filter_QIE8
Definition: DigiTask.h:85
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:958
hcaldqm::Container2D::initialize
virtual void initialize(std::string const &folder, hashfunctions::HashType, quantity::Quantity *, quantity::Quantity *, quantity::Quantity *qz=new quantity::ValueQuantity(quantity::fN), int debug=0)
Definition: Container2D.cc:32
DigiTask::_xNChsNominal
hcaldqm::ContainerXXX< uint32_t > _xNChsNominal
Definition: DigiTask.h:163
DigiTask::_cSumQ_SubdetPM
hcaldqm::Container1D _cSumQ_SubdetPM
Definition: DigiTask.h:92
DigiTask::_cOccupancyvsiphi_SubdetPM
hcaldqm::Container1D _cOccupancyvsiphi_SubdetPM
Definition: DigiTask.h:135
hcaldqm::ContainerProf2D::book
void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="") override
Definition: ContainerProf2D.cc:51
DigiTask::_cTimingCutvsieta_Subdet
hcaldqm::ContainerProf1D _cTimingCutvsieta_Subdet
Definition: DigiTask.h:121
DigiTask::_cLETDCTime_SubdetPM
hcaldqm::Container1D _cLETDCTime_SubdetPM
Definition: DigiTask.h:170
DigiTask::_dhashmap
hcaldqm::electronicsmap::ElectronicsMap _dhashmap
Definition: DigiTask.h:78
hcaldqm::ContainerProf2D::initialize
void initialize(std::string const &folder, hashfunctions::HashType, quantity::Quantity *, quantity::Quantity *, quantity::Quantity *qz=new quantity::ValueQuantity(quantity::fEnergy), int debug=0) override
Definition: ContainerProf2D.cc:27
DigiTask::_vflags
std::vector< hcaldqm::flag::Flag > _vflags
Definition: DigiTask.h:59
hcaldqm::Container2D::setBinContent
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:144
DigiTask::_cTimingCut_ElectronicsVME
hcaldqm::ContainerProf2D _cTimingCut_ElectronicsVME
Definition: DigiTask.h:115
DigiTask::_cLETDCvsADC_SubdetPM
hcaldqm::Container2D _cLETDCvsADC_SubdetPM
Definition: DigiTask.h:168
hcaldqm::ContainerXXX::get
virtual STDTYPE & get(HcalDetId const &)
Definition: ContainerXXX.h:247
hcaldqm::DQTask::_dbService
edm::ESHandle< HcalDbService > _dbService
Definition: DQTask.h:72
DigiTask::_cBadTDCvsBX_SubdetPM
hcaldqm::Container1D _cBadTDCvsBX_SubdetPM
Definition: DigiTask.h:175
hcaldqm::Container1D::setLumiFlag
virtual void setLumiFlag()
Definition: Container1D.cc:1361
hcaldqm::utilities::getFEDList
std::vector< int > getFEDList(HcalElectronicsMap const *)
Definition: Utilities.cc:95
DigiTask::_cOccupancy_Crate
hcaldqm::Container2D _cOccupancy_Crate
Definition: DigiTask.h:132
edm::EventSetup::get
T get() const
Definition: EventSetup.h:63
hcaldqm::utilities::getFEDuTCAList
std::vector< int > getFEDuTCAList(HcalElectronicsMap const *)
Definition: Utilities.cc:145
hcaldqm::Container2D::fill
void fill(HcalDetId const &) override
Definition: Container2D.cc:59
HcalDbService::getHcalMapping
const HcalElectronicsMap * getHcalMapping() const
Definition: HcalDbService.cc:309
hcaldqm::filter::HashFilter::initialize
virtual void initialize(FilterType ftype, hashfunctions::HashType htype, std::vector< uint32_t > const &)
Definition: HashFilter.cc:28
hcaldqm::DQModule::_subsystem
std::string _subsystem
Definition: DQModule.h:64
DigiTask::_filter_uTCA
hcaldqm::filter::HashFilter _filter_uTCA
Definition: DigiTask.h:82
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:159
DigiTask::_xNChs
hcaldqm::ContainerXXX< uint32_t > _xNChs
Definition: DigiTask.h:162
hcaldqm::ContainerSingle2D::fill
virtual void fill(int, int)
Definition: ContainerSingle2D.cc:133
DigiTask::_ehashmap
hcaldqm::electronicsmap::ElectronicsMap _ehashmap
Definition: DigiTask.h:77
DigiTask::_cOccupancyCutvsLS_Subdet
hcaldqm::ContainerProf1D _cOccupancyCutvsLS_Subdet
Definition: DigiTask.h:153
hcaldqm::ContainerXXX::initialize
virtual void initialize(hashfunctions::HashType, int debug=0)
Definition: ContainerXXX.h:89
DigiTask::_filter_QIE1011
hcaldqm::filter::HashFilter _filter_QIE1011
Definition: DigiTask.h:84
hcaldqm::Container2D::book
void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="") override
Definition: Container2D.cc:895
DigiTask::_cBadTDCCount_depth
hcaldqm::Container2D _cBadTDCCount_depth
Definition: DigiTask.h:177
DigiTask::_cLETDCTime_depth
hcaldqm::ContainerProf2D _cLETDCTime_depth
Definition: DigiTask.h:171
DigiTask::_vhashFEDs
std::vector< uint32_t > _vhashFEDs
Definition: DigiTask.h:72
QIE11DataFrame::samples
int samples() const
total number of samples in the digi
Definition: QIE11DataFrame.h:57
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:973
DigiTask::_cOccupancyvsieta_Subdet
hcaldqm::Container1D _cOccupancyvsieta_Subdet
Definition: DigiTask.h:136
HcalElectronicsId
Readout chain identification for Hcal.
Definition: HcalElectronicsId.h:32
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:37
nanoDQM_cff.Flag
Flag
Definition: nanoDQM_cff.py:62
hcaldqm::utilities::isFEDHF
bool isFEDHF(HcalElectronicsId const &)
Definition: Utilities.cc:187
DigiTask::_cTimingCut_ElectronicsuTCA
hcaldqm::ContainerProf2D _cTimingCut_ElectronicsuTCA
Definition: DigiTask.h:116
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::_cSumQ_SubdetPM_QIE1011
hcaldqm::Container1D _cSumQ_SubdetPM_QIE1011
Definition: DigiTask.h:100
DigiTask::_cTimingCut_SubdetPM
hcaldqm::Container1D _cTimingCut_SubdetPM
Definition: DigiTask.h:112
edm::Run
Definition: Run.h:44
DigiTask::_unknownIdsPresent
bool _unknownIdsPresent
Definition: DigiTask.h:194
DigiTask::_cOccupancy_FEDuTCA
hcaldqm::Container2D _cOccupancy_FEDuTCA
Definition: DigiTask.h:129
DigiTask::_cCapidMinusBXmod4_SubdetPM
hcaldqm::Container1D _cCapidMinusBXmod4_SubdetPM
Definition: DigiTask.h:184
cuy.ib
ib
Definition: cuy.py:661
DigiTask::_xUniHF
hcaldqm::ContainerXXX< uint32_t > _xUniHF
Definition: DigiTask.h:161
hcaldqm::utilities::sumQ
double sumQ(DIGI const &digi, double ped, int i=0, int j=3)
Definition: Utilities.h:121
DigiTask::_cfC_SubdetPM
hcaldqm::Container1D _cfC_SubdetPM
Definition: DigiTask.h:91
DigiTask::_cOccupancy_ElectronicsuTCA
hcaldqm::Container2D _cOccupancy_ElectronicsuTCA
Definition: DigiTask.h:131
Generated for CMSSW Reference Manual by doxygen 1.8.11