#include <EELaserTask.h>

Inheritance diagram for EELaserTask:

Public Member Functions
	EELaserTask (const edm::ParameterSet &ps)
	Constructor. More...

virtual	~EELaserTask ()
	Destructor. More...

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
	Analyze. More...

void	beginJob (void)
	BeginJob. More...

void	beginRun (const edm::Run &r, const edm::EventSetup &c)
	BeginRun. More...

void	cleanup (void)
	Cleanup. More...

void	endJob (void)
	EndJob. More...

void	endRun (const edm::Run &r, const edm::EventSetup &c)
	EndRun. More...

void	reset (void)
	Reset. More...

void	setup (void)
	Setup. More...

Protected Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Attributes
DQMStore *	dqmStore_

edm::InputTag	EcalPnDiodeDigiCollection_

edm::InputTag	EcalRawDataCollection_

edm::InputTag	EcalUncalibratedRecHitCollection_

edm::InputTag	EEDigiCollection_

bool	enableCleanup_

int	ievt_

bool	init_

std::vector< int >	laserWavelengths_

MonitorElement *	meAmplMapL1_ [18]

MonitorElement *	meAmplMapL2_ [18]

MonitorElement *	meAmplMapL3_ [18]

MonitorElement *	meAmplMapL4_ [18]

MonitorElement *	meAmplPNMapL1_ [18]

MonitorElement *	meAmplPNMapL2_ [18]

MonitorElement *	meAmplPNMapL3_ [18]

MonitorElement *	meAmplPNMapL4_ [18]

MonitorElement *	meAmplSummaryMapL1_ [2]

MonitorElement *	meAmplSummaryMapL2_ [2]

MonitorElement *	meAmplSummaryMapL3_ [2]

MonitorElement *	meAmplSummaryMapL4_ [2]

MonitorElement *	mePnAmplMapG01L1_ [18]

MonitorElement *	mePnAmplMapG01L2_ [18]

MonitorElement *	mePnAmplMapG01L3_ [18]

MonitorElement *	mePnAmplMapG01L4_ [18]

MonitorElement *	mePnAmplMapG16L1_ [18]

MonitorElement *	mePnAmplMapG16L2_ [18]

MonitorElement *	mePnAmplMapG16L3_ [18]

MonitorElement *	mePnAmplMapG16L4_ [18]

MonitorElement *	mePnPedMapG01L1_ [18]

MonitorElement *	mePnPedMapG01L2_ [18]

MonitorElement *	mePnPedMapG01L3_ [18]

MonitorElement *	mePnPedMapG01L4_ [18]

MonitorElement *	mePnPedMapG16L1_ [18]

MonitorElement *	mePnPedMapG16L2_ [18]

MonitorElement *	mePnPedMapG16L3_ [18]

MonitorElement *	mePnPedMapG16L4_ [18]

bool	mergeRuns_

MonitorElement *	meShapeMapL1_ [18]

MonitorElement *	meShapeMapL2_ [18]

MonitorElement *	meShapeMapL3_ [18]

MonitorElement *	meShapeMapL4_ [18]

MonitorElement *	meTimeMapL1_ [18]

MonitorElement *	meTimeMapL2_ [18]

MonitorElement *	meTimeMapL3_ [18]

MonitorElement *	meTimeMapL4_ [18]

std::string	prefixME_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

Definition at line 20 of file EELaserTask.h.

Constructor & Destructor Documentation

EELaserTask::EELaserTask ( const edm::ParameterSet & ps )

Constructor.

Definition at line 33 of file EELaserTask.cc.

References dqmStore_, EcalPnDiodeDigiCollection_, EcalRawDataCollection_, EcalUncalibratedRecHitCollection_, EEDigiCollection_, enableCleanup_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, init_, laserWavelengths_, meAmplMapL1_, meAmplMapL2_, meAmplMapL3_, meAmplMapL4_, meAmplPNMapL1_, meAmplPNMapL2_, meAmplPNMapL3_, meAmplPNMapL4_, meAmplSummaryMapL1_, meAmplSummaryMapL2_, meAmplSummaryMapL3_, meAmplSummaryMapL4_, mePnAmplMapG01L1_, mePnAmplMapG01L2_, mePnAmplMapG01L3_, mePnAmplMapG01L4_, mePnAmplMapG16L1_, mePnAmplMapG16L2_, mePnAmplMapG16L3_, mePnAmplMapG16L4_, mePnPedMapG01L1_, mePnPedMapG01L2_, mePnPedMapG01L3_, mePnPedMapG01L4_, mePnPedMapG16L1_, mePnPedMapG16L2_, mePnPedMapG16L3_, mePnPedMapG16L4_, mergeRuns_, meShapeMapL1_, meShapeMapL2_, meShapeMapL3_, meShapeMapL4_, meTimeMapL1_, meTimeMapL2_, meTimeMapL3_, meTimeMapL4_, cppFunctionSkipper::operator, prefixME_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                  {
 
   init_ = false;
 
   dqmStore_ = edm::Service<DQMStore>().operator->();
 
   prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
 
   enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
 
   mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
 
   EcalRawDataCollection_ = ps.getParameter<edm::InputTag>("EcalRawDataCollection");
   EEDigiCollection_ = ps.getParameter<edm::InputTag>("EEDigiCollection");
   EcalPnDiodeDigiCollection_ = ps.getParameter<edm::InputTag>("EcalPnDiodeDigiCollection");
   EcalUncalibratedRecHitCollection_ = ps.getParameter<edm::InputTag>("EcalUncalibratedRecHitCollection");
 
   // vector of enabled wavelengths (Default to all 4)
   laserWavelengths_.reserve(4);
   for ( unsigned int i = 1; i <= 4; i++ ) laserWavelengths_.push_back(i);
   laserWavelengths_ = ps.getUntrackedParameter<std::vector<int> >("laserWavelengths", laserWavelengths_);
 
   for (int i = 0; i < 18; i++) {
     meShapeMapL1_[i] = 0;
     meAmplMapL1_[i] = 0;
     meTimeMapL1_[i] = 0;
     meAmplPNMapL1_[i] = 0;
     mePnAmplMapG01L1_[i] = 0;
     mePnPedMapG01L1_[i] = 0;
     mePnAmplMapG16L1_[i] = 0;
     mePnPedMapG16L1_[i] = 0;
 
     meShapeMapL2_[i] = 0;
     meAmplMapL2_[i] = 0;
     meTimeMapL2_[i] = 0;
     meAmplPNMapL2_[i] = 0;
     mePnAmplMapG01L2_[i] = 0;
     mePnPedMapG01L2_[i] = 0;
     mePnAmplMapG16L2_[i] = 0;
     mePnPedMapG16L2_[i] = 0;
 
     meShapeMapL3_[i] = 0;
     meAmplMapL3_[i] = 0;
     meTimeMapL3_[i] = 0;
     meAmplPNMapL3_[i] = 0;
     mePnAmplMapG01L3_[i] = 0;
     mePnPedMapG01L3_[i] = 0;
     mePnAmplMapG16L3_[i] = 0;
     mePnPedMapG16L3_[i] = 0;
 
     meShapeMapL4_[i] = 0;
     meAmplMapL4_[i] = 0;
     meTimeMapL4_[i] = 0;
     meAmplPNMapL4_[i] = 0;
     mePnAmplMapG01L4_[i] = 0;
     mePnPedMapG01L4_[i] = 0;
     mePnAmplMapG16L4_[i] = 0;
     mePnPedMapG16L4_[i] = 0;
   }
 
   for(int i=0; i<2; i++){
     meAmplSummaryMapL1_[i] = 0;
     meAmplSummaryMapL2_[i] = 0;
     meAmplSummaryMapL3_[i] = 0;
     meAmplSummaryMapL4_[i] = 0;
   }
 
 }

