Detector.cc

Go to the documentation of this file.

#include "mutils/nmath.h"

#include "Detector.h"
#include "Graybody.h"

namespace mu = mutils;

namespace manticore {

namespace detector {

extern "C" {

double fWidth(double nu, void *d)
{ return static_cast<Detector*>(d)->response(nu)/nu; }

}

}  // namescape detector


void Detector::init()
{
  const char *const fn = "Detector::init";

  // Effective band limits.
  freqRange_ = {resp_.front().first, resp_.back().first};
  MU_DEBUG(fn, "%s response range is [%.2f, %.2f] microns", name(),
           1e4*c_light/freqRange_.second, 1e4*c_light/freqRange_.first);

  // Initialize response interpolator.
  spline_ = initSpline(resp_);

  // Compute and set reference bandwidth.
  calcWidth();
}


double Detector::response(double nu, gsl_interp_accel *acc) const
{
  if (nu >= freqRange_.first && nu <= freqRange_.second) {
    return gsl_spline_eval(spline_, nu, acc);
  } else {
    // Truncate respone at specified band edges.
    return 0.0;
  }
}


void Detector::calcWidth()
{
  const char *const fn = "Detector::calcWidth";
  width_ = center_ * mu::kgauss_r(detector::fWidth, this,
                                  freqRange_.first, freqRange_.second, 1e-5);
  MU_DEBUG(fn, "%s reference bandwidth: %.2f microns, %.3e Hz",
           name(), width_*1e4*c_light/(center_*center_), width_);
}


double Detector::ccf(const Graybody &gray, double T, double Sigma) const
{ return gray.F(T, Sigma, this)/(width_*gray.Inu(center_, T, Sigma)); }

}  // namespace manticore