FlattenedMoffatComponent.cpp

Go to the documentation of this file.

 
/*
 * FlattenedMoffatComponent.cpp
 *
 *  Created on: 2019 M01 7
 *      Author: mschefer
 */
 
#include <math.h>
 
#include "ModelFitting/Models/FlattenedMoffatComponent.h"
 
namespace ModelFitting {
 
FlattenedMoffatComponent::FlattenedMoffatComponent(std::shared_ptr<BasicParameter>  max_intensity, std::shared_ptr<BasicParameter> moffat_index,
    std::shared_ptr<BasicParameter> minkowski_distance_param, std::shared_ptr<BasicParameter> flat_top_offset) :
      m_max_intensity { max_intensity },
      m_moffat_index { moffat_index },
      m_minkowski_distance_param { minkowski_distance_param },
      m_flat_top_offset { flat_top_offset }
{
}
 
double FlattenedMoffatComponent::getValue(double x, double y) {
 
  auto z = pow(pow(fabs(x), m_minkowski_distance_param->getValue()) + pow(fabs(y), m_minkowski_distance_param->getValue()),
      1 / m_minkowski_distance_param->getValue()) -  m_flat_top_offset->getValue();
 
  if (z < 0) {
    return m_max_intensity->getValue();
  } else {
    return m_max_intensity->getValue() * pow(1+z*z, -m_moffat_index->getValue());
  }
}
 
void FlattenedMoffatComponent::updateRasterizationInfo(double, double) {
}
 
std::vector<ModelComponent::ModelSample> FlattenedMoffatComponent::getSharpSampling() {
  return {};
}
 
bool FlattenedMoffatComponent::insideSharpRegion(double, double) {
  return false;
}
 
}