#include <TimingTask.h>

Inheritance diagram for ecaldqm::TimingTask:

Public Member Functions
bool	analyze (void const *, Collections) override

bool	filterRunType (short const *) override

void	runOnRecHits (EcalRecHitCollection const &, Collections)

void	runOnUncalibRecHits (EcalUncalibratedRecHitCollection const &)

	TimingTask ()

	~TimingTask () override

Public Member Functions inherited from ecaldqm::DQWorkerTask
virtual void	addDependencies (DependencySet &)

	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	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

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

virtual void	bookMEs (DQMStore::IBooker &)

	DQWorker ()

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

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	beginEvent (edm::Event const &, edm::EventSetup const &, bool const &, bool &) override

void	setParams (edm::ParameterSet const &) override

Private Attributes
double	bxBin_

std::vector< int >	bxBinEdges_

float	chi2ThresholdEB_

float	chi2ThresholdEE_

float	energyThresholdEB_

float	energyThresholdEE_

float	energyThresholdEEFwd_

MESet *	meTimeMapByLS

bool	splashSwitch_

float	timeErrorThreshold_

float	timingVsBXThreshold_

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 10 of file TimingTask.h.

Constructor & Destructor Documentation

ecaldqm::TimingTask::TimingTask ( )

Definition at line 13 of file TimingTask.cc.

       : DQWorkerTask(),
         bxBinEdges_(),
         bxBin_(0.),
         chi2ThresholdEB_(0.),
         chi2ThresholdEE_(0.),
         energyThresholdEB_(0.),
         energyThresholdEE_(0.),
         energyThresholdEEFwd_(0.),
         timingVsBXThreshold_(0.),
         timeErrorThreshold_(0.),
         meTimeMapByLS(nullptr) {}

ecaldqm::TimingTask::~TimingTask ( )

inlineoverride

Definition at line 13 of file TimingTask.h.

13 {}

Member Function Documentation

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

inlineoverridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 41 of file TimingTask.h.

References ecaldqm::kEBRecHit, ecaldqm::kEBUncalibRecHit, ecaldqm::kEERecHit, ecaldqm::kEEUncalibRecHit, runOnRecHits(), and runOnUncalibRecHits().

                                                                          {
     switch (_collection) {
       case kEBRecHit:
       case kEERecHit:
         if (_p)
           runOnRecHits(*static_cast<EcalRecHitCollection const*>(_p), _collection);
         return true;
         break;
       case kEBUncalibRecHit:
       case kEEUncalibRecHit:
         if (_p)
           runOnUncalibRecHits(*static_cast<EcalUncalibratedRecHitCollection const*>(_p));
         return true;
         break;
       default:
         break;
     }
     return false;
   }

void ecaldqm::TimingTask::beginEvent	(	edm::Event const &	_evt,
		edm::EventSetup const &	_es,
		bool const &	ByLumiResetSwitch,
		bool &
	)

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 51 of file TimingTask.cc.

References ecaldqm::MESetCollection::at(), edm::EventBase::bunchCrossing(), bxBin_, bxBinEdges_, ecaldqm::DQWorker::GetElectronicsMap(), ecaldqm::DQWorker::Timestamp::iLumi, ecaldqm::DQWorker::MEs_, meTimeMapByLS, ecaldqm::MESet::reset(), ecaldqm::DQWorker::timestamp_, and cuda_std::upper_bound().

                                                                                                                   {
     using namespace std;
     std::vector<int>::iterator pBin = std::upper_bound(bxBinEdges_.begin(), bxBinEdges_.end(), _evt.bunchCrossing());
     bxBin_ = static_cast<int>(pBin - bxBinEdges_.begin()) - 0.5;
     if (ByLumiResetSwitch) {
       meTimeMapByLS = &MEs_.at("TimeMapByLS");
       if (timestamp_.iLumi % 10 == 0)
         meTimeMapByLS->reset(GetElectronicsMap());
     }
   }

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

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 39 of file TimingTask.cc.

