GyotoKappaDistributionSynchrotronSpectrum.h

Go to the documentation of this file.

 
/*
    Copyright 2018 Frederic Vincent, Thibaut Paumard
 
    This file is part of Gyoto.
 
    Gyoto is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
 
    Gyoto is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.
 
    You should have received a copy of the GNU General Public License
    along with Gyoto.  If not, see <http://www.gnu.org/licenses/>.
 */
 
#ifndef __GyotoKappaDistributionSynchrotronSpectrum_H_ 
#define __GyotoKappaDistributionSynchrotronSpectrum_H_ 
#include "GyotoSpectrum.h"
#include <GyotoBlackBodySpectrum.h>
 
namespace Gyoto {
  namespace Spectrum {
    class KappaDistributionSynchrotron;
  }
}
 
class Gyoto::Spectrum::KappaDistributionSynchrotron : public Gyoto::Spectrum::Generic {
  friend class Gyoto::SmartPointer<Gyoto::Spectrum::KappaDistributionSynchrotron>;
 protected:
  SmartPointer<Spectrum::BlackBody> spectrumBB_; 
  double numberdensityCGS_; 
  double angle_B_pem_; 
  double cyclotron_freq_; 
  double thetae_; 
  double kappaindex_; 
  double hypergeometric_; 
  bool angle_averaged_; 
  double gamma_max_; 
 
 public:
  GYOTO_OBJECT;
 
  KappaDistributionSynchrotron();
  KappaDistributionSynchrotron(const KappaDistributionSynchrotron &);
  virtual KappaDistributionSynchrotron * clone() const; 
 
  using Gyoto::Spectrum::Generic::operator();
  virtual double operator()(double nu) const;
#ifndef GYOTO_SWIGIMPORTED
  virtual double operator()(double nu,double ,double ds) const;
#endif
  // NB: the second argument, opacity in the Spectrum API
  // is useless here
 
  double numberdensityCGS() const;
  void numberdensityCGS(double rho);
  double angle_B_pem() const;
  void angle_B_pem(double aa);
  double cyclotron_freq() const;
  void cyclotron_freq(double ff);
  double thetae() const;
  void thetae(double th);
  double kappaindex() const;
  void kappaindex(double ind);
  double hypergeometric() const;
  void hypergeometric(double hh);
  bool angle_averaged() const;
  void angle_averaged(bool ang);
  double gamma_max() const;
  void gamma_max(double gammaMax);
  
  double jnuCGS(double nu) const;
  double jQnuCGS(double nu) const;
  double jUnuCGS(double nu) const;
  double jVnuCGS(double nu) const;
 
  double alphanuCGS(double nu) const;
  double alphaQnuCGS(double nu) const;
  double alphaUnuCGS(double nu) const;
  double alphaVnuCGS(double nu) const;
 
  double rQnuCGS(double nu) const;
  double rUnuCGS(double nu) const;
  double rVnuCGS(double nu) const;
 
  void radiativeQ(double jnu[], // output
                  double anu[], // output
                  double const nu_ems[],
                  size_t nbnu
                  ) ;
  void radiativeQ(double jInu[], double jQnu[], double jUnu[], double jVnu[],
        double aInu[], double aQnu[], double aUnu[], double aVnu[],
        double rQnu[], double rUnu[], double rVnu[],
        double const nu_ems[],
        size_t nbnu );
};
 
#endif