################################ # Diagram information ################################ name = "Bhabha-npl-dbox" edges = [[1, 2], [2, 5], [3, 5], [4, 5], [4, 6], [1, 6], [3, 6]] nodes = [1, 2, 3, 4] internal_masses = [0, Mb2, Mb2, 0, Ma2, Ma2, 0] external_masses = [Ma2, Mb2, Mb2, Ma2] U = x[1]*x[3] + x[1]*x[4] + x[1]*x[5] + x[1]*x[7] + x[2]*x[3] + x[2]*x[4] + x[2]*x[5] + x[2]*x[7] + x[3]*x[4] + x[3]*x[5] + x[3]*x[6] + x[4]*x[6] + x[4]*x[7] + x[5]*x[6] + x[5]*x[7] + x[6]*x[7] F = -Mb2*x[1]*x[3]^2 + (Ma2 + Mb2 - s - t)*x[1]*x[3]*x[5] + t*x[1]*x[4]*x[7] - Ma2*x[1]*x[5]^2 - Mb2*x[2]^2*x[3] - Mb2*x[2]^2*x[4] - Mb2*x[2]^2*x[5] - Mb2*x[2]^2*x[7] - Mb2*x[2]*x[3]^2 - 2*Mb2*x[2]*x[3]*x[4] - 2*Mb2*x[2]*x[3]*x[5] + (-Ma2 - Mb2 + s)*x[2]*x[3]*x[6] + (-Ma2 - Mb2 + s)*x[2]*x[4]*x[6] - Ma2*x[2]*x[5]^2 + (-Ma2 - Mb2 + s)*x[2]*x[5]*x[6] + (-Ma2 - Mb2 + s)*x[2]*x[5]*x[7] + (-Ma2 - Mb2 + s)*x[2]*x[6]*x[7] - Mb2*x[3]^2*x[4] - Mb2*x[3]^2*x[5] - Mb2*x[3]^2*x[6] + (-Ma2 - Mb2 + s)*x[3]*x[4]*x[6] - Ma2*x[3]*x[5]^2 - 2*Ma2*x[3]*x[5]*x[6] - Ma2*x[3]*x[6]^2 - Ma2*x[4]*x[6]^2 - Ma2*x[5]^2*x[6] - Ma2*x[5]^2*x[7] - Ma2*x[5]*x[6]^2 - 2*Ma2*x[5]*x[6]*x[7] - Ma2*x[6]^2*x[7] parameters = [Ma2, Mb2, s, t] variables = [x[1], x[2], x[3], x[4], x[5], x[6], x[7]] χ_generic = 111 f_vector = [35, 184, 402, 457, 291, 103, 18] ################################ # Component 1 ################################ D[1] = Ma2 χ[1] = 51 weights[1] = [[-6, -1, 0, -4, -5, -1, 0], [2, 0, 2, -2, -3, 2, 2], [2, 2, 0, -2, -3, 2, 2], [-2, -1, 0, 2, 2, -3, 2], [-4, -4, -1, -8, -2, -6, -2], [2, 2, 2, -2, -3, 0, 2], [0, -1, -3, -6, -6, -1, -10], [-2, -1, 2, 2, 0, -3, 2], [-4, -4, -3, -2, -1, -6, -10], [0, 0, 2, -2, -3, 2, 2], [0, 2, 0, -2, -3, 2, 2], [-6, -2, -1, -6, -4, -4, 0], [0, 2, 2, -2, -3, 0, 2], [-2, 1, 1, -2, -4, 1, -2], [-4, -1, 0, 0, -2, -3, 0], [-6, -2, -1, 0, -3, -4, -4], [2, 0, 0, -2, -3, 2, 2], [2, 0, 2, -2, -3, 0, 2], [2, 0, 2, -2, -3, 2, 0], [2, 2, 0, -2, -3, 0, 2], [2, 2, 0, -2, -3, 2, 0], [-2, -1, 0, 0, 2, -3, 2], [-2, -1, 0, 2, 0, -3, 2], [-2, -1, 0, 2, 2, -3, 0], [0, 0, 0, -4, -3, -2, 0], [0, -2, -2, -8, -5, -4, -6], [-2, -1, 2, 0, 0, -3, 2], [-2, -3, 0, -4, 0, -5, -4], [2, 2, 2, -2, -3, 0, 0], [0, 0, 0, 0, -1, -1, 0], [-2, -3, -5, -4, -5, -5, -12], [-2, -1, 2, 2, 0, -3, 0], [-4, -3, -3, -6, -4, -1, 0], [0, 0, 0, -2, -3, 2, 2], [0, 0, 2, -2, -3, 0, 2], [0, 0, 2, -2, -3, 2, 0], [0, 2, 0, -2, -3, 0, 2], [0, 2, 0, -2, -3, 2, 0], [-2, 0, 0, -4, -3, -2, 0], [-4, -1, 0, -2, -2, -3, 0], [0, 2, 2, -2, -3, 0, 0], [-1, 0, 0, 0, -1, -1, 0], [-4, -1, -1, 0, -3, -3, -4], [-4, -1, 0, 0, -2, -3, -2], [0, -2, -2, -4, -3, 0, 0], [0, -4, -4, -8, -5, -2, -6], [2, 0, 0, -2, -3, 0, 2], [2, 0, 0, -2, -3, 2, 0], [2, 0, 2, -2, -3, 0, 0], [2, 2, 0, -2, -3, 0, 0], [-2, -1, 0, 0, 0, -3, 2], [-2, -1, 0, 0, 2, -3, 0], [-2, -1, 0, 2, 0, -3, 0], [0, 0, 0, -4, -3, -2, -2], [0, 0, 0, -1, -1, -1, 0], [0, -1, -1, -3, -2, -2, -3], [-2, -1, 2, 0, 0, -3, 0], [0, 0, 0, 0, -1, -1, -1], [0, 1, 0, 1, 1, 0, 1], [0, 1, -1, 1, 0, 1, -2], [0, 1, -1, 1, 1, 0, -2], [-2, -1, 2, 0, 0, 1, 2], [-2, -1, 2, 2, 0, 1, 0], [-4, -3, 0, -2, -2, -3, 0], [-2, -2, -2, -4, -3, 0, 0], [-4, -3, 0, 0, -2, -3, -2], [0, 0, 0, -2, -3, 0, 2], [0, 0, 0, -2, -3, 2, 0], [0, 0, 2, -2, -3, 0, 0], [0, 2, 0, -2, -3, 0, 0], [-1, 0, 0, -1, -1, -1, 0], [-1, 0, 0, 0, -1, -1, -1], [1, 0, 0, 1, 0, 1, 1], [1, 0, 1, 0, 0, 1, 1], [0, -1, 0, -1, -1, -1, 0], [0, -2, -2, -4, -3, 0, -2], [2, 0, 0, -2, -3, 0, 0], [1, 1, 0, 0, 0, 1, 1], [1, 1, -1, 0, 0, 1, -2], [1, 1, -1, 0, 1, 0, -2], [-2, -1, 0, 0, 0, -3, 0], [0, 0, 0, -1, -1, -1, -1]] computed_with[1] = ["PLD_sym"] ################################ # Component 2 ################################ D[2] = Ma2 + Mb2 - 1//2*s - 1//2*t χ[2] = 106 weights[2] = [[-2, 1, 0, 2, 0, 1, 2], [-4, -3, -2, 0, -2, -1, 0], [-1, 0, -1, 0, -1, 0, 0], [-4, -1, -2, 0, -2, -3, 0], [-2, -1, 0, 2, 0, 1, 2], [0, 1, 0, 1, 0, 1, 1], [-2, 1, 0, 2, 0, -1, 2], [-4, -3, -2, 0, -2, -3, 0]] computed_with[2] = ["PLD_sym"] ################################ # Component 3 ################################ D[3] = Ma2 - 1//4*t χ[3] = 102 weights[3] = [[-1, 0, -2, -1, -1, 0, -3], [-1, 0, -2, -1, 0, -1, -3], [-2, 1, 2, 0, 0, 1, 0], [0, 2, 2, -2, 1, 0, 0], [-4, -3, 0, -2, -2, -1, -2], [-2, 0, 0, -4, -3, -2, -2], [-4, -1, 0, -2, -2, -3, -2], [0, 1, 0, 0, 0, 1, -2], [0, 1, 0, 0, 1, 0, -2]] computed_with[3] = ["PLD_sym"] ################################ # Component 4 ################################ D[4] = Ma2 - Mb2 χ[4] = 107 weights[4] = [[0, -1, -1, 0, -1, 0, 0], [0, 0, -1, 0, -1, -1, 0], [0, -1, -1, 0, -1, -1, 0]] computed_with[4] = ["PLD_sym", "PLD_num"] ################################ # Component 5 ################################ D[5] = Ma2^2 - 2*Ma2*Mb2 - 2*Ma2*s + Mb2^2 - 2*Mb2*s + s^2 χ[5] = 63 weights[5] = [[-2, -3, 0, 2, 2, -3, 2], [-4, -6, -1, -8, -2, -6, -2], [-2, -3, 2, 2, 0, -3, 2], [-2, -5, -5, -4, -5, -3, -12], [-4, -6, -3, -2, -1, -6, -10], [-2, -3, -5, -10, -4, -5, -4], [-2, -3, 0, 0, 2, -3, 2], [-2, -3, 0, 2, 0, -3, 2], [-2, -3, 0, 2, 2, -3, 0], [-2, -3, 2, 0, 0, -3, 2], [-2, -5, 0, -4, 0, -5, -4], [-2, -5, -5, -10, -4, -5, -4], [-2, -3, 2, 2, 0, -3, 0], [-2, -5, -5, -4, -5, -5, -12], [-2, -3, 0, 0, 0, -3, 2], [-2, -3, 0, 0, 2, -3, 0], [0, -1, -1, 0, -1, 0, -1], [-2, -3, 0, 2, 0, -3, 0], [0, -1, -1, 0, -1, -1, 0], [-2, -3, 2, 0, 0, -3, 0], [0, -2, -2, -3, -2, -1, -3], [0, 0, -1, -1, -1, -1, 0], [0, -1, -2, -3, -2, -2, -3], [0, 1, 0, 0, 1, 0, 1], [1, 0, 0, 0, 0, 0, 0]] computed_with[5] = ["PLD_sym", "PLD_num"] ################################ # Component 6 ################################ D[6] = Ma2^2 - 2*Ma2*Mb2 - 2*Ma2*s + Mb2^2 - 2*Mb2*s + s^2 + s*t χ[6] = 80 weights[6] = [[-2, 1, 0, 0, 0, 1, 0], [0, 1, 0, 0, 0, 0, 0], [0, 0, 0, 0, 1, 0, 0], [0, 0, 0, 0, 0, 0, 1], [0, 0, 0, 1, 0, 0, 0], [-1, -1, -1, -1, -1, -1, -1], [0, 0, -1, 0, 0, 0, -2], [0, 0, 1, 0, 0, 0, 0], [0, 0, 0, -2, -1, 0, 0], [0, 0, 0, 0, 0, 0, 0]] computed_with[6] = ["PLD_sym", "PLD_num"] ################################ # Component 7 ################################ D[7] = Ma2^2 - 2*Ma2*Mb2 - 2*Ma2*s - 2*Ma2*t + Mb2^2 - 2*Mb2*s - 2*Mb2*t + s^2 + 2*s*t + t^2 χ[7] = 89 weights[7] = [[-8, -3, -4, -2, -4, -3, -2], [-1, -1, -1, 0, -1, 0, 0], [-4, -3, -3, -4, -3, -1, 0], [-1, 0, -1, 0, -1, -1, 0], [-4, -1, -3, 0, -3, -3, -4], [-2, -1, 0, 0, 0, 1, 2], [0, 1, 0, 0, 0, 1, 1], [-2, 1, 0, 0, 0, -1, 2], [-2, -1, 0, 2, 0, 1, 0], [0, 1, 0, 1, 0, 1, 0], [-2, 1, 0, 2, 0, -1, 0]] computed_with[7] = ["PLD_sym", "PLD_num"] ################################ # Component 8 ################################ D[8] = Mb2 χ[8] = 51 weights[8] = [[-6, -1, -6, 0, 0, -1, -6], [-2, -3, 0, 2, 2, -1, 2], [-4, -6, -1, -8, -2, -4, -2], [-2, -3, 2, 2, 0, -1, 2], [-4, -6, -3, -2, -1, -4, -10], [1, 0, -2, 1, 1, 1, -2], [0, -1, -7, -8, -2, -1, -8], [1, 1, -2, 1, 0, 1, -2], [1, 1, -2, 1, 1, 0, -2], [-4, -3, -2, 0, 0, -1, 0], [-6, -4, -3, -4, -1, -2, 0], [-6, -4, -5, 0, -1, -2, -8], [0, 0, -2, 1, 1, 1, -2], [-2, 1, -5, -2, 1, 1, -4], [0, 1, -2, 1, 0, 1, -2], [0, 1, -2, 1, 1, 0, -2], [-2, -3, 0, 0, 2, -1, 2], [-2, -3, 0, 2, 0, -1, 2], [-2, -3, 0, 2, 2, -1, 0], [0, -1, -1, 0, 0, 0, 0], [-2, -3, 2, 0, 0, -1, 2], [-2, -5, 0, -4, 0, -3, -4], [-2, -5, -5, -10, -4, -3, -4], [-2, -3, 2, 2, 0, -1, 0], [0, -1, -2, 0, 0, 0, -3], [0, -2, -3, -3, -1, -1, -5], [1, 0, -2, 0, 1, 1, -2], [1, 0, -2, 1, 0, 1, -2], [1, 0, -2, 1, 1, 0, -2], [1, 1, -2, 0, 0, 1, -2], [1, 1, -2, 0, 1, 0, -2], [1, 1, -2, 1, 0, 0, -2], [-4, -3, -2, -2, 0, -1, 0], [-4, -3, -2, 0, 0, -1, -2], [-1, -1, -1, 0, 0, 0, 0], [-4, -3, -3, -4, -1, -1, 0], [-1, -1, -2, 0, 0, 0, -3], [0, 0, -2, 0, 1, 1, -2], [0, 0, -2, 1, 0, 1, -2], [0, 0, -2, 1, 1, 0, -2], [0, 1, -2, 0, 0, 1, -2], [0, 1, -2, 0, 1, 0, -2], [0, 1, -2, 1, 0, 0, -2], [-4, -1, -5, 0, -3, -3, -8], [-2, -3, 0, 0, 0, -1, 2], [-2, -3, 0, 0, 2, -1, 0], [0, -1, -1, -1, 0, 0, 0], [-2, -3, 0, 2, 0, -1, 0], [0, -1, -1, 0, 0, 0, -1], [-2, -3, 2, 0, 0, -1, 0], [0, -2, -2, -3, -1, -1, -3], [0, -1, -2, -1, 0, 0, -3], [1, 0, -2, 0, 0, 1, -2], [1, 0, -2, 0, 1, 0, -2], [1, 0, -2, 1, 0, 0, -2], [1, 1, -2, 0, 0, 0, -2], [0, -1, -3, -3, -2, -2, -5], [0, 0, -2, 0, -1, -1, -3], [-2, 1, 0, 0, 2, -1, 2], [-2, 1, 0, 2, 2, -1, 0], [-1, -1, -1, -1, 0, 0, 0], [-4, -3, -2, -2, 0, -3, 0], [-1, -1, -1, 0, 0, 0, -1], [-4, -3, -2, 0, 0, -3, -2], [0, 0, -2, 0, 0, 1, -2], [0, 0, -2, 0, 1, 0, -2], [0, 0, -2, 1, 0, 0, -2], [0, 1, -2, 0, 0, 0, -2], [-1, 0, -2, 0, -1, -1, -3], [0, 0, 2, -2, -1, 2, 2], [0, 2, 0, -2, -1, 2, 2], [0, 0, 1, 1, 0, 1, 1], [-2, -3, 0, 0, 0, -1, 0], [0, -1, -1, -1, 0, 0, -1], [0, -1, -1, 0, 0, -1, -1], [2, 0, 2, -2, -1, 2, 0], [1, 0, -2, 0, 0, 0, -2], [0, 0, -2, -1, -1, -1, -3], [2, 2, 0, -2, -1, 2, 0], [1, 1, 0, 1, 0, 0, 1], [1, 1, 0, 1, 0, 1, 0], [1, 1, 0, 1, 1, 0, 0], [-2, -1, 0, 0, 2, -1, 2]] computed_with[8] = ["PLD_sym", "PLD_num"] ################################ # Component 9 ################################ D[9] = Mb2 - 1//4*t χ[9] = 102 weights[9] = [[-2, 1, 0, 0, 2, 1, 0], [-4, -3, -2, -2, 0, -1, -2], [-1, -1, -2, -1, 0, 0, -3], [-4, -1, -2, -2, 0, -3, -2], [0, 0, 2, -2, 1, 2, 0], [0, 2, 0, -2, 1, 2, 0], [-2, -2, 0, -4, -3, 0, -2], [-2, 0, -2, -4, -3, 0, -2], [0, 0, 0, 0, 1, 1, -2], [-2, -1, 0, 0, 2, 1, 0], [0, 1, 0, 0, 1, 1, 0], [-2, 1, 0, 0, 2, -1, 0]] computed_with[9] = ["PLD_sym", "PLD_num"] ################################ # Component 10 ################################ D[10] = s χ[10] = 98 weights[10] = [[0, -1, -1, 0, 0, -1, 0], [1, 0, 1, 2, 2, 0, 2], [2, 0, 4, 0, 3, 0, 4], [0, -2, 0, -4, -3, -2, 0], [0, -1, 0, 0, -1, 0, -1], [0, -1, 0, 0, -1, -1, 0], [1, 0, 2, 2, 1, 0, 2], [0, -1, -2, 0, -1, 0, -3], [0, -1, -2, 0, 0, -1, -3], [1, 0, 0, 1, 0, 1, -2], [1, 0, 1, 2, 2, 0, -1], [0, 0, -1, -1, 0, -1, 0], [2, 2, 0, -2, 1, 0, 2], [0, 0, -2, -4, -3, -2, 0], [0, -1, -1, -1, 0, -1, 0], [1, 0, 1, 1, 2, 0, 2], [1, 0, 1, 2, 1, 0, 2], [0, -1, -1, 0, 0, -1, -1], [1, 0, 1, 2, 2, 0, 1], [1, 0, 0, 1, 1, 0, 1], [1, 0, 2, 1, 1, 0, 2], [0, -1, 0, -1, -1, -1, 0], [1, 0, 2, 1, 2, 0, 1], [2, 0, 2, -2, -1, 0, 2], [2, 0, 0, -2, 1, 0, 2], [0, -2, -2, -4, -3, -2, 0], [1, 0, 1, 1, 0, 1, 0], [0, -1, 0, 0, -1, -1, -1], [1, 0, 2, 2, 1, 0, 1], [1, 0, 1, 1, 0, 0, 1], [1, 0, -1, 1, 0, 1, -2], [1, 0, -1, 1, 1, 0, -2], [0, -1, -2, 0, -1, -1, -3], [1, 0, 0, 1, 0, 0, -2], [1, 1, 0, 0, 1, 0, 1], [2, 2, 0, -2, -1, 0, 2]] computed_with[10] = ["PLD_sym", "PLD_num"] ################################ # Component 11 ################################ D[11] = t χ[11] = 55 weights[11] = [[-8, -5, -6, -10, -5, -5, -12], [-4, -5, -5, -10, -4, -3, -6], [-4, -3, -5, -6, -5, -5, -12], [0, 1, -1, 0, 1, 1, -2], [0, 2, 2, -2, -1, 2, 0], [-2, -1, 2, 0, 2, 1, 0], [-2, 1, 2, 0, 2, -1, 0], [-4, -5, -2, -6, -2, -5, -6], [0, 1, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 1], [0, 0, 0, 1, 0, 0, 0], [-1, -1, -1, -1, -1, -1, -1], [0, 0, 0, 0, 0, 0, 0]] computed_with[11] = ["PLD_sym", "PLD_num"]