BeamWavelengthItem.cpp

Go to the documentation of this file.

//  ************************************************************************************************
//
//  BornAgain: simulate and fit reflection and scattering
//
//! @file      GUI/Model/Device/BeamWavelengthItem.cpp
//! @brief     Implements class BeamWavelengthItem
//!
//! @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/Model/Device/BeamWavelengthItem.h"
#include "GUI/Model/CatDevice/DistributionItemCatalog.h"
#include "GUI/Support/XML/Serialize.h"
 
namespace {
 
const double default_wl = 0.1;
 
void initDistribution(DistributionItem* newDistribution,
                      const DistributionItem* currentDistribution)
{
    newDistribution->setLimits(RealLimits::positive());
    newDistribution->setUnit(Unit::nanometer);
 
    double meanVal = default_wl;
    if (auto* cd = dynamic_cast<const SymmetricResolutionItem*>(currentDistribution))
        meanVal = cd->mean();
 
    if (auto* nd = dynamic_cast<SymmetricResolutionItem*>(newDistribution)) {
        nd->setMean(meanVal);
        nd->setMeanDecimals(4);
    }
 
    newDistribution->initDistribution(meanVal);
};
 
} // namespace
 
BeamWavelengthItem::BeamWavelengthItem()
    : BeamWavelengthItem(false)
{
}
 
BeamWavelengthItem::BeamWavelengthItem(bool onlySymmetricDistributions)
{
    m_distribution.initWithInitializer<DistributionItemCatalog>(
        "Distribution", "", "distribution",
        onlySymmetricDistributions ? DistributionItemCatalog::symmetricTypes()
                                   : DistributionItemCatalog::types(),
        initDistribution);
}
 
void BeamWavelengthItem::serialize(Streamer& s)
{
    s.assertVersion(0);
    Serialize::rwSelected<DistributionItemCatalog>(s, m_distribution);
 
    if (s.xmlReader()) {
        m_distribution->setLimits(RealLimits::positive());
        m_distribution->setUnit(Unit::nanometer);
    }
}
 
double BeamWavelengthItem::wavelength() const
{
    return BeamDistributionItem::meanValue();
}
 
// ------------------------------------------------------------------------------------------------
 
SpecularBeamWavelengthItem::SpecularBeamWavelengthItem()
    : BeamWavelengthItem(true)
{
}
 
void SpecularBeamWavelengthItem::setToRange(const RealLimits& limits)
{
    if (auto* symmetricDistribution = dynamic_cast<SymmetricResolutionItem*>(distribution())) {
        if (!limits.isInRange(wavelength())) {
            const double new_value =
                limits.isLimited() ? (limits.upperLimit() - limits.lowerLimit()) / 2. : default_wl;
            symmetricDistribution->setMean(new_value);
        }
        symmetricDistribution->setLimits(limits);
    }
}