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