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