Reference Manual

Go to the documentation of this file.
00001 /*
00002  * \file EBLaserTask.cc
00003  *
00004  * $Date: 2011/08/30 09:30:32 $
00005  * $Revision: 1.134 $
00006  * \author G. Della Ricca
00007  *
00008 */
00009 
00010 #include <iostream>
00011 #include <fstream>
00012 #include <vector>
00013 #include <algorithm>
00014 
00015 #include "FWCore/ServiceRegistry/interface/Service.h"
00016 #include "FWCore/MessageLogger/interface/MessageLogger.h"
00017 
00018 #include "DQMServices/Core/interface/MonitorElement.h"
00019 
00020 #include "DQMServices/Core/interface/DQMStore.h"
00021 
00022 #include "DataFormats/EcalRawData/interface/EcalRawDataCollections.h"
00023 #include "DataFormats/EcalDetId/interface/EBDetId.h"
00024 #include "DataFormats/EcalDigi/interface/EBDataFrame.h"
00025 #include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
00026 #include "DataFormats/EcalRecHit/interface/EcalUncalibratedRecHit.h"
00027 #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
00028 
00029 #include "DQM/EcalCommon/interface/Numbers.h"
00030 
00031 #include "DQM/EcalBarrelMonitorTasks/interface/EBLaserTask.h"
00032 
00033 EBLaserTask::EBLaserTask(const edm::ParameterSet& ps){
00034 
00035   init_ = false;
00036 
00037   dqmStore_ = edm::Service<DQMStore>().operator->();
00038 
00039   prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
00040 
00041   enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
00042 
00043   mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
00044 
00045   EcalRawDataCollection_ = ps.getParameter<edm::InputTag>("EcalRawDataCollection");
00046   EBDigiCollection_ = ps.getParameter<edm::InputTag>("EBDigiCollection");
00047   EcalPnDiodeDigiCollection_ = ps.getParameter<edm::InputTag>("EcalPnDiodeDigiCollection");
00048   EcalUncalibratedRecHitCollection_ = ps.getParameter<edm::InputTag>("EcalUncalibratedRecHitCollection");
00049 
00050   // vector of enabled wavelengths (Default to all 4)
00051   laserWavelengths_.reserve(4);
00052   for ( unsigned int i = 1; i <= 4; i++ ) laserWavelengths_.push_back(i);
00053   laserWavelengths_ = ps.getUntrackedParameter<std::vector<int> >("laserWavelengths", laserWavelengths_);
00054 
00055   for (int i = 0; i < 36; i++) {
00056     meShapeMapL1_[i] = 0;
00057     meAmplMapL1_[i] = 0;
00058     meTimeMapL1_[i] = 0;
00059     meAmplPNMapL1_[i] = 0;
00060     mePnAmplMapG01L1_[i] = 0;
00061     mePnPedMapG01L1_[i] = 0;
00062     mePnAmplMapG16L1_[i] = 0;
00063     mePnPedMapG16L1_[i] = 0;
00064 
00065     meShapeMapL2_[i] = 0;
00066     meAmplMapL2_[i] = 0;
00067     meTimeMapL2_[i] = 0;
00068     meAmplPNMapL2_[i] = 0;
00069     mePnAmplMapG01L2_[i] = 0;
00070     mePnPedMapG01L2_[i] = 0;
00071     mePnAmplMapG16L2_[i] = 0;
00072     mePnPedMapG16L2_[i] = 0;
00073 
00074     meShapeMapL3_[i] = 0;
00075     meAmplMapL3_[i] = 0;
00076     meTimeMapL3_[i] = 0;
00077     meAmplPNMapL3_[i] = 0;
00078     mePnAmplMapG01L3_[i] = 0;
00079     mePnPedMapG01L3_[i] = 0;
00080     mePnAmplMapG16L3_[i] = 0;
00081     mePnPedMapG16L3_[i] = 0;
00082 
00083     meShapeMapL4_[i] = 0;
00084     meAmplMapL4_[i] = 0;
00085     meTimeMapL4_[i] = 0;
00086     meAmplPNMapL4_[i] = 0;
00087     mePnAmplMapG01L4_[i] = 0;
00088     mePnPedMapG01L4_[i] = 0;
00089     mePnAmplMapG16L4_[i] = 0;
00090     mePnPedMapG16L4_[i] = 0;
00091   }
00092 
00093   meAmplSummaryMapL1_ = 0;
00094   meAmplSummaryMapL2_ = 0;
00095   meAmplSummaryMapL3_ = 0;
00096   meAmplSummaryMapL4_ = 0;
00097 
00098 }
00099 
00100 EBLaserTask::~EBLaserTask(){
00101 
00102 }
00103 
00104 void EBLaserTask::beginJob(void){
00105 
00106   ievt_ = 0;
00107 
00108   if ( dqmStore_ ) {
00109     dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask");
00110     dqmStore_->rmdir(prefixME_ + "/EBLaserTask");
00111   }
00112 
00113 }
00114 
00115 void EBLaserTask::beginRun(const edm::Run& r, const edm::EventSetup& c) {
00116 
00117   Numbers::initGeometry(c, false);
00118 
00119   if ( ! mergeRuns_ ) this->reset();
00120 
00121 }
00122 
00123 void EBLaserTask::endRun(const edm::Run& r, const edm::EventSetup& c) {
00124 
00125 }
00126 
00127 void EBLaserTask::reset(void) {
00128 
00129   for (int i = 0; i < 36; i++) {
00130     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
00131       if ( meShapeMapL1_[i] )  meShapeMapL1_[i]->Reset();
00132       if ( meAmplMapL1_[i] ) meAmplMapL1_[i]->Reset();
00133       if ( meTimeMapL1_[i] ) meTimeMapL1_[i]->Reset();
00134       if ( meAmplPNMapL1_[i] ) meAmplPNMapL1_[i]->Reset();
00135     }
00136 
00137     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
00138       if ( meShapeMapL2_[i] )  meShapeMapL2_[i]->Reset();
00139       if ( meAmplMapL2_[i] ) meAmplMapL2_[i]->Reset();
00140       if ( meTimeMapL2_[i] ) meTimeMapL2_[i]->Reset();
00141       if ( meAmplPNMapL2_[i] ) meAmplPNMapL2_[i]->Reset();
00142     }
00143 
00144     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
00145       if ( meShapeMapL3_[i] )  meShapeMapL3_[i]->Reset();
00146       if ( meAmplMapL3_[i] ) meAmplMapL3_[i]->Reset();
00147       if ( meTimeMapL3_[i] ) meTimeMapL3_[i]->Reset();
00148       if ( meAmplPNMapL3_[i] ) meAmplPNMapL3_[i]->Reset();
00149     }
00150 
00151     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
00152       if ( meShapeMapL4_[i] )  meShapeMapL4_[i]->Reset();
00153       if ( meAmplMapL4_[i] ) meAmplMapL4_[i]->Reset();
00154       if ( meTimeMapL4_[i] ) meTimeMapL4_[i]->Reset();
00155       if ( meAmplPNMapL4_[i] ) meAmplPNMapL4_[i]->Reset();
00156     }
00157 
00158     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
00159       if ( mePnAmplMapG01L1_[i] ) mePnAmplMapG01L1_[i]->Reset();
00160       if ( mePnPedMapG01L1_[i] ) mePnPedMapG01L1_[i]->Reset();
00161 
00162       if ( mePnAmplMapG16L1_[i] ) mePnAmplMapG16L1_[i]->Reset();
00163       if ( mePnPedMapG16L1_[i] ) mePnPedMapG16L1_[i]->Reset();
00164     }
00165 
00166     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
00167       if ( mePnAmplMapG01L2_[i] ) mePnAmplMapG01L2_[i]->Reset();
00168       if ( mePnPedMapG01L2_[i] ) mePnPedMapG01L2_[i]->Reset();
00169 
00170       if ( mePnAmplMapG16L2_[i] ) mePnAmplMapG16L2_[i]->Reset();
00171       if ( mePnPedMapG16L2_[i] ) mePnPedMapG16L2_[i]->Reset();
00172     }
00173 
00174     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
00175       if ( mePnAmplMapG01L3_[i] ) mePnAmplMapG01L3_[i]->Reset();
00176       if ( mePnPedMapG01L3_[i] ) mePnPedMapG01L3_[i]->Reset();
00177 
00178       if ( mePnAmplMapG16L3_[i] ) mePnAmplMapG16L3_[i]->Reset();
00179       if ( mePnPedMapG16L3_[i] ) mePnPedMapG16L3_[i]->Reset();
00180     }
00181 
00182     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
00183       if ( mePnAmplMapG01L4_[i] ) mePnAmplMapG01L4_[i]->Reset();
00184       if ( mePnPedMapG01L4_[i] ) mePnPedMapG01L4_[i]->Reset();
00185 
00186       if ( mePnAmplMapG16L4_[i] ) mePnAmplMapG16L4_[i]->Reset();
00187       if ( mePnPedMapG16L4_[i] ) mePnPedMapG16L4_[i]->Reset();
00188     }
00189   }
00190 
00191   if( meAmplSummaryMapL1_ ) meAmplSummaryMapL1_->Reset();
00192   if( meAmplSummaryMapL2_ ) meAmplSummaryMapL2_->Reset();
00193   if( meAmplSummaryMapL3_ ) meAmplSummaryMapL3_->Reset();
00194   if( meAmplSummaryMapL4_ ) meAmplSummaryMapL4_->Reset();
00195 
00196 }
00197 
00198 void EBLaserTask::setup(void){
00199 
00200   init_ = true;
00201 
00202   std::string name;
00203   std::stringstream LaserN, LN;
00204 
00205   if ( dqmStore_ ) {
00206     dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask");
00207 
00208     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
00209 
00210       LaserN.str("");
00211       LaserN << "Laser" << 1;
00212       LN.str("");
00213       LN << "L" << 1;
00214 
00215       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str());
00216       for (int i = 0; i < 36; i++) {
00217         name = "EBLT shape " + Numbers::sEB(i+1) + " " + LN.str();
00218         meShapeMapL1_[i] = dqmStore_->bookProfile2D(name, name, 1700, 0., 1700., 10, 0., 10., 4096, 0., 4096., "s");
00219         meShapeMapL1_[i]->setAxisTitle("channel", 1);
00220         meShapeMapL1_[i]->setAxisTitle("sample", 2);
00221         meShapeMapL1_[i]->setAxisTitle("amplitude", 3);
00222         dqmStore_->tag(meShapeMapL1_[i], i+1);
00223         name = "EBLT amplitude " + Numbers::sEB(i+1) + " " + LN.str();
00224         meAmplMapL1_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 4096, 0., 4096.*12., "s");
00225         meAmplMapL1_[i]->setAxisTitle("ieta", 1);
00226         meAmplMapL1_[i]->setAxisTitle("iphi", 2);
00227         dqmStore_->tag(meAmplMapL1_[i], i+1);
00228         name = "EBLT timing " + Numbers::sEB(i+1) + " " + LN.str();
00229         meTimeMapL1_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 250, 0., 10., "s");
00230         meTimeMapL1_[i]->setAxisTitle("ieta", 1);
00231         meTimeMapL1_[i]->setAxisTitle("iphi", 2);
00232         dqmStore_->tag(meTimeMapL1_[i], i+1);
00233         name = "EBLT amplitude over PN " + Numbers::sEB(i+1) + " " + LN.str();
00234         meAmplPNMapL1_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 4096, 0., 4096.*12., "s");
00235         meAmplPNMapL1_[i]->setAxisTitle("ieta", 1);
00236         meAmplPNMapL1_[i]->setAxisTitle("iphi", 2);
00237         dqmStore_->tag(meAmplPNMapL1_[i], i+1);
00238       }
00239 
00240       name = "EBLT amplitude map " + LN.str();
00241       meAmplSummaryMapL1_ = dqmStore_->bookProfile2D(name, name, 72, 0., 360., 34, -85., 85., 0., 4096.);
00242       meAmplSummaryMapL1_->setAxisTitle("jphi", 1);
00243       meAmplSummaryMapL1_->setAxisTitle("jeta", 2);
00244 
00245       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str() + "/PN");
00246 
00247       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str() + "/PN/Gain01");
00248       for (int i = 0; i < 36; i++) {
00249         name = "EBLT PNs amplitude " + Numbers::sEB(i+1) + " G01 " + LN.str(); 
00250         mePnAmplMapG01L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00251         mePnAmplMapG01L1_[i]->setAxisTitle("channel", 1);
00252         mePnAmplMapG01L1_[i]->setAxisTitle("amplitude", 2);
00253         dqmStore_->tag(mePnAmplMapG01L1_[i], i+1);
00254 
00255         name = "EBLT PNs pedestal " + Numbers::sEB(i+1) + " G01 " + LN.str(); 
00256         mePnPedMapG01L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00257         mePnPedMapG01L1_[i]->setAxisTitle("channel", 1);
00258         mePnPedMapG01L1_[i]->setAxisTitle("pedestal", 2);
00259         dqmStore_->tag(mePnPedMapG01L1_[i], i+1);
00260       }
00261 
00262       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str() + "/PN/Gain16");
00263       for (int i = 0; i < 36; i++) {
00264         name = "EBLT PNs amplitude " + Numbers::sEB(i+1) + " G16 " + LN.str(), 
00265         mePnAmplMapG16L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00266         mePnAmplMapG16L1_[i]->setAxisTitle("channel", 1);
00267         mePnAmplMapG16L1_[i]->setAxisTitle("amplitude", 2);
00268         dqmStore_->tag(mePnAmplMapG16L1_[i], i+1);
00269 
00270         name = "EBLT PNs pedestal " + Numbers::sEB(i+1) + " G16 " + LN.str(); 
00271         mePnPedMapG16L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00272         mePnPedMapG16L1_[i]->setAxisTitle("channel", 1);
00273         mePnPedMapG16L1_[i]->setAxisTitle("pedestal", 2);
00274         dqmStore_->tag(mePnPedMapG16L1_[i], i+1);
00275       }
00276 
00277 
00278     }
00279 
00280     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
00281 
00282       LaserN.str("");
00283       LaserN << "Laser" << 2;
00284       LN.str("");
00285       LN << "L" << 2;
00286 
00287       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str());
00288       for (int i = 0; i < 36; i++) {
00289         name = "EBLT shape " + Numbers::sEB(i+1) + " " + LN.str();
00290         meShapeMapL2_[i] = dqmStore_->bookProfile2D(name, name, 1700, 0., 1700., 10, 0., 10., 4096, 0., 4096., "s");
00291         meShapeMapL2_[i]->setAxisTitle("channel", 1);
00292         meShapeMapL2_[i]->setAxisTitle("sample", 2);
00293         meShapeMapL2_[i]->setAxisTitle("amplitude", 3);
00294         dqmStore_->tag(meShapeMapL2_[i], i+1);
00295         name = "EBLT amplitude " + Numbers::sEB(i+1) + " " + LN.str();
00296         meAmplMapL2_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 4096, 0., 4096.*12., "s");
00297         meAmplMapL2_[i]->setAxisTitle("ieta", 1);
00298         meAmplMapL2_[i]->setAxisTitle("iphi", 2);
00299         dqmStore_->tag(meAmplMapL2_[i], i+1);
00300         name = "EBLT timing " + Numbers::sEB(i+1) + " " + LN.str();
00301         meTimeMapL2_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 250, 0., 10., "s");
00302         meTimeMapL2_[i]->setAxisTitle("ieta", 1);
00303         meTimeMapL2_[i]->setAxisTitle("iphi", 2);
00304         dqmStore_->tag(meTimeMapL2_[i], i+1);
00305         name = "EBLT amplitude over PN " + Numbers::sEB(i+1) + " " + LN.str();
00306         meAmplPNMapL2_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 4096, 0., 4096.*12., "s");
00307         meAmplPNMapL2_[i]->setAxisTitle("ieta", 1);
00308         meAmplPNMapL2_[i]->setAxisTitle("iphi", 2);
00309         dqmStore_->tag(meAmplPNMapL2_[i], i+1);
00310       }
00311 
00312       name = "EBLT amplitude map " + LN.str();
00313       meAmplSummaryMapL2_ = dqmStore_->bookProfile2D(name, name, 72, 0., 360., 34, -85., 85., 0., 4096.);
00314       meAmplSummaryMapL2_->setAxisTitle("jphi", 1);
00315       meAmplSummaryMapL2_->setAxisTitle("jeta", 2);
00316 
00317       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str() + "/PN");
00318 
00319       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str() + "/PN/Gain01");
00320       for (int i = 0; i < 36; i++) {
00321         name = "EBLT PNs amplitude " + Numbers::sEB(i+1) + " G01 " + LN.str(); 
00322         mePnAmplMapG01L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00323         mePnAmplMapG01L2_[i]->setAxisTitle("channel", 1);
00324         mePnAmplMapG01L2_[i]->setAxisTitle("amplitude", 2);
00325         dqmStore_->tag(mePnAmplMapG01L2_[i], i+1);
00326 
00327         name = "EBLT PNs pedestal " + Numbers::sEB(i+1) + " G01 " + LN.str(); 
00328         mePnPedMapG01L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00329         mePnPedMapG01L2_[i]->setAxisTitle("channel", 1);
00330         mePnPedMapG01L2_[i]->setAxisTitle("pedestal", 2);
00331         dqmStore_->tag(mePnPedMapG01L2_[i], i+1);
00332       }
00333 
00334       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str() + "/PN/Gain16");
00335       for (int i = 0; i < 36; i++) {
00336         name = "EBLT PNs amplitude " + Numbers::sEB(i+1) + " G16 " + LN.str(), 
00337         mePnAmplMapG16L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00338         mePnAmplMapG16L2_[i]->setAxisTitle("channel", 1);
00339         mePnAmplMapG16L2_[i]->setAxisTitle("amplitude", 2);
00340         dqmStore_->tag(mePnAmplMapG16L2_[i], i+1);
00341 
00342         name = "EBLT PNs pedestal " + Numbers::sEB(i+1) + " G16 " + LN.str(); 
00343         mePnPedMapG16L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00344         mePnPedMapG16L2_[i]->setAxisTitle("channel", 1);
00345         mePnPedMapG16L2_[i]->setAxisTitle("pedestal", 2);
00346         dqmStore_->tag(mePnPedMapG16L2_[i], i+1);
00347       }
00348 
00349 
00350     }
00351 
00352     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
00353 
00354       LaserN.str("");
00355       LaserN << "Laser" << 3;
00356       LN.str("");
00357       LN << "L" << 3;
00358 
00359       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str());
00360       for (int i = 0; i < 36; i++) {
00361         name = "EBLT shape " + Numbers::sEB(i+1) + " " + LN.str();
00362         meShapeMapL3_[i] = dqmStore_->bookProfile2D(name, name, 1700, 0., 1700., 10, 0., 10., 4096, 0., 4096., "s");
00363         meShapeMapL3_[i]->setAxisTitle("channel", 1);
00364         meShapeMapL3_[i]->setAxisTitle("sample", 2);
00365         meShapeMapL3_[i]->setAxisTitle("amplitude", 3);
00366         dqmStore_->tag(meShapeMapL3_[i], i+1);
00367         name = "EBLT amplitude " + Numbers::sEB(i+1) + " " + LN.str();
00368         meAmplMapL3_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 4096, 0., 4096.*12., "s");
00369         meAmplMapL3_[i]->setAxisTitle("ieta", 1);
00370         meAmplMapL3_[i]->setAxisTitle("iphi", 2);
00371         dqmStore_->tag(meAmplMapL3_[i], i+1);
00372         name = "EBLT timing " + Numbers::sEB(i+1) + " " + LN.str();
00373         meTimeMapL3_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 250, 0., 10., "s");
00374         meTimeMapL3_[i]->setAxisTitle("ieta", 1);
00375         meTimeMapL3_[i]->setAxisTitle("iphi", 2);
00376         dqmStore_->tag(meTimeMapL3_[i], i+1);
00377         name = "EBLT amplitude over PN " + Numbers::sEB(i+1) + " " + LN.str();
00378         meAmplPNMapL3_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 4096, 0., 4096.*12., "s");
00379         meAmplPNMapL3_[i]->setAxisTitle("ieta", 1);
00380         meAmplPNMapL3_[i]->setAxisTitle("iphi", 2);
00381         dqmStore_->tag(meAmplPNMapL3_[i], i+1);
00382       }
00383 
00384       name = "EBLT amplitude map " + LN.str();
00385       meAmplSummaryMapL3_ = dqmStore_->bookProfile2D(name, name, 72, 0., 360., 34, -85., 85., 0., 4096.);
00386       meAmplSummaryMapL3_->setAxisTitle("jphi", 1);
00387       meAmplSummaryMapL3_->setAxisTitle("jeta", 2);
00388 
00389       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str() + "/PN");
00390 
00391       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str() + "/PN/Gain01");
00392       for (int i = 0; i < 36; i++) {
00393         name = "EBLT PNs amplitude " + Numbers::sEB(i+1) + " G01 " + LN.str(); 
00394         mePnAmplMapG01L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00395         mePnAmplMapG01L3_[i]->setAxisTitle("channel", 1);
00396         mePnAmplMapG01L3_[i]->setAxisTitle("amplitude", 2);
00397         dqmStore_->tag(mePnAmplMapG01L3_[i], i+1);
00398 
00399         name = "EBLT PNs pedestal " + Numbers::sEB(i+1) + " G01 " + LN.str(); 
00400         mePnPedMapG01L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00401         mePnPedMapG01L3_[i]->setAxisTitle("channel", 1);
00402         mePnPedMapG01L3_[i]->setAxisTitle("pedestal", 2);
00403         dqmStore_->tag(mePnPedMapG01L3_[i], i+1);
00404       }
00405 
00406       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str() + "/PN/Gain16");
00407       for (int i = 0; i < 36; i++) {
00408         name = "EBLT PNs amplitude " + Numbers::sEB(i+1) + " G16 " + LN.str(), 
00409         mePnAmplMapG16L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00410         mePnAmplMapG16L3_[i]->setAxisTitle("channel", 1);
00411         mePnAmplMapG16L3_[i]->setAxisTitle("amplitude", 2);
00412         dqmStore_->tag(mePnAmplMapG16L3_[i], i+1);
00413 
00414         name = "EBLT PNs pedestal " + Numbers::sEB(i+1) + " G16 " + LN.str(); 
00415         mePnPedMapG16L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00416         mePnPedMapG16L3_[i]->setAxisTitle("channel", 1);
00417         mePnPedMapG16L3_[i]->setAxisTitle("pedestal", 2);
00418         dqmStore_->tag(mePnPedMapG16L3_[i], i+1);
00419       }
00420 
00421     }
00422 
00423     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
00424 
00425       LaserN.str("");
00426       LaserN << "Laser" << 4;
00427       LN.str("");
00428       LN << "L" << 4;
00429 
00430       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str());
00431       for (int i = 0; i < 36; i++) {
00432         name = "EBLT shape " + Numbers::sEB(i+1) + " " + LN.str();
00433         meShapeMapL4_[i] = dqmStore_->bookProfile2D(name, name, 1700, 0., 1700., 10, 0., 10., 4096, 0., 4096., "s");
00434         meShapeMapL4_[i]->setAxisTitle("channel", 1);
00435         meShapeMapL4_[i]->setAxisTitle("sample", 2);
00436         meShapeMapL4_[i]->setAxisTitle("amplitude", 3);
00437         dqmStore_->tag(meShapeMapL4_[i], i+1);
00438         name = "EBLT amplitude " + Numbers::sEB(i+1) + " " + LN.str();
00439         meAmplMapL4_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 4096, 0., 4096.*12., "s");
00440         meAmplMapL4_[i]->setAxisTitle("ieta", 1);
00441         meAmplMapL4_[i]->setAxisTitle("iphi", 2);
00442         dqmStore_->tag(meAmplMapL4_[i], i+1);
00443         name = "EBLT timing " + Numbers::sEB(i+1) + " " + LN.str();
00444         meTimeMapL4_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 250, 0., 10., "s");
00445         meTimeMapL4_[i]->setAxisTitle("ieta", 1);
00446         meTimeMapL4_[i]->setAxisTitle("iphi", 2);
00447         dqmStore_->tag(meTimeMapL4_[i], i+1);
00448         name = "EBLT amplitude over PN " + Numbers::sEB(i+1) + " " + LN.str();
00449         meAmplPNMapL4_[i] = dqmStore_->bookProfile2D(name, name, 85, 0., 85., 20, 0., 20., 4096, 0., 4096.*12., "s");
00450         meAmplPNMapL4_[i]->setAxisTitle("ieta", 1);
00451         meAmplPNMapL4_[i]->setAxisTitle("iphi", 2);
00452         dqmStore_->tag(meAmplPNMapL4_[i], i+1);
00453       }
00454 
00455       name = "EBLT amplitude map " + LN.str();
00456       meAmplSummaryMapL4_ = dqmStore_->bookProfile2D(name, name, 72, 0., 360., 34, -85., 85., 0., 4096.);
00457       meAmplSummaryMapL4_->setAxisTitle("jphi", 1);
00458       meAmplSummaryMapL4_->setAxisTitle("jeta", 2);
00459 
00460       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str() + "/PN");
00461 
00462       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str() + "/PN/Gain01");
00463       for (int i = 0; i < 36; i++) {
00464         name = "EBLT PNs amplitude " + Numbers::sEB(i+1) + " G01 " + LN.str(); 
00465         mePnAmplMapG01L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00466         mePnAmplMapG01L4_[i]->setAxisTitle("channel", 1);
00467         mePnAmplMapG01L4_[i]->setAxisTitle("amplitude", 2);
00468         dqmStore_->tag(mePnAmplMapG01L4_[i], i+1);
00469 
00470         name = "EBLT PNs pedestal " + Numbers::sEB(i+1) + " G01 " + LN.str(); 
00471         mePnPedMapG01L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00472         mePnPedMapG01L4_[i]->setAxisTitle("channel", 1);
00473         mePnPedMapG01L4_[i]->setAxisTitle("pedestal", 2);
00474         dqmStore_->tag(mePnPedMapG01L4_[i], i+1);
00475       }
00476 
00477       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/" + LaserN.str() + "/PN/Gain16");
00478       for (int i = 0; i < 36; i++) {
00479         name = "EBLT PNs amplitude " + Numbers::sEB(i+1) + " G16 " + LN.str(), 
00480         mePnAmplMapG16L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00481         mePnAmplMapG16L4_[i]->setAxisTitle("channel", 1);
00482         mePnAmplMapG16L4_[i]->setAxisTitle("amplitude", 2);
00483         dqmStore_->tag(mePnAmplMapG16L4_[i], i+1);
00484 
00485         name = "EBLT PNs pedestal " + Numbers::sEB(i+1) + " G16 " + LN.str(); 
00486         mePnPedMapG16L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
00487         mePnPedMapG16L4_[i]->setAxisTitle("channel", 1);
00488         mePnPedMapG16L4_[i]->setAxisTitle("pedestal", 2);
00489         dqmStore_->tag(mePnPedMapG16L4_[i], i+1);
00490       }
00491 
00492     }
00493 
00494   }
00495 
00496 }
00497 
00498 void EBLaserTask::cleanup(void){
00499 
00500   if ( ! init_ ) return;
00501 
00502   if ( dqmStore_ ) {
00503     dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask");
00504 
00505     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
00506       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser1");
00507       for (int i = 0; i < 36; i++) {
00508         if ( meShapeMapL1_[i] )  dqmStore_->removeElement( meShapeMapL1_[i]->getName() );
00509         meShapeMapL1_[i] = 0;
00510         if ( meAmplMapL1_[i] ) dqmStore_->removeElement( meAmplMapL1_[i]->getName() );
00511         meAmplMapL1_[i] = 0;
00512         if ( meTimeMapL1_[i] ) dqmStore_->removeElement( meTimeMapL1_[i]->getName() );
00513         meTimeMapL1_[i] = 0;
00514         if ( meAmplPNMapL1_[i] ) dqmStore_->removeElement( meAmplPNMapL1_[i]->getName() );
00515         meAmplPNMapL1_[i] = 0;
00516       }
00517       if( meAmplSummaryMapL1_ ) dqmStore_->removeElement( meAmplSummaryMapL1_->getName() );
00518       meAmplSummaryMapL1_ = 0;
00519     }
00520 
00521     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
00522       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser2");
00523       for (int i = 0; i < 36; i++) {
00524         if ( meShapeMapL2_[i] )  dqmStore_->removeElement( meShapeMapL2_[i]->getName() );
00525         meShapeMapL2_[i] = 0;
00526         if ( meAmplMapL2_[i] ) dqmStore_->removeElement( meAmplMapL2_[i]->getName() );
00527         meAmplMapL2_[i] = 0;
00528         if ( meTimeMapL2_[i] ) dqmStore_->removeElement( meTimeMapL2_[i]->getName() );
00529         meTimeMapL2_[i] = 0;
00530         if ( meAmplPNMapL2_[i] ) dqmStore_->removeElement( meAmplPNMapL2_[i]->getName() );
00531         meAmplPNMapL2_[i] = 0;
00532       }
00533       if( meAmplSummaryMapL2_ ) dqmStore_->removeElement( meAmplSummaryMapL2_->getName() );
00534       meAmplSummaryMapL2_ = 0;
00535     }
00536 
00537     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
00538       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser3");
00539       for (int i = 0; i < 36; i++) {
00540         if ( meShapeMapL3_[i] )  dqmStore_->removeElement( meShapeMapL3_[i]->getName() );
00541         meShapeMapL3_[i] = 0;
00542         if ( meAmplMapL3_[i] ) dqmStore_->removeElement( meAmplMapL3_[i]->getName() );
00543         meAmplMapL3_[i] = 0;
00544         if ( meTimeMapL3_[i] ) dqmStore_->removeElement( meTimeMapL3_[i]->getName() );
00545         meTimeMapL3_[i] = 0;
00546         if ( meAmplPNMapL3_[i] ) dqmStore_->removeElement( meAmplPNMapL3_[i]->getName() );
00547         meAmplPNMapL3_[i] = 0;
00548       }
00549       if( meAmplSummaryMapL3_ ) dqmStore_->removeElement( meAmplSummaryMapL3_->getName() );
00550       meAmplSummaryMapL3_ = 0;
00551     }
00552 
00553     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
00554       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser4");
00555       for (int i = 0; i < 36; i++) {
00556         if ( meShapeMapL4_[i] )  dqmStore_->removeElement( meShapeMapL4_[i]->getName() );
00557         meShapeMapL4_[i] = 0;
00558         if ( meAmplMapL4_[i] ) dqmStore_->removeElement( meAmplMapL4_[i]->getName() );
00559         meAmplMapL4_[i] = 0;
00560         if ( meTimeMapL4_[i] ) dqmStore_->removeElement( meTimeMapL4_[i]->getName() );
00561         meTimeMapL4_[i] = 0;
00562         if ( meAmplPNMapL4_[i] ) dqmStore_->removeElement( meAmplPNMapL4_[i]->getName() );
00563         meAmplPNMapL4_[i] = 0;
00564       }
00565       if( meAmplSummaryMapL4_ ) dqmStore_->removeElement( meAmplSummaryMapL4_->getName() );
00566       meAmplSummaryMapL4_ = 0;
00567     }
00568 
00569     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
00570       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser1/PN");
00571 
00572       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser1/PN/Gain01");
00573       for (int i = 0; i < 36; i++) {
00574         if ( mePnAmplMapG01L1_[i] ) dqmStore_->removeElement( mePnAmplMapG01L1_[i]->getName() );
00575         mePnAmplMapG01L1_[i] = 0;
00576         if ( mePnPedMapG01L1_[i] ) dqmStore_->removeElement( mePnPedMapG01L1_[i]->getName() );
00577         mePnPedMapG01L1_[i] = 0;
00578       }
00579 
00580       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser1/PN/Gain16");
00581       for (int i = 0; i < 36; i++) {
00582         if ( mePnAmplMapG16L1_[i] ) dqmStore_->removeElement( mePnAmplMapG16L1_[i]->getName() );
00583         mePnAmplMapG16L1_[i] = 0;
00584         if ( mePnPedMapG16L1_[i] ) dqmStore_->removeElement( mePnPedMapG16L1_[i]->getName() );
00585         mePnPedMapG16L1_[i] = 0;
00586       }
00587     }
00588 
00589     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
00590       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser2/PN");
00591 
00592       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser2/PN/Gain01");
00593       for (int i = 0; i < 36; i++) {
00594         if ( mePnAmplMapG01L2_[i] ) dqmStore_->removeElement( mePnAmplMapG01L2_[i]->getName() );
00595         mePnAmplMapG01L2_[i] = 0;
00596         if ( mePnPedMapG01L2_[i] ) dqmStore_->removeElement( mePnPedMapG01L2_[i]->getName() );
00597         mePnPedMapG01L2_[i] = 0;
00598       }
00599 
00600       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser2/PN/Gain16");
00601       for (int i = 0; i < 36; i++) {
00602         if ( mePnAmplMapG16L2_[i] ) dqmStore_->removeElement( mePnAmplMapG16L2_[i]->getName() );
00603         mePnAmplMapG16L2_[i] = 0;
00604         if ( mePnPedMapG16L2_[i] ) dqmStore_->removeElement( mePnPedMapG16L2_[i]->getName() );
00605         mePnPedMapG16L2_[i] = 0;
00606       }
00607     }
00608 
00609     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
00610       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser3/PN");
00611 
00612       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser3/PN/Gain01");
00613       for (int i = 0; i < 36; i++) {
00614         if ( mePnAmplMapG01L3_[i] ) dqmStore_->removeElement( mePnAmplMapG01L3_[i]->getName() );
00615         mePnAmplMapG01L3_[i] = 0;
00616         if ( mePnPedMapG01L3_[i] ) dqmStore_->removeElement( mePnPedMapG01L3_[i]->getName() );
00617         mePnPedMapG01L3_[i] = 0;
00618       }
00619 
00620       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser3/PN/Gain16");
00621       for (int i = 0; i < 36; i++) {
00622         if ( mePnAmplMapG16L3_[i] ) dqmStore_->removeElement( mePnAmplMapG16L3_[i]->getName() );
00623         mePnAmplMapG16L3_[i] = 0;
00624         if ( mePnPedMapG16L3_[i] ) dqmStore_->removeElement( mePnPedMapG16L3_[i]->getName() );
00625         mePnPedMapG16L3_[i] = 0;
00626       }
00627     }
00628 
00629     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
00630       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser4/PN");
00631 
00632       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser4/PN/Gain01");
00633       for (int i = 0; i < 36; i++) {
00634         if ( mePnAmplMapG01L4_[i] ) dqmStore_->removeElement( mePnAmplMapG01L4_[i]->getName() );
00635         mePnAmplMapG01L4_[i] = 0;
00636         if ( mePnPedMapG01L4_[i] ) dqmStore_->removeElement( mePnPedMapG01L4_[i]->getName() );
00637         mePnPedMapG01L4_[i] = 0;
00638       }
00639 
00640       dqmStore_->setCurrentFolder(prefixME_ + "/EBLaserTask/Laser4/PN/Gain16");
00641       for (int i = 0; i < 36; i++) {
00642         if ( mePnAmplMapG16L4_[i] ) dqmStore_->removeElement( mePnAmplMapG16L4_[i]->getName() );
00643         mePnAmplMapG16L4_[i] = 0;
00644         if ( mePnPedMapG16L4_[i] ) dqmStore_->removeElement( mePnPedMapG16L4_[i]->getName() );
00645         mePnPedMapG16L4_[i] = 0;
00646       }
00647 
00648     }
00649 
00650   }
00651 
00652   init_ = false;
00653 
00654 }
00655 
00656 void EBLaserTask::endJob(void){
00657 
00658   edm::LogInfo("EBLaserTask") << "analyzed " << ievt_ << " events";
00659 
00660   if ( enableCleanup_ ) this->cleanup();
00661 
00662 }
00663 
00664 void EBLaserTask::analyze(const edm::Event& e, const edm::EventSetup& c){
00665 
00666   bool enable = false;
00667   int runType[36];
00668   for (int i=0; i<36; i++) runType[i] = -1;
00669   int rtHalf[36];
00670   for (int i=0; i<36; i++) rtHalf[i] = -1;
00671   int waveLength[36];
00672   for (int i=0; i<36; i++) waveLength[i] = -1;
00673 
00674   edm::Handle<EcalRawDataCollection> dcchs;
00675 
00676   if ( e.getByLabel(EcalRawDataCollection_, dcchs) ) {
00677 
00678     for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
00679 
00680       if ( Numbers::subDet( *dcchItr ) != EcalBarrel ) continue;
00681 
00682       int ism = Numbers::iSM( *dcchItr, EcalBarrel );
00683 
00684       runType[ism-1] = dcchItr->getRunType();
00685       rtHalf[ism-1] = dcchItr->getRtHalf();
00686       waveLength[ism-1] = dcchItr->getEventSettings().wavelength;
00687 
00688       if ( dcchItr->getRunType() == EcalDCCHeaderBlock::LASER_STD ||
00689            dcchItr->getRunType() == EcalDCCHeaderBlock::LASER_GAP ) enable = true;
00690 
00691     }
00692 
00693   } else {
00694 
00695     edm::LogWarning("EBLaserTask") << EcalRawDataCollection_ << " not available";
00696 
00697   }
00698 
00699   if ( ! enable ) return;
00700 
00701   if ( ! init_ ) this->setup();
00702 
00703   ievt_++;
00704 
00705   edm::Handle<EBDigiCollection> digis;
00706 
00707   if ( e.getByLabel(EBDigiCollection_, digis) ) {
00708 
00709     int nebd = digis->size();
00710     LogDebug("EBLaserTask") << "event " << ievt_ << " digi collection size " << nebd;
00711 
00712     for ( EBDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr ) {
00713 
00714       EBDetId id = digiItr->id();
00715 
00716       int ic = id.ic();
00717 
00718       int ism = Numbers::iSM( id );
00719 
00720       if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
00721                runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
00722 
00723       if ( rtHalf[ism-1] != Numbers::RtHalf(id) ) continue;
00724 
00725       EBDataFrame dataframe = (*digiItr);
00726 
00727       for (int i = 0; i < 10; i++) {
00728 
00729         int adc = dataframe.sample(i).adc();
00730         float gain = 1.;
00731 
00732         MonitorElement* meShapeMap = 0;
00733 
00734         if ( dataframe.sample(i).gainId() == 1 ) gain = 1./12.;
00735         if ( dataframe.sample(i).gainId() == 2 ) gain = 1./ 6.;
00736         if ( dataframe.sample(i).gainId() == 3 ) gain = 1./ 1.;
00737 
00738         if ( rtHalf[ism-1] == 0 || rtHalf[ism-1] == 1 ) {
00739 
00740           if ( waveLength[ism-1] == 0 ) meShapeMap = meShapeMapL1_[ism-1];
00741           if ( waveLength[ism-1] == 1 ) meShapeMap = meShapeMapL2_[ism-1];
00742           if ( waveLength[ism-1] == 2 ) meShapeMap = meShapeMapL3_[ism-1];
00743           if ( waveLength[ism-1] == 3 ) meShapeMap = meShapeMapL4_[ism-1];
00744 
00745         } else {
00746 
00747           edm::LogWarning("EBLaserTask") << " RtHalf = " << rtHalf[ism-1];
00748 
00749         }
00750 
00751 //        float xval = float(adc) * gain;
00752         float xval = float(adc);
00753 
00754         if ( meShapeMap ) meShapeMap->Fill(ic - 0.5, i + 0.5, xval);
00755 
00756       }
00757 
00758     }
00759 
00760   } else {
00761 
00762     edm::LogWarning("EBLaserTask") << EBDigiCollection_ << " not available";
00763 
00764   }
00765 
00766   float adcA[36];
00767   float adcB[36];
00768 
00769   for ( int i = 0; i < 36; i++ ) {
00770     adcA[i] = 0.;
00771     adcB[i] = 0.;
00772   }
00773 
00774   edm::Handle<EcalPnDiodeDigiCollection> pns;
00775 
00776   if ( e.getByLabel(EcalPnDiodeDigiCollection_, pns) ) {
00777 
00778     int nep = pns->size();
00779     LogDebug("EBLaserTask") << "event " << ievt_ << " pns collection size " << nep;
00780 
00781     for ( EcalPnDiodeDigiCollection::const_iterator pnItr = pns->begin(); pnItr != pns->end(); ++pnItr ) {
00782 
00783       if ( Numbers::subDet( pnItr->id() ) != EcalBarrel ) continue;
00784 
00785       int ism = Numbers::iSM( pnItr->id() );
00786 
00787       int num = pnItr->id().iPnId();
00788 
00789       if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
00790                runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
00791 
00792       float xvalped = 0.;
00793 
00794       for (int i = 0; i < 4; i++) {
00795 
00796         int adc = pnItr->sample(i).adc();
00797 
00798         MonitorElement* mePNPed = 0;
00799 
00800         if ( pnItr->sample(i).gainId() == 0 ) {
00801           if ( waveLength[ism-1] == 0 ) mePNPed = mePnPedMapG01L1_[ism-1];
00802           if ( waveLength[ism-1] == 1 ) mePNPed = mePnPedMapG01L2_[ism-1];
00803           if ( waveLength[ism-1] == 2 ) mePNPed = mePnPedMapG01L3_[ism-1];
00804           if ( waveLength[ism-1] == 3 ) mePNPed = mePnPedMapG01L4_[ism-1];
00805         }
00806         if ( pnItr->sample(i).gainId() == 1 ) {
00807           if ( waveLength[ism-1] == 0 ) mePNPed = mePnPedMapG16L1_[ism-1];
00808           if ( waveLength[ism-1] == 1 ) mePNPed = mePnPedMapG16L2_[ism-1];
00809           if ( waveLength[ism-1] == 2 ) mePNPed = mePnPedMapG16L3_[ism-1];
00810           if ( waveLength[ism-1] == 3 ) mePNPed = mePnPedMapG16L4_[ism-1];
00811         }
00812 
00813         float xval = float(adc);
00814 
00815         if ( mePNPed ) mePNPed->Fill(num - 0.5, xval);
00816 
00817         xvalped = xvalped + xval;
00818 
00819       }
00820 
00821       xvalped = xvalped / 4;
00822 
00823       float xvalmax = 0.;
00824 
00825       MonitorElement* mePN = 0;
00826 
00827       for (int i = 0; i < 50; i++) {
00828 
00829         int adc = pnItr->sample(i).adc();
00830 
00831         float xval = float(adc);
00832 
00833         if ( xval >= xvalmax ) xvalmax = xval;
00834 
00835       }
00836 
00837       xvalmax = xvalmax - xvalped;
00838 
00839       if ( pnItr->sample(0).gainId() == 0 ) {
00840         if ( waveLength[ism-1] == 0 ) mePN = mePnAmplMapG01L1_[ism-1];
00841         if ( waveLength[ism-1] == 1 ) mePN = mePnAmplMapG01L2_[ism-1];
00842         if ( waveLength[ism-1] == 2 ) mePN = mePnAmplMapG01L3_[ism-1];
00843         if ( waveLength[ism-1] == 3 ) mePN = mePnAmplMapG01L4_[ism-1];
00844       }
00845       if ( pnItr->sample(0).gainId() == 1 ) {
00846         if ( waveLength[ism-1] == 0 ) mePN = mePnAmplMapG16L1_[ism-1];
00847         if ( waveLength[ism-1] == 1 ) mePN = mePnAmplMapG16L2_[ism-1];
00848         if ( waveLength[ism-1] == 2 ) mePN = mePnAmplMapG16L3_[ism-1];
00849         if ( waveLength[ism-1] == 3 ) mePN = mePnAmplMapG16L4_[ism-1];
00850       }
00851 
00852       if ( mePN ) mePN->Fill(num - 0.5, xvalmax);
00853 
00854       if ( num == 1 ) adcA[ism-1] = xvalmax;
00855       if ( num == 6 ) adcB[ism-1] = xvalmax;
00856 
00857     }
00858 
00859   } else {
00860 
00861     edm::LogWarning("EBLaserTask") << EcalPnDiodeDigiCollection_ << " not available";
00862 
00863   }
00864 
00865   edm::Handle<EcalUncalibratedRecHitCollection> hits;
00866 
00867   if ( e.getByLabel(EcalUncalibratedRecHitCollection_, hits) ) {
00868 
00869     int neh = hits->size();
00870     LogDebug("EBLaserTask") << "event " << ievt_ << " hits collection size " << neh;
00871 
00872     for ( EcalUncalibratedRecHitCollection::const_iterator hitItr = hits->begin(); hitItr != hits->end(); ++hitItr ) {
00873 
00874       EBDetId id = hitItr->id();
00875 
00876       int ic = id.ic();
00877       int ie = (ic-1)/20 + 1;
00878       int ip = (ic-1)%20 + 1;
00879 
00880       int ism = Numbers::iSM( id );
00881 
00882       float xie = ie - 0.5;
00883       float xip = ip - 0.5;
00884 
00885       if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
00886                runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
00887 
00888       if ( rtHalf[ism-1] != Numbers::RtHalf(id) ) continue;
00889 
00890       MonitorElement* meAmplMap = 0;
00891       MonitorElement* meTimeMap = 0;
00892       MonitorElement* meAmplPNMap = 0;
00893       MonitorElement* meAmplSummaryMap = 0;
00894 
00895       if ( rtHalf[ism-1] == 0 || rtHalf[ism-1] == 1 ) {
00896 
00897         if ( waveLength[ism-1] == 0 ) {
00898           meAmplMap = meAmplMapL1_[ism-1];
00899           meTimeMap = meTimeMapL1_[ism-1];
00900           meAmplPNMap = meAmplPNMapL1_[ism-1];
00901           meAmplSummaryMap = meAmplSummaryMapL1_;
00902         }
00903         if ( waveLength[ism-1] == 1 ) {
00904           meAmplMap = meAmplMapL2_[ism-1];
00905           meTimeMap = meTimeMapL2_[ism-1];
00906           meAmplPNMap = meAmplPNMapL2_[ism-1];
00907           meAmplSummaryMap = meAmplSummaryMapL2_;
00908         }
00909         if ( waveLength[ism-1] == 2 ) {
00910           meAmplMap = meAmplMapL3_[ism-1];
00911           meTimeMap = meTimeMapL3_[ism-1];
00912           meAmplPNMap = meAmplPNMapL3_[ism-1];
00913           meAmplSummaryMap = meAmplSummaryMapL3_;
00914         }
00915         if ( waveLength[ism-1] == 3 ) {
00916           meAmplMap = meAmplMapL4_[ism-1];
00917           meTimeMap = meTimeMapL4_[ism-1];
00918           meAmplPNMap = meAmplPNMapL4_[ism-1];
00919           meAmplSummaryMap = meAmplSummaryMapL4_;
00920         }
00921 
00922       } else {
00923 
00924         edm::LogWarning("EBLaserTask") << " RtHalf = " << rtHalf[ism-1];
00925 
00926       }
00927 
00928       float xval = hitItr->amplitude();
00929       if ( xval <= 0. ) xval = 0.0;
00930       float yval = hitItr->jitter() + 5.0;
00931       if ( yval <= 0. ) yval = 0.0;
00932       float zval = hitItr->pedestal();
00933       if ( zval <= 0. ) zval = 0.0;
00934 
00935       if ( meAmplMap ) meAmplMap->Fill(xie, xip, xval);
00936 
00937       if ( xval > 12. ) {
00938         if ( meTimeMap ) meTimeMap->Fill(xie, xip, yval);
00939       }
00940 
00941       float wval = 0.;
00942 
00943       if ( rtHalf[ism-1] == 0 ) {
00944 
00945         if ( adcA[ism-1] != 0. ) wval = xval / adcA[ism-1];
00946 
00947       } else if ( rtHalf[ism-1] == 1 ) {
00948 
00949         if ( adcB[ism-1] != 0. ) wval = xval / adcB[ism-1];
00950 
00951       } else {
00952 
00953         edm::LogWarning("EBLaserTask") << " RtHalf = " << rtHalf[ism-1];
00954 
00955       }
00956 
00957       if ( meAmplPNMap ) meAmplPNMap->Fill(xie, xip, wval);
00958 
00959       float xjp = id.iphi() - 0.5;
00960       float xje = id.ieta() - 0.5 * id.zside();
00961 
00962       if( meAmplSummaryMap ) meAmplSummaryMap->Fill(xjp, xje, xval);
00963 
00964     }
00965 
00966   } else {
00967 
00968     edm::LogWarning("EBLaserTask") << EcalUncalibratedRecHitCollection_ << " not available";
00969 
00970   }
00971 
00972 }
00973
CMSSW_4_4_3_patch1/src/DQM/EcalBarrelMonitorTasks/src/EBLaserTask.cc
