+ debmsg

Resample/Coherence/SumDWBA.cpp

@@ -19,6 +19,7 @@
 #include "Resample/Flux/MatrixFlux.h"
 #include "Resample/Flux/ScalarFlux.h"
 #include "Sample/Scattering/IFormFactor.h"
+#include <iostream>
 
 namespace {
 
@@ -86,10 +87,14 @@ complex_t SumDWBA::coherentFF(const DiffuseElement& ele) const
     // Note that the order of multiplication matters:
     // If a prefactor is 0, then theFF() won't be called.
 
-    const complex_t term_S = T_in * T_out * m_ff->theFF(q_TT);
+    std::cout <<  "SumDWBA::coherentFF call theFF with q=" << q_TT.getQ() << std::endl;
+    const complex_t result = m_ff->theFF(q_TT);
+    std::cout <<  "SumDWBA::coherentFF obtain result=" << result << std::endl << std::endl;
+    const complex_t term_S = T_in * T_out * result;
     const complex_t term_RS = R_in * T_out * m_ff->theFF(q_RT);
     const complex_t term_SR = T_in * R_out * m_ff->theFF(q_TR);
     const complex_t term_RSR = R_in * R_out * m_ff->theFF(q_RR);
+    std::cout <<  "SumDWBA::coherentFF done with the 3 other terms" << std::endl << std::endl;
 
     return term_S + term_RS + term_SR + term_RSR;
 }