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