Fix CA cell hardCurvCut, to apply absolute value

Author: mmasciov

RecoTracker/PixelSeeding/plugins/GPUCACell.h

@@ -208,7 +208,7 @@ class GPUCACellT {
 
     CircleEq<float> eq(x1, y1, x2, y2, x3, y3);
 
-    if (eq.curvature() > maxCurv)
+    if (std::abs(eq.curvature()) > maxCurv)
       return false;
 
     return std::abs(eq.dca0()) < region_origin_radius_plus_tolerance * std::abs(eq.curvature());
@@ -224,7 +224,7 @@ class GPUCACellT {
                                                  const float maxCurv) {
     CircleEq<float> eq(x1, y1, x2, y2, x3, y3);
 
-    if (eq.curvature() > maxCurv)
+    if (std::abs(eq.curvature()) > maxCurv)
       return false;
 
     return std::abs(eq.dca0()) < region_origin_radius_plus_tolerance * std::abs(eq.curvature());

RecoTracker/PixelSeeding/plugins/alpaka/CACell.h

@@ -213,7 +213,7 @@ namespace ALPAKA_ACCELERATOR_NAMESPACE {
 
       CircleEq<float> eq(x1, y1, x2, y2, x3, y3);
 
-      if (eq.curvature() > maxCurv)
+      if (std::abs(eq.curvature()) > maxCurv)
         return false;
 
       return std::abs(eq.dca0()) < region_origin_radius_plus_tolerance * std::abs(eq.curvature());
@@ -230,7 +230,7 @@ namespace ALPAKA_ACCELERATOR_NAMESPACE {
         const float maxCurv) {
       CircleEq<float> eq(x1, y1, x2, y2, x3, y3);
 
-      if (eq.curvature() > maxCurv)
+      if (std::abs(eq.curvature()) > maxCurv)
         return false;
 
       return std::abs(eq.dca0()) < region_origin_radius_plus_tolerance * std::abs(eq.curvature());