#include <TotemTimingConversions.h>

Public Member Functions
float	timeOfFirstSample (const TotemTimingDigi &digi) const

float	timePrecision (const TotemTimingDigi &digi) const

std::vector< float >	timeSamples (const TotemTimingDigi &digi) const

	TotemTimingConversions (bool mergeTimePeaks, const PPSTimingCalibration &calibration)

float	triggerTime (const TotemTimingDigi &digi) const

std::vector< float >	voltSamples (const TotemTimingDigi &digi) const

Private Attributes
PPSTimingCalibration	calibration_

reco::FormulaEvaluator	calibrationFunction_

bool	mergeTimePeaks_

Static Private Attributes
static constexpr int	ACCEPTED_TIME_RADIUS = 4

static constexpr unsigned long	CELL0_MASK = 0xfffffff000

static constexpr float	SAMPIC_ADC_V = 1. / 256

static constexpr int	SAMPIC_DEFAULT_OFFSET = 30

static constexpr int	SAMPIC_MAX_NUMBER_OF_SAMPLES = 64

static constexpr float	SAMPIC_SAMPLING_PERIOD_NS = 1. / 7.695

Detailed Description

Definition at line 22 of file TotemTimingConversions.h.

Constructor & Destructor Documentation

TotemTimingConversions::TotemTimingConversions	(	bool	mergeTimePeaks,
		const PPSTimingCalibration &	calibration
	)

Definition at line 16 of file TotemTimingConversions.cc.

                                                                                                            :
   calibration_(calibration), mergeTimePeaks_(mergeTimePeaks),
   calibrationFunction_(calibration_.formula())
 {}

Member Function Documentation

float TotemTimingConversions::timeOfFirstSample ( const TotemTimingDigi & digi ) const

Definition at line 24 of file TotemTimingConversions.cc.

Referenced by timeSamples(), and triggerTime().

 {
   unsigned int offsetOfSamples = digi.getEventInfo().getOffsetOfSamples();
   if (offsetOfSamples == 0)
     offsetOfSamples = SAMPIC_DEFAULT_OFFSET; // FW 0 is not sending this, FW > 0 yes
 
   unsigned int timestamp = (digi.getCellInfo() <= SAMPIC_MAX_NUMBER_OF_SAMPLES/2)
     ? digi.getTimestampA()
     : digi.getTimestampB();
 
   int cell0TimeClock = timestamp +
     ((digi.getFPGATimestamp()-timestamp) & CELL0_MASK)
     - digi.getEventInfo().getL1ATimestamp()
     + digi.getEventInfo().getL1ALatency();
 
   // time of first cell
   float cell0TimeInstant = SAMPIC_MAX_NUMBER_OF_SAMPLES * SAMPIC_SAMPLING_PERIOD_NS * cell0TimeClock;
 
   // time of triggered cell
   float firstCellTimeInstant = (digi.getCellInfo() < offsetOfSamples)
     ? cell0TimeInstant + digi.getCellInfo() * SAMPIC_SAMPLING_PERIOD_NS
     : cell0TimeInstant - (SAMPIC_MAX_NUMBER_OF_SAMPLES - digi.getCellInfo())*SAMPIC_SAMPLING_PERIOD_NS;
 
   int db = digi.getHardwareBoardId();
   int sampic = digi.getHardwareSampicId();
   int channel = digi.getHardwareChannelId();
   float t = firstCellTimeInstant + calibration_.timeOffset(db, sampic, channel);
   //NOTE: If no time offset is set, timeOffset returns 0
 
   if (mergeTimePeaks_) {
     if (t < -ACCEPTED_TIME_RADIUS)
       t += SAMPIC_MAX_NUMBER_OF_SAMPLES * SAMPIC_SAMPLING_PERIOD_NS;
     if (t > ACCEPTED_TIME_RADIUS)
       t -= SAMPIC_MAX_NUMBER_OF_SAMPLES * SAMPIC_SAMPLING_PERIOD_NS;
   }
   return t;
 }

float TotemTimingConversions::timePrecision ( const TotemTimingDigi & digi ) const

