Detector.h

Go to the documentation of this file.

#ifndef DETECTOR_H
#define DETECTOR_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) { }

  virtual ~Detector() {
    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_; }

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

  bool isExtended() const noexcept { return extend_; }

  virtual double response(double nu, gsl_interp_accel *acc = nullptr) const;

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

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

  void calcWidth();

  std::pair<double,double> freqRange_;

  std::pair<double,double> sysError_;

  tableType resp_;

  bool  extend_;

  double center_ = 0.0,  
         width_  = 0.0;  

  gsl_spline *spline_ = nullptr;
};

}  // namespace manticore

#endif  // DETECTOR_H