migrating from ParticleLayoutContribution to ReLayout

Resample/Processed/ReLayout.cpp

@@ -14,16 +14,34 @@
 
 #include "Resample/Processed/ReLayout.h"
 #include "Resample/Coherence/CoheringSubparticles.h"
+#include "Resample/Interparticle/DecouplingApproximationStrategy.h"
 #include "Resample/Interparticle/IInterparticleStrategy.h"
+#include "Resample/Interparticle/SSCAStrategy.h"
 #include "Sample/Aggregate/IInterference.h"
+#include "Sample/Aggregate/InterferenceRadialParacrystal.h"
 
 namespace {
 
+std::unique_ptr<IInterparticleStrategy> processedInterference(const ReLayout& re_layout,
+                                                              const SimulationOptions& sim_params,
+                                                              bool polarized)
+{
+    const IInterference* iff = re_layout.interferenceFunction();
+
+    const OwningVector<const CoheringSubparticles>& weighted_formfactors = re_layout.subparticles();
+
+    const auto* radial_para = dynamic_cast<const InterferenceRadialParacrystal*>(iff);
+    if (radial_para)
+        if (double kappa = radial_para->kappa(); kappa > 0.0)
+            return std::make_unique<SSCAStrategy>(weighted_formfactors, radial_para, sim_params,
+                                                  polarized, kappa);
+
+    return std::make_unique<DecouplingApproximationStrategy>(weighted_formfactors, iff, sim_params,
+                                                             polarized);
+}
+
 } // namespace
 
-//  ************************************************************************************************
-//  class implementation
-//  ************************************************************************************************
 
 ReLayout::ReLayout(double surface_density, OwningVector<const CoheringSubparticles>&& particles,
                    const IInterference* iff)