EELaserTask::~EELaserTask ( )

virtual

Destructor.

Definition at line 102 of file EELaserTask.cc.

                          {
 
 }

Member Function Documentation

void EELaserTask::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 726 of file EELaserTask.cc.

References ecalMGPA::adc(), EcalMGPASample::adc(), EcalEndcap, EcalPnDiodeDigiCollection_, EcalRawDataCollection_, EcalUncalibratedRecHitCollection_, EEDigiCollection_, MonitorElement::Fill(), edm::Event::getByLabel(), NumbersPn::getPNs(), i, Numbers::icEE(), ievt_, init_, NumbersPn::ipnEE(), ecalpyutils::ism(), Numbers::iSM(), EEDetId::ix(), EcalDCCHeaderBlock::LASER_GAP, EcalDCCHeaderBlock::LASER_STD, LogDebug, max(), meAmplMapL1_, meAmplMapL2_, meAmplMapL3_, meAmplMapL4_, meAmplPNMapL1_, meAmplPNMapL2_, meAmplPNMapL3_, meAmplPNMapL4_, meAmplSummaryMapL1_, meAmplSummaryMapL2_, meAmplSummaryMapL3_, meAmplSummaryMapL4_, mePnAmplMapG01L1_, mePnAmplMapG01L2_, mePnAmplMapG01L3_, mePnAmplMapG01L4_, mePnAmplMapG16L1_, mePnAmplMapG16L2_, mePnAmplMapG16L3_, mePnAmplMapG16L4_, mePnPedMapG01L1_, mePnPedMapG01L2_, mePnPedMapG01L3_, mePnPedMapG01L4_, mePnPedMapG16L1_, mePnPedMapG16L2_, mePnPedMapG16L3_, mePnPedMapG16L4_, meShapeMapL1_, meShapeMapL2_, meShapeMapL3_, meShapeMapL4_, meTimeMapL1_, meTimeMapL2_, meTimeMapL3_, meTimeMapL4_, min, Numbers::RtHalf(), environment_file_cff::runType, EcalDataFrame::sample(), setup(), Numbers::subDet(), and dtDQMClient_cfg::threshold.

                                                                   {
 
   bool enable = false;
   int runType[18];
   for (int i=0; i<18; i++) runType[i] = -1;
   unsigned rtHalf[18];
   for (int i=0; i<18; i++) rtHalf[i] = -1;
   int waveLength[18];
   for (int i=0; i<18; i++) waveLength[i] = -1;
 
   edm::Handle<EcalRawDataCollection> dcchs;
 
   if ( e.getByLabel(EcalRawDataCollection_, dcchs) ) {
 
     for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
 
       if ( Numbers::subDet( *dcchItr ) != EcalEndcap ) continue;
 
       int ism = Numbers::iSM( *dcchItr, EcalEndcap );
 
       runType[ism-1] = dcchItr->getRunType();
       rtHalf[ism-1] = dcchItr->getRtHalf();
       waveLength[ism-1] = dcchItr->getEventSettings().wavelength;
 
       if ( dcchItr->getRunType() == EcalDCCHeaderBlock::LASER_STD ||
            dcchItr->getRunType() == EcalDCCHeaderBlock::LASER_GAP ) enable = true;
 
     }
 
   } else {
 
     edm::LogWarning("EELaserTask") << EcalRawDataCollection_ << " not available";
 
   }
 
   if ( ! enable ) return;
 
   if ( ! init_ ) this->setup();
 
   ievt_++;
 
   bool numPN[80];
   float adcPN[80];
   for ( int i = 0; i < 80; i++ ) {
     numPN[i] = false;
     adcPN[i] = 0.;
   }
 
   std::vector<int> PNs;
   PNs.reserve(12);
 
   edm::Handle<EEDigiCollection> digis;
 
   if ( e.getByLabel(EEDigiCollection_, digis) ) {
 
     int maxpos[10];
     for(int i(0); i < 10; i++)
       maxpos[i] = 0;
     int nReadouts(0);
 
     for ( EEDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr ) {
 
       EEDetId id = digiItr->id();
 
       int ism = Numbers::iSM( id );
 
       if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
                runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
 
       if ( rtHalf[ism-1] != Numbers::RtHalf(id) ) continue;
 
       nReadouts++;
 
       EEDataFrame dataframe = (*digiItr);
 
       int iMax(-1);
       float max(0.);
       float min(4096.);
       for (int i = 0; i < 10; i++) {
         int adc = dataframe.sample(i).adc();
     if(adc > max){
       max = adc;
       iMax = i;
     }
     if(adc < min)
       min = adc;
       }
       if(iMax >= 0 && max - min > 20.)
     maxpos[iMax] += 1;
 
     }
 
     int threshold(nReadouts / 2);
     enable = false;
     for(int i(0); i < 10; i++){
       if(maxpos[i] > threshold){
     enable = true;
     break;
       }
     }
 
     if(!enable) return;
 
     int need = digis->size();
     LogDebug("EELaserTask") << "event " << ievt_ << " digi collection size " << need;
 
     for ( EEDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr ) {
 
       EEDetId id = digiItr->id();
 
       int ix = id.ix();
       int iy = id.iy();
 
       int ism = Numbers::iSM( id );
 
       if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
                runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
 
       if ( rtHalf[ism-1] != Numbers::RtHalf(id) ) continue;
 
       int ic = Numbers::icEE(ism, ix, iy);
 
       EEDataFrame dataframe = (*digiItr);
 
       for (int i = 0; i < 10; i++) {
 
         int adc = dataframe.sample(i).adc();
 
         MonitorElement* meShapeMap = 0;
 
         if ( rtHalf[ism-1] == 0 || rtHalf[ism-1] == 1 ) {
 
           if ( waveLength[ism-1] == 0 ) meShapeMap = meShapeMapL1_[ism-1];
           if ( waveLength[ism-1] == 1 ) meShapeMap = meShapeMapL2_[ism-1];
           if ( waveLength[ism-1] == 2 ) meShapeMap = meShapeMapL3_[ism-1];
           if ( waveLength[ism-1] == 3 ) meShapeMap = meShapeMapL4_[ism-1];
 
         } else {
 
           edm::LogWarning("EELaserTask") << " RtHalf = " << rtHalf[ism-1];
 
         }
 
         float xval = float(adc);
 
         if ( meShapeMap ) meShapeMap->Fill(ic - 0.5, i + 0.5, xval);
 
       }
 
       NumbersPn::getPNs( ism, ix, iy, PNs );
 
       for (unsigned int i=0; i<PNs.size(); i++) {
         int ipn = PNs[i];
         if ( ipn >= 0 && ipn < 80 ) numPN[ipn] = true;
       }
 
     }
 
   } else {
 
     edm::LogWarning("EELaserTask") << EEDigiCollection_ << " not available";
 
   }
 
   edm::Handle<EcalPnDiodeDigiCollection> pns;
 
   if ( e.getByLabel(EcalPnDiodeDigiCollection_, pns) ) {
 
     int nep = pns->size();
     LogDebug("EELaserTask") << "event " << ievt_ << " pns collection size " << nep;
 
     for ( EcalPnDiodeDigiCollection::const_iterator pnItr = pns->begin(); pnItr != pns->end(); ++pnItr ) {
 
       if ( Numbers::subDet( pnItr->id() ) != EcalEndcap ) continue;
 
       int ism = Numbers::iSM( pnItr->id() );
 
       int num = pnItr->id().iPnId();
 
       if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
                runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
 
       int ipn = NumbersPn::ipnEE( ism, num );
 
       if ( ipn >= 0 && ipn < 80 && numPN[ipn] == false ) continue;
 
       float xvalped = 0.;
 
       for (int i = 0; i < 4; i++) {
 
         int adc = pnItr->sample(i).adc();
 
         MonitorElement* mePNPed = 0;
 
         if ( pnItr->sample(i).gainId() == 0 ) {
           if ( waveLength[ism-1] == 0 ) mePNPed = mePnPedMapG01L1_[ism-1];
           if ( waveLength[ism-1] == 1 ) mePNPed = mePnPedMapG01L2_[ism-1];
           if ( waveLength[ism-1] == 2 ) mePNPed = mePnPedMapG01L3_[ism-1];
           if ( waveLength[ism-1] == 3 ) mePNPed = mePnPedMapG01L4_[ism-1];
         }
         if ( pnItr->sample(i).gainId() == 1 ) {
           if ( waveLength[ism-1] == 0 ) mePNPed = mePnPedMapG16L1_[ism-1];
           if ( waveLength[ism-1] == 1 ) mePNPed = mePnPedMapG16L2_[ism-1];
           if ( waveLength[ism-1] == 2 ) mePNPed = mePnPedMapG16L3_[ism-1];
           if ( waveLength[ism-1] == 3 ) mePNPed = mePnPedMapG16L4_[ism-1];
         }
 
         float xval = float(adc);
 
         if ( mePNPed ) mePNPed->Fill(num - 0.5, xval);
 
         xvalped = xvalped + xval;
 
       }
 
       xvalped = xvalped / 4;
 
       float xvalmax = 0.;
 
       MonitorElement* mePN = 0;
 
       for (int i = 0; i < 50; i++) {
 
         int adc = pnItr->sample(i).adc();
 
         float xval = float(adc);
 
         if ( xval >= xvalmax ) xvalmax = xval;
 
       }
 
       xvalmax = xvalmax - xvalped;
 
       if ( pnItr->sample(0).gainId() == 0 ) {
         if ( waveLength[ism-1] == 0 ) mePN = mePnAmplMapG01L1_[ism-1];
         if ( waveLength[ism-1] == 1 ) mePN = mePnAmplMapG01L2_[ism-1];
         if ( waveLength[ism-1] == 2 ) mePN = mePnAmplMapG01L3_[ism-1];
         if ( waveLength[ism-1] == 3 ) mePN = mePnAmplMapG01L4_[ism-1];
       }
       if ( pnItr->sample(0).gainId() == 1 ) {
         if ( waveLength[ism-1] == 0 ) mePN = mePnAmplMapG16L1_[ism-1];
         if ( waveLength[ism-1] == 1 ) mePN = mePnAmplMapG16L2_[ism-1];
         if ( waveLength[ism-1] == 2 ) mePN = mePnAmplMapG16L3_[ism-1];
         if ( waveLength[ism-1] == 3 ) mePN = mePnAmplMapG16L4_[ism-1];
       }
 
       if ( mePN ) mePN->Fill(num - 0.5, xvalmax);
 
       if ( ipn >= 0 && ipn < 80 ) adcPN[ipn] = xvalmax;
 
     }
 
   } else {
 
     edm::LogWarning("EELaserTask") << EcalPnDiodeDigiCollection_ << " not available";
 
   }
 
   edm::Handle<EcalUncalibratedRecHitCollection> hits;
 
   if ( e.getByLabel(EcalUncalibratedRecHitCollection_, hits) ) {
 
     int neh = hits->size();
     LogDebug("EELaserTask") << "event " << ievt_ << " hits collection size " << neh;
 
     for ( EcalUncalibratedRecHitCollection::const_iterator hitItr = hits->begin(); hitItr != hits->end(); ++hitItr ) {
 
       EEDetId id = hitItr->id();
 
       int ix = id.ix();
       int iy = id.iy();
       int iz;
 
       int ism = Numbers::iSM( id );
 
       if ( ism >= 1 && ism <= 9 ){
     ix = 101 - ix;
     iz = 0;
       }else{
     iz = 1;
       }
 
       float xix = ix - 0.5;
       float xiy = iy - 0.5;
 
       if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
                runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
 
       if ( rtHalf[ism-1] != Numbers::RtHalf(id) ) continue;
 
       MonitorElement* meAmplMap = 0;
       MonitorElement* meTimeMap = 0;
       MonitorElement* meAmplPNMap = 0;
       MonitorElement* meAmplSummaryMap = 0;
 
       if ( rtHalf[ism-1] == 0 || rtHalf[ism-1] == 1 ) {
 
         if ( waveLength[ism-1] == 0 ) {
           meAmplMap = meAmplMapL1_[ism-1];
           meTimeMap = meTimeMapL1_[ism-1];
           meAmplPNMap = meAmplPNMapL1_[ism-1];
       meAmplSummaryMap = meAmplSummaryMapL1_[iz];
         }
         if ( waveLength[ism-1] == 1 ) {
           meAmplMap = meAmplMapL2_[ism-1];
           meTimeMap = meTimeMapL2_[ism-1];
           meAmplPNMap = meAmplPNMapL2_[ism-1];
       meAmplSummaryMap = meAmplSummaryMapL2_[iz];
         }
         if ( waveLength[ism-1] == 2 ) {
           meAmplMap = meAmplMapL3_[ism-1];
           meTimeMap = meTimeMapL3_[ism-1];
           meAmplPNMap = meAmplPNMapL3_[ism-1];
       meAmplSummaryMap = meAmplSummaryMapL3_[iz];
         }
         if ( waveLength[ism-1] == 3 ) {
           meAmplMap = meAmplMapL4_[ism-1];
           meTimeMap = meTimeMapL4_[ism-1];
           meAmplPNMap = meAmplPNMapL4_[ism-1];
       meAmplSummaryMap = meAmplSummaryMapL4_[iz];
         }
 
       } else {
 
         edm::LogWarning("EELaserTask") << " RtHalf = " << rtHalf[ism-1];
 
       }
 
       float xval = hitItr->amplitude();
       if ( xval <= 0. ) xval = 0.0;
       float yval = hitItr->jitter() + 5.0;
       if ( yval <= 0. ) yval = 0.0;
       float zval = hitItr->pedestal();
       if ( zval <= 0. ) zval = 0.0;
 
       if ( meAmplMap ) meAmplMap->Fill(xix, xiy, xval);
 
       if ( xval > 16. ) {
         if ( meTimeMap ) meTimeMap->Fill(xix, xiy, yval);
       }
 
       float wval = 0.;
 
       NumbersPn::getPNs( ism, ix, iy, PNs );
 
       if ( PNs.size() > 0 ) {
         int ipn = PNs[0];
         if ( ipn >= 0 && ipn < 80 ) {
           if ( adcPN[ipn] != 0. ) wval = xval / adcPN[ipn];
         }
       }
 
       if ( meAmplPNMap ) meAmplPNMap->Fill(xix, xiy, wval);
 
       if( meAmplSummaryMap ) meAmplSummaryMap->Fill(xix, xiy, xval);
 
     }
 
   } else {
 
     edm::LogWarning("EELaserTask") << EcalUncalibratedRecHitCollection_ << " not available";
 
   }
 
 }

