IsoRotDiff.cpp

Go to the documentation of this file.

// Mantid Repository : https://github.com/mantidproject/mantid
//
// Copyright © 2018 ISIS Rutherford Appleton Laboratory UKRI,
//   NScD Oak Ridge National Laboratory, European Spallation Source,
//   Institut Laue - Langevin & CSNS, Institute of High Energy Physics, CAS
// SPDX - License - Identifier: GPL - 3.0 +
#include "MantidCurveFitting/Functions/IsoRotDiff.h"
#include "MantidAPI/FunctionFactory.h"
#include "MantidAPI/IFunction.h"
#include "MantidAPI/ImmutableCompositeFunction.h"
#include "MantidCurveFitting/Functions/ElasticIsoRotDiff.h"
#include "MantidCurveFitting/Functions/InelasticIsoRotDiff.h"
#include "MantidKernel/Logger.h"
#include <boost/math/special_functions/bessel.hpp>
#include <cmath>
#include <limits>
 
using Mantid::API::IFunction;
using Mantid::API::ImmutableCompositeFunction;
using Attribute = Mantid::API::IFunction::Attribute;
 
namespace {
Mantid::Kernel::Logger g_log("IsoRotDiff");
}
 
namespace Mantid::CurveFitting::Functions {
 
DECLARE_FUNCTION(IsoRotDiff)
 
 
void IsoRotDiff::trickleDownAttribute(const std::string &name) {
  for (size_t iFun = 0; iFun < this->nFunctions(); iFun++) {
    auto fun = this->getFunction(iFun);
    if (fun->hasAttribute(name)) {
      fun->setAttribute(name, this->getAttribute(name));
    }
  }
}
 
void IsoRotDiff::declareAttribute(const std::string &name, const Attribute &defaultValue) {
  ImmutableCompositeFunction::declareAttribute(name, defaultValue);
  this->trickleDownAttribute(name);
}
 
void IsoRotDiff::setAttribute(const std::string &name, const Attribute &att) {
  ImmutableCompositeFunction::setAttribute(name, att);
  this->trickleDownAttribute(name);
}
 
void IsoRotDiff::init() {
  m_elastic = std::dynamic_pointer_cast<ElasticIsoRotDiff>(
      API::FunctionFactory::Instance().createFunction("ElasticIsoRotDiff"));
  this->addFunction(m_elastic);
  m_inelastic = std::dynamic_pointer_cast<InelasticIsoRotDiff>(
      API::FunctionFactory::Instance().createFunction("InelasticIsoRotDiff"));
  this->addFunction(m_inelastic);
 
  this->declareAttribute("Q", API::IFunction::Attribute(0.3));
  this->declareAttribute("N", API::IFunction::Attribute(25));
 
  // Set the ties between Elastic and Inelastic parameters
  this->addDefaultTies("f1.Height=f0.Height,f1.Radius=f0.Radius,f1.Centre=f0.Centre");
  this->applyTies();
}
 
} // namespace Mantid::CurveFitting::Functions