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