Powerlaw.hpp

Go to the documentation of this file.

#pragma once
 
#include <iostream>
 
#include <gsl/gsl_spline.h>
 
#include "Particles.hpp"
 
namespace kariba {
 
class Powerlaw : public Particles {
  protected:
    double pspec, plnorm;
    double pmin, pcut; // was pmax, has been adjusted in the code 
    bool isEfficient;    
 
  public:
    Powerlaw(size_t size);
 
    virtual void set_p(double min, double ucom, double bfield, double betaeff, double r,
                       double fsc);
    virtual void set_p(double min, double gmax);
    virtual void set_ndens();
    virtual void set_pspec(double s1);
    virtual void set_norm(double n);
 
    virtual bool get_Efficiency() const { return isEfficient; }
 
    virtual void cooling_steadystate(double ucom, double n0, double bfield, double r,
                                     double tshift);
    virtual double max_p(double ucom, double bfield, double betaeff, double r, double fsc);
 
    // // some functions for protons
    // virtual void set_energy(double gpmin, double fsc, double f_beta, double bfield, double r_g,
    //                         double z, double r, int infosw, double protdens, double nwind,
    //                         double Uradjet, const std::string& outputConfiguration,
    //                         const std::string& source);
    // virtual void check_secondary_charged_syn(double bfield, double gpmax);
    // virtual void ProtonTimescales(double& logdgp, double fsc, double f_beta, double bfield,
    //                               double gpmin, double& gpmax, double r_g, double z, double r,
    //                               int infosw, double nwind, double Uradjet,
    //                               const std::string& outputConfiguration,
    //                               const std::string& source);
    // virtual double sigma_pp(double Ep);
    // virtual double set_normprot(double nprot);
    // virtual void set_gdens(double& plfrac_p, double Up, double protdens);
    // virtual void set_gdens(double r, double protdens, double nwind, double bfield, double plfrac,
    //                        double Urad);
    // virtual void set_gdens_pdens(double r, double beta, double Ljet, double ep, double pspec,
    //                              double& protdens);
 
    // // secondary electrons from pp
    // virtual void set_pp_elecs(gsl_interp_accel* acc_Jp, gsl_spline* spline_Jp, double ntot_prot,
    //                           double nwind, double plfrac, double gammap_min, double gammap_max,
    //                           double bfield, double r);
    // // convert secondary electrons from pg from Neutrinos units proper units for
    // // synchrotron radiation
    // virtual void set_pg_electrons(const std::vector<double>& energy,
    //                               const std::vector<double>& density, double f_beta, double r,
    //                               double vol, double B);
 
    // // function for electrons from γγ annihilation
    // virtual void Qggeefunction(double r, double vol, double bfield, size_t phot_number,
    //                            const std::vector<double>& en_perseg,
    //                            const std::vector<double>& lum_perseg, double gmax);
 
    virtual void test();
};
 
}    // namespace kariba