SpecularComputation.cpp

Go to the documentation of this file.

//  ************************************************************************************************
//
//  BornAgain: simulate and fit reflection and scattering
//
//! @file      Sim/Computation/SpecularComputation.cpp
//! @brief     Implements class SpecularComputation.
//!
//! @homepage  http://www.bornagainproject.org
//! @license   GNU General Public License v3 or higher (see COPYING)
//! @copyright Forschungszentrum Jülich GmbH 2018
//! @authors   Scientific Computing Group at MLZ (see CITATION, AUTHORS)
//
//  ************************************************************************************************
 
#include "Sim/Computation/SpecularComputation.h"
#include "Base/Progress/ProgressHandler.h"
#include "Resample/Element/SpecularElement.h"
#include "Resample/Flux/IFlux.h"
#include "Resample/Processed/ReSample.h"
#include "Resample/Specular/ComputeFluxMagnetic.h"
#include "Resample/Specular/ComputeFluxScalar.h"
 
SpecularComputation::SpecularComputation(const reSample& re_sample,
                                         const SimulationOptions& options,
                                         ProgressHandler& progress, SpecularElementIter begin_it,
                                         SpecularElementIter end_it)
    : IComputation(std::move(re_sample), options, progress)
    , m_begin_it(begin_it)
    , m_end_it(end_it)
{
}
 
SpecularComputation::~SpecularComputation() = default;
 
void SpecularComputation::runProtected()
{
    const SliceStack& slices = m_re_sample.averageSlices();
 
    for (auto it = m_begin_it; it != m_end_it; ++it) {
        SpecularElement& ele = *it;
        if (!ele.isCalculated())
            continue;
        std::vector<complex_t> kz_stack = ele.produceKz(slices);
        if (m_re_sample.polarizing()) {
            const SpinMatrix R = Compute::SpecularMagnetic::topLayerR(slices, kz_stack, true);
 
            const auto& polMatrix = ele.polMatrices().polarizerMatrix();
            const auto& anaMatrix = ele.polMatrices().analyzerMatrix();
 
            const complex_t trace = (polMatrix * R.adjoint() * anaMatrix * R).trace();
 
            ele.setIntensity(std::abs(trace));
        } else {
            const complex_t R = Compute::SpecularScalar::topLayerR(slices, kz_stack);
            ele.setIntensity(std::norm(R));
        }
 
        stepProgress();
    }
}