CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EEPedestalTask.cc
Go to the documentation of this file.
1 /*
2  * \file EEPedestalTask.cc
3  *
4  * \author G. Della Ricca
5  *
6 */
7 
8 #include <iostream>
9 #include <sstream>
10 #include <iomanip>
11 #include <vector>
12 
13 #include "FWCore/ServiceRegistry/interface/Service.h"
14 #include "FWCore/MessageLogger/interface/MessageLogger.h"
15 
16 #include "DQMServices/Core/interface/MonitorElement.h"
17 
18 #include "DQMServices/Core/interface/DQMStore.h"
19 
20 #include "DataFormats/EcalDetId/interface/EEDetId.h"
21 #include "DataFormats/EcalDigi/interface/EEDataFrame.h"
22 
23 #include "DQM/EcalCommon/interface/Numbers.h"
24 
25 #include "DQM/EcalEndcapMonitorTasks/interface/EEPedestalTask.h"
26 
27 // #define COMMON_NOISE_ANALYSIS
28 
29 EEPedestalTask::EEPedestalTask(const edm::ParameterSet& ps){
30 
31  init_ = false;
32 
33  dqmStore_ = edm::Service<DQMStore>().operator->();
34 
35  prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
36 
37  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
38 
39  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
40 
41  EcalRawDataCollection_ = consumes<EcalRawDataCollection>(ps.getParameter<edm::InputTag>("EcalRawDataCollection"));
42  EEDigiCollection_ = consumes<EEDigiCollection>(ps.getParameter<edm::InputTag>("EEDigiCollection"));
43  EcalPnDiodeDigiCollection_ = consumes<EcalPnDiodeDigiCollection>(ps.getParameter<edm::InputTag>("EcalPnDiodeDigiCollection"));
44 
45  MGPAGains_.reserve(3);
46  for ( unsigned int i = 1; i <= 3; i++ ) MGPAGains_.push_back(i);
47  MGPAGains_ = ps.getUntrackedParameter<std::vector<int> >("MGPAGains", MGPAGains_);
48 
49  MGPAGainsPN_.reserve(2);
50  for ( unsigned int i = 1; i <= 3; i++ ) MGPAGainsPN_.push_back(i);
51  MGPAGainsPN_ = ps.getUntrackedParameter<std::vector<int> >("MGPAGainsPN", MGPAGainsPN_);
52 
53  for (int i = 0; i < 18; i++) {
54  mePedMapG01_[i] = 0;
55  mePedMapG06_[i] = 0;
56  mePedMapG12_[i] = 0;
57 #ifdef COMMON_NOISE_ANALYSIS
58  mePed3SumMapG01_[i] = 0;
59  mePed3SumMapG06_[i] = 0;
60  mePed3SumMapG12_[i] = 0;
61  mePed5SumMapG01_[i] = 0;
62  mePed5SumMapG06_[i] = 0;
63  mePed5SumMapG12_[i] = 0;
64 #endif
65  mePnPedMapG01_[i] = 0;
66  mePnPedMapG16_[i] = 0;
67  }
68 
69 }
70 
71 EEPedestalTask::~EEPedestalTask(){
72 
73 }
74 
75 void EEPedestalTask::beginJob(void){
76 
77  ievt_ = 0;
78 
79  if ( dqmStore_ ) {
80  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask");
81  dqmStore_->rmdir(prefixME_ + "/EEPedestalTask");
82  }
83 
84 }
85 
86 void EEPedestalTask::beginRun(const edm::Run& r, const edm::EventSetup& c) {
87 
88  Numbers::initGeometry(c, false);
89 
90  if ( ! mergeRuns_ ) this->reset();
91 
92 }
93 
94 void EEPedestalTask::endRun(const edm::Run& r, const edm::EventSetup& c) {
95 
96 }
97 
98 void EEPedestalTask::reset(void) {
99 
100  for (int i = 0; i < 18; i++) {
101  if (find(MGPAGains_.begin(), MGPAGains_.end(), 1) != MGPAGains_.end() ) {
102  if ( mePedMapG01_[i] ) mePedMapG01_[i]->Reset();
103  }
104  if (find(MGPAGains_.begin(), MGPAGains_.end(), 6) != MGPAGains_.end() ) {
105  if ( mePedMapG06_[i] ) mePedMapG06_[i]->Reset();
106  }
107  if (find(MGPAGains_.begin(), MGPAGains_.end(), 12) != MGPAGains_.end() ) {
108  if ( mePedMapG12_[i] ) mePedMapG12_[i]->Reset();
109  }
110 #ifdef COMMON_NOISE_ANALYSIS
111  if ( mePed3SumMapG01_[i] ) mePed3SumMapG01_[i]->Reset();
112  if ( mePed3SumMapG06_[i] ) mePed3SumMapG06_[i]->Reset();
113  if ( mePed3SumMapG12_[i] ) mePed3SumMapG12_[i]->Reset();
114  if ( mePed5SumMapG01_[i] ) mePed5SumMapG01_[i]->Reset();
115  if ( mePed5SumMapG06_[i] ) mePed5SumMapG06_[i]->Reset();
116  if ( mePed5SumMapG12_[i] ) mePed5SumMapG12_[i]->Reset();
117 #endif
118  if (find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), 1) != MGPAGainsPN_.end() ) {
119  if ( mePnPedMapG01_[i] ) mePnPedMapG01_[i]->Reset();
120  }
121  if (find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), 16) != MGPAGainsPN_.end() ) {
122  if ( mePnPedMapG16_[i] ) mePnPedMapG16_[i]->Reset();
123  }
124  }
125 
126 }
127 
128 void EEPedestalTask::setup(void){
129 
130  init_ = true;
131 
132  std::string name;
133  std::stringstream GainN, GN;
134 
135  if ( dqmStore_ ) {
136  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask");
137 
138  if (find(MGPAGains_.begin(), MGPAGains_.end(), 1) != MGPAGains_.end() ) {
139 
140  GainN.str("");
141  GainN << "Gain" << std::setw(2) << std::setfill('0') << 1;
142  GN.str("");
143  GN << "G" << std::setw(2) << std::setfill('0') << 1;
144 
145  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask/" + GainN.str());
146  for (int i = 0; i < 18; i++) {
147  name = "EEPT pedestal " + Numbers::sEE(i+1) + " " + GN.str();
148  mePedMapG01_[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");
149  mePedMapG01_[i]->setAxisTitle("ix", 1);
150  if ( i+1 >= 1 && i+1 <= 9 ) mePedMapG01_[i]->setAxisTitle("101-ix", 1);
151  mePedMapG01_[i]->setAxisTitle("iy", 2);
152  dqmStore_->tag(mePedMapG01_[i], i+1);
153 
154 #ifdef COMMON_NOISE_ANALYSIS
155  name = "EEPT pedestal 3sum " + Numbers::sEE(i+1) + " " + GN.str();
156  mePed3SumMapG01_[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");
157  mePed3SumMapG01_[i]->setAxisTitle("ix", 1);
158  if ( i+1 >= 1 && i+1 <= 9 ) mePed3SumMapG01_[i]->setAxisTitle("101-ix", 1);
159  mePed3SumMapG01_[i]->setAxisTitle("iy", 2);
160  dqmStore_->tag(mePed3SumMapG01_[i], i+1);
161 
162  name = "EEPT pedestal 5sum " + Numbers::sEE(i+1) + " " + GN.str();
163  mePed5SumMapG01_[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");
164  mePed5SumMapG01_[i]->setAxisTitle("ix", 1);
165  if ( i+1 >= 1 && i+1 <= 9 ) mePed5SumMapG01_[i]->setAxisTitle("101-ix", 1);
166  mePed5SumMapG01_[i]->setAxisTitle("iy", 2);
167  dqmStore_->tag(mePed5SumMapG01_[i], i+1);
168 #endif
169  }
170 
171  }
172 
173  if (find(MGPAGains_.begin(), MGPAGains_.end(), 6) != MGPAGains_.end() ) {
174 
175  GainN.str("");
176  GainN << "Gain" << std::setw(2) << std::setfill('0') << 6;
177  GN.str("");
178  GN << "G" << std::setw(2) << std::setfill('0') << 6;
179 
180  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask/" + GainN.str());
181  for (int i = 0; i < 18; i++) {
182  name = "EEPT pedestal " + Numbers::sEE(i+1) + " " + GN.str();
183  mePedMapG06_[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");
184  mePedMapG06_[i]->setAxisTitle("ix", 1);
185  if ( i+1 >= 1 && i+1 <= 9 ) mePedMapG06_[i]->setAxisTitle("101-ix", 1);
186  mePedMapG06_[i]->setAxisTitle("iy", 2);
187  dqmStore_->tag(mePedMapG06_[i], i+1);
188 
189 #ifdef COMMON_NOISE_ANALYSIS
190  name = "EEPT pedestal 3sum " + Numbers::sEE(i+1) + " " + GN.str();
191  mePed3SumMapG06_[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");
192  mePed3SumMapG06_[i]->setAxisTitle("ix", 1);
193  if ( i+1 >= 1 && i+1 <= 9 ) mePed3SumMapG06_[i]->setAxisTitle("101-ix", 1);
194  mePed3SumMapG06_[i]->setAxisTitle("iy", 2);
195  dqmStore_->tag(mePed3SumMapG06_[i], i+1);
196 
197  name = "EEPT pedestal 5sum " + Numbers::sEE(i+1) + " " + GN.str();
198  mePed5SumMapG06_[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");
199  mePed5SumMapG06_[i]->setAxisTitle("ix", 1);
200  if ( i+1 >= 1 && i+1 <= 9 ) mePed5SumMapG06_[i]->setAxisTitle("101-ix", 1);
201  mePed5SumMapG06_[i]->setAxisTitle("iy", 2);
202  dqmStore_->tag(mePed5SumMapG06_[i], i+1);
203 #endif
204  }
205 
206  }
207 
208  if (find(MGPAGains_.begin(), MGPAGains_.end(), 12) != MGPAGains_.end() ) {
209 
210  GainN.str("");
211  GainN << "Gain" << std::setw(2) << std::setfill('0') << 12;
212  GN.str("");
213  GN << "G" << std::setw(2) << std::setfill('0') << 12;
214 
215  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask/" + GainN.str());
216  for (int i = 0; i < 18; i++) {
217  name = "EEPT pedestal " + Numbers::sEE(i+1) + " " + GN.str();
218  mePedMapG12_[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");
219  mePedMapG12_[i]->setAxisTitle("ix", 1);
220  if ( i+1 >= 1 && i+1 <= 9 ) mePedMapG12_[i]->setAxisTitle("101-ix", 1);
221  mePedMapG12_[i]->setAxisTitle("iy", 2);
222  dqmStore_->tag(mePedMapG12_[i], i+1);
223 
224 #ifdef COMMON_NOISE_ANALYSIS
225  name = "EEPT pedestal 3sum " + Numbers::sEE(i+1) + " " + GN.str();
226  mePed3SumMapG12_[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");
227  mePed3SumMapG12_[i]->setAxisTitle("ix", 1);
228  if ( i+1 >= 1 && i+1 <= 9 ) mePed3SumMapG12_[i]->setAxisTitle("101-ix", 1);
229  mePed3SumMapG12_[i]->setAxisTitle("iy", 2);
230  dqmStore_->tag(mePed3SumMapG12_[i], i+1);
231 
232  name = "EEPT pedestal 5sum " + Numbers::sEE(i+1) + " " + GN.str();
233  mePed5SumMapG12_[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");
234  mePed5SumMapG12_[i]->setAxisTitle("ix", 1);
235  if ( i+1 >= 1 && i+1 <= 9 ) mePed5SumMapG12_[i]->setAxisTitle("101-ix", 1);
236  mePed5SumMapG12_[i]->setAxisTitle("iy", 2);
237  dqmStore_->tag(mePed5SumMapG12_[i], i+1);
238 #endif
239  }
240 
241  }
242 
243  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask/PN");
244 
245  if (find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), 1) != MGPAGainsPN_.end() ) {
246 
247  GainN.str("");
248  GainN << "Gain" << std::setw(2) << std::setfill('0') << 1;
249  GN.str("");
250  GN << "G" << std::setw(2) << std::setfill('0') << 1;
251 
252  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask/PN/" + GainN.str());
253  for (int i = 0; i < 18; i++) {
254  name = "EEPDT PNs pedestal " + Numbers::sEE(i+1) + " " + GN.str();
255  mePnPedMapG01_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
256  mePnPedMapG01_[i]->setAxisTitle("channel", 1);
257  mePnPedMapG01_[i]->setAxisTitle("pedestal", 2);
258  dqmStore_->tag(mePnPedMapG01_[i], i+1);
259  }
260 
261  }
262 
263  if (find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), 16) != MGPAGainsPN_.end() ) {
264 
265  GainN.str("");
266  GainN << "Gain" << std::setw(2) << std::setfill('0') << 16;
267  GN.str("");
268  GN << "G" << std::setw(2) << std::setfill('0') << 16;
269 
270  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask/PN/" + GainN.str());
271  for (int i = 0; i < 18; i++) {
272  name = "EEPDT PNs pedestal " + Numbers::sEE(i+1) + " " + GN.str();
273  mePnPedMapG16_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
274  mePnPedMapG16_[i]->setAxisTitle("channel", 1);
275  mePnPedMapG16_[i]->setAxisTitle("pedestal", 2);
276  dqmStore_->tag(mePnPedMapG16_[i], i+1);
277  }
278 
279  }
280 
281  }
282 
283 }
284 
285 void EEPedestalTask::cleanup(void){
286 
287  if ( ! init_ ) return;
288 
289  if ( dqmStore_ ) {
290  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask");
291 
292  if (find(MGPAGains_.begin(), MGPAGains_.end(), 1) != MGPAGains_.end() ) {
293 
294  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask/Gain01");
295  for ( int i = 0; i < 18; i++ ) {
296  if ( mePedMapG01_[i] ) dqmStore_->removeElement( mePedMapG01_[i]->getName() );
297  mePedMapG01_[i] = 0;
298 #ifdef COMMON_NOISE_ANALYSIS
299  if ( mePed3SumMapG01_[i] ) dqmStore_->removeElement( mePed3SumMapG01_[i]->getName() );
300  mePed3SumMapG01_[i] = 0;
301  if ( mePed5SumMapG01_[i] ) dqmStore_->removeElement( mePed5SumMapG01_[i]->getName() );
302  mePed5SumMapG01_[i] = 0;
303 #endif
304  }
305 
306  }
307 
308  if (find(MGPAGains_.begin(), MGPAGains_.end(), 6) != MGPAGains_.end() ) {
309 
310  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask/Gain06");
311  for ( int i = 0; i < 18; i++ ) {
312  if ( mePedMapG06_[i] ) dqmStore_->removeElement( mePedMapG06_[i]->getName() );
313  mePedMapG06_[i] = 0;
314 #ifdef COMMON_NOISE_ANALYSIS
315  if ( mePed3SumMapG06_[i] ) dqmStore_->removeElement( mePed3SumMapG06_[i]->getName() );
316  mePed3SumMapG06_[i] = 0;
317  if ( mePed5SumMapG06_[i] ) dqmStore_->removeElement( mePed5SumMapG06_[i]->getName() );
318  mePed5SumMapG06_[i] = 0;
319 #endif
320  }
321 
322  }
323 
324  if (find(MGPAGains_.begin(), MGPAGains_.end(), 12) != MGPAGains_.end() ) {
325 
326  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask/Gain12");
327  for ( int i = 0; i < 18; i++ ) {
328  if ( mePedMapG12_[i] ) dqmStore_->removeElement( mePedMapG12_[i]->getName() );
329  mePedMapG12_[i] = 0;
330 #ifdef COMMON_NOISE_ANALYSIS
331  if ( mePed3SumMapG12_[i] ) dqmStore_->removeElement( mePed3SumMapG12_[i]->getName() );
332  mePed3SumMapG12_[i] = 0;
333  if ( mePed5SumMapG12_[i] ) dqmStore_->removeElement( mePed5SumMapG12_[i]->getName() );
334  mePed5SumMapG12_[i] = 0;
335 #endif
336  }
337 
338  }
339 
340  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask/PN");
341 
342  if (find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), 1) != MGPAGainsPN_.end() ) {
343 
344  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask/PN/Gain01");
345  for ( int i = 0; i < 18; i++ ) {
346  if ( mePnPedMapG01_[i]) dqmStore_->removeElement( mePnPedMapG01_[i]->getName() );
347  mePnPedMapG01_[i] = 0;
348  }
349 
350  }
351 
352  if (find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), 16) != MGPAGainsPN_.end() ) {
353 
354  dqmStore_->setCurrentFolder(prefixME_ + "/EEPedestalTask/PN/Gain16");
355  for ( int i = 0; i < 18; i++ ) {
356  if ( mePnPedMapG16_[i]) dqmStore_->removeElement( mePnPedMapG16_[i]->getName() );
357  mePnPedMapG16_[i] = 0;
358  }
359 
360  }
361 
362  }
363 
364  init_ = false;
365 
366 }
367 
368 void EEPedestalTask::endJob(void){
369 
370  edm::LogInfo("EEPedestalTask") << "analyzed " << ievt_ << " events";
371 
372  if ( enableCleanup_ ) this->cleanup();
373 
374 }
375 
376 void EEPedestalTask::analyze(const edm::Event& e, const edm::EventSetup& c){
377 
378  bool enable = false;
379  int runType[18];
380  for (int i=0; i<18; i++) runType[i] = -1;
381 
382  edm::Handle<EcalRawDataCollection> dcchs;
383 
384  if ( e.getByToken(EcalRawDataCollection_, dcchs) ) {
385 
386  for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
387 
388  if ( Numbers::subDet( *dcchItr ) != EcalEndcap ) continue;
389 
390  int ism = Numbers::iSM( *dcchItr, EcalEndcap );
391 
392  runType[ism-1] = dcchItr->getRunType();
393 
394  if ( dcchItr->getRunType() == EcalDCCHeaderBlock::PEDESTAL_STD ||
395  dcchItr->getRunType() == EcalDCCHeaderBlock::PEDESTAL_GAP ) enable = true;
396 
397  }
398 
399  } else {
400 
401  edm::LogWarning("EEPedestalTask") << "EcalRawDataCollection not available";
402 
403  }
404 
405  if ( ! enable ) return;
406 
407  if ( ! init_ ) this->setup();
408 
409  ievt_++;
410 
411  edm::Handle<EEDigiCollection> digis;
412 
413  if ( e.getByToken(EEDigiCollection_, digis) ) {
414 
415  int need = digis->size();
416  LogDebug("EEPedestalTask") << "event " << ievt_ << " digi collection size " << need;
417 
418  float xmap01[18][50][50];
419  float xmap06[18][50][50];
420  float xmap12[18][50][50];
421 
422  for ( int ism = 1; ism <= 18; ism++ ) {
423  for ( int ix = 1; ix <= 50; ix++ ) {
424  for ( int iy = 1; iy <= 50; iy++ ) {
425 
426  xmap01[ism-1][ix-1][iy-1] = 0.;
427  xmap06[ism-1][ix-1][iy-1] = 0.;
428  xmap12[ism-1][ix-1][iy-1] = 0.;
429 
430  }
431  }
432  }
433 
434  for ( EEDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr ) {
435 
436  EEDetId id = digiItr->id();
437 
438  int ix = id.ix();
439  int iy = id.iy();
440 
441  int ism = Numbers::iSM( id );
442 
443  if ( ism >= 1 && ism <= 9 ) ix = 101 - ix;
444 
445  float xix = ix - 0.5;
446  float xiy = iy - 0.5;
447 
448  if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::PEDESTAL_STD ||
449  runType[ism-1] == EcalDCCHeaderBlock::PEDESTAL_GAP ) ) continue;
450 
451  EEDataFrame dataframe = (*digiItr);
452 
453  for (int i = 0; i < 10; i++) {
454 
455  int adc = dataframe.sample(i).adc();
456 
457  MonitorElement* mePedMap = 0;
458 
459  if ( dataframe.sample(i).gainId() == 1 ) mePedMap = mePedMapG12_[ism-1];
460  if ( dataframe.sample(i).gainId() == 2 ) mePedMap = mePedMapG06_[ism-1];
461  if ( dataframe.sample(i).gainId() == 3 ) mePedMap = mePedMapG01_[ism-1];
462 
463  float xval = float(adc);
464 
465  if ( mePedMap ) mePedMap->Fill(xix, xiy, xval);
466 
467  if ( dataframe.sample(i).gainId() == 1 ) xmap12[ism-1][ix-1-Numbers::ix0EE(ism)][iy-1-Numbers::iy0EE(ism)] = xmap12[ism-1][ix-1-Numbers::ix0EE(ism)][iy-1-Numbers::iy0EE(ism)] + xval;
468  if ( dataframe.sample(i).gainId() == 2 ) xmap06[ism-1][ix-1-Numbers::ix0EE(ism)][iy-1-Numbers::iy0EE(ism)] = xmap06[ism-1][ix-1-Numbers::ix0EE(ism)][iy-1-Numbers::iy0EE(ism)] + xval;
469  if ( dataframe.sample(i).gainId() == 3 ) xmap01[ism-1][ix-1-Numbers::ix0EE(ism)][iy-1-Numbers::iy0EE(ism)] = xmap01[ism-1][ix-1-Numbers::ix0EE(ism)][iy-1-Numbers::iy0EE(ism)] + xval;
470 
471  }
472 
473  xmap12[ism-1][ix-1-Numbers::ix0EE(ism)][iy-1-Numbers::iy0EE(ism)]=xmap12[ism-1][ix-1-Numbers::ix0EE(ism)][iy-1-Numbers::iy0EE(ism)]/10.;
474  xmap06[ism-1][ix-1-Numbers::ix0EE(ism)][iy-1-Numbers::iy0EE(ism)]=xmap06[ism-1][ix-1-Numbers::ix0EE(ism)][iy-1-Numbers::iy0EE(ism)]/10.;
475  xmap01[ism-1][ix-1-Numbers::ix0EE(ism)][iy-1-Numbers::iy0EE(ism)]=xmap01[ism-1][ix-1-Numbers::ix0EE(ism)][iy-1-Numbers::iy0EE(ism)]/10.;
476 
477  }
478 
479  // to be re-done using the 3x3 & 5x5 Selectors (if faster)
480 
481 #ifdef COMMON_NOISE_ANALYSIS
482  for ( int ism = 1; ism <= 18; ism++ ) {
483  for ( int ix = 1; ix <= 50; ix++ ) {
484  for ( int iy = 1; iy <= 50; iy++ ) {
485 
486  int xix = ix + Numbers::ix0EE(ism);
487  int xiy = iy + Numbers::iy0EE(ism);
488 
489  if ( ism >= 1 && ism <= 9 ) xix = 101 - xix;
490 
491  if ( ! Numbers::validEE(ism, xix, xiy) ) continue;
492 
493  float xix = ix - 0.5;
494  float xiy = iy - 0.5;
495 
496  float x3val01;
497  float x3val06;
498  float x3val12;
499 
500  if ( ix >= 2 && ix <= 49 && iy >= 2 && iy <= 49 ) {
501 
502  x3val01 = 0.;
503  x3val06 = 0.;
504  x3val12 = 0.;
505  for ( int i = -1; i <= +1; i++ ) {
506  for ( int j = -1; j <= +1; j++ ) {
507 
508  x3val01 = x3val01 + xmap01[ism-1][ix-1+i][iy-1+j];
509  x3val06 = x3val06 + xmap06[ism-1][ix-1+i][iy-1+j];
510  x3val12 = x3val12 + xmap12[ism-1][ix-1+i][iy-1+j];
511 
512  }
513  }
514  x3val01 = x3val01 / 9.;
515  x3val06 = x3val06 / 9.;
516  x3val12 = x3val12 / 9.;
517  if ( mePed3SumMapG01_[ism-1] && x3val01 != 0. ) mePed3SumMapG01_[ism-1]->Fill(xix+Numbers::ix0EE(ism), xiy+Numbers::iy0EE(ism), x3val01);
518  if ( mePed3SumMapG06_[ism-1] && x3val06 != 0. ) mePed3SumMapG06_[ism-1]->Fill(xix+Numbers::ix0EE(ism), xiy+Numbers::iy0EE(ism), x3val06);
519  if ( mePed3SumMapG12_[ism-1] && x3val12 != 0. ) mePed3SumMapG12_[ism-1]->Fill(xix+Numbers::ix0EE(ism), xiy+Numbers::iy0EE(ism), x3val12);
520 
521  }
522 
523  float x5val01;
524  float x5val06;
525  float x5val12;
526 
527  if ( ix >= 3 && ix <= 48 && iy >= 3 && iy <= 48 ) {
528 
529  x5val01 = 0.;
530  x5val06 = 0.;
531  x5val12 = 0.;
532  for ( int i = -2; i <= +2; i++ ) {
533  for ( int j = -2; j <= +2; j++ ) {
534 
535  x5val01 = x5val01 + xmap01[ism-1][ix-1+i][iy-1+j];
536  x5val06 = x5val06 + xmap06[ism-1][ix-1+i][iy-1+j];
537  x5val12 = x5val12 + xmap12[ism-1][ix-1+i][iy-1+j];
538 
539  }
540  }
541  x5val01 = x5val01 / 25.;
542  x5val06 = x5val06 / 25.;
543  x5val12 = x5val12 / 25.;
544  if ( mePed5SumMapG01_[ism-1] && x5val01 != 0. ) mePed5SumMapG01_[ism-1]->Fill(xix+Numbers::ix0EE(ism), xiy+Numbers::iy0EE(ism), x5val01);
545  if ( mePed5SumMapG06_[ism-1] && x5val06 != 0. ) mePed5SumMapG06_[ism-1]->Fill(xix+Numbers::ix0EE(ism), xiy+Numbers::iy0EE(ism), x5val06);
546  if ( mePed5SumMapG12_[ism-1] && x5val12 != 0. ) mePed5SumMapG12_[ism-1]->Fill(xix+Numbers::ix0EE(ism), xiy+Numbers::iy0EE(ism), x5val12);
547 
548  }
549 
550  }
551  }
552  }
553 #endif
554 
555  } else {
556 
557  edm::LogWarning("EEPedestalTask") << "EEDigiCollection not available";
558 
559  }
560 
561  edm::Handle<EcalPnDiodeDigiCollection> pns;
562 
563  if ( e.getByToken(EcalPnDiodeDigiCollection_, pns) ) {
564 
565  int nep = pns->size();
566  LogDebug("EEPedestalTask") << "event " << ievt_ << " pns collection size " << nep;
567 
568  for ( EcalPnDiodeDigiCollection::const_iterator pnItr = pns->begin(); pnItr != pns->end(); ++pnItr ) {
569 
570  if ( Numbers::subDet( pnItr->id() ) != EcalEndcap ) continue;
571 
572  int ism = Numbers::iSM( pnItr->id() );
573 
574  int num = pnItr->id().iPnId();
575 
576  if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::PEDESTAL_STD ||
577  runType[ism-1] == EcalDCCHeaderBlock::PEDESTAL_GAP ) ) continue;
578 
579  for (int i = 0; i < 50; i++) {
580 
581  int adc = pnItr->sample(i).adc();
582 
583  MonitorElement* mePNPed = 0;
584 
585  if ( pnItr->sample(i).gainId() == 0 ) mePNPed = mePnPedMapG01_[ism-1];
586  if ( pnItr->sample(i).gainId() == 1 ) mePNPed = mePnPedMapG16_[ism-1];
587 
588  float xval = float(adc);
589 
590  if ( mePNPed ) mePNPed->Fill(num - 0.5, xval);
591 
592  }
593 
594  }
595 
596  } else {
597 
598  edm::LogWarning("EEPedestalTask") << "EcalPnDiodeDigiCollection not available";
599 
600  }
601 
602 }
603 
ecalMGPA::adc
int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:11
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
EEPedestalTask::mePedMapG06_
MonitorElement * mePedMapG06_[18]
Definition: EEPedestalTask.h:76
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
Numbers::sEE
static std::string sEE(const unsigned ism)
Definition: Numbers.cc:223
i
int i
Definition: DBlmapReader.cc:9
Numbers::validEE
static bool validEE(const unsigned ism, const unsigned ix, const unsigned iy)
Definition: Numbers.cc:829
EEPedestalTask::mePedMapG12_
MonitorElement * mePedMapG12_[18]
Definition: EEPedestalTask.h:77
EEDetId::ix
int ix() const
Definition: EEDetId.h:76
EEPedestalTask::cleanup
void cleanup(void)
Cleanup.
Definition: EEPedestalTask.cc:285
edm::DataFrameContainer::const_iterator
boost::transform_iterator< IterHelp, boost::counting_iterator< int > > const_iterator
Definition: DataFrameContainer.h:51
DQMStore::rmdir
void rmdir(const std::string &fullpath)
Definition: DQMStore.cc:2730
EEPedestalTask::mePedMapG01_
MonitorElement * mePedMapG01_[18]
Definition: EEPedestalTask.h:75
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434
Numbers.h
Some &quot;id&quot; conversions.
EEPedestalTask::setup
void setup(void)
Setup.
Definition: EEPedestalTask.cc:128
edm::SortedCollection< EcalDCCHeaderBlock >::const_iterator
std::vector< EcalDCCHeaderBlock >::const_iterator const_iterator
Definition: SortedCollection.h:80
EcalDataFrame::sample
EcalMGPASample sample(int i) const
Definition: EcalDataFrame.h:29
Numbers::ix0EE
static int ix0EE(const unsigned ism)
Definition: Numbers.cc:770
EEPedestalTask::mePnPedMapG16_
MonitorElement * mePnPedMapG16_[18]
Definition: EEPedestalTask.h:88
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
EcalMGPASample::gainId
int gainId() const
get the gainId (2 bits)
Definition: EcalMGPASample.h:36
EEPedestalTask::reset
void reset(void)
Reset.
Definition: EEPedestalTask.cc:98
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
EEPedestalTask::EEPedestalTask
EEPedestalTask(const edm::ParameterSet &ps)
Constructor.
Definition: EEPedestalTask.cc:29
EEPedestalTask::endJob
void endJob(void)
EndJob.
Definition: EEPedestalTask.cc:368
EEPedestalTask::beginJob
void beginJob(void)
BeginJob.
Definition: EEPedestalTask.cc:75
EEPedestalTask::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Analyze.
Definition: EEPedestalTask.cc:376
EEPedestalTask::endRun
void endRun(const edm::Run &r, const edm::EventSetup &c)
EndRun.
Definition: EEPedestalTask.cc:94
EEPedestalTask::~EEPedestalTask
virtual ~EEPedestalTask()
Destructor.
Definition: EEPedestalTask.cc:71
DQMStore::removeElement
void removeElement(const std::string &name)
Definition: DQMStore.cc:2772
EEPedestalTask::mePed3SumMapG06_
MonitorElement * mePed3SumMapG06_[18]
Definition: EEPedestalTask.h:80
j
int j
Definition: DBlmapReader.cc:9
EEPedestalTask::mePed5SumMapG06_
MonitorElement * mePed5SumMapG06_[18]
Definition: EEPedestalTask.h:84
EEPedestalTask::mePed3SumMapG01_
MonitorElement * mePed3SumMapG01_[18]
Definition: EEPedestalTask.h:79
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
EEPedestalTask::MGPAGains_
std::vector< int > MGPAGains_
Definition: EEPedestalTask.h:72
EEPedestalTask::prefixME_
std::string prefixME_
Definition: EEPedestalTask.h:63
environment_file_cff.runType
tuple runType
Definition: environment_file_cff.py:30
EEPedestalTask::mePed5SumMapG12_
MonitorElement * mePed5SumMapG12_[18]
Definition: EEPedestalTask.h:85
Numbers::initGeometry
static void initGeometry(const edm::EventSetup &setup, bool verbose=false)
Definition: Numbers.cc:47
EEPedestalTask::mePnPedMapG01_
MonitorElement * mePnPedMapG01_[18]
Definition: EEPedestalTask.h:87
EEPedestalTask::MGPAGainsPN_
std::vector< int > MGPAGainsPN_
Definition: EEPedestalTask.h:73
EEPedestalTask::mePed5SumMapG01_
MonitorElement * mePed5SumMapG01_[18]
Definition: EEPedestalTask.h:83
EEPedestalTask::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
BeginRun.
Definition: EEPedestalTask.cc:86
Numbers::iSM
static unsigned iSM(const unsigned ism, const EcalSubdetector subdet)
Definition: Numbers.cc:243
EEPedestalTask::EcalRawDataCollection_
edm::EDGetTokenT< EcalRawDataCollection > EcalRawDataCollection_
Definition: EEPedestalTask.h:69
EEPedestalTask::EEDigiCollection_
edm::EDGetTokenT< EEDigiCollection > EEDigiCollection_
Definition: EEPedestalTask.h:70
Numbers::subDet
static EcalSubdetector subDet(const EBDetId &id)
Definition: Numbers.cc:142
EEPedestalTask::dqmStore_
DQMStore * dqmStore_
Definition: EEPedestalTask.h:61
ecalpyutils::ism
int ism(int ieta, int iphi)
Definition: EcalPyUtils.cc:56
EEPedestalTask::mePed3SumMapG12_
MonitorElement * mePed3SumMapG12_[18]
Definition: EEPedestalTask.h:81
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
EEPedestalTask::EcalPnDiodeDigiCollection_
edm::EDGetTokenT< EcalPnDiodeDigiCollection > EcalPnDiodeDigiCollection_
Definition: EEPedestalTask.h:71
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:584
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