################################
# Diagram information
################################

name = "B4"

edges = [[1, 2], [1, 2], [1, 2], [1, 2]]
nodes = [1, 1, 2, 2]
internal_masses = [m2, m2, m2, m2]
external_masses = [M[1], M[2], M[3], M[4]]

U = x[1]*x[2]*x[3] + x[1]*x[2]*x[4] + x[1]*x[3]*x[4] + x[2]*x[3]*x[4]
F = -m2*x[1]^2*x[2]*x[3] - m2*x[1]^2*x[2]*x[4] - m2*x[1]^2*x[3]*x[4] - m2*x[1]*x[2]^2*x[3] - m2*x[1]*x[2]^2*x[4] - m2*x[1]*x[2]*x[3]^2 + (-4*m2 + s)*x[1]*x[2]*x[3]*x[4] - m2*x[1]*x[2]*x[4]^2 - m2*x[1]*x[3]^2*x[4] - m2*x[1]*x[3]*x[4]^2 - m2*x[2]^2*x[3]*x[4] - m2*x[2]*x[3]^2*x[4] - m2*x[2]*x[3]*x[4]^2
parameters = [m2, s]
variables = [x[1], x[2], x[3], x[4]]

χ_generic = 15
f_vector = [16, 42, 42, 16]

################################
# Component 1
################################

D[1] = m2
χ[1] = 1
weights[1] = [[-1, 1, 1, 3], [-1, 1, 3, 1], [-1, 3, 1, 1], [1, -1, 1, 3], [1, -1, 3, 1], [1, 1, -1, 3], [1, 1, 3, -1], [1, 3, -1, 1], [1, 3, 1, -1], [3, -1, 1, 1], [3, 1, -1, 1], [3, 1, 1, -1], [-1, 0, 1, 1], [-1, 1, 0, 1], [-1, -1, 0, 1], [-1, 0, -1, 1], [-1, 1, 1, 0], [-1, -1, 1, 0], [-1, 0, 1, -1], [-1, 1, -1, 0], [-1, 1, 0, -1], [0, -1, 1, 1], [0, -1, -1, 1], [1, -1, 0, 1], [0, -1, 1, -1], [1, -1, 1, 0], [0, 1, -1, 1], [1, 0, -1, 1], [0, 1, 1, -1], [1, 0, 1, -1], [0, 1, -1, -1], [1, 1, -1, 0], [1, 1, 0, -1], [1, -1, -1, 0], [1, -1, 0, -1], [1, 0, -1, -1], [-1, 0, 0, 0], [-1, -1, 0, 0], [-1, 0, -1, 0], [-1, -1, -1, 0], [-1, 0, 0, -1], [-1, -1, 0, -1], [-1, 0, -1, -1], [0, -1, 0, 0], [0, -1, -1, 0], [0, -1, 0, -1], [0, 0, -1, 0], [0, 0, 0, -1], [0, 0, -1, -1], [0, -1, -1, -1]]
computed_with[1] = ["PLD_sym"]

################################
# Component 2
################################

D[2] = m2 - 1//16*s
χ[2] = 14
weights[2] = [[-1, -1, -1, -1]]
computed_with[2] = ["PLD_sym", "HyperInt"]

################################
# Component 3
################################

D[3] = m2 - 1//4*s
χ[3] = 11
weights[3] = [[-1, -1, -1, -1], [-1, -1, -1, -1], [-1, -1, -1, -1], [-1, -1, -1, -1]]
computed_with[3] = ["PLD_sym", "HyperInt"]

################################
# Component 4
################################

D[4] = s
χ[4] = 5
weights[4] = [[-1, -1, -1, -1], [-1, -1, -1, -1], [-1, -1, -1, -1], [0, 0, 0, 0]]
computed_with[4] = ["PLD_sym"]