From ce0336b98af55478977116155fb8822a1281e9e5 Mon Sep 17 00:00:00 2001 From: Alexander Voigt Date: Sat, 13 Jan 2024 16:32:39 +0100 Subject: [PATCH] remove unused functions --- .../mh2_eft/ThresholdLoopFunctions.cpp | 68 ------------------- 1 file changed, 68 deletions(-) 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::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::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::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 {