CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EEIntegrityTask.cc
Go to the documentation of this file.
1 /*
2  * \file EEIntegrityTask.cc
3  *
4  * $Date: 2012/04/27 13:46:15 $
5  * $Revision: 1.61 $
6  * \author G. Della Ricca
7  *
8  */
9 
10 #include <iostream>
11 #include <fstream>
12 #include <vector>
13 
14 #include "FWCore/ServiceRegistry/interface/Service.h"
15 #include "FWCore/MessageLogger/interface/MessageLogger.h"
16 
17 #include "DQMServices/Core/interface/MonitorElement.h"
18 
19 #include "DQMServices/Core/interface/DQMStore.h"
20 
21 #include "DataFormats/EcalDetId/interface/EEDetId.h"
22 
23 #include "DataFormats/EcalDetId/interface/EcalDetIdCollections.h"
24 
25 #include "DQM/EcalCommon/interface/Numbers.h"
26 
27 #include "DQM/EcalEndcapMonitorTasks/interface/EEIntegrityTask.h"
28 
29 EEIntegrityTask::EEIntegrityTask(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  subfolder_ = ps.getUntrackedParameter<std::string>("subfolder", "");
38 
39  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
40 
41  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
42 
43  EEDetIdCollection0_ = ps.getParameter<edm::InputTag>("EEDetIdCollection0");
44  EEDetIdCollection1_ = ps.getParameter<edm::InputTag>("EEDetIdCollection1");
45  EEDetIdCollection2_ = ps.getParameter<edm::InputTag>("EEDetIdCollection2");
46  EEDetIdCollection3_ = ps.getParameter<edm::InputTag>("EEDetIdCollection3");
47  EcalElectronicsIdCollection1_ = ps.getParameter<edm::InputTag>("EcalElectronicsIdCollection1");
48  EcalElectronicsIdCollection2_ = ps.getParameter<edm::InputTag>("EcalElectronicsIdCollection2");
49  EcalElectronicsIdCollection3_ = ps.getParameter<edm::InputTag>("EcalElectronicsIdCollection3");
50  EcalElectronicsIdCollection4_ = ps.getParameter<edm::InputTag>("EcalElectronicsIdCollection4");
51  EcalElectronicsIdCollection5_ = ps.getParameter<edm::InputTag>("EcalElectronicsIdCollection5");
52  EcalElectronicsIdCollection6_ = ps.getParameter<edm::InputTag>("EcalElectronicsIdCollection6");
53 
54  meIntegrityDCCSize = 0;
55  for (int i = 0; i < 18; i++) {
56  meIntegrityGain[i] = 0;
57  meIntegrityChId[i] = 0;
58  meIntegrityGainSwitch[i] = 0;
59  meIntegrityTTId[i] = 0;
60  meIntegrityTTBlockSize[i] = 0;
61  meIntegrityMemChId[i] = 0;
62  meIntegrityMemGain[i] = 0;
63  meIntegrityMemTTId[i] = 0;
64  meIntegrityMemTTBlockSize[i] = 0;
65  }
66  meIntegrityErrorsByLumi = 0;
67 }
68 
69 
70 EEIntegrityTask::~EEIntegrityTask(){
71 
72 }
73 
74 void EEIntegrityTask::beginJob(void){
75 
76  ievt_ = 0;
77 
78  if ( dqmStore_ ) {
79  dqmStore_->setCurrentFolder(prefixME_ + "/EEIntegrityTask");
80  if(subfolder_.size())
81  dqmStore_->setCurrentFolder(prefixME_ + "/EEIntegrityTask/" + subfolder_);
82  dqmStore_->rmdir(prefixME_ + "/EEIntegrityTask");
83  }
84 
85 }
86 
87 void EEIntegrityTask::beginLuminosityBlock(const edm::LuminosityBlock& lumiBlock, const edm::EventSetup& iSetup) {
88 
89  if ( meIntegrityErrorsByLumi ) meIntegrityErrorsByLumi->Reset();
90 
91 }
92 
93 void EEIntegrityTask::endLuminosityBlock(const edm::LuminosityBlock& lumiBlock, const edm::EventSetup& iSetup) {
94 }
95 
96 void EEIntegrityTask::beginRun(const edm::Run& r, const edm::EventSetup& c) {
97 
98  Numbers::initGeometry(c, false);
99 
100  if ( ! mergeRuns_ ) this->reset();
101 
102 }
103 
104 void EEIntegrityTask::endRun(const edm::Run& r, const edm::EventSetup& c) {
105 
106 }
107 
108 void EEIntegrityTask::reset(void) {
109 
110  if ( meIntegrityDCCSize ) meIntegrityDCCSize->Reset();
111  for (int i = 0; i < 18; i++) {
112  if ( meIntegrityGain[i] ) meIntegrityGain[i]->Reset();
113  if ( meIntegrityChId[i] ) meIntegrityChId[i]->Reset();
114  if ( meIntegrityGainSwitch[i] ) meIntegrityGainSwitch[i]->Reset();
115  if ( meIntegrityTTId[i] ) meIntegrityTTId[i]->Reset();
116  if ( meIntegrityTTBlockSize[i] ) meIntegrityTTBlockSize[i]->Reset();
117  if ( meIntegrityMemChId[i] ) meIntegrityMemChId[i]->Reset();
118  if ( meIntegrityMemGain[i] ) meIntegrityMemGain[i]->Reset();
119  if ( meIntegrityMemTTId[i] ) meIntegrityMemTTId[i]->Reset();
120  if ( meIntegrityMemTTBlockSize[i] ) meIntegrityMemTTBlockSize[i]->Reset();
121  }
122  if ( meIntegrityErrorsByLumi ) meIntegrityErrorsByLumi->Reset();
123 
124 }
125 
126 void EEIntegrityTask::setup(void){
127 
128  init_ = true;
129 
130  std::string name;
131  std::string dir;
132 
133  if ( dqmStore_ ) {
134 
135  dir = prefixME_ + "/EEIntegrityTask";
136  if(subfolder_.size())
137  dir += "/" + subfolder_;
138 
139  dqmStore_->setCurrentFolder(dir);
140 
141  // checking when number of towers in data different than expected from header
142  name = "EEIT DCC size error";
143  meIntegrityDCCSize = dqmStore_->book1D(name, name, 18, 1., 19.);
144  for (int i = 0; i < 18; i++) {
145  meIntegrityDCCSize->setBinLabel(i+1, Numbers::sEE(i+1).c_str(), 1);
146  }
147 
148  // checking the number of integrity errors in each DCC for each lumi
149  // crystal integrity error is weighted by 1/850
150  // tower integrity error is weighted by 1/34
151  // bin 0 contains the number of processed events in the lumi (for normalization)
152  name = "EEIT weighted integrity errors by lumi";
153  meIntegrityErrorsByLumi = dqmStore_->book1D(name, name, 18, 1., 19.);
154  meIntegrityErrorsByLumi->setLumiFlag();
155  for (int i = 0; i < 18; i++) {
156  meIntegrityErrorsByLumi->setBinLabel(i+1, Numbers::sEE(i+1).c_str(), 1);
157  }
158 
159  // checking when the gain is 0
160  dqmStore_->setCurrentFolder(dir + "/Gain");
161  for (int i = 0; i < 18; i++) {
162  name = "EEIT gain " + Numbers::sEE(i+1);
163  meIntegrityGain[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.);
164  meIntegrityGain[i]->setAxisTitle("ix", 1);
165  if ( i+1 >= 1 && i+1 <= 9 ) meIntegrityGain[i]->setAxisTitle("101-ix", 1);
166  meIntegrityGain[i]->setAxisTitle("iy", 2);
167  dqmStore_->tag(meIntegrityGain[i], i+1);
168  }
169 
170  // checking when channel has unexpected or invalid ID
171  dqmStore_->setCurrentFolder(dir + "/ChId");
172  for (int i = 0; i < 18; i++) {
173  name = "EEIT ChId " + Numbers::sEE(i+1);
174  meIntegrityChId[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.);
175  meIntegrityChId[i]->setAxisTitle("ix", 1);
176  if ( i+1 >= 1 && i+1 <= 9 ) meIntegrityChId[i]->setAxisTitle("101-ix", 1);
177  meIntegrityChId[i]->setAxisTitle("iy", 2);
178  dqmStore_->tag(meIntegrityChId[i], i+1);
179  }
180 
181  // checking when channel has unexpected or invalid ID
182  dqmStore_->setCurrentFolder(dir + "/GainSwitch");
183  for (int i = 0; i < 18; i++) {
184  name = "EEIT gain switch " + Numbers::sEE(i+1);
185  meIntegrityGainSwitch[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.);
186  meIntegrityGainSwitch[i]->setAxisTitle("ix", 1);
187  if ( i+1 >= 1 && i+1 <= 9 ) meIntegrityGainSwitch[i]->setAxisTitle("101-ix", 1);
188  meIntegrityGainSwitch[i]->setAxisTitle("iy", 2);
189  dqmStore_->tag(meIntegrityGainSwitch[i], i+1);
190  }
191 
192  // checking when trigger tower has unexpected or invalid ID
193  dqmStore_->setCurrentFolder(dir + "/TTId");
194  for (int i = 0; i < 18; i++) {
195  name = "EEIT TTId " + Numbers::sEE(i+1);
196  meIntegrityTTId[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.);
197  meIntegrityTTId[i]->setAxisTitle("ix", 1);
198  if ( i+1 >= 1 && i+1 <= 9 ) meIntegrityTTId[i]->setAxisTitle("101-ix", 1);
199  meIntegrityTTId[i]->setAxisTitle("iy", 2);
200  dqmStore_->tag(meIntegrityTTId[i], i+1);
201  }
202 
203  // checking when trigger tower has unexpected or invalid size
204  dqmStore_->setCurrentFolder(dir + "/TTBlockSize");
205  for (int i = 0; i < 18; i++) {
206  name = "EEIT TTBlockSize " + Numbers::sEE(i+1);
207  meIntegrityTTBlockSize[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.);
208  meIntegrityTTBlockSize[i]->setAxisTitle("ix", 1);
209  if ( i+1 >= 1 && i+1 <= 9 ) meIntegrityTTBlockSize[i]->setAxisTitle("101-ix", 1);
210  meIntegrityTTBlockSize[i]->setAxisTitle("iy", 2);
211  dqmStore_->tag(meIntegrityTTBlockSize[i], i+1);
212  }
213 
214  // checking when mem channels have unexpected ID
215  dqmStore_->setCurrentFolder(dir + "/MemChId");
216  for (int i = 0; i < 18; i++) {
217  name = "EEIT MemChId " + Numbers::sEE(i+1);
218  meIntegrityMemChId[i] = dqmStore_->book2D(name, name, 10, 0., 10., 5, 0., 5.);
219  meIntegrityMemChId[i]->setAxisTitle("pseudo-strip", 1);
220  meIntegrityMemChId[i]->setAxisTitle("channel", 2);
221  dqmStore_->tag(meIntegrityMemChId[i], i+1);
222  }
223 
224  // checking when mem samples have second bit encoding the gain different from 0
225  // note: strictly speaking, this does not corrupt the mem sample gain value (since only first bit is considered)
226  // but indicates that data are not completely correct
227  dqmStore_->setCurrentFolder(dir + "/MemGain");
228  for (int i = 0; i < 18; i++) {
229  name = "EEIT MemGain " + Numbers::sEE(i+1);
230  meIntegrityMemGain[i] = dqmStore_->book2D(name, name, 10, 0., 10., 5, 0., 5.);
231  meIntegrityMemGain[i]->setAxisTitle("pseudo-strip", 1);
232  meIntegrityMemGain[i]->setAxisTitle("channel", 2);
233  dqmStore_->tag(meIntegrityMemGain[i], i+1);
234  }
235 
236  // checking when mem tower block has unexpected ID
237  dqmStore_->setCurrentFolder(dir + "/MemTTId");
238  for (int i = 0; i < 18; i++) {
239  name = "EEIT MemTTId " + Numbers::sEE(i+1);
240  meIntegrityMemTTId[i] = dqmStore_->book2D(name, name, 2, 0., 2., 1, 0., 1.);
241  meIntegrityMemTTId[i]->setAxisTitle("pseudo-strip", 1);
242  meIntegrityMemTTId[i]->setAxisTitle("channel", 2);
243  dqmStore_->tag(meIntegrityMemTTId[i], i+1);
244  }
245 
246  // checking when mem tower block has invalid size
247  dqmStore_->setCurrentFolder(dir + "/MemSize");
248  for (int i = 0; i < 18; i++) {
249  name = "EEIT MemSize " + Numbers::sEE(i+1);
250  meIntegrityMemTTBlockSize[i] = dqmStore_->book2D(name, name, 2, 0., 2., 1, 0., 1.);
251  meIntegrityMemTTBlockSize[i]->setAxisTitle("pseudo-strip", 1);
252  meIntegrityMemTTBlockSize[i]->setAxisTitle("channel", 2);
253  dqmStore_->tag(meIntegrityMemTTBlockSize[i], i+1);
254  }
255 
256  }
257 
258 }
259 
260 void EEIntegrityTask::cleanup(void){
261 
262  if ( ! init_ ) return;
263 
264  if ( dqmStore_ ) {
265 
266  std::string dir;
267 
268  dir = prefixME_ + "/EEIntegrityTask";
269  if(subfolder_.size())
270  dir += "/" + subfolder_;
271 
272  dqmStore_->setCurrentFolder(dir + "");
273 
274  if ( meIntegrityDCCSize ) dqmStore_->removeElement( meIntegrityDCCSize->getName() );
275  meIntegrityDCCSize = 0;
276 
277  if ( meIntegrityErrorsByLumi ) dqmStore_->removeElement( meIntegrityErrorsByLumi->getName() );
278  meIntegrityErrorsByLumi = 0;
279 
280  dqmStore_->setCurrentFolder(dir + "/Gain");
281  for (int i = 0; i < 18; i++) {
282  if ( meIntegrityGain[i] ) dqmStore_->removeElement( meIntegrityGain[i]->getName() );
283  meIntegrityGain[i] = 0;
284  }
285 
286  dqmStore_->setCurrentFolder(dir + "/ChId");
287  for (int i = 0; i < 18; i++) {
288  if ( meIntegrityChId[i] ) dqmStore_->removeElement( meIntegrityChId[i]->getName() );
289  meIntegrityChId[i] = 0;
290  }
291 
292  dqmStore_->setCurrentFolder(dir + "/GainSwitch");
293  for (int i = 0; i < 18; i++) {
294  if ( meIntegrityGainSwitch[i] ) dqmStore_->removeElement( meIntegrityGainSwitch[i]->getName() );
295  meIntegrityGainSwitch[i] = 0;
296  }
297 
298  dqmStore_->setCurrentFolder(dir + "/TTId");
299  for (int i = 0; i < 18; i++) {
300  if ( meIntegrityTTId[i] ) dqmStore_->removeElement( meIntegrityTTId[i]->getName() );
301  meIntegrityTTId[i] = 0;
302  }
303 
304  dqmStore_->setCurrentFolder(dir + "/TTBlockSize");
305  for (int i = 0; i < 18; i++) {
306  if ( meIntegrityTTBlockSize[i] ) dqmStore_->removeElement( meIntegrityTTBlockSize[i]->getName() );
307  meIntegrityTTBlockSize[i] = 0;
308  }
309 
310  dqmStore_->setCurrentFolder(dir + "/MemChId");
311  for (int i = 0; i < 18; i++) {
312  if ( meIntegrityMemChId[i] ) dqmStore_->removeElement( meIntegrityMemChId[i]->getName() );
313  meIntegrityMemChId[i] = 0;
314  }
315 
316  dqmStore_->setCurrentFolder(dir + "/MemGain");
317  for (int i = 0; i < 18; i++) {
318  if ( meIntegrityMemGain[i] ) dqmStore_->removeElement( meIntegrityMemGain[i]->getName() );
319  meIntegrityMemGain[i] = 0;
320  }
321 
322  dqmStore_->setCurrentFolder(dir + "/MemTTId");
323  for (int i = 0; i < 18; i++) {
324  if ( meIntegrityMemTTId[i] ) dqmStore_->removeElement( meIntegrityMemTTId[i]->getName() );
325  meIntegrityMemTTId[i] = 0;
326  }
327 
328  dqmStore_->setCurrentFolder(dir + "/MemSize");
329  for (int i = 0; i < 18; i++) {
330  if ( meIntegrityMemTTBlockSize[i] ) dqmStore_->removeElement( meIntegrityMemTTBlockSize[i]->getName() );
331  meIntegrityMemTTBlockSize[i] = 0;
332  }
333 
334  }
335 
336  init_ = false;
337 
338 }
339 
340 void EEIntegrityTask::endJob(void){
341 
342  edm::LogInfo("EEIntegrityTask") << "analyzed " << ievt_ << " events";
343 
344  if ( enableCleanup_ ) this->cleanup();
345 
346 }
347 
348 void EEIntegrityTask::analyze(const edm::Event& e, const edm::EventSetup& c){
349 
350  if ( ! init_ ) this->setup();
351 
352  ievt_++;
353 
354  // fill bin 0 with number of events in the lumi
355  if ( meIntegrityErrorsByLumi ) meIntegrityErrorsByLumi->Fill(0.);
356 
357  edm::Handle<EEDetIdCollection> ids0;
358 
359  if ( e.getByLabel(EEDetIdCollection0_, ids0) ) {
360 
361  for ( EEDetIdCollection::const_iterator idItr = ids0->begin(); idItr != ids0->end(); ++idItr ) {
362 
363  int ism = Numbers::iSM( *idItr );
364 
365  float xism = ism + 0.5;
366 
367  if ( meIntegrityDCCSize ) meIntegrityDCCSize->Fill(xism);
368 
369  }
370 
371  } else {
372 
373 // edm::LogWarning("EEIntegrityTask") << EEDetIdCollection0_ << " not available";
374 
375  }
376 
377  edm::Handle<EEDetIdCollection> ids1;
378 
379  if ( e.getByLabel(EEDetIdCollection1_, ids1) ) {
380 
381  for ( EEDetIdCollection::const_iterator idItr = ids1->begin(); idItr != ids1->end(); ++idItr ) {
382 
383  EEDetId id = (*idItr);
384 
385  int ix = id.ix();
386  int iy = id.iy();
387 
388  int ism = Numbers::iSM( id );
389  float xism = ism + 0.5;
390 
391  if ( ism >= 1 && ism <= 9 ) ix = 101 - ix;
392 
393  float xix = ix - 0.5;
394  float xiy = iy - 0.5;
395 
396  if ( meIntegrityGain[ism-1] ) meIntegrityGain[ism-1]->Fill(xix, xiy);
397  if ( meIntegrityErrorsByLumi ) meIntegrityErrorsByLumi->Fill(xism, 1./850.);
398 
399  }
400 
401  } else {
402 
403  edm::LogWarning("EEIntegrityTask") << EEDetIdCollection1_ << " not available";
404 
405  }
406 
407  edm::Handle<EEDetIdCollection> ids2;
408 
409  if ( e.getByLabel(EEDetIdCollection2_, ids2) ) {
410 
411  for ( EEDetIdCollection::const_iterator idItr = ids2->begin(); idItr != ids2->end(); ++idItr ) {
412 
413  EEDetId id = (*idItr);
414 
415  int ix = id.ix();
416  int iy = id.iy();
417 
418  int ism = Numbers::iSM( id );
419  float xism = ism + 0.5;
420 
421  if ( ism >= 1 && ism <= 9 ) ix = 101 - ix;
422 
423  float xix = ix - 0.5;
424  float xiy = iy - 0.5;
425 
426  if ( meIntegrityChId[ism-1] ) meIntegrityChId[ism-1]->Fill(xix, xiy);
427  if ( meIntegrityErrorsByLumi ) meIntegrityErrorsByLumi->Fill(xism, 1./850.);
428 
429  }
430 
431  } else {
432 
433  edm::LogWarning("EEIntegrityTask") << EEDetIdCollection2_ << " not available";
434 
435  }
436 
437  edm::Handle<EEDetIdCollection> ids3;
438 
439  if ( e.getByLabel(EEDetIdCollection3_, ids3) ) {
440 
441  for ( EEDetIdCollection::const_iterator idItr = ids3->begin(); idItr != ids3->end(); ++idItr ) {
442 
443  EEDetId id = (*idItr);
444 
445  int ix = id.ix();
446  int iy = id.iy();
447 
448  int ism = Numbers::iSM( id );
449  float xism = ism + 0.5;
450 
451  if ( ism >= 1 && ism <= 9 ) ix = 101 - ix;
452 
453  float xix = ix - 0.5;
454  float xiy = iy - 0.5;
455 
456  if ( meIntegrityGainSwitch[ism-1] ) meIntegrityGainSwitch[ism-1]->Fill(xix, xiy);
457  if ( meIntegrityErrorsByLumi ) meIntegrityErrorsByLumi->Fill(xism, 1./850.);
458 
459  }
460 
461  } else {
462 
463  edm::LogWarning("EEIntegrityTask") << EEDetIdCollection3_ << " not available";
464 
465  }
466 
467  edm::Handle<EcalElectronicsIdCollection> ids4;
468 
469  if ( e.getByLabel(EcalElectronicsIdCollection1_, ids4) ) {
470 
471  for ( EcalElectronicsIdCollection::const_iterator idItr = ids4->begin(); idItr != ids4->end(); ++idItr ) {
472 
473  if ( Numbers::subDet( *idItr ) != EcalEndcap ) continue;
474 
475  int ism = Numbers::iSM( *idItr );
476  float xism = ism + 0.5;
477 
478  std::vector<DetId>* crystals = Numbers::crystals( *idItr );
479 
480  for ( unsigned int i=0; i<crystals->size(); i++ ) {
481 
482  EEDetId id = (*crystals)[i];
483 
484  int ix = id.ix();
485  int iy = id.iy();
486 
487  if ( ism >= 1 && ism <= 9 ) ix = 101 - ix;
488 
489  float xix = ix - 0.5;
490  float xiy = iy - 0.5;
491 
492  if ( meIntegrityTTId[ism-1] ) meIntegrityTTId[ism-1]->Fill(xix, xiy);
493  if ( meIntegrityErrorsByLumi ) meIntegrityErrorsByLumi->Fill(xism, 1./34./crystals->size());
494 
495  }
496 
497  }
498 
499  } else {
500 
501  edm::LogWarning("EEIntegrityTask") << EcalElectronicsIdCollection1_ << " not available";
502 
503  }
504 
505  edm::Handle<EcalElectronicsIdCollection> ids5;
506 
507  if ( e.getByLabel(EcalElectronicsIdCollection2_, ids5) ) {
508 
509  for ( EcalElectronicsIdCollection::const_iterator idItr = ids5->begin(); idItr != ids5->end(); ++idItr ) {
510 
511  if ( Numbers::subDet( *idItr ) != EcalEndcap ) continue;
512 
513  int ism = Numbers::iSM( *idItr );
514  float xism = ism + 0.5;
515 
516  std::vector<DetId>* crystals = Numbers::crystals( *idItr );
517 
518  for ( unsigned int i=0; i<crystals->size(); i++ ) {
519 
520  EEDetId id = (*crystals)[i];
521 
522  int ix = id.ix();
523  int iy = id.iy();
524 
525  if ( ism >= 1 && ism <= 9 ) ix = 101 - ix;
526 
527  float xix = ix - 0.5;
528  float xiy = iy - 0.5;
529 
530  if ( meIntegrityTTBlockSize[ism-1] ) meIntegrityTTBlockSize[ism-1]->Fill(xix, xiy);
531  if ( meIntegrityErrorsByLumi ) meIntegrityErrorsByLumi->Fill(xism, 1./34./crystals->size());
532 
533  }
534 
535  }
536 
537  } else {
538 
539  edm::LogWarning("EEIntegrityTask") << EcalElectronicsIdCollection2_ << " not available";
540 
541  }
542 
543  edm::Handle<EcalElectronicsIdCollection> ids6;
544 
545  if ( e.getByLabel(EcalElectronicsIdCollection3_, ids6) ) {
546 
547  for ( EcalElectronicsIdCollection::const_iterator idItr = ids6->begin(); idItr != ids6->end(); ++idItr ) {
548 
549  if ( Numbers::subDet( *idItr ) != EcalEndcap ) continue;
550 
551  int ism = Numbers::iSM( *idItr );
552 
553  int itt = idItr->towerId();
554  float iTt = itt + 0.5 - 69;
555 
556  if ( meIntegrityMemTTId[ism-1] ) meIntegrityMemTTId[ism-1]->Fill(iTt,0);
557 
558  }
559 
560  } else {
561 
562  edm::LogWarning("EEIntegrityTask") << EcalElectronicsIdCollection3_ << " not available";
563 
564  }
565 
566  edm::Handle<EcalElectronicsIdCollection> ids7;
567 
568  if ( e.getByLabel(EcalElectronicsIdCollection4_, ids7) ) {
569 
570  for ( EcalElectronicsIdCollection::const_iterator idItr = ids7->begin(); idItr != ids7->end(); ++idItr ) {
571 
572  if ( Numbers::subDet( *idItr ) != EcalEndcap ) continue;
573 
574  int ism = Numbers::iSM( *idItr );
575 
576  int itt = idItr->towerId();
577  float iTt = itt + 0.5 - 69;
578 
579  if ( meIntegrityMemTTBlockSize[ism-1] ) meIntegrityMemTTBlockSize[ism-1]->Fill(iTt,0);
580 
581  }
582 
583  } else {
584 
585  edm::LogWarning("EEIntegrityTask") << EcalElectronicsIdCollection4_ << " not available";
586 
587  }
588 
589  edm::Handle<EcalElectronicsIdCollection> ids8;
590 
591  if ( e.getByLabel(EcalElectronicsIdCollection5_, ids8) ) {
592 
593  for ( EcalElectronicsIdCollection::const_iterator idItr = ids8->begin(); idItr != ids8->end(); ++idItr ) {
594 
595  if ( Numbers::subDet( *idItr ) != EcalEndcap ) continue;
596 
597  int ism = Numbers::iSM( *idItr );
598 
599  int chid = idItr->channelId();
600  int ie = EEIntegrityTask::chMemAbscissa[chid-1];
601  int ip = EEIntegrityTask::chMemOrdinate[chid-1];
602 
603  int itt = idItr->towerId();
604  ie += (itt-69)*5;
605 
606  float xix = ie - 0.5;
607  float xiy = ip - 0.5;
608 
609  if ( meIntegrityMemChId[ism-1] ) meIntegrityMemChId[ism-1]->Fill(xix,xiy);
610 
611  }
612 
613  } else {
614 
615  edm::LogWarning("EEIntegrityTask") << EcalElectronicsIdCollection5_ << " not available";
616 
617  }
618 
619  edm::Handle<EcalElectronicsIdCollection> ids9;
620 
621  if ( e.getByLabel(EcalElectronicsIdCollection6_, ids9) ) {
622 
623  for ( EcalElectronicsIdCollection::const_iterator idItr = ids9->begin(); idItr != ids9->end(); ++idItr ) {
624 
625  if ( Numbers::subDet( *idItr ) != EcalEndcap ) continue;
626 
627  int ism = Numbers::iSM( *idItr );
628 
629  int chid = idItr->channelId();
630  int ie = EEIntegrityTask::chMemAbscissa[chid-1];
631  int ip = EEIntegrityTask::chMemOrdinate[chid-1];
632 
633  int itt = idItr->towerId();
634  ie += (itt-69)*5;
635 
636  float xix = ie - 0.5;
637  float xiy = ip - 0.5;
638 
639  if ( meIntegrityMemGain[ism-1] ) meIntegrityMemGain[ism-1]->Fill(xix,xiy);
640 
641  }
642 
643  } else {
644 
645  edm::LogWarning("EEIntegrityTask") << EcalElectronicsIdCollection6_ << " not available";
646 
647  }
648 
649 }// end analyze
650 
651 const int EEIntegrityTask::chMemAbscissa [25] = {
652  1, 1, 1, 1, 1,
653  2, 2, 2, 2, 2,
654  3, 3, 3, 3, 3,
655  4, 4, 4, 4, 4,
656  5, 5, 5, 5, 5
657 };
658 
659 const int EEIntegrityTask::chMemOrdinate [25] = {
660  1, 2, 3, 4, 5,
661  5, 4, 3, 2, 1,
662  1, 2, 3, 4, 5,
663  5, 4, 3, 2, 1,
664  1, 2, 3, 4, 5
665 };
666 
EEIntegrityTask::EEDetIdCollection3_
edm::InputTag EEDetIdCollection3_
Definition: EEIntegrityTask.h:80
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
MonitorElement::getName
const std::string & getName(void) const
get name of ME
Definition: MonitorElement.h:99
Numbers::sEE
static std::string sEE(const unsigned ism)
Definition: Numbers.cc:226
i
int i
Definition: DBlmapReader.cc:9
EEIntegrityTask::EEDetIdCollection2_
edm::InputTag EEDetIdCollection2_
Definition: EEIntegrityTask.h:79
EEIntegrityTask::meIntegrityTTId
MonitorElement * meIntegrityTTId[18]
Definition: EEIntegrityTask.h:92
EEDetId::ix
int ix() const
Definition: EEDetId.h:71
EEIntegrityTask::EEDetIdCollection0_
edm::InputTag EEDetIdCollection0_
Definition: EEIntegrityTask.h:77
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
Numbers.h
Some &quot;id&quot; conversions.
EEIntegrityTask::meIntegrityErrorsByLumi
MonitorElement * meIntegrityErrorsByLumi
Definition: EEIntegrityTask.h:99
EEIntegrityTask::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Analyze.
Definition: EEIntegrityTask.cc:348
MonitorElement::setBinLabel
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:889
EEIntegrityTask::meIntegrityTTBlockSize
MonitorElement * meIntegrityTTBlockSize[18]
Definition: EEIntegrityTask.h:93
EEIntegrityTask::cleanup
void cleanup(void)
Cleanup.
Definition: EEIntegrityTask.cc:260
EEIntegrityTask::meIntegrityDCCSize
MonitorElement * meIntegrityDCCSize
Definition: EEIntegrityTask.h:98
EEIntegrityTask::beginLuminosityBlock
void beginLuminosityBlock(const edm::LuminosityBlock &lumiBlock, const edm::EventSetup &iSetup)
BeginLuminosityBlock.
Definition: EEIntegrityTask.cc:87
Numbers::ix0EE
static int ix0EE(const unsigned ism)
Definition: Numbers.cc:773
Numbers::iy0EE
static int iy0EE(const unsigned ism)
Definition: Numbers.cc:812
EEIntegrityTask::EcalElectronicsIdCollection1_
edm::InputTag EcalElectronicsIdCollection1_
Definition: EEIntegrityTask.h:82
EEIntegrityTask::EcalElectronicsIdCollection3_
edm::InputTag EcalElectronicsIdCollection3_
Definition: EEIntegrityTask.h:84
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:152
EEIntegrityTask::setup
void setup(void)
Setup.
Definition: EEIntegrityTask.cc:126
EEIntegrityTask::EcalElectronicsIdCollection5_
edm::InputTag EcalElectronicsIdCollection5_
Definition: EEIntegrityTask.h:86
DQMStore::tag
void tag(MonitorElement *me, unsigned int myTag)
Definition: DQMStore.cc:1354
EEIntegrityTask::prefixME_
std::string prefixME_
Definition: EEIntegrityTask.h:69
Numbers::crystals
static std::vector< DetId > * crystals(const EcalTrigTowerDetId &id)
Definition: Numbers.cc:585
EEIntegrityTask::meIntegrityGainSwitch
MonitorElement * meIntegrityGainSwitch[18]
Definition: EEIntegrityTask.h:91
EEIntegrityTask::beginJob
void beginJob(void)
BeginJob.
Definition: EEIntegrityTask.cc:74
EEIntegrityTask::meIntegrityMemChId
MonitorElement * meIntegrityMemChId[18]
Definition: EEIntegrityTask.h:94
EEIntegrityTask::EEDetIdCollection1_
edm::InputTag EEDetIdCollection1_
Definition: EEIntegrityTask.h:78
EEIntegrityTask::endRun
void endRun(const edm::Run &r, const edm::EventSetup &c)
EndRun.
Definition: EEIntegrityTask.cc:104
DQMStore::removeElement
void removeElement(const std::string &name)
Definition: DQMStore.cc:2572
EEIntegrityTask::reset
void reset(void)
Reset.
Definition: EEIntegrityTask.cc:108
edm::getName
std::string getName(Reflex::Type &cc)
Definition: ClassFiller.cc:18
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:356
EEIntegrityTask::chMemAbscissa
static const int chMemAbscissa[25]
Definition: EEIntegrityTask.h:103
EEIntegrityTask::meIntegrityGain
MonitorElement * meIntegrityGain[18]
Definition: EEIntegrityTask.h:90
EEIntegrityTask::EcalElectronicsIdCollection4_
edm::InputTag EcalElectronicsIdCollection4_
Definition: EEIntegrityTask.h:85
EEIntegrityTask::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
BeginRun.
Definition: EEIntegrityTask.cc:96
EEIntegrityTask::meIntegrityMemTTId
MonitorElement * meIntegrityMemTTId[18]
Definition: EEIntegrityTask.h:96
Numbers::initGeometry
static void initGeometry(const edm::EventSetup &setup, bool verbose=false)
Definition: Numbers.cc:50
EEIntegrityTask::chMemOrdinate
static const int chMemOrdinate[25]
Definition: EEIntegrityTask.h:104
EEIntegrityTask::dqmStore_
DQMStore * dqmStore_
Definition: EEIntegrityTask.h:67
EEIntegrityTask::EcalElectronicsIdCollection6_
edm::InputTag EcalElectronicsIdCollection6_
Definition: EEIntegrityTask.h:87
EEIntegrityTask::EEIntegrityTask
EEIntegrityTask(const edm::ParameterSet &ps)
Constructor.
Definition: EEIntegrityTask.cc:29
EEIntegrityTask::endJob
void endJob(void)
EndJob.
Definition: EEIntegrityTask.cc:340
Numbers::iSM
static unsigned iSM(const unsigned ism, const EcalSubdetector subdet)
Definition: Numbers.cc:246
EEIntegrityTask::subfolder_
std::string subfolder_
Definition: EEIntegrityTask.h:71
EEIntegrityTask::endLuminosityBlock
void endLuminosityBlock(const edm::LuminosityBlock &lumiBlock, const edm::EventSetup &iSetup)
EndLuminosityBlock.
Definition: EEIntegrityTask.cc:93
EEIntegrityTask::meIntegrityMemGain
MonitorElement * meIntegrityMemGain[18]
Definition: EEIntegrityTask.h:95
dir
dbl *** dir
Definition: mlp_gen.cc:35
Numbers::subDet
static EcalSubdetector subDet(const EBDetId &id)
Definition: Numbers.cc:145
MonitorElement::setLumiFlag
void setLumiFlag(void)
this ME is meant to be stored for each luminosity section
Definition: MonitorElement.h:136
EEIntegrityTask::meIntegrityChId
MonitorElement * meIntegrityChId[18]
Definition: EEIntegrityTask.h:89
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:845
ecalpyutils::ism
int ism(int ieta, int iphi)
Definition: EcalPyUtils.cc:56
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:916
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:422
EEIntegrityTask::~EEIntegrityTask
virtual ~EEIntegrityTask()
Destructor.
Definition: EEIntegrityTask.cc:70
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:429
EEIntegrityTask::EcalElectronicsIdCollection2_
edm::InputTag EcalElectronicsIdCollection2_
Definition: EEIntegrityTask.h:83
edm::Run
Definition: Run.h:33
EEIntegrityTask::meIntegrityMemTTBlockSize
MonitorElement * meIntegrityMemTTBlockSize[18]
Definition: EEIntegrityTask.h:97
edm::EDCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: EDCollection.h:20
Generated for CMSSW Reference Manual by doxygen 1.8.5