References EcalDCCHeaderBlock::COSMIC, EcalDCCHeaderBlock::COSMICS_GLOBAL, EcalDCCHeaderBlock::COSMICS_LOCAL, EcalDCCHeaderBlock::MTCC, ecaldqm::nDCC, EcalDCCHeaderBlock::PHYSICS_GLOBAL, and EcalDCCHeaderBlock::PHYSICS_LOCAL.

                                                       {
     for (int iFED(0); iFED < nDCC; iFED++) {
       if (_runType[iFED] == EcalDCCHeaderBlock::COSMIC || _runType[iFED] == EcalDCCHeaderBlock::MTCC ||
           _runType[iFED] == EcalDCCHeaderBlock::COSMICS_GLOBAL ||
           _runType[iFED] == EcalDCCHeaderBlock::PHYSICS_GLOBAL || _runType[iFED] == EcalDCCHeaderBlock::COSMICS_LOCAL ||
           _runType[iFED] == EcalDCCHeaderBlock::PHYSICS_LOCAL)
         return true;
     }
 
     return false;
   }

void ecaldqm::TimingTask::runOnRecHits	(	EcalRecHitCollection const &	_hits,
		Collections	_collection
	)

Definition at line 62 of file TimingTask.cc.

References ecaldqm::MESetCollection::at(), edm::SortedCollection< T, SORT >::begin(), bxBin_, chi2ThresholdEB_, chi2ThresholdEE_, EcalBarrel, EcalEndcap, edm::SortedCollection< T, SORT >::end(), relval_parameters_module::energy, ClusterTask_cfi::energyThreshold, energyThresholdEB_, energyThresholdEE_, energyThresholdEEFwd_, ecaldqm::MESet::fill(), ecaldqm::DQWorker::getEcalDQMSetupObjects(), gpuClustering::id, hit::id, ecaldqm::isForward(), EcalRecHit::kGood, EcalRecHit::kOutOfTime, ecaldqm::DQWorker::MEs_, meTimeMapByLS, ecaldqm::DQWorker::onlineMode_, splashSwitch_, timeErrorThreshold_, and timingVsBXThreshold_.

Referenced by analyze().

                                                                                           {
     MESet& meTimeAmp(MEs_.at("TimeAmp"));
     MESet& meTimeAmpAll(MEs_.at("TimeAmpAll"));
     MESet& meTimingVsBX(onlineMode_ ? MEs_.at("BarrelTimingVsBX") : MEs_.at("BarrelTimingVsBXFineBinned"));
     MESet& meTimeAll(MEs_.at("TimeAll"));
     MESet& meTimeAllMap(MEs_.at("TimeAllMap"));
     MESet& meTimeMap(MEs_.at("TimeMap"));  // contains cumulative run stats => not suitable for Trend plots
     MESet& meTime1D(MEs_.at("Time1D"));
     MESet& meChi2(MEs_.at("Chi2"));
 
     uint32_t mask(~((0x1 << EcalRecHit::kGood) | (0x1 << EcalRecHit::kOutOfTime)));
     int signedSubdet;
 
     std::for_each(_hits.begin(), _hits.end(), [&](EcalRecHitCollection::value_type const& hit) {
       if (hit.checkFlagMask(mask))
         return;
 
       DetId id(hit.id());
 
       float time(hit.time());
       float energy(hit.energy());
 
       float energyThreshold;
       if (id.subdetId() == EcalBarrel) {
         energyThreshold = energyThresholdEB_;
         signedSubdet = EcalBarrel;
       } else {
         energyThreshold = (isForward(id)) ? energyThresholdEEFwd_ : energyThresholdEE_;
         EEDetId eeId(hit.id());
         if (eeId.zside() < 0)
           signedSubdet = -EcalEndcap;
         else
           signedSubdet = EcalEndcap;
       }
 
       if (energy > energyThreshold)
         meChi2.fill(getEcalDQMSetupObjects(), signedSubdet, hit.chi2());
 
       if (!splashSwitch_) {  //Not applied for splash events
         float chi2Threshold;
         if (id.subdetId() == EcalBarrel)
           chi2Threshold = chi2ThresholdEB_;
         else
           chi2Threshold = chi2ThresholdEE_;
 
         //Apply cut on chi2 of pulse shape fit
         if (hit.chi2() > chi2Threshold)
           return;
       }
       // Apply cut based on timing error of rechit
       if (hit.timeError() > timeErrorThreshold_)
         return;
 
       meTimeAmp.fill(getEcalDQMSetupObjects(), id, energy, time);
       meTimeAmpAll.fill(getEcalDQMSetupObjects(), id, energy, time);
 
       if (energy > timingVsBXThreshold_ && signedSubdet == EcalBarrel)
         meTimingVsBX.fill(getEcalDQMSetupObjects(), bxBin_, time);
 
       if (energy > energyThreshold) {
         meTimeAll.fill(getEcalDQMSetupObjects(), id, time);
         meTimeMap.fill(getEcalDQMSetupObjects(), id, time);
         meTimeMapByLS->fill(getEcalDQMSetupObjects(), id, time);
         meTime1D.fill(getEcalDQMSetupObjects(), id, time);
         meTimeAllMap.fill(getEcalDQMSetupObjects(), id, time);
       }
     });
   }

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

