MaskResultsPresenter.cpp

Go to the documentation of this file.

//  ************************************************************************************************
//
//  BornAgain: simulate and fit reflection and scattering
//
//! @file      GUI/coregui/Views/MaskWidgets/MaskResultsPresenter.cpp
//! @brief     Implements class MaskResultsPresenter
//!
//! @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 "GUI/coregui/Views/MaskWidgets/MaskResultsPresenter.h"
#include "Device/Detector/DetectorMask.h"
#include "Device/Detector/RegionOfInterest.h"
#include "GUI/coregui/Models/IntensityDataItem.h"
#include "GUI/coregui/Models/MaskItems.h"
#include "GUI/coregui/Models/SessionModel.h"
#include <QVBoxLayout>
 
MaskResultsPresenter::MaskResultsPresenter(QWidget* parent)
    : QObject(parent), m_interpolation_flag_backup(false)
{
}
 
void MaskResultsPresenter::setMaskContext(SessionModel* maskModel,
                                          const QModelIndex& maskContainerIndex,
                                          IntensityDataItem* intensityItem)
{
    m_maskModel = maskModel;
    m_maskContainerIndex = maskContainerIndex;
    m_intensityDataItem = intensityItem;
}
 
void MaskResultsPresenter::resetContext()
{
    setMaskContext(nullptr, QModelIndex(), nullptr);
}
 
void MaskResultsPresenter::updatePresenter(MaskEditorFlags::PresentationType presentationType)
{
    if (!m_maskContainerIndex.isValid())
        return;
 
    if (presentationType == MaskEditorFlags::MASK_PRESENTER) {
        setShowMaskMode();
    } else if (presentationType == MaskEditorFlags::MASK_EDITOR) {
        setOriginalMode();
    }
}
 
//! Update IntensityDataItem in SessionModel to represent masked areas. Corresponding
//! bins of OutputData will be put to zero.
 
void MaskResultsPresenter::setShowMaskMode()
{
    if (OutputData<double>* maskedData = createMaskPresentation()) {
        backup_data();
        m_intensityDataItem->setOutputData(maskedData);
        m_intensityDataItem->setItemValue(IntensityDataItem::P_IS_INTERPOLATED, false);
    } else {
        m_dataBackup.reset();
    }
}
 
//! Restores original state of IntensityDataItem
 
void MaskResultsPresenter::setOriginalMode()
{
    if (m_dataBackup) {
        m_intensityDataItem->setOutputData(m_dataBackup->clone());
        m_intensityDataItem->setItemValue(IntensityDataItem::P_IS_INTERPOLATED,
                                          m_interpolation_flag_backup);
    }
}
 
void MaskResultsPresenter::backup_data()
{
    m_interpolation_flag_backup =
        m_intensityDataItem->getItemValue(IntensityDataItem::P_IS_INTERPOLATED).toBool();
    m_dataBackup.reset(m_intensityDataItem->getOutputData()->clone());
}
 
//! Constructs OutputData which contains original intensity data except masked areas,
//! and areas outside of ROI, where bin content is set to zero.
 
OutputData<double>* MaskResultsPresenter::createMaskPresentation() const
{
    // Requesting mask information
    std::unique_ptr<RegionOfInterest> roi;
    DetectorMask detectorMask;
    for (int i_row = m_maskModel->rowCount(m_maskContainerIndex); i_row > 0; --i_row) {
        QModelIndex itemIndex = m_maskModel->index(i_row - 1, 0, m_maskContainerIndex);
        if (MaskItem* maskItem = dynamic_cast<MaskItem*>(m_maskModel->itemForIndex(itemIndex))) {
            if (maskItem->modelType() == "RegionOfInterest") {
                double xlow = maskItem->getItemValue(RectangleItem::P_XLOW).toDouble();
                double ylow = maskItem->getItemValue(RectangleItem::P_YLOW).toDouble();
                double xup = maskItem->getItemValue(RectangleItem::P_XUP).toDouble();
                double yup = maskItem->getItemValue(RectangleItem::P_YUP).toDouble();
                roi.reset(new RegionOfInterest(*m_intensityDataItem->getOutputData(), xlow, ylow,
                                               xup, yup));
            } else {
                std::unique_ptr<IShape2D> shape(maskItem->createShape());
                bool mask_value = maskItem->getItemValue(MaskItem::P_MASK_VALUE).toBool();
                detectorMask.addMask(*shape.get(), mask_value);
            }
        }
    }
 
    if (!detectorMask.hasMasks() && !roi)
        return 0;
 
    OutputData<double>* result = m_intensityDataItem->getOutputData()->clone();
    detectorMask.initMaskData(*result);
 
    for (size_t i = 0; i < result->getAllocatedSize(); ++i) {
        if (detectorMask.isMasked(i))
            (*result)[i] = 0.0;
        if (roi && !roi->isInROI(i))
            (*result)[i] = 0.0;
    }
 
    return result;
}