Compton.hpp

Go to the documentation of this file.

#pragma once
 
#include <gsl/gsl_spline2d.h>
 
#include "Radiation.hpp"
 
namespace kariba {
 
class Compton : public Radiation {
  protected:
    size_t Niter;          
    double tau, ypar;      
    double rphot;          
    double escape_corr;    
 
    std::vector<double> seed_energ;    
    std::vector<double> seed_urad;     
    std::vector<double>
        iter_urad;    
 
    gsl_spline* seed_ph;           
    gsl_interp_accel* acc_seed;    
 
    gsl_spline* iter_ph;           
    gsl_interp_accel* acc_iter;    
 
    gsl_spline2d* esc_p_sph;    
    gsl_spline2d* esc_p_cyl;    
    gsl_interp_accel* acc_tau;    
    gsl_interp_accel* acc_Te;     
 
  public:
    ~Compton();
    Compton(size_t size, size_t seed_size);
 
    virtual double comintegral(size_t it, double blim, double ulim, double nu, double numin,
                               double numax, gsl_spline* eldis, gsl_interp_accel* acc_eldis);
    virtual void compton_spectrum(double gmin, double gmax, gsl_spline* eldis,
                                  gsl_interp_accel* acc_eldis);
 
    virtual void cyclosyn_seed(const std::vector<double>& seed_arr,
                               const std::vector<double>& seed_lum);
    virtual void bb_seed_k(const std::vector<double>& seed_arr, double Urad, double Tbb);
    virtual void bb_seed_kev(const std::vector<double>& seed_energ, double Urad, double Tbb);
    virtual void shsdisk_seed(const std::vector<double>& seed_arr, double tin, double rin,
                              double rout, double h, double z);
 
    virtual void set_frequency(double numin, double numax);
    virtual void set_tau(double n, double gam);
    virtual void set_tau(double _tau);
    virtual void set_escape(double escape);
    virtual void set_niter(double nu0, double Te);
    virtual void set_niter(size_t n);
    virtual void seed_freq_array(const std::vector<double>& seed_energ);
 
    virtual double get_tau() const { return tau; };
 
    virtual double get_ypar() const { return ypar; };
 
    virtual void reset();
    virtual void urad_test();
    virtual void test();
 
    friend double comfnc(double ein, void* p);
    friend double comint(double gam, void* p);
    friend double disk_integral(double alfa, void* p);
};
 
}    // namespace kariba