Loading...
Searching...
No Matches
CommonUtils.hh
//----------------------------------------------------------------------
// 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