From 6780a0a4379f1a297b117d5ab239f7dca703ae02 Mon Sep 17 00:00:00 2001 From: Alexander Voigt Date: Fri, 29 Sep 2017 09:40:37 +0200 Subject: [PATCH] use complex logarithm where the argument can be negative --- meta/SM/HSSUSY_corrections.m | 49 +++++++++++++++++++----------------- 1 file changed, 26 insertions(+), 23 deletions(-) diff --git a/meta/SM/HSSUSY_corrections.m b/meta/SM/HSSUSY_corrections.m index 3e0993485..6c6aab1a8 100644 --- a/meta/SM/HSSUSY_corrections.m +++ b/meta/SM/HSSUSY_corrections.m @@ -2679,10 +2679,11 @@ shiftDROS\[Alpha]t\[Alpha]Sgeneral[mtilde_, mQ3_, mU3_, QXt_, M3_] := Re@ Block[{x1 = (mQ3)^2/M3^2, x2 = (mU3)^2/M3^2, xt = QXt/M3, xs = mtilde^2/M3^2}, + ( 1/(64 \[Pi]^4) (-(( - 2 ((-1 + x1)^2 Log[1 - x1] + + 2 ((-1 + x1)^2 ComplexLog[1 - x1] + x1 (3 + x1 - (-2 + x1) Log[x1] + 2 Log[xs/x1])))/x1^2) - ( - 2 ((-1 + x2)^2 Log[1 - x2] + + 2 ((-1 + x2)^2 ComplexLog[1 - x2] + x2 (3 + x2 - (-2 + x2) Log[x2] + 2 Log[xs/x2])))/x2^2 + 1/(x1^2 x2^2) xt (1/(x1 - @@ -2692,34 +2693,34 @@ x2]))/((-1 + x1) (x1 - x2) (-1 + x2))) + 1/(x1 - x2) xt (-7 x1 + x1/(-1 + x1) + x1/(-1 + x2) + 7 x2 + x2/( - 1 - x1) + x2/(1 - x2) - (4 (-1 + x1)^2 Log[1 - x1])/ + 1 - x1) + x2/(1 - x2) - (4 (-1 + x1)^2 ComplexLog[1 - x1])/ x1 + ((-x1 + x2) Log[x1])/(-1 + x1)^2 - - 8 Log[1 - x2] + (4 Log[1 - x2])/x2 + - 4 x2 Log[1 - x2] - (x1 Log[x2])/(-1 + x2)^2 + ( + 8 ComplexLog[1 - x2] + (4 ComplexLog[1 - x2])/x2 + + 4 x2 ComplexLog[1 - x2] - (x1 Log[x2])/(-1 + x2)^2 + ( x2 Log[x2])/(-1 + x2)^2 + 9 x1 Log[x1/xs] - x2 Log[x1/xs] + x1 Log[x2/xs] - 9 x2 Log[x2/xs] + 4 x1 Log[xs] - 4 x2 Log[xs]) - ( - 4 ((-1 + x1) x2 Log[1 - x1] + - x1 ((-1 + x2) Log[1 - x2] - 2 x2 (2 + Log[xs]))))/( + 4 ((-1 + x1) x2 ComplexLog[1 - x1] + + x1 ((-1 + x2) ComplexLog[1 - x2] - 2 x2 (2 + Log[xs]))))/( x1 x2) + ( - xt ((-1 + x1)^2 Log[1 - x1] + + xt ((-1 + x1)^2 ComplexLog[1 - x1] + x1 (3 + x1 - (-2 + x1) Log[x1] + 2 Log[xs/x1])))/ x1^2 + ( - xt ((-1 + x2)^2 Log[1 - x2] + + xt ((-1 + x2)^2 ComplexLog[1 - x2] + x2 (3 + x2 - (-2 + x2) Log[x2] + 2 Log[xs/x2])))/ x2^2) + xt (1/(x1 - x2)^2 2 (-2 x1 + - 2 x2 + (x1 + x2) Log[x1/x2]) ((-1 + x1)^2 x2^2 Log[ + 2 x2 + (x1 + x2) Log[x1/x2]) ((-1 + x1)^2 x2^2 ComplexLog[ 1 - x1] - x1 ((-2 + x1) x2^2 Log[x1] + - x1 (-1 + x2)^2 Log[1 - x2] + + x1 (-1 + x2)^2 ComplexLog[1 - x2] + x2 (3 x1 - 3 x2 - x1 (-2 + x2) Log[x2] - 2 x2 Log[xs/x1] + 2 x1 Log[xs/x2]))) + x1^2 xt (1/(x1^2 (x1 - x2)^3) 2 xt (-x1 + x2 + - x1 Log[x1/x2]) ((-1 + x1)^2 x2^2 Log[1 - x1] - + x1 Log[x1/x2]) ((-1 + x1)^2 x2^2 ComplexLog[1 - x1] - x1 ((-2 + x1) x2^2 Log[x1] + - x1 (-1 + x2)^2 Log[1 - x2] + + x1 (-1 + x2)^2 ComplexLog[1 - x2] + x2 (3 x1 - 3 x2 - x1 (-2 + x2) Log[x2] - 2 x2 Log[xs/x1] + 2 x1 Log[xs/x2]))) - 1/((-1 + x1/x2)^4 x2) (-2 x1 + @@ -2731,28 +2732,30 @@ 1/(x1 - x2) xt (-7 x1 + x1/(-1 + x1) + x1/(-1 + x2) + 7 x2 + x2/(1 - x1) + x2/(1 - x2) - ( - 4 (-1 + x1)^2 Log[1 - x1])/ + 4 (-1 + x1)^2 ComplexLog[1 - x1])/ x1 + ((-x1 + x2) Log[x1])/(-1 + x1)^2 - - 8 Log[1 - x2] + (4 Log[1 - x2])/x2 + - 4 x2 Log[1 - x2] - (x1 Log[x2])/(-1 + x2)^2 + ( + 8 ComplexLog[1 - x2] + (4 ComplexLog[1 - x2])/x2 + + 4 x2 ComplexLog[1 - x2] - (x1 Log[x2])/(-1 + x2)^2 + ( x2 Log[x2])/(-1 + x2)^2 + 9 x1 Log[x1/xs] - x2 Log[x1/xs] + x1 Log[x2/xs] - 9 x2 Log[x2/xs] + 4 x1 Log[xs] - 4 x2 Log[xs]) - ( - 4 ((-1 + x1) x2 Log[1 - x1] + - x1 ((-1 + x2) Log[1 - x2] - + 4 ((-1 + x1) x2 ComplexLog[1 - x1] + + x1 ((-1 + x2) ComplexLog[1 - x2] - 2 x2 (2 + Log[xs]))))/(x1 x2) - ( - xt ((-1 + x1)^2 Log[1 - x1] + + xt ((-1 + x1)^2 ComplexLog[1 - x1] + x1 (3 + x1 - (-2 + x1) Log[x1] + 2 Log[xs/x1])))/ x1^2 + ( - 3 xt ((-1 + x2)^2 Log[1 - x2] + + 3 xt ((-1 + x2)^2 ComplexLog[1 - x2] + x2 (3 + x2 - (-2 + x2) Log[x2] + 2 Log[xs/x2])))/ x2^2) + 1/(x1 x2^3) - 6 xt ((-1 + x1)^2 x2^2 Log[1 - x1] - + 6 xt ((-1 + x1)^2 x2^2 ComplexLog[1 - x1] - x1 ((-2 + x1) x2^2 Log[x1] + - x1 (-1 + x2)^2 Log[1 - x2] + + x1 (-1 + x2)^2 ComplexLog[1 - x2] + x2 (3 x1 - 3 x2 - x1 (-2 + x2) Log[x2] - - 2 x2 Log[xs/x1] + 2 x1 Log[xs/x2]))))))))]; + 2 x2 Log[xs/x1] + 2 x1 Log[xs/x2])))))))) + ) +]; (* shift DR -> OS at O(at*as), degenerate masses, from SUSYHD 1.1 *) shiftDROS\[Alpha]t\[Alpha]Sdegenerate[mtilde_, M3_, Xt_] :=