Skip to content
Snippets Groups Projects
Commit 8ac72bac authored by Van Herck, Walter's avatar Van Herck, Walter
Browse files

Use the dummy DW factor in the SSCA approximation

parent d38f6841
No related branches found
No related tags found
No related merge requests found
......@@ -14,6 +14,7 @@
#include "SSCApproximationStrategy.h"
#include "FormFactorCoherentSum.h"
#include "IInterferenceFunction.h"
#include "InterferenceFunctionUtils.h"
#include "SimulationElement.h"
......@@ -49,7 +50,8 @@ double SSCApproximationStrategy::scalarCalculation(const SimulationElement& sim_
complex_t p2kappa = m_helper.getCharacteristicSizeCoupling(qp, m_formfactor_wrappers);
complex_t omega = m_helper.getCharacteristicDistribution(qp, mP_iff.get());
double interference_intensity = 2.0 * (mean_ff_norm * omega / (1.0 - p2kappa * omega)).real();
return diffuse_intensity + interference_intensity;
double dw_factor = mP_iff->DWfactor(sim_element.getMeanQ());
return diffuse_intensity + dw_factor * interference_intensity;
}
//! This is the polarized version
......@@ -75,5 +77,6 @@ double SSCApproximationStrategy::polarizedCalculation(const SimulationElement& s
Eigen::Matrix2cd diffuse_matrix2 = polarization_handler.getAnalyzerOperator() * diffuse_matrix;
double interference_trace = std::abs(interference_matrix.trace());
double diffuse_trace = std::abs(diffuse_matrix2.trace());
return diffuse_trace + interference_trace;
double dw_factor = mP_iff->DWfactor(sim_element.getMeanQ());
return diffuse_trace + dw_factor * interference_trace;
}
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment