Radiation.hpp

Go to the documentation of this file.

#pragma once
 
#include <string>
#include <vector>
 
#include <gsl/gsl_spline.h>
 
namespace kariba {
 
struct CyclosynEmisParams {
    double nu;
    double b;
    gsl_spline* syn;
    gsl_interp_accel* acc_syn;
    gsl_spline* eldis;
    gsl_interp_accel* acc_eldis;
};
 
struct CyclosynAbsParams {
    double nu;
    double b;
    gsl_spline* syn;
    gsl_interp_accel* acc_syn;
    gsl_spline* derivs;
    gsl_interp_accel* acc_derivs;
};
 
struct ComintParams {
    double eph;
    double ephmin;
    double ephmax;
    gsl_spline* eldis;
    gsl_interp_accel* acc_eldis;
    gsl_spline* phodis;
    gsl_interp_accel* acc_phodis;
};
 
struct ComfncParams {
    double game;
    double e1;
    gsl_spline* phodis;
    gsl_interp_accel* acc_phodis;
};
 
struct DiskObsParams {
    double tin;
    double rin;
    double nu;
};
 
struct DiskIcParams {
    double gamma;
    double beta;
    double tin;
    double rin;
    double rout;
    double h;
    double z;
    double nu;
};
 
class Radiation {
  protected:
    std::vector<double> en_phot;         
    std::vector<double> num_phot;        
    std::vector<double> en_phot_obs;     
    std::vector<double> num_phot_obs;    
 
    double r, z;             
    double vol;              
    double beta;             
    double dopfac, angle;    
    double dopnum;           
    bool counterjet;         
    std::string geometry;    
 
  public:
    Radiation(size_t size);
 
    const std::vector<double>& get_energy() const { return en_phot; }
 
    const std::vector<double>& get_nphot() const { return num_phot; }
 
    const std::vector<double>& get_energy_obs() const { return en_phot_obs; }
 
    const std::vector<double>& get_nphot_obs() const { return num_phot_obs; }
 
    size_t get_size() const { return en_phot.size(); }
 
    double get_volume() const { return vol; }
 
    double integrated_luminosity(double numin, double numax);
 
    void set_beaming(double theta, double speed, double doppler);
    void set_inclination(double theta);
    void set_geometry(const std::string& geom, double l1, double l2);
    void set_geometry(const std::string& geom, double l1);
 
    void set_counterjet(bool flag);
    void test_arrays();
};
}    // namespace kariba