DUNE-DAQ
DUNE Trigger and Data Acquisition software
Loading...
Searching...
No Matches
dunedaq::utilities::TimestampEstimator Class Reference

TimestampEstimator is an implementation of TimestampEstimatorBase that uses TimeSync messages from an input queue to estimate the current timestamp. More...

#include <TimestampEstimator.hpp>

Inheritance diagram for dunedaq::utilities::TimestampEstimator:
[legend]
Collaboration diagram for dunedaq::utilities::TimestampEstimator:
[legend]

Classes

struct  TimeSyncPoint
 

Public Member Functions

 TimestampEstimator (uint32_t run_number, uint64_t clock_frequency_hz)
 
 TimestampEstimator (uint64_t clock_frequency_hz)
 
virtual ~TimestampEstimator ()
 
uint64_t get_timestamp_estimate () const override
 Returns the current timestamp estimate or a special value if no valid timestamp is available.
 
std::chrono::microseconds get_wait_estimate (uint64_t ts) const override
 
void add_timestamp_datapoint (uint64_t daq_time, uint64_t system_time)
 
template<class T >
void timesync_callback (const T &tsync)
 
uint64_t get_received_timesync_count () const
 
- Public Member Functions inherited from dunedaq::utilities::TimestampEstimatorBase
virtual ~TimestampEstimatorBase ()=default
 
WaitStatus wait_for_valid_timestamp (std::atomic< bool > &continue_flag)
 Waits for a valid timestamp to become available. Returns a status value that indicates whether a valid timestamp is available or not.
 
WaitStatus wait_for_timestamp (uint64_t ts, std::atomic< bool > &continue_flag)
 

Private Attributes

std::atomic< TimeSyncPointm_current_timestamp_estimate
 
uint64_t m_clock_frequency_hz
 
uint64_t m_most_recent_daq_time
 
uint64_t m_most_recent_system_time
 
std::mutex m_datapoint_mutex
 
uint32_t m_run_number {0}
 
std::atomic< uint64_t > m_received_timesync_count
 
uint32_t m_current_process_id
 

Additional Inherited Members

- Public Types inherited from dunedaq::utilities::TimestampEstimatorBase
enum  WaitStatus { kFinished , kInterrupted }
 

Detailed Description

TimestampEstimator is an implementation of TimestampEstimatorBase that uses TimeSync messages from an input queue to estimate the current timestamp.

Definition at line 27 of file TimestampEstimator.hpp.

Constructor & Destructor Documentation

◆ TimestampEstimator() [1/2]

dunedaq::utilities::TimestampEstimator::TimestampEstimator ( uint32_t run_number,
uint64_t clock_frequency_hz )

Definition at line 21 of file TimestampEstimator.cpp.

22 : TimestampEstimator(clock_frequency_hz)
23{
24 m_run_number = run_number;
25}
TimestampEstimator(uint32_t run_number, uint64_t clock_frequency_hz)

◆ TimestampEstimator() [2/2]

dunedaq::utilities::TimestampEstimator::TimestampEstimator ( uint64_t clock_frequency_hz)
explicit

Definition at line 27 of file TimestampEstimator.cpp.

29 TimeSyncPoint{ std::numeric_limits<uint64_t>::max(), std::chrono::time_point<std::chrono::steady_clock>() })
30 , m_clock_frequency_hz(clock_frequency_hz)
33 , m_run_number(0)
35{
36 m_current_process_id = static_cast<uint32_t>(getpid());
37}
std::atomic< TimeSyncPoint > m_current_timestamp_estimate
std::atomic< uint64_t > m_received_timesync_count

◆ ~TimestampEstimator()

dunedaq::utilities::TimestampEstimator::~TimestampEstimator ( )
virtual

Definition at line 39 of file TimestampEstimator.cpp.

39{}

Member Function Documentation

◆ add_timestamp_datapoint()

void dunedaq::utilities::TimestampEstimator::add_timestamp_datapoint ( uint64_t daq_time,
uint64_t system_time )

Definition at line 74 of file TimestampEstimator.cpp.

