CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EBTestPulseTask.cc
Go to the documentation of this file.
1 /*
2  * \file EBTestPulseTask.cc
3  *
4  * $Date: 2011/08/23 00:30:59 $
5  * $Revision: 1.116.4.2 $
6  * \author G. Della Ricca
7  * \author G. Franzoni
8  *
9 */
10 
11 #include <iostream>
12 #include <sstream>
13 #include <iomanip>
14 #include <vector>
15 
16 #include "FWCore/ServiceRegistry/interface/Service.h"
17 #include "FWCore/MessageLogger/interface/MessageLogger.h"
18 
19 #include "DQMServices/Core/interface/MonitorElement.h"
20 
21 #include "DQMServices/Core/interface/DQMStore.h"
22 
23 #include "DataFormats/EcalRawData/interface/EcalRawDataCollections.h"
24 #include "DataFormats/EcalDetId/interface/EBDetId.h"
25 #include "DataFormats/EcalDigi/interface/EBDataFrame.h"
26 #include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
27 #include "DataFormats/EcalRecHit/interface/EcalUncalibratedRecHit.h"
28 #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
29 
30 #include "DQM/EcalCommon/interface/Numbers.h"
31 
32 #include "DQM/EcalBarrelMonitorTasks/interface/EBTestPulseTask.h"
33 
34 EBTestPulseTask::EBTestPulseTask(const edm::ParameterSet& ps){
35 
36  init_ = false;
37 
38  dqmStore_ = edm::Service<DQMStore>().operator->();
39 
40  prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
41 
42  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
43 
44  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
45 
46  EcalRawDataCollection_ = ps.getParameter<edm::InputTag>("EcalRawDataCollection");
47  EBDigiCollection_ = ps.getParameter<edm::InputTag>("EBDigiCollection");
48  EcalPnDiodeDigiCollection_ = ps.getParameter<edm::InputTag>("EcalPnDiodeDigiCollection");
49  EcalUncalibratedRecHitCollection_ = ps.getParameter<edm::InputTag>("EcalUncalibratedRecHitCollection");
50 
51  MGPAGains_.reserve(3);
52  for ( unsigned int i = 1; i <= 3; i++ ) MGPAGains_.push_back(i);
53  MGPAGains_ = ps.getUntrackedParameter<std::vector<int> >("MGPAGains", MGPAGains_);
54 
55  MGPAGainsPN_.reserve(2);
56  for ( unsigned int i = 1; i <= 3; i++ ) MGPAGainsPN_.push_back(i);
57  MGPAGainsPN_ = ps.getUntrackedParameter<std::vector<int> >("MGPAGainsPN", MGPAGainsPN_);
58 
59  for (int i = 0; i < 36; i++) {
60  meShapeMapG01_[i] = 0;
61  meAmplMapG01_[i] = 0;
62  meShapeMapG06_[i] = 0;
63  meAmplMapG06_[i] = 0;
64  meShapeMapG12_[i] = 0;
65  meAmplMapG12_[i] = 0;
66  mePnAmplMapG01_[i] = 0;
67  mePnPedMapG01_[i] = 0;
68  mePnAmplMapG16_[i] = 0;
69  mePnPedMapG16_[i] = 0;
70  }
71 
72 
73 }
74 
75 EBTestPulseTask::~EBTestPulseTask(){
76 
77 }
78 
79 void EBTestPulseTask::beginJob(void){
80 
81  ievt_ = 0;
82 
83  if ( dqmStore_ ) {
84  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask");
85  dqmStore_->rmdir(prefixME_ + "/EBTestPulseTask");
86  }
87 
88 }
89 
90 void EBTestPulseTask::beginRun(const edm::Run& r, const edm::EventSetup& c) {
91 
92  Numbers::initGeometry(c, false);
93 
94  if ( ! mergeRuns_ ) this->reset();
95 
96 }
97 
98 void EBTestPulseTask::endRun(const edm::Run& r, const edm::EventSetup& c) {
99 
100 }
101 
102 void EBTestPulseTask::reset(void) {
103 
104  for (int i = 0; i < 36; i++) {
105  if (find(MGPAGains_.begin(), MGPAGains_.end(), 1) != MGPAGains_.end() ) {
106  if ( meShapeMapG01_[i] ) meShapeMapG01_[i]->Reset();
107  if ( meAmplMapG01_[i] ) meAmplMapG01_[i]->Reset();
108  }
109  if (find(MGPAGains_.begin(), MGPAGains_.end(), 6) != MGPAGains_.end() ) {
110  if ( meShapeMapG06_[i] ) meShapeMapG06_[i]->Reset();
111  if ( meAmplMapG06_[i] ) meAmplMapG06_[i]->Reset();
112  }
113  if (find(MGPAGains_.begin(), MGPAGains_.end(), 12) != MGPAGains_.end() ) {
114  if ( meShapeMapG12_[i] ) meShapeMapG12_[i]->Reset();
115  if ( meAmplMapG12_[i] ) meAmplMapG12_[i]->Reset();
116  }
117  if (find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), 1) != MGPAGainsPN_.end() ) {
118  if ( mePnAmplMapG01_[i] ) mePnAmplMapG01_[i]->Reset();
119  if ( mePnPedMapG01_[i] ) mePnPedMapG01_[i]->Reset();
120  }
121  if (find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), 16) != MGPAGainsPN_.end() ) {
122  if ( mePnAmplMapG16_[i] ) mePnAmplMapG16_[i]->Reset();
123  if ( mePnPedMapG16_[i] ) mePnPedMapG16_[i]->Reset();
124  }
125  }
126 
127 }
128 
129 void EBTestPulseTask::setup(void){
130 
131  init_ = true;
132 
133  std::string name;
134  std::stringstream GainN, GN;
135 
136  if ( dqmStore_ ) {
137  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask");
138 
139  if (find(MGPAGains_.begin(), MGPAGains_.end(), 1) != MGPAGains_.end() ) {
140 
141  GainN.str("");
142  GainN << "Gain" << std::setw(2) << std::setfill('0') << 1;
143  GN.str("");
144  GN << "G" << std::setw(2) << std::setfill('0') << 1;
145 
146  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask/" + GainN.str());
147  for (int i = 0; i < 36; i++) {
148  name = "EBTPT shape " + Numbers::sEB(i+1) + " " + GN.str();
149  meShapeMapG01_[i] = dqmStore_->bookProfile2D(name, name, 1700, 0., 1700., 10, 0., 10., 4096, 0., 4096., "s");
150  meShapeMapG01_[i]->setAxisTitle("channel", 1);
151  meShapeMapG01_[i]->setAxisTitle("sample", 2);
152  meShapeMapG01_[i]->setAxisTitle("amplitude", 3);
153  dqmStore_->tag(meShapeMapG01_[i], i+1);
154 
155  name = "EBTPT amplitude " + Numbers::sEB(i+1) + " " + GN.str(),
156  meAmplMapG01_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 4096, 0., 4096.*12., "s");
157  meAmplMapG01_[i]->setAxisTitle("ieta", 1);
158  meAmplMapG01_[i]->setAxisTitle("iphi", 2);
159  dqmStore_->tag(meAmplMapG01_[i], i+1);
160  }
161 
162  }
163 
164  if (find(MGPAGains_.begin(), MGPAGains_.end(), 6) != MGPAGains_.end() ) {
165 
166  GainN.str("");
167  GainN << "Gain" << std::setw(2) << std::setfill('0') << 6;
168  GN.str("");
169  GN << "G" << std::setw(2) << std::setfill('0') << 6;
170 
171  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask/" + GainN.str());
172  for (int i = 0; i < 36; i++) {
173  name = "EBTPT shape " + Numbers::sEB(i+1) + " " + GN.str();
174  meShapeMapG06_[i] = dqmStore_->bookProfile2D(name, name, 1700, 0., 1700., 10, 0., 10., 4096, 0., 4096., "s");
175  meShapeMapG06_[i]->setAxisTitle("channel", 1);
176  meShapeMapG06_[i]->setAxisTitle("sample", 2);
177  meShapeMapG06_[i]->setAxisTitle("amplitude", 3);
178  dqmStore_->tag(meShapeMapG06_[i], i+1);
179 
180  name = "EBTPT amplitude " + Numbers::sEB(i+1) + " " + GN.str(),
181  meAmplMapG06_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 4096, 0., 4096.*12., "s");
182  meAmplMapG06_[i]->setAxisTitle("ieta", 1);
183  meAmplMapG06_[i]->setAxisTitle("iphi", 2);
184  dqmStore_->tag(meAmplMapG06_[i], i+1);
185  }
186 
187  }
188 
189  if (find(MGPAGains_.begin(), MGPAGains_.end(), 12) != MGPAGains_.end() ) {
190 
191  GainN.str("");
192  GainN << "Gain" << std::setw(2) << std::setfill('0') << 12;
193  GN.str("");
194  GN << "G" << std::setw(2) << std::setfill('0') << 12;
195 
196  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask/" + GainN.str());
197  for (int i = 0; i < 36; i++) {
198  name = "EBTPT shape " + Numbers::sEB(i+1) + " " + GN.str();
199  meShapeMapG12_[i] = dqmStore_->bookProfile2D(name, name, 1700, 0., 1700., 10, 0., 10., 4096, 0., 4096., "s");
200  meShapeMapG12_[i]->setAxisTitle("channel", 1);
201  meShapeMapG12_[i]->setAxisTitle("sample", 2);
202  meShapeMapG12_[i]->setAxisTitle("amplitude", 3);
203  dqmStore_->tag(meShapeMapG12_[i], i+1);
204 
205  name = "EBTPT amplitude " + Numbers::sEB(i+1) + " " + GN.str(),
206  meAmplMapG12_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 4096, 0., 4096.*12., "s");
207  meAmplMapG12_[i]->setAxisTitle("ieta", 1);
208  meAmplMapG12_[i]->setAxisTitle("iphi", 2);
209  dqmStore_->tag(meAmplMapG12_[i], i+1);
210  }
211 
212  }
213 
214  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask/PN");
215 
216  if (find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), 1) != MGPAGainsPN_.end() ) {
217 
218  GainN.str("");
219  GainN << "Gain" << std::setw(2) << std::setfill('0') << 1;
220  GN.str("");
221  GN << "G" << std::setw(2) << std::setfill('0') << 1;
222 
223  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask/PN/" + GainN.str());
224  for (int i = 0; i < 36; i++) {
225  name = "EBTPT PNs amplitude " + Numbers::sEB(i+1) + " " + GN.str();
226  mePnAmplMapG01_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
227  mePnAmplMapG01_[i]->setAxisTitle("channel", 1);
228  mePnAmplMapG01_[i]->setAxisTitle("amplitude", 2);
229  dqmStore_->tag(mePnAmplMapG01_[i], i+1);
230 
231  name = "EBTPT PNs pedestal " + Numbers::sEB(i+1) + " " + GN.str();
232  mePnPedMapG01_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
233  mePnPedMapG01_[i]->setAxisTitle("channel", 1);
234  mePnPedMapG01_[i]->setAxisTitle("pedestal", 2);
235  dqmStore_->tag(mePnPedMapG01_[i], i+1);
236  }
237 
238  }
239 
240  if (find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), 16) != MGPAGainsPN_.end() ) {
241 
242  GainN.str("");
243  GainN << "Gain" << std::setw(2) << std::setfill('0') << 16;
244  GN.str("");
245  GN << "G" << std::setw(2) << std::setfill('0') << 16;
246 
247  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask/PN/" + GainN.str());
248  for (int i = 0; i < 36; i++) {
249  name = "EBTPT PNs amplitude " + Numbers::sEB(i+1) + " " + GN.str();
250  mePnAmplMapG16_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
251  mePnAmplMapG16_[i]->setAxisTitle("channel", 1);
252  mePnAmplMapG16_[i]->setAxisTitle("amplitude", 2);
253  dqmStore_->tag(mePnAmplMapG16_[i], i+1);
254 
255  name = "EBTPT PNs pedestal " + Numbers::sEB(i+1) + " " + GN.str();
256  mePnPedMapG16_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
257  mePnPedMapG16_[i]->setAxisTitle("channel", 1);
258  mePnPedMapG16_[i]->setAxisTitle("pedestal", 2);
259  dqmStore_->tag(mePnPedMapG16_[i], i+1);
260  }
261 
262  }
263 
264  }
265 
266 }
267 
268 void EBTestPulseTask::cleanup(void){
269 
270  if ( ! init_ ) return;
271 
272  if ( dqmStore_ ) {
273  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask");
274 
275  if (find(MGPAGains_.begin(), MGPAGains_.end(), 1) != MGPAGains_.end() ) {
276 
277  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask/Gain01");
278  for (int i = 0; i < 36; i++) {
279  if ( meShapeMapG01_[i] ) dqmStore_->removeElement( meShapeMapG01_[i]->getName() );
280  meShapeMapG01_[i] = 0;
281  if ( meAmplMapG01_[i] ) dqmStore_->removeElement( meAmplMapG01_[i]->getName() );
282  meAmplMapG01_[i] = 0;
283  }
284 
285  }
286 
287  if (find(MGPAGains_.begin(), MGPAGains_.end(), 6) != MGPAGains_.end() ) {
288 
289  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask/Gain06");
290  for (int i = 0; i < 36; i++) {
291  if ( meShapeMapG06_[i] ) dqmStore_->removeElement( meShapeMapG06_[i]->getName() );
292  meShapeMapG06_[i] = 0;
293  if ( meAmplMapG06_[i] ) dqmStore_->removeElement( meAmplMapG06_[i]->getName() );
294  meAmplMapG06_[i] = 0;
295  }
296 
297  }
298 
299  if (find(MGPAGains_.begin(), MGPAGains_.end(), 12) != MGPAGains_.end() ) {
300 
301  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask/Gain12");
302  for (int i = 0; i < 36; i++) {
303  if ( meShapeMapG12_[i] ) dqmStore_->removeElement( meShapeMapG12_[i]->getName() );
304  meShapeMapG12_[i] = 0;
305  if ( meAmplMapG12_[i] ) dqmStore_->removeElement( meAmplMapG12_[i]->getName() );
306  meAmplMapG12_[i] = 0;
307  }
308 
309  }
310 
311  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask/PN");
312 
313  if (find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), 1) != MGPAGainsPN_.end() ) {
314 
315  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask/PN/Gain01");
316  for (int i = 0; i < 36; i++) {
317  if ( mePnAmplMapG01_[i] ) dqmStore_->removeElement( mePnAmplMapG01_[i]->getName() );
318  mePnAmplMapG01_[i] = 0;
319  if ( mePnPedMapG01_[i] ) dqmStore_->removeElement( mePnPedMapG01_[i]->getName() );
320  mePnPedMapG01_[i] = 0;
321  }
322 
323  }
324 
325 
326  if (find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), 16) != MGPAGainsPN_.end() ) {
327 
328  dqmStore_->setCurrentFolder(prefixME_ + "/EBTestPulseTask/PN/Gain16");
329  for (int i = 0; i < 36; i++) {
330  if ( mePnAmplMapG16_[i] ) dqmStore_->removeElement( mePnAmplMapG16_[i]->getName() );
331  mePnAmplMapG16_[i] = 0;
332  if ( mePnPedMapG16_[i] ) dqmStore_->removeElement( mePnPedMapG16_[i]->getName() );
333  mePnPedMapG16_[i] = 0;
334  }
335 
336  }
337 
338  }
339 
340  init_ = false;
341 
342 }
343 
344 void EBTestPulseTask::endJob(void){
345 
346  edm::LogInfo("EBTestPulseTask") << "analyzed " << ievt_ << " events";
347 
348  if ( enableCleanup_ ) this->cleanup();
349 
350 }
351 
352 void EBTestPulseTask::analyze(const edm::Event& e, const edm::EventSetup& c){
353 
354  bool enable = false;
355  int runType[36];
356  for (int i=0; i<36; i++) runType[i] = -1;
357  int mgpaGain[36];
358  for (int i=0; i<36; i++) mgpaGain[i] = -1;
359 
360  edm::Handle<EcalRawDataCollection> dcchs;
361 
362  if ( e.getByLabel(EcalRawDataCollection_, dcchs) ) {
363 
364  for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
365 
366  if ( Numbers::subDet( *dcchItr ) != EcalBarrel ) continue;
367 
368  int ism = Numbers::iSM( *dcchItr, EcalBarrel );
369 
370  runType[ism-1] = dcchItr->getRunType();
371  mgpaGain[ism-1] = dcchItr->getMgpaGain();
372 
373  if ( dcchItr->getRunType() == EcalDCCHeaderBlock::TESTPULSE_MGPA ||
374  dcchItr->getRunType() == EcalDCCHeaderBlock::TESTPULSE_GAP ) enable = true;
375 
376  }
377 
378  } else {
379 
380  edm::LogWarning("EBTestPulseTask") << EcalRawDataCollection_ << " not available";
381 
382  }
383 
384  if ( ! enable ) return;
385 
386  if ( ! init_ ) this->setup();
387 
388  ievt_++;
389 
390  edm::Handle<EBDigiCollection> digis;
391 
392  if ( e.getByLabel(EBDigiCollection_, digis) ) {
393 
394  int nebd = digis->size();
395  LogDebug("EBTestPulseTask") << "event " << ievt_ << " digi collection size " << nebd;
396 
397  for ( EBDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr ) {
398 
399  EBDetId id = digiItr->id();
400 
401  int ic = id.ic();
402 
403  int ism = Numbers::iSM( id );
404 
405  if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::TESTPULSE_MGPA ||
406  runType[ism-1] == EcalDCCHeaderBlock::TESTPULSE_GAP ) ) continue;
407 
408  EBDataFrame dataframe = (*digiItr);
409 
410  for (int i = 0; i < 10; i++) {
411 
412  int adc = dataframe.sample(i).adc();
413  float gain = 1.;
414 
415  MonitorElement* meShapeMap = 0;
416 
417  if ( dataframe.sample(i).gainId() == 1 ) gain = 1./12.;
418  if ( dataframe.sample(i).gainId() == 2 ) gain = 1./ 6.;
419  if ( dataframe.sample(i).gainId() == 3 ) gain = 1./ 1.;
420 
421  if ( mgpaGain[ism-1] == 3 ) meShapeMap = meShapeMapG01_[ism-1];
422  if ( mgpaGain[ism-1] == 2 ) meShapeMap = meShapeMapG06_[ism-1];
423  if ( mgpaGain[ism-1] == 1 ) meShapeMap = meShapeMapG12_[ism-1];
424 
425 // float xval = float(adc) * gain;
426  float xval = float(adc);
427 
428  if ( meShapeMap ) meShapeMap->Fill(ic - 0.5, i + 0.5, xval);
429 
430  }
431 
432  }
433 
434  } else {
435 
436  edm::LogWarning("EBTestPulseTask") << EBDigiCollection_ << " not available";
437 
438  }
439 
440  edm::Handle<EcalUncalibratedRecHitCollection> hits;
441 
442  if ( e.getByLabel(EcalUncalibratedRecHitCollection_, hits) ) {
443 
444  int neh = hits->size();
445  LogDebug("EBTestPulseTask") << "event " << ievt_ << " hits collection size " << neh;
446 
447  for ( EcalUncalibratedRecHitCollection::const_iterator hitItr = hits->begin(); hitItr != hits->end(); ++hitItr ) {
448 
449  EBDetId id = hitItr->id();
450 
451  int ic = id.ic();
452  int ie = (ic-1)/20 + 1;
453  int ip = (ic-1)%20 + 1;
454 
455  int ism = Numbers::iSM( id );
456 
457  float xie = ie - 0.5;
458  float xip = ip - 0.5;
459 
460  if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::TESTPULSE_MGPA ||
461  runType[ism-1] == EcalDCCHeaderBlock::TESTPULSE_GAP ) ) continue;
462 
463  MonitorElement* meAmplMap = 0;
464 
465  if ( mgpaGain[ism-1] == 3 ) meAmplMap = meAmplMapG01_[ism-1];
466  if ( mgpaGain[ism-1] == 2 ) meAmplMap = meAmplMapG06_[ism-1];
467  if ( mgpaGain[ism-1] == 1 ) meAmplMap = meAmplMapG12_[ism-1];
468 
469  float xval = hitItr->amplitude();
470  if ( xval <= 0. ) xval = 0.0;
471 
472 // if ( mgpaGain[ism-1] == 3 ) xval = xval * 1./12.;
473 // if ( mgpaGain[ism-1] == 2 ) xval = xval * 1./ 2.;
474 // if ( mgpaGain[ism-1] == 1 ) xval = xval * 1./ 1.;
475 
476  if ( meAmplMap ) meAmplMap->Fill(xie, xip, xval);
477 
478  }
479 
480  } else {
481 
482  edm::LogWarning("EBTestPulseTask") << EcalUncalibratedRecHitCollection_ << " not available";
483 
484  }
485 
486  edm::Handle<EcalPnDiodeDigiCollection> pns;
487 
488  if ( e.getByLabel(EcalPnDiodeDigiCollection_, pns) ) {
489 
490  int nep = pns->size();
491  LogDebug("EBTestPulseTask") << "event " << ievt_ << " pns collection size " << nep;
492 
493  for ( EcalPnDiodeDigiCollection::const_iterator pnItr = pns->begin(); pnItr != pns->end(); ++pnItr ) {
494 
495  if ( Numbers::subDet( pnItr->id() ) != EcalBarrel ) continue;
496 
497  int ism = Numbers::iSM( pnItr->id() );
498 
499  int num = pnItr->id().iPnId();
500 
501  if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::TESTPULSE_MGPA ||
502  runType[ism-1] == EcalDCCHeaderBlock::TESTPULSE_GAP ) ) continue;
503 
504  float xvalped = 0.;
505 
506  for (int i = 0; i < 4; i++) {
507 
508  int adc = pnItr->sample(i).adc();
509 
510  MonitorElement* mePNPed = 0;
511 
512  if ( pnItr->sample(i).gainId() == 0 ) mePNPed = mePnPedMapG01_[ism-1];
513  if ( pnItr->sample(i).gainId() == 1 ) mePNPed = mePnPedMapG16_[ism-1];
514 
515  float xval = float(adc);
516 
517  if ( mePNPed ) mePNPed->Fill(num - 0.5, xval);
518 
519  xvalped = xvalped + xval;
520 
521  }
522 
523  xvalped = xvalped / 4;
524 
525  float xvalmax = 0.;
526 
527  MonitorElement* mePN = 0;
528 
529  for (int i = 0; i < 50; i++) {
530 
531  int adc = pnItr->sample(i).adc();
532 
533  float xval = float(adc);
534 
535  if ( xval >= xvalmax ) xvalmax = xval;
536 
537  }
538 
539  xvalmax = xvalmax - xvalped;
540 
541  if ( pnItr->sample(0).gainId() == 0 ) mePN = mePnAmplMapG01_[ism-1];
542  if ( pnItr->sample(0).gainId() == 1 ) mePN = mePnAmplMapG16_[ism-1];
543 
544  if ( mePN ) mePN->Fill(num - 0.5, xvalmax);
545 
546  }
547 
548  } else {
549 
550  edm::LogWarning("EBTestPulseTask") << EcalPnDiodeDigiCollection_ << " not available";
551 
552  }
553 
554 }
555 
ecalMGPA::adc
int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:11
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:516
EBTestPulseTask::endJob
void endJob(void)
EndJob.
Definition: EBTestPulseTask.cc:344
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
EBTestPulseTask::meAmplMapG06_
MonitorElement * meAmplMapG06_[36]
Definition: EBTestPulseTask.h:81
i
int i
Definition: DBlmapReader.cc:9
Numbers::sEB
static std::string sEB(const int ism)
Definition: Numbers.cc:87
EBTestPulseTask::EBTestPulseTask
EBTestPulseTask(const edm::ParameterSet &ps)
Constructor.
Definition: EBTestPulseTask.cc:34
EBTestPulseTask::dqmStore_
DQMStore * dqmStore_
Definition: EBTestPulseTask.h:61
EBTestPulseTask::meShapeMapG06_
MonitorElement * meShapeMapG06_[36]
Definition: EBTestPulseTask.h:77
edm::DataFrameContainer::const_iterator
boost::transform_iterator< IterHelp, boost::counting_iterator< int > > const_iterator
Definition: DataFrameContainer.h:51
EBTestPulseTask::mePnAmplMapG01_
MonitorElement * mePnAmplMapG01_[36]
Definition: EBTestPulseTask.h:84
DQMStore::rmdir
void rmdir(const std::string &fullpath)
Definition: DQMStore.cc:2311
EBTestPulseTask::meAmplMapG01_
MonitorElement * meAmplMapG01_[36]
Definition: EBTestPulseTask.h:80
EBTestPulseTask::reset
void reset(void)
Reset.
Definition: EBTestPulseTask.cc:102
Numbers.h
Some &quot;id&quot; conversions.
EBTestPulseTask::mePnAmplMapG16_
MonitorElement * mePnAmplMapG16_[36]
Definition: EBTestPulseTask.h:85
EBTestPulseTask::cleanup
void cleanup(void)
Cleanup.
Definition: EBTestPulseTask.cc:268
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:81
EcalDataFrame::sample
EcalMGPASample sample(int i) const
Definition: EcalDataFrame.h:28
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
EcalMGPASample::gainId
int gainId() const
get the gainId (2 bits)
Definition: EcalMGPASample.h:37
EBTestPulseTask::beginJob
void beginJob(void)
BeginJob.
Definition: EBTestPulseTask.cc:79
EBTestPulseTask::setup
void setup(void)
Setup.
Definition: EBTestPulseTask.cc:129
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:147
EBTestPulseTask::meShapeMapG12_
MonitorElement * meShapeMapG12_[36]
Definition: EBTestPulseTask.h:78
DQMStore::tag
void tag(MonitorElement *me, unsigned int myTag)
Definition: DQMStore.cc:1151
Numbers::iSM
static int iSM(const int ism, const EcalSubdetector subdet)
Definition: Numbers.cc:222
DQMStore::removeElement
void removeElement(const std::string &name)
Definition: DQMStore.cc:2353
EBTestPulseTask::meShapeMapG01_
MonitorElement * meShapeMapG01_[36]
Definition: EBTestPulseTask.h:76
EBTestPulseTask::EcalRawDataCollection_
edm::InputTag EcalRawDataCollection_
Definition: EBTestPulseTask.h:69
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::getName
std::string getName(Reflex::Type &cc)
Definition: ClassFiller.cc:18
EBTestPulseTask::meAmplMapG12_
MonitorElement * meAmplMapG12_[36]
Definition: EBTestPulseTask.h:82
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
EBTestPulseTask::mePnPedMapG01_
MonitorElement * mePnPedMapG01_[36]
Definition: EBTestPulseTask.h:87
EBTestPulseTask::~EBTestPulseTask
virtual ~EBTestPulseTask()
Destructor.
Definition: EBTestPulseTask.cc:75
Numbers::initGeometry
static void initGeometry(const edm::EventSetup &setup, bool verbose=false)
Definition: Numbers.cc:45
evf::utils::num
long long int num
Definition: procUtils.cc:71
EBTestPulseTask::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
BeginRun.
Definition: EBTestPulseTask.cc:90
EBTestPulseTask::mePnPedMapG16_
MonitorElement * mePnPedMapG16_[36]
Definition: EBTestPulseTask.h:88
EBTestPulseTask::MGPAGainsPN_
std::vector< int > MGPAGainsPN_
Definition: EBTestPulseTask.h:74
EBTestPulseTask::EBDigiCollection_
edm::InputTag EBDigiCollection_
Definition: EBTestPulseTask.h:70
EBTestPulseTask::EcalUncalibratedRecHitCollection_
edm::InputTag EcalUncalibratedRecHitCollection_
Definition: EBTestPulseTask.h:72
Numbers::subDet
static EcalSubdetector subDet(const EBDetId &id)
Definition: Numbers.cc:136
EBTestPulseTask::EcalPnDiodeDigiCollection_
edm::InputTag EcalPnDiodeDigiCollection_
Definition: EBTestPulseTask.h:71
ecalpyutils::ism
int ism(int ieta, int iphi)
Definition: EcalPyUtils.cc:56
EBTestPulseTask::MGPAGains_
std::vector< int > MGPAGains_
Definition: EBTestPulseTask.h:73
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
EBTestPulseTask::endRun
void endRun(const edm::Run &r, const edm::EventSetup &c)
EndRun.
Definition: EBTestPulseTask.cc:98
edm::Run
Definition: Run.h:32
EcalMGPASample::adc
int adc() const
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:35
EBTestPulseTask::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Analyze.
Definition: EBTestPulseTask.cc:352
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
EBTestPulseTask::prefixME_
std::string prefixME_
Definition: EBTestPulseTask.h:63
Generated for CMSSW Reference Manual by doxygen 1.8.5