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

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

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

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

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

void	setParams (edm::ParameterSet const &) override

Private Attributes
float	chi2ThresholdEB_

float	chi2ThresholdEE_

float	energyThresholdEB_

float	energyThresholdEE_

MESet *	meTimeMapByLS

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

Constructor & Destructor Documentation

ecaldqm::TimingTask::TimingTask ( )

Definition at line 11 of file TimingTask.cc.

                          :
     DQWorkerTask(),
     chi2ThresholdEB_(0.),
     chi2ThresholdEE_(0.),
     energyThresholdEB_(0.),
     energyThresholdEE_(0.),
     meTimeMapByLS(0)
   {
   }

ecaldqm::TimingTask::~TimingTask ( )

inline

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 46 of file TimingTask.cc.

References ecaldqm::DQWorker::Timestamp::iLumi, ecaldqm::DQWorker::MEs_, meTimeMapByLS, ecaldqm::MESet::reset(), and ecaldqm::DQWorker::timestamp_.

   {
     // Fill separate MEs with only 10 LSs worth of stats
     // Used to correctly fill Presample Trend plots:
     // 1 pt:10 LS in Trend plots
     meTimeMapByLS = &MEs_.at("TimeMapByLS");
     if ( timestamp_.iLumi % 10 == 0 )
       meTimeMapByLS->reset();
   }

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

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 31 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 57 of file TimingTask.cc.

References edm::SortedCollection< T, SORT >::begin(), chi2ThresholdEB_, chi2ThresholdEE_, EcalBarrel, EcalEndcap, edm::SortedCollection< T, SORT >::end(), relval_parameters_module::energy, energyThresholdEB_, energyThresholdEE_, ecaldqm::MESet::fill(), hit::id, ecaldqm::kEBRecHit, EcalRecHit::kGood, EcalRecHit::kOutOfTime, ecaldqm::DQWorker::MEs_, meTimeMapByLS, and dtDQMClient_cfg::threshold.

Referenced by analyze().

   {
     MESet& meTimeAmp(MEs_.at("TimeAmp"));
     MESet& meTimeAmpAll(MEs_.at("TimeAmpAll"));
     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)));
     float threshold(_collection == kEBRecHit ? energyThresholdEB_ : energyThresholdEE_);
     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());
 
                     // Apply cut on chi2 of pulse shape fit
                     float chi2Threshold = ( id.subdetId() == EcalBarrel ) ? chi2ThresholdEB_ : chi2ThresholdEE_;
                     if ( id.subdetId() == EcalBarrel )
                       signedSubdet = EcalBarrel;
                     else {
                       EEDetId eeId( hit.id() );
                       if ( eeId.zside() < 0 )
                         signedSubdet = -EcalEndcap;
                       else
                         signedSubdet =  EcalEndcap;
                     }
                     if ( energy > threshold )
                       meChi2.fill(signedSubdet, hit.chi2());
                     if ( hit.chi2() > chi2Threshold ) return;
 
                     meTimeAmp.fill(id, energy, time);
                     meTimeAmpAll.fill(id, energy, time);
 
                     if(energy > threshold){
                       meTimeAll.fill(id, time);
                       meTimeMap.fill(id, time);
                       meTimeMapByLS->fill(id, time); 
                       meTime1D.fill(id, time);
                       meTimeAllMap.fill(id, time);
                     }
                   });
   }

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

Definition at line 110 of file TimingTask.cc.

References funct::abs(), edm::SortedCollection< T, SORT >::begin(), chi2ThresholdEB_, chi2ThresholdEE_, EcalBarrel, edm::SortedCollection< T, SORT >::end(), energyThresholdEB_, energyThresholdEE_, 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 = ( id.subdetId() == EcalBarrel ) ? chi2ThresholdEB_ : chi2ThresholdEE_;
       if( uhitItr->chi2() > chi2Threshold ) continue;
 
       // Apply amplitude cut based on approx rechit energy
       float amp( uhitItr->amplitude() );
       float ampThreshold( id.subdetId() == EcalBarrel ? energyThresholdEB_*20. : energyThresholdEE_*5. ); // 1 GeV ~ ( EB:20, EE:5 ) ADC
       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( abs(hitTime) >= 5. ) continue;
 
       // Fill MEs
       meTimeAmpBXm.fill( id,amp,uhitItr->outOfTimeAmplitude(4) ); // BX-1
       meTimeAmpBXp.fill( id,amp,uhitItr->outOfTimeAmplitude(6) ); // BX+1
 
     }
   }

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 22 of file TimingTask.cc.

References chi2ThresholdEB_, chi2ThresholdEE_, energyThresholdEB_, energyThresholdEE_, and edm::ParameterSet::getUntrackedParameter().

   {
     chi2ThresholdEB_   = _params.getUntrackedParameter<double>("chi2ThresholdEB");
     chi2ThresholdEE_   = _params.getUntrackedParameter<double>("chi2ThresholdEE");
     energyThresholdEB_ = _params.getUntrackedParameter<double>("energyThresholdEB");
     energyThresholdEE_ = _params.getUntrackedParameter<double>("energyThresholdEE");
   }