diff --git a/source/himalaya/mh2_eft/ThresholdLoopFunctions.cpp b/source/himalaya/mh2_eft/ThresholdLoopFunctions.cpp
index f0ee137..e53cc78 100644
--- a/source/himalaya/mh2_eft/ThresholdLoopFunctions.cpp
+++ b/source/himalaya/mh2_eft/ThresholdLoopFunctions.cpp
@@ -1231,74 +1231,6 @@ double f8(double r1, double r2) noexcept
    return 1.5 * result;
 }
 
-double fth1(double y) noexcept
-{
-   const double eps = 10.*std::numeric_limits<double>::epsilon();
-
-   if (is_zero(y, eps)) {
-      return 0.;
-   }
-
-   if (is_equal(std::abs(y), 1.0, eps)) {
-      return -1.;
-   }
-
-   if (!is_zero(y, eps) && !is_equal(std::abs(y), 1.0, eps)) {
-      const double y2 = sqr(y);
-
-      return y2*std::log(y2) / (1. - y2);
-   }
-
-   return 0.;
-}
-
-double fth2(double y) noexcept
-{
-   const double eps = 10.*std::numeric_limits<double>::epsilon();
-
-   if (is_zero(y, eps)) {
-      return 0.;
-   }
-
-   if (is_equal(std::abs(y), 1.0, eps)) {
-      return 0.5;
-   }
-
-   if (!is_zero(y, eps) && !is_equal(std::abs(y), 1.0, eps)) {
-      const double y2 = sqr(y);
-
-      return (1. + y2*std::log(y2) / (1. - y2)) / (1 - y2);
-   }
-
-   return 0.;
-}
-
-double fth3(double y) noexcept
-{
-   const double eps = 10.*std::numeric_limits<double>::epsilon();
-   const double z2 = 1.644934066848226;
-
-   if (is_zero(y, eps)) {
-      return z2;
-   }
-
-   if (is_equal(std::abs(y), 1.0, eps)) {
-      return -9./4.;
-   }
-
-   if (!is_zero(y, eps) && !is_equal(std::abs(y), 1.0, eps)) {
-      const double y2 = sqr(y);
-      const double y4 = sqr(y2);
-      const double ly = std::log(y2);
-
-      return (-1. + 2*y2 + 2*y4)
-         *(-z2 - y2*ly + ly*std::log(1. - y2) + dilog(y2))
-         / sqr(1 - y2);
-   }
-
-   return 0.;
-}
-
 /// Delta function from hep-ph/0907.47682v1
 double delta_xyz(double x, double y, double z) noexcept
 {