source.h

Go to the documentation of this file.

#ifndef SOURCE_H
#define SOURCE_H

#include <vector>
#include <map>
class Source
{
public:
    /* Constructor */
    Source();

    /* Destructor - it's important that this is virtual ! */
    virtual ~Source();

    virtual std::vector<double> get_spectral_density(std::vector<double> wavelength);

    virtual double get_spectral_density(double wavelength);

    virtual std::vector<double> get_spectrum(std::vector<double> wavelength);

    void set_doppler_shift(double shift);

    void set_integration_steps(int n);

private:
    double integral_s(double a, double b, int n);

    double shift; 
    int integration_steps; 
};

class Constant : public Source{
public:
    /* Constructor */
    Constant();
    Constant(double value);
    double get_spectral_density(double wavelength);
private:
    /* Constant spectral density value */
    double value;
};

class IdealEtalon : public Source{
public:
    IdealEtalon(double d, double n, double theta, double R);
    static double coefficient_of_finesse(double R);

    double FSR();

    double F();

    static double T(double wl, double theta, double d, double n, double cF);
    double get_spectral_density(double wavelength);
private:
    double d;
    double n;
    double theta;
    double R;
    double cF;
};


class Blackbody : public Source{
public:
    Blackbody(double T);
    static double planck(const double& T, const double& wavelength);
    double get_spectral_density(double wavelength);
private:
    double T; 
};

class PhoenixSpectrum : public Source{
public:
    PhoenixSpectrum(std::string spectrum_file, std::string wavelength_file, const double& min_wavelength, const double& max_wavelength);
    void read_spectrum(std::string spectrum_file, std::string wavelength_file, const double& min_wavelength, const double& max_wavelength);
    double get_spectral_density(double wavelength);
private:
    std::map<double, double> data;

};

class LineList : public Source{
public:
    LineList(std::string linelist);
    void read_spectrum(std::string linelist);
    double get_spectral_density(double wavelength);
    std::vector<double> get_spectrum(std::vector<double> wavelength);
    std::vector<double> get_wavelength();
private:
    std::map<double, double> data;
};
#endif // SOURCE_H