IsophotalFluxTask.cpp

Go to the documentation of this file.

/*
 * IsophotalTask.cpp
 *
 *  Created on: Sep 28, 2016
 *      Author: mschefer
 */
 
#include "SEImplementation/Property/PixelCoordinateList.h"
#include "SEImplementation/Plugin/DetectionFramePixelValues/DetectionFramePixelValues.h"
#include "SEImplementation/Plugin/IsophotalFlux/IsophotalFlux.h"
#include "SEImplementation/Plugin/DetectionFrameInfo/DetectionFrameInfo.h"
 
#include "SEImplementation/Plugin/IsophotalFlux/IsophotalFluxTask.h"
 
namespace SourceXtractor {
 
void IsophotalFluxTask::computeProperties(SourceInterface& source) const {
  const auto& detection_frame_info = source.getProperty<DetectionFrameInfo>();
  const auto& pixel_values = source.getProperty<DetectionFramePixelValues>().getValues();
  const auto& pixel_variances = source.getProperty<DetectionFramePixelValues>().getVariances();
 
  SeFloat total_flux = 0.0;
  SeFloat total_variance = 0.0;
 
  auto variance_iter = pixel_variances.begin();
  for (auto value : pixel_values) {
    auto variance = *variance_iter++;
 
    total_flux += value;
    total_variance += variance;
  }
 
  // Add variance from gain
  SeFloat gain = detection_frame_info.getGain();
  if (gain > 0.0) {
    total_variance += total_flux / gain;
  }
 
  auto flux_error = sqrt(total_variance);
 
  auto mag = total_flux > 0.0 ? -2.5*log10(total_flux) + m_magnitude_zero_point : std::numeric_limits<SeFloat>::quiet_NaN();
  auto mag_error = 1.0857 * flux_error / total_flux;
 
  source.setProperty<IsophotalFlux>(total_flux, flux_error, mag, mag_error);
}
 
}