source: trunk/FACT++/src/Time.cc@ 11209

Last change on this file since 11209 was 11208, checked in by tbretz, 14 years ago
Added a new constructor which takes special_values as argument.
File size: 8.5 KB
Line 
1// **************************************************************************
2/** @class Time
3
4@brief Adds some functionality to boost::posix_time::ptime for our needs
5
6This is basically a wrapper around boost::posix_time::ptime which is made
7to adapt the functionality to our needs. Time can store the current
8data and time with a precision up to nanoseconds if provided by the
9undrlaying system, otherwise microsecond precision is used.
10
11It main purpose is to provide needed constructors and simplyfy the
12conversion of dates and times from and to a string/stream.
13
14Note that posix_time (as Posix times have) has a limited range. You cannot
15use it for example for very early years of the last century.
16
17@section Examples
18
19 - An example can be found in \ref time.cc
20
21@section References
22
23 - <A HREF="http://www.boost.org/doc/libs/1_45_0/doc/html/date_time.html">BOOST++ date_time (V1.45.0)</A>
24
25**/
26// **************************************************************************
27#include "Time.h"
28
29using namespace std;
30using namespace boost::posix_time;
31
32const boost::gregorian::date Time::fUnixOffset(1970, 1, 1);
33
34const Time Time::None(Time::none);
35
36// strftime
37const _time_format Time::reset = 0;
38const _time_format Time::def = "%c";
39const _time_format Time::std = "%x %X%F";
40const _time_format Time::sql = "%Y-%m-%d %H:%M:%S.%f";
41const _time_format Time::ssql = "%Y-%m-%d %H:%M:%S";
42const _time_format Time::iso = "%Y%m%dT%H%M%S%F%q";
43const _time_format Time::magic = "%Y %m %d %H %M %S %f";
44const _time_format Time::smagic = "%Y %m %d %H %M %S";
45
46// --------------------------------------------------------------------------
47//
48//! Construct a Time object with either UTC or local time, or without any
49//! particular time.
50//!
51//! @param typ
52//! enum as defined in Time::init_t
53//
54Time::Time(enum init_t typ)
55{
56 switch (typ)
57 {
58 case utc:
59 *this = microsec_clock::universal_time();
60 break;
61 case local:
62 *this = microsec_clock::local_time();
63 break;
64 case none:
65 break;
66 }
67}
68
69
70// --------------------------------------------------------------------------
71//
72//! Construct a Time object with a date_time::special_value, e.g.
73//!
74//! - neg_infin
75//! - pos_infin
76//! - not_a_date_time
77//! - max_date_time
78//! - min_date_time
79//!
80//!
81//! @param val
82//! date_time::special_value
83//
84Time::Time(const boost::date_time::special_values &val) : ptime(val)
85{
86}
87
88// --------------------------------------------------------------------------
89//
90//! Construct a Time object from seconds since 1970/1/1 and number of
91//! milliseconds, as for example returned by gettimeofday()
92//!
93//! @param tm
94//! seconds since 1970/1/1
95//!
96//! @param millisec
97//! number of milliseconds
98//
99Time::Time(const time_t &tm, const int &millisec)
100: ptime(fUnixOffset, time_duration(0, 0, tm, millisec))
101{
102}
103
104// --------------------------------------------------------------------------
105//
106//! Construct a Time from a date and time.
107//!
108//! @param year, month, day, hh, mm, ss, microsec
109//! A full date and time down to microsecond precision. From the end
110//! arguments can be omitted.
111//!
112Time::Time(short year, unsigned char month, unsigned char day,
113 unsigned char hh, unsigned char mm, unsigned char ss, unsigned int microsec)
114// Last argument is fractional_seconds ( correct with num_fractional_digits() )
115: ptime(boost::gregorian::date(year, month, day),
116 time_duration(hh, mm, ss, microsec*pow(10, time_of_day().num_fractional_digits()-6)))
117{
118}
119
120// --------------------------------------------------------------------------
121//
122//! Set the Time object to a given MJD. Note that this involves
123//! conversion from double. So converting forth and back many many
124//! times might results in drifts.
125//!
126//! @param mjd
127//! Modified Julian Date
128//!
129void Time::Mjd(double mjd)
130{
131 // Convert MJD to seconds since offset
132 mjd -= 40587;
133 mjd *= 24*60*60;
134
135 const int exp = time_of_day().num_fractional_digits();
136 const double frac = fmod(mjd, 1)*pow(10, exp);
137
138 *this = ptime(fUnixOffset, time_duration(0, 0, mjd, frac));
139}
140
141// --------------------------------------------------------------------------
142//
143//! @returns the seconds of the day including the fractional seconds.
144//!
145double Time::SecondsOfDay() const
146{
147 const time_duration tod = time_of_day();
148
149 const int exp = tod.num_fractional_digits();
150
151 const double frac = tod.fractional_seconds()/pow(10, exp);
152 const double sec = tod.total_seconds()+frac;
153
154 return sec;
155}
156
157// --------------------------------------------------------------------------
158//
159//! Get the current MJD. Note that this involves
160//! conversion to double. So converting forth and back many many
161//! times might results in drifts.
162//!
163//! @returns
164//! Modified Julian Date
165//!
166double Time::Mjd() const
167{
168 return date().modjulian_day()+SecondsOfDay()/(24*60*60);
169
170 /*
171 const time_duration mjd = *this - ptime(fUnixOffset);
172 const double sec = mjd.total_seconds()+mjd.fractional_seconds()/1e6;
173 return sec/(24*60*60)+40587;
174 */
175}
176
177// --------------------------------------------------------------------------
178//
179// @returns seconds since 1970/1/1
180//
181double Time::UnixTime() const
182{
183 return (date().modjulian_day()-40587)*24*60*60 + SecondsOfDay();
184}
185
186// --------------------------------------------------------------------------
187//
188//! @returns the time in a format needed for root's TAxis
189//!
190double Time::RootTime() const
191{
192 return (date().modjulian_day()-49718)*24*60*60 + SecondsOfDay();
193}
194
195// --------------------------------------------------------------------------
196//
197//! Returns a string with the contents of the Time object formated
198//! as defined in format.
199//!
200//! @param format
201//! format description of the string to be returned. For details
202//! see the boost documentation or the man page of strftime
203//!
204//! @returns
205//! A string with the time formatted as requested. Note some special
206//! strings might be returned in case the time is invalid.
207//
208string Time::GetAsStr(const char *format) const
209{
210 stringstream out;
211 out << Time::fmt(format) << *this;
212 return out.str();
213}
214
215// --------------------------------------------------------------------------
216//
217//! Sets the time of the Time object to a time corresponding to
218//! the one given as argument. It is evaluated according to the given
219//! format.
220//!
221//! @param str
222//! The time as a string which should be converted to the Time object
223//!
224//! @param format
225//! format description of the string to be returned. For details
226//! see the boost documentation or the man page of strftime
227//!
228void Time::SetFromStr(const string &str, const char *format)
229{
230 // FIXME: exception handline
231 stringstream stream;
232 stream << str;
233 stream >> Time::fmt(format) >> *this;
234}
235
236// --------------------------------------------------------------------------
237//
238//! A stream manipulator which sets the streams Time output format
239//! as defined in the argument.
240//!
241//! @param format
242//! format description of the manipulator be returned. For details
243//! see the boost documentation or the man page of strftime
244//!
245//! @returns
246//! a stream manipulator for the given format
247//!
248const _time_format Time::fmt(const char *format)
249{
250 return format;
251}
252
253// --------------------------------------------------------------------------
254//
255//! Sets the locale discription of the stream (the way how a time is
256//! output) to the format defined by the given manipulator.
257//!
258//! Example:
259//! \code
260//! Time t();
261//! cout << Time::fmt("%Y:%m:%d %H:%M:%S.%f") << t << endl;
262//! \endcode
263//!
264//! @param out
265//! Reference to the stream
266//!
267//! @param f
268//! Time format described by a manipulator
269//!
270//! @returns
271//! A reference to the stream
272//!
273ostream &operator<<(ostream &out, const _time_format &f)
274{
275 const locale loc(locale::classic(),
276 f.ptr==0 ? 0 : new time_facet(f.ptr));
277
278 out.imbue(loc);
279
280 return out;
281}
282
283// --------------------------------------------------------------------------
284//
285//! Sets the locale discription of the stream (the way how a time is
286//! input) to the format defined by the given manipulator.
287//!
288//! Example:
289//! \code
290//! stringstream s;
291//! s << "09.09.1974 21:59";
292//!
293//! Time t;
294//! s >> Time::fmt("%d.%m.%Y %H:%M") >> t;
295//! \endcode
296//!
297//! @param in
298//! Reference to the stream
299//!
300//! @param f
301//! Time format described by a manipulator
302//!
303//! @returns
304//! A reference to the stream
305//!
306istream &operator>>(istream &in, const _time_format &f)
307{
308 const locale loc(locale::classic(),
309 f.ptr==0 ? 0 : new time_input_facet(f.ptr));
310
311 in.imbue(loc);
312
313 return in;
314}
315
Note: See TracBrowser for help on using the repository browser.