#include <LedTask.h>

Inheritance diagram for ecaldqm::LedTask:

Public Member Functions
void	addDependencies (DependencySet &) override

bool	analyze (void const *, Collections) override

void	beginEvent (edm::Event const &, edm::EventSetup const &, bool const &, bool &) override

void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	beginRun (edm::Run const &, edm::EventSetup const &) override

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

bool	filterRunType (short const *) override

	LedTask ()

void	runOnDigis (EEDigiCollection const &)

void	runOnPnDigis (EcalPnDiodeDigiCollection const &)

void	runOnRawData (EcalRawDataCollection const &)

void	runOnUncalibRecHits (EcalUncalibratedRecHitCollection const &)

	~LedTask () override

Public Member Functions inherited from ecaldqm::DQWorkerTask
	DQWorkerTask ()

virtual void	endEvent (edm::Event const &, edm::EventSetup const &)

virtual bool	filterTrigger (edm::TriggerResultsByName const &)

virtual void	setTokens (edm::ConsumesCollector &)

	~DQWorkerTask () override

Public Member Functions inherited from ecaldqm::DQWorker
virtual void	bookMEs (DQMStore::IBooker &)

	DQWorker ()

virtual void	endRun (edm::Run const &, edm::EventSetup const &)

EcalDQMSetupObjects const	getEcalDQMSetupObjects ()

EcalElectronicsMapping const *	GetElectronicsMap ()

CaloGeometry const *	GetGeometry ()

std::string const &	getName () const

CaloTopology const *	GetTopology ()

EcalTrigTowerConstituentsMap const *	GetTrigTowerMap ()

bool	onlineMode () const

virtual void	releaseMEs ()

void	setEventNumber (edm::EventNumber_t _e)

void	setLumiNumber (edm::LuminosityBlockNumber_t _l)

void	setRunNumber (edm::RunNumber_t _r)

void	setSetupObjects (edm::EventSetup const &)

void	setSetupObjectsEndLumi (edm::EventSetup const &)

void	setTime (time_t _t)

void	setTokens (edm::ConsumesCollector &)

virtual	~DQWorker () noexcept(false)

Private Member Functions
void	setParams (edm::ParameterSet const &) override

Private Attributes
int	emptyLS_

int	emptyLSLimit_

bool	enable_ [nEEDCC]

int	isemptyLS

std::map< unsigned, float >	pnAmp_

unsigned	rtHalf_ [nEEDCC]

unsigned	wavelength_ [nEEDCC]

std::map< int, unsigned >	wlToME_

Additional Inherited Members
Public Types inherited from ecaldqm::DQWorkerTask
typedef EcalDCCHeaderBlock::EcalDCCEventSettings	EventSettings

Static Public Member Functions inherited from ecaldqm::DQWorkerTask
static void	fillDescriptions (edm::ParameterSetDescription &)

Static Public Member Functions inherited from ecaldqm::DQWorker
static void	fillDescriptions (edm::ParameterSetDescription &_desc)

Public Attributes inherited from ecaldqm::DQWorker
edm::ESGetToken < EcalElectronicsMapping, EcalMappingRcd >	elecMapHandle

edm::ESGetToken < EcalElectronicsMapping, EcalMappingRcd >	elecMapHandleEndLumi

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geomHandle

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geomHandleEndLumi

edm::ESGetToken< CaloTopology, CaloTopologyRecord >	topoHandle

edm::ESGetToken< CaloTopology, CaloTopologyRecord >	topoHandleEndLumi

edm::ESGetToken < EcalTrigTowerConstituentsMap, IdealGeometryRecord >	ttMapHandle

edm::ESGetToken < EcalTrigTowerConstituentsMap, IdealGeometryRecord >	ttMapHandleEndLumi

Protected Types inherited from ecaldqm::DQWorker
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Protected Member Functions inherited from ecaldqm::DQWorkerTask
void	setME (edm::ParameterSet const &) final

