CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EELaserTask.cc
Go to the documentation of this file.
1 /*
2  * \file EELaserTask.cc
3  *
4  * \author G. Della Ricca
5  *
6 */
7 
8 #include <iostream>
9 #include <sstream>
10 #include <vector>
11 
12 #include "FWCore/ServiceRegistry/interface/Service.h"
13 #include "FWCore/MessageLogger/interface/MessageLogger.h"
14 
15 #include "DQMServices/Core/interface/MonitorElement.h"
16 
17 #include "DQMServices/Core/interface/DQMStore.h"
18 
19 #include "DataFormats/EcalDetId/interface/EEDetId.h"
20 #include "DataFormats/EcalDigi/interface/EEDataFrame.h"
21 #include "DataFormats/EcalRecHit/interface/EcalUncalibratedRecHit.h"
22 
23 #include "DQM/EcalCommon/interface/Numbers.h"
24 #include "DQM/EcalCommon/interface/NumbersPn.h"
25 
26 #include "DQM/EcalEndcapMonitorTasks/interface/EELaserTask.h"
27 
28 EELaserTask::EELaserTask(const edm::ParameterSet& ps){
29 
30  init_ = false;
31 
32  dqmStore_ = edm::Service<DQMStore>().operator->();
33 
34  prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
35 
36  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
37 
38  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
39 
40  EcalRawDataCollection_ = consumes<EcalRawDataCollection>(ps.getParameter<edm::InputTag>("EcalRawDataCollection"));
41  EEDigiCollection_ = consumes<EEDigiCollection>(ps.getParameter<edm::InputTag>("EEDigiCollection"));
42  EcalPnDiodeDigiCollection_ = consumes<EcalPnDiodeDigiCollection>(ps.getParameter<edm::InputTag>("EcalPnDiodeDigiCollection"));
43  EcalUncalibratedRecHitCollection_ = consumes<EcalUncalibratedRecHitCollection>(ps.getParameter<edm::InputTag>("EcalUncalibratedRecHitCollection"));
44 
45  // vector of enabled wavelengths (Default to all 4)
46  laserWavelengths_.reserve(4);
47  for ( unsigned int i = 1; i <= 4; i++ ) laserWavelengths_.push_back(i);
48  laserWavelengths_ = ps.getUntrackedParameter<std::vector<int> >("laserWavelengths", laserWavelengths_);
49 
50  for (int i = 0; i < 18; i++) {
51  meShapeMapL1_[i] = 0;
52  meAmplMapL1_[i] = 0;
53  meTimeMapL1_[i] = 0;
54  meAmplPNMapL1_[i] = 0;
55  mePnAmplMapG01L1_[i] = 0;
56  mePnPedMapG01L1_[i] = 0;
57  mePnAmplMapG16L1_[i] = 0;
58  mePnPedMapG16L1_[i] = 0;
59 
60  meShapeMapL2_[i] = 0;
61  meAmplMapL2_[i] = 0;
62  meTimeMapL2_[i] = 0;
63  meAmplPNMapL2_[i] = 0;
64  mePnAmplMapG01L2_[i] = 0;
65  mePnPedMapG01L2_[i] = 0;
66  mePnAmplMapG16L2_[i] = 0;
67  mePnPedMapG16L2_[i] = 0;
68 
69  meShapeMapL3_[i] = 0;
70  meAmplMapL3_[i] = 0;
71  meTimeMapL3_[i] = 0;
72  meAmplPNMapL3_[i] = 0;
73  mePnAmplMapG01L3_[i] = 0;
74  mePnPedMapG01L3_[i] = 0;
75  mePnAmplMapG16L3_[i] = 0;
76  mePnPedMapG16L3_[i] = 0;
77 
78  meShapeMapL4_[i] = 0;
79  meAmplMapL4_[i] = 0;
80  meTimeMapL4_[i] = 0;
81  meAmplPNMapL4_[i] = 0;
82  mePnAmplMapG01L4_[i] = 0;
83  mePnPedMapG01L4_[i] = 0;
84  mePnAmplMapG16L4_[i] = 0;
85  mePnPedMapG16L4_[i] = 0;
86  }
87 
88  for(int i=0; i<2; i++){
89  meAmplSummaryMapL1_[i] = 0;
90  meAmplSummaryMapL2_[i] = 0;
91  meAmplSummaryMapL3_[i] = 0;
92  meAmplSummaryMapL4_[i] = 0;
93  }
94 
95 }
96 
97 EELaserTask::~EELaserTask(){
98 
99 }
100 
101 void EELaserTask::beginJob(void){
102 
103  ievt_ = 0;
104 
105  if ( dqmStore_ ) {
106  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask");
107  dqmStore_->rmdir(prefixME_ + "/EELaserTask");
108  }
109 
110 }
111 
112 void EELaserTask::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 EELaserTask::endRun(const edm::Run& r, const edm::EventSetup& c) {
121 
122 }
123 
124 void EELaserTask::reset(void) {
125 
126  for (int i = 0; i < 18; i++) {
127  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
128  if ( meShapeMapL1_[i] ) meShapeMapL1_[i]->Reset();
129  if ( meAmplMapL1_[i] ) meAmplMapL1_[i]->Reset();
130  if ( meTimeMapL1_[i] ) meTimeMapL1_[i]->Reset();
131  if ( meAmplPNMapL1_[i] ) meAmplPNMapL1_[i]->Reset();
132  }
133 
134  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
135  if ( meShapeMapL2_[i] ) meShapeMapL2_[i]->Reset();
136  if ( meAmplMapL2_[i] ) meAmplMapL2_[i]->Reset();
137  if ( meTimeMapL2_[i] ) meTimeMapL2_[i]->Reset();
138  if ( meAmplPNMapL2_[i] ) meAmplPNMapL2_[i]->Reset();
139  }
140 
141  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
142  if ( meShapeMapL3_[i] ) meShapeMapL3_[i]->Reset();
143  if ( meAmplMapL3_[i] ) meAmplMapL3_[i]->Reset();
144  if ( meTimeMapL3_[i] ) meTimeMapL3_[i]->Reset();
145  if ( meAmplPNMapL3_[i] ) meAmplPNMapL3_[i]->Reset();
146  }
147 
148  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
149  if ( meShapeMapL4_[i] ) meShapeMapL4_[i]->Reset();
150  if ( meAmplMapL4_[i] ) meAmplMapL4_[i]->Reset();
151  if ( meTimeMapL4_[i] ) meTimeMapL4_[i]->Reset();
152  if ( meAmplPNMapL4_[i] ) meAmplPNMapL4_[i]->Reset();
153  }
154 
155  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
156  if ( mePnAmplMapG01L1_[i] ) mePnAmplMapG01L1_[i]->Reset();
157  if ( mePnPedMapG01L1_[i] ) mePnPedMapG01L1_[i]->Reset();
158 
159  if ( mePnAmplMapG16L1_[i] ) mePnAmplMapG16L1_[i]->Reset();
160  if ( mePnPedMapG16L1_[i] ) mePnPedMapG16L1_[i]->Reset();
161  }
162 
163 
164  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
165  if ( mePnAmplMapG01L2_[i] ) mePnAmplMapG01L2_[i]->Reset();
166  if ( mePnPedMapG01L2_[i] ) mePnPedMapG01L2_[i]->Reset();
167 
168  if ( mePnAmplMapG16L2_[i] ) mePnAmplMapG16L2_[i]->Reset();
169  if ( mePnPedMapG16L2_[i] ) mePnPedMapG16L2_[i]->Reset();
170  }
171 
172  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
173  if ( mePnAmplMapG01L3_[i] ) mePnAmplMapG01L3_[i]->Reset();
174  if ( mePnPedMapG01L3_[i] ) mePnPedMapG01L3_[i]->Reset();
175 
176  if ( mePnAmplMapG16L3_[i] ) mePnAmplMapG16L3_[i]->Reset();
177  if ( mePnPedMapG16L3_[i] ) mePnPedMapG16L3_[i]->Reset();
178  }
179 
180  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
181  if ( mePnAmplMapG01L4_[i] ) mePnAmplMapG01L4_[i]->Reset();
182  if ( mePnPedMapG01L4_[i] ) mePnPedMapG01L4_[i]->Reset();
183 
184  if ( mePnAmplMapG16L4_[i] ) mePnAmplMapG16L4_[i]->Reset();
185  if ( mePnPedMapG16L4_[i] ) mePnPedMapG16L4_[i]->Reset();
186  }
187 
188  }
189 
190  for(int i=0; i<2; i++){
191  if( meAmplSummaryMapL1_[i] ) meAmplSummaryMapL1_[i]->Reset();
192  if( meAmplSummaryMapL2_[i] ) meAmplSummaryMapL2_[i]->Reset();
193  if( meAmplSummaryMapL3_[i] ) meAmplSummaryMapL3_[i]->Reset();
194  if( meAmplSummaryMapL4_[i] ) meAmplSummaryMapL4_[i]->Reset();
195  }
196 
197 }
198 
199 void EELaserTask::setup(void){
200 
201  init_ = true;
202 
203  std::string name;
204  std::stringstream LaserN, LN;
205 
206  if ( dqmStore_ ) {
207  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask");
208 
209  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
210 
211  LaserN.str("");
212  LaserN << "Laser" << 1;
213  LN.str("");
214  LN << "L" << 1;
215 
216  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str());
217  for (int i = 0; i < 18; i++) {
218  name = "EELT shape " + Numbers::sEE(i+1) + " " + LN.str();
219  meShapeMapL1_[i] = dqmStore_->bookProfile2D(name, name, 850, 0., 850., 10, 0., 10., 4096, 0., 4096., "s");
220  meShapeMapL1_[i]->setAxisTitle("channel", 1);
221  meShapeMapL1_[i]->setAxisTitle("sample", 2);
222  meShapeMapL1_[i]->setAxisTitle("amplitude", 3);
223  dqmStore_->tag(meShapeMapL1_[i], i+1);
224 
225  name = "EELT amplitude " + Numbers::sEE(i+1) + " " + LN.str();
226  meAmplMapL1_[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.*12., "s");
227  meAmplMapL1_[i]->setAxisTitle("ix", 1);
228  if ( i+1 >= 1 && i+1 <= 9 ) meAmplMapL1_[i]->setAxisTitle("101-ix", 1);
229  meAmplMapL1_[i]->setAxisTitle("iy", 2);
230  dqmStore_->tag(meAmplMapL1_[i], i+1);
231 
232  name = "EELT timing " + Numbers::sEE(i+1) + " " + LN.str();
233  meTimeMapL1_[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., 250, 0., 10., "s");
234  meTimeMapL1_[i]->setAxisTitle("ix", 1);
235  if ( i+1 >= 1 && i+1 <= 9 ) meTimeMapL1_[i]->setAxisTitle("101-ix", 1);
236  meTimeMapL1_[i]->setAxisTitle("iy", 2);
237  dqmStore_->tag(meTimeMapL1_[i], i+1);
238 
239  name = "EELT amplitude over PN " + Numbers::sEE(i+1) + " " + LN.str();
240  meAmplPNMapL1_[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.*12., "s");
241  meAmplPNMapL1_[i]->setAxisTitle("ix", 1);
242  if ( i+1 >= 1 && i+1 <= 9 ) meAmplPNMapL1_[i]->setAxisTitle("101-ix", 1);
243  meAmplPNMapL1_[i]->setAxisTitle("iy", 2);
244  dqmStore_->tag(meAmplPNMapL1_[i], i+1);
245  }
246 
247  name = "EELT amplitude map " + LN.str() + " EE -";
248  meAmplSummaryMapL1_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
249  meAmplSummaryMapL1_[0]->setAxisTitle("ix", 1);
250  meAmplSummaryMapL1_[0]->setAxisTitle("iy", 2);
251 
252  name = "EELT amplitude map " + LN.str() + " EE +";
253  meAmplSummaryMapL1_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
254  meAmplSummaryMapL1_[1]->setAxisTitle("ix", 1);
255  meAmplSummaryMapL1_[1]->setAxisTitle("iy", 2);
256 
257  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN");
258 
259  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain01");
260 
261  for (int i = 0; i < 18; i++) {
262  name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G01 " + LN.str();
263  mePnAmplMapG01L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
264  mePnAmplMapG01L1_[i]->setAxisTitle("channel", 1);
265  mePnAmplMapG01L1_[i]->setAxisTitle("amplitude", 2);
266  dqmStore_->tag(mePnAmplMapG01L1_[i], i+1);
267 
268  name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G01 " + LN.str();
269  mePnPedMapG01L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
270  mePnPedMapG01L1_[i]->setAxisTitle("channel", 1);
271  mePnPedMapG01L1_[i]->setAxisTitle("pedestal", 2);
272  dqmStore_->tag(mePnPedMapG01L1_[i], i+1);
273  }
274 
275  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain16");
276 
277  for (int i = 0; i < 18; i++) {
278  name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G16 " + LN.str();
279  mePnAmplMapG16L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
280  mePnAmplMapG16L1_[i]->setAxisTitle("channel", 1);
281  mePnAmplMapG16L1_[i]->setAxisTitle("amplitude", 2);
282  dqmStore_->tag(mePnAmplMapG16L1_[i], i+1);
283 
284  name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G16 " + LN.str();
285  mePnPedMapG16L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
286  mePnPedMapG16L1_[i]->setAxisTitle("channel", 1);
287  mePnPedMapG16L1_[i]->setAxisTitle("pedestal", 2);
288  dqmStore_->tag(mePnPedMapG16L1_[i], i+1);
289  }
290 
291  }
292 
293  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
294 
295  LaserN.str("");
296  LaserN << "Laser" << 2;
297  LN.str("");
298  LN << "L" << 2;
299 
300  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str());
301  for (int i = 0; i < 18; i++) {
302  name = "EELT shape " + Numbers::sEE(i+1) + " " + LN.str();
303  meShapeMapL2_[i] = dqmStore_->bookProfile2D(name, name, 850, 0., 850., 10, 0., 10., 4096, 0., 4096., "s");
304  meShapeMapL2_[i]->setAxisTitle("channel", 1);
305  meShapeMapL2_[i]->setAxisTitle("sample", 2);
306  meShapeMapL2_[i]->setAxisTitle("amplitude", 3);
307  dqmStore_->tag(meShapeMapL2_[i], i+1);
308 
309  name = "EELT amplitude " + Numbers::sEE(i+1) + " " + LN.str();
310  meAmplMapL2_[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.*12., "s");
311  meAmplMapL2_[i]->setAxisTitle("ix", 1);
312  if ( i+1 >= 1 && i+1 <= 9 ) meAmplMapL2_[i]->setAxisTitle("101-ix", 1);
313  meAmplMapL2_[i]->setAxisTitle("iy", 2);
314  dqmStore_->tag(meAmplMapL2_[i], i+1);
315 
316  name = "EELT timing " + Numbers::sEE(i+1) + " " + LN.str();
317  meTimeMapL2_[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., 250, 0., 10., "s");
318  meTimeMapL2_[i]->setAxisTitle("ix", 1);
319  if ( i+1 >= 1 && i+1 <= 9 ) meTimeMapL2_[i]->setAxisTitle("101-ix", 1);
320  meTimeMapL2_[i]->setAxisTitle("iy", 2);
321  dqmStore_->tag(meTimeMapL2_[i], i+1);
322 
323  name = "EELT amplitude over PN " + Numbers::sEE(i+1) + " " + LN.str();
324  meAmplPNMapL2_[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.*12., "s");
325  meAmplPNMapL2_[i]->setAxisTitle("ix", 1);
326  if ( i+1 >= 1 && i+1 <= 9 ) meAmplPNMapL2_[i]->setAxisTitle("101-ix", 1);
327  meAmplPNMapL2_[i]->setAxisTitle("iy", 2);
328  dqmStore_->tag(meAmplPNMapL2_[i], i+1);
329  }
330 
331  name = "EELT amplitude map " + LN.str() + " EE -";
332  meAmplSummaryMapL2_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
333  meAmplSummaryMapL2_[0]->setAxisTitle("ix", 1);
334  meAmplSummaryMapL2_[0]->setAxisTitle("iy", 2);
335 
336  name = "EELT amplitude map " + LN.str() + " EE +";
337  meAmplSummaryMapL2_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
338  meAmplSummaryMapL2_[1]->setAxisTitle("ix", 1);
339  meAmplSummaryMapL2_[1]->setAxisTitle("iy", 2);
340 
341  for (int i = 0; i < 18; i++) {
342  name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G01 " + LN.str();
343  mePnAmplMapG01L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
344  mePnAmplMapG01L2_[i]->setAxisTitle("channel", 1);
345  mePnAmplMapG01L2_[i]->setAxisTitle("amplitude", 2);
346  dqmStore_->tag(mePnAmplMapG01L2_[i], i+1);
347 
348  name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G01 " + LN.str();
349  mePnPedMapG01L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
350  mePnPedMapG01L2_[i]->setAxisTitle("channel", 1);
351  mePnPedMapG01L2_[i]->setAxisTitle("pedestal", 2);
352  dqmStore_->tag(mePnPedMapG01L2_[i], i+1);
353  }
354 
355  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain16");
356 
357  for (int i = 0; i < 18; i++) {
358  name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G16 " + LN.str();
359  mePnAmplMapG16L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
360  mePnAmplMapG16L2_[i]->setAxisTitle("channel", 1);
361  mePnAmplMapG16L2_[i]->setAxisTitle("amplitude", 2);
362  dqmStore_->tag(mePnAmplMapG16L2_[i], i+1);
363 
364  name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G16 " + LN.str();
365  mePnPedMapG16L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
366  mePnPedMapG16L2_[i]->setAxisTitle("channel", 1);
367  mePnPedMapG16L2_[i]->setAxisTitle("pedestal", 2);
368  dqmStore_->tag(mePnPedMapG16L2_[i], i+1);
369  }
370 
371  }
372 
373  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
374 
375  LaserN.str("");
376  LaserN << "Laser" << 3;
377  LN.str("");
378  LN << "L" << 3;
379 
380  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str());
381  for (int i = 0; i < 18; i++) {
382  name = "EELT shape " + Numbers::sEE(i+1) + " " + LN.str();
383  meShapeMapL3_[i] = dqmStore_->bookProfile2D(name, name, 850, 0., 850., 10, 0., 10., 4096, 0., 4096., "s");
384  meShapeMapL3_[i]->setAxisTitle("channel", 1);
385  meShapeMapL3_[i]->setAxisTitle("sample", 2);
386  meShapeMapL3_[i]->setAxisTitle("amplitude", 3);
387  dqmStore_->tag(meShapeMapL3_[i], i+1);
388 
389  name = "EELT amplitude " + Numbers::sEE(i+1) + " " + LN.str();
390  meAmplMapL3_[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.*12., "s");
391  meAmplMapL3_[i]->setAxisTitle("ix", 1);
392  if ( i+1 >= 1 && i+1 <= 9 ) meAmplMapL3_[i]->setAxisTitle("101-ix", 1);
393  meAmplMapL3_[i]->setAxisTitle("iy", 2);
394  dqmStore_->tag(meAmplMapL3_[i], i+1);
395 
396  name = "EELT timing " + Numbers::sEE(i+1) + " " + LN.str();
397  meTimeMapL3_[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., 250, 0., 10., "s");
398  meTimeMapL3_[i]->setAxisTitle("ix", 1);
399  if ( i+1 >= 1 && i+1 <= 9 ) meTimeMapL3_[i]->setAxisTitle("101-ix", 1);
400  meTimeMapL3_[i]->setAxisTitle("iy", 2);
401  dqmStore_->tag(meTimeMapL3_[i], i+1);
402 
403  name = "EELT amplitude over PN " + Numbers::sEE(i+1) + " " + LN.str();
404  meAmplPNMapL3_[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.*12., "s");
405  meAmplPNMapL3_[i]->setAxisTitle("ix", 1);
406  if ( i+1 >= 1 && i+1 <= 9 ) meAmplPNMapL3_[i]->setAxisTitle("101-ix", 1);
407  meAmplPNMapL3_[i]->setAxisTitle("iy", 2);
408  dqmStore_->tag(meAmplPNMapL3_[i], i+1);
409  }
410 
411  name = "EELT amplitude map " + LN.str() + " EE -";
412  meAmplSummaryMapL3_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
413  meAmplSummaryMapL3_[0]->setAxisTitle("ix", 1);
414  meAmplSummaryMapL3_[0]->setAxisTitle("iy", 2);
415 
416  name = "EELT amplitude map " + LN.str() + " EE +";
417  meAmplSummaryMapL3_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
418  meAmplSummaryMapL3_[1]->setAxisTitle("ix", 1);
419  meAmplSummaryMapL3_[1]->setAxisTitle("iy", 2);
420 
421  for (int i = 0; i < 18; i++) {
422  name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G01 " + LN.str();
423  mePnAmplMapG01L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
424  mePnAmplMapG01L3_[i]->setAxisTitle("channel", 1);
425  mePnAmplMapG01L3_[i]->setAxisTitle("amplitude", 2);
426  dqmStore_->tag(mePnAmplMapG01L3_[i], i+1);
427 
428  name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G01 " + LN.str();
429  mePnPedMapG01L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
430  mePnPedMapG01L3_[i]->setAxisTitle("channel", 1);
431  mePnPedMapG01L3_[i]->setAxisTitle("pedestal", 2);
432  dqmStore_->tag(mePnPedMapG01L3_[i], i+1);
433  }
434 
435  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain16");
436 
437  for (int i = 0; i < 18; i++) {
438  name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G16 " + LN.str();
439  mePnAmplMapG16L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
440  mePnAmplMapG16L3_[i]->setAxisTitle("channel", 1);
441  mePnAmplMapG16L3_[i]->setAxisTitle("amplitude", 2);
442  dqmStore_->tag(mePnAmplMapG16L3_[i], i+1);
443 
444  name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G16 " + LN.str();
445  mePnPedMapG16L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
446  mePnPedMapG16L3_[i]->setAxisTitle("channel", 1);
447  mePnPedMapG16L3_[i]->setAxisTitle("pedestal", 2);
448  dqmStore_->tag(mePnPedMapG16L3_[i], i+1);
449  }
450 
451  }
452 
453  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
454 
455  LaserN.str("");
456  LaserN << "Laser" << 4;
457  LN.str("");
458  LN << "L" << 4;
459 
460  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str());
461  for (int i = 0; i < 18; i++) {
462  name = "EELT shape " + Numbers::sEE(i+1) + " " + LN.str();
463  meShapeMapL4_[i] = dqmStore_->bookProfile2D(name, name, 850, 0., 850., 10, 0., 10., 4096, 0., 4096., "s");
464  meShapeMapL4_[i]->setAxisTitle("channel", 1);
465  meShapeMapL4_[i]->setAxisTitle("sample", 2);
466  meShapeMapL4_[i]->setAxisTitle("amplitude", 3);
467  dqmStore_->tag(meShapeMapL4_[i], i+1);
468 
469  name = "EELT amplitude " + Numbers::sEE(i+1) + " " + LN.str();
470  meAmplMapL4_[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.*12., "s");
471  meAmplMapL4_[i]->setAxisTitle("ix", 1);
472  if ( i+1 >= 1 && i+1 <= 9 ) meAmplMapL4_[i]->setAxisTitle("101-ix", 1);
473  meAmplMapL4_[i]->setAxisTitle("iy", 2);
474  dqmStore_->tag(meAmplMapL4_[i], i+1);
475 
476  name = "EELT timing " + Numbers::sEE(i+1) + " " + LN.str();
477  meTimeMapL4_[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., 250, 0., 10., "s");
478  meTimeMapL4_[i]->setAxisTitle("ix", 1);
479  if ( i+1 >= 1 && i+1 <= 9 ) meTimeMapL4_[i]->setAxisTitle("101-ix", 1);
480  meTimeMapL4_[i]->setAxisTitle("iy", 2);
481  dqmStore_->tag(meTimeMapL4_[i], i+1);
482 
483  name = "EELT amplitude over PN " + Numbers::sEE(i+1) + " " + LN.str();
484  meAmplPNMapL4_[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.*12., "s");
485  meAmplPNMapL4_[i]->setAxisTitle("ix", 1);
486  if ( i+1 >= 1 && i+1 <= 9 ) meAmplPNMapL4_[i]->setAxisTitle("101-ix", 1);
487  meAmplPNMapL4_[i]->setAxisTitle("iy", 2);
488  dqmStore_->tag(meAmplPNMapL4_[i], i+1);
489  }
490 
491  name = "EELT amplitude map " + LN.str() + " EE -";
492  meAmplSummaryMapL4_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
493  meAmplSummaryMapL4_[0]->setAxisTitle("ix", 1);
494  meAmplSummaryMapL4_[0]->setAxisTitle("iy", 2);
495 
496  name = "EELT amplitude map " + LN.str() + " EE +";
497  meAmplSummaryMapL4_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
498  meAmplSummaryMapL4_[1]->setAxisTitle("ix", 1);
499  meAmplSummaryMapL4_[1]->setAxisTitle("iy", 2);
500 
501  for (int i = 0; i < 18; i++) {
502  name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G01 " + LN.str();
503  mePnAmplMapG01L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
504  mePnAmplMapG01L4_[i]->setAxisTitle("channel", 1);
505  mePnAmplMapG01L4_[i]->setAxisTitle("amplitude", 2);
506  dqmStore_->tag(mePnAmplMapG01L4_[i], i+1);
507 
508  name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G01 " + LN.str();
509  mePnPedMapG01L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
510  mePnPedMapG01L4_[i]->setAxisTitle("channel", 1);
511  mePnPedMapG01L4_[i]->setAxisTitle("pedestal", 2);
512  dqmStore_->tag(mePnPedMapG01L4_[i], i+1);
513  }
514 
515  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain16");
516 
517  for (int i = 0; i < 18; i++) {
518  name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G16 " + LN.str();
519  mePnAmplMapG16L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
520  mePnAmplMapG16L4_[i]->setAxisTitle("channel", 1);
521  mePnAmplMapG16L4_[i]->setAxisTitle("amplitude", 2);
522  dqmStore_->tag(mePnAmplMapG16L4_[i], i+1);
523 
524  name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G16 " + LN.str();
525  mePnPedMapG16L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
526  mePnPedMapG16L4_[i]->setAxisTitle("channel", 1);
527  mePnPedMapG16L4_[i]->setAxisTitle("pedestal", 2);
528  dqmStore_->tag(mePnPedMapG16L4_[i], i+1);
529  }
530 
531  }
532 
533  }
534 
535 }
536 
537 void EELaserTask::cleanup(void){
538 
539  if ( ! init_ ) return;
540 
541  if ( dqmStore_ ) {
542  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask");
543 
544  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
545 
546  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser1");
547  for (int i = 0; i < 18; i++) {
548  if ( meShapeMapL1_[i] ) dqmStore_->removeElement( meShapeMapL1_[i]->getName() );
549  meShapeMapL1_[i] = 0;
550  if ( meAmplMapL1_[i] ) dqmStore_->removeElement( meAmplMapL1_[i]->getName() );
551  meAmplMapL1_[i] = 0;
552  if ( meTimeMapL1_[i] ) dqmStore_->removeElement( meTimeMapL1_[i]->getName() );
553  meTimeMapL1_[i] = 0;
554  if ( meAmplPNMapL1_[i] ) dqmStore_->removeElement( meAmplPNMapL1_[i]->getName() );
555  meAmplPNMapL1_[i] = 0;
556  }
557 
558  for(int i=0; i<2; i++)
559  if( meAmplSummaryMapL1_[i] ) dqmStore_->removeElement( meAmplSummaryMapL1_[i]->getName() );
560  }
561 
562  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
563 
564  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser2");
565  for (int i = 0; i < 18; i++) {
566  if ( meShapeMapL2_[i] ) dqmStore_->removeElement( meShapeMapL2_[i]->getName() );
567  meShapeMapL2_[i] = 0;
568  if ( meAmplMapL2_[i] ) dqmStore_->removeElement( meAmplMapL2_[i]->getName() );
569  meAmplMapL2_[i] = 0;
570  if ( meTimeMapL2_[i] ) dqmStore_->removeElement( meTimeMapL2_[i]->getName() );
571  meTimeMapL2_[i] = 0;
572  if ( meAmplPNMapL2_[i] ) dqmStore_->removeElement( meAmplPNMapL2_[i]->getName() );
573  meAmplPNMapL2_[i] = 0;
574  }
575 
576  for(int i=0; i<2; i++)
577  if( meAmplSummaryMapL2_[i] ) dqmStore_->removeElement( meAmplSummaryMapL2_[i]->getName() );
578 
579  }
580 
581  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
582 
583  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser3");
584  for (int i = 0; i < 18; i++) {
585  if ( meShapeMapL3_[i] ) dqmStore_->removeElement( meShapeMapL3_[i]->getName() );
586  meShapeMapL3_[i] = 0;
587  if ( meAmplMapL3_[i] ) dqmStore_->removeElement( meAmplMapL3_[i]->getName() );
588  meAmplMapL3_[i] = 0;
589  if ( meTimeMapL3_[i] ) dqmStore_->removeElement( meTimeMapL3_[i]->getName() );
590  meTimeMapL3_[i] = 0;
591  if ( meAmplPNMapL3_[i] ) dqmStore_->removeElement( meAmplPNMapL3_[i]->getName() );
592  meAmplPNMapL3_[i] = 0;
593  }
594 
595  for(int i=0; i<2; i++)
596  if( meAmplSummaryMapL3_[i] ) dqmStore_->removeElement( meAmplSummaryMapL3_[i]->getName() );
597 
598  }
599 
600  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
601 
602  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser4");
603  for (int i = 0; i < 18; i++) {
604  if ( meShapeMapL4_[i] ) dqmStore_->removeElement( meShapeMapL4_[i]->getName() );
605  meShapeMapL4_[i] = 0;
606  if ( meAmplMapL4_[i] ) dqmStore_->removeElement( meAmplMapL4_[i]->getName() );
607  meAmplMapL4_[i] = 0;
608  if ( meTimeMapL4_[i] ) dqmStore_->removeElement( meTimeMapL4_[i]->getName() );
609  meTimeMapL4_[i] = 0;
610  if ( meAmplPNMapL4_[i] ) dqmStore_->removeElement( meAmplPNMapL4_[i]->getName() );
611  meAmplPNMapL4_[i] = 0;
612  }
613 
614  for(int i=0; i<2; i++)
615  if( meAmplSummaryMapL4_[i] ) dqmStore_->removeElement( meAmplSummaryMapL4_[i]->getName() );
616 
617  }
618 
619  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
620 
621  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser1/PN");
622 
623  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser1/PN/Gain01");
624  for (int i = 0; i < 18; i++) {
625  if ( mePnAmplMapG01L1_[i] ) dqmStore_->removeElement( mePnAmplMapG01L1_[i]->getName() );
626  mePnAmplMapG01L1_[i] = 0;
627  if ( mePnPedMapG01L1_[i] ) dqmStore_->removeElement( mePnPedMapG01L1_[i]->getName() );
628  mePnPedMapG01L1_[i] = 0;
629  }
630 
631  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser1/PN/Gain16");
632  for (int i = 0; i < 18; i++) {
633  if ( mePnAmplMapG16L1_[i] ) dqmStore_->removeElement( mePnAmplMapG16L1_[i]->getName() );
634  mePnAmplMapG16L1_[i] = 0;
635  if ( mePnPedMapG16L1_[i] ) dqmStore_->removeElement( mePnPedMapG16L1_[i]->getName() );
636  mePnPedMapG16L1_[i] = 0;
637  }
638 
639  }
640 
641  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
642 
643  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser2/PN");
644 
645  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser2/PN/Gain01");
646  for (int i = 0; i < 18; i++) {
647  if ( mePnAmplMapG01L2_[i] ) dqmStore_->removeElement( mePnAmplMapG01L2_[i]->getName() );
648  mePnAmplMapG01L2_[i] = 0;
649  if ( mePnPedMapG01L2_[i] ) dqmStore_->removeElement( mePnPedMapG01L2_[i]->getName() );
650  mePnPedMapG01L2_[i] = 0;
651  }
652 
653  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser2/PN/Gain16");
654  for (int i = 0; i < 18; i++) {
655  if ( mePnAmplMapG16L2_[i] ) dqmStore_->removeElement( mePnAmplMapG16L2_[i]->getName() );
656  mePnAmplMapG16L2_[i] = 0;
657  if ( mePnPedMapG16L2_[i] ) dqmStore_->removeElement( mePnPedMapG16L2_[i]->getName() );
658  mePnPedMapG16L2_[i] = 0;
659  }
660 
661  }
662 
663  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
664 
665  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser3/PN");
666 
667  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser3/PN/Gain01");
668  for (int i = 0; i < 18; i++) {
669  if ( mePnAmplMapG01L3_[i] ) dqmStore_->removeElement( mePnAmplMapG01L3_[i]->getName() );
670  mePnAmplMapG01L3_[i] = 0;
671  if ( mePnPedMapG01L3_[i] ) dqmStore_->removeElement( mePnPedMapG01L3_[i]->getName() );
672  mePnPedMapG01L3_[i] = 0;
673  }
674 
675  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser3/PN/Gain16");
676  for (int i = 0; i < 18; i++) {
677  if ( mePnAmplMapG16L3_[i] ) dqmStore_->removeElement( mePnAmplMapG16L3_[i]->getName() );
678  mePnAmplMapG16L3_[i] = 0;
679  if ( mePnPedMapG16L3_[i] ) dqmStore_->removeElement( mePnPedMapG16L3_[i]->getName() );
680  mePnPedMapG16L3_[i] = 0;
681  }
682 
683  }
684 
685  if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
686 
687  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser4/PN");
688 
689  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser4/PN/Gain01");
690  for (int i = 0; i < 18; i++) {
691  if ( mePnAmplMapG01L4_[i] ) dqmStore_->removeElement( mePnAmplMapG01L4_[i]->getName() );
692  mePnAmplMapG01L4_[i] = 0;
693  if ( mePnPedMapG01L4_[i] ) dqmStore_->removeElement( mePnPedMapG01L4_[i]->getName() );
694  mePnPedMapG01L4_[i] = 0;
695  }
696 
697  dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser4/PN/Gain16");
698  for (int i = 0; i < 18; i++) {
699  if ( mePnAmplMapG16L4_[i] ) dqmStore_->removeElement( mePnAmplMapG16L4_[i]->getName() );
700  mePnAmplMapG16L4_[i] = 0;
701  if ( mePnPedMapG16L4_[i] ) dqmStore_->removeElement( mePnPedMapG16L4_[i]->getName() );
702  mePnPedMapG16L4_[i] = 0;
703  }
704 
705  }
706 
707  }
708 
709  init_ = false;
710 
711 }
712 
713 void EELaserTask::endJob(void){
714 
715  edm::LogInfo("EELaserTask") << "analyzed " << ievt_ << " events";
716 
717  if ( enableCleanup_ ) this->cleanup();
718 
719 }
720 
721 void EELaserTask::analyze(const edm::Event& e, const edm::EventSetup& c){
722 
723  bool enable = false;
724  int runType[18];
725  for (int i=0; i<18; i++) runType[i] = -1;
726  unsigned rtHalf[18];
727  for (int i=0; i<18; i++) rtHalf[i] = -1;
728  int waveLength[18];
729  for (int i=0; i<18; i++) waveLength[i] = -1;
730 
731  edm::Handle<EcalRawDataCollection> dcchs;
732 
733  if ( e.getByToken(EcalRawDataCollection_, dcchs) ) {
734 
735  for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
736 
737  if ( Numbers::subDet( *dcchItr ) != EcalEndcap ) continue;
738 
739  int ism = Numbers::iSM( *dcchItr, EcalEndcap );
740 
741  runType[ism-1] = dcchItr->getRunType();
742  rtHalf[ism-1] = dcchItr->getRtHalf();
743  waveLength[ism-1] = dcchItr->getEventSettings().wavelength;
744 
745  if ( dcchItr->getRunType() == EcalDCCHeaderBlock::LASER_STD ||
746  dcchItr->getRunType() == EcalDCCHeaderBlock::LASER_GAP ) enable = true;
747 
748  }
749 
750  } else {
751 
752  edm::LogWarning("EELaserTask") << "EcalRawDataCollection not available";
753 
754  }
755 
756  if ( ! enable ) return;
757 
758  if ( ! init_ ) this->setup();
759 
760  ievt_++;
761 
762  bool numPN[80];
763  float adcPN[80];
764  for ( int i = 0; i < 80; i++ ) {
765  numPN[i] = false;
766  adcPN[i] = 0.;
767  }
768 
769  std::vector<int> PNs;
770  PNs.reserve(12);
771 
772  edm::Handle<EEDigiCollection> digis;
773 
774  if ( e.getByToken(EEDigiCollection_, digis) ) {
775 
776  int maxpos[10];
777  for(int i(0); i < 10; i++)
778  maxpos[i] = 0;
779  int nReadouts(0);
780 
781  for ( EEDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr ) {
782 
783  EEDetId id = digiItr->id();
784 
785  int ism = Numbers::iSM( id );
786 
787  if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
788  runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
789 
790  if ( rtHalf[ism-1] != Numbers::RtHalf(id) ) continue;
791 
792  nReadouts++;
793 
794  EEDataFrame dataframe = (*digiItr);
795 
796  int iMax(-1);
797  float max(0.);
798  float min(4096.);
799  for (int i = 0; i < 10; i++) {
800  int adc = dataframe.sample(i).adc();
801  if(adc > max){
802  max = adc;
803  iMax = i;
804  }
805  if(adc < min)
806  min = adc;
807  }
808  if(iMax >= 0 && max - min > 20.)
809  maxpos[iMax] += 1;
810 
811  }
812 
813  int threshold(nReadouts / 2);
814  enable = false;
815  for(int i(0); i < 10; i++){
816  if(maxpos[i] > threshold){
817  enable = true;
818  break;
819  }
820  }
821 
822  if(!enable) return;
823 
824  int need = digis->size();
825  LogDebug("EELaserTask") << "event " << ievt_ << " digi collection size " << need;
826 
827  for ( EEDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr ) {
828 
829  EEDetId id = digiItr->id();
830 
831  int ix = id.ix();
832  int iy = id.iy();
833 
834  int ism = Numbers::iSM( id );
835 
836  if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
837  runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
838 
839  if ( rtHalf[ism-1] != Numbers::RtHalf(id) ) continue;
840 
841  int ic = Numbers::icEE(ism, ix, iy);
842 
843  EEDataFrame dataframe = (*digiItr);
844 
845  for (int i = 0; i < 10; i++) {
846 
847  int adc = dataframe.sample(i).adc();
848 
849  MonitorElement* meShapeMap = 0;
850 
851  if ( rtHalf[ism-1] == 0 || rtHalf[ism-1] == 1 ) {
852 
853  if ( waveLength[ism-1] == 0 ) meShapeMap = meShapeMapL1_[ism-1];
854  if ( waveLength[ism-1] == 1 ) meShapeMap = meShapeMapL2_[ism-1];
855  if ( waveLength[ism-1] == 2 ) meShapeMap = meShapeMapL3_[ism-1];
856  if ( waveLength[ism-1] == 3 ) meShapeMap = meShapeMapL4_[ism-1];
857 
858  } else {
859 
860  edm::LogWarning("EELaserTask") << " RtHalf = " << rtHalf[ism-1];
861 
862  }
863 
864  float xval = float(adc);
865 
866  if ( meShapeMap ) meShapeMap->Fill(ic - 0.5, i + 0.5, xval);
867 
868  }
869 
870  NumbersPn::getPNs( ism, ix, iy, PNs );
871 
872  for (unsigned int i=0; i<PNs.size(); i++) {
873  int ipn = PNs[i];
874  if ( ipn >= 0 && ipn < 80 ) numPN[ipn] = true;
875  }
876 
877  }
878 
879  } else {
880 
881  edm::LogWarning("EELaserTask") << "EEDigiCollection not available";
882 
883  }
884 
885  edm::Handle<EcalPnDiodeDigiCollection> pns;
886 
887  if ( e.getByToken(EcalPnDiodeDigiCollection_, pns) ) {
888 
889  int nep = pns->size();
890  LogDebug("EELaserTask") << "event " << ievt_ << " pns collection size " << nep;
891 
892  for ( EcalPnDiodeDigiCollection::const_iterator pnItr = pns->begin(); pnItr != pns->end(); ++pnItr ) {
893 
894  if ( Numbers::subDet( pnItr->id() ) != EcalEndcap ) continue;
895 
896  int ism = Numbers::iSM( pnItr->id() );
897 
898  int num = pnItr->id().iPnId();
899 
900  if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
901  runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
902 
903  int ipn = NumbersPn::ipnEE( ism, num );
904 
905  if ( ipn >= 0 && ipn < 80 && numPN[ipn] == false ) continue;
906 
907  float xvalped = 0.;
908 
909  for (int i = 0; i < 4; i++) {
910 
911  int adc = pnItr->sample(i).adc();
912 
913  MonitorElement* mePNPed = 0;
914 
915  if ( pnItr->sample(i).gainId() == 0 ) {
916  if ( waveLength[ism-1] == 0 ) mePNPed = mePnPedMapG01L1_[ism-1];
917  if ( waveLength[ism-1] == 1 ) mePNPed = mePnPedMapG01L2_[ism-1];
918  if ( waveLength[ism-1] == 2 ) mePNPed = mePnPedMapG01L3_[ism-1];
919  if ( waveLength[ism-1] == 3 ) mePNPed = mePnPedMapG01L4_[ism-1];
920  }
921  if ( pnItr->sample(i).gainId() == 1 ) {
922  if ( waveLength[ism-1] == 0 ) mePNPed = mePnPedMapG16L1_[ism-1];
923  if ( waveLength[ism-1] == 1 ) mePNPed = mePnPedMapG16L2_[ism-1];
924  if ( waveLength[ism-1] == 2 ) mePNPed = mePnPedMapG16L3_[ism-1];
925  if ( waveLength[ism-1] == 3 ) mePNPed = mePnPedMapG16L4_[ism-1];
926  }
927 
928  float xval = float(adc);
929 
930  if ( mePNPed ) mePNPed->Fill(num - 0.5, xval);
931 
932  xvalped = xvalped + xval;
933 
934  }
935 
936  xvalped = xvalped / 4;
937 
938  float xvalmax = 0.;
939 
940  MonitorElement* mePN = 0;
941 
942  for (int i = 0; i < 50; i++) {
943 
944  int adc = pnItr->sample(i).adc();
945 
946  float xval = float(adc);
947 
948  if ( xval >= xvalmax ) xvalmax = xval;
949 
950  }
951 
952  xvalmax = xvalmax - xvalped;
953 
954  if ( pnItr->sample(0).gainId() == 0 ) {
955  if ( waveLength[ism-1] == 0 ) mePN = mePnAmplMapG01L1_[ism-1];
956  if ( waveLength[ism-1] == 1 ) mePN = mePnAmplMapG01L2_[ism-1];
957  if ( waveLength[ism-1] == 2 ) mePN = mePnAmplMapG01L3_[ism-1];
958  if ( waveLength[ism-1] == 3 ) mePN = mePnAmplMapG01L4_[ism-1];
959  }
960  if ( pnItr->sample(0).gainId() == 1 ) {
961  if ( waveLength[ism-1] == 0 ) mePN = mePnAmplMapG16L1_[ism-1];
962  if ( waveLength[ism-1] == 1 ) mePN = mePnAmplMapG16L2_[ism-1];
963  if ( waveLength[ism-1] == 2 ) mePN = mePnAmplMapG16L3_[ism-1];
964  if ( waveLength[ism-1] == 3 ) mePN = mePnAmplMapG16L4_[ism-1];
965  }
966 
967  if ( mePN ) mePN->Fill(num - 0.5, xvalmax);
968 
969  if ( ipn >= 0 && ipn < 80 ) adcPN[ipn] = xvalmax;
970 
971  }
972 
973  } else {
974 
975  edm::LogWarning("EELaserTask") << "EcalPnDiodeDigiCollection not available";
976 
977  }
978 
979  edm::Handle<EcalUncalibratedRecHitCollection> hits;
980 
981  if ( e.getByToken(EcalUncalibratedRecHitCollection_, hits) ) {
982 
983  int neh = hits->size();
984  LogDebug("EELaserTask") << "event " << ievt_ << " hits collection size " << neh;
985 
986  for ( EcalUncalibratedRecHitCollection::const_iterator hitItr = hits->begin(); hitItr != hits->end(); ++hitItr ) {
987 
988  EEDetId id = hitItr->id();
989 
990  int ix = id.ix();
991  int iy = id.iy();
992  int iz;
993 
994  int ism = Numbers::iSM( id );
995 
996  if ( ism >= 1 && ism <= 9 ){
997  ix = 101 - ix;
998  iz = 0;
999  }else{
1000  iz = 1;
1001  }
1002 
1003  float xix = ix - 0.5;
1004  float xiy = iy - 0.5;
1005 
1006  if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
1007  runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
1008 
1009  if ( rtHalf[ism-1] != Numbers::RtHalf(id) ) continue;
1010 
1011  MonitorElement* meAmplMap = 0;
1012  MonitorElement* meTimeMap = 0;
1013  MonitorElement* meAmplPNMap = 0;
1014  MonitorElement* meAmplSummaryMap = 0;
1015 
1016  if ( rtHalf[ism-1] == 0 || rtHalf[ism-1] == 1 ) {
1017 
1018  if ( waveLength[ism-1] == 0 ) {
1019  meAmplMap = meAmplMapL1_[ism-1];
1020  meTimeMap = meTimeMapL1_[ism-1];
1021  meAmplPNMap = meAmplPNMapL1_[ism-1];
1022  meAmplSummaryMap = meAmplSummaryMapL1_[iz];
1023  }
1024  if ( waveLength[ism-1] == 1 ) {
1025  meAmplMap = meAmplMapL2_[ism-1];
1026  meTimeMap = meTimeMapL2_[ism-1];
1027  meAmplPNMap = meAmplPNMapL2_[ism-1];
1028  meAmplSummaryMap = meAmplSummaryMapL2_[iz];
1029  }
1030  if ( waveLength[ism-1] == 2 ) {
1031  meAmplMap = meAmplMapL3_[ism-1];
1032  meTimeMap = meTimeMapL3_[ism-1];
1033  meAmplPNMap = meAmplPNMapL3_[ism-1];
1034  meAmplSummaryMap = meAmplSummaryMapL3_[iz];
1035  }
1036  if ( waveLength[ism-1] == 3 ) {
1037  meAmplMap = meAmplMapL4_[ism-1];
1038  meTimeMap = meTimeMapL4_[ism-1];
1039  meAmplPNMap = meAmplPNMapL4_[ism-1];
1040  meAmplSummaryMap = meAmplSummaryMapL4_[iz];
1041  }
1042 
1043  } else {
1044 
1045  edm::LogWarning("EELaserTask") << " RtHalf = " << rtHalf[ism-1];
1046 
1047  }
1048 
1049  float xval = hitItr->amplitude();
1050  if ( xval <= 0. ) xval = 0.0;
1051  float yval = hitItr->jitter() + 5.0;
1052  if ( yval <= 0. ) yval = 0.0;
1053  float zval = hitItr->pedestal();
1054  if ( zval <= 0. ) zval = 0.0;
1055 
1056  if ( meAmplMap ) meAmplMap->Fill(xix, xiy, xval);
1057 
1058  if ( xval > 16. ) {
1059  if ( meTimeMap ) meTimeMap->Fill(xix, xiy, yval);
1060  }
1061 
1062  float wval = 0.;
1063 
1064  NumbersPn::getPNs( ism, ix, iy, PNs );
1065 
1066  if ( PNs.size() > 0 ) {
1067  int ipn = PNs[0];
1068  if ( ipn >= 0 && ipn < 80 ) {
1069  if ( adcPN[ipn] != 0. ) wval = xval / adcPN[ipn];
1070  }
1071  }
1072 
1073  if ( meAmplPNMap ) meAmplPNMap->Fill(xix, xiy, wval);
1074 
1075  if( meAmplSummaryMap ) meAmplSummaryMap->Fill(xix, xiy, xval);
1076 
1077  }
1078 
1079  } else {
1080 
1081  edm::LogWarning("EELaserTask") << "EcalUncalibratedRecHitCollection not available";
1082 
1083  }
1084 
1085 }
1086 
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
ecalMGPA::adc
int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:11
EELaserTask::mePnAmplMapG01L2_
MonitorElement * mePnAmplMapG01L2_[18]
Definition: EELaserTask.h:89
NumbersPn.h
Some &quot;id&quot; conversions.
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
Numbers::sEE
static std::string sEE(const unsigned ism)
Definition: Numbers.cc:223
EELaserTask::meTimeMapL4_
MonitorElement * meTimeMapL4_[18]
Definition: EELaserTask.h:105
i
int i
Definition: DBlmapReader.cc:9
NumbersPn::ipnEE
static int ipnEE(const int ism, const int ipnid)
Definition: NumbersPn.cc:18
EELaserTask::mePnAmplMapG16L2_
MonitorElement * mePnAmplMapG16L2_[18]
Definition: EELaserTask.h:91
EEDetId::ix
int ix() const
Definition: EEDetId.h:76
EELaserTask::EEDigiCollection_
edm::EDGetTokenT< EEDigiCollection > EEDigiCollection_
Definition: EELaserTask.h:71
edm::DataFrameContainer::const_iterator
boost::transform_iterator< IterHelp, boost::counting_iterator< int > > const_iterator
Definition: DataFrameContainer.h:51
EELaserTask::meTimeMapL3_
MonitorElement * meTimeMapL3_[18]
Definition: EELaserTask.h:96
Numbers::icEE
static unsigned icEE(const unsigned ism, const unsigned ix, const unsigned iy)
Definition: Numbers.cc:736
EELaserTask::mePnAmplMapG16L1_
MonitorElement * mePnAmplMapG16L1_[18]
Definition: EELaserTask.h:82
DQMStore::rmdir
void rmdir(const std::string &fullpath)
Definition: DQMStore.cc:2730
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434
EELaserTask::prefixME_
std::string prefixME_
Definition: EELaserTask.h:64
Numbers.h
Some &quot;id&quot; conversions.
EELaserTask::mePnPedMapG01L4_
MonitorElement * mePnPedMapG01L4_[18]
Definition: EELaserTask.h:108
EELaserTask::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Analyze.
Definition: EELaserTask.cc:721
EELaserTask::mePnPedMapG01L2_
MonitorElement * mePnPedMapG01L2_[18]
Definition: EELaserTask.h:90
edm::SortedCollection< EcalDCCHeaderBlock >::const_iterator
std::vector< EcalDCCHeaderBlock >::const_iterator const_iterator
Definition: SortedCollection.h:80
EELaserTask::meAmplMapL2_
MonitorElement * meAmplMapL2_[18]
Definition: EELaserTask.h:86
EcalDataFrame::sample
EcalMGPASample sample(int i) const
Definition: EcalDataFrame.h:29
EELaserTask::mePnAmplMapG01L1_
MonitorElement * mePnAmplMapG01L1_[18]
Definition: EELaserTask.h:80
Numbers::ix0EE
static int ix0EE(const unsigned ism)
Definition: Numbers.cc:770
EELaserTask::mePnPedMapG01L3_
MonitorElement * mePnPedMapG01L3_[18]
Definition: EELaserTask.h:99
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
Numbers::iy0EE
static int iy0EE(const unsigned ism)
Definition: Numbers.cc:809
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:161
DQMStore::tag
void tag(MonitorElement *me, unsigned int myTag)
Definition: DQMStore.cc:1509
EELaserTask::mePnPedMapG16L2_
MonitorElement * mePnPedMapG16L2_[18]
Definition: EELaserTask.h:92
EELaserTask::mePnAmplMapG01L4_
MonitorElement * mePnAmplMapG01L4_[18]
Definition: EELaserTask.h:107
EELaserTask::~EELaserTask
virtual ~EELaserTask()
Destructor.
Definition: EELaserTask.cc:97
EELaserTask::meTimeMapL2_
MonitorElement * meTimeMapL2_[18]
Definition: EELaserTask.h:87
EELaserTask::mePnPedMapG16L4_
MonitorElement * mePnPedMapG16L4_[18]
Definition: EELaserTask.h:110
EELaserTask::meShapeMapL3_
MonitorElement * meShapeMapL3_[18]
Definition: EELaserTask.h:94
max
const T & max(const T &a, const T &b)
Definition: MaterialBudgetTrackerHistos.cc:4
EELaserTask::meAmplMapL4_
MonitorElement * meAmplMapL4_[18]
Definition: EELaserTask.h:104
EELaserTask::mePnPedMapG16L1_
MonitorElement * mePnPedMapG16L1_[18]
Definition: EELaserTask.h:83
EELaserTask::meTimeMapL1_
MonitorElement * meTimeMapL1_[18]
Definition: EELaserTask.h:78
DQMStore::removeElement
void removeElement(const std::string &name)
Definition: DQMStore.cc:2772
EELaserTask::EcalPnDiodeDigiCollection_
edm::EDGetTokenT< EcalPnDiodeDigiCollection > EcalPnDiodeDigiCollection_
Definition: EELaserTask.h:72
EELaserTask::meAmplMapL3_
MonitorElement * meAmplMapL3_[18]
Definition: EELaserTask.h:95
EELaserTask::meAmplSummaryMapL1_
MonitorElement * meAmplSummaryMapL1_[2]
Definition: EELaserTask.h:112
dtDQMClient_cfg.threshold
tuple threshold
Definition: dtDQMClient_cfg.py:16
EELaserTask::meAmplSummaryMapL4_
MonitorElement * meAmplSummaryMapL4_[2]
Definition: EELaserTask.h:115
DQMStore::bookProfile
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
Definition: DQMStore.cc:1186
EELaserTask::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
BeginRun.
Definition: EELaserTask.cc:112
EELaserTask::reset
void reset(void)
Reset.
Definition: EELaserTask.cc:124
EELaserTask::dqmStore_
DQMStore * dqmStore_
Definition: EELaserTask.h:62
EELaserTask::meAmplPNMapL3_
MonitorElement * meAmplPNMapL3_[18]
Definition: EELaserTask.h:97
EELaserTask::laserWavelengths_
std::vector< int > laserWavelengths_
Definition: EELaserTask.h:74
EELaserTask::mePnPedMapG01L1_
MonitorElement * mePnPedMapG01L1_[18]
Definition: EELaserTask.h:81
EELaserTask::EcalRawDataCollection_
edm::EDGetTokenT< EcalRawDataCollection > EcalRawDataCollection_
Definition: EELaserTask.h:70
EELaserTask::meAmplMapL1_
MonitorElement * meAmplMapL1_[18]
Definition: EELaserTask.h:77
EELaserTask::EELaserTask
EELaserTask(const edm::ParameterSet &ps)
Constructor.
Definition: EELaserTask.cc:28
environment_file_cff.runType
tuple runType
Definition: environment_file_cff.py:30
EELaserTask::EcalUncalibratedRecHitCollection_
edm::EDGetTokenT< EcalUncalibratedRecHitCollection > EcalUncalibratedRecHitCollection_
Definition: EELaserTask.h:73
Numbers::initGeometry
static void initGeometry(const edm::EventSetup &setup, bool verbose=false)
Definition: Numbers.cc:47
EELaserTask::meAmplSummaryMapL3_
MonitorElement * meAmplSummaryMapL3_[2]
Definition: EELaserTask.h:114
EELaserTask::mePnAmplMapG16L3_
MonitorElement * mePnAmplMapG16L3_[18]
Definition: EELaserTask.h:100
EELaserTask::mePnAmplMapG16L4_
MonitorElement * mePnAmplMapG16L4_[18]
Definition: EELaserTask.h:109
EELaserTask::setup
void setup(void)
Setup.
Definition: EELaserTask.cc:199
EELaserTask::endJob
void endJob(void)
EndJob.
Definition: EELaserTask.cc:713
EELaserTask::mePnAmplMapG01L3_
MonitorElement * mePnAmplMapG01L3_[18]
Definition: EELaserTask.h:98
EELaserTask::cleanup
void cleanup(void)
Cleanup.
Definition: EELaserTask.cc:537
EELaserTask::mePnPedMapG16L3_
MonitorElement * mePnPedMapG16L3_[18]
Definition: EELaserTask.h:101
EELaserTask::meAmplPNMapL4_
MonitorElement * meAmplPNMapL4_[18]
Definition: EELaserTask.h:106
EELaserTask::meAmplPNMapL1_
MonitorElement * meAmplPNMapL1_[18]
Definition: EELaserTask.h:79
Numbers::iSM
static unsigned iSM(const unsigned ism, const EcalSubdetector subdet)
Definition: Numbers.cc:243
EELaserTask::meAmplSummaryMapL2_
MonitorElement * meAmplSummaryMapL2_[2]
Definition: EELaserTask.h:113
NumbersPn::getPNs
static void getPNs(const int ism, const int ix, const int iy, std::vector< int > &PNsInLM)
Definition: NumbersPn.cc:48
Numbers::RtHalf
static unsigned RtHalf(const EBDetId &id)
Definition: Numbers.cc:605
EELaserTask::mergeRuns_
bool mergeRuns_
Definition: EELaserTask.h:68
EELaserTask::endRun
void endRun(const edm::Run &r, const edm::EventSetup &c)
EndRun.
Definition: EELaserTask.cc:120
Numbers::subDet
static EcalSubdetector subDet(const EBDetId &id)
Definition: Numbers.cc:142
EELaserTask::enableCleanup_
bool enableCleanup_
Definition: EELaserTask.h:66
ecalpyutils::ism
int ism(int ieta, int iphi)
Definition: EcalPyUtils.cc:56
EELaserTask::meAmplPNMapL2_
MonitorElement * meAmplPNMapL2_[18]
Definition: EELaserTask.h:88
EELaserTask::meShapeMapL1_
MonitorElement * meShapeMapL1_[18]
Definition: EELaserTask.h:76
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
EELaserTask::beginJob
void beginJob(void)
BeginJob.
Definition: EELaserTask.cc:101
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:584
EELaserTask::meShapeMapL4_
MonitorElement * meShapeMapL4_[18]
Definition: EELaserTask.h:103
EELaserTask::meShapeMapL2_
MonitorElement * meShapeMapL2_[18]
Definition: EELaserTask.h:85
edm::Run
Definition: Run.h:41
EcalMGPASample::adc
int adc() const
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:34
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