CHANGE PHYSICS: use m_refined_stack instead of m_slices

Resample/Processed/ProcessedSample.cpp

@@ -253,12 +253,12 @@ const IFresnelMap* ProcessedSample::fresnelMap() const
 
 double ProcessedSample::sliceTopZ(size_t i) const
 {
-    return m_slices.at(i).zTop();
+    return m_refined_stack.at(i).zTop();
 }
 
 double ProcessedSample::sliceBottomZ(size_t i) const
 {
-    return m_slices.at(i).zBottom();
+    return m_refined_stack.at(i).zBottom();
 }
 
 bool ProcessedSample::containsMagneticMaterial() const
@@ -273,7 +273,7 @@ bool ProcessedSample::polarizing() const
 
 bool ProcessedSample::hasRoughness() const
 {
-    for (const auto& slice : m_slices)
+    for (const auto& slice : m_refined_stack)
         if (slice.topRoughness())
             return true;
     return false;
@@ -286,8 +286,8 @@ double ProcessedSample::crossCorrSpectralFun(const kvector_t kvec, size_t j, siz
         return 0.0;
     const double z_j = sliceBottomZ(j);
     const double z_k = sliceBottomZ(k);
-    const LayerRoughness* rough_j = m_slices.at(j + 1).topRoughness();
-    const LayerRoughness* rough_k = m_slices.at(k + 1).topRoughness();
+    const LayerRoughness* rough_j = m_refined_stack.at(j + 1).topRoughness();
+    const LayerRoughness* rough_k = m_refined_stack.at(k + 1).topRoughness();
     if (!rough_j || !rough_k)
         return 0.0;
     const double sigma_j = rough_j->getSigma();

Resample/Processed/ProcessedSample.h

@@ -60,7 +60,7 @@ public:
 private:
     const MultiLayer& m_sample;
     const bool m_polarized;
-    const SliceStack m_slices;
+    const SliceStack m_slices; //!< used only during construction
     const std::vector<ProcessedLayout> m_layouts;
     const SliceStack m_refined_stack;
     const std::unique_ptr<const IFresnelMap> m_fresnel_map;