vector Delta is now defined with length * cos and sin(theta_lDelta)

src/diffractivecrosssection.cpp

@@ -114,14 +114,10 @@ std::complex<double> DiffractiveCrossSection::SumAmplitudes()
 
     //TODO: loop over all values
     double length_l = 1./(100*lattice_spacing); //[GeV] This is the magnitude of the vector ... whose direction is chosen to be in the x direction (1,0) For starters << a (lattice const)
-    double vector_l[2] = {length_l, 0.}; //this is the vector l 
+    double vector_l[2] = {length_l, 0.}; 
     double theta_lDelta = M_PI;
     double length_Delta = 0.5; //[GeV] This is the magnitude of the vector ... For starters <1 Gev
-    // given the first vector's (l) components and the inner product, find the second vector's (Delta) components
-    // sign of y component unknown, choose +
-    // TODO: this is nonsense, fix!
-    double vector_Delta[2] = {length_Delta*length_l*cos(theta_lDelta), pow(length_Delta,2) - pow(length_Delta*length_l*cos(theta_lDelta),2)};
-
+    double vector_Delta[2] = {length_Delta*cos(theta_lDelta), length_Delta*sin(theta_lDelta)};
 
     // quark flavors
     const double elementary_charge = sqrt(4.*M_PI/137.036);