CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EETriggerTowerTask.cc
Go to the documentation of this file.
1 /*
2  * \file EETriggerTowerTask.cc
3  *
4  * $Date: 2011/06/08 18:15:04 $
5  * $Revision: 1.75.8.1 $
6  * \author G. Della Ricca
7  * \author E. Di Marco
8  *
9 */
10 
11 #include <vector>
12 
13 #include "FWCore/ServiceRegistry/interface/Service.h"
14 #include "FWCore/MessageLogger/interface/MessageLogger.h"
15 
16 #include "DataFormats/Common/interface/TriggerResults.h"
17 #include "DataFormats/EcalDetId/interface/EEDetId.h"
18 #include "DQMServices/Core/interface/MonitorElement.h"
19 #include "DQMServices/Core/interface/DQMStore.h"
20 
21 #include "DQM/EcalCommon/interface/Numbers.h"
22 
23 #include "DQM/EcalEndcapMonitorTasks/interface/EETriggerTowerTask.h"
24 #include "FWCore/Common/interface/TriggerNames.h"
25 
26 #include <iostream>
27 
28 const int EETriggerTowerTask::nTTEta = 20;
29 const int EETriggerTowerTask::nTTPhi = 20;
30 const int EETriggerTowerTask::nSM = 18;
31 
32 EETriggerTowerTask::EETriggerTowerTask(const edm::ParameterSet& ps) {
33 
34  init_ = false;
35 
36  dqmStore_ = edm::Service<DQMStore>().operator->();
37 
38  prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
39 
40  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
41 
42  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
43 
44  meEtSpectrumReal_[0] = 0;
45  meEtSpectrumReal_[1] = 0;
46  meEtSpectrumEmul_[0] = 0;
47  meEtSpectrumEmul_[1] = 0;
48  meEtSpectrumEmulMax_[0] = 0;
49  meEtSpectrumEmulMax_[1] = 0;
50  meEtBxReal_[0] = 0;
51  meEtBxReal_[1] = 0;
52  meOccupancyBxReal_[0] = 0;
53  meOccupancyBxReal_[1] = 0;
54  meTCCTimingCalo_[0] = 0;
55  meTCCTimingCalo_[1] = 0;
56  meTCCTimingMuon_[0] = 0;
57  meTCCTimingMuon_[1] = 0;
58  meEmulMatchIndex1D_[0] = 0;
59  meEmulMatchIndex1D_[1] = 0;
60  meEmulMatchMaxIndex1D_[0] = 0;
61  meEmulMatchMaxIndex1D_[1] = 0;
62 
63  reserveArray(meEtMapReal_);
64  reserveArray(meVetoReal_);
65  reserveArray(meEtMapEmul_);
66  reserveArray(meVetoEmul_);
67  reserveArray(meEmulError_);
68  reserveArray(meEmulMatch_);
69  reserveArray(meVetoEmulError_);
70 
71  realCollection_ = ps.getParameter<edm::InputTag>("EcalTrigPrimDigiCollectionReal");
72  emulCollection_ = ps.getParameter<edm::InputTag>("EcalTrigPrimDigiCollectionEmul");
73  EEDigiCollection_ = ps.getParameter<edm::InputTag>("EEDigiCollection");
74  HLTResultsCollection_ = ps.getParameter<edm::InputTag>("HLTResultsCollection");
75 
76  HLTCaloHLTBit_ = ps.getUntrackedParameter<std::string>("HLTCaloHLTBit", "");
77  HLTMuonHLTBit_ = ps.getUntrackedParameter<std::string>("HLTMuonHLTBit", "");
78 
79  outputFile_ = ps.getUntrackedParameter<std::string>("OutputRootFile", "");
80 
81  LogDebug("EETriggerTowerTask") << "REAL digis: " << realCollection_;
82  LogDebug("EETriggerTowerTask") << "EMULATED digis: " << emulCollection_;
83 
84 }
85 
86 EETriggerTowerTask::~EETriggerTowerTask(){
87 
88 }
89 
90 void EETriggerTowerTask::reserveArray( array1& array ) {
91 
92  array.reserve( nSM );
93  array.resize( nSM, static_cast<MonitorElement*>(0) );
94 
95 }
96 
97 void EETriggerTowerTask::beginJob(void){
98 
99  ievt_ = 0;
100 
101  if ( dqmStore_ ) {
102  dqmStore_->setCurrentFolder(prefixME_ + "/EETriggerTowerTask");
103  dqmStore_->rmdir(prefixME_ + "/EETriggerTowerTask");
104  }
105 
106 }
107 
108 void EETriggerTowerTask::beginRun(const edm::Run& r, const edm::EventSetup& c) {
109 
110  Numbers::initGeometry(c, false);
111 
112  if ( ! mergeRuns_ ) this->reset();
113 
114 }
115 
116 void EETriggerTowerTask::endRun(const edm::Run& r, const edm::EventSetup& c) {
117 
118 }
119 
120 void EETriggerTowerTask::reset(void) {
121 
122  for (int iside = 0; iside < 2; iside++) {
123  if ( meEtSpectrumReal_[iside] ) meEtSpectrumReal_[iside]->Reset();
124  if ( meEtSpectrumEmul_[iside] ) meEtSpectrumEmul_[iside]->Reset();
125  if ( meEtSpectrumEmulMax_[iside] ) meEtSpectrumEmulMax_[iside]->Reset();
126  if ( meEtBxReal_[iside] ) meEtBxReal_[iside]->Reset();
127  if ( meOccupancyBxReal_[iside] ) meOccupancyBxReal_[iside]->Reset();
128  if ( meTCCTimingCalo_[iside] ) meTCCTimingCalo_[iside]->Reset();
129  if ( meTCCTimingMuon_[iside] ) meTCCTimingMuon_[iside]->Reset();
130  if ( meEmulMatchIndex1D_[iside] ) meEmulMatchIndex1D_[iside]->Reset();
131  if ( meEmulMatchMaxIndex1D_[iside] ) meEmulMatchMaxIndex1D_[iside]->Reset();
132  }
133 
134  for (int i = 0; i < 18; i++) {
135 
136  if ( meEtMapReal_[i] ) meEtMapReal_[i]->Reset();
137  if ( meVetoReal_[i] ) meVetoReal_[i]->Reset();
138  if ( meEtMapEmul_[i] ) meEtMapEmul_[i]->Reset();
139  if ( meVetoEmul_[i] ) meVetoEmul_[i]->Reset();
140  if ( meEmulError_[i] ) meEmulError_[i]->Reset();
141  if ( meEmulMatch_[i] ) meEmulMatch_[i]->Reset();
142  if ( meVetoEmulError_[i] ) meVetoEmulError_[i]->Reset();
143 
144  }
145 
146 }
147 
148 void EETriggerTowerTask::setup(void){
149 
150  init_ = true;
151 
152  if ( dqmStore_ ) {
153  setup( "Real Digis",
154  (prefixME_ + "/EETriggerTowerTask").c_str(), false );
155 
156  setup( "Emulated Digis",
157  (prefixME_ + "/EETriggerTowerTask/Emulated").c_str(), true);
158  }
159  else {
160  edm::LogError("EETriggerTowerTask") << "Bad DQMStore, cannot book MonitorElements.";
161  }
162 }
163 
164 void EETriggerTowerTask::setup( std::string const &nameext,
165  std::string const &folder,
166  bool emulated ) {
167 
168  array1* meEtMap = &meEtMapReal_;
169  array1* meVeto = &meVetoReal_;
170 
171  if ( emulated ) {
172  meEtMap = &meEtMapEmul_;
173  meVeto = &meVetoEmul_;
174  }
175 
176  dqmStore_->setCurrentFolder(folder);
177 
178  std::string name;
179 
180  if (!emulated) {
181  name = "EETTT Et spectrum " + nameext + " EE -";
182  meEtSpectrumReal_[0] = dqmStore_->book1D(name, name, 256, 0., 256.);
183  meEtSpectrumReal_[0]->setAxisTitle("energy (ADC)", 1);
184 
185  name = "EETTT Et spectrum " + nameext + " EE +";
186  meEtSpectrumReal_[1] = dqmStore_->book1D(name, name, 256, 0., 256.);
187  meEtSpectrumReal_[1]->setAxisTitle("energy (ADC)", 1);
188 
189  name = "EETTT TP matching index EE -";
190  meEmulMatchIndex1D_[0] = dqmStore_->book1D(name, name, 7, -1., 6.);
191  meEmulMatchIndex1D_[0]->setAxisTitle("TP data matching emulator", 1);
192 
193  name = "EETTT TP matching index EE +";
194  meEmulMatchIndex1D_[1] = dqmStore_->book1D(name, name, 7, -1., 6.);
195  meEmulMatchIndex1D_[1]->setAxisTitle("TP data matching emulator", 1);
196 
197  name = "EETTT max TP matching index EE -";
198  meEmulMatchMaxIndex1D_[0] = dqmStore_->book1D(name, name, 7, -1., 6.);
199  meEmulMatchMaxIndex1D_[0]->setAxisTitle("Max TP data matching emulator", 1);
200 
201  name = "EETTT max TP matching index EE +";
202  meEmulMatchMaxIndex1D_[1] = dqmStore_->book1D(name, name, 7, -1., 6.);
203  meEmulMatchMaxIndex1D_[1]->setAxisTitle("Max TP data matching emulator", 1);
204 
205  double xbins[51];
206  for ( int i=0; i<=11; i++ ) xbins[i] = i-1; // begin of orbit
207  // abort gap in presence of calibration: [3381-3500]
208  // abort gap in absence of calibration: [3444-3500]
209  // using the wider abort gap always, start finer binning at bx=3371
210  for ( int i=12; i<=22; i++) xbins[i] = 3371+i-12;
211  // use 29 bins for the abort gap
212  for ( int i=23; i<=50; i++) xbins[i] = 3382+(i-23)*6;
213 
214  name = "EETTT Et vs bx " + nameext + " EE -";
215  meEtBxReal_[0] = dqmStore_->bookProfile(name, name, 50, xbins, 256, 0, 256);
216  meEtBxReal_[0]->setAxisTitle("bunch crossing", 1);
217  meEtBxReal_[0]->setAxisTitle("energy (ADC)", 2);
218 
219  name = "EETTT Et vs bx " + nameext + " EE +";
220  meEtBxReal_[1] = dqmStore_->bookProfile(name, name, 50, xbins, 256, 0, 256);
221  meEtBxReal_[1]->setAxisTitle("bunch crossing", 1);
222  meEtBxReal_[1]->setAxisTitle("energy (ADC)", 2);
223 
224  name = "EETTT TP occupancy vs bx " + nameext + " EE -";
225  meOccupancyBxReal_[0] = dqmStore_->bookProfile(name, name, 50, xbins, 2448, 0, 2448);
226  meOccupancyBxReal_[0]->setAxisTitle("bunch crossing", 1);
227  meOccupancyBxReal_[0]->setAxisTitle("TP number", 2);
228 
229  name = "EETTT TP occupancy vs bx " + nameext + " EE +";
230  meOccupancyBxReal_[1] = dqmStore_->bookProfile(name, name, 50, xbins, 2448, 0, 2448);
231  meOccupancyBxReal_[1]->setAxisTitle("bunch crossing", 1);
232  meOccupancyBxReal_[1]->setAxisTitle("TP number", 2);
233 
234  if ( HLTCaloHLTBit_ != "" ) {
235  name = "EETTT TCC timing calo triggers " + nameext + " EE -";
236  meTCCTimingCalo_[0] = dqmStore_->book2D(name, name, 36, 1, 37, 7, -1., 6.);
237  meTCCTimingCalo_[0]->setAxisTitle("nTCC", 1);
238  meTCCTimingCalo_[0]->setAxisTitle("TP data matching emulator", 2);
239 
240  name = "EETTT TCC timing calo triggers " + nameext + " EE +";
241  meTCCTimingCalo_[1] = dqmStore_->book2D(name, name, 36, 73, 109, 7, -1., 6.);
242  meTCCTimingCalo_[1]->setAxisTitle("nTCC", 1);
243  meTCCTimingCalo_[1]->setAxisTitle("TP data matching emulator", 2);
244  }
245 
246  if ( HLTMuonHLTBit_ != "" ) {
247  name = "EETTT TCC timing muon triggers " + nameext + " EE -";
248  meTCCTimingMuon_[0] = dqmStore_->book2D(name, name, 36, 1, 37, 7, -1., 6.);
249  meTCCTimingMuon_[0]->setAxisTitle("nTCC", 1);
250  meTCCTimingMuon_[0]->setAxisTitle("TP data matching emulator", 2);
251 
252  name = "EETTT TCC timing muon triggers " + nameext + " EE +";
253  meTCCTimingMuon_[1] = dqmStore_->book2D(name, name, 36, 73, 109, 7, -1., 6.);
254  meTCCTimingMuon_[1]->setAxisTitle("nTCC", 1);
255  meTCCTimingMuon_[1]->setAxisTitle("TP data matching emulator", 2);
256  }
257 
258  } else {
259  name = "EETTT Et spectrum " + nameext + " EE -";
260  meEtSpectrumEmul_[0] = dqmStore_->book1D(name, name, 256, 0., 256.);
261  meEtSpectrumEmul_[0]->setAxisTitle("energy (ADC)", 1);
262 
263  name = "EETTT Et spectrum " + nameext + " EE +";
264  meEtSpectrumEmul_[1] = dqmStore_->book1D(name, name, 256, 0., 256.);
265  meEtSpectrumEmul_[1]->setAxisTitle("energy (ADC)", 1);
266 
267  name = "EETTT Et spectrum " + nameext + " max EE -";
268  meEtSpectrumEmulMax_[0] = dqmStore_->book1D(name, name, 256, 0., 256.);
269  meEtSpectrumEmulMax_[0]->setAxisTitle("energy (ADC)", 1);
270 
271  name = "EETTT Et spectrum " + nameext + " max EE +";
272  meEtSpectrumEmulMax_[1] = dqmStore_->book1D(name, name, 256, 0., 256.);
273  meEtSpectrumEmulMax_[1]->setAxisTitle("energy (ADC)", 1);
274  }
275 
276  for (int i = 0; i < 18; i++) {
277 
278  name = "EETTT Et map " + nameext + " " + Numbers::sEE(i+1);
279  (*meEtMap)[i] = dqmStore_->bookProfile2D(name, name,
280  50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50.,
281  50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50.,
282  256, 0, 256.);
283  (*meEtMap)[i]->setAxisTitle("ix", 1);
284  if ( i+1 >= 1 && i+1 <= 9 ) (*meEtMap)[i]->setAxisTitle("101-ix", 1);
285  (*meEtMap)[i]->setAxisTitle("iy", 2);
286  dqmStore_->tag((*meEtMap)[i], i+1);
287 
288  if (!emulated) {
289 
290  name = "EETTT EmulError " + Numbers::sEE(i+1);
291  meEmulError_[i] = dqmStore_->book2D(name, name,
292  50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50.,
293  50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50. );
294  meEmulError_[i]->setAxisTitle("ix", 1);
295  if ( i+1 >= 1 && i+1 <= 9 ) meEmulError_[i]->setAxisTitle("101-ix", 1);
296  meEmulError_[i]->setAxisTitle("iy", 2);
297  dqmStore_->tag(meEmulError_[i], i+1);
298 
299  name = "EETTT EmulMatch " + Numbers::sEE(i+1);
300  meEmulMatch_[i] = dqmStore_->book3D(name, name,
301  50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50.,
302  50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50.,
303  6, 0., 6.);
304  meEmulMatch_[i]->setAxisTitle("ix", 1);
305  if ( i+1 >= 1 && i+1 <= 9 ) meEmulMatch_[i]->setAxisTitle("101-ix", 1);
306  meEmulMatch_[i]->setAxisTitle("iy", 2);
307  dqmStore_->tag(meEmulMatch_[i], i+1);
308 
309  name = "EETTT EmulFineGrainVetoError " + Numbers::sEE(i+1);
310  meVetoEmulError_[i] = dqmStore_->book2D(name, name,
311  50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50.,
312  50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50.);
313  meVetoEmulError_[i]->setAxisTitle("ix", 1);
314  if ( i+1 >= 1 && i+1 <= 9 ) meVetoEmulError_[i]->setAxisTitle("101-ix", 1);
315  meVetoEmulError_[i]->setAxisTitle("iy", 2);
316  dqmStore_->tag(meVetoEmulError_[i], i+1);
317 
318  }
319  }
320 
321 }
322 
323 void EETriggerTowerTask::cleanup(void) {
324 
325  if ( ! init_ ) return;
326 
327  if ( dqmStore_ ) {
328 
329  if ( !outputFile_.empty() ) dqmStore_->save( outputFile_.c_str() );
330 
331  dqmStore_->rmdir( prefixME_ + "/EETriggerTowerTask" );
332 
333  }
334 
335  init_ = false;
336 
337 }
338 
339 void EETriggerTowerTask::endJob(void){
340 
341  edm::LogInfo("EETriggerTowerTask") << "analyzed " << ievt_ << " events";
342 
343  if ( enableCleanup_ ) this->cleanup();
344 
345 }
346 
347 void EETriggerTowerTask::analyze(const edm::Event& e, const edm::EventSetup& c){
348 
349  if ( ! init_ ) this->setup();
350 
351  ievt_++;
352 
353  edm::Handle<EcalTrigPrimDigiCollection> realDigis;
354 
355  if ( e.getByLabel(realCollection_, realDigis) ) {
356 
357  int neetpd = realDigis->size();
358  LogDebug("EETriggerTowerTask") << "event " << ievt_ << " trigger primitive digi collection size: " << neetpd;
359 
360  processDigis( e,
361  realDigis,
362  meEtMapReal_,
363  meVetoReal_);
364 
365  } else {
366  edm::LogWarning("EETriggerTowerTask") << realCollection_ << " not available";
367  }
368 
369  edm::Handle<EcalTrigPrimDigiCollection> emulDigis;
370 
371  if ( e.getByLabel(emulCollection_, emulDigis) ) {
372 
373  edm::Handle<edm::TriggerResults> hltResults;
374 
375  if ( !e.getByLabel(HLTResultsCollection_, hltResults) ) {
376  HLTResultsCollection_ = edm::InputTag(HLTResultsCollection_.label(), HLTResultsCollection_.instance(), "HLT");
377  }
378 
379  if ( !e.getByLabel(HLTResultsCollection_, hltResults) ) {
380  HLTResultsCollection_ = edm::InputTag(HLTResultsCollection_.label(), HLTResultsCollection_.instance(), "FU");
381  }
382 
383  if ( e.getByLabel(HLTResultsCollection_, hltResults) ) {
384 
385  processDigis( e,
386  emulDigis,
387  meEtMapEmul_,
388  meVetoEmul_,
389  realDigis,
390  hltResults);
391 
392  } else {
393  edm::LogWarning("EETriggerTowerTask") << HLTResultsCollection_ << " not available";
394  }
395 
396  } else {
397  edm::LogInfo("EETriggerTowerTask") << emulCollection_ << " not available";
398  }
399 
400 }
401 
402 void
403 EETriggerTowerTask::processDigis( const edm::Event& e, const edm::Handle<EcalTrigPrimDigiCollection>& digis,
404  array1& meEtMap,
405  array1& meVeto,
406  const edm::Handle<EcalTrigPrimDigiCollection>& compDigis,
407  const edm::Handle<edm::TriggerResults> & hltResults ) {
408 
409  int bx = e.bunchCrossing();
410  int nTP[2];
411  nTP[0] = nTP[1] = 0;
412 
413  // indexes are: readoutCrystalsInTower[TCCId][iTT]
414  int readoutCrystalsInTower[108][41];
415  for (int itcc = 0; itcc < 108; itcc++) {
416  for (int itt = 0; itt < 41; itt++) readoutCrystalsInTower[itcc][itt] = 0;
417  }
418 
419  if ( compDigis.isValid() ) {
420 
421  edm::Handle<EEDigiCollection> crystalDigis;
422 
423  if ( e.getByLabel(EEDigiCollection_, crystalDigis) ) {
424 
425  for ( EEDigiCollection::const_iterator cDigiItr = crystalDigis->begin(); cDigiItr != crystalDigis->end(); ++cDigiItr ) {
426 
427  EEDetId id = cDigiItr->id();
428 
429  int ix = id.ix();
430  int iy = id.iy();
431  int ism = Numbers::iSM( id );
432  int itcc = Numbers::iTCC( ism, EcalEndcap, ix, iy );
433  int itt = Numbers::iTT( ism, EcalEndcap, ix, iy );
434 
435  readoutCrystalsInTower[itcc-1][itt-1]++;
436 
437  }
438 
439  } else {
440  edm::LogWarning("EETriggerTowerTask") << EEDigiCollection_ << " not available";
441  }
442 
443  }
444 
445  bool caloTrg = false;
446  bool muonTrg = false;
447 
448  if ( hltResults.isValid() ) {
449 
450  int ntrigs = hltResults->size();
451  if ( ntrigs!=0 ) {
452 
453  const edm::TriggerNames & triggerNames = e.triggerNames(*hltResults);
454 
455  for ( int itrig = 0; itrig != ntrigs; ++itrig ) {
456  std::string trigName = triggerNames.triggerName(itrig);
457  bool accept = hltResults->accept(itrig);
458 
459  if ( trigName == HLTCaloHLTBit_ ) caloTrg = accept;
460 
461  if ( trigName == HLTMuonHLTBit_ ) muonTrg = accept;
462 
463  }
464 
465  } else {
466  edm::LogWarning("EBTriggerTowerTask") << " zero size trigger names in input TriggerResults";
467  }
468 
469  }
470 
471  for ( EcalTrigPrimDigiCollection::const_iterator tpdigiItr = digis->begin(); tpdigiItr != digis->end(); ++tpdigiItr ) {
472 
473  if ( Numbers::subDet( tpdigiItr->id() ) != EcalEndcap ) continue;
474 
475  int ismt = Numbers::iSM( tpdigiItr->id() );
476  int itt = Numbers::iTT( tpdigiItr->id() );
477  int itcc = Numbers::iTCC( tpdigiItr->id() );
478 
479  std::vector<DetId>* crystals = Numbers::crystals( tpdigiItr->id() );
480 
481  float xvalEt = tpdigiItr->compressedEt();
482  float xvalVeto = 0.5 + tpdigiItr->fineGrain();
483 
484  bool good = true;
485  bool goodVeto = true;
486 
487  int compDigiInterest = -1;
488 
489  bool matchSample[6];
490  for (int j=0; j<6; j++) matchSample[j]=false;
491 
492  if ( compDigis.isValid() ) {
493 
494  if ( ismt >= 1 && ismt <= 9 ) {
495  if ( meEtSpectrumEmul_[0] ) meEtSpectrumEmul_[0]->Fill( xvalEt );
496  } else {
497  if ( meEtSpectrumEmul_[1] ) meEtSpectrumEmul_[1]->Fill( xvalEt );
498  }
499 
500  float maxEt = 0;
501  int maxTPIndex = -1;
502  for (int j=0; j<5; j++) {
503  float EtTP = (*tpdigiItr)[j].compressedEt();
504  if ( EtTP > maxEt ) {
505  maxEt = EtTP;
506  maxTPIndex = j+1;
507  }
508  }
509 
510  if ( ismt >= 1 && ismt <= 9 ) {
511  if ( meEtSpectrumEmulMax_[0] ) meEtSpectrumEmulMax_[0]->Fill( maxEt );
512  if ( meEmulMatchMaxIndex1D_[0] && maxEt > 0 ) meEmulMatchMaxIndex1D_[0]->Fill( maxTPIndex );
513  } else {
514  if ( meEtSpectrumEmulMax_[1] ) meEtSpectrumEmulMax_[1]->Fill( maxEt );
515  if ( meEmulMatchMaxIndex1D_[1] && maxEt > 0 ) meEmulMatchMaxIndex1D_[1]->Fill( maxTPIndex );
516  }
517 
518  EcalTrigPrimDigiCollection::const_iterator compDigiItr = compDigis->find( tpdigiItr->id().rawId() );
519  if ( compDigiItr != compDigis->end() ) {
520  int compDigiEt = compDigiItr->compressedEt();
521  compDigiInterest = (compDigiItr->ttFlag() & 0x3);
522 
523  if ( ismt >= 1 && ismt <= 9 ) {
524  if ( compDigiEt > 0 ) nTP[0]++;
525  if ( meEtSpectrumReal_[0] ) meEtSpectrumReal_[0]->Fill( compDigiEt );
526  if ( meEtBxReal_[0] && compDigiEt > 0 ) meEtBxReal_[0]->Fill( bx, compDigiEt );
527  } else {
528  if ( compDigiEt > 0 ) nTP[1]++;
529  if ( meEtSpectrumReal_[1] ) meEtSpectrumReal_[1]->Fill( compDigiEt );
530  if ( meEtBxReal_[1] && compDigiEt > 0 ) meEtBxReal_[1]->Fill( bx, compDigiEt );
531  }
532 
533  // compare the 5 TPs with different time-windows
534  // sample 0 means no match, 1-5: sample of the TP that matches
535  matchSample[0]=false;
536  bool matchedAny=false;
537 
538  for (int j=0; j<5; j++) {
539  if ((*tpdigiItr)[j].compressedEt() == compDigiEt ) {
540  matchSample[j+1]=true;
541  matchedAny=true;
542  } else {
543  matchSample[j+1]=false;
544  }
545  }
546 
547  if (!matchedAny) matchSample[0]=true;
548 
549  // check if the tower has been readout completely and if it is medium or high interest
550  if (readoutCrystalsInTower[itcc-1][itt-1] == int(crystals->size()) &&
551  (compDigiInterest == 1 || compDigiInterest == 3) && compDigiEt > 0) {
552 
553  if ( tpdigiItr->compressedEt() != compDigiEt ) {
554  good = false;
555  }
556  if ( tpdigiItr->fineGrain() != compDigiItr->fineGrain() ) {
557  goodVeto = false;
558  }
559 
560  for (int j=0; j<6; j++) {
561  if (matchSample[j]) {
562 
563  int index = ( j==0 ) ? -1 : j;
564 
565  if ( ismt >= 1 && ismt <= 9 ) {
566  meEmulMatchIndex1D_[0]->Fill(index+0.5);
567  } else {
568  meEmulMatchIndex1D_[1]->Fill(index+0.5);
569  }
570 
571  for ( unsigned int i=0; i<crystals->size(); i++ ) {
572 
573  EEDetId id = (*crystals)[i];
574 
575  int ix = id.ix();
576  int iy = id.iy();
577 
578  if ( ismt >= 1 && ismt <= 9 ) ix = 101 - ix;
579 
580  float xix = ix-0.5;
581  float xiy = iy-0.5;
582 
583  meEmulMatch_[ismt-1]->Fill(xix, xiy, j+0.5);
584  if ( ismt >= 1 && ismt <= 9 ) {
585  if ( meTCCTimingCalo_[0] && caloTrg ) meTCCTimingCalo_[0]->Fill( itcc, index+0.5 );
586  if ( meTCCTimingMuon_[0] && muonTrg ) meTCCTimingMuon_[0]->Fill( itcc, index+0.5 );
587  } else {
588  if ( meTCCTimingCalo_[1] && caloTrg ) meTCCTimingCalo_[1]->Fill( itcc, index+0.5 );
589  if ( meTCCTimingMuon_[1] && muonTrg ) meTCCTimingMuon_[1]->Fill( itcc, index+0.5 );
590  }
591 
592  } // loop on crystals
593 
594  }
595  }
596 
597  } // check readout
598 
599  } else {
600  good = false;
601  goodVeto = false;
602  }
603 
604  for ( unsigned int i=0; i<crystals->size(); i++ ) {
605 
606  EEDetId id = (*crystals)[i];
607 
608  int ix = id.ix();
609  int iy = id.iy();
610 
611  if ( ismt >= 1 && ismt <= 9 ) ix = 101 - ix;
612 
613  float xix = ix-0.5;
614  float xiy = iy-0.5;
615 
616  if (!good ) {
617  if ( meEmulError_[ismt-1] ) meEmulError_[ismt-1]->Fill(xix, xiy);
618  }
619  if (!goodVeto) {
620  if ( meVetoEmulError_[ismt-1] ) meVetoEmulError_[ismt-1]->Fill(xix, xiy);
621  }
622 
623  } // loop on crystals
624 
625  } // compDigis.isValid
626 
627  for ( unsigned int i=0; i<crystals->size(); i++ ) {
628 
629  EEDetId id = (*crystals)[i];
630 
631  int ix = id.ix();
632  int iy = id.iy();
633 
634  if ( ismt >= 1 && ismt <= 9 ) ix = 101 - ix;
635 
636  float xix = ix-0.5;
637  float xiy = iy-0.5;
638 
639  if ( meEtMap[ismt-1] ) meEtMap[ismt-1]->Fill(xix, xiy, xvalEt);
640  if ( meVeto[ismt-1] ) meVeto[ismt-1]->Fill(xix, xiy, xvalVeto);
641 
642  } // loop on crystals
643 
644  } // loop on TP
645 
646  if ( meOccupancyBxReal_[0] ) meOccupancyBxReal_[0]->Fill( bx, nTP[0] );
647  if ( meOccupancyBxReal_[1] ) meOccupancyBxReal_[1]->Fill( bx, nTP[1] );
648 
649 }
650 
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:516
EETriggerTowerTask::enableCleanup_
bool enableCleanup_
remove MEs
Definition: EETriggerTowerTask.h:126
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
EETriggerTowerTask::endJob
void endJob(void)
EndJob.
Definition: EETriggerTowerTask.cc:339
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
edm::Event::triggerNames
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
Definition: Event.cc:208
EETriggerTowerTask::~EETriggerTowerTask
virtual ~EETriggerTowerTask()
Destructor.
Definition: EETriggerTowerTask.cc:86
EEDetId::ix
int ix() const
Definition: EEDetId.h:71
EETriggerTowerTask::meVetoEmul_
array1 meVetoEmul_
Emulated fine grain veto vs ix vs iy, for each SM.
Definition: EETriggerTowerTask.h:107
edm::DataFrameContainer::const_iterator
boost::transform_iterator< IterHelp, boost::counting_iterator< int > > const_iterator
Definition: DataFrameContainer.h:51
fw3dlego::xbins
const double xbins[]
Definition: fw3dlego_xbins.cc:18
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:514
EETriggerTowerTask::meEtBxReal_
MonitorElement * meEtBxReal_[2]
number and ET average of TP vs bx for the whole EE+/-
Definition: EETriggerTowerTask.h:158
DQMStore::rmdir
void rmdir(const std::string &fullpath)
Definition: DQMStore.cc:2311
Numbers::iy0EE
static int iy0EE(const int ism)
Definition: Numbers.cc:987
DQMStore::book3D
MonitorElement * book3D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, int nchZ, double lowZ, double highZ)
Book 3D histogram.
Definition: DQMStore.cc:776
Numbers.h
Some &quot;id&quot; conversions.
EETriggerTowerTask::outputFile_
std::string outputFile_
debug output root file. if empty, no output file created.
Definition: EETriggerTowerTask.h:146
Numbers::sEE
static std::string sEE(const int ism)
Definition: Numbers.cc:205
DQMStore::save
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:1898
EETriggerTowerTask::nSM
static const int nSM
number of supermodules
Definition: EETriggerTowerTask.h:40
EETriggerTowerTask::setup
void setup(void)
Setup.
Definition: EETriggerTowerTask.cc:148
EETriggerTowerTask::nTTPhi
static const int nTTPhi
number of trigger towers in phi
Definition: EETriggerTowerTask.h:37
EETriggerTowerTask::meOccupancyBxReal_
MonitorElement * meOccupancyBxReal_[2]
Definition: EETriggerTowerTask.h:159
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:81
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:62
EETriggerTowerTask::reserveArray
void reserveArray(array1 &array)
reserve an array to hold one histogram per supermodule
Definition: EETriggerTowerTask.cc:90
EETriggerTowerTask::meEtSpectrumEmul_
MonitorElement * meEtSpectrumEmul_[2]
Definition: EETriggerTowerTask.h:154
EETriggerTowerTask::emulCollection_
edm::InputTag emulCollection_
to find the input collection of emulated digis
Definition: EETriggerTowerTask.h:135
EETriggerTowerTask::beginJob
void beginJob(void)
BeginJob.
Definition: EETriggerTowerTask.cc:97
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:21
Numbers::iTCC
static int iTCC(const int ism, const EcalSubdetector subdet, const int i1, const int i2)
Definition: Numbers.cc:621
EETriggerTowerTask::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
BeginRun.
Definition: EETriggerTowerTask.cc:108
Numbers::iTT
static int iTT(const int ism, const EcalSubdetector subdet, const int i1, const int i2)
Definition: Numbers.cc:526
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:147
DQMStore::tag
void tag(MonitorElement *me, unsigned int myTag)
Definition: DQMStore.cc:1151
Numbers::crystals
static std::vector< DetId > * crystals(const EcalTrigTowerDetId &id)
Definition: Numbers.cc:740
EETriggerTowerTask::HLTResultsCollection_
edm::InputTag HLTResultsCollection_
to find the input collection of HLT bits
Definition: EETriggerTowerTask.h:141
EETriggerTowerTask::array1
std::vector< MonitorElement * > array1
1D array
Definition: EETriggerTowerTask.h:72
EETriggerTowerTask::processDigis
void processDigis(const edm::Event &e, const edm::Handle< EcalTrigPrimDigiCollection > &digis, array1 &meEtMap, array1 &meVeto, const edm::Handle< EcalTrigPrimDigiCollection > &compDigis=edm::Handle< EcalTrigPrimDigiCollection >(), const edm::Handle< edm::TriggerResults > &hltResults=edm::Handle< edm::TriggerResults >())
process a collection of digis, either real or emulated
Definition: EETriggerTowerTask.cc:403
Numbers::iSM
static int iSM(const int ism, const EcalSubdetector subdet)
Definition: Numbers.cc:222
EETriggerTowerTask::nTTEta
static const int nTTEta
number of trigger towers in eta
Definition: EETriggerTowerTask.h:34
EETriggerTowerTask::dqmStore_
DQMStore * dqmStore_
DQM back-end interface.
Definition: EETriggerTowerTask.h:120
EETriggerTowerTask::meEmulMatchMaxIndex1D_
MonitorElement * meEmulMatchMaxIndex1D_[2]
Definition: EETriggerTowerTask.h:150
EETriggerTowerTask::meEtMapEmul_
array1 meEtMapEmul_
Emulated Et vs ix vs iy, for each SM.
Definition: EETriggerTowerTask.h:104
JetPDSkim_cfg.hltResults
tuple hltResults
Definition: JetPDSkim_cfg.py:151
EETriggerTowerTask::meVetoReal_
array1 meVetoReal_
fine grain veto vs ix vs iy, for each SM
Definition: EETriggerTowerTask.h:101
j
int j
Definition: DBlmapReader.cc:9
DQMStore::bookProfile
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
Definition: DQMStore.cc:828
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:355
EETriggerTowerTask::meEtMapReal_
array1 meEtMapReal_
Et vs ix vs iy, for each SM.
Definition: EETriggerTowerTask.h:98
EETriggerTowerTask::prefixME_
std::string prefixME_
path to MEs
Definition: EETriggerTowerTask.h:123
EETriggerTowerTask::EETriggerTowerTask
EETriggerTowerTask(const edm::ParameterSet &ps)
Constructor.
Definition: EETriggerTowerTask.cc:32
EETriggerTowerTask::meEtSpectrumEmulMax_
MonitorElement * meEtSpectrumEmulMax_[2]
Definition: EETriggerTowerTask.h:155
EETriggerTowerTask::mergeRuns_
bool mergeRuns_
merge MEs across runs
Definition: EETriggerTowerTask.h:129
EETriggerTowerTask::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Analyze.
Definition: EETriggerTowerTask.cc:347
Numbers::initGeometry
static void initGeometry(const edm::EventSetup &setup, bool verbose=false)
Definition: Numbers.cc:45
edm::TriggerNames::triggerName
std::string const & triggerName(unsigned int index) const
Definition: TriggerNames.cc:27
EETriggerTowerTask::endRun
void endRun(const edm::Run &r, const edm::EventSetup &c)
EndRun.
Definition: EETriggerTowerTask.cc:116
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:25
EETriggerTowerTask::reset
void reset(void)
Reset.
Definition: EETriggerTowerTask.cc:120
Numbers::ix0EE
static int ix0EE(const int ism)
Definition: Numbers.cc:952
EETriggerTowerTask::meEtSpectrumReal_
MonitorElement * meEtSpectrumReal_[2]
ET spectrums for the whole EE+/-.
Definition: EETriggerTowerTask.h:153
EETriggerTowerTask::realCollection_
edm::InputTag realCollection_
to find the input collection of real digis
Definition: EETriggerTowerTask.h:132
Numbers::subDet
static EcalSubdetector subDet(const EBDetId &id)
Definition: Numbers.cc: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
EETriggerTowerTask::meTCCTimingCalo_
MonitorElement * meTCCTimingCalo_[2]
TCC timing.
Definition: EETriggerTowerTask.h:162
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
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:26
EETriggerTowerTask::meEmulMatchIndex1D_
MonitorElement * meEmulMatchIndex1D_[2]
1D emulator match 1D
Definition: EETriggerTowerTask.h:149
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:232
EETriggerTowerTask::cleanup
void cleanup(void)
Cleanup.
Definition: EETriggerTowerTask.cc:323
EETriggerTowerTask::EEDigiCollection_
edm::InputTag EEDigiCollection_
to find the input collection of crystal digis
Definition: EETriggerTowerTask.h:138
edm::Run
Definition: Run.h:32
EETriggerTowerTask::meTCCTimingMuon_
MonitorElement * meTCCTimingMuon_[2]
Definition: EETriggerTowerTask.h:163
EETriggerTowerTask::ievt_
int ievt_
local event counter
Definition: EETriggerTowerTask.h:95
DQMStore::bookProfile2D
MonitorElement * bookProfile2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, int nchZ, double lowZ, double highZ, const char *option="s")
Definition: DQMStore.cc:972
Generated for CMSSW Reference Manual by doxygen 1.8.5