Definition at line 77 of file TotemTimingConversions.cc.

References calibration_, TotemTimingDigi::getHardwareBoardId(), TotemTimingDigi::getHardwareChannelId(), TotemTimingDigi::getHardwareSampicId(), and PPSTimingCalibration::timePrecision().

 {
   int db = digi.getHardwareBoardId();
   int sampic = digi.getHardwareSampicId();
   int channel = digi.getHardwareChannelId();
   return calibration_.timePrecision(db, sampic, channel);
 }

std::vector< float > TotemTimingConversions::timeSamples ( const TotemTimingDigi & digi ) const

Definition at line 88 of file TotemTimingConversions.cc.

References TotemTimingDigi::getNumberOfSamples(), mps_fire::i, SAMPIC_SAMPLING_PERIOD_NS, ntuplemaker::time, and timeOfFirstSample().

 {
   std::vector<float> time(digi.getNumberOfSamples());
   for (unsigned int i = 0; i < time.size(); ++i)
     time.at(i) = timeOfFirstSample(digi) + i * SAMPIC_SAMPLING_PERIOD_NS;
   return time;
 }

float TotemTimingConversions::triggerTime ( const TotemTimingDigi & digi ) const

Definition at line 65 of file TotemTimingConversions.cc.

References TotemTimingDigi::getEventInfo(), TotemTimingEventInfo::getOffsetOfSamples(), SAMPIC_MAX_NUMBER_OF_SAMPLES, SAMPIC_SAMPLING_PERIOD_NS, and timeOfFirstSample().

 {
   unsigned int offsetOfSamples = digi.getEventInfo().getOffsetOfSamples();
   if (offsetOfSamples == 0)
     offsetOfSamples = 30; // FW 0 is not sending this, FW > 0 yes
 
   return timeOfFirstSample(digi) + (SAMPIC_MAX_NUMBER_OF_SAMPLES - offsetOfSamples) * SAMPIC_SAMPLING_PERIOD_NS;
 }

std::vector< float > TotemTimingConversions::voltSamples ( const TotemTimingDigi & digi ) const

Definition at line 100 of file TotemTimingConversions.cc.

References calibration_, calibrationFunction_, data, reco::FormulaEvaluator::evaluate(), TotemTimingDigi::getCellInfo(), TotemTimingDigi::getHardwareBoardId(), TotemTimingDigi::getHardwareChannelId(), TotemTimingDigi::getHardwareSampicId(), TotemTimingDigi::getSamples(), reco::FormulaEvaluator::numberOfParameters(), reco::FormulaEvaluator::numberOfVariables(), PPSTimingCalibration::parameters(), SAMPIC_ADC_V, SAMPIC_MAX_NUMBER_OF_SAMPLES, and simplePhotonAnalyzer_cfi::sample.

 {
   std::vector<float> data;
   if (calibrationFunction_.numberOfVariables() != 1)
     for (const auto& sample : digi.getSamples())
       data.emplace_back(SAMPIC_ADC_V * sample);
   else {
     unsigned int db = digi.getHardwareBoardId();
     unsigned int sampic = digi.getHardwareSampicId();
     unsigned int channel = digi.getHardwareChannelId();
     unsigned int cell = digi.getCellInfo();
     for (const auto& sample : digi.getSamples()) {
       // ring buffer on Sampic, so accounting for samples register boundary
       const unsigned short sample_cell = (cell++) % SAMPIC_MAX_NUMBER_OF_SAMPLES;
       auto parameters = calibration_.parameters(db, sampic, channel, sample_cell);
       if (parameters.empty() || parameters.size() != calibrationFunction_.numberOfParameters())
         throw cms::Exception("TotemTimingConversions:voltSamples")
           << "Invalid calibrations retrieved for Sampic digi"
           << " (" << db << ", " << sampic << ", " << channel << ", " << sample_cell << ")!";
       data.emplace_back(calibrationFunction_.evaluate(std::vector<double>{(double)sample}, parameters));
     }
   }
   return data;
 }