Protected Member Functions inherited from ecaldqm::DQWorker
void	initialize (std::string const &_name, edm::ParameterSet const &)

void	print_ (std::string const &, int=0) const

virtual void	setSource (edm::ParameterSet const &)

void	setVerbosity (int _verbosity)

Protected Attributes inherited from ecaldqm::DQWorker
bool	booked_

MESetCollection	MEs_

std::string	name_

bool	onlineMode_

Timestamp	timestamp_

int	verbosity_

bool	willConvertToEDM_

Detailed Description

Definition at line 14 of file LedTask.h.

Constructor & Destructor Documentation

ecaldqm::LedTask::LedTask ( )

Definition at line 9 of file LedTask.cc.

References enable_, ecaldqm::nEEDCC, rtHalf_, and wavelength_.

                    : DQWorkerTask(), wlToME_(), pnAmp_(), emptyLS_(0), emptyLSLimit_(0) {
     std::fill_n(enable_, nEEDCC, false);
     std::fill_n(wavelength_, nEEDCC, 0);
     std::fill_n(rtHalf_, nEEDCC, 0);
   }

ecaldqm::LedTask::~LedTask ( )

inlineoverride

Definition at line 17 of file LedTask.h.

17 {}

Member Function Documentation

void ecaldqm::LedTask::addDependencies ( DependencySet & _dependencies )

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 33 of file LedTask.cc.

References ecaldqm::kEcalRawData, ecaldqm::kEEDigi, ecaldqm::kEELaserLedUncalibRecHit, ecaldqm::kPnDiodeDigi, and ecaldqm::DependencySet::push_back().

                                                             {
     _dependencies.push_back(Dependency(kEEDigi, kEcalRawData));
     _dependencies.push_back(Dependency(kPnDiodeDigi, kEEDigi, kEcalRawData));
     _dependencies.push_back(Dependency(kEELaserLedUncalibRecHit, kPnDiodeDigi, kEEDigi, kEcalRawData));
   }

bool ecaldqm::LedTask::analyze	(	void const *	_p,
		Collections	_collection
	)

inlineoverridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 51 of file LedTask.h.

References ecaldqm::kEcalRawData, ecaldqm::kEEDigi, ecaldqm::kEELaserLedUncalibRecHit, ecaldqm::kPnDiodeDigi, runOnDigis(), runOnPnDigis(), runOnRawData(), and runOnUncalibRecHits().

                                                                       {
     switch (_collection) {
       case kEcalRawData:
         if (_p)
           runOnRawData(*static_cast<EcalRawDataCollection const*>(_p));
         return true;
         break;
       case kEEDigi:
         if (_p)
           runOnDigis(*static_cast<EEDigiCollection const*>(_p));
         return true;
         break;
       case kPnDiodeDigi:
         if (_p)
           runOnPnDigis(*static_cast<EcalPnDiodeDigiCollection const*>(_p));
         return true;
         break;
       case kEELaserLedUncalibRecHit:
         if (_p)
           runOnUncalibRecHits(*static_cast<EcalUncalibratedRecHitCollection const*>(_p));
         return true;
         break;
       default:
         break;
     }
 
     return false;
   }

void ecaldqm::LedTask::beginEvent	(	edm::Event const &	,
		edm::EventSetup const &	,
		bool const &	,
		bool &
	)

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 64 of file LedTask.cc.

References pnAmp_.

64 { pnAmp_.clear(); }

ecaldqm::LedTask::pnAmp_

std::map< unsigned, float > pnAmp_

Definition: LedTask.h:44

void ecaldqm::LedTask::beginLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

overridevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 58 of file LedTask.cc.

References emptyLS_, emptyLSLimit_, and isemptyLS.

                                                                                     {
     isemptyLS = 0;
     if (emptyLS_ + 1 > emptyLSLimit_)
       emptyLS_ = -1;
   }

void ecaldqm::LedTask::beginRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 56 of file LedTask.cc.

References emptyLS_.

56 { emptyLS_ = 0; }

ecaldqm::LedTask::emptyLS_

