#include <TotemT2RecHitProducerAlgorithm.h>

Inheritance diagram for TotemT2RecHitProducerAlgorithm:

Public Member Functions
void	build (const TotemGeometry &, const edmNew::DetSetVector< TotemT2Digi > &, edmNew::DetSetVector< TotemT2RecHit > &) override

	TimingRecHitProducerAlgorithm (const edm::ParameterSet &iConfig)

Public Member Functions inherited from TimingRecHitProducerAlgorithm< TotemGeometry, edmNew::DetSetVector< TotemT2Digi >, edmNew::DetSetVector< TotemT2RecHit > >
void	setCalibration (const PPSTimingCalibration &calib, const PPSTimingCalibrationLUT &calibLUT)

	TimingRecHitProducerAlgorithm (const edm::ParameterSet &iConfig)

virtual	~TimingRecHitProducerAlgorithm ()=default

Additional Inherited Members
Protected Attributes inherited from TimingRecHitProducerAlgorithm< TotemGeometry, edmNew::DetSetVector< TotemT2Digi >, edmNew::DetSetVector< TotemT2RecHit > >
bool	apply_calib_
	Switch on/off the timing calibration. More...

const PPSTimingCalibration *	calib_
	DB-loaded calibration object. More...

std::unique_ptr< reco::FormulaEvaluator >	calib_fct_
	Timing correction formula. More...

const PPSTimingCalibrationLUT *	calibLUT_

double	ts_to_ns_
	Conversion constant between time slice and absolute time (in ns) More...

Detailed Description

Definition at line 20 of file TotemT2RecHitProducerAlgorithm.h.

Member Function Documentation

◆ build()

void TotemT2RecHitProducerAlgorithm::build	(	const TotemGeometry &	geom,
		const edmNew::DetSetVector< TotemT2Digi > &	input,
		edmNew::DetSetVector< TotemT2RecHit > &	output
	)

overridevirtual

Implements TimingRecHitProducerAlgorithm< TotemGeometry, edmNew::DetSetVector< TotemT2Digi >, edmNew::DetSetVector< TotemT2RecHit > >.

Definition at line 15 of file TotemT2RecHitProducerAlgorithm.cc.

Referenced by TotemT2RecHitProducer::produce().

                                                                                       {
   for (const auto& vec : input) {
     const TotemT2DetId detid(vec.detId());
     const int sector = detid.arm(), plane = detid.plane(), channel = detid.channel();
 
     // retrieve the timing calibration part for this channel
     const auto& ch_params = (apply_calib_) ? calib_->parameters(sector, 0, plane, channel) : std::vector<double>{};
     // default values for offset + time precision if calibration object not found
     const double ch_t_offset = (apply_calib_) ? calib_->timeOffset(sector, 0, plane, channel) : 0.;
     const double ch_t_precis = (apply_calib_) ? calib_->timePrecision(sector, 0, plane, channel) : 0.;
 
     // prepare the output collection filler
     edmNew::DetSetVector<TotemT2RecHit>::FastFiller filler(output, detid);
 
     for (const auto& digi : vec) {
       const int t_lead = digi.leadingEdge(), t_trail = digi.trailingEdge();
       if (t_lead == 0 && t_trail == 0)  // skip invalid digis
         continue;
       double tot = -1., ch_t_twc = 0.;
       if (t_lead != 0 && t_trail != 0) {
         tot = (t_trail - t_lead) * ts_to_ns_;  // in ns
         if (calib_fct_ && apply_calib_) {      // compute the time-walk correction
           ch_t_twc = calib_fct_->evaluate(std::vector<double>{tot}, ch_params);
           if (edm::isNotFinite(ch_t_twc))
             ch_t_twc = 0.;
         }
       }
 
       // retrieve the geometry element associated to this DetID
       const auto& tile = geom.tile(detid);
 
       // store to the output collection
       filler.emplace_back(tile.centre(), t_lead * ts_to_ns_ - ch_t_offset - ch_t_twc, ch_t_precis, tot);
     }
   }
 }