WIP, see TODO in SumAmplitudes

src/ipglasma.cpp

@@ -645,11 +645,19 @@ double IPGlasma::SumAmplitudes()
     
                     int q2[2] = {k,l};
                     double r2 = DistanceToOrigin(q2);
-                }
-                   r = sqrt(pow((X(i)-X(k)),2)+pow(Y(j)-Y(l)),2);
-
-                   double WF = WFGammaqqbarLongitudinal(double r, double z, double Q, struct* quark_flavor); 
-                   sum += Amplitude(0, quark, q2)*WF*WF;
+                    double r = sqrt(pow((X(i)-X(k)),2)+pow(Y(j)-Y(l)),2);
+
+                    double WF = WFGammaqqbarLongitudinal(double r, double z, double Q, struct* quark_flavor); 
+                    sum += Amplitude(0, quark, q2)*WF;
+
+                    //TODO: Finish eq (12) from notes
+                    //          Finish amplitude
+                    //              Add phase factor
+                    //          Add z0 z1 integrations
+                    //          Add l Delta integrations
+                    //          Add delta functions
+                    //          Add prefactor and integration measure
+                    //      
                 }
             }
         }