int emptyLS_

Definition: LedTask.h:46

void ecaldqm::LedTask::endLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

overridevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 336 of file LedTask.cc.

References emptyLS_, and isemptyLS.

                                                                                   {
     if (isemptyLS == 1)
       emptyLS_ += 1;
   }

bool ecaldqm::LedTask::filterRunType ( short const * _runType )

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 39 of file LedTask.cc.

References CreateSqliteForCondDB_cfg::enable, enable_, ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::kEEmHigh, EcalDCCHeaderBlock::LED_GAP, EcalDCCHeaderBlock::LED_STD, ecaldqm::nDCC, and ecaldqm::nEBDCC.

                                                    {
     bool enable(false);
 
     for (unsigned iDCC(0); iDCC != nDCC; iDCC++) {
       if (iDCC >= kEBmLow && iDCC <= kEBpHigh)
         continue;
       unsigned index(iDCC <= kEEmHigh ? iDCC : iDCC - nEBDCC);
       if (_runType[iDCC] == EcalDCCHeaderBlock::LED_STD || _runType[iDCC] == EcalDCCHeaderBlock::LED_GAP) {
         enable = true;
         enable_[index] = true;
       } else
         enable_[index] = false;
     }
 
     return enable;
   }

void ecaldqm::LedTask::runOnDigis ( EEDigiCollection const & _digis )

Definition at line 112 of file LedTask.cc.

References gpuClustering::adc, EcalMGPASample::adc(), ecaldqm::MESetCollection::at(), edm::DataFrameContainer::begin(), ecaldqm::dccId(), emptyLS_, CreateSqliteForCondDB_cfg::enable, enable_, edm::DataFrameContainer::end(), ecaldqm::MESetMulti::fill(), ecaldqm::DQWorker::getEcalDQMSetupObjects(), ecaldqm::DQWorker::GetElectronicsMap(), mps_fire::i, gpuClustering::id, isemptyLS, ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::kEEmHigh, SiStripPI::max, ecaldqm::DQWorker::MEs_, min(), ecaldqm::nEBDCC, ecaldqm::nEEDCC, DetId::null(), pnAmp_, ecaldqm::pnForCrystal(), DetId::rawId(), ecaldqm::rtHalf(), rtHalf_, EcalDataFrame::sample(), dtDQMClient_cfg::threshold, wavelength_, and wlToME_.