Definition at line 134 of file TimingTask.cc.

References funct::abs(), ecaldqm::MESetCollection::at(), edm::SortedCollection< T, SORT >::begin(), chi2ThresholdEB_, chi2ThresholdEE_, EcalBarrel, edm::SortedCollection< T, SORT >::end(), energyThresholdEB_, energyThresholdEE_, energyThresholdEEFwd_, ecaldqm::DQWorker::getEcalDQMSetupObjects(), gpuClustering::id, ecaldqm::isForward(), EcalUncalibratedRecHit::kGood, and ecaldqm::DQWorker::MEs_.

Referenced by analyze().

                                                                                      {
     MESet& meTimeAmpBXm(MEs_.at("TimeAmpBXm"));
     MESet& meTimeAmpBXp(MEs_.at("TimeAmpBXp"));
 
     for (EcalUncalibratedRecHitCollection::const_iterator uhitItr(_uhits.begin()); uhitItr != _uhits.end(); ++uhitItr) {
       // Apply reconstruction quality cuts
       if (!uhitItr->checkFlag(EcalUncalibratedRecHit::kGood))
         continue;
       DetId id(uhitItr->id());
       float chi2Threshold = 0.;
       float ampThreshold = 0.;
       if (id.subdetId() == EcalBarrel) {
         chi2Threshold = chi2ThresholdEB_;
         ampThreshold = 20. * energyThresholdEB_;  // 1 GeV ~ 20 ADC in EB
       } else {
         chi2Threshold = chi2ThresholdEE_;
         ampThreshold = 5. * ((isForward(id)) ? energyThresholdEEFwd_ : energyThresholdEE_);  // 1 GeV ~ 5 ADC in EE
       }
 
       if (uhitItr->chi2() > chi2Threshold)
         continue;
 
       // Apply amplitude cut based on approx rechit energy
       float amp(uhitItr->amplitude());
       if (amp < ampThreshold)
         continue;
 
       // Apply jitter timing cut based on approx rechit timing
       float timeOff(id.subdetId() == EcalBarrel ? 0.4 : 1.8);
       float hitTime(uhitItr->jitter() * 25. + timeOff);  // 1 jitter ~ 25 ns
       if (std::abs(hitTime) >= 5.)
         continue;
 
       // Fill MEs
       meTimeAmpBXm.fill(getEcalDQMSetupObjects(), id, amp, uhitItr->outOfTimeAmplitude(4));  // BX-1
       meTimeAmpBXp.fill(getEcalDQMSetupObjects(), id, amp, uhitItr->outOfTimeAmplitude(6));  // BX+1
     }
   }

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 26 of file TimingTask.cc.

References bxBinEdges_, chi2ThresholdEB_, chi2ThresholdEE_, energyThresholdEB_, energyThresholdEE_, energyThresholdEEFwd_, edm::ParameterSet::getUntrackedParameter(), ecaldqm::DQWorker::onlineMode_, splashSwitch_, timeErrorThreshold_, and timingVsBXThreshold_.

                                                            {
     bxBinEdges_ = onlineMode_ ? _params.getUntrackedParameter<std::vector<int> >("bxBins")
                               : _params.getUntrackedParameter<std::vector<int> >("bxBinsFine");
     chi2ThresholdEB_ = _params.getUntrackedParameter<double>("chi2ThresholdEB");
     chi2ThresholdEE_ = _params.getUntrackedParameter<double>("chi2ThresholdEE");
     energyThresholdEB_ = _params.getUntrackedParameter<double>("energyThresholdEB");
     energyThresholdEE_ = _params.getUntrackedParameter<double>("energyThresholdEE");
     energyThresholdEEFwd_ = _params.getUntrackedParameter<double>("energyThresholdEEFwd");
     timingVsBXThreshold_ = _params.getUntrackedParameter<double>("timingVsBXThreshold");
     timeErrorThreshold_ = _params.getUntrackedParameter<double>("timeErrorThreshold");
     splashSwitch_ = _params.getUntrackedParameter<bool>("splashSwitch", false);
   }