1 | //---------------------------------------------------------------------- |
---|
2 | // This file distributed with FastJet has been obtained from SpartyJet |
---|
3 | // v2.20.0 by Pierre-Antoine Delsart, Kurtis L. Geerlings, Joey |
---|
4 | // Huston, Brian T. Martin and Chris Vermilion |
---|
5 | // For details, see http://www.pa.msu.edu/~huston/SpartyJet/ |
---|
6 | // http://projects.hepforge.org/spartyjet/ |
---|
7 | // |
---|
8 | // Changes from the original file are listed below. |
---|
9 | //---------------------------------------------------------------------- |
---|
10 | |
---|
11 | // History of changes from the original CommonUtils.hh file in |
---|
12 | // SpartyJet v2.20 |
---|
13 | // |
---|
14 | // 2009-01-15 Gregory Soyez <soyez@fastjet.fr> |
---|
15 | // |
---|
16 | // * put the code in the fastjet::atlas namespace |
---|
17 | |
---|
18 | #ifndef _JETCOMMONUTILS_HH_ |
---|
19 | #define _JETCOMMONUTILS_HH_ |
---|
20 | |
---|
21 | #include <ctime> |
---|
22 | #include <algorithm> |
---|
23 | #include <cmath> |
---|
24 | |
---|
25 | #include <fastjet/internal/base.hh> |
---|
26 | |
---|
27 | FASTJET_BEGIN_NAMESPACE |
---|
28 | |
---|
29 | namespace atlas{ |
---|
30 | |
---|
31 | // ************************************************************** |
---|
32 | // phi conversions |
---|
33 | // ************************************************************** |
---|
34 | inline float to_minusPI_PI(float phi){ |
---|
35 | while(phi < -M_PI) phi += 2*M_PI; |
---|
36 | while(phi >= M_PI) phi -= 2*M_PI; |
---|
37 | return phi; |
---|
38 | } |
---|
39 | inline float to_zero_2PI(float phi){ |
---|
40 | while(phi < 0) phi += 2*M_PI; |
---|
41 | while(phi >= 2*M_PI) phi -= 2*M_PI; |
---|
42 | return phi; |
---|
43 | } |
---|
44 | |
---|
45 | |
---|
46 | |
---|
47 | // ************************************************************** |
---|
48 | // List utils |
---|
49 | // ************************************************************** |
---|
50 | // Destroy all pointers in a container, and clear it |
---|
51 | // T must be a container of pointers ex. list<T2*> |
---|
52 | template<class T> |
---|
53 | void clear_list(T & list){ |
---|
54 | typedef typename T::iterator it_t; |
---|
55 | it_t it = list.begin(); |
---|
56 | it_t itE = list.end(); |
---|
57 | for(; it != itE; ++it){ |
---|
58 | delete *it; |
---|
59 | } |
---|
60 | list.clear(); |
---|
61 | } |
---|
62 | |
---|
63 | |
---|
64 | |
---|
65 | |
---|
66 | |
---|
67 | // ************************************************************** |
---|
68 | // timing |
---|
69 | // ************************************************************** |
---|
70 | class stopwatch { |
---|
71 | public : |
---|
72 | stopwatch() : m_total(0){}; |
---|
73 | void start(){m_last = std::clock();}; |
---|
74 | void resume(){m_last = std::clock();}; |
---|
75 | float pause() { |
---|
76 | std::clock_t now=std::clock(); |
---|
77 | m_total = m_total + now - m_last; |
---|
78 | m_last = now; |
---|
79 | return convert(); |
---|
80 | } |
---|
81 | float stop(){float t=pause(); m_total = std::clock_t(0); return t;} |
---|
82 | protected: |
---|
83 | std::clock_t m_last; |
---|
84 | std::clock_t m_total; |
---|
85 | |
---|
86 | float convert(){ return float(m_total)*1000/CLOCKS_PER_SEC;} |
---|
87 | }; |
---|
88 | |
---|
89 | |
---|
90 | } // namespace atlas |
---|
91 | |
---|
92 | FASTJET_END_NAMESPACE |
---|
93 | #endif |
---|