TrkVKalUtils.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TRKVKALVRTCORE_TRKVKALUTILS_H
#define TRKVKALVRTCORE_TRKVKALUTILS_H
#include "boost/core/noinit_adaptor.hpp"
#include <math.h>
#include <algorithm>
#include <vector>
 
namespace Trk {
 
  struct Vect3DF{
    double X; double Y; double Z;
    Vect3DF() noexcept  { X=Y=Z=0.; }
    void Set(double arr[3]) noexcept { X=arr[0]; Y=arr[1]; Z=arr[2]; }
    double Dist3D( Vect3DF & v) const noexcept  { return sqrt(pow(X-v.X,2)+pow(Y-v.Y,2)+pow(Z-v.Z,2)); }
  };
 
  struct VectMOM { double Px; double Py; double Pz; double E; };
 
  inline void cfdcopy(double *source, double *target, int n)
    {  std::copy(source, source + n, target); }
 
  inline double cfddist3D(double *V1, double *V2)
    { return sqrt( (V2[0]-V1[0])*(V2[0]-V1[0]) + (V2[1]-V1[1])*(V2[1]-V1[1]) + (V2[2]-V1[2])*(V2[2]-V1[2]) ); }
 
 
//
//  Implementation of fortran function SIGN()
//------------------------------------------------
  inline double d_sign(double value, double sign)
    {
      if( value >= 0){
        if (sign >=0) return  value;
        if (sign < 0) return -value;
      }
      if( value < 0){
        if (sign >=0) return -value;
        if (sign < 0) return  value;
      }
      return value;
   }
 
 
   template <class T>
   using noinit_vector = std::vector<T, boost::noinit_adaptor<std::allocator<T> > >;
 
 
} /* End namespace */
#endif