Referenced by analyze().

                                                          {
     MESet& meOccupancy(MEs_.at("Occupancy"));
     MESet& meShape(MEs_.at("Shape"));
     MESet& meSignalRate(MEs_.at("SignalRate"));
 
     int nReadouts[nEEDCC];
     int maxpos[nEEDCC][10];
     for (unsigned index(0); index < nEEDCC; ++index) {
       nReadouts[index] = 0;
       for (int i(0); i < 10; i++)
         maxpos[index][i] = 0;
     }
 
     for (EEDigiCollection::const_iterator digiItr(_digis.begin()); digiItr != _digis.end(); ++digiItr) {
       const DetId& id(digiItr->id());
 
       unsigned iDCC(dccId(id, GetElectronicsMap()) - 1);
       if (iDCC >= kEBmLow && iDCC <= kEBpHigh)
         continue;
       unsigned index(iDCC <= kEEmHigh ? iDCC : iDCC - nEBDCC);
 
       if (!enable_[index])
         continue;
       if (rtHalf(id, GetElectronicsMap()) != rtHalf_[index])
         continue;
 
       meOccupancy.fill(getEcalDQMSetupObjects(), id);
 
       ++nReadouts[index];
 
       EcalDataFrame dataFrame(*digiItr);
 
       int iMax(-1);
       int max(0);
       int 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 > 3)  // normal RMS of pedestal is ~2.5
         maxpos[index][iMax] += 1;
     }
 
     // signal existence check
     bool enable(false);
     bool ledOnExpected(emptyLS_ >= 0);
 
     unsigned iME(-1);
 
     for (int index(0); index < nEEDCC; ++index) {
       if (nReadouts[index] == 0) {
         enable_[index] = false;
         continue;
       }
 
       int threshold(nReadouts[index] / 3);
       if (ledOnExpected)
         enable_[index] = false;
 
       for (int i(0); i < 10; i++) {
         if (maxpos[index][i] > threshold) {
           enable = true;
           enable_[index] = true;
           break;
         }
       }
 
       if (iME != wlToME_[wavelength_[index]]) {
         iME = wlToME_[wavelength_[index]];
         static_cast<MESetMulti&>(meSignalRate).use(iME);
       }
 
       meSignalRate.fill(
           getEcalDQMSetupObjects(), (index <= kEEmHigh ? index : index + nEBDCC) + 1, enable_[index] ? 1 : 0);
     }
 
     if (!enable && isemptyLS >= 0)
       isemptyLS = 1;
     else if (enable)
       isemptyLS = -1;
 
     if (enable)
       emptyLS_ = 0;
     else if (ledOnExpected)
       return;
 
     iME = -1;
 
     for (EEDigiCollection::const_iterator digiItr(_digis.begin()); digiItr != _digis.end(); ++digiItr) {
       const DetId& id(digiItr->id());
 
       unsigned iDCC(dccId(id, GetElectronicsMap()) - 1);
       if (iDCC >= kEBmLow && iDCC <= kEBpHigh)
         continue;
       unsigned index(iDCC <= kEEmHigh ? iDCC : iDCC - nEBDCC);
 
       if (!enable_[index])
         continue;
       if (rtHalf(id, GetElectronicsMap()) != rtHalf_[index])
         continue;
 
       if (iME != wlToME_[wavelength_[index]]) {
         iME = wlToME_[wavelength_[index]];
         static_cast<MESetMulti&>(meShape).use(iME);
       }
 
       // EcalDataFrame is not a derived class of edm::DataFrame, but can take edm::DataFrame in the constructor
       EcalDataFrame dataFrame(*digiItr);
 
       for (int iSample(0); iSample < 10; iSample++)
         meShape.fill(getEcalDQMSetupObjects(), id, iSample + 0.5, float(dataFrame.sample(iSample).adc()));
 
       EcalPnDiodeDetId pnidA(pnForCrystal(id, 'a', GetElectronicsMap()));
       EcalPnDiodeDetId pnidB(pnForCrystal(id, 'b', GetElectronicsMap()));
       if (pnidA.null() || pnidB.null())
         continue;
       pnAmp_.insert(std::make_pair(pnidA.rawId(), 0.));
       pnAmp_.insert(std::make_pair(pnidB.rawId(), 0.));
     }
   }

void ecaldqm::LedTask::runOnPnDigis ( EcalPnDiodeDigiCollection const & _digis )

Definition at line 238 of file LedTask.cc.

References ecaldqm::MESetCollection::at(), edm::SortedCollection< T, SORT >::begin(), ecaldqm::dccId(), edm::SortedCollection< T, SORT >::end(), ecaldqm::MESetMulti::fill(), ecaldqm::DQWorker::getEcalDQMSetupObjects(), ecaldqm::DQWorker::GetElectronicsMap(), gpuClustering::id, ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::kEEmHigh, SiStripPI::max, ecaldqm::DQWorker::MEs_, ecaldqm::nEBDCC, EcalCondDBWriter_cfi::pedestal, pnAmp_, wavelength_, and wlToME_.

Referenced by analyze().

                                                                     {
     MESet& mePNAmplitude(MEs_.at("PNAmplitude"));
 
     unsigned iME(-1);
 
     for (EcalPnDiodeDigiCollection::const_iterator digiItr(_digis.begin()); digiItr != _digis.end(); ++digiItr) {
       if (digiItr->sample(0).gainId() != 0 && digiItr->sample(0).gainId() != 1)
         continue;
 
       const EcalPnDiodeDetId& id(digiItr->id());
 
       std::map<uint32_t, float>::iterator ampItr(pnAmp_.find(id.rawId()));
       if (ampItr == pnAmp_.end())
         continue;
 
       unsigned iDCC(dccId(id, GetElectronicsMap()) - 1);
       if (iDCC >= kEBmLow && iDCC <= kEBpHigh)
         continue;
       unsigned index(iDCC <= kEEmHigh ? iDCC : iDCC - nEBDCC);
 
       float pedestal(0.);
       for (int iSample(0); iSample < 4; iSample++)
         pedestal += digiItr->sample(iSample).adc();
       pedestal /= 4.;
 
       float max(0.);
       for (int iSample(0); iSample < 50; iSample++) {
         float amp(digiItr->sample(iSample).adc() - pedestal);
         if (amp > max)
           max = amp;
       }
 
       if (iME != wlToME_[wavelength_[index]]) {
         iME = wlToME_[wavelength_[index]];
         static_cast<MESetMulti&>(mePNAmplitude).use(iME);
       }
 
       mePNAmplitude.fill(getEcalDQMSetupObjects(), id, max);
 
       ampItr->second = max;
     }
   }

void ecaldqm::LedTask::runOnRawData ( EcalRawDataCollection const & _rawData )

Definition at line 66 of file LedTask.cc.

References ecaldqm::MESetCollection::at(), edm::SortedCollection< T, SORT >::begin(), enable_, edm::SortedCollection< T, SORT >::end(), ecaldqm::DQWorker::getEcalDQMSetupObjects(), ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::kEEmHigh, ecaldqm::DQWorker::MEs_, ecaldqm::nEBDCC, ecaldqm::nEEDCC, rtHalf_, wavelength_, and wlToME_.

Referenced by analyze().

                                                                   {
     MESet& meCalibStatus(MEs_.at("CalibStatus"));
     for (EcalRawDataCollection::const_iterator rItr(_rawData.begin()); rItr != _rawData.end(); ++rItr) {
       unsigned iDCC(rItr->id() - 1);
       if (iDCC >= kEBmLow && iDCC <= kEBpHigh)
         continue;
       unsigned index(iDCC <= kEEmHigh ? iDCC : iDCC - nEBDCC);
 
       if (!enable_[index]) {
         wavelength_[index] = -1;
         rtHalf_[index] = -1;
         continue;
       }
       if (rItr->getEventSettings().wavelength == 0)
         wavelength_[index] = 1;
       else if (rItr->getEventSettings().wavelength == 2)
         wavelength_[index] = 2;
       else
         wavelength_[index] = -1;
 
       if (wlToME_.find(wavelength_[index]) == wlToME_.end())
         enable_[index] = false;
 
       rtHalf_[index] = rItr->getRtHalf();
     }
     bool LedStatus[2];
     for (unsigned iW(0); iW < 2; iW++) {
       LedStatus[iW] = false;
     }
     for (unsigned index(0); index < nEEDCC; ++index) {
       switch (wavelength_[index]) {
         case 1:
           LedStatus[0] = true;
           break;
         case 2:
           LedStatus[1] = true;
           break;
         default:
           break;
       }
     }
     for (unsigned iWL(0); iWL < 2; iWL++) {
       meCalibStatus.fill(getEcalDQMSetupObjects(), double(iWL + 3), LedStatus[iWL] ? 1 : 0);
     }
   }

void ecaldqm::LedTask::runOnUncalibRecHits ( EcalUncalibratedRecHitCollection const & _uhits )

Definition at line 281 of file LedTask.cc.

References ecaldqm::MESetCollection::at(), edm::SortedCollection< T, SORT >::begin(), ecaldqm::dccId(), enable_, edm::SortedCollection< T, SORT >::end(), ecaldqm::MESetMulti::fill(), ecaldqm::DQWorker::getEcalDQMSetupObjects(), ecaldqm::DQWorker::GetElectronicsMap(), gpuClustering::id, if(), ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::kEEmHigh, SiStripPI::max, ecaldqm::DQWorker::MEs_, ecaldqm::nEBDCC, pnAmp_, ecaldqm::pnForCrystal(), ecaldqm::rtHalf(), rtHalf_, wavelength_, and wlToME_.

Referenced by analyze().

                                                                                   {
     using namespace std;
 
     MESet& meAmplitude(MEs_.at("Amplitude"));
     MESet& meAmplitudeSummary(MEs_.at("AmplitudeSummary"));
     MESet& meTiming(MEs_.at("Timing"));
     MESet& meAOverP(MEs_.at("AOverP"));
 
     unsigned iME(-1);
 
     for (EcalUncalibratedRecHitCollection::const_iterator uhitItr(_uhits.begin()); uhitItr != _uhits.end(); ++uhitItr) {
       EEDetId id(uhitItr->id());
 
       unsigned iDCC(dccId(id, GetElectronicsMap()) - 1);
       if (iDCC >= kEBmLow && iDCC <= kEBpHigh)
         continue;
       unsigned index(iDCC <= kEEmHigh ? iDCC : iDCC - nEBDCC);
 
       if (!enable_[index])
         continue;
       if (rtHalf(id, GetElectronicsMap()) != rtHalf_[index])
         continue;
 
       if (iME != wlToME_[wavelength_[index]]) {
         iME = wlToME_[wavelength_[index]];
         static_cast<MESetMulti&>(meAmplitude).use(iME);
         static_cast<MESetMulti&>(meAmplitudeSummary).use(iME);
         static_cast<MESetMulti&>(meTiming).use(iME);
         static_cast<MESetMulti&>(meAOverP).use(iME);
       }
 
       float amp(max((double)uhitItr->amplitude(), 0.));
       float jitter(max((double)uhitItr->jitter() + 5.0, 0.));
 
       meAmplitude.fill(getEcalDQMSetupObjects(), id, amp);
       meAmplitudeSummary.fill(getEcalDQMSetupObjects(), id, amp);
       meTiming.fill(getEcalDQMSetupObjects(), id, jitter);
 
       float aop(0.);
 
       map<uint32_t, float>::iterator ampItrA(pnAmp_.find(pnForCrystal(id, 'a', GetElectronicsMap())));
       map<uint32_t, float>::iterator ampItrB(pnAmp_.find(pnForCrystal(id, 'b', GetElectronicsMap())));
       if (ampItrA == pnAmp_.end() && ampItrB == pnAmp_.end())
         continue;
       else if (ampItrB == pnAmp_.end())
         aop = amp / ampItrA->second;
       else if (ampItrA == pnAmp_.end())
         aop = amp / ampItrB->second;
       else
         aop = amp / (ampItrA->second + ampItrB->second) * 2.;
 
       meAOverP.fill(getEcalDQMSetupObjects(), id, aop);
     }
   }

void ecaldqm::LedTask::setParams ( edm::ParameterSet const & _params )

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 15 of file LedTask.cc.

References ecaldqm::MESetCollection::at(), emptyLSLimit_, Exception, ecaldqm::MESetMulti::getIndex(), edm::ParameterSet::getUntrackedParameter(), ecaldqm::DQWorker::MEs_, and wlToME_.

                                                         {
     emptyLSLimit_ = _params.getUntrackedParameter<int>("emptyLSLimit");
 
     std::vector<int> ledWavelengths(_params.getUntrackedParameter<std::vector<int> >("ledWavelengths"));
 
     MESet::PathReplacements repl;
 
     MESetMulti& amplitude(static_cast<MESetMulti&>(MEs_.at("Amplitude")));
     unsigned nWL(ledWavelengths.size());
     for (unsigned iWL(0); iWL != nWL; ++iWL) {
       int wl(ledWavelengths[iWL]);
       if (wl != 1 && wl != 2)
         throw cms::Exception("InvalidConfiguration") << "Led Wavelength";
       repl["wl"] = std::to_string(wl);
       wlToME_[wl] = amplitude.getIndex(repl);
     }
   }