MoffatModelEvaluator.cpp

Go to the documentation of this file.

/*
 * MoffatModelFittingUtils.cpp
 *
 *  Created on: 2019 M02 20
 *      Author: mschefer
 */
 
 
#include "AlexandriaKernel/memory_tools.h"
#include "ModelFitting/Parameters/ManualParameter.h"
#include "ModelFitting/Models/ExtendedModel.h"
#include "ModelFitting/Models/FlattenedMoffatComponent.h"
 
#include "SEImplementation/Image/ImageInterfaceTraits.h"
 
#include "SEImplementation/Plugin/MoffatModelFitting/MoffatModelFitting.h"
#include "SEImplementation/Plugin/MoffatModelFitting/MoffatModelEvaluator.h"
 
namespace SourceXtractor {
 
using namespace ModelFitting;
using Euclid::make_unique;
 
MoffatModelEvaluator::MoffatModelEvaluator(const MoffatModelFitting& model) {
  m_iterations = model.getIterations();
 
  auto x = std::make_shared<ManualParameter>(model.getX());
  auto y = std::make_shared<ManualParameter>(model.getY());
  auto moffat_i0 = std::make_shared<ManualParameter>(model.getMoffatI0());
  auto moffat_index = std::make_shared<ManualParameter>(model.getMoffatIndex());
  auto minkowski_exponent = std::make_shared<ManualParameter>(model.getMinkowksiExponent());
  auto flat_top_offset = std::make_shared<ManualParameter>(model.getTopOffset());
  double size = model.getSize();
  auto x_scale = std::make_shared<ManualParameter>(model.getXScale());
  auto y_scale= std::make_shared<ManualParameter>(model.getYScale());
  auto moffat_rotation= std::make_shared<ManualParameter>(model.getMoffatRotation());
 
  // Moffat model
  std::vector<std::unique_ptr<ModelComponent>> component_list {};
  auto moff = Euclid::make_unique<FlattenedMoffatComponent>(moffat_i0, moffat_index, minkowski_exponent, flat_top_offset);
  component_list.clear();
  component_list.emplace_back(std::move(moff));
 
  m_model = std::make_shared<ExtendedModel<ImageInterfaceTypePtr>>(
      std::move(component_list), x_scale, y_scale, moffat_rotation, size, size, x, y);
}
 
 
}