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:
    const double log_floor = -230; // about 1e-100 
    size_t Niter;          
    double tau, ypar;      
    double rphot;          
    double escape_corr;    
 
    std::vector<double> target_energy;    
    std::vector<double> log_target_diff_spec;     
    std::vector<double>
        log_target_diff_spec_iter;    
 
 
    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;     
 
    gsl_integration_workspace* w1;
    gsl_integration_workspace* w2;
 
  public:
    ~Compton();
    Compton(size_t size, size_t target_size);
 
    friend double comfnc(double logein, void* pars);
    friend double comint(double gam, void* pars);
    friend double disk_integral(double alfa, void* p);
 
    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 set_target_energy_array(const std::vector<double>& new_target_energy);
    virtual void set_target_frequency_array(const std::vector<double>& new_target_frequency);
 
    virtual void add_target_diff_spec(const std::vector<double>& new_target_diff_spec);
    virtual void add_target_energy_density(
        const std::vector<double>& new_target_energy,
        const std::vector<double>& new_target_energy_density);
    virtual void add_target_number_density(
        const std::vector<double>& new_target_energy,
        const std::vector<double>& new_target_number_density);
    virtual void add_target_number_density_on_freq_grid(
        const std::vector<double>& new_target_frequency,
        const std::vector<double>& new_target_number_density);
 
    virtual void cyclosyn_seed(const std::vector<double>& syn_energies,
                               const std::vector<double>& syn_number_rates);
    virtual void bb_seed_k(double Urad, double Tbb);
    virtual void bb_seed_kev(double Urad, double Tbb);
    virtual void shsdisk_seed(double tin, double rin, double rout, double h, double z);
 
    virtual void set_niter(double nu0, double Te);
    virtual void set_niter(size_t n);
    virtual void set_tau(double n, double gam);
    virtual void set_tau(double _tau);
    virtual void set_frequency(double numin, double numax);
    virtual void set_escape(double escape);
 
    virtual std::vector<double> get_target_energy() const;
    virtual std::vector<double> get_target_diff_spec() const;
 
    virtual double get_tau() const { return tau; };
    virtual double get_ypar() const { return ypar; };
 
    virtual void reset();
    virtual void urad_test();
    virtual void test();
 
    // -------------------------------------------------------------------------
    // Legacy API — backwards compatibility wrappers
    // The seed_arr parameter sets the energy grid, but it is cleaner to set it
    // explicitly with set_target_energy_array() / set_target_frequency_array().
    // -------------------------------------------------------------------------
 
 
    virtual void bb_seed_k(const std::vector<double>& seed_arr, double Urad, double Tbb) {
        set_target_energy_array(seed_arr);
        bb_seed_k(Urad, Tbb);
    }
 
    virtual void bb_seed_kev(const std::vector<double>& seed_arr, double Urad, double Tbb) {
        set_target_energy_array(seed_arr);
        bb_seed_kev(Urad, Tbb);
    }
 
    virtual void shsdisk_seed(const std::vector<double>& seed_arr, double tin, double rin,
                              double rout, double h, double z) {
        set_target_energy_array(seed_arr);
        shsdisk_seed(tin, rin, rout, h, z);
    }
};
 
}    // namespace kariba