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