void EELaserTask::beginJob ( void )

protectedvirtual

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 106 of file EELaserTask.cc.

References dqmStore_, ievt_, prefixME_, DQMStore::rmdir(), and DQMStore::setCurrentFolder().

                               {
 
   ievt_ = 0;
 
   if ( dqmStore_ ) {
     dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask");
     dqmStore_->rmdir(prefixME_ + "/EELaserTask");
   }
 
 }

void EELaserTask::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

protectedvirtual

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 117 of file EELaserTask.cc.

References Numbers::initGeometry(), mergeRuns_, and reset().

                                                                   {
 
   Numbers::initGeometry(c, false);
 
   if ( ! mergeRuns_ ) this->reset();
 
 }

void EELaserTask::cleanup ( void )

protected

Cleanup.

Definition at line 542 of file EELaserTask.cc.

References dqmStore_, spr::find(), i, init_, laserWavelengths_, meAmplMapL1_, meAmplMapL2_, meAmplMapL3_, meAmplMapL4_, meAmplPNMapL1_, meAmplPNMapL2_, meAmplPNMapL3_, meAmplPNMapL4_, meAmplSummaryMapL1_, meAmplSummaryMapL2_, meAmplSummaryMapL3_, meAmplSummaryMapL4_, mePnAmplMapG01L1_, mePnAmplMapG01L2_, mePnAmplMapG01L3_, mePnAmplMapG01L4_, mePnAmplMapG16L1_, mePnAmplMapG16L2_, mePnAmplMapG16L3_, mePnAmplMapG16L4_, mePnPedMapG01L1_, mePnPedMapG01L2_, mePnPedMapG01L3_, mePnPedMapG01L4_, mePnPedMapG16L1_, mePnPedMapG16L2_, mePnPedMapG16L3_, mePnPedMapG16L4_, meShapeMapL1_, meShapeMapL2_, meShapeMapL3_, meShapeMapL4_, meTimeMapL1_, meTimeMapL2_, meTimeMapL3_, meTimeMapL4_, prefixME_, DQMStore::removeElement(), and DQMStore::setCurrentFolder().

