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