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

name = "par"

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

U = x[1]*x[3] + x[1]*x[4] + x[2]*x[3] + x[2]*x[4] + x[3]*x[4]
F = s*x[1]*x[2]*x[3] + s*x[1]*x[2]*x[4] + M[4]*x[1]*x[3]*x[4] + M[3]*x[2]*x[3]*x[4]
parameters = [M[3], M[4], s]
variables = [x[1], x[2], x[3], x[4]]

χ_generic = 4
f_vector = [9, 24, 24, 9]

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

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

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

D[2] = M[3]^2 - 2*M[3]*M[4] - 2*M[3]*s + M[4]^2 - 2*M[4]*s + s^2
χ[2] = 3
weights[2] = [[0, 0, -1, 0], [0, 0, 0, -1], [0, 0, 0, 0]]
computed_with[2] = ["PLD_sym", "HyperInt"]

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

D[3] = M[4]
χ[3] = 2
weights[3] = [[-2, 0, -2, -2], [-1, 0, 0, 0]]
computed_with[3] = ["PLD_sym"]

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

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