Referenced by endJob().

                              {
 
   if ( ! init_ ) return;
 
   if ( dqmStore_ ) {
     dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask");
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser1");
       for (int i = 0; i < 18; i++) {
         if ( meShapeMapL1_[i] )  dqmStore_->removeElement( meShapeMapL1_[i]->getName() );
         meShapeMapL1_[i] = 0;
         if ( meAmplMapL1_[i] ) dqmStore_->removeElement( meAmplMapL1_[i]->getName() );
         meAmplMapL1_[i] = 0;
         if ( meTimeMapL1_[i] ) dqmStore_->removeElement( meTimeMapL1_[i]->getName() );
         meTimeMapL1_[i] = 0;
         if ( meAmplPNMapL1_[i] ) dqmStore_->removeElement( meAmplPNMapL1_[i]->getName() );
         meAmplPNMapL1_[i] = 0;
       }
 
       for(int i=0; i<2; i++)
     if( meAmplSummaryMapL1_[i] ) dqmStore_->removeElement( meAmplSummaryMapL1_[i]->getName() );
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser2");
       for (int i = 0; i < 18; i++) {
         if ( meShapeMapL2_[i] )  dqmStore_->removeElement( meShapeMapL2_[i]->getName() );
         meShapeMapL2_[i] = 0;
         if ( meAmplMapL2_[i] ) dqmStore_->removeElement( meAmplMapL2_[i]->getName() );
         meAmplMapL2_[i] = 0;
         if ( meTimeMapL2_[i] ) dqmStore_->removeElement( meTimeMapL2_[i]->getName() );
         meTimeMapL2_[i] = 0;
         if ( meAmplPNMapL2_[i] ) dqmStore_->removeElement( meAmplPNMapL2_[i]->getName() );
         meAmplPNMapL2_[i] = 0;
       }
 
       for(int i=0; i<2; i++)
     if( meAmplSummaryMapL2_[i] ) dqmStore_->removeElement( meAmplSummaryMapL2_[i]->getName() );
 
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser3");
       for (int i = 0; i < 18; i++) {
         if ( meShapeMapL3_[i] )  dqmStore_->removeElement( meShapeMapL3_[i]->getName() );
         meShapeMapL3_[i] = 0;
         if ( meAmplMapL3_[i] ) dqmStore_->removeElement( meAmplMapL3_[i]->getName() );
         meAmplMapL3_[i] = 0;
         if ( meTimeMapL3_[i] ) dqmStore_->removeElement( meTimeMapL3_[i]->getName() );
         meTimeMapL3_[i] = 0;
         if ( meAmplPNMapL3_[i] ) dqmStore_->removeElement( meAmplPNMapL3_[i]->getName() );
         meAmplPNMapL3_[i] = 0;
       }
 
       for(int i=0; i<2; i++)
     if( meAmplSummaryMapL3_[i] ) dqmStore_->removeElement( meAmplSummaryMapL3_[i]->getName() );
 
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser4");
       for (int i = 0; i < 18; i++) {
         if ( meShapeMapL4_[i] )  dqmStore_->removeElement( meShapeMapL4_[i]->getName() );
         meShapeMapL4_[i] = 0;
         if ( meAmplMapL4_[i] ) dqmStore_->removeElement( meAmplMapL4_[i]->getName() );
         meAmplMapL4_[i] = 0;
         if ( meTimeMapL4_[i] ) dqmStore_->removeElement( meTimeMapL4_[i]->getName() );
         meTimeMapL4_[i] = 0;
         if ( meAmplPNMapL4_[i] ) dqmStore_->removeElement( meAmplPNMapL4_[i]->getName() );
         meAmplPNMapL4_[i] = 0;
       }
 
       for(int i=0; i<2; i++)
     if( meAmplSummaryMapL4_[i] ) dqmStore_->removeElement( meAmplSummaryMapL4_[i]->getName() );
 
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser1/PN");
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser1/PN/Gain01");
       for (int i = 0; i < 18; i++) {
         if ( mePnAmplMapG01L1_[i] ) dqmStore_->removeElement( mePnAmplMapG01L1_[i]->getName() );
         mePnAmplMapG01L1_[i] = 0;
         if ( mePnPedMapG01L1_[i] ) dqmStore_->removeElement( mePnPedMapG01L1_[i]->getName() );
         mePnPedMapG01L1_[i] = 0;
       }
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser1/PN/Gain16");
       for (int i = 0; i < 18; i++) {
         if ( mePnAmplMapG16L1_[i] ) dqmStore_->removeElement( mePnAmplMapG16L1_[i]->getName() );
         mePnAmplMapG16L1_[i] = 0;
         if ( mePnPedMapG16L1_[i] ) dqmStore_->removeElement( mePnPedMapG16L1_[i]->getName() );
         mePnPedMapG16L1_[i] = 0;
       }
 
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser2/PN");
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser2/PN/Gain01");
       for (int i = 0; i < 18; i++) {
         if ( mePnAmplMapG01L2_[i] ) dqmStore_->removeElement( mePnAmplMapG01L2_[i]->getName() );
         mePnAmplMapG01L2_[i] = 0;
         if ( mePnPedMapG01L2_[i] ) dqmStore_->removeElement( mePnPedMapG01L2_[i]->getName() );
         mePnPedMapG01L2_[i] = 0;
       }
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser2/PN/Gain16");
       for (int i = 0; i < 18; i++) {
         if ( mePnAmplMapG16L2_[i] ) dqmStore_->removeElement( mePnAmplMapG16L2_[i]->getName() );
         mePnAmplMapG16L2_[i] = 0;
         if ( mePnPedMapG16L2_[i] ) dqmStore_->removeElement( mePnPedMapG16L2_[i]->getName() );
         mePnPedMapG16L2_[i] = 0;
       }
 
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser3/PN");
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser3/PN/Gain01");
       for (int i = 0; i < 18; i++) {
         if ( mePnAmplMapG01L3_[i] ) dqmStore_->removeElement( mePnAmplMapG01L3_[i]->getName() );
         mePnAmplMapG01L3_[i] = 0;
         if ( mePnPedMapG01L3_[i] ) dqmStore_->removeElement( mePnPedMapG01L3_[i]->getName() );
         mePnPedMapG01L3_[i] = 0;
       }
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser3/PN/Gain16");
       for (int i = 0; i < 18; i++) {
         if ( mePnAmplMapG16L3_[i] ) dqmStore_->removeElement( mePnAmplMapG16L3_[i]->getName() );
         mePnAmplMapG16L3_[i] = 0;
         if ( mePnPedMapG16L3_[i] ) dqmStore_->removeElement( mePnPedMapG16L3_[i]->getName() );
         mePnPedMapG16L3_[i] = 0;
       }
 
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser4/PN");
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser4/PN/Gain01");
       for (int i = 0; i < 18; i++) {
         if ( mePnAmplMapG01L4_[i] ) dqmStore_->removeElement( mePnAmplMapG01L4_[i]->getName() );
         mePnAmplMapG01L4_[i] = 0;
         if ( mePnPedMapG01L4_[i] ) dqmStore_->removeElement( mePnPedMapG01L4_[i]->getName() );
         mePnPedMapG01L4_[i] = 0;
       }
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser4/PN/Gain16");
       for (int i = 0; i < 18; i++) {
         if ( mePnAmplMapG16L4_[i] ) dqmStore_->removeElement( mePnAmplMapG16L4_[i]->getName() );
         mePnAmplMapG16L4_[i] = 0;
         if ( mePnPedMapG16L4_[i] ) dqmStore_->removeElement( mePnPedMapG16L4_[i]->getName() );
         mePnPedMapG16L4_[i] = 0;
       }
 
     }
 
   }
 
   init_ = false;
 
 }

