VectorUtil.cxx

Go to the documentation of this file.

// @(#)root/mathcore:$Id$
// Authors: W. Brown, M. Fischler, L. Moneta    2006
 
 /**********************************************************************
  *                                                                    *
  * Copyright (c) 2006 , LCG ROOT MathLib Team                         *
  *                                                                    *
  *                                                                    *
  **********************************************************************/
 
// Implementation of VectorUtil functions
//
// Created by: Lorenzo Moneta 22 Aug 2006
//
#include "Math/GenVector/VectorUtil.h"
 
namespace ROOT {
 
namespace Math {
 
 
double VectorUtil::Phi_0_2pi(double angle) {
   // returns phi angle in the interval (0,2*PI]
   if ( angle <= 2.*M_PI && angle > 0 ) return angle;
 
   if ( angle > 0 ) {
      int n = static_cast<int>( angle/(2.*M_PI) );
      angle -= 2.*M_PI*n;
   } else {
      int n = static_cast<int>( -(angle)/(2.*M_PI) );
      angle += 2.*M_PI*(n+1);
   }
   return angle;
}
 
double VectorUtil::Phi_mpi_pi(double angle) {
   // returns phi angle in the interval (-PI,PI]
 
   if ( angle <= M_PI && angle > -M_PI ) return angle;
 
   if ( angle > 0 ) {
      int n = static_cast<int>( (angle+M_PI)/(2.*M_PI) );
      angle -= 2*M_PI*n;
   } else {
      int n = static_cast<int>( -(angle-M_PI)/(2.*M_PI) );
      angle += 2*M_PI*n;
   }
   return angle;
}
 
 
 
} //namespace Math
} //namespace ROOT