Distance3DNormed.h

/**************************************************************************
* BASF2 (Belle Analysis Framework 2)                                     *
* Copyright(C) 2014 - Belle II Collaboration                             *
*                                                                        *
* Author: The Belle II Collaboration                                     *
* Contributors: Jakob Lettenbichler (jakob.lettenbichler@oeaw.ac.at)     *
*                                                                        *
* This software is provided "as is" without any warranty.                *
**************************************************************************/
 
#pragma once
 
#include <tracking/trackFindingVXD/filterMap/filterFramework/SelectionVariable.h>
#include <cmath>
 
#define DISTANCE3DNORMED_NAME Distance3DNormed
namespace Belle2 {
  template <typename PointType >
  class DISTANCE3DNORMED_NAME : public SelectionVariable< PointType , 2, double > {
  public:
 
    PUT_NAME_FUNCTION(DISTANCE3DNORMED_NAME);
 
 
 
    static double value(const PointType& outerHit, const PointType& innerHit)
    {
 
      double result =
        (std::pow(outerHit.X() - innerHit.X() , 2)
         + std::pow(outerHit.Y() - innerHit.Y() , 2))
        /
        (std::pow(outerHit.X() - innerHit.X() , 2)
         + std::pow(outerHit.Y() - innerHit.Y() , 2)
         + std::pow(outerHit.Z() - innerHit.Z() , 2));
 
      return
        (std::isnan(result) || std::isinf(result)) ? 0 : result;
    }
 
 
  };
 
}