CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EBOccupancyTask.cc
Go to the documentation of this file.
1 /*
2  * \file EBOccupancyTask.cc
3  *
4  * \author G. Della Ricca
5  * \author G. Franzoni
6  *
7 */
8 
9 #include <iostream>
10 #include <fstream>
11 
12 #include "FWCore/ServiceRegistry/interface/Service.h"
13 #include "FWCore/MessageLogger/interface/MessageLogger.h"
14 #include "FWCore/Framework/interface/ESHandle.h"
15 
16 #include "DQMServices/Core/interface/MonitorElement.h"
17 
18 #include "DQMServices/Core/interface/DQMStore.h"
19 
20 #include "DataFormats/EcalRawData/interface/EcalRawDataCollections.h"
21 #include "DataFormats/EcalDetId/interface/EBDetId.h"
22 #include "DataFormats/EcalDigi/interface/EBDataFrame.h"
23 #include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
24 #include "DataFormats/EcalRecHit/interface/EcalUncalibratedRecHit.h"
25 #include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
26 #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
27 
28 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h"
29 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgoRcd.h"
30 
31 #include "DQM/EcalCommon/interface/Numbers.h"
32 
33 #include "DQM/EcalBarrelMonitorTasks/interface/EBOccupancyTask.h"
34 
35 EBOccupancyTask::EBOccupancyTask(const edm::ParameterSet& ps){
36 
37  init_ = false;
38 
39  dqmStore_ = edm::Service<DQMStore>().operator->();
40 
41  prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
42 
43  subfolder_ = ps.getUntrackedParameter<std::string>("subfolder", "");
44 
45  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
46 
47  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
48 
49  EcalRawDataCollection_ = consumes<EcalRawDataCollection>(ps.getParameter<edm::InputTag>("EcalRawDataCollection"));
50  EBDigiCollection_ = consumes<EBDigiCollection>(ps.getParameter<edm::InputTag>("EBDigiCollection"));
51  EcalPnDiodeDigiCollection_ = consumes<EcalPnDiodeDigiCollection>(ps.getParameter<edm::InputTag>("EcalPnDiodeDigiCollection"));
52  EcalRecHitCollection_ = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("EcalRecHitCollection"));
53  EcalTrigPrimDigiCollection_ = consumes<EcalTrigPrimDigiCollection>(ps.getParameter<edm::InputTag>("EcalTrigPrimDigiCollection"));
54 
55  for (int i = 0; i < 36; i++) {
56  meOccupancy_[i] = 0;
57  meOccupancyMem_[i] = 0;
58  meEBRecHitEnergy_[i] = 0;
59  meSpectrum_[i] = 0;
60  }
61 
62  meEBRecHitSpectrum_ = 0;
63 
64  meEBDigiOccupancy_ = 0;
65  meEBDigiOccupancyProjEta_ = 0;
66  meEBDigiOccupancyProjPhi_ = 0;
67 
68  meEBRecHitOccupancy_ = 0;
69  meEBRecHitOccupancyProjEta_ = 0;
70  meEBRecHitOccupancyProjPhi_ = 0;
71 
72  meEBRecHitOccupancyThr_ = 0;
73  meEBRecHitOccupancyProjEtaThr_ = 0;
74  meEBRecHitOccupancyProjPhiThr_ = 0;
75 
76  meEBTrigPrimDigiOccupancy_ = 0;
77  meEBTrigPrimDigiOccupancyProjEta_ = 0;
78  meEBTrigPrimDigiOccupancyProjPhi_ = 0;
79 
80  meEBTrigPrimDigiOccupancyThr_ = 0;
81  meEBTrigPrimDigiOccupancyProjEtaThr_ = 0;
82  meEBTrigPrimDigiOccupancyProjPhiThr_ = 0;
83 
84  meEBTestPulseDigiOccupancy_ = 0;
85  meEBLaserDigiOccupancy_ = 0;
86  meEBPedestalDigiOccupancy_ = 0;
87 
88  recHitEnergyMin_ = 0.300; // GeV
89  trigPrimEtMin_ = 4.; // 4 ADCs == 1 GeV
90 
91 }
92 
93 EBOccupancyTask::~EBOccupancyTask(){
94 
95 }
96 
97 void EBOccupancyTask::beginJob(void){
98 
99  ievt_ = 0;
100 
101  if ( dqmStore_ ) {
102  dqmStore_->setCurrentFolder(prefixME_ + "/EBOccupancyTask");
103  if(subfolder_.size())
104  dqmStore_->setCurrentFolder(prefixME_ + "/EBOccupancyTask/" + subfolder_);
105  dqmStore_->rmdir(prefixME_ + "/EBOccupancyTask");
106  }
107 
108 }
109 
110 void EBOccupancyTask::beginRun(const edm::Run& r, const edm::EventSetup& c) {
111 
112  Numbers::initGeometry(c, false);
113 
114  if ( ! mergeRuns_ ) this->reset();
115 
116 }
117 
118 void EBOccupancyTask::endRun(const edm::Run& r, const edm::EventSetup& c) {
119 
120 }
121 
122 void EBOccupancyTask::reset(void) {
123 
124  for (int i = 0; i < 36; i++) {
125  if ( meOccupancy_[i] ) meOccupancy_[i]->Reset();
126  if ( meOccupancyMem_[i] ) meOccupancyMem_[i]->Reset();
127  if ( meEBRecHitEnergy_[i] ) meEBRecHitEnergy_[i]->Reset();
128  if ( meSpectrum_[i] ) meSpectrum_[i]->Reset();
129  }
130 
131  if ( meEBRecHitSpectrum_ ) meEBRecHitSpectrum_->Reset();
132 
133  if ( meEBDigiOccupancy_ ) meEBDigiOccupancy_->Reset();
134  if ( meEBDigiOccupancyProjEta_ ) meEBDigiOccupancyProjEta_->Reset();
135  if ( meEBDigiOccupancyProjPhi_ ) meEBDigiOccupancyProjPhi_->Reset();
136 
137  if ( meEBRecHitOccupancy_ ) meEBRecHitOccupancy_->Reset();
138  if ( meEBRecHitOccupancyProjEta_ ) meEBRecHitOccupancyProjEta_->Reset();
139  if ( meEBRecHitOccupancyProjPhi_ ) meEBRecHitOccupancyProjPhi_->Reset();
140 
141  if ( meEBRecHitOccupancyThr_ ) meEBRecHitOccupancyThr_->Reset();
142  if ( meEBRecHitOccupancyProjEtaThr_ ) meEBRecHitOccupancyProjEtaThr_->Reset();
143  if ( meEBRecHitOccupancyProjPhiThr_ ) meEBRecHitOccupancyProjPhiThr_->Reset();
144 
145  if ( meEBTrigPrimDigiOccupancy_ ) meEBTrigPrimDigiOccupancy_->Reset();
146  if ( meEBTrigPrimDigiOccupancyProjEta_ ) meEBTrigPrimDigiOccupancyProjEta_->Reset();
147  if ( meEBTrigPrimDigiOccupancyProjPhi_ ) meEBTrigPrimDigiOccupancyProjPhi_->Reset();
148 
149  if ( meEBTrigPrimDigiOccupancyThr_ ) meEBTrigPrimDigiOccupancyThr_->Reset();
150  if ( meEBTrigPrimDigiOccupancyProjEtaThr_ ) meEBTrigPrimDigiOccupancyProjEtaThr_->Reset();
151  if ( meEBTrigPrimDigiOccupancyProjPhiThr_ ) meEBTrigPrimDigiOccupancyProjPhiThr_->Reset();
152 
153  if ( meEBTestPulseDigiOccupancy_ ) meEBTestPulseDigiOccupancy_->Reset();
154  if ( meEBLaserDigiOccupancy_ ) meEBLaserDigiOccupancy_->Reset();
155  if ( meEBPedestalDigiOccupancy_ ) meEBPedestalDigiOccupancy_->Reset();
156 
157 }
158 
159 void EBOccupancyTask::setup(void){
160 
161  init_ = true;
162 
163  std::string name;
164 
165  if ( dqmStore_ ) {
166  dqmStore_->setCurrentFolder(prefixME_ + "/EBOccupancyTask");
167  if(subfolder_.size())
168  dqmStore_->setCurrentFolder(prefixME_ + "/EBOccupancyTask/" + subfolder_);
169 
170  for (int i = 0; i < 36; i++) {
171  name = "EBOT digi occupancy " + Numbers::sEB(i+1);
172  meOccupancy_[i] = dqmStore_->book2D(name, name, 85, 0., 85., 20, 0., 20.);
173  meOccupancy_[i]->setAxisTitle("ieta", 1);
174  meOccupancy_[i]->setAxisTitle("iphi", 2);
175  dqmStore_->tag(meOccupancy_[i], i+1);
176 
177  name = "EBOT MEM digi occupancy " + Numbers::sEB(i+1);
178  meOccupancyMem_[i] = dqmStore_->book2D(name, name, 10, 0., 10., 5, 0., 5.);
179  meOccupancyMem_[i]->setAxisTitle("pseudo-strip", 1);
180  meOccupancyMem_[i]->setAxisTitle("channel", 2);
181  dqmStore_->tag(meOccupancyMem_[i], i+1);
182 
183  name = "EBOT rec hit energy " + Numbers::sEB(i+1);
184  meEBRecHitEnergy_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 100, 0., 100., "s");
185  meEBRecHitEnergy_[i]->setAxisTitle("ieta", 1);
186  meEBRecHitEnergy_[i]->setAxisTitle("iphi", 2);
187  meEBRecHitEnergy_[i]->setAxisTitle("energy (GeV)", 3);
188  dqmStore_->tag(meEBRecHitEnergy_[i], i+1);
189 
190  name = "EBOT energy spectrum " + Numbers::sEB(i+1);
191  meSpectrum_[i] = dqmStore_->book1D(name, name, 100, 0., 1.5);
192  meSpectrum_[i]->setAxisTitle("energy (GeV)", 1);
193  dqmStore_->tag(meSpectrum_[i], i+1);
194  }
195 
196  name = "EBOT rec hit spectrum";
197  meEBRecHitSpectrum_ = dqmStore_->book1D(name, name, 100, 0., 10.);
198  meEBRecHitSpectrum_->setAxisTitle("energy (GeV)", 1);
199 
200  name = "EBOT digi occupancy";
201  meEBDigiOccupancy_ = dqmStore_->book2D(name, name, 72, 0., 360., 34, -85., 85.);
202  meEBDigiOccupancy_->setAxisTitle("jphi", 1);
203  meEBDigiOccupancy_->setAxisTitle("jeta", 2);
204  name = "EBOT digi occupancy projection eta";
205  meEBDigiOccupancyProjEta_ = dqmStore_->book1DD(name, name, 34, -85., 85.);
206  meEBDigiOccupancyProjEta_->setAxisTitle("jeta", 1);
207  meEBDigiOccupancyProjEta_->setAxisTitle("number of digis", 2);
208  name = "EBOT digi occupancy projection phi";
209  meEBDigiOccupancyProjPhi_ = dqmStore_->book1DD(name, name, 72, 0., 360.);
210  meEBDigiOccupancyProjPhi_->setAxisTitle("jphi", 1);
211  meEBDigiOccupancyProjPhi_->setAxisTitle("number of digis", 2);
212 
213  name = "EBOT rec hit occupancy";
214  meEBRecHitOccupancy_ = dqmStore_->book2D(name, name, 72, 0., 360., 34, -85., 85.);
215  meEBRecHitOccupancy_->setAxisTitle("jphi", 1);
216  meEBRecHitOccupancy_->setAxisTitle("jeta", 2);
217  name = "EBOT rec hit occupancy projection eta";
218  meEBRecHitOccupancyProjEta_ = dqmStore_->book1DD(name, name, 34, -85., 85.);
219  meEBRecHitOccupancyProjEta_->setAxisTitle("jeta", 1);
220  meEBRecHitOccupancyProjEta_->setAxisTitle("number of hits", 2);
221  name = "EBOT rec hit occupancy projection phi";
222  meEBRecHitOccupancyProjPhi_ = dqmStore_->book1DD(name, name, 72, 0., 360.);
223  meEBRecHitOccupancyProjPhi_->setAxisTitle("jphi", 1);
224  meEBRecHitOccupancyProjPhi_->setAxisTitle("number of hits", 2);
225 
226  name = "EBOT rec hit thr occupancy";
227  meEBRecHitOccupancyThr_ = dqmStore_->book2D(name, name, 72, 0., 360., 34, -85., 85.);
228  meEBRecHitOccupancyThr_->setAxisTitle("jphi", 1);
229  meEBRecHitOccupancyThr_->setAxisTitle("jeta", 2);
230  name = "EBOT rec hit thr occupancy projection eta";
231  meEBRecHitOccupancyProjEtaThr_ = dqmStore_->book1DD(name, name, 34, -85., 85.);
232  meEBRecHitOccupancyProjEtaThr_->setAxisTitle("jeta", 1);
233  meEBRecHitOccupancyProjEtaThr_->setAxisTitle("number of hits", 2);
234  name = "EBOT rec hit thr occupancy projection phi";
235  meEBRecHitOccupancyProjPhiThr_ = dqmStore_->book1DD(name, name, 72, 0., 360.);
236  meEBRecHitOccupancyProjPhiThr_->setAxisTitle("jphi", 1);
237  meEBRecHitOccupancyProjPhiThr_->setAxisTitle("number of hits", 2);
238 
239  name = "EBOT TP digi occupancy";
240  meEBTrigPrimDigiOccupancy_ = dqmStore_->book2D(name, name, 72, 0., 72., 34, -17., 17.);
241  meEBTrigPrimDigiOccupancy_->setAxisTitle("jphi'", 1);
242  meEBTrigPrimDigiOccupancy_->setAxisTitle("jeta'", 2);
243  name = "EBOT TP digi occupancy projection eta";
244  meEBTrigPrimDigiOccupancyProjEta_ = dqmStore_->book1DD(name, name, 34, -17., 17.);
245  meEBTrigPrimDigiOccupancyProjEta_->setAxisTitle("jeta'", 1);
246  meEBTrigPrimDigiOccupancyProjEta_->setAxisTitle("number of TP digis", 2);
247  name = "EBOT TP digi occupancy projection phi";
248  meEBTrigPrimDigiOccupancyProjPhi_ = dqmStore_->book1DD(name, name, 72, 0., 72.);
249  meEBTrigPrimDigiOccupancyProjPhi_->setAxisTitle("jphi'", 1);
250  meEBTrigPrimDigiOccupancyProjPhi_->setAxisTitle("number of TP digis", 2);
251 
252  name = "EBOT TP digi thr occupancy";
253  meEBTrigPrimDigiOccupancyThr_ = dqmStore_->book2D(name, name, 72, 0., 72., 34, -17., 17.);
254  meEBTrigPrimDigiOccupancyThr_->setAxisTitle("jphi'", 1);
255  meEBTrigPrimDigiOccupancyThr_->setAxisTitle("jeta'", 2);
256  name = "EBOT TP digi thr occupancy projection eta";
257  meEBTrigPrimDigiOccupancyProjEtaThr_ = dqmStore_->book1DD(name, name, 34, -17., 17.);
258  meEBTrigPrimDigiOccupancyProjEtaThr_->setAxisTitle("jeta'", 1);
259  meEBTrigPrimDigiOccupancyProjEtaThr_->setAxisTitle("number of TP digis", 2);
260  name = "EBOT TP digi thr occupancy projection phi";
261  meEBTrigPrimDigiOccupancyProjPhiThr_ = dqmStore_->book1DD(name, name, 72, 0., 72.);
262  meEBTrigPrimDigiOccupancyProjPhiThr_->setAxisTitle("jphi'", 1);
263  meEBTrigPrimDigiOccupancyProjPhiThr_->setAxisTitle("number of TP digis", 2);
264 
265  name = "EBOT test pulse digi occupancy";
266  meEBTestPulseDigiOccupancy_ = dqmStore_->book2D(name, name, 72, 0., 360., 34, -85., 85.);
267  meEBTestPulseDigiOccupancy_->setAxisTitle("jphi'", 1);
268  meEBTestPulseDigiOccupancy_->setAxisTitle("jeta'", 2);
269 
270  name = "EBOT laser digi occupancy";
271  meEBLaserDigiOccupancy_ = dqmStore_->book2D(name, name, 72, 0., 360., 34, -85., 85.);
272  meEBLaserDigiOccupancy_->setAxisTitle("jphi'", 1);
273  meEBLaserDigiOccupancy_->setAxisTitle("jeta'", 2);
274 
275  name = "EBOT pedestal digi occupancy";
276  meEBPedestalDigiOccupancy_ = dqmStore_->book2D(name, name, 72, 0., 360., 34, -85., 85.);
277  meEBPedestalDigiOccupancy_->setAxisTitle("jphi'", 1);
278  meEBPedestalDigiOccupancy_->setAxisTitle("jeta'", 2);
279 
280  }
281 
282 }
283 
284 void EBOccupancyTask::cleanup(void){
285 
286  if ( ! init_ ) return;
287 
288  if ( dqmStore_ ) {
289  dqmStore_->setCurrentFolder(prefixME_ + "/EBOccupancyTask");
290  if(subfolder_.size())
291  dqmStore_->setCurrentFolder(prefixME_ + "/EBOccupancyTask/" + subfolder_);
292 
293  for (int i = 0; i < 36; i++) {
294  if ( meOccupancy_[i] ) dqmStore_->removeElement( meOccupancy_[i]->getName() );
295  meOccupancy_[i] = 0;
296  if ( meOccupancyMem_[i] ) dqmStore_->removeElement( meOccupancyMem_[i]->getName() );
297  meOccupancyMem_[i] = 0;
298  if ( meEBRecHitEnergy_[i] ) dqmStore_->removeElement( meEBRecHitEnergy_[i]->getName() );
299  meEBRecHitEnergy_[i] = 0;
300  if ( meSpectrum_[i] ) dqmStore_->removeElement( meSpectrum_[i]->getName() );
301  meSpectrum_[i] = 0;
302  }
303 
304  if ( meEBRecHitSpectrum_ ) dqmStore_->removeElement( meEBRecHitSpectrum_->getName() );
305  meEBRecHitSpectrum_ = 0;
306 
307  if ( meEBDigiOccupancy_ ) dqmStore_->removeElement( meEBDigiOccupancy_->getName() );
308  meEBDigiOccupancy_ = 0;
309  if ( meEBDigiOccupancyProjEta_ ) dqmStore_->removeElement( meEBDigiOccupancyProjEta_->getName() );
310  meEBDigiOccupancyProjEta_ = 0;
311  if ( meEBDigiOccupancyProjPhi_ ) dqmStore_->removeElement( meEBDigiOccupancyProjPhi_->getName() );
312  meEBDigiOccupancyProjPhi_ = 0;
313 
314  if ( meEBRecHitOccupancy_ ) dqmStore_->removeElement( meEBRecHitOccupancy_->getName() );
315  meEBRecHitOccupancy_ = 0;
316  if ( meEBRecHitOccupancyProjEta_ ) dqmStore_->removeElement( meEBRecHitOccupancyProjEta_->getName() );
317  meEBRecHitOccupancyProjEta_ = 0;
318  if ( meEBRecHitOccupancyProjPhi_ ) dqmStore_->removeElement( meEBRecHitOccupancyProjPhi_->getName() );
319  meEBRecHitOccupancyProjPhi_ = 0;
320 
321  if ( meEBRecHitOccupancyThr_ ) dqmStore_->removeElement( meEBRecHitOccupancyThr_->getName() );
322  meEBRecHitOccupancyThr_ = 0;
323  if ( meEBRecHitOccupancyProjEtaThr_ ) dqmStore_->removeElement( meEBRecHitOccupancyProjEtaThr_->getName() );
324  meEBRecHitOccupancyProjEtaThr_ = 0;
325  if ( meEBRecHitOccupancyProjPhiThr_ ) dqmStore_->removeElement( meEBRecHitOccupancyProjPhiThr_->getName() );
326  meEBRecHitOccupancyProjPhiThr_ = 0;
327 
328  if ( meEBTrigPrimDigiOccupancy_ ) dqmStore_->removeElement( meEBTrigPrimDigiOccupancy_->getName() );
329  meEBTrigPrimDigiOccupancy_ = 0;
330  if ( meEBTrigPrimDigiOccupancyProjEta_ ) dqmStore_->removeElement( meEBTrigPrimDigiOccupancyProjEta_->getName() );
331  meEBTrigPrimDigiOccupancyProjEta_ = 0;
332  if ( meEBTrigPrimDigiOccupancyProjPhi_ ) dqmStore_->removeElement( meEBTrigPrimDigiOccupancyProjPhi_->getName() );
333  meEBTrigPrimDigiOccupancyProjPhi_ = 0;
334 
335  if ( meEBTrigPrimDigiOccupancyThr_ ) dqmStore_->removeElement( meEBTrigPrimDigiOccupancyThr_->getName() );
336  meEBTrigPrimDigiOccupancyThr_ = 0;
337  if ( meEBTrigPrimDigiOccupancyProjEtaThr_ ) dqmStore_->removeElement( meEBTrigPrimDigiOccupancyProjEtaThr_->getName() );
338  meEBTrigPrimDigiOccupancyProjEtaThr_ = 0;
339  if ( meEBTrigPrimDigiOccupancyProjPhiThr_ ) dqmStore_->removeElement( meEBTrigPrimDigiOccupancyProjPhiThr_->getName() );
340  meEBTrigPrimDigiOccupancyProjPhiThr_ = 0;
341 
342  if ( meEBTestPulseDigiOccupancy_ ) dqmStore_->removeElement( meEBTestPulseDigiOccupancy_->getName() );
343  meEBTestPulseDigiOccupancy_ = 0;
344 
345  if ( meEBLaserDigiOccupancy_ ) dqmStore_->removeElement( meEBLaserDigiOccupancy_->getName() );
346  meEBLaserDigiOccupancy_ = 0;
347 
348  if ( meEBPedestalDigiOccupancy_ ) dqmStore_->removeElement( meEBPedestalDigiOccupancy_->getName() );
349  meEBPedestalDigiOccupancy_ = 0;
350 
351  }
352 
353  init_ = false;
354 
355 }
356 
357 void EBOccupancyTask::endJob(void) {
358 
359  edm::LogInfo("EBOccupancyTask") << "analyzed " << ievt_ << " events";
360 
361  if ( enableCleanup_ ) this->cleanup();
362 
363 }
364 
365 void EBOccupancyTask::analyze(const edm::Event& e, const edm::EventSetup& c){
366 
367  if ( ! init_ ) this->setup();
368 
369  ievt_++;
370 
371  int runType[36] = { notdata };
372 
373  edm::Handle<EcalRawDataCollection> dcchs;
374 
375  if ( e.getByToken(EcalRawDataCollection_, dcchs) ) {
376 
377  for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
378 
379  if ( Numbers::subDet( *dcchItr ) != EcalBarrel ) continue;
380 
381  int ism = Numbers::iSM( *dcchItr, EcalBarrel );
382 
383  int runtype = dcchItr->getRunType();
384 
385  if ( runtype == EcalDCCHeaderBlock::COSMIC ||
386  runtype == EcalDCCHeaderBlock::MTCC ||
387  runtype == EcalDCCHeaderBlock::COSMICS_GLOBAL ||
388  runtype == EcalDCCHeaderBlock::PHYSICS_GLOBAL ||
389  runtype == EcalDCCHeaderBlock::COSMICS_LOCAL ||
390  runtype == EcalDCCHeaderBlock::PHYSICS_LOCAL ) runType[ism-1] = physics;
391  if ( runtype == EcalDCCHeaderBlock::TESTPULSE_MGPA ||
392  runtype == EcalDCCHeaderBlock::TESTPULSE_GAP ) runType[ism-1] = testpulse;
393  if ( runtype == EcalDCCHeaderBlock::LASER_STD ||
394  runtype == EcalDCCHeaderBlock::LASER_GAP ) runType[ism-1] = laser;
395  if ( runtype == EcalDCCHeaderBlock::PEDESTAL_STD ||
396  runtype == EcalDCCHeaderBlock::PEDESTAL_GAP ) runType[ism-1] = pedestal;
397 
398  }
399 
400  } else {
401  edm::LogWarning("EBOccupancyTask") << "EcalRawDataCollection not available";
402  }
403 
404  edm::Handle<EBDigiCollection> digis;
405 
406  if ( e.getByToken(EBDigiCollection_, digis) ) {
407 
408  int nebd = digis->size();
409  LogDebug("EBOccupancyTask") << "event " << ievt_ << " digi collection size " << nebd;
410 
411  for ( EBDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr ) {
412 
413  EBDetId id = digiItr->id();
414 
415  int ic = id.ic();
416  int ie = (ic-1)/20 + 1;
417  int ip = (ic-1)%20 + 1;
418 
419  int ism = Numbers::iSM( id );
420 
421  float xie = ie - 0.5;
422  float xip = ip - 0.5;
423 
424  if ( xie <= 0. || xie >= 85. || xip <= 0. || xip >= 20. ) {
425  edm::LogWarning("EBOccupancyTask") << " det id = " << id;
426  edm::LogWarning("EBOccupancyTask") << " sm, ieta, iphi " << ism << " " << ie << " " << ip;
427  edm::LogWarning("EBOccupancyTask") << " xie, xip " << xie << " " << xip;
428  }
429 
430  if ( meOccupancy_[ism-1] ) meOccupancy_[ism-1]->Fill(xie, xip);
431 
432  int ebeta = id.ieta();
433  int ebphi = id.iphi();
434 
435  float xebeta = ebeta - 0.5*id.zside();
436  float xebphi = ebphi - 0.5;
437 
438  if ( runType[ism-1] == physics || runType[ism-1] == notdata ) {
439 
440  if ( meEBDigiOccupancy_ ) meEBDigiOccupancy_->Fill( xebphi, xebeta );
441  if ( meEBDigiOccupancyProjEta_ ) meEBDigiOccupancyProjEta_->Fill( xebeta );
442  if ( meEBDigiOccupancyProjPhi_ ) meEBDigiOccupancyProjPhi_->Fill( xebphi );
443 
444  }
445 
446  if ( runType[ism-1] == testpulse ) {
447 
448  if ( meEBTestPulseDigiOccupancy_ ) meEBTestPulseDigiOccupancy_->Fill( xebphi, xebeta );
449 
450  }
451 
452  if ( runType[ism-1] == laser ) {
453 
454  if ( meEBLaserDigiOccupancy_ ) meEBLaserDigiOccupancy_->Fill( xebphi, xebeta );
455 
456  }
457 
458  if ( runType[ism-1] == pedestal ) {
459 
460  if ( meEBPedestalDigiOccupancy_ ) meEBPedestalDigiOccupancy_->Fill( xebphi, xebeta );
461 
462  }
463 
464  }
465 
466  } else {
467 
468  edm::LogWarning("EBOccupancyTask") << "EBDigiCollection not available";
469 
470  }
471 
472  edm::Handle<EcalPnDiodeDigiCollection> PNs;
473 
474  if ( e.getByToken(EcalPnDiodeDigiCollection_, PNs) ) {
475 
476  // filling mem occupancy only for the 5 channels belonging
477  // to a fully reconstructed PN's
478 
479  for ( EcalPnDiodeDigiCollection::const_iterator pnItr = PNs->begin(); pnItr != PNs->end(); ++pnItr ) {
480 
481  if ( Numbers::subDet( pnItr->id() ) != EcalBarrel ) continue;
482 
483  int ism = Numbers::iSM( pnItr->id() );
484 
485  float PnId = pnItr->id().iPnId();
486 
487  PnId = PnId - 0.5;
488  float st = 0.0;
489 
490  for (int chInStrip = 1; chInStrip <= 5; chInStrip++){
491  if ( meOccupancyMem_[ism-1] ) {
492  st = chInStrip - 0.5;
493  meOccupancyMem_[ism-1]->Fill(PnId, st);
494  }
495  }
496 
497  }
498 
499  } else {
500 
501  edm::LogWarning("EBOccupancyTask") << "EcalPnDiodeDigiCollection not available";
502 
503  }
504 
505  edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
506  c.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
507 
508  edm::Handle<EcalRecHitCollection> rechits;
509 
510  if ( e.getByToken(EcalRecHitCollection_, rechits) ) {
511 
512  int nebrh = rechits->size();
513  LogDebug("EBOccupancyTask") << "event " << ievt_ << " rec hits collection size " << nebrh;
514 
515  for ( EcalRecHitCollection::const_iterator rechitItr = rechits->begin(); rechitItr != rechits->end(); ++rechitItr ) {
516 
517  EBDetId id = rechitItr->id();
518 
519  // global coordinates
520  int ebeta = id.ieta();
521  int ebphi = id.iphi();
522 
523  float xebeta = ebeta - 0.5*id.zside();
524  float xebphi = ebphi - 0.5;
525 
526  int ism = Numbers::iSM( id );
527 
528  // local coordinates
529  int ic = id.ic();
530  int ie = (ic-1)/20 + 1;
531  int ip = (ic-1)%20 + 1;
532 
533  float xie = ie - 0.5;
534  float xip = ip - 0.5;
535 
536  if ( runType[ism-1] == physics || runType[ism-1] == notdata ) {
537 
538  if ( meEBRecHitOccupancy_ ) meEBRecHitOccupancy_->Fill( xebphi, xebeta );
539  if ( meEBRecHitOccupancyProjEta_ ) meEBRecHitOccupancyProjEta_->Fill( xebeta );
540  if ( meEBRecHitOccupancyProjPhi_ ) meEBRecHitOccupancyProjPhi_->Fill( xebphi );
541 
542  uint32_t flag = rechitItr->recoFlag();
543 
544  uint32_t sev = sevlv->severityLevel( id, *rechits);
545 
546  if ( rechitItr->energy() > recHitEnergyMin_ && flag == EcalRecHit::kGood && sev == EcalSeverityLevel::kGood ) {
547 
548  if ( meEBRecHitOccupancyThr_ ) meEBRecHitOccupancyThr_->Fill( xebphi, xebeta );
549  if ( meEBRecHitOccupancyProjEtaThr_ ) meEBRecHitOccupancyProjEtaThr_->Fill( xebeta );
550  if ( meEBRecHitOccupancyProjPhiThr_ ) meEBRecHitOccupancyProjPhiThr_->Fill( xebphi );
551 
552  }
553 
554  if ( flag == EcalRecHit::kGood && sev == EcalSeverityLevel::kGood ) {
555  if ( meEBRecHitEnergy_[ism-1] ) meEBRecHitEnergy_[ism-1]->Fill( xie, xip, rechitItr->energy() );
556  if ( meSpectrum_[ism-1] ) meSpectrum_[ism-1]->Fill( rechitItr->energy() );
557  if ( meEBRecHitSpectrum_ ) meEBRecHitSpectrum_->Fill( rechitItr->energy() );
558  }
559 
560  }
561  }
562 
563  } else {
564 
565  edm::LogWarning("EBOccupancyTask") << "EcalRecHitCollection not available";
566 
567  }
568 
569  edm::Handle<EcalTrigPrimDigiCollection> trigPrimDigis;
570 
571  if ( e.getByToken(EcalTrigPrimDigiCollection_, trigPrimDigis) ) {
572 
573  int nebtpg = trigPrimDigis->size();
574  LogDebug("EBOccupancyTask") << "event " << ievt_ << " trigger primitives digis collection size " << nebtpg;
575 
576  for ( EcalTrigPrimDigiCollection::const_iterator tpdigiItr = trigPrimDigis->begin(); tpdigiItr != trigPrimDigis->end(); ++tpdigiItr ) {
577 
578  if ( Numbers::subDet( tpdigiItr->id() ) != EcalBarrel ) continue;
579 
580  int ebeta = tpdigiItr->id().ieta();
581  int ebphi = tpdigiItr->id().iphi();
582 
583  // phi_tower: change the range from global to SM-local
584  // phi==0 is in the middle of a SM
585  ebphi = ebphi + 2;
586  if ( ebphi > 72 ) ebphi = ebphi - 72;
587 
588  float xebeta = ebeta-0.5*tpdigiItr->id().zside();
589  float xebphi = ebphi-0.5;
590 
591  int ism = Numbers::iSM( tpdigiItr->id() );
592 
593  if ( runType[ism-1] == physics || runType[ism-1] == notdata ) {
594 
595  if ( meEBTrigPrimDigiOccupancy_ ) meEBTrigPrimDigiOccupancy_->Fill( xebphi, xebeta );
596  if ( meEBTrigPrimDigiOccupancyProjEta_ ) meEBTrigPrimDigiOccupancyProjEta_->Fill( xebeta );
597  if ( meEBTrigPrimDigiOccupancyProjPhi_ ) meEBTrigPrimDigiOccupancyProjPhi_->Fill( xebphi );
598 
599  if ( tpdigiItr->compressedEt() > trigPrimEtMin_ ) {
600 
601  if ( meEBTrigPrimDigiOccupancyThr_ ) meEBTrigPrimDigiOccupancyThr_->Fill( xebphi, xebeta );
602  if ( meEBTrigPrimDigiOccupancyProjEtaThr_ ) meEBTrigPrimDigiOccupancyProjEtaThr_->Fill( xebeta );
603  if ( meEBTrigPrimDigiOccupancyProjPhiThr_ ) meEBTrigPrimDigiOccupancyProjPhiThr_->Fill( xebphi );
604 
605  }
606  }
607  }
608 
609  } else {
610 
611  edm::LogWarning("EBOccupancyTask") << "EcalTrigPrimDigiCollection not available";
612 
613  }
614 
615 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
EBOccupancyTask::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Analyze.
Definition: EBOccupancyTask.cc:365
EBOccupancyTask::~EBOccupancyTask
virtual ~EBOccupancyTask()
Destructor.
Definition: EBOccupancyTask.cc:93
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
MonitorElement::getName
const std::string & getName(void) const
get name of ME
Definition: MonitorElement.h:108
i
int i
Definition: DBlmapReader.cc:9
EBOccupancyTask::meEBTestPulseDigiOccupancy_
MonitorElement * meEBTestPulseDigiOccupancy_
Definition: EBOccupancyTask.h:102
EBOccupancyTask::meEBTrigPrimDigiOccupancyProjEtaThr_
MonitorElement * meEBTrigPrimDigiOccupancyProjEtaThr_
Definition: EBOccupancyTask.h:100
EBOccupancyTask::meEBTrigPrimDigiOccupancy_
MonitorElement * meEBTrigPrimDigiOccupancy_
Definition: EBOccupancyTask.h:96
EBOccupancyTask::meEBRecHitOccupancyProjEta_
MonitorElement * meEBRecHitOccupancyProjEta_
Definition: EBOccupancyTask.h:91
EBOccupancyTask::EcalRecHitCollection_
edm::EDGetTokenT< EcalRecHitCollection > EcalRecHitCollection_
Definition: EBOccupancyTask.h:75
edm::DataFrameContainer::const_iterator
boost::transform_iterator< IterHelp, boost::counting_iterator< int > > const_iterator
Definition: DataFrameContainer.h:51
EBOccupancyTask::meOccupancy_
MonitorElement * meOccupancy_[36]
Definition: EBOccupancyTask.h:81
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:872
EBOccupancyTask::prefixME_
std::string prefixME_
Definition: EBOccupancyTask.h:64
DQMStore::rmdir
void rmdir(const std::string &fullpath)
Definition: DQMStore.cc:2730
EBOccupancyTask::meEBTrigPrimDigiOccupancyProjPhi_
MonitorElement * meEBTrigPrimDigiOccupancyProjPhi_
Definition: EBOccupancyTask.h:98
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434
Numbers.h
Some &quot;id&quot; conversions.
EBOccupancyTask::beginJob
void beginJob(void)
BeginJob.
Definition: EBOccupancyTask.cc:97
EBOccupancyTask::meEBDigiOccupancy_
MonitorElement * meEBDigiOccupancy_
Definition: EBOccupancyTask.h:87
Numbers::sEB
static std::string sEB(const unsigned ism)
Definition: Numbers.cc:91
EBOccupancyTask::EcalRawDataCollection_
edm::EDGetTokenT< EcalRawDataCollection > EcalRawDataCollection_
Definition: EBOccupancyTask.h:72
edm::SortedCollection< EcalDCCHeaderBlock >::const_iterator
std::vector< EcalDCCHeaderBlock >::const_iterator const_iterator
Definition: SortedCollection.h:80
EBOccupancyTask::dqmStore_
DQMStore * dqmStore_
Definition: EBOccupancyTask.h:62
EBOccupancyTask::meEBLaserDigiOccupancy_
MonitorElement * meEBLaserDigiOccupancy_
Definition: EBOccupancyTask.h:103
DQMStore::book1DD
MonitorElement * book1DD(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1S histogram.
Definition: DQMStore.cc:904
EBOccupancyTask::meEBRecHitOccupancyThr_
MonitorElement * meEBRecHitOccupancyThr_
Definition: EBOccupancyTask.h:93
EBOccupancyTask::meEBRecHitEnergy_
MonitorElement * meEBRecHitEnergy_[36]
Definition: EBOccupancyTask.h:83
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
EBOccupancyTask::meEBRecHitSpectrum_
MonitorElement * meEBRecHitSpectrum_
Definition: EBOccupancyTask.h:86
EBOccupancyTask::meOccupancyMem_
MonitorElement * meOccupancyMem_[36]
Definition: EBOccupancyTask.h:82
EBOccupancyTask::meEBTrigPrimDigiOccupancyProjEta_
MonitorElement * meEBTrigPrimDigiOccupancyProjEta_
Definition: EBOccupancyTask.h:97
EBOccupancyTask::meEBTrigPrimDigiOccupancyThr_
MonitorElement * meEBTrigPrimDigiOccupancyThr_
Definition: EBOccupancyTask.h:99
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:161
DQMStore::tag
void tag(MonitorElement *me, unsigned int myTag)
Definition: DQMStore.cc:1509
EBOccupancyTask::reset
void reset(void)
Reset.
Definition: EBOccupancyTask.cc:122
EBOccupancyTask::cleanup
void cleanup(void)
Cleanup.
Definition: EBOccupancyTask.cc:284
EBOccupancyTask::EBDigiCollection_
edm::EDGetTokenT< EBDigiCollection > EBDigiCollection_
Definition: EBOccupancyTask.h:73
EBOccupancyTask::meEBRecHitOccupancyProjPhi_
MonitorElement * meEBRecHitOccupancyProjPhi_
Definition: EBOccupancyTask.h:92
EBOccupancyTask::subfolder_
std::string subfolder_
Definition: EBOccupancyTask.h:66
DQMStore::removeElement
void removeElement(const std::string &name)
Definition: DQMStore.cc:2772
EBOccupancyTask::EcalTrigPrimDigiCollection_
edm::EDGetTokenT< EcalTrigPrimDigiCollection > EcalTrigPrimDigiCollection_
Definition: EBOccupancyTask.h:76
EBOccupancyTask::meSpectrum_
MonitorElement * meSpectrum_[36]
Definition: EBOccupancyTask.h:84
EBDetId::ieta
int ieta() const
get the crystal ieta
Definition: EBDetId.h:51
EBOccupancyTask::meEBRecHitOccupancyProjPhiThr_
MonitorElement * meEBRecHitOccupancyProjPhiThr_
Definition: EBOccupancyTask.h:95
EBOccupancyTask::endRun
void endRun(const edm::Run &r, const edm::EventSetup &c)
EndRun.
Definition: EBOccupancyTask.cc:118
EBDetId::ic
int ic() const
get ECAL/crystal number inside SM
Definition: EBDetId.cc:46
environment_file_cff.runType
tuple runType
Definition: environment_file_cff.py:30
EBOccupancyTask::meEBDigiOccupancyProjEta_
MonitorElement * meEBDigiOccupancyProjEta_
Definition: EBOccupancyTask.h:88
Numbers::initGeometry
static void initGeometry(const edm::EventSetup &setup, bool verbose=false)
Definition: Numbers.cc:47
EBOccupancyTask::meEBPedestalDigiOccupancy_
MonitorElement * meEBPedestalDigiOccupancy_
Definition: EBOccupancyTask.h:104
EBOccupancyTask::EcalPnDiodeDigiCollection_
edm::EDGetTokenT< EcalPnDiodeDigiCollection > EcalPnDiodeDigiCollection_
Definition: EBOccupancyTask.h:74
EBOccupancyTask::meEBDigiOccupancyProjPhi_
MonitorElement * meEBDigiOccupancyProjPhi_
Definition: EBOccupancyTask.h:89
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
EBOccupancyTask::meEBTrigPrimDigiOccupancyProjPhiThr_
MonitorElement * meEBTrigPrimDigiOccupancyProjPhiThr_
Definition: EBOccupancyTask.h:101
HI_PhotonSkim_cff.rechits
tuple rechits
Definition: HI_PhotonSkim_cff.py:76
Numbers::iSM
static unsigned iSM(const unsigned ism, const EcalSubdetector subdet)
Definition: Numbers.cc:243
EBOccupancyTask::meEBRecHitOccupancyProjEtaThr_
MonitorElement * meEBRecHitOccupancyProjEtaThr_
Definition: EBOccupancyTask.h:94
Numbers::subDet
static EcalSubdetector subDet(const EBDetId &id)
Definition: Numbers.cc:142
EBOccupancyTask::endJob
void endJob(void)
EndJob.
Definition: EBOccupancyTask.cc:357
EBOccupancyTask::EBOccupancyTask
EBOccupancyTask(const edm::ParameterSet &ps)
Constructor.
Definition: EBOccupancyTask.cc:35
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:1000
ecalpyutils::ism
int ism(int ieta, int iphi)
Definition: EcalPyUtils.cc:56
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:928
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:434
EBOccupancyTask::meEBRecHitOccupancy_
MonitorElement * meEBRecHitOccupancy_
Definition: EBOccupancyTask.h:90
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:584
edm::Run
Definition: Run.h:41
EBOccupancyTask::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
BeginRun.
Definition: EBOccupancyTask.cc:110
EBOccupancyTask::setup
void setup(void)
Setup.
Definition: EBOccupancyTask.cc:159
DQMStore::bookProfile2D
MonitorElement * bookProfile2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, int nchZ, double lowZ, double highZ, const char *option="s")
Definition: DQMStore.cc:1330
Generated for CMSSW Reference Manual by doxygen 1.8.5