Detector.h

Go to the documentation of this file.

#ifndef DETECTOR_H
#define DETECTOR_H

#include <gsl/gsl_integration.h>
#include <gsl/gsl_spline.h>

#include <mutils/util.h>

namespace manticore {

// Forward reference.
class Graybody;


class Detector {
public:
  Detector(bool extend = true) noexcept : extend_(extend) { }

  ~Detector() {
    if (integ_)  { gsl_integration_workspace_free(integ_); }
    if (spline_) { gsl_spline_free(spline_); }
  }

  virtual std::string name() const = 0;

  virtual void init();

  virtual double freqBand() const noexcept { return center_; }

  virtual double freqWidth() const noexcept { return width_; }

  virtual std::pair<double, double> freqRange() const noexcept
  { return freqRange_; }

  bool isExtended() const noexcept { return extend_; }

  virtual double response(double nu) const;

  virtual double ccf(const Graybody &gray, double T, double Sigma = 1.0) const;

protected:
  static constexpr size_t subLimit = 512;

  static constexpr int intKey = GSL_INTEG_GAUSS41;

  using tableType = std::vector<std::pair<double, double>>;

  void calcWidth();

  std::pair<double,double> freqRange_;

  tableType resp_;

  bool  extend_;

  double center_ = 0.0,  
         width_  = 0.0;  

  gsl_integration_workspace *integ_ = nullptr;

  gsl_spline *spline_ = nullptr;
};

}  // namespace manticore

#endif  // DETECTOR_H