CMS 3D CMS Logo
RecHitTask.cc
Go to the documentation of this file.
1 #include "DQM/HcalTasks/interface/RecHitTask.h"
2 #include <math.h>
3 
4 using namespace hcaldqm;
5 using namespace hcaldqm::constants;
6 using namespace hcaldqm::filter;
7 
8 RecHitTask::RecHitTask(edm::ParameterSet const& ps):
9  DQTask(ps)
10 {
11  _tagHBHE = ps.getUntrackedParameter<edm::InputTag>("tagHBHE",
12  edm::InputTag("hbhereco"));
13  _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO",
14  edm::InputTag("horeco"));
15  _tagHF = ps.getUntrackedParameter<edm::InputTag>("tagHF",
16  edm::InputTag("hfreco"));
17  _tagPreHF = ps.getUntrackedParameter<edm::InputTag>("tagPreHF",
18  edm::InputTag(""));
19  _hfPreRecHitsAvailable = ps.getUntrackedParameter<bool>("hfPreRecHitsAvailable", false);
20 
21  _tokHBHE = consumes<HBHERecHitCollection>(_tagHBHE);
22  _tokHO = consumes<HORecHitCollection>(_tagHO);
23  _tokHF = consumes<HFRecHitCollection>(_tagHF);
24  _tokPreHF = consumes<HFPreRecHitCollection>(_tagPreHF);
25 
26  _cutE_HBHE = ps.getUntrackedParameter<double>("cutE_HBHE", 5);
27  _cutE_HO = ps.getUntrackedParameter<double>("cutE_HO", 5);
28  _cutE_HF = ps.getUntrackedParameter<double>("cutE_HF", 5);
29  _thresh_unihf = ps.getUntrackedParameter<double>("thresh_unihf", 0.2);
30 
31  // order must be the same as in RecoFlag enum
32  _vflags.resize(nRecoFlag);
33  _vflags[fUni]=flag::Flag("UniSlotHF");
34  _vflags[fTCDS]=flag::Flag("TCDS");
35  _vflags[fUnknownIds] = flag::Flag("UnknownIds");
36 }
37 
38 /* virtual */ void RecHitTask::bookHistograms(DQMStore::IBooker& ib,
39  edm::Run const& r, edm::EventSetup const& es)
40 {
41  DQTask::bookHistograms(ib,r,es);
42 
43  // GET WHAT YOU NEED
44  edm::ESHandle<HcalDbService> dbs;
45  es.get<HcalDbRecord>().get(dbs);
46  _emap = dbs->getHcalMapping();
47 
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  std::vector<uint32_t> vhashHF;
59  vhashHF.push_back(hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdet](HcalDetId(HcalForward, 29,1,1)));
60  _filter_HF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fSubdet, vhashHF);
61 
62  // INITIALIZE FIRST
63  // Energy
64  _cEnergy_Subdet.initialize(_name, "Energy", hcaldqm::hashfunctions::fSubdet,
65  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy),
66  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
67  _cEnergy_depth.initialize(_name, "Energy", hcaldqm::hashfunctions::fdepth,
68  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
69  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
70  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy, true),0);
71 
72  // Timing
73  _cTimingCut_SubdetPM.initialize(_name, "TimingCut",
74  hcaldqm::hashfunctions::fSubdetPM,
75  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
76  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
77  _cTimingvsEnergy_SubdetPM.initialize(_name, "TimingvsEnergy",
78  hcaldqm::hashfunctions::fSubdetPM,
79  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy, true),
80  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
81  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
82 
83  _cTimingCut_HBHEPartition.initialize(_name, "TimingCut",
84  hcaldqm::hashfunctions::fHBHEPartition,
85  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
86  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
87  _cTimingCut_depth.initialize(_name, "TimingCut",
88  hcaldqm::hashfunctions::fdepth,
89  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
90  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
91  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
92  _cTimingCutvsLS_SubdetPM.initialize(_name, "TimingCutvsLS",
93  hcaldqm::hashfunctions::fSubdetPM,
94  new hcaldqm::quantity::LumiSection(_maxLS),
95  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
96 
97  // Occupancy
98  _cOccupancy_depth.initialize(_name, "Occupancy",
99  hcaldqm::hashfunctions::fdepth,
100  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
101  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
102  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
103 
104  _cOccupancyvsLS_Subdet.initialize(_name, "OccupancyvsLS",
105  hcaldqm::hashfunctions::fSubdet,
106  new hcaldqm::quantity::LumiSection(_maxLS),
107  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),0);
108 
109  _cOccupancyCut_depth.initialize(_name, "OccupancyCut",
110  hcaldqm::hashfunctions::fdepth,
111  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
112  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
113  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
114 
115  if (_hfPreRecHitsAvailable) {
116  _cDAAsymmetryVsCharge_SubdetPM.initialize(_name, "ChargeVsAsymmetry",
117  hcaldqm::hashfunctions::fSubdetPM,
118  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDualAnodeAsymmetry),
119  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000),
120  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
121  _cDAAsymmetryMean_cut_depth.initialize(_name, "AsymmetryMean",
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::fDualAnodeAsymmetry),0);
126  _cDAAsymmetry_cut_SubdetPM.initialize(_name, "Asymmetry",
127  hcaldqm::hashfunctions::fSubdetPM,
128  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDualAnodeAsymmetry),
129  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN), 0);
130  }
131 
132  // INITIALIZE HISTOGRAMS to be used only in Online
133  if (_ptype==fOnline)
134  {
135  _cEnergyvsieta_Subdet.initialize(_name, "Energyvsieta",
136  hcaldqm::hashfunctions::fSubdet,
137  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
138  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),0);
139  _cEnergyvsiphi_SubdetPM.initialize(_name, "Energyvsiphi",
140  hcaldqm::hashfunctions::fSubdetPM,
141  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
142  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),0);
143  _cEnergyvsLS_SubdetPM.initialize(_name, "EnergyvsLS",
144  hcaldqm::hashfunctions::fSubdetPM,
145  new hcaldqm::quantity::LumiSection(_maxLS),
146  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),0);
147  _cEnergyvsBX_SubdetPM.initialize(_name, "EnergyvsBX",
148  hcaldqm::hashfunctions::fSubdetPM,
149  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
150  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),0);
151  _cTimingCutvsieta_Subdet.initialize(_name, "TimingCutvsieta",
152  hcaldqm::hashfunctions::fSubdet,
153  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
154  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
155  _cTimingCutvsiphi_SubdetPM.initialize(_name, "TimingCutvsiphi",
156  hcaldqm::hashfunctions::fSubdetPM,
157  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
158  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
159  _cTimingCutvsBX_SubdetPM.initialize(_name, "TimingCutvsBX",
160  hcaldqm::hashfunctions::fSubdetPM,
161  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
162  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
163  _cOccupancyvsiphi_SubdetPM.initialize(_name, "Occupancyvsiphi",
164  hcaldqm::hashfunctions::fSubdetPM,
165  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
166  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
167  _cOccupancyvsieta_Subdet.initialize(_name, "Occupancyvsieta",
168  hcaldqm::hashfunctions::fSubdet,
169  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
170  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
171  _cOccupancyCutvsiphi_SubdetPM.initialize(_name, "OccupancyCutvsiphi",
172  hcaldqm::hashfunctions::fSubdetPM,
173  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
174  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
175  _cOccupancyCutvsieta_Subdet.initialize(_name, "OccupancyCutvsieta",
176  hcaldqm::hashfunctions::fSubdet,
177  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
178  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
179  _cOccupancyCutvsBX_Subdet.initialize(_name, "OccupancyCutvsBX",
180  hcaldqm::hashfunctions::fSubdet,
181  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
182  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
183  _cOccupancyCutvsiphivsLS_SubdetPM.initialize(_name,
184  "OccupancyCutvsiphivsLS", hcaldqm::hashfunctions::fSubdetPM,
185  new hcaldqm::quantity::LumiSection(_maxLS),
186  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
187  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
188  _cOccupancyCutvsLS_Subdet.initialize(_name, "OccupancyCutvsLS",
189  hcaldqm::hashfunctions::fSubdet,
190  new hcaldqm::quantity::LumiSection(_maxLS),
191  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),0);
192  }
193 
194  // FED-based plots
195  if (_ptype != fOffline) { // hidefed2crate
196  std::vector<int> vFEDs = hcaldqm::utilities::getFEDList(_emap);
197  std::vector<int> vFEDsVME = hcaldqm::utilities::getFEDVMEList(_emap);
198  std::vector<int> vFEDsuTCA = hcaldqm::utilities::getFEDuTCAList(_emap);
199 
200  // push the rawIds of each fed into the vector
201  for (std::vector<int>::const_iterator it=vFEDsVME.begin();
202  it!=vFEDsVME.end(); ++it)
203  _vhashFEDs.push_back(HcalElectronicsId(
204  FIBERCH_MIN, FIBER_VME_MIN, SPIGOT_MIN, (*it)-FED_VME_MIN).rawId());
205  for (std::vector<int>::const_iterator it=vFEDsuTCA.begin();
206  it!=vFEDsuTCA.end(); ++it)
207  {
208  std::pair<uint16_t, uint16_t> cspair = utilities::fed2crate(*it);
209  _vhashFEDs.push_back(HcalElectronicsId(
210  cspair.first, cspair.second, FIBER_uTCA_MIN1,
211  FIBERCH_MIN, false).rawId());
212  }
213 
214  _cTimingCut_FEDVME.initialize(_name, "TimingCut",
215  hcaldqm::hashfunctions::fFED,
216  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
217  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
218  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
219  _cTimingCut_FEDuTCA.initialize(_name, "TimingCut",
220  hcaldqm::hashfunctions::fFED,
221  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
222  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
223  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
224  _cTimingCutvsLS_FED.initialize(_name, "TimingCutvsLS",
225  hcaldqm::hashfunctions::fFED,
226  new hcaldqm::quantity::LumiSection(_maxLS),
227  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
228  _cTimingCut_ElectronicsVME.initialize(_name, "TimingCut",
229  hcaldqm::hashfunctions::fElectronics,
230  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
231  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
232  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
233  _cTimingCut_ElectronicsuTCA.initialize(_name, "TimingCut",
234  hcaldqm::hashfunctions::fElectronics,
235  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
236  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
237  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
238 
239  _cOccupancy_FEDVME.initialize(_name, "Occupancy",
240  hcaldqm::hashfunctions::fFED,
241  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
242  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
243  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
244  _cOccupancy_FEDuTCA.initialize(_name, "Occupancy",
245  hcaldqm::hashfunctions::fFED,
246  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
247  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
248  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
249  _cOccupancy_ElectronicsVME.initialize(_name, "Occupancy",
250  hcaldqm::hashfunctions::fElectronics,
251  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
252  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
253  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
254  _cOccupancy_ElectronicsuTCA.initialize(_name, "Occupancy",
255  hcaldqm::hashfunctions::fElectronics,
256  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
257  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
258  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
259 
260  _cOccupancyCut_FEDVME.initialize(_name, "OccupancyCut",
261  hcaldqm::hashfunctions::fFED,
262  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
263  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
264  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
265  _cOccupancyCut_FEDuTCA.initialize(_name, "OccupancyCut",
266  hcaldqm::hashfunctions::fFED,
267  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
268  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
269  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
270  _cOccupancyCut_ElectronicsVME.initialize(_name, "OccupancyCut",
271  hcaldqm::hashfunctions::fElectronics,
272  new hcaldqm::quantity::FEDQuantity(vFEDsVME),
273  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
274  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
275  _cOccupancyCut_ElectronicsuTCA.initialize(_name, "OccupancyCut",
276  hcaldqm::hashfunctions::fElectronics,
277  new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
278  new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
279  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
280  if (_ptype == fOnline) {
281  _cSummaryvsLS_FED.initialize(_name, "SummaryvsLS",
282  hcaldqm::hashfunctions::fFED,
283  new hcaldqm::quantity::LumiSection(_maxLS),
284  new hcaldqm::quantity::FlagQuantity(_vflags),
285  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),0);
286  _cSummaryvsLS.initialize(_name, "SummaryvsLS",
287  new hcaldqm::quantity::LumiSection(_maxLS),
288  new hcaldqm::quantity::FEDQuantity(vFEDs),
289  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),0);
290 
291  _xUniHF.initialize(hcaldqm::hashfunctions::fFEDSlot);
292  _xUni.initialize(hcaldqm::hashfunctions::fFED);
293  }
294  }
295 
296  // BOOK HISTOGRAMS
297  char cutstr[200];
298  sprintf(cutstr, "_EHBHE%dHO%dHF%d", int(_cutE_HBHE),
299  int(_cutE_HO), int(_cutE_HF));
300  char cutstr2[200];
301  sprintf(cutstr2, "_EHF%d", int(_cutE_HF));
302 
303  // Energy
304  _cEnergy_Subdet.book(ib, _emap, _subsystem);
305  _cEnergy_depth.book(ib, _emap, _subsystem);
306 
307  // Timing
308  _cTimingCut_SubdetPM.book(ib, _emap, _subsystem);
309  _cTimingvsEnergy_SubdetPM.book(ib, _emap, _subsystem);
310  if (_ptype != fOffline) { // hidefed2crate
311  _cTimingCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
312  _cTimingCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
313  _cTimingCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
314  _cTimingCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
315  }
316  _cTimingCut_HBHEPartition.book(ib, _emap, _subsystem);
317  _cTimingCut_depth.book(ib, _emap, _subsystem);
318  _cTimingCutvsLS_SubdetPM.book(ib, _emap, _subsystem);
319  if (_ptype != fOffline) { // hidefed2crate
320  _cTimingCutvsLS_FED.book(ib, _emap, _subsystem);
321  }
322 
323  // Occupancy
324  _cOccupancy_depth.book(ib, _emap, _subsystem);
325  if (_ptype != fOffline) { // hidefed2crate
326  _cOccupancy_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
327  _cOccupancy_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
328  _cOccupancy_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
329  _cOccupancy_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
330  }
331  _cOccupancyvsLS_Subdet.book(ib, _emap, _subsystem);
332  _cOccupancyCut_depth.book(ib, _emap, _subsystem);
333  if (_ptype != fOffline) { // hidefed2crate
334  _cOccupancyCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
335  _cOccupancyCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
336  _cOccupancyCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
337  _cOccupancyCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
338  }
339 
340  if (_hfPreRecHitsAvailable) {
341  _cDAAsymmetryVsCharge_SubdetPM.book(ib, _emap, _filter_HF, _subsystem);
342  _cDAAsymmetryMean_cut_depth.book(ib, _emap, _filter_HF, _subsystem);
343  _cDAAsymmetry_cut_SubdetPM.book(ib, _emap, _filter_HF, _subsystem);
344  }
345 
346  // BOOK HISTOGRAMS to be used only in Online
347  if (_ptype==fOnline)
348  {
349  _cEnergyvsLS_SubdetPM.book(ib, _emap, _subsystem);
350  _cEnergyvsieta_Subdet.book(ib, _emap, _subsystem);
351  _cEnergyvsiphi_SubdetPM.book(ib, _emap, _subsystem);
352  _cEnergyvsBX_SubdetPM.book(ib, _emap, _subsystem);
353  _cTimingCutvsieta_Subdet.book(ib, _emap, _subsystem);
354  _cTimingCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
355  _cTimingCutvsBX_SubdetPM.book(ib, _emap, _subsystem);
356  _cOccupancyvsiphi_SubdetPM.book(ib, _emap, _subsystem);
357  _cOccupancyvsieta_Subdet.book(ib, _emap, _subsystem);
358  _cOccupancyCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
359  _cOccupancyCutvsieta_Subdet.book(ib, _emap, _subsystem);
360  _cOccupancyCutvsBX_Subdet.book(ib, _emap, _subsystem);
361  _cOccupancyCutvsiphivsLS_SubdetPM.book(ib, _emap, _subsystem);
362  _cOccupancyCutvsLS_Subdet.book(ib, _emap, _subsystem);
363  _cSummaryvsLS_FED.book(ib, _emap, _subsystem);
364  _cSummaryvsLS.book(ib, _subsystem);
365 
366  std::vector<uint32_t> vhashFEDHF;
367  vhashFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN,
368  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
369  vhashFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN,
370  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
371  vhashFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN,
372  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
373  vhashFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN+6,
374  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
375  vhashFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN+6,
376  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
377  vhashFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN+6,
378  FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
379  HashFilter filter_FEDHF;
380  filter_FEDHF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fFED,
381  vhashFEDHF);
382 
383  _gids = _emap->allPrecisionId();
384  if (_ptype != fOffline) { // hidefed2crate
385  _xUniHF.book(_emap, filter_FEDHF);
386  _xUni.book(_emap);
387  }
388  }
389 
390  // initialize hash map
391  _ehashmap.initialize(_emap, hcaldqm::electronicsmap::fD2EHashMap);
392 
393  // book some mes...
394  ib.setCurrentFolder(_subsystem+"/"+_name);
395  meUnknownIds1LS = ib.book1D("UnknownIds", "UnknownIds",
396  1, 0, 1);
397  _unknownIdsPresent = false;
398  meUnknownIds1LS->setLumiFlag();
399 }
400 
401 /* virtual */ void RecHitTask::_resetMonitors(hcaldqm::UpdateFreq uf)
402 {
403  switch(uf)
404  {
405  case hcaldqm::f1LS:
406  _unknownIdsPresent = false;
407  break;
408  default:
409  break;
410  }
411 
412  DQTask::_resetMonitors(uf);
413 }
414 
415 /* virtual */ void RecHitTask::_process(edm::Event const& e,
416  edm::EventSetup const&)
417 {
418  edm::Handle<HBHERecHitCollection> chbhe;
419  edm::Handle<HORecHitCollection> cho;
420  edm::Handle<HFRecHitCollection> chf;
421 
422  if (!(e.getByToken(_tokHBHE, chbhe)))
423  _logger.dqmthrow("Collection HBHERecHitCollection not available "
424  + _tagHBHE.label() + " " + _tagHBHE.instance());
425  if (!(e.getByToken(_tokHO, cho)))
426  _logger.dqmthrow("Collection HORecHitCollection not available "
427  + _tagHO.label() + " " + _tagHO.instance());
428  if (!(e.getByToken(_tokHF, chf)))
429  _logger.dqmthrow("Collection HFRecHitCollection not available "
430  + _tagHF.label() + " " + _tagHF.instance());
431 
432  edm::Handle<HFPreRecHitCollection> cprehf;
433  if (_hfPreRecHitsAvailable) {
434  if (!(e.getByToken(_tokPreHF, cprehf)))
435  _logger.dqmthrow("Collection HFPreRecHitCollection not available "
436  + _tagPreHF.label() + " " + _tagPreHF.instance());
437  }
438 
439  // extract some info per event
440  int bx = e.bunchCrossing();
441 
442  // To fill histograms outside of the loop, you need to determine if there were
443  // any valid det ids first
444  uint32_t rawidValid = 0;
445  uint32_t rawidHBValid = 0;
446  uint32_t rawidHEValid = 0;
447 
448  double ehbm = 0; double ehbp = 0;
449  double ehem = 0; double ehep = 0;
450  int nChsHB = 0; int nChsHE = 0;
451  int nChsHBCut = 0; int nChsHECut = 0;
452  for (HBHERecHitCollection::const_iterator it=chbhe->begin();
453  it!=chbhe->end(); ++it)
454  {
455  double energy = it->energy();
456  double timing = it->time();
457 
458  // Explicit check on the DetIds present in the Collection
459  HcalDetId did = it->id();
460  uint32_t rawid = _ehashmap.lookup(did);
461  /*
462  * Needs to be removed as DetIds that belong to the HEP17 after combination
463  * are not present in the emap
464  * Removed until further notice!
465  *
466  */
467  //if (rawid==0)
468  //{meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
469 
470  HcalElectronicsId const& eid(rawid);
471  rawidValid = did.rawId();
472  if (did.subdet()==HcalBarrel)
473  rawidHBValid = did.rawId();
474  else if (did.subdet()==HcalEndcap)
475  rawidHEValid = did.rawId();
476 
477  _cEnergy_Subdet.fill(did, energy);
478  _cTimingvsEnergy_SubdetPM.fill(did, energy, timing);
479  _cOccupancy_depth.fill(did);
480  did.subdet()==HcalBarrel?did.ieta()>0?ehbp+=energy:ehbm+=energy:
481  did.ieta()>0?ehep+=energy:ehem+=energy;
482 
483  // ONLINE ONLY!
484  if (_ptype==fOnline)
485  {
486  _cOccupancyvsiphi_SubdetPM.fill(did);
487  _cOccupancyvsieta_Subdet.fill(did);
488  }
489  // ^^^ONLINE ONLY!
490  //
491  if (_ptype != fOffline) { // hidefed2crate
492  // Also, for these electronics plots, require that the channel was found in the emap.
493  if (rawid != 0) {
494  if (eid.isVMEid())
495  {
496  _cOccupancy_FEDVME.fill(eid);
497  _cOccupancy_ElectronicsVME.fill(eid);
498  } else if (eid.isUTCAid()) {
499  _cOccupancy_FEDuTCA.fill(eid);
500  _cOccupancy_ElectronicsuTCA.fill(eid);
501  }
502  }
503  }
504 
505  if (energy>_cutE_HBHE)
506  {
507  // ONLINE ONLY!
508  if (_ptype==fOnline)
509  {
510  _cEnergyvsLS_SubdetPM.fill(did, _currentLS, energy);
511  _cEnergyvsBX_SubdetPM.fill(did, bx, energy);
512  _cEnergyvsieta_Subdet.fill(did, energy);
513  _cEnergyvsiphi_SubdetPM.fill(did, energy);
514  _cTimingCutvsieta_Subdet.fill(did, timing);
515  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
516  _cTimingCutvsBX_SubdetPM.fill(did, bx, timing);
517  _cOccupancyCutvsiphi_SubdetPM.fill(did);
518  _cOccupancyCutvsieta_Subdet.fill(did);
519  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
520  }
521  // ^^^ONLINE ONLY!
522  _cEnergy_depth.fill(did, energy);
523  _cTimingCut_SubdetPM.fill(did, timing);
524  _cTimingCut_HBHEPartition.fill(did, timing);
525  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
526 
527  // ONLINE
528  if (_ptype==fOnline) {
529  if (rawid != 0) {
530  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
531  }
532  _cTimingCut_depth.fill(did, timing);
533  // ^^^ONLINE
534  } else {
535  if (_ptype != fOffline) { // hidefed2crate
536  if (rawid != 0) {
537  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
538  }
539  }
540  _cTimingCut_depth.fill(did, timing);
541  }
542  _cOccupancyCut_depth.fill(did);
543  if (_ptype != fOffline) { // hidefed2crate
544  if (rawid != 0) {
545  if (eid.isVMEid()){
546  // ONLINE
547  if (_ptype==fOnline)
548  {
549  _cTimingCut_FEDVME.fill(eid, timing);
550  _cTimingCut_ElectronicsVME.fill(eid, timing);
551  } // ^^^ ONLINE
552  else
553  {
554  _cTimingCut_FEDVME.fill(eid, timing);
555  _cTimingCut_ElectronicsVME.fill(eid, timing);
556  }
557  // ^^^ONLINE
558 
559  _cOccupancyCut_FEDVME.fill(eid);
560  _cOccupancyCut_ElectronicsVME.fill(eid);
561  } else if (eid.isUTCAid()) {
562  if (_ptype==fOnline)
563  {
564  // time constraints are explicit!
565  _cTimingCut_FEDuTCA.fill(eid, timing);
566  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
567  }
568  else
569  {
570  _cTimingCut_FEDuTCA.fill(eid, timing);
571  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
572  }
573  _cOccupancyCut_FEDuTCA.fill(eid);
574  _cOccupancyCut_ElectronicsuTCA.fill(eid);
575  }
576  }
577  }
578  did.subdet()==HcalBarrel?nChsHBCut++:nChsHECut++;
579  }
580  did.subdet()==HcalBarrel?nChsHB++:nChsHE++;
581  }
582 
583  if (rawidHBValid!=0 && rawidHEValid!=0)
584  {
585  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHBValid), _currentLS,
586  nChsHB);
587  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHEValid), _currentLS,
588  nChsHE);
589 
590  // ONLINE ONLY!
591  if (_ptype==fOnline)
592  {
593  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHBValid),
594  bx, nChsHBCut);
595  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHEValid),
596  bx, nChsHECut);
597  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHBValid),
598  _currentLS, nChsHBCut);
599  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHEValid),
600  _currentLS, nChsHECut);
601  }
602  // ^^^ONLINE ONLY!
603  }
604 
605  // reset
606  rawidValid = 0;
607 
608  int nChsHO = 0; int nChsHOCut = 0;
609  double ehop = 0; double ehom = 0;
610  for (HORecHitCollection::const_iterator it=cho->begin();
611  it!=cho->end(); ++it)
612  {
613  double energy = it->energy();
614  double timing = it->time();
615 
616  // Explicit check on the DetIds present in the Collection
617  HcalDetId did = it->id();
618  uint32_t rawid = _ehashmap.lookup(did);
619  if (rawid==0)
620  {meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
621  HcalElectronicsId const& eid(rawid);
622  if (did.subdet()==HcalOuter)
623  rawidValid = did.rawId();
624 
625  _cEnergy_Subdet.fill(did, energy);
626  _cTimingvsEnergy_SubdetPM.fill(did, energy, timing);
627  _cOccupancy_depth.fill(did);
628  did.ieta()>0?ehop+=energy:ehom+=energy;
629 
630  // IMPORTANT: ONLINE ONLY!
631  if (_ptype==fOnline)
632  {
633  _cOccupancyvsiphi_SubdetPM.fill(did);
634  _cOccupancyvsieta_Subdet.fill(did);
635  }
636  // ONLINE ONLY!
637 
638  if (_ptype != fOffline) { // hidefed2crate
639  if (eid.isVMEid())
640  {
641  _cOccupancy_FEDVME.fill(eid);
642  _cOccupancy_ElectronicsVME.fill(eid);
643  }
644  else
645  {
646  _cOccupancy_FEDuTCA.fill(eid);
647  _cOccupancy_ElectronicsuTCA.fill(eid);
648  }
649  }
650 
651  if (energy>_cutE_HO)
652  {
653  // ONLINE ONLY!
654  if (_ptype==fOnline)
655  {
656  _cEnergyvsLS_SubdetPM.fill(did, _currentLS, energy);
657  _cEnergyvsBX_SubdetPM.fill(did, bx, energy);
658  _cEnergyvsieta_Subdet.fill(did, energy);
659  _cEnergyvsiphi_SubdetPM.fill(did, energy);
660  _cTimingCutvsieta_Subdet.fill(did, timing);
661  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
662  _cTimingCutvsBX_SubdetPM.fill(did, bx, timing);
663  _cOccupancyCutvsiphi_SubdetPM.fill(did);
664  _cOccupancyCutvsieta_Subdet.fill(did);
665  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
666  }
667  // ^^^ONLINE ONLY!
668 
669  _cEnergy_depth.fill(did, energy);
670  _cTimingCut_SubdetPM.fill(did, timing);
671  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
672  if (_ptype != fOffline) { // hidefed2crate
673  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
674  }
675  _cOccupancyCut_depth.fill(did);
676  _cTimingCut_depth.fill(did, timing);
677  if (_ptype != fOffline) { // hidefed2crate
678  if (eid.isVMEid())
679  {
680  _cTimingCut_FEDVME.fill(eid, timing);
681  _cTimingCut_ElectronicsVME.fill(eid, timing);
682  _cOccupancyCut_FEDVME.fill(eid);
683  _cOccupancyCut_ElectronicsVME.fill(eid);
684  }
685  else
686  {
687  _cTimingCut_FEDuTCA.fill(eid, timing);
688  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
689  _cOccupancyCut_FEDuTCA.fill(eid);
690  _cOccupancyCut_ElectronicsuTCA.fill(eid);
691  }
692  }
693  nChsHOCut++;
694  }
695  nChsHO++;
696  }
697 
698  if (rawidValid!=0)
699  {
700  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS,
701  nChsHO);
702  // ONLINE ONLY!
703  if (_ptype==fOnline)
704  {
705  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid),
706  bx, nChsHOCut);
707  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid),
708  _currentLS, nChsHOCut);
709  }
710  // ^^^ONLINE ONLY!
711  }
712 
713  //reset
714  rawidValid = 0;
715 
716  int nChsHF = 0; int nChsHFCut = 0;
717  double ehfp = 0; double ehfm = 0;
718  for (HFRecHitCollection::const_iterator it=chf->begin();
719  it!=chf->end(); ++it)
720  {
721  double energy = it->energy();
722  double timing = it->time();
723 
724  // Explicit check on the DetIds present in the Collection
725  HcalDetId did = it->id();
726  uint32_t rawid = _ehashmap.lookup(did);
727  if (rawid==0)
728  {meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
729  HcalElectronicsId const& eid(rawid);
730  if (did.subdet()==HcalForward)
731  rawidValid = did.rawId();
732 
733  _cEnergy_Subdet.fill(did, energy);
734  _cTimingvsEnergy_SubdetPM.fill(did, energy, timing);
735  _cOccupancy_depth.fill(did);
736  did.ieta()>0?ehfp+=energy:ehfm+=energy;
737 
738  // IMPORTANT:
739  // only for Online Processing
740  //
741  if (_ptype==fOnline)
742  {
743  _cOccupancyvsiphi_SubdetPM.fill(did);
744  _cOccupancyvsieta_Subdet.fill(did);
745  }
746  // ONLINE ONLY!
747 
748  if (_ptype != fOffline) { // hidefed2crate
749  if (eid.isVMEid())
750  {
751  _cOccupancy_FEDVME.fill(eid);
752  _cOccupancy_ElectronicsVME.fill(eid);
753  }
754  else
755  {
756  _cOccupancy_FEDuTCA.fill(eid);
757  _cOccupancy_ElectronicsuTCA.fill(eid);
758  }
759  }
760 
761  if (energy>_cutE_HF)
762  {
763  // ONLINE ONLY!
764  if (_ptype==fOnline)
765  {
766  _cEnergyvsLS_SubdetPM.fill(did, _currentLS, energy);
767  _cEnergyvsBX_SubdetPM.fill(did, bx, energy);
768  _cEnergyvsieta_Subdet.fill(did, energy);
769  _cEnergyvsiphi_SubdetPM.fill(did, energy);
770  _cTimingCutvsieta_Subdet.fill(did, timing);
771  _cTimingCutvsiphi_SubdetPM.fill(did, timing);
772  _cTimingCutvsBX_SubdetPM.fill(did, bx, timing);
773  _cOccupancyCutvsiphi_SubdetPM.fill(did);
774  _cOccupancyCutvsieta_Subdet.fill(did);
775  _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
776  if (_ptype != fOffline) { // hidefed2crate
777  _xUniHF.get(eid)++;
778  }
779  }
780  // ^^^ONLINE ONLY!
781  _cEnergy_depth.fill(did, energy);
782  _cTimingCut_SubdetPM.fill(did, timing);
783  _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
784  if (_ptype != fOffline) { // hidefed2crate
785  _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
786  }
787  _cOccupancyCut_depth.fill(did);
788  _cTimingCut_depth.fill(did, timing);
789  if (_ptype != fOffline) { // hidefed2crate
790  if (eid.isVMEid())
791  {
792  _cTimingCut_FEDVME.fill(eid, timing);
793  _cTimingCut_ElectronicsVME.fill(eid, timing);
794  _cOccupancyCut_FEDVME.fill(eid);
795  _cOccupancyCut_ElectronicsVME.fill(eid);
796  }
797  else
798  {
799  _cTimingCut_FEDuTCA.fill(eid, timing);
800  _cTimingCut_ElectronicsuTCA.fill(eid, timing);
801  _cOccupancyCut_FEDuTCA.fill(eid);
802  _cOccupancyCut_ElectronicsuTCA.fill(eid);
803  }
804  }
805  nChsHFCut++;
806  }
807  nChsHF++;
808  }
809 
810  if (rawidValid!=0)
811  {
812  _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS,
813  nChsHF);
814  // ONLINE ONLY!
815  if (_ptype==fOnline)
816  {
817  _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid),
818  bx, nChsHFCut);
819  _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid),
820  _currentLS, nChsHFCut);
821  }
822  // ^^^ONLINE ONLY!
823  }
824 
825  // Loop over HFPreRecHits to get charge and charge asymmetry
826  if (_hfPreRecHitsAvailable) {
827  for (HFPreRecHitCollection::const_iterator it=cprehf->begin();
828  it!=cprehf->end(); ++it)
829  {
830  HcalDetId did = it->id();
831  if (_filter_HF.filter(did)) {
832  continue;
833  }
834  std::pair<float, bool> chargeAsymmetry = it->chargeAsymmetry(0.);
835  std::pair<float, bool> chargeAsymmetryCut = it->chargeAsymmetry(20.);
836 
837  if (chargeAsymmetry.second) {
838  _cDAAsymmetryVsCharge_SubdetPM.fill(did, chargeAsymmetry.first, it->charge());
839  }
840  if (chargeAsymmetryCut.second) {
841  _cDAAsymmetryMean_cut_depth.fill(did, chargeAsymmetryCut.first);
842  _cDAAsymmetry_cut_SubdetPM.fill(did, chargeAsymmetryCut.first);
843  }
844  }
845  }
846 }
847 
848 /* virtual */ void RecHitTask::beginLuminosityBlock(edm::LuminosityBlock const&
849  lb, edm::EventSetup const& es)
850 {
851  DQTask::beginLuminosityBlock(lb, es);
852 }
853 
854 /* virtual */ void RecHitTask::endLuminosityBlock(edm::LuminosityBlock const&
855  lb, edm::EventSetup const& es)
856 {
857  if (_ptype!=fOnline)
858  return;
859 
860  //
861  // GENERATE STATUS ONLY FOR ONLINE
862  //
863 // for (std::vector<HcalGenericDetId>::const_iterator it=gids.begin();
864 // it!=gids.end(); ++it)
865 // {}
866 
867  for (uintCompactMap::const_iterator it=_xUniHF.begin();
868  it!=_xUniHF.end(); ++it)
869  {
870  uint32_t hash1 = it->first;
871  HcalElectronicsId eid1(hash1);
872  double x1 = it->second;
873 
874  for (uintCompactMap::const_iterator jt=_xUniHF.begin();
875  jt!=_xUniHF.end(); ++jt)
876  {
877  if (jt==it)
878  continue;
879  double x2 = jt->second;
880  if (x2==0)
881  continue;
882  if (x1/x2<_thresh_unihf)
883  _xUni.get(eid1)++;
884  }
885  }
886 
887  if (_ptype != fOffline) { // hidefed2crate
888  for (std::vector<uint32_t>::const_iterator it=_vhashFEDs.begin();
889  it!=_vhashFEDs.end(); ++it)
890  {
891  flag::Flag fSum("RECO");
892  HcalElectronicsId eid = HcalElectronicsId(*it);
893 
894  std::vector<uint32_t>::const_iterator cit=std::find(
895  _vcdaqEids.begin(), _vcdaqEids.end(), *it);
896  if (cit==_vcdaqEids.end())
897  {
898  // not @cDAQ
899  for (uint32_t iflag=0; iflag<_vflags.size(); iflag++)
900  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag),
901  int(flag::fNCDAQ));
902  _cSummaryvsLS.setBinContent(eid, _currentLS, int(flag::fNCDAQ));
903  continue;
904  }
905 
906  // FED is @cDAQ
907  if (hcaldqm::utilities::isFEDHF(eid) && (_runkeyVal==0 || _runkeyVal==4))
908  {
909  if (_xUni.get(eid)>0)
910  _vflags[fUni]._state = flag::fPROBLEMATIC;
911  else
912  _vflags[fUni]._state = flag::fGOOD;
913  }
914 
915  if (_unknownIdsPresent)
916  _vflags[fUnknownIds]._state = flag::fBAD;
917  else
918  _vflags[fUnknownIds]._state = flag::fGOOD;
919 
920  int iflag=0;
921  for (std::vector<flag::Flag>::iterator ft=_vflags.begin();
922  ft!=_vflags.end(); ++ft)
923  {
924  _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag),
925  int(ft->_state));
926  fSum+=(*ft);
927  iflag++;
928 
929  // reset after using
930  ft->reset();
931  }
932  _cSummaryvsLS.setBinContent(eid, _currentLS, fSum._state);
933  }
934  _xUniHF.reset(); _xUni.reset();
935  }
936 
937  // in the end always do the DQTask::endLumi
938  DQTask::endLuminosityBlock(lb, es);
939 }
940 
941 DEFINE_FWK_MODULE(RecHitTask);
942 
RecHitTask::_xUniHF
hcaldqm::ContainerXXX< uint32_t > _xUniHF
Definition: RecHitTask.h:125
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
RecHitTask::_cOccupancyCut_depth
hcaldqm::Container2D _cOccupancyCut_depth
Definition: RecHitTask.h:120
RecHitTask::_cSummaryvsLS_FED
hcaldqm::Container2D _cSummaryvsLS_FED
Definition: RecHitTask.h:137
RecHitTask::_vhashFEDs
std::vector< uint32_t > _vhashFEDs
Definition: RecHitTask.h:57
RecHitTask::_cOccupancyvsLS_Subdet
hcaldqm::ContainerProf1D _cOccupancyvsLS_Subdet
Definition: RecHitTask.h:110
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
RecHitTask::_tokHBHE
edm::EDGetTokenT< HBHERecHitCollection > _tokHBHE
Definition: RecHitTask.h:48
RecHitTask::_cEnergy_depth
hcaldqm::ContainerProf2D _cEnergy_depth
Definition: RecHitTask.h:83
RecHitTask::_cTimingvsEnergy_SubdetPM
hcaldqm::Container2D _cTimingvsEnergy_SubdetPM
Definition: RecHitTask.h:88
hcaldqm::constants::CRATE_VME_MIN
int const CRATE_VME_MIN
Definition: Constants.h:111
RecHitTask::_cOccupancyCut_FEDVME
hcaldqm::Container2D _cOccupancyCut_FEDVME
Definition: RecHitTask.h:115
RecHitTask::_cOccupancy_FEDVME
hcaldqm::Container2D _cOccupancy_FEDVME
Definition: RecHitTask.h:106
HcalDetId::subdet
HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:49
RecHitTask::_cutE_HO
double _cutE_HO
Definition: RecHitTask.h:53
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
RecHitTask::_tokHF
edm::EDGetTokenT< HFRecHitCollection > _tokHF
Definition: RecHitTask.h:50
RecHitTask::_xUni
hcaldqm::ContainerXXX< uint32_t > _xUni
Definition: RecHitTask.h:125
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
RecHitTask::_cOccupancyCutvsiphi_SubdetPM
hcaldqm::Container1D _cOccupancyCutvsiphi_SubdetPM
Definition: RecHitTask.h:121
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
RecHitTask::_cEnergyvsieta_Subdet
hcaldqm::ContainerProf1D _cEnergyvsieta_Subdet
Definition: RecHitTask.h:81
RecHitTask::_cEnergyvsiphi_SubdetPM
hcaldqm::ContainerProf1D _cEnergyvsiphi_SubdetPM
Definition: RecHitTask.h:82
RecHitTask::_cOccupancyCutvsiphivsLS_SubdetPM
hcaldqm::Container2D _cOccupancyCutvsiphivsLS_SubdetPM
Definition: RecHitTask.h:124
RecHitTask::_gids
std::vector< HcalGenericDetId > _gids
Definition: RecHitTask.h:136
hcaldqm::ContainerXXX::begin
virtual CompactMap::const_iterator begin()
Definition: ContainerXXX.h:75
hcaldqm::UpdateFreq
UpdateFreq
Definition: DQTask.h:17
hcaldqm::flag::Flag::_state
State _state
Definition: Flag.h:79
edm::SortedCollection< HBHERecHit >::const_iterator
std::vector< HBHERecHit >::const_iterator const_iterator
Definition: SortedCollection.h:82
RecHitTask::_cTimingCut_SubdetPM
hcaldqm::Container1D _cTimingCut_SubdetPM
Definition: RecHitTask.h:91
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:66
RecHitTask::RecHitTask
RecHitTask(edm::ParameterSet const &)
Definition: RecHitTask.cc:8
hcaldqm::constants::SPIGOT_MIN
int const SPIGOT_MIN
Definition: Constants.h:136
RecHitTask::endLuminosityBlock
virtual void endLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &)
Definition: RecHitTask.cc:854
hcaldqm::ContainerXXX::reset
virtual void reset()
Definition: ContainerXXX.h:366
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
RecHitTask::_cTimingCutvsieta_Subdet
hcaldqm::ContainerProf1D _cTimingCutvsieta_Subdet
Definition: RecHitTask.h:100
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:20
hcaldqm::constants::FIBER_VME_MIN
int const FIBER_VME_MIN
Definition: Constants.h:141
RecHitTask::_cEnergyvsBX_SubdetPM
hcaldqm::ContainerProf1D _cEnergyvsBX_SubdetPM
Definition: RecHitTask.h:85
bookHistograms
void bookHistograms(fwlite::EventContainer &eventCont)
Definition: secvtxMassTemplates.cc:293
RecHitTask::_cOccupancyCut_ElectronicsVME
hcaldqm::Container2D _cOccupancyCut_ElectronicsVME
Definition: RecHitTask.h:117
hcaldqm::constants::FIBERCH_MIN
int const FIBERCH_MIN
Definition: Constants.h:151
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
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
RecHitTask::_cTimingCut_ElectronicsuTCA
hcaldqm::ContainerProf2D _cTimingCut_ElectronicsuTCA
Definition: RecHitTask.h:96
hcaldqm::ContainerXXX::end
virtual CompactMap::const_iterator end()
Definition: ContainerXXX.h:77
RecHitTask::_vflags
std::vector< hcaldqm::flag::Flag > _vflags
Definition: RecHitTask.h:60
RecHitTask::_cutE_HBHE
double _cutE_HBHE
Definition: RecHitTask.h:53
hcaldqm::DQModule::_ptype
ProcessingType _ptype
Definition: DQModule.h:59
hcaldqm::Container1D::fill
virtual void fill(uint32_t)
Definition: Container1D.cc:82
RecHitTask::_cEnergy_Subdet
hcaldqm::Container1D _cEnergy_Subdet
Definition: RecHitTask.h:80
hcaldqm::ContainerProf1D::initialize
virtual void initialize(std::string const &folder, hashfunctions::HashType, Quantity *, Quantity *, int debug=0)
Definition: ContainerProf1D.cc:26
RecHitTask::_cTimingCutvsLS_FED
hcaldqm::ContainerProf1D _cTimingCutvsLS_FED
Definition: RecHitTask.h:98
RecHitTask::_cTimingCut_FEDVME
hcaldqm::ContainerProf2D _cTimingCut_FEDVME
Definition: RecHitTask.h:93
RecHitTask::_cTimingCutvsLS_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsLS_SubdetPM
Definition: RecHitTask.h:99
hcaldqm::flag::Flag::reset
void reset()
Definition: Flag.h:76
RecHitTask::_cOccupancyCut_ElectronicsuTCA
hcaldqm::Container2D _cOccupancyCut_ElectronicsuTCA
Definition: RecHitTask.h:118
hcaldqm::ContainerSingle2D::setBinContent
virtual void setBinContent(int, int, int)
Definition: ContainerSingle2D.cc:209
RecHitTask::_cOccupancyCutvsBX_Subdet
hcaldqm::ContainerProf1D _cOccupancyCutvsBX_Subdet
Definition: RecHitTask.h:123
Logger::dqmthrow
void dqmthrow(std::string const &msg) const
Definition: Logger.h:15
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
RecHitTask::_tagPreHF
edm::InputTag _tagPreHF
Definition: RecHitTask.h:46
HcalDetId::ieta
int ieta() const
get the cell ieta
Definition: HcalDetId.h:56
RecHitTask::_cOccupancy_ElectronicsVME
hcaldqm::Container2D _cOccupancy_ElectronicsVME
Definition: RecHitTask.h:108
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
RecHitTask::_cutE_HF
double _cutE_HF
Definition: RecHitTask.h:53
hcaldqm::hashfunctions::hash_did
hash_function_did const hash_did[nHashType_did]
Definition: HashFunctions.h:167
HcalElectronicsMap::allPrecisionId
std::vector< HcalGenericDetId > allPrecisionId() const
Definition: HcalElectronicsMap.cc:139
hcaldqm::DQModule::_logger
Logger _logger
Definition: DQModule.h:70
RecHitTask::_cOccupancyCut_FEDuTCA
hcaldqm::Container2D _cOccupancyCut_FEDuTCA
Definition: RecHitTask.h:116
RecHitTask::_ehashmap
hcaldqm::electronicsmap::ElectronicsMap _ehashmap
Definition: RecHitTask.h:70
RecHitTask::meUnknownIds1LS
MonitorElement * meUnknownIds1LS
Definition: RecHitTask.h:133
RecHitTask::_tokPreHF
edm::EDGetTokenT< HFPreRecHitCollection > _tokPreHF
Definition: RecHitTask.h:51
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
RecHitTask::_cSummaryvsLS
hcaldqm::ContainerSingle2D _cSummaryvsLS
Definition: RecHitTask.h:138
RecHitTask::_thresh_unihf
double _thresh_unihf
Definition: RecHitTask.h:54
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
HcalElectronicsId::isVMEid
bool isVMEid() const
Definition: HcalElectronicsId.h:49
RecHitTask::_cOccupancy_ElectronicsuTCA
hcaldqm::Container2D _cOccupancy_ElectronicsuTCA
Definition: RecHitTask.h:109
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:277
hcaldqm::ContainerXXX::book
virtual void book(HcalElectronicsMap const *)
Definition: ContainerXXX.h:99
RecHitTask::_tagHO
edm::InputTag _tagHO
Definition: RecHitTask.h:44
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
hcaldqm::Container1D::book
virtual void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="")
Definition: Container1D.cc:957
RecHitTask::_filter_HF
hcaldqm::filter::HashFilter _filter_HF
Definition: RecHitTask.h:77
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
RecHitTask::_cOccupancy_FEDuTCA
hcaldqm::Container2D _cOccupancy_FEDuTCA
Definition: RecHitTask.h:107
RecHitTask::_tagHF
edm::InputTag _tagHF
Definition: RecHitTask.h:45
RecHitTask::_process
virtual void _process(edm::Event const &, edm::EventSetup const &)
Definition: RecHitTask.cc:415
RecHitTask::_cDAAsymmetryVsCharge_SubdetPM
hcaldqm::Container2D _cDAAsymmetryVsCharge_SubdetPM
Definition: RecHitTask.h:128
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
RecHitTask::_cTimingCut_HBHEPartition
hcaldqm::Container1D _cTimingCut_HBHEPartition
Definition: RecHitTask.h:92
RecHitTask::_filter_uTCA
hcaldqm::filter::HashFilter _filter_uTCA
Definition: RecHitTask.h:74
RecHitTask::_cOccupancyCutvsLS_Subdet
hcaldqm::ContainerProf1D _cOccupancyCutvsLS_Subdet
Definition: RecHitTask.h:119
hcaldqm::ContainerXXX::get
virtual STDTYPE & get(HcalDetId const &)
Definition: ContainerXXX.h:249
RecHitTask::_cOccupancyvsieta_Subdet
hcaldqm::Container1D _cOccupancyvsieta_Subdet
Definition: RecHitTask.h:112
hcaldqm::utilities::getFEDList
std::vector< int > getFEDList(HcalElectronicsMap const *)
Definition: Utilities.cc:95
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
RecHitTask::_hfPreRecHitsAvailable
bool _hfPreRecHitsAvailable
Definition: RecHitTask.h:47
hcaldqm::DQModule::_subsystem
std::string _subsystem
Definition: DQModule.h:64
RecHitTask::bookHistograms
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
Definition: RecHitTask.cc:38
RecHitTask::_cTimingCutvsBX_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsBX_SubdetPM
Definition: RecHitTask.h:102
RecHitTask::_tokHO
edm::EDGetTokenT< HORecHitCollection > _tokHO
Definition: RecHitTask.h:49
hcaldqm::DQTask::_vcdaqEids
std::vector< uint32_t > _vcdaqEids
Definition: DQTask.h:65
RecHitTask::_tagHBHE
edm::InputTag _tagHBHE
Definition: RecHitTask.h:43
RecHitTask::_cTimingCut_FEDuTCA
hcaldqm::ContainerProf2D _cTimingCut_FEDuTCA
Definition: RecHitTask.h:94
hcaldqm::ContainerXXX::initialize
virtual void initialize(hashfunctions::HashType, int debug=0)
Definition: ContainerXXX.h:91
RecHitTask::beginLuminosityBlock
virtual void beginLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &)
Definition: RecHitTask.cc:848
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
RecHitTask::_resetMonitors
virtual void _resetMonitors(hcaldqm::UpdateFreq)
Definition: RecHitTask.cc:401
hcaldqm::ContainerProf2D::book
virtual void book(DQMStore::IBooker &, HcalElectronicsMap const *, std::string subsystem="Hcal", std::string aux="")
Definition: ContainerProf2D.cc:51
RecHitTask::_cOccupancy_depth
hcaldqm::Container2D _cOccupancy_depth
Definition: RecHitTask.h:105
RecHitTask::_cDAAsymmetryMean_cut_depth
hcaldqm::ContainerProf2D _cDAAsymmetryMean_cut_depth
Definition: RecHitTask.h:129
RecHitTask::_cTimingCutvsiphi_SubdetPM
hcaldqm::ContainerProf1D _cTimingCutvsiphi_SubdetPM
Definition: RecHitTask.h:101
RecHitTask::_cTimingCut_depth
hcaldqm::ContainerProf2D _cTimingCut_depth
Definition: RecHitTask.h:97
HcalElectronicsId
Readout chain identification for Hcal.
Definition: HcalElectronicsId.h:32
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:37
RecHitTask::_cTimingCut_ElectronicsVME
hcaldqm::ContainerProf2D _cTimingCut_ElectronicsVME
Definition: RecHitTask.h:95
RecHitTask::_filter_VME
hcaldqm::filter::HashFilter _filter_VME
Definition: RecHitTask.h:73
RecHitTask::_cOccupancyCutvsieta_Subdet
hcaldqm::Container1D _cOccupancyCutvsieta_Subdet
Definition: RecHitTask.h:122
hcaldqm::utilities::isFEDHF
bool isFEDHF(HcalElectronicsId const &)
Definition: Utilities.cc:187
RecHitTask::_cOccupancyvsiphi_SubdetPM
hcaldqm::Container1D _cOccupancyvsiphi_SubdetPM
Definition: RecHitTask.h:111
hcaldqm::constants::SLOT_uTCA_MIN
int const SLOT_uTCA_MIN
Definition: Constants.h:123
HcalElectronicsId::isUTCAid
bool isUTCAid() const
Definition: HcalElectronicsId.h:50
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:303
edm::Run
Definition: Run.h:42
RecHitTask::_unknownIdsPresent
bool _unknownIdsPresent
Definition: RecHitTask.h:134
cuy.ib
ib
Definition: cuy.py:660
RecHitTask::_cDAAsymmetry_cut_SubdetPM
hcaldqm::Container1D _cDAAsymmetry_cut_SubdetPM
Definition: RecHitTask.h:130
RecHitTask::_cEnergyvsLS_SubdetPM
hcaldqm::ContainerProf1D _cEnergyvsLS_SubdetPM
Definition: RecHitTask.h:84
Generated for CMSSW Reference Manual by doxygen 1.8.11