*   Nterms = 28

L   [C(1,[p],2,[p],3,[p],4)] =
       + k1.e2*k1.e3*e1.e4*ss23
       - k1.e2*k1.e4*e1.e3*ss23
       + k1.e2*k2.e3*e1.e4*ss12
       + k1.e2*k2.e3*e1.e4*ss23
       - k1.e2*k2.e4*e1.e3*ss23
       - k1.e2*k3.e1*e3.e4*ss23
       - k1.e3*k2.e1*e2.e4*ss23
       + k1.e3*k2.e4*e1.e2*ss23
       - k1.e3*k3.e2*e1.e4*ss12
       + k1.e4*k2.e1*e2.e3*ss12
       + k1.e4*k2.e1*e2.e3*ss23
       - k1.e4*k2.e3*e1.e2*ss12
       - k1.e4*k2.e3*e1.e2*ss23
       + k1.e4*k3.e2*e1.e3*ss12
       - k2.e1*k2.e3*e2.e4*ss12
       - k2.e1*k2.e3*e2.e4*ss23
       + k2.e1*k2.e4*e2.e3*ss12
       + k2.e1*k2.e4*e2.e3*ss23
       + k2.e1*k3.e2*e3.e4*ss12
       + k2.e1*k3.e2*e3.e4*ss23
       - k2.e3*k3.e1*e2.e4*ss12
       + k2.e4*k3.e1*e2.e3*ss12
       + k3.e1*k3.e2*e3.e4*ss12
       - e1.e2*e3.e4*ss12*ss23
       - e1.e2*e3.e4*ss23^2
       + e1.e3*e2.e4*ss12*ss23
       - e1.e4*e2.e3*ss12^2
       - e1.e4*e2.e3*ss12*ss23
      ;