75{
76 using namespace std::chrono;
77
78 std::scoped_lock<std::mutex> lk(m_datapoint_mutex);
79
80 // First, update the latest timestamp
81 TimeSyncPoint estimate = m_current_timestamp_estimate.load();
82 int64_t diff = estimate.daq_time - daq_time;
83 TLOG_DEBUG(TLVL_TIME_SYNC_PROPERTIES) << "Got a TimeSync timestamp = " << daq_time
84 << ", system time = " << system_time << " when current timestamp estimate was "
85 << estimate.daq_time << ". diff=" << diff;
86
87 if (m_most_recent_daq_time == std::numeric_limits<uint64_t>::max() || daq_time > m_most_recent_daq_time) {
88 m_most_recent_daq_time = daq_time;
89 m_most_recent_system_time = system_time;
90 }
91
92 if (m_most_recent_daq_time != std::numeric_limits<uint64_t>::max()) {
93 // Update the current timestamp estimate, based on the most recently-read TimeSync
94 using namespace std::chrono;
95
96 auto time_now =
97 static_cast<uint64_t>(duration_cast<microseconds>(system_clock::now().time_since_epoch()).count()); // NOLINT
98 auto steady_time_now = steady_clock::now();
99
100 // (PAR 2021-07-22) We only want to _increase_ our timestamp
101 // estimate, not _decrease_ it, so we only attempt the update if
102 // our system time is later than the latest time sync's system
103 // time. We can get TimeSync messages from the "future" if
104 // they're coming from another host whose clock is not exactly
105 // synchronized with ours: that's fine, but if the discrepancy
106 // is large, then badness could happen, so emit a warning
107
108 if (time_now < m_most_recent_system_time - 10000) {
109 ers::warning(EarlyTimeSync(ERS_HERE, m_most_recent_system_time - time_now));
110 }
111
112 if (time_now > m_most_recent_system_time) {
113
114 auto delta_time = time_now - m_most_recent_system_time;
116 << "Time diff between current system and latest TimeSync system time [us]: " << delta_time;
117
118 // Warn user if current system time is more than 1s ahead of latest TimeSync system time. This could be a sign of
119 // an issue, e.g. machine times out of sync
120 if (delta_time > 1e6)
121 ers::warning(LateTimeSync(ERS_HERE, delta_time));
122
123 const uint64_t new_timestamp = m_most_recent_daq_time + delta_time * m_clock_frequency_hz / 1000000;
124
125 // Don't ever decrease the timestamp; just wait until enough
126 // time passes that we want to increase it
127 if (estimate.daq_time == std::numeric_limits<uint64_t>::max() || new_timestamp >= estimate.daq_time) {
129 << "Storing new timestamp estimate of " << new_timestamp << " ticks (..." << std::fixed
130 << std::setprecision(8)
131 << (static_cast<double>(new_timestamp % (m_clock_frequency_hz * 1000)) /
132 static_cast<double>(m_clock_frequency_hz))
133 << " sec), mrt.daq_time is " << m_most_recent_daq_time << " ticks (..."
134 << (static_cast<double>(m_most_recent_daq_time % (m_clock_frequency_hz * 1000)) /
135 static_cast<double>(m_clock_frequency_hz))
136 << " sec), delta_time is " << delta_time << " usec, clock_freq is " << m_clock_frequency_hz << " Hz";
137 m_current_timestamp_estimate.store(TimeSyncPoint{ new_timestamp, steady_time_now });
138 } else {
139 TLOG_DEBUG(TLVL_TIME_SYNC_NOTES) << "Not updating timestamp estimate backwards from "
140 << m_current_timestamp_estimate.load().daq_time << " to " << new_timestamp;
141 }
142 }
143 }
144}
#define ERS_HERE
#define TLOG_DEBUG(lvl,...)
Definition Logging.hpp:112
void warning(const Issue &issue)
Definition ers.hpp:115
#define TLVL_TIME_SYNC_PROPERTIES
Definition Issues.hpp:16
#define TLVL_TIME_SYNC_NEW_ESTIMATE
Definition Issues.hpp:18
#define TLVL_TIME_SYNC_NOTES
Definition Issues.hpp:17

◆ get_received_timesync_count()

uint64_t dunedaq::utilities::TimestampEstimator::get_received_timesync_count ( ) const
inline

Definition at line 45 of file TimestampEstimator.hpp.

45{ return m_received_timesync_count.load(); }

◆ get_timestamp_estimate()

uint64_t dunedaq::utilities::TimestampEstimator::get_timestamp_estimate ( ) const
overridevirtual