void EELaserTask::endJob ( void )

protectedvirtual

EndJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 718 of file EELaserTask.cc.

References cleanup(), enableCleanup_, and ievt_.

                             {
 
   edm::LogInfo("EELaserTask") << "analyzed " << ievt_ << " events";
 
   if ( enableCleanup_ ) this->cleanup();
 
 }

void EELaserTask::endRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

protectedvirtual

EndRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 125 of file EELaserTask.cc.

                                                                 {
 
 }

void EELaserTask::reset ( void )

protected

Reset.

Definition at line 129 of file EELaserTask.cc.

References spr::find(), i, laserWavelengths_, meAmplMapL1_, meAmplMapL2_, meAmplMapL3_, meAmplMapL4_, meAmplPNMapL1_, meAmplPNMapL2_, meAmplPNMapL3_, meAmplPNMapL4_, meAmplSummaryMapL1_, meAmplSummaryMapL2_, meAmplSummaryMapL3_, meAmplSummaryMapL4_, mePnAmplMapG01L1_, mePnAmplMapG01L2_, mePnAmplMapG01L3_, mePnAmplMapG01L4_, mePnAmplMapG16L1_, mePnAmplMapG16L2_, mePnAmplMapG16L3_, mePnAmplMapG16L4_, mePnPedMapG01L1_, mePnPedMapG01L2_, mePnPedMapG01L3_, mePnPedMapG01L4_, mePnPedMapG16L1_, mePnPedMapG16L2_, mePnPedMapG16L3_, mePnPedMapG16L4_, meShapeMapL1_, meShapeMapL2_, meShapeMapL3_, meShapeMapL4_, meTimeMapL1_, meTimeMapL2_, meTimeMapL3_, meTimeMapL4_, and MonitorElement::Reset().

