CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EEOccupancyTask.cc
Go to the documentation of this file.
1 /*
2  * \file EEOccupancyTask.cc
3  *
4  * $Date: 2012/04/27 13:46:15 $
5  * $Revision: 1.94 $
6  * \author G. Della Ricca
7  * \author G. Franzoni
8  *
9 */
10 
11 #include <iostream>
12 #include <fstream>
13 #include <vector>
14 
15 #include "FWCore/ServiceRegistry/interface/Service.h"
16 #include "FWCore/MessageLogger/interface/MessageLogger.h"
17 #include "FWCore/Framework/interface/ESHandle.h"
18 
19 #include "DQMServices/Core/interface/MonitorElement.h"
20 
21 #include "DQMServices/Core/interface/DQMStore.h"
22 
23 #include "DataFormats/EcalRawData/interface/EcalRawDataCollections.h"
24 #include "DataFormats/EcalDetId/interface/EEDetId.h"
25 #include "DataFormats/EcalDigi/interface/EEDataFrame.h"
26 #include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
27 #include "DataFormats/EcalRecHit/interface/EcalUncalibratedRecHit.h"
28 #include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
29 #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
30 
31 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h"
32 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgoRcd.h"
33 
34 #include "DQM/EcalCommon/interface/Numbers.h"
35 
36 #include "DQM/EcalEndcapMonitorTasks/interface/EEOccupancyTask.h"
37 
38 EEOccupancyTask::EEOccupancyTask(const edm::ParameterSet& ps){
39 
40  init_ = false;
41 
42  initCaloGeometry_ = false;
43 
44  dqmStore_ = edm::Service<DQMStore>().operator->();
45 
46  prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
47 
48  subfolder_ = ps.getUntrackedParameter<std::string>("subfolder", "");
49 
50  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
51 
52  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
53 
54  EcalRawDataCollection_ = ps.getParameter<edm::InputTag>("EcalRawDataCollection");
55  EEDigiCollection_ = ps.getParameter<edm::InputTag>("EEDigiCollection");
56  EcalPnDiodeDigiCollection_ = ps.getParameter<edm::InputTag>("EcalPnDiodeDigiCollection");
57  EcalRecHitCollection_ = ps.getParameter<edm::InputTag>("EcalRecHitCollection");
58  EcalTrigPrimDigiCollection_ = ps.getParameter<edm::InputTag>("EcalTrigPrimDigiCollection");
59 
60  for (int i = 0; i < 18; i++) {
61  meOccupancy_[i] = 0;
62  meOccupancyMem_[i] = 0;
63  meEERecHitEnergy_[i] = 0;
64  meSpectrum_[i] = 0;
65  }
66 
67  meEERecHitSpectrum_[0] = 0;
68  meEERecHitSpectrum_[1] = 0;
69 
70  meEEDigiOccupancy_[0] = 0;
71  meEEDigiOccupancyProEta_[0] = 0;
72  meEEDigiOccupancyProPhi_[0] = 0;
73  meEEDigiOccupancy_[1] = 0;
74  meEEDigiOccupancyProEta_[1] = 0;
75  meEEDigiOccupancyProPhi_[1] = 0;
76 
77  meEERecHitOccupancy_[0] = 0;
78  meEERecHitOccupancyProEta_[0] = 0;
79  meEERecHitOccupancyProPhi_[0] = 0;
80  meEERecHitOccupancy_[1] = 0;
81  meEERecHitOccupancyProEta_[1] = 0;
82  meEERecHitOccupancyProPhi_[1] = 0;
83 
84  meEERecHitOccupancyThr_[0] = 0;
85  meEERecHitOccupancyProEtaThr_[0] = 0;
86  meEERecHitOccupancyProPhiThr_[0] = 0;
87  meEERecHitOccupancyThr_[1] = 0;
88  meEERecHitOccupancyProEtaThr_[1] = 0;
89  meEERecHitOccupancyProPhiThr_[1] = 0;
90 
91  meEETrigPrimDigiOccupancy_[0] = 0;
92  meEETrigPrimDigiOccupancyProEta_[0] = 0;
93  meEETrigPrimDigiOccupancyProPhi_[0] = 0;
94  meEETrigPrimDigiOccupancy_[1] = 0;
95  meEETrigPrimDigiOccupancyProEta_[1] = 0;
96  meEETrigPrimDigiOccupancyProPhi_[1] = 0;
97 
98  meEETrigPrimDigiOccupancyThr_[0] = 0;
99  meEETrigPrimDigiOccupancyProEtaThr_[0] = 0;
100  meEETrigPrimDigiOccupancyProPhiThr_[0] = 0;
101  meEETrigPrimDigiOccupancyThr_[1] = 0;
102  meEETrigPrimDigiOccupancyProEtaThr_[1] = 0;
103  meEETrigPrimDigiOccupancyProPhiThr_[1] = 0;
104 
105  meEETestPulseDigiOccupancy_[0] = 0;
106  meEETestPulseDigiOccupancy_[1] = 0;
107 
108  meEELaserDigiOccupancy_[0] = 0;
109  meEELaserDigiOccupancy_[1] = 0;
110 
111  meEELedDigiOccupancy_[0] = 0;
112  meEELedDigiOccupancy_[1] = 0;
113 
114  meEEPedestalDigiOccupancy_[0] = 0;
115  meEEPedestalDigiOccupancy_[1] = 0;
116 
117  recHitEnergyMin_ = 0.500; // GeV
118  trigPrimEtMin_ = 4.; // 4 ADCs == 1 GeV
119 
120  for (int i = 0; i < EEDetId::kSizeForDenseIndexing; i++) {
121  geometryEE[i][0] = 0;
122  geometryEE[i][1] = 0;
123  }
124 
125 }
126 
127 EEOccupancyTask::~EEOccupancyTask(){
128 
129 }
130 
131 void EEOccupancyTask::beginJob(void){
132 
133  ievt_ = 0;
134 
135  if ( dqmStore_ ) {
136  dqmStore_->setCurrentFolder(prefixME_ + "/EEOccupancyTask");
137  if(subfolder_.size())
138  dqmStore_->setCurrentFolder(prefixME_ + "/EEOccupancyTask/" + subfolder_);
139  dqmStore_->rmdir(prefixME_ + "/EEOccupancyTask");
140  }
141 
142 }
143 
144 void EEOccupancyTask::beginRun(const edm::Run& r, const edm::EventSetup& c) {
145 
146  Numbers::initGeometry(c, false);
147 
148  if( !initCaloGeometry_ ) {
149  c.get<CaloGeometryRecord>().get(pGeometry_);
150  initCaloGeometry_ = true;
151  }
152 
153  if ( ! mergeRuns_ ) this->reset();
154 
155 }
156 
157 void EEOccupancyTask::endRun(const edm::Run& r, const edm::EventSetup& c) {
158 
159 }
160 
161 void EEOccupancyTask::reset(void) {
162 
163  for (int i = 0; i < 18; i++) {
164  if ( meOccupancy_[i] ) meOccupancy_[i]->Reset();
165  if ( meOccupancyMem_[i] ) meOccupancyMem_[i]->Reset();
166  if ( meEERecHitEnergy_[i] ) meEERecHitEnergy_[i]->Reset();
167  if ( meSpectrum_[i] ) meSpectrum_[i]->Reset();
168  }
169 
170  if ( meEERecHitSpectrum_[0] ) meEERecHitSpectrum_[0]->Reset();
171  if ( meEERecHitSpectrum_[1] ) meEERecHitSpectrum_[1]->Reset();
172 
173  if ( meEEDigiOccupancy_[0] ) meEEDigiOccupancy_[0]->Reset();
174  if ( meEEDigiOccupancyProEta_[0] ) meEEDigiOccupancyProEta_[0]->Reset();
175  if ( meEEDigiOccupancyProPhi_[0] ) meEEDigiOccupancyProPhi_[0]->Reset();
176  if ( meEEDigiOccupancy_[1] ) meEEDigiOccupancy_[1]->Reset();
177  if ( meEEDigiOccupancyProEta_[1] ) meEEDigiOccupancyProEta_[1]->Reset();
178  if ( meEEDigiOccupancyProPhi_[1] ) meEEDigiOccupancyProPhi_[1]->Reset();
179 
180  if ( meEERecHitOccupancy_[0] ) meEERecHitOccupancy_[0]->Reset();
181  if ( meEERecHitOccupancyProEta_[0] ) meEERecHitOccupancyProEta_[0]->Reset();
182  if ( meEERecHitOccupancyProPhi_[0] ) meEERecHitOccupancyProPhi_[0]->Reset();
183  if ( meEERecHitOccupancy_[1] ) meEERecHitOccupancy_[1]->Reset();
184  if ( meEERecHitOccupancyProEta_[1] ) meEERecHitOccupancyProEta_[1]->Reset();
185  if ( meEERecHitOccupancyProPhi_[1] ) meEERecHitOccupancyProPhi_[1]->Reset();
186 
187  if ( meEERecHitOccupancyThr_[0] ) meEERecHitOccupancyThr_[0]->Reset();
188  if ( meEERecHitOccupancyProEtaThr_[0] ) meEERecHitOccupancyProEtaThr_[0]->Reset();
189  if ( meEERecHitOccupancyProPhiThr_[0] ) meEERecHitOccupancyProPhiThr_[0]->Reset();
190  if ( meEERecHitOccupancyThr_[1] ) meEERecHitOccupancyThr_[1]->Reset();
191  if ( meEERecHitOccupancyProEtaThr_[1] ) meEERecHitOccupancyProEtaThr_[1]->Reset();
192  if ( meEERecHitOccupancyProPhiThr_[1] ) meEERecHitOccupancyProPhiThr_[1]->Reset();
193 
194  if ( meEETrigPrimDigiOccupancy_[0] ) meEETrigPrimDigiOccupancy_[0]->Reset();
195  if ( meEETrigPrimDigiOccupancyProEta_[0] ) meEETrigPrimDigiOccupancyProEta_[0]->Reset();
196  if ( meEETrigPrimDigiOccupancyProPhi_[0] ) meEETrigPrimDigiOccupancyProPhi_[0]->Reset();
197  if ( meEETrigPrimDigiOccupancy_[1] ) meEETrigPrimDigiOccupancy_[1]->Reset();
198  if ( meEETrigPrimDigiOccupancyProEta_[1] ) meEETrigPrimDigiOccupancyProEta_[1]->Reset();
199  if ( meEETrigPrimDigiOccupancyProPhi_[1] ) meEETrigPrimDigiOccupancyProPhi_[1]->Reset();
200 
201  if ( meEETrigPrimDigiOccupancyThr_[0] ) meEETrigPrimDigiOccupancyThr_[0]->Reset();
202  if ( meEETrigPrimDigiOccupancyProEtaThr_[0] ) meEETrigPrimDigiOccupancyProEtaThr_[0]->Reset();
203  if ( meEETrigPrimDigiOccupancyProPhiThr_[0] ) meEETrigPrimDigiOccupancyProPhiThr_[0]->Reset();
204  if ( meEETrigPrimDigiOccupancyThr_[1] ) meEETrigPrimDigiOccupancyThr_[1]->Reset();
205  if ( meEETrigPrimDigiOccupancyProEtaThr_[1] ) meEETrigPrimDigiOccupancyProEtaThr_[1]->Reset();
206  if ( meEETrigPrimDigiOccupancyProPhiThr_[1] ) meEETrigPrimDigiOccupancyProPhiThr_[1]->Reset();
207 
208  if ( meEETestPulseDigiOccupancy_[0] ) meEETestPulseDigiOccupancy_[0]->Reset();
209  if ( meEETestPulseDigiOccupancy_[1] ) meEETestPulseDigiOccupancy_[1]->Reset();
210 
211  if ( meEELaserDigiOccupancy_[0] ) meEELaserDigiOccupancy_[0]->Reset();
212  if ( meEELaserDigiOccupancy_[1] ) meEELaserDigiOccupancy_[1]->Reset();
213 
214  if ( meEELedDigiOccupancy_[0] ) meEELedDigiOccupancy_[0]->Reset();
215  if ( meEELedDigiOccupancy_[1] ) meEELedDigiOccupancy_[1]->Reset();
216 
217  if ( meEEPedestalDigiOccupancy_[0] ) meEEPedestalDigiOccupancy_[0]->Reset();
218  if ( meEEPedestalDigiOccupancy_[1] ) meEEPedestalDigiOccupancy_[1]->Reset();
219 
220 }
221 
222 void EEOccupancyTask::setup(void){
223 
224  init_ = true;
225 
226  std::string name;
227 
228  if ( dqmStore_ ) {
229  dqmStore_->setCurrentFolder(prefixME_ + "/EEOccupancyTask");
230  if(subfolder_.size())
231  dqmStore_->setCurrentFolder(prefixME_ + "/EEOccupancyTask/" + subfolder_);
232 
233  for (int i = 0; i < 18; i++) {
234  name = "EEOT digi occupancy " + Numbers::sEE(i+1);
235  meOccupancy_[i] = dqmStore_->book2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50.);
236  meOccupancy_[i]->setAxisTitle("ix", 1);
237  if ( i+1 >= 1 && i+1 <= 9 ) meOccupancy_[i]->setAxisTitle("101-ix", 1);
238  meOccupancy_[i]->setAxisTitle("iy", 2);
239  dqmStore_->tag(meOccupancy_[i], i+1);
240 
241  name = "EEOT MEM digi occupancy " + Numbers::sEE(i+1);
242  meOccupancyMem_[i] = dqmStore_->book2D(name, name, 10, 0., 10., 5, 0., 5.);
243  meOccupancyMem_[i]->setAxisTitle("pseudo-strip", 1);
244  meOccupancyMem_[i]->setAxisTitle("channel", 2);
245  dqmStore_->tag(meOccupancyMem_[i], i+1);
246 
247  name = "EEOT rec hit energy " + Numbers::sEE(i+1);
248  meEERecHitEnergy_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096., "s");
249  meEERecHitEnergy_[i]->setAxisTitle("ix", 1);
250  if ( i+1 >= 1 && i+1 <= 9 ) meEERecHitEnergy_[i]->setAxisTitle("101-ix", 1);
251  meEERecHitEnergy_[i]->setAxisTitle("iy", 2);
252  meEERecHitEnergy_[i]->setAxisTitle("energy (GeV)", 3);
253  dqmStore_->tag(meEERecHitEnergy_[i], i+1);
254 
255  name = "EEOT energy spectrum " + Numbers::sEE(i+1);
256  meSpectrum_[i] = dqmStore_->book1D(name, name, 100, 0., 1.5);
257  meSpectrum_[i]->setAxisTitle("energy (GeV)", 1);
258  dqmStore_->tag(meSpectrum_[i], i+1);
259  }
260 
261  name = "EEOT rec hit spectrum EE -";
262  meEERecHitSpectrum_[0] = dqmStore_->book1D(name, name, 100, 0., 10.);
263  meEERecHitSpectrum_[0]->setAxisTitle("energy (GeV)", 1);
264 
265  name = "EEOT rec hit spectrum EE +";
266  meEERecHitSpectrum_[1] = dqmStore_->book1D(name, name, 100, 0., 10.);
267  meEERecHitSpectrum_[1]->setAxisTitle("energy (GeV)", 1);
268 
269  name = "EEOT digi occupancy EE -";
270  meEEDigiOccupancy_[0] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
271  meEEDigiOccupancy_[0]->setAxisTitle("jx", 1);
272  meEEDigiOccupancy_[0]->setAxisTitle("jy", 2);
273  name = "EEOT digi occupancy EE - projection eta";
274  meEEDigiOccupancyProEta_[0] = dqmStore_->book1DD(name, name, 22, -3.0, -1.479);
275  meEEDigiOccupancyProEta_[0]->setAxisTitle("eta", 1);
276  meEEDigiOccupancyProEta_[0]->setAxisTitle("number of digis", 2);
277  name = "EEOT digi occupancy EE - projection phi";
278  meEEDigiOccupancyProPhi_[0] = dqmStore_->book1DD(name, name, 50, -M_PI, M_PI);
279  meEEDigiOccupancyProPhi_[0]->setAxisTitle("phi", 1);
280  meEEDigiOccupancyProPhi_[0]->setAxisTitle("number of digis", 2);
281 
282  name = "EEOT digi occupancy EE +";
283  meEEDigiOccupancy_[1] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
284  meEEDigiOccupancy_[1]->setAxisTitle("jx", 1);
285  meEEDigiOccupancy_[1]->setAxisTitle("jy", 2);
286  name = "EEOT digi occupancy EE + projection eta";
287  meEEDigiOccupancyProEta_[1] = dqmStore_->book1DD(name, name, 22, 1.479, 3.0);
288  meEEDigiOccupancyProEta_[1]->setAxisTitle("eta", 1);
289  meEEDigiOccupancyProEta_[1]->setAxisTitle("number of digis", 2);
290  name = "EEOT digi occupancy EE + projection phi";
291  meEEDigiOccupancyProPhi_[1] = dqmStore_->book1DD(name, name, 50, -M_PI, M_PI);
292  meEEDigiOccupancyProPhi_[1]->setAxisTitle("phi", 1);
293  meEEDigiOccupancyProPhi_[1]->setAxisTitle("number of digis", 2);
294 
295  name = "EEOT rec hit occupancy EE -";
296  meEERecHitOccupancy_[0] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
297  meEERecHitOccupancy_[0]->setAxisTitle("jx", 1);
298  meEERecHitOccupancy_[0]->setAxisTitle("jy", 2);
299  name = "EEOT rec hit occupancy EE - projection eta";
300  meEERecHitOccupancyProEta_[0] = dqmStore_->book1DD(name, name, 22, -3.0, -1.479);
301  meEERecHitOccupancyProEta_[0]->setAxisTitle("eta", 1);
302  meEERecHitOccupancyProEta_[0]->setAxisTitle("number of hits", 2);
303  name = "EEOT rec hit occupancy EE - projection phi";
304  meEERecHitOccupancyProPhi_[0] = dqmStore_->book1DD(name, name, 50, -M_PI, M_PI);
305  meEERecHitOccupancyProPhi_[0]->setAxisTitle("phi", 1);
306  meEERecHitOccupancyProPhi_[0]->setAxisTitle("number of hits", 2);
307 
308  name = "EEOT rec hit occupancy EE +";
309  meEERecHitOccupancy_[1] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
310  meEERecHitOccupancy_[1]->setAxisTitle("jx", 1);
311  meEERecHitOccupancy_[1]->setAxisTitle("jy", 2);
312  name = "EEOT rec hit occupancy EE + projection eta";
313  meEERecHitOccupancyProEta_[1] = dqmStore_->book1DD(name, name, 22, 1.479, 3.0);
314  meEERecHitOccupancyProEta_[1]->setAxisTitle("eta", 1);
315  meEERecHitOccupancyProEta_[1]->setAxisTitle("number of hits", 2);
316  name = "EEOT rec hit occupancy EE + projection phi";
317  meEERecHitOccupancyProPhi_[1] = dqmStore_->book1DD(name, name, 50, -M_PI, M_PI);
318  meEERecHitOccupancyProPhi_[1]->setAxisTitle("phi", 1);
319  meEERecHitOccupancyProPhi_[1]->setAxisTitle("number of hits", 2);
320 
321  name = "EEOT rec hit thr occupancy EE -";
322  meEERecHitOccupancyThr_[0] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
323  meEERecHitOccupancyThr_[0]->setAxisTitle("jx", 1);
324  meEERecHitOccupancyThr_[0]->setAxisTitle("jy", 2);
325  name = "EEOT rec hit thr occupancy EE - projection eta";
326  meEERecHitOccupancyProEtaThr_[0] = dqmStore_->book1DD(name, name, 22, -3.0, -1.479);
327  meEERecHitOccupancyProEtaThr_[0]->setAxisTitle("eta", 1);
328  meEERecHitOccupancyProEtaThr_[0]->setAxisTitle("number of hits", 2);
329  name = "EEOT rec hit thr occupancy EE - projection phi";
330  meEERecHitOccupancyProPhiThr_[0] = dqmStore_->book1DD(name, name, 50, -M_PI, M_PI);
331  meEERecHitOccupancyProPhiThr_[0]->setAxisTitle("phi", 1);
332  meEERecHitOccupancyProPhiThr_[0]->setAxisTitle("number of hits", 2);
333 
334  name = "EEOT rec hit thr occupancy EE +";
335  meEERecHitOccupancyThr_[1] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
336  meEERecHitOccupancyThr_[1]->setAxisTitle("jx", 1);
337  meEERecHitOccupancyThr_[1]->setAxisTitle("jy", 2);
338  name = "EEOT rec hit thr occupancy EE + projection eta";
339  meEERecHitOccupancyProEtaThr_[1] = dqmStore_->book1DD(name, name, 22, 1.479, 3.0);
340  meEERecHitOccupancyProEtaThr_[1]->setAxisTitle("eta", 1);
341  meEERecHitOccupancyProEtaThr_[1]->setAxisTitle("number of hits", 2);
342  name = "EEOT rec hit thr occupancy EE + projection phi";
343  meEERecHitOccupancyProPhiThr_[1] = dqmStore_->book1DD(name, name, 50, -M_PI, M_PI);
344  meEERecHitOccupancyProPhiThr_[1]->setAxisTitle("phi", 1);
345  meEERecHitOccupancyProPhiThr_[1]->setAxisTitle("number of hits", 2);
346 
347  name = "EEOT TP digi occupancy EE -";
348  meEETrigPrimDigiOccupancy_[0] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
349  meEETrigPrimDigiOccupancy_[0]->setAxisTitle("jx", 1);
350  meEETrigPrimDigiOccupancy_[0]->setAxisTitle("jy", 2);
351  name = "EEOT TP digi occupancy EE - projection eta";
352  meEETrigPrimDigiOccupancyProEta_[0] = dqmStore_->book1DD(name, name, 22, -3.0, -1.479);
353  meEETrigPrimDigiOccupancyProEta_[0]->setAxisTitle("eta", 1);
354  meEETrigPrimDigiOccupancyProEta_[0]->setAxisTitle("number of TP digis", 2);
355  name = "EEOT TP digi occupancy EE - projection phi";
356  meEETrigPrimDigiOccupancyProPhi_[0] = dqmStore_->book1DD(name, name, 50, -M_PI, M_PI);
357  meEETrigPrimDigiOccupancyProPhi_[0]->setAxisTitle("phi", 1);
358  meEETrigPrimDigiOccupancyProPhi_[0]->setAxisTitle("number of TP digis", 2);
359 
360  name = "EEOT TP digi occupancy EE +";
361  meEETrigPrimDigiOccupancy_[1] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
362  meEETrigPrimDigiOccupancy_[1]->setAxisTitle("jx", 1);
363  meEETrigPrimDigiOccupancy_[1]->setAxisTitle("jy", 2);
364  name = "EEOT TP digi occupancy EE + projection eta";
365  meEETrigPrimDigiOccupancyProEta_[1] = dqmStore_->book1DD(name, name, 22, 1.479, 3.0);
366  meEETrigPrimDigiOccupancyProEta_[1]->setAxisTitle("eta", 1);
367  meEETrigPrimDigiOccupancyProEta_[1]->setAxisTitle("number of TP digis", 2);
368  name = "EEOT TP digi occupancy EE + projection phi";
369  meEETrigPrimDigiOccupancyProPhi_[1] = dqmStore_->book1DD(name, name, 50, -M_PI, M_PI);
370  meEETrigPrimDigiOccupancyProPhi_[1]->setAxisTitle("phi", 1);
371  meEETrigPrimDigiOccupancyProPhi_[1]->setAxisTitle("number of TP digis", 2);
372 
373  name = "EEOT TP digi thr occupancy EE -";
374  meEETrigPrimDigiOccupancyThr_[0] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
375  meEETrigPrimDigiOccupancyThr_[0]->setAxisTitle("jx", 1);
376  meEETrigPrimDigiOccupancyThr_[0]->setAxisTitle("jy", 2);
377  name = "EEOT TP digi thr occupancy EE - projection eta";
378  meEETrigPrimDigiOccupancyProEtaThr_[0] = dqmStore_->book1DD(name, name, 22, -3.0, -1.479);
379  meEETrigPrimDigiOccupancyProEtaThr_[0]->setAxisTitle("eta", 1);
380  meEETrigPrimDigiOccupancyProEtaThr_[0]->setAxisTitle("number of TP digis", 2);
381  name = "EEOT TP digi thr occupancy EE - projection phi";
382  meEETrigPrimDigiOccupancyProPhiThr_[0] = dqmStore_->book1DD(name, name, 50, -M_PI, M_PI);
383  meEETrigPrimDigiOccupancyProPhiThr_[0]->setAxisTitle("phi", 1);
384  meEETrigPrimDigiOccupancyProPhiThr_[0]->setAxisTitle("number of TP digis", 2);
385 
386  name = "EEOT TP digi thr occupancy EE +";
387  meEETrigPrimDigiOccupancyThr_[1] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
388  meEETrigPrimDigiOccupancyThr_[1]->setAxisTitle("jx", 1);
389  meEETrigPrimDigiOccupancyThr_[1]->setAxisTitle("jy", 2);
390  name = "EEOT TP digi thr occupancy EE + projection eta";
391  meEETrigPrimDigiOccupancyProEtaThr_[1] = dqmStore_->book1DD(name, name, 22, 1.479, 3.0);
392  meEETrigPrimDigiOccupancyProEtaThr_[1]->setAxisTitle("eta", 1);
393  meEETrigPrimDigiOccupancyProEtaThr_[1]->setAxisTitle("number of TP digis", 2);
394  name = "EEOT TP digi thr occupancy EE + projection phi";
395  meEETrigPrimDigiOccupancyProPhiThr_[1] = dqmStore_->book1DD(name, name, 50, -M_PI, M_PI);
396  meEETrigPrimDigiOccupancyProPhiThr_[1]->setAxisTitle("phi", 1);
397  meEETrigPrimDigiOccupancyProPhiThr_[1]->setAxisTitle("number of TP digis", 2);
398 
399  name = "EEOT test pulse digi occupancy EE -";
400  meEETestPulseDigiOccupancy_[0] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
401  meEETestPulseDigiOccupancy_[0]->setAxisTitle("jx", 1);
402  meEETestPulseDigiOccupancy_[0]->setAxisTitle("jy", 2);
403 
404  name = "EEOT test pulse digi occupancy EE +";
405  meEETestPulseDigiOccupancy_[1] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
406  meEETestPulseDigiOccupancy_[1]->setAxisTitle("jx", 1);
407  meEETestPulseDigiOccupancy_[1]->setAxisTitle("jy", 2);
408 
409  name = "EEOT led digi occupancy EE -";
410  meEELedDigiOccupancy_[0] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
411  meEELedDigiOccupancy_[0]->setAxisTitle("jx", 1);
412  meEELedDigiOccupancy_[0]->setAxisTitle("jy", 2);
413 
414  name = "EEOT led digi occupancy EE +";
415  meEELedDigiOccupancy_[1] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
416  meEELedDigiOccupancy_[1]->setAxisTitle("jx", 1);
417  meEELedDigiOccupancy_[1]->setAxisTitle("jy", 2);
418 
419  name = "EEOT laser digi occupancy EE -";
420  meEELaserDigiOccupancy_[0] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
421  meEELaserDigiOccupancy_[0]->setAxisTitle("jx", 1);
422  meEELaserDigiOccupancy_[0]->setAxisTitle("jy", 2);
423 
424  name = "EEOT laser digi occupancy EE +";
425  meEELaserDigiOccupancy_[1] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
426  meEELaserDigiOccupancy_[1]->setAxisTitle("jx", 1);
427  meEELaserDigiOccupancy_[1]->setAxisTitle("jy", 2);
428 
429  name = "EEOT pedestal digi occupancy EE -";
430  meEEPedestalDigiOccupancy_[0] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
431  meEEPedestalDigiOccupancy_[0]->setAxisTitle("jx", 1);
432  meEEPedestalDigiOccupancy_[0]->setAxisTitle("jy", 2);
433 
434  name = "EEOT pedestal digi occupancy EE +";
435  meEEPedestalDigiOccupancy_[1] = dqmStore_->book2D(name, name, 100, 0., 100., 100, 0., 100.);
436  meEEPedestalDigiOccupancy_[1]->setAxisTitle("jx", 1);
437  meEEPedestalDigiOccupancy_[1]->setAxisTitle("jy", 2);
438 
439  }
440 
441 }
442 
443 void EEOccupancyTask::cleanup(void){
444 
445  if ( ! init_ ) return;
446 
447  if ( dqmStore_ ) {
448  dqmStore_->setCurrentFolder(prefixME_ + "/EEOccupancyTask");
449  if(subfolder_.size())
450  dqmStore_->setCurrentFolder(prefixME_ + "/EEOccupancyTask/" + subfolder_);
451 
452  for (int i = 0; i < 18; i++) {
453  if ( meOccupancy_[i] ) dqmStore_->removeElement( meOccupancy_[i]->getName() );
454  meOccupancy_[i] = 0;
455  if ( meOccupancyMem_[i] ) dqmStore_->removeElement( meOccupancyMem_[i]->getName() );
456  meOccupancyMem_[i] = 0;
457  if ( meEERecHitEnergy_[i] ) dqmStore_->removeElement( meEERecHitEnergy_[i]->getName() );
458  meEERecHitEnergy_[i] = 0;
459  if ( meSpectrum_[i] ) dqmStore_->removeElement( meSpectrum_[i]->getName() );
460  meSpectrum_[i] = 0;
461  }
462 
463  if ( meEERecHitSpectrum_[0] ) dqmStore_->removeElement( meEERecHitSpectrum_[0]->getName() );
464  meEERecHitSpectrum_[0] = 0;
465  if ( meEERecHitSpectrum_[1] ) dqmStore_->removeElement( meEERecHitSpectrum_[1]->getName() );
466  meEERecHitSpectrum_[1] = 0;
467 
468  if ( meEEDigiOccupancy_[0] ) dqmStore_->removeElement( meEEDigiOccupancy_[0]->getName() );
469  meEEDigiOccupancy_[0] = 0;
470  if ( meEEDigiOccupancyProEta_[0] ) dqmStore_->removeElement( meEEDigiOccupancyProEta_[0]->getName() );
471  meEEDigiOccupancyProEta_[0] = 0;
472  if ( meEEDigiOccupancyProPhi_[0] ) dqmStore_->removeElement( meEEDigiOccupancyProPhi_[0]->getName() );
473  meEEDigiOccupancyProPhi_[0] = 0;
474 
475  if ( meEEDigiOccupancy_[1] ) dqmStore_->removeElement( meEEDigiOccupancy_[1]->getName() );
476  meEEDigiOccupancy_[1] = 0;
477  if ( meEEDigiOccupancyProEta_[1] ) dqmStore_->removeElement( meEEDigiOccupancyProEta_[1]->getName() );
478  meEEDigiOccupancyProEta_[1] = 0;
479  if ( meEEDigiOccupancyProPhi_[1] ) dqmStore_->removeElement( meEEDigiOccupancyProPhi_[1]->getName() );
480  meEEDigiOccupancyProPhi_[1] = 0;
481 
482  if ( meEERecHitOccupancy_[0] ) dqmStore_->removeElement( meEERecHitOccupancy_[0]->getName() );
483  meEERecHitOccupancy_[0] = 0;
484  if ( meEERecHitOccupancyProEta_[0] ) dqmStore_->removeElement( meEERecHitOccupancyProEta_[0]->getName() );
485  meEERecHitOccupancyProEta_[0] = 0;
486  if ( meEERecHitOccupancyProPhi_[0] ) dqmStore_->removeElement( meEERecHitOccupancyProPhi_[0]->getName() );
487  meEERecHitOccupancyProPhi_[0] = 0;
488 
489  if ( meEERecHitOccupancy_[1] ) dqmStore_->removeElement( meEERecHitOccupancy_[1]->getName() );
490  meEERecHitOccupancy_[1] = 0;
491  if ( meEERecHitOccupancyProEta_[1] ) dqmStore_->removeElement( meEERecHitOccupancyProEta_[1]->getName() );
492  meEERecHitOccupancyProEta_[1] = 0;
493  if ( meEERecHitOccupancyProPhi_[1] ) dqmStore_->removeElement( meEERecHitOccupancyProPhi_[1]->getName() );
494  meEERecHitOccupancyProPhi_[1] = 0;
495 
496  if ( meEERecHitOccupancyThr_[0] ) dqmStore_->removeElement( meEERecHitOccupancyThr_[0]->getName() );
497  meEERecHitOccupancyThr_[0] = 0;
498  if ( meEERecHitOccupancyProEtaThr_[0] ) dqmStore_->removeElement( meEERecHitOccupancyProEtaThr_[0]->getName() );
499  meEERecHitOccupancyProEtaThr_[0] = 0;
500  if ( meEERecHitOccupancyProPhiThr_[0] ) dqmStore_->removeElement( meEERecHitOccupancyProPhiThr_[0]->getName() );
501  meEERecHitOccupancyProPhiThr_[0] = 0;
502 
503  if ( meEERecHitOccupancyThr_[1] ) dqmStore_->removeElement( meEERecHitOccupancyThr_[1]->getName() );
504  meEERecHitOccupancyThr_[1] = 0;
505  if ( meEERecHitOccupancyProEtaThr_[1] ) dqmStore_->removeElement( meEERecHitOccupancyProEtaThr_[1]->getName() );
506  meEERecHitOccupancyProEtaThr_[1] = 0;
507  if ( meEERecHitOccupancyProPhiThr_[1] ) dqmStore_->removeElement( meEERecHitOccupancyProPhiThr_[1]->getName() );
508  meEERecHitOccupancyProPhiThr_[1] = 0;
509 
510  if ( meEETrigPrimDigiOccupancy_[0] ) dqmStore_->removeElement( meEETrigPrimDigiOccupancy_[0]->getName() );
511  meEETrigPrimDigiOccupancy_[0] = 0;
512  if ( meEETrigPrimDigiOccupancyProEta_[0] ) dqmStore_->removeElement( meEETrigPrimDigiOccupancyProEta_[0]->getName() );
513  meEETrigPrimDigiOccupancyProEta_[0] = 0;
514  if ( meEETrigPrimDigiOccupancyProPhi_[0] ) dqmStore_->removeElement( meEETrigPrimDigiOccupancyProPhi_[0]->getName() );
515  meEETrigPrimDigiOccupancyProPhi_[0] = 0;
516 
517  if ( meEETrigPrimDigiOccupancy_[1] ) dqmStore_->removeElement( meEETrigPrimDigiOccupancy_[1]->getName() );
518  meEETrigPrimDigiOccupancy_[1] = 0;
519  if ( meEETrigPrimDigiOccupancyProEta_[1] ) dqmStore_->removeElement( meEETrigPrimDigiOccupancyProEta_[1]->getName() );
520  meEETrigPrimDigiOccupancyProEta_[1] = 0;
521  if ( meEETrigPrimDigiOccupancyProPhi_[1] ) dqmStore_->removeElement( meEETrigPrimDigiOccupancyProPhi_[1]->getName() );
522  meEETrigPrimDigiOccupancyProPhi_[1] = 0;
523 
524  if ( meEETrigPrimDigiOccupancyThr_[0] ) dqmStore_->removeElement( meEETrigPrimDigiOccupancyThr_[0]->getName() );
525  meEETrigPrimDigiOccupancyThr_[0] = 0;
526  if ( meEETrigPrimDigiOccupancyProEtaThr_[0] ) dqmStore_->removeElement( meEETrigPrimDigiOccupancyProEtaThr_[0]->getName() );
527  meEETrigPrimDigiOccupancyProEtaThr_[0] = 0;
528  if ( meEETrigPrimDigiOccupancyProPhiThr_[0] ) dqmStore_->removeElement( meEETrigPrimDigiOccupancyProPhiThr_[0]->getName() );
529  meEETrigPrimDigiOccupancyProPhiThr_[0] = 0;
530 
531  if ( meEETrigPrimDigiOccupancyThr_[1] ) dqmStore_->removeElement( meEETrigPrimDigiOccupancyThr_[1]->getName() );
532  meEETrigPrimDigiOccupancyThr_[1] = 0;
533  if ( meEETrigPrimDigiOccupancyProEtaThr_[1] ) dqmStore_->removeElement( meEETrigPrimDigiOccupancyProEtaThr_[1]->getName() );
534  meEETrigPrimDigiOccupancyProEtaThr_[1] = 0;
535  if ( meEETrigPrimDigiOccupancyProPhiThr_[1] ) dqmStore_->removeElement( meEETrigPrimDigiOccupancyProPhiThr_[1]->getName() );
536  meEETrigPrimDigiOccupancyProPhiThr_[1] = 0;
537 
538  if ( meEETestPulseDigiOccupancy_[0] ) dqmStore_->removeElement( meEETestPulseDigiOccupancy_[0]->getName() );
539  meEETestPulseDigiOccupancy_[0] = 0;
540  if ( meEETestPulseDigiOccupancy_[1] ) dqmStore_->removeElement( meEETestPulseDigiOccupancy_[1]->getName() );
541  meEETestPulseDigiOccupancy_[1] = 0;
542 
543  if ( meEELaserDigiOccupancy_[0] ) dqmStore_->removeElement( meEELaserDigiOccupancy_[0]->getName() );
544  meEELaserDigiOccupancy_[0] = 0;
545  if ( meEELaserDigiOccupancy_[1] ) dqmStore_->removeElement( meEELaserDigiOccupancy_[1]->getName() );
546  meEELaserDigiOccupancy_[1] = 0;
547 
548  if ( meEELedDigiOccupancy_[0] ) dqmStore_->removeElement( meEELedDigiOccupancy_[0]->getName() );
549  meEELedDigiOccupancy_[0] = 0;
550  if ( meEELedDigiOccupancy_[1] ) dqmStore_->removeElement( meEELedDigiOccupancy_[1]->getName() );
551  meEELedDigiOccupancy_[1] = 0;
552 
553  if ( meEEPedestalDigiOccupancy_[0] ) dqmStore_->removeElement( meEEPedestalDigiOccupancy_[0]->getName() );
554  meEEPedestalDigiOccupancy_[0] = 0;
555  if ( meEEPedestalDigiOccupancy_[1] ) dqmStore_->removeElement( meEEPedestalDigiOccupancy_[1]->getName() );
556  meEEPedestalDigiOccupancy_[1] = 0;
557 
558  }
559 
560  init_ = false;
561 
562 }
563 
564 void EEOccupancyTask::endJob(void) {
565 
566  edm::LogInfo("EEOccupancyTask") << "analyzed " << ievt_ << " events";
567 
568  if ( enableCleanup_ ) this->cleanup();
569 
570 }
571 
572 void EEOccupancyTask::analyze(const edm::Event& e, const edm::EventSetup& c){
573 
574  if ( ! init_ ) this->setup();
575 
576  ievt_++;
577 
578  int runType[18] = { notdata };
579 
580  edm::Handle<EcalRawDataCollection> dcchs;
581 
582  if ( e.getByLabel(EcalRawDataCollection_, dcchs) ) {
583 
584  for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
585 
586  if ( Numbers::subDet( *dcchItr ) != EcalEndcap ) continue;
587 
588  int ism = Numbers::iSM( *dcchItr, EcalEndcap );
589 
590  int runtype = dcchItr->getRunType();
591 
592  if ( runtype == EcalDCCHeaderBlock::COSMIC ||
593  runtype == EcalDCCHeaderBlock::MTCC ||
594  runtype == EcalDCCHeaderBlock::COSMICS_GLOBAL ||
595  runtype == EcalDCCHeaderBlock::PHYSICS_GLOBAL ||
596  runtype == EcalDCCHeaderBlock::COSMICS_LOCAL ||
597  runtype == EcalDCCHeaderBlock::PHYSICS_LOCAL ) runType[ism-1] = physics;
598  if ( runtype == EcalDCCHeaderBlock::TESTPULSE_MGPA ||
599  runtype == EcalDCCHeaderBlock::TESTPULSE_GAP ) runType[ism-1] = testpulse;
600  if ( runtype == EcalDCCHeaderBlock::LASER_STD ||
601  runtype == EcalDCCHeaderBlock::LASER_GAP ) runType[ism-1] = laser;
602  if ( runtype == EcalDCCHeaderBlock::LED_STD ||
603  runtype == EcalDCCHeaderBlock::LED_GAP ) runType[ism-1] = led;
604  if ( runtype == EcalDCCHeaderBlock::PEDESTAL_STD ||
605  runtype == EcalDCCHeaderBlock::PEDESTAL_GAP ) runType[ism-1] = pedestal;
606 
607  }
608 
609  } else {
610  edm::LogWarning("EEOccupancyTask") << EcalRawDataCollection_ << " not available";
611  }
612 
613  edm::Handle<EEDigiCollection> digis;
614 
615  if ( e.getByLabel(EEDigiCollection_, digis) ) {
616 
617  int need = digis->size();
618  LogDebug("EEOccupancyTask") << "event " << ievt_ << " digi collection size " << need;
619 
620  for ( EEDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr ) {
621 
622  EEDetId id = digiItr->id();
623 
624  int ix = id.ix();
625  int iy = id.iy();
626 
627  int hi = id.hashedIndex();
628 
629  if ( geometryEE[hi][0] == 0 ) {
630  const GlobalPoint& pos = pGeometry_->getGeometry(id)->getPosition();
631  geometryEE[hi][0] = pos.eta();
632  geometryEE[hi][1] = pos.phi();
633  }
634 
635  float eta = geometryEE[hi][0];
636  float phi = geometryEE[hi][1];
637 
638  int ism = Numbers::iSM( id );
639 
640  if ( ism >= 1 && ism <= 9 ) ix = 101 - ix;
641 
642  float xix = ix - 0.5;
643  float xiy = iy - 0.5;
644 
645  if ( xix <= 0. || xix >= 100. || xiy <= 0. || xiy >= 100. ) {
646  edm::LogWarning("EEOccupancyTask") << " det id = " << id;
647  edm::LogWarning("EEOccupancyTask") << " sm, ix, iw " << ism << " " << ix << " " << iy;
648  edm::LogWarning("EEOccupancyTask") << " xix, xiy " << xix << " " << xiy;
649  }
650 
651  if ( meOccupancy_[ism-1] ) meOccupancy_[ism-1]->Fill( xix, xiy );
652 
653  int eex = id.ix();
654  int eey = id.iy();
655 
656  float xeex = eex - 0.5;
657  float xeey = eey - 0.5;
658 
659  if ( runType[ism-1] == physics || runType[ism-1] == notdata ) {
660 
661  if ( ism >=1 && ism <= 9 ) {
662  if ( meEEDigiOccupancy_[0] ) meEEDigiOccupancy_[0]->Fill( xeex, xeey );
663  if ( meEEDigiOccupancyProEta_[0] ) meEEDigiOccupancyProEta_[0]->Fill( eta );
664  if ( meEEDigiOccupancyProPhi_[0] ) meEEDigiOccupancyProPhi_[0]->Fill( phi );
665  } else {
666  if ( meEEDigiOccupancy_[1] ) meEEDigiOccupancy_[1]->Fill( xeex, xeey );
667  if ( meEEDigiOccupancyProEta_[1] ) meEEDigiOccupancyProEta_[1]->Fill( eta );
668  if ( meEEDigiOccupancyProPhi_[1] ) meEEDigiOccupancyProPhi_[1]->Fill( phi );
669  }
670 
671  }
672 
673  if ( runType[ism-1] == testpulse ) {
674 
675  if ( ism >=1 && ism <= 9 ) {
676  if ( meEETestPulseDigiOccupancy_[0] ) meEETestPulseDigiOccupancy_[0]->Fill( xeex, xeey );
677  } else {
678  if ( meEETestPulseDigiOccupancy_[1] ) meEETestPulseDigiOccupancy_[1]->Fill( xeex, xeey );
679  }
680 
681  }
682 
683  if ( runType[ism-1] == laser ) {
684 
685  if ( ism >=1 && ism <= 9 ) {
686  if ( meEELaserDigiOccupancy_[0] ) meEELaserDigiOccupancy_[0]->Fill( xeex, xeey );
687  } else {
688  if ( meEELaserDigiOccupancy_[1] ) meEELaserDigiOccupancy_[1]->Fill( xeex, xeey );
689  }
690 
691  }
692 
693  if ( runType[ism-1] == led ) {
694 
695  if ( ism >=1 && ism <= 9 ) {
696  if ( meEELedDigiOccupancy_[0] ) meEELedDigiOccupancy_[0]->Fill( xeex, xeey );
697  } else {
698  if ( meEELedDigiOccupancy_[1] ) meEELedDigiOccupancy_[1]->Fill( xeex, xeey );
699  }
700 
701  }
702 
703  if ( runType[ism-1] == pedestal ) {
704 
705  if ( ism >=1 && ism <= 9 ) {
706  if ( meEEPedestalDigiOccupancy_[0] ) meEEPedestalDigiOccupancy_[0]->Fill( xeex, xeey );
707  } else {
708  if ( meEEPedestalDigiOccupancy_[1] ) meEEPedestalDigiOccupancy_[1]->Fill( xeex, xeey );
709  }
710 
711  }
712 
713  }
714 
715  } else {
716 
717  edm::LogWarning("EEOccupancyTask") << EEDigiCollection_ << " not available";
718 
719  }
720 
721  edm::Handle<EcalPnDiodeDigiCollection> PNs;
722 
723  if ( e.getByLabel(EcalPnDiodeDigiCollection_, PNs) ) {
724 
725  // filling mem occupancy only for the 5 channels belonging
726  // to a fully reconstructed PN's
727 
728  for ( EcalPnDiodeDigiCollection::const_iterator pnItr = PNs->begin(); pnItr != PNs->end(); ++pnItr ) {
729 
730  if ( Numbers::subDet( pnItr->id() ) != EcalEndcap ) continue;
731 
732  int ism = Numbers::iSM( pnItr->id() );
733 
734  float PnId = pnItr->id().iPnId();
735 
736  PnId = PnId - 0.5;
737  float st = 0.0;
738 
739  for (int chInStrip = 1; chInStrip <= 5; chInStrip++){
740  if ( meOccupancyMem_[ism-1] ) {
741  st = chInStrip - 0.5;
742  meOccupancyMem_[ism-1]->Fill(PnId, st);
743  }
744  }
745 
746  }
747 
748  } else {
749 
750  edm::LogWarning("EEOccupancyTask") << EcalPnDiodeDigiCollection_ << " not available";
751 
752  }
753 
754  edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
755  c.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
756 
757  edm::Handle<EcalRecHitCollection> rechits;
758 
759  if ( e.getByLabel(EcalRecHitCollection_, rechits) ) {
760 
761  int nebrh = rechits->size();
762  LogDebug("EEOccupancyTask") << "event " << ievt_ << " rec hits collection size " << nebrh;
763 
764  for ( EcalRecHitCollection::const_iterator rechitItr = rechits->begin(); rechitItr != rechits->end(); ++rechitItr ) {
765 
766  EEDetId id = rechitItr->id();
767 
768  int eex = id.ix();
769  int eey = id.iy();
770 
771  int hi = id.hashedIndex();
772 
773  if ( geometryEE[hi][0] == 0 ) {
774  const GlobalPoint& pos = pGeometry_->getGeometry(id)->getPosition();
775  geometryEE[hi][0] = pos.eta();
776  geometryEE[hi][1] = pos.phi();
777  }
778 
779  float eta = geometryEE[hi][0];
780  float phi = geometryEE[hi][1];
781 
782  int ism = Numbers::iSM( id );
783 
784  // sector view (from electronics)
785  float xix = ( ism >= 1 && ism <= 9 ) ? 101 - eex - 0.5 : eex - 0.5;
786  float xiy = eey - 0.5;
787 
788  // physics view (from IP)
789  float xeex = eex - 0.5;
790  float xeey = eey - 0.5;
791 
792  if ( runType[ism-1] == physics || runType[ism-1] == notdata ) {
793 
794  if ( ism >= 1 && ism <= 9 ) {
795  if ( meEERecHitOccupancy_[0] ) meEERecHitOccupancy_[0]->Fill( xeex, xeey );
796  if ( meEERecHitOccupancyProEta_[0] ) meEERecHitOccupancyProEta_[0]->Fill( eta );
797  if ( meEERecHitOccupancyProPhi_[0] ) meEERecHitOccupancyProPhi_[0]->Fill( phi );
798  } else {
799  if ( meEERecHitOccupancy_[1] ) meEERecHitOccupancy_[1]->Fill( xeex, xeey );
800  if ( meEERecHitOccupancyProEta_[1] ) meEERecHitOccupancyProEta_[1]->Fill( eta );
801  if ( meEERecHitOccupancyProPhi_[1] ) meEERecHitOccupancyProPhi_[1]->Fill( phi );
802  }
803 
804  uint32_t flag = rechitItr->recoFlag();
805 
806  uint32_t sev = sevlv->severityLevel(id, *rechits);
807 
808  if ( rechitItr->energy() > recHitEnergyMin_ && flag == EcalRecHit::kGood && sev == EcalSeverityLevel::kGood ) {
809 
810  if ( ism >= 1 && ism <= 9 ) {
811  if ( meEERecHitOccupancyThr_[0] ) meEERecHitOccupancyThr_[0]->Fill( xeex, xeey );
812  if ( meEERecHitOccupancyProEtaThr_[0] ) meEERecHitOccupancyProEtaThr_[0]->Fill( eta );
813  if ( meEERecHitOccupancyProPhiThr_[0] ) meEERecHitOccupancyProPhiThr_[0]->Fill( phi );
814  } else {
815  if ( meEERecHitOccupancyThr_[1] ) meEERecHitOccupancyThr_[1]->Fill( xeex, xeey );
816  if ( meEERecHitOccupancyProEtaThr_[1] ) meEERecHitOccupancyProEtaThr_[1]->Fill( eta );
817  if ( meEERecHitOccupancyProPhiThr_[1] ) meEERecHitOccupancyProPhiThr_[1]->Fill( phi );
818  }
819 
820  }
821 
822  if ( flag == EcalRecHit::kGood && sev == EcalSeverityLevel::kGood ) {
823  if ( meEERecHitEnergy_[ism-1] ) meEERecHitEnergy_[ism-1]->Fill( xix, xiy, rechitItr->energy() );
824  if ( meSpectrum_[ism-1] ) meSpectrum_[ism-1]->Fill( rechitItr->energy() );
825  if ( ism >= 1 && ism <= 9 ) meEERecHitSpectrum_[0]->Fill( rechitItr->energy() );
826  else meEERecHitSpectrum_[1]->Fill( rechitItr->energy() );
827  }
828 
829  }
830  }
831 
832  } else {
833 
834  edm::LogWarning("EEOccupancyTask") << EcalRecHitCollection_ << " not available";
835 
836  }
837 
838  edm::Handle<EcalTrigPrimDigiCollection> trigPrimDigis;
839 
840  if ( e.getByLabel(EcalTrigPrimDigiCollection_, trigPrimDigis) ) {
841 
842  int nebtpg = trigPrimDigis->size();
843  LogDebug("EEOccupancyTask") << "event " << ievt_ << " trigger primitives digis collection size " << nebtpg;
844 
845  for ( EcalTrigPrimDigiCollection::const_iterator tpdigiItr = trigPrimDigis->begin(); tpdigiItr != trigPrimDigis->end(); ++tpdigiItr ) {
846 
847  if ( Numbers::subDet( tpdigiItr->id() ) != EcalEndcap ) continue;
848 
849  int ism = Numbers::iSM( tpdigiItr->id() );
850 
851  std::vector<DetId>* crystals = Numbers::crystals( tpdigiItr->id() );
852 
853  for ( unsigned int i=0; i<crystals->size(); i++ ) {
854 
855  EEDetId id = (*crystals)[i];
856 
857  int eex = id.ix();
858  int eey = id.iy();
859 
860  int hi = id.hashedIndex();
861 
862  if ( geometryEE[hi][0] == 0 ) {
863  const GlobalPoint& pos = pGeometry_->getGeometry(id)->getPosition();
864  geometryEE[hi][0] = pos.eta();
865  geometryEE[hi][1] = pos.phi();
866  }
867 
868  float eta = geometryEE[hi][0];
869  float phi = geometryEE[hi][1];
870 
871  float xeex = eex - 0.5;
872  float xeey = eey - 0.5;
873 
874  if ( runType[ism-1] == physics || runType[ism-1] == notdata ) {
875 
876  if ( ism >= 1 && ism <= 9 ) {
877  if ( meEETrigPrimDigiOccupancy_[0] ) meEETrigPrimDigiOccupancy_[0]->Fill( xeex, xeey );
878  if ( meEETrigPrimDigiOccupancyProEta_[0] ) meEETrigPrimDigiOccupancyProEta_[0]->Fill( eta );
879  if ( meEETrigPrimDigiOccupancyProPhi_[0] ) meEETrigPrimDigiOccupancyProPhi_[0]->Fill( phi );
880  } else {
881  if ( meEETrigPrimDigiOccupancy_[1] ) meEETrigPrimDigiOccupancy_[1]->Fill( xeex, xeey );
882  if ( meEETrigPrimDigiOccupancyProEta_[1] ) meEETrigPrimDigiOccupancyProEta_[1]->Fill( eta );
883  if ( meEETrigPrimDigiOccupancyProPhi_[1] ) meEETrigPrimDigiOccupancyProPhi_[1]->Fill( phi );
884  }
885 
886  if ( tpdigiItr->compressedEt() > trigPrimEtMin_ ) {
887 
888  if ( ism >= 1 && ism <= 9 ) {
889  if ( meEETrigPrimDigiOccupancyThr_[0] ) meEETrigPrimDigiOccupancyThr_[0]->Fill( xeex, xeey );
890  if ( meEETrigPrimDigiOccupancyProEtaThr_[0] ) meEETrigPrimDigiOccupancyProEtaThr_[0]->Fill( eta );
891  if ( meEETrigPrimDigiOccupancyProPhiThr_[0] ) meEETrigPrimDigiOccupancyProPhiThr_[0]->Fill( phi );
892  } else {
893  if ( meEETrigPrimDigiOccupancyThr_[1] ) meEETrigPrimDigiOccupancyThr_[1]->Fill( xeex, xeey );
894  if ( meEETrigPrimDigiOccupancyProEtaThr_[1] ) meEETrigPrimDigiOccupancyProEtaThr_[1]->Fill( eta );
895  if ( meEETrigPrimDigiOccupancyProPhiThr_[1] ) meEETrigPrimDigiOccupancyProPhiThr_[1]->Fill( phi );
896  }
897 
898  }
899 
900  }
901  }
902  }
903 
904  } else {
905 
906  edm::LogWarning("EEOccupancyTask") << EcalTrigPrimDigiCollection_ << " not available";
907 
908  }
909 
910 }
911 
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
EEOccupancyTask::meEEPedestalDigiOccupancy_
MonitorElement * meEEPedestalDigiOccupancy_[2]
Definition: EEOccupancyTask.h:111
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
EEOccupancyTask::meEETrigPrimDigiOccupancyProEtaThr_
MonitorElement * meEETrigPrimDigiOccupancyProEtaThr_[2]
Definition: EEOccupancyTask.h:106
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
Numbers::sEE
static std::string sEE(const unsigned ism)
Definition: Numbers.cc:226
i
int i
Definition: DBlmapReader.cc:9
EEOccupancyTask::beginJob
void beginJob(void)
BeginJob.
Definition: EEOccupancyTask.cc:131
flag
long int flag
Definition: mlp_lapack.h:47
EEOccupancyTask::meEETestPulseDigiOccupancy_
MonitorElement * meEETestPulseDigiOccupancy_[2]
Definition: EEOccupancyTask.h:108
EEDetId::ix
int ix() const
Definition: EEDetId.h:71
edm::DataFrameContainer::const_iterator
boost::transform_iterator< IterHelp, boost::counting_iterator< int > > const_iterator
Definition: DataFrameContainer.h:51
EEOccupancyTask::meEERecHitEnergy_
MonitorElement * meEERecHitEnergy_[18]
Definition: EEOccupancyTask.h:89
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:717
DQMStore::rmdir
void rmdir(const std::string &fullpath)
Definition: DQMStore.cc:2530
EEOccupancyTask::setup
void setup(void)
Setup.
Definition: EEOccupancyTask.cc:222
Numbers.h
Some &quot;id&quot; conversions.
EEOccupancyTask::meEETrigPrimDigiOccupancyProPhi_
MonitorElement * meEETrigPrimDigiOccupancyProPhi_[2]
Definition: EEOccupancyTask.h:104
EEOccupancyTask::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
BeginRun.
Definition: EEOccupancyTask.cc:144
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:68
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:80
EEOccupancyTask::meEELaserDigiOccupancy_
MonitorElement * meEELaserDigiOccupancy_[2]
Definition: EEOccupancyTask.h:109
EEOccupancyTask::meEERecHitOccupancyProEtaThr_
MonitorElement * meEERecHitOccupancyProEtaThr_[2]
Definition: EEOccupancyTask.h:100
EEOccupancyTask::EcalPnDiodeDigiCollection_
edm::InputTag EcalPnDiodeDigiCollection_
Definition: EEOccupancyTask.h:80
EEOccupancyTask::meSpectrum_
MonitorElement * meSpectrum_[18]
Definition: EEOccupancyTask.h:90
EEOccupancyTask::endRun
void endRun(const edm::Run &r, const edm::EventSetup &c)
EndRun.
Definition: EEOccupancyTask.cc:157
EEOccupancyTask::meEEDigiOccupancyProEta_
MonitorElement * meEEDigiOccupancyProEta_[2]
Definition: EEOccupancyTask.h:94
Numbers::ix0EE
static int ix0EE(const unsigned ism)
Definition: Numbers.cc:773
DQMStore::book1DD
MonitorElement * book1DD(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1S histogram.
Definition: DQMStore.cc:749
EEOccupancyTask::pGeometry_
edm::ESHandle< CaloGeometry > pGeometry_
Definition: EEOccupancyTask.h:116
eta
T eta() const
Definition: Basic3DVectorLD.h:172
EEOccupancyTask::meEERecHitOccupancyProPhi_
MonitorElement * meEERecHitOccupancyProPhi_[2]
Definition: EEOccupancyTask.h:98
Numbers::iy0EE
static int iy0EE(const unsigned ism)
Definition: Numbers.cc:812
EEOccupancyTask::~EEOccupancyTask
virtual ~EEOccupancyTask()
Destructor.
Definition: EEOccupancyTask.cc:127
EEOccupancyTask::geometryEE
float geometryEE[EEDetId::kSizeForDenseIndexing][2]
Definition: EEOccupancyTask.h:118
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:152
EEOccupancyTask::endJob
void endJob(void)
EndJob.
Definition: EEOccupancyTask.cc:564
DQMStore::tag
void tag(MonitorElement *me, unsigned int myTag)
Definition: DQMStore.cc:1354
EEOccupancyTask::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Analyze.
Definition: EEOccupancyTask.cc:572
EEOccupancyTask::meOccupancyMem_
MonitorElement * meOccupancyMem_[18]
Definition: EEOccupancyTask.h:88
EEOccupancyTask::meEETrigPrimDigiOccupancyThr_
MonitorElement * meEETrigPrimDigiOccupancyThr_[2]
Definition: EEOccupancyTask.h:105
errorMatrix2Lands_multiChannel.id
tuple id
Definition: errorMatrix2Lands_multiChannel.py:151
EEOccupancyTask::EcalTrigPrimDigiCollection_
edm::InputTag EcalTrigPrimDigiCollection_
Definition: EEOccupancyTask.h:82
Numbers::crystals
static std::vector< DetId > * crystals(const EcalTrigTowerDetId &id)
Definition: Numbers.cc:585
EEOccupancyTask::EcalRawDataCollection_
edm::InputTag EcalRawDataCollection_
Definition: EEOccupancyTask.h:78
DQMStore::removeElement
void removeElement(const std::string &name)
Definition: DQMStore.cc:2572
edm::getName
std::string getName(Reflex::Type &cc)
Definition: ClassFiller.cc:18
EEOccupancyTask::meEERecHitOccupancyProPhiThr_
MonitorElement * meEERecHitOccupancyProPhiThr_[2]
Definition: EEOccupancyTask.h:101
EEOccupancyTask::meEERecHitOccupancy_
MonitorElement * meEERecHitOccupancy_[2]
Definition: EEOccupancyTask.h:96
EEOccupancyTask::meEEDigiOccupancy_
MonitorElement * meEEDigiOccupancy_[2]
Definition: EEOccupancyTask.h:93
EEOccupancyTask::meEETrigPrimDigiOccupancy_
MonitorElement * meEETrigPrimDigiOccupancy_[2]
Definition: EEOccupancyTask.h:102
EEOccupancyTask::EEOccupancyTask
EEOccupancyTask(const edm::ParameterSet &ps)
Constructor.
Definition: EEOccupancyTask.cc:38
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:356
EEOccupancyTask::meEERecHitOccupancyThr_
MonitorElement * meEERecHitOccupancyThr_[2]
Definition: EEOccupancyTask.h:99
EEOccupancyTask::meEEDigiOccupancyProPhi_
MonitorElement * meEEDigiOccupancyProPhi_[2]
Definition: EEOccupancyTask.h:95
environment_file_cff.runType
tuple runType
Definition: environment_file_cff.py:30
M_PI
#define M_PI
Definition: BFit3D.cc:3
Numbers::initGeometry
static void initGeometry(const edm::EventSetup &setup, bool verbose=false)
Definition: Numbers.cc:50
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
EEOccupancyTask::reset
void reset(void)
Reset.
Definition: EEOccupancyTask.cc:161
EEOccupancyTask::dqmStore_
DQMStore * dqmStore_
Definition: EEOccupancyTask.h:68
EEOccupancyTask::meEELedDigiOccupancy_
MonitorElement * meEELedDigiOccupancy_[2]
Definition: EEOccupancyTask.h:110
EEOccupancyTask::meEETrigPrimDigiOccupancyProEta_
MonitorElement * meEETrigPrimDigiOccupancyProEta_[2]
Definition: EEOccupancyTask.h:103
EEOccupancyTask::EcalRecHitCollection_
edm::InputTag EcalRecHitCollection_
Definition: EEOccupancyTask.h:81
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:75
EEOccupancyTask::meOccupancy_
MonitorElement * meOccupancy_[18]
Definition: EEOccupancyTask.h:87
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
EEOccupancyTask::meEERecHitOccupancyProEta_
MonitorElement * meEERecHitOccupancyProEta_[2]
Definition: EEOccupancyTask.h:97
EEOccupancyTask::prefixME_
std::string prefixME_
Definition: EEOccupancyTask.h:70
EEOccupancyTask::subfolder_
std::string subfolder_
Definition: EEOccupancyTask.h:72
EEOccupancyTask::EEDigiCollection_
edm::InputTag EEDigiCollection_
Definition: EEOccupancyTask.h:79
Numbers::subDet
static EcalSubdetector subDet(const EBDetId &id)
Definition: Numbers.cc:145
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
EEOccupancyTask::cleanup
void cleanup(void)
Cleanup.
Definition: EEOccupancyTask.cc:443
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
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:429
edm::Run
Definition: Run.h:33
EEOccupancyTask::meEETrigPrimDigiOccupancyProPhiThr_
MonitorElement * meEETrigPrimDigiOccupancyProPhiThr_[2]
Definition: EEOccupancyTask.h:107
EEOccupancyTask::meEERecHitSpectrum_
MonitorElement * meEERecHitSpectrum_[2]
Definition: EEOccupancyTask.h:92
phi
Definition: DDAxes.h:10
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