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
|
---|