Referenced by beginRun().

                             {
 
   for (int i = 0; i < 18; i++) {
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
       if ( meShapeMapL1_[i] )  meShapeMapL1_[i]->Reset();
       if ( meAmplMapL1_[i] ) meAmplMapL1_[i]->Reset();
       if ( meTimeMapL1_[i] ) meTimeMapL1_[i]->Reset();
       if ( meAmplPNMapL1_[i] ) meAmplPNMapL1_[i]->Reset();
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
       if ( meShapeMapL2_[i] )  meShapeMapL2_[i]->Reset();
       if ( meAmplMapL2_[i] ) meAmplMapL2_[i]->Reset();
       if ( meTimeMapL2_[i] ) meTimeMapL2_[i]->Reset();
       if ( meAmplPNMapL2_[i] ) meAmplPNMapL2_[i]->Reset();
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
       if ( meShapeMapL3_[i] )  meShapeMapL3_[i]->Reset();
       if ( meAmplMapL3_[i] ) meAmplMapL3_[i]->Reset();
       if ( meTimeMapL3_[i] ) meTimeMapL3_[i]->Reset();
       if ( meAmplPNMapL3_[i] ) meAmplPNMapL3_[i]->Reset();
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
       if ( meShapeMapL4_[i] )  meShapeMapL4_[i]->Reset();
       if ( meAmplMapL4_[i] ) meAmplMapL4_[i]->Reset();
       if ( meTimeMapL4_[i] ) meTimeMapL4_[i]->Reset();
       if ( meAmplPNMapL4_[i] ) meAmplPNMapL4_[i]->Reset();
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
       if ( mePnAmplMapG01L1_[i] ) mePnAmplMapG01L1_[i]->Reset();
       if ( mePnPedMapG01L1_[i] ) mePnPedMapG01L1_[i]->Reset();
 
       if ( mePnAmplMapG16L1_[i] ) mePnAmplMapG16L1_[i]->Reset();
       if ( mePnPedMapG16L1_[i] ) mePnPedMapG16L1_[i]->Reset();
     }
 
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
       if ( mePnAmplMapG01L2_[i] ) mePnAmplMapG01L2_[i]->Reset();
       if ( mePnPedMapG01L2_[i] ) mePnPedMapG01L2_[i]->Reset();
 
       if ( mePnAmplMapG16L2_[i] ) mePnAmplMapG16L2_[i]->Reset();
       if ( mePnPedMapG16L2_[i] ) mePnPedMapG16L2_[i]->Reset();
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
       if ( mePnAmplMapG01L3_[i] ) mePnAmplMapG01L3_[i]->Reset();
       if ( mePnPedMapG01L3_[i] ) mePnPedMapG01L3_[i]->Reset();
 
       if ( mePnAmplMapG16L3_[i] ) mePnAmplMapG16L3_[i]->Reset();
       if ( mePnPedMapG16L3_[i] ) mePnPedMapG16L3_[i]->Reset();
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
       if ( mePnAmplMapG01L4_[i] ) mePnAmplMapG01L4_[i]->Reset();
       if ( mePnPedMapG01L4_[i] ) mePnPedMapG01L4_[i]->Reset();
 
       if ( mePnAmplMapG16L4_[i] ) mePnAmplMapG16L4_[i]->Reset();
       if ( mePnPedMapG16L4_[i] ) mePnPedMapG16L4_[i]->Reset();
     }
 
   }
 
   for(int i=0; i<2; i++){
     if( meAmplSummaryMapL1_[i] ) meAmplSummaryMapL1_[i]->Reset();
     if( meAmplSummaryMapL2_[i] ) meAmplSummaryMapL2_[i]->Reset();
     if( meAmplSummaryMapL3_[i] ) meAmplSummaryMapL3_[i]->Reset();
     if( meAmplSummaryMapL4_[i] ) meAmplSummaryMapL4_[i]->Reset();
   }
 
 }

void EELaserTask::setup ( void )

protected

Setup.

Definition at line 204 of file EELaserTask.cc.

References DQMStore::bookProfile(), DQMStore::bookProfile2D(), dqmStore_, spr::find(), i, init_, Numbers::ix0EE(), Numbers::iy0EE(), laserWavelengths_, meAmplMapL1_, meAmplMapL2_, meAmplMapL3_, meAmplMapL4_, meAmplPNMapL1_, meAmplPNMapL2_, meAmplPNMapL3_, meAmplPNMapL4_, meAmplSummaryMapL1_, meAmplSummaryMapL2_, meAmplSummaryMapL3_, meAmplSummaryMapL4_, mePnAmplMapG01L1_, mePnAmplMapG01L2_, mePnAmplMapG01L3_, mePnAmplMapG01L4_, mePnAmplMapG16L1_, mePnAmplMapG16L2_, mePnAmplMapG16L3_, mePnAmplMapG16L4_, mePnPedMapG01L1_, mePnPedMapG01L2_, mePnPedMapG01L3_, mePnPedMapG01L4_, mePnPedMapG16L1_, mePnPedMapG16L2_, mePnPedMapG16L3_, mePnPedMapG16L4_, meShapeMapL1_, meShapeMapL2_, meShapeMapL3_, meShapeMapL4_, meTimeMapL1_, meTimeMapL2_, meTimeMapL3_, meTimeMapL4_, mergeVDriftHistosByStation::name, prefixME_, Numbers::sEE(), MonitorElement::setAxisTitle(), DQMStore::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and DQMStore::tag().

Referenced by analyze().

                            {
 
   init_ = true;
 
   std::string name;
   std::stringstream LaserN, LN;
 
   if ( dqmStore_ ) {
     dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask");
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
 
       LaserN.str("");
       LaserN << "Laser" << 1;
       LN.str("");
       LN << "L" << 1;
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str());
       for (int i = 0; i < 18; i++) {
     name = "EELT shape " + Numbers::sEE(i+1) + " " + LN.str(); 
         meShapeMapL1_[i] = dqmStore_->bookProfile2D(name, name, 850, 0., 850., 10, 0., 10., 4096, 0., 4096., "s");
         meShapeMapL1_[i]->setAxisTitle("channel", 1);
         meShapeMapL1_[i]->setAxisTitle("sample", 2);
         meShapeMapL1_[i]->setAxisTitle("amplitude", 3);
         dqmStore_->tag(meShapeMapL1_[i], i+1);
 
     name = "EELT amplitude " + Numbers::sEE(i+1) + " " + LN.str(); 
         meAmplMapL1_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
         meAmplMapL1_[i]->setAxisTitle("ix", 1);
         if ( i+1 >= 1 && i+1 <= 9 ) meAmplMapL1_[i]->setAxisTitle("101-ix", 1);
         meAmplMapL1_[i]->setAxisTitle("iy", 2);
         dqmStore_->tag(meAmplMapL1_[i], i+1);
 
     name = "EELT timing " + Numbers::sEE(i+1) + " " + LN.str(); 
         meTimeMapL1_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 250, 0., 10., "s");
         meTimeMapL1_[i]->setAxisTitle("ix", 1);
         if ( i+1 >= 1 && i+1 <= 9 ) meTimeMapL1_[i]->setAxisTitle("101-ix", 1);
         meTimeMapL1_[i]->setAxisTitle("iy", 2);
         dqmStore_->tag(meTimeMapL1_[i], i+1);
 
     name = "EELT amplitude over PN " + Numbers::sEE(i+1) + " " + LN.str();
         meAmplPNMapL1_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
         meAmplPNMapL1_[i]->setAxisTitle("ix", 1);
         if ( i+1 >= 1 && i+1 <= 9 ) meAmplPNMapL1_[i]->setAxisTitle("101-ix", 1);
         meAmplPNMapL1_[i]->setAxisTitle("iy", 2);
         dqmStore_->tag(meAmplPNMapL1_[i], i+1);
       }
 
       name = "EELT amplitude map " + LN.str() + " EE -";
       meAmplSummaryMapL1_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
       meAmplSummaryMapL1_[0]->setAxisTitle("ix", 1);
       meAmplSummaryMapL1_[0]->setAxisTitle("iy", 2);
 
       name = "EELT amplitude map " + LN.str() + " EE +";
       meAmplSummaryMapL1_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
       meAmplSummaryMapL1_[1]->setAxisTitle("ix", 1);
       meAmplSummaryMapL1_[1]->setAxisTitle("iy", 2);
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN");
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain01");
 
       for (int i = 0; i < 18; i++) {
     name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G01 " + LN.str(); 
         mePnAmplMapG01L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnAmplMapG01L1_[i]->setAxisTitle("channel", 1);
         mePnAmplMapG01L1_[i]->setAxisTitle("amplitude", 2);
         dqmStore_->tag(mePnAmplMapG01L1_[i], i+1);
 
     name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G01 " + LN.str(); 
         mePnPedMapG01L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnPedMapG01L1_[i]->setAxisTitle("channel", 1);
         mePnPedMapG01L1_[i]->setAxisTitle("pedestal", 2);
         dqmStore_->tag(mePnPedMapG01L1_[i], i+1);
       }
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain16");
 
       for (int i = 0; i < 18; i++) {
     name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G16 " + LN.str();
         mePnAmplMapG16L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnAmplMapG16L1_[i]->setAxisTitle("channel", 1);
         mePnAmplMapG16L1_[i]->setAxisTitle("amplitude", 2);
         dqmStore_->tag(mePnAmplMapG16L1_[i], i+1);
 
     name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G16 " + LN.str();
         mePnPedMapG16L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnPedMapG16L1_[i]->setAxisTitle("channel", 1);
         mePnPedMapG16L1_[i]->setAxisTitle("pedestal", 2);
         dqmStore_->tag(mePnPedMapG16L1_[i], i+1);
       }
 
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
 
       LaserN.str("");
       LaserN << "Laser" << 2;
       LN.str("");
       LN << "L" << 2;
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str());
       for (int i = 0; i < 18; i++) {
     name = "EELT shape " + Numbers::sEE(i+1) + " " + LN.str(); 
         meShapeMapL2_[i] = dqmStore_->bookProfile2D(name, name, 850, 0., 850., 10, 0., 10., 4096, 0., 4096., "s");
         meShapeMapL2_[i]->setAxisTitle("channel", 1);
         meShapeMapL2_[i]->setAxisTitle("sample", 2);
         meShapeMapL2_[i]->setAxisTitle("amplitude", 3);
         dqmStore_->tag(meShapeMapL2_[i], i+1);
 
     name = "EELT amplitude " + Numbers::sEE(i+1) + " " + LN.str(); 
         meAmplMapL2_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
         meAmplMapL2_[i]->setAxisTitle("ix", 1);
         if ( i+1 >= 1 && i+1 <= 9 ) meAmplMapL2_[i]->setAxisTitle("101-ix", 1);
         meAmplMapL2_[i]->setAxisTitle("iy", 2);
         dqmStore_->tag(meAmplMapL2_[i], i+1);
 
     name = "EELT timing " + Numbers::sEE(i+1) + " " + LN.str(); 
         meTimeMapL2_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 250, 0., 10., "s");
         meTimeMapL2_[i]->setAxisTitle("ix", 1);
         if ( i+1 >= 1 && i+1 <= 9 ) meTimeMapL2_[i]->setAxisTitle("101-ix", 1);
         meTimeMapL2_[i]->setAxisTitle("iy", 2);
         dqmStore_->tag(meTimeMapL2_[i], i+1);
 
     name = "EELT amplitude over PN " + Numbers::sEE(i+1) + " " + LN.str();
         meAmplPNMapL2_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
         meAmplPNMapL2_[i]->setAxisTitle("ix", 1);
         if ( i+1 >= 1 && i+1 <= 9 ) meAmplPNMapL2_[i]->setAxisTitle("101-ix", 1);
         meAmplPNMapL2_[i]->setAxisTitle("iy", 2);
         dqmStore_->tag(meAmplPNMapL2_[i], i+1);
       }
 
       name = "EELT amplitude map " + LN.str() + " EE -";
       meAmplSummaryMapL2_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
       meAmplSummaryMapL2_[0]->setAxisTitle("ix", 1);
       meAmplSummaryMapL2_[0]->setAxisTitle("iy", 2);
 
       name = "EELT amplitude map " + LN.str() + " EE +";
       meAmplSummaryMapL2_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
       meAmplSummaryMapL2_[1]->setAxisTitle("ix", 1);
       meAmplSummaryMapL2_[1]->setAxisTitle("iy", 2);
 
       for (int i = 0; i < 18; i++) {
     name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G01 " + LN.str(); 
         mePnAmplMapG01L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnAmplMapG01L2_[i]->setAxisTitle("channel", 1);
         mePnAmplMapG01L2_[i]->setAxisTitle("amplitude", 2);
         dqmStore_->tag(mePnAmplMapG01L2_[i], i+1);
 
     name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G01 " + LN.str(); 
         mePnPedMapG01L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnPedMapG01L2_[i]->setAxisTitle("channel", 1);
         mePnPedMapG01L2_[i]->setAxisTitle("pedestal", 2);
         dqmStore_->tag(mePnPedMapG01L2_[i], i+1);
       }
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain16");
 
       for (int i = 0; i < 18; i++) {
     name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G16 " + LN.str();
         mePnAmplMapG16L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnAmplMapG16L2_[i]->setAxisTitle("channel", 1);
         mePnAmplMapG16L2_[i]->setAxisTitle("amplitude", 2);
         dqmStore_->tag(mePnAmplMapG16L2_[i], i+1);
 
     name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G16 " + LN.str();
         mePnPedMapG16L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnPedMapG16L2_[i]->setAxisTitle("channel", 1);
         mePnPedMapG16L2_[i]->setAxisTitle("pedestal", 2);
         dqmStore_->tag(mePnPedMapG16L2_[i], i+1);
       }
 
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
 
       LaserN.str("");
       LaserN << "Laser" << 3;
       LN.str("");
       LN << "L" << 3;
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str());
       for (int i = 0; i < 18; i++) {
     name = "EELT shape " + Numbers::sEE(i+1) + " " + LN.str(); 
         meShapeMapL3_[i] = dqmStore_->bookProfile2D(name, name, 850, 0., 850., 10, 0., 10., 4096, 0., 4096., "s");
         meShapeMapL3_[i]->setAxisTitle("channel", 1);
         meShapeMapL3_[i]->setAxisTitle("sample", 2);
         meShapeMapL3_[i]->setAxisTitle("amplitude", 3);
         dqmStore_->tag(meShapeMapL3_[i], i+1);
 
     name = "EELT amplitude " + Numbers::sEE(i+1) + " " + LN.str(); 
         meAmplMapL3_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
         meAmplMapL3_[i]->setAxisTitle("ix", 1);
         if ( i+1 >= 1 && i+1 <= 9 ) meAmplMapL3_[i]->setAxisTitle("101-ix", 1);
         meAmplMapL3_[i]->setAxisTitle("iy", 2);
         dqmStore_->tag(meAmplMapL3_[i], i+1);
 
     name = "EELT timing " + Numbers::sEE(i+1) + " " + LN.str(); 
         meTimeMapL3_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 250, 0., 10., "s");
         meTimeMapL3_[i]->setAxisTitle("ix", 1);
         if ( i+1 >= 1 && i+1 <= 9 ) meTimeMapL3_[i]->setAxisTitle("101-ix", 1);
         meTimeMapL3_[i]->setAxisTitle("iy", 2);
         dqmStore_->tag(meTimeMapL3_[i], i+1);
 
     name = "EELT amplitude over PN " + Numbers::sEE(i+1) + " " + LN.str();
         meAmplPNMapL3_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
         meAmplPNMapL3_[i]->setAxisTitle("ix", 1);
         if ( i+1 >= 1 && i+1 <= 9 ) meAmplPNMapL3_[i]->setAxisTitle("101-ix", 1);
         meAmplPNMapL3_[i]->setAxisTitle("iy", 2);
         dqmStore_->tag(meAmplPNMapL3_[i], i+1);
       }
 
       name = "EELT amplitude map " + LN.str() + " EE -";
       meAmplSummaryMapL3_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
       meAmplSummaryMapL3_[0]->setAxisTitle("ix", 1);
       meAmplSummaryMapL3_[0]->setAxisTitle("iy", 2);
 
       name = "EELT amplitude map " + LN.str() + " EE +";
       meAmplSummaryMapL3_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
       meAmplSummaryMapL3_[1]->setAxisTitle("ix", 1);
       meAmplSummaryMapL3_[1]->setAxisTitle("iy", 2);
 
       for (int i = 0; i < 18; i++) {
     name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G01 " + LN.str(); 
         mePnAmplMapG01L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnAmplMapG01L3_[i]->setAxisTitle("channel", 1);
         mePnAmplMapG01L3_[i]->setAxisTitle("amplitude", 2);
         dqmStore_->tag(mePnAmplMapG01L3_[i], i+1);
 
     name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G01 " + LN.str(); 
         mePnPedMapG01L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnPedMapG01L3_[i]->setAxisTitle("channel", 1);
         mePnPedMapG01L3_[i]->setAxisTitle("pedestal", 2);
         dqmStore_->tag(mePnPedMapG01L3_[i], i+1);
       }
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain16");
 
       for (int i = 0; i < 18; i++) {
     name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G16 " + LN.str();
         mePnAmplMapG16L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnAmplMapG16L3_[i]->setAxisTitle("channel", 1);
         mePnAmplMapG16L3_[i]->setAxisTitle("amplitude", 2);
         dqmStore_->tag(mePnAmplMapG16L3_[i], i+1);
 
     name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G16 " + LN.str();
         mePnPedMapG16L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnPedMapG16L3_[i]->setAxisTitle("channel", 1);
         mePnPedMapG16L3_[i]->setAxisTitle("pedestal", 2);
         dqmStore_->tag(mePnPedMapG16L3_[i], i+1);
       }
 
     }
 
     if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
 
       LaserN.str("");
       LaserN << "Laser" << 4;
       LN.str("");
       LN << "L" << 4;
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str());
       for (int i = 0; i < 18; i++) {
     name = "EELT shape " + Numbers::sEE(i+1) + " " + LN.str(); 
         meShapeMapL4_[i] = dqmStore_->bookProfile2D(name, name, 850, 0., 850., 10, 0., 10., 4096, 0., 4096., "s");
         meShapeMapL4_[i]->setAxisTitle("channel", 1);
         meShapeMapL4_[i]->setAxisTitle("sample", 2);
         meShapeMapL4_[i]->setAxisTitle("amplitude", 3);
         dqmStore_->tag(meShapeMapL4_[i], i+1);
 
     name = "EELT amplitude " + Numbers::sEE(i+1) + " " + LN.str(); 
         meAmplMapL4_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
         meAmplMapL4_[i]->setAxisTitle("ix", 1);
         if ( i+1 >= 1 && i+1 <= 9 ) meAmplMapL4_[i]->setAxisTitle("101-ix", 1);
         meAmplMapL4_[i]->setAxisTitle("iy", 2);
         dqmStore_->tag(meAmplMapL4_[i], i+1);
 
     name = "EELT timing " + Numbers::sEE(i+1) + " " + LN.str(); 
         meTimeMapL4_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 250, 0., 10., "s");
         meTimeMapL4_[i]->setAxisTitle("ix", 1);
         if ( i+1 >= 1 && i+1 <= 9 ) meTimeMapL4_[i]->setAxisTitle("101-ix", 1);
         meTimeMapL4_[i]->setAxisTitle("iy", 2);
         dqmStore_->tag(meTimeMapL4_[i], i+1);
 
     name = "EELT amplitude over PN " + Numbers::sEE(i+1) + " " + LN.str();
         meAmplPNMapL4_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
         meAmplPNMapL4_[i]->setAxisTitle("ix", 1);
         if ( i+1 >= 1 && i+1 <= 9 ) meAmplPNMapL4_[i]->setAxisTitle("101-ix", 1);
         meAmplPNMapL4_[i]->setAxisTitle("iy", 2);
         dqmStore_->tag(meAmplPNMapL4_[i], i+1);
       }
 
       name = "EELT amplitude map " + LN.str() + " EE -";
       meAmplSummaryMapL4_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
       meAmplSummaryMapL4_[0]->setAxisTitle("ix", 1);
       meAmplSummaryMapL4_[0]->setAxisTitle("iy", 2);
 
       name = "EELT amplitude map " + LN.str() + " EE +";
       meAmplSummaryMapL4_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
       meAmplSummaryMapL4_[1]->setAxisTitle("ix", 1);
       meAmplSummaryMapL4_[1]->setAxisTitle("iy", 2);
 
       for (int i = 0; i < 18; i++) {
     name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G01 " + LN.str(); 
         mePnAmplMapG01L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnAmplMapG01L4_[i]->setAxisTitle("channel", 1);
         mePnAmplMapG01L4_[i]->setAxisTitle("amplitude", 2);
         dqmStore_->tag(mePnAmplMapG01L4_[i], i+1);
 
     name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G01 " + LN.str(); 
         mePnPedMapG01L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnPedMapG01L4_[i]->setAxisTitle("channel", 1);
         mePnPedMapG01L4_[i]->setAxisTitle("pedestal", 2);
         dqmStore_->tag(mePnPedMapG01L4_[i], i+1);
       }
 
       dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain16");
 
       for (int i = 0; i < 18; i++) {
     name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G16 " + LN.str();
         mePnAmplMapG16L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnAmplMapG16L4_[i]->setAxisTitle("channel", 1);
         mePnAmplMapG16L4_[i]->setAxisTitle("amplitude", 2);
         dqmStore_->tag(mePnAmplMapG16L4_[i], i+1);
 
     name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G16 " + LN.str();
         mePnPedMapG16L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
         mePnPedMapG16L4_[i]->setAxisTitle("channel", 1);
         mePnPedMapG16L4_[i]->setAxisTitle("pedestal", 2);
         dqmStore_->tag(mePnPedMapG16L4_[i], i+1);
       }
 
     }
 
   }
 
 }