Returns the current timestamp estimate or a special value if no valid timestamp is available.

Returns
the current estimated timestamp (in units of DUNE Timing System ticks) or std::numeric_limits<uint64_t>::max() if no valid timestamp is currently available

Implements dunedaq::utilities::TimestampEstimatorBase.

Definition at line 47 of file TimestampEstimator.cpp.

48{
49 using namespace std::chrono;
50
51 TimeSyncPoint estimate = m_current_timestamp_estimate.load();
52 // 27-May-2025, KAB: added check if a valid timestamp is available and, if not, return early
53 // with the special value that indicates that none is available.
54 if (estimate.daq_time == std::numeric_limits<uint64_t>::max()) {return estimate.daq_time;}
55
56 auto delta_time_us = duration_cast<microseconds>(steady_clock::now() - estimate.system_time).count();
57
58 const uint64_t new_timestamp = estimate.daq_time + delta_time_us * m_clock_frequency_hz / 1000000;
59
60 return new_timestamp;
61}

◆ get_wait_estimate()

std::chrono::microseconds dunedaq::utilities::TimestampEstimator::get_wait_estimate ( uint64_t ts) const
overridevirtual

Implements dunedaq::utilities::TimestampEstimatorBase.

Definition at line 64 of file TimestampEstimator.cpp.

65{
67 if (now > ts)
68 return std::chrono::microseconds(0);
69 auto diff = ts - now;
70 return std::chrono::microseconds(static_cast<long>(diff * 1000000. / m_clock_frequency_hz));
71}
uint64_t get_timestamp_estimate() const override
Returns the current timestamp estimate or a special value if no valid timestamp is available.
static int64_t now()
PDS Frame with unphysical timestamp detected with ts

◆ timesync_callback()

template<class T >
void dunedaq::utilities::TimestampEstimator::timesync_callback ( const T & tsync)

Definition at line 7 of file TimestampEstimator.hxx.

8{
10 TLOG_DEBUG(TLVL_TIME_SYNC_PROPERTIES) << "Got a TimeSync run=" << tsync.run_number << " local run=" << m_run_number
11 << " seqno=" << tsync.sequence_number
12 << " source_pid=" << tsync.source_pid;
13 if (tsync.run_number == m_run_number && tsync.source_pid != m_current_process_id) {
14 add_timestamp_datapoint(tsync.daq_time, tsync.system_time);
15 } else {
16 TLOG_DEBUG(0) << "Discarded TimeSync message from run " << tsync.run_number << " during run "
17 << m_run_number << " with pid " << tsync.source_pid << " and timestamp " << tsync.daq_time;
18 }
19}
void add_timestamp_datapoint(uint64_t daq_time, uint64_t system_time)

Member Data Documentation

◆ m_clock_frequency_hz

uint64_t dunedaq::utilities::TimestampEstimator::m_clock_frequency_hz
private

Definition at line 56 of file TimestampEstimator.hpp.

◆ m_current_process_id

uint32_t dunedaq::utilities::TimestampEstimator::m_current_process_id
private

Definition at line 62 of file TimestampEstimator.hpp.

◆ m_current_timestamp_estimate

std::atomic<TimeSyncPoint> dunedaq::utilities::TimestampEstimator::m_current_timestamp_estimate
private

Definition at line 53 of file TimestampEstimator.hpp.

◆ m_datapoint_mutex

std::mutex dunedaq::utilities::TimestampEstimator::m_datapoint_mutex
private

Definition at line 59 of file TimestampEstimator.hpp.

◆ m_most_recent_daq_time

uint64_t dunedaq::utilities::TimestampEstimator::m_most_recent_daq_time
private

Definition at line 57 of file TimestampEstimator.hpp.

◆ m_most_recent_system_time

uint64_t dunedaq::utilities::TimestampEstimator::m_most_recent_system_time
private

Definition at line 58 of file TimestampEstimator.hpp.

◆ m_received_timesync_count

std::atomic<uint64_t> dunedaq::utilities::TimestampEstimator::m_received_timesync_count
private

Definition at line 61 of file TimestampEstimator.hpp.

◆ m_run_number

uint32_t dunedaq::utilities::TimestampEstimator::m_run_number {0}
private

Definition at line 60 of file TimestampEstimator.hpp.

60{0};

The documentation for this class was generated from the following files: