#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 &) 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 ()

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

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

virtual bool	filterTrigger (edm::TriggerResultsByName const &)

void	recoverStats ()

virtual void	setTokens (edm::ConsumesCollector &)

void	softReset ()

virtual	~DQWorkerTask ()

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

	DQWorker ()

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

std::string const &	getName () const

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	setTime (time_t _t)

virtual	~DQWorker ()

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

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)

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::DQWorkerTask
std::set< std::string >	resettable_

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 ( )

inline

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 41 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 50 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 &
	)

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 82 of file LedTask.cc.

References pnAmp_.

   {
     pnAmp_.clear();
   }

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

overridevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 75 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 69 of file LedTask.cc.

References emptyLS_.

   {
     emptyLS_ = 0;
   }

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

overridevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 340 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 49 of file LedTask.cc.

References enable_, cmsHarvester::index, 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 136 of file LedTask.cc.

References ecalMGPA::adc(), EcalMGPASample::adc(), edm::DataFrameContainer::begin(), ecaldqm::dccId(), emptyLS_, enable_, edm::DataFrameContainer::end(), ecaldqm::MESetMulti::fill(), i, cmsHarvester::index, isemptyLS, ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::kEEmHigh, bookConverter::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) - 1);
       if(iDCC >= kEBmLow && iDCC <= kEBpHigh) continue;
       unsigned index(iDCC <= kEEmHigh ? iDCC : iDCC - nEBDCC);
 
       if(!enable_[index]) continue;
       if(rtHalf(id) != rtHalf_[index]) continue;
 
       meOccupancy.fill(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((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) - 1);
       if(iDCC >= kEBmLow && iDCC <= kEBpHigh) continue;
       unsigned index(iDCC <= kEEmHigh ? iDCC : iDCC - nEBDCC);
 
       if(!enable_[index]) continue;
       if(rtHalf(id) != 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(id, iSample + 0.5, float(dataFrame.sample(iSample).adc()));
 
       EcalPnDiodeDetId pnidA(pnForCrystal(id, 'a'));
       EcalPnDiodeDetId pnidB(pnForCrystal(id, 'b'));
       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 249 of file LedTask.cc.

References edm::SortedCollection< T, SORT >::begin(), ecaldqm::dccId(), edm::SortedCollection< T, SORT >::end(), ecaldqm::MESetMulti::fill(), cmsHarvester::index, ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::kEEmHigh, bookConverter::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) - 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(id, max);
 
       ampItr->second = max;
     }
   }

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

Definition at line 88 of file LedTask.cc.

References edm::SortedCollection< T, SORT >::begin(), enable_, edm::SortedCollection< T, SORT >::end(), cmsHarvester::index, 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(double(iWL+3), LedStatus[iWL]? 1:0);
     }  
   }

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

Definition at line 290 of file LedTask.cc.

References edm::SortedCollection< T, SORT >::begin(), ecaldqm::dccId(), enable_, edm::SortedCollection< T, SORT >::end(), ecaldqm::MESetMulti::fill(), cmsHarvester::index, ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::kEEmHigh, bookConverter::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) - 1);
       if(iDCC >= kEBmLow && iDCC <= kEBpHigh) continue;
       unsigned index(iDCC <= kEEmHigh ? iDCC : iDCC - nEBDCC);
 
       if(!enable_[index]) continue;
       if(rtHalf(id) != 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(id, amp);
       meAmplitudeSummary.fill(id, amp);
       meTiming.fill(id, jitter);
 
       float aop(0.);
 
       map<uint32_t, float>::iterator ampItrA(pnAmp_.find(pnForCrystal(id, 'a')));
       map<uint32_t, float>::iterator ampItrB(pnAmp_.find(pnForCrystal(id, 'b')));
       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(id, aop);
     }
   }

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

privatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 22 of file LedTask.cc.

References 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);
     }
   }