Gyoto
GyotoKappaDistributionSynchrotronSpectrum.h
Go to the documentation of this file.
1 
7 /*
8  Copyright 2018 Frederic Vincent, Thibaut Paumard
9 
10  This file is part of Gyoto.
11 
12  Gyoto is free software: you can redistribute it and/or modify
13  it under the terms of the GNU General Public License as published by
14  the Free Software Foundation, either version 3 of the License, or
15  (at your option) any later version.
16 
17  Gyoto is distributed in the hope that it will be useful,
18  but WITHOUT ANY WARRANTY; without even the implied warranty of
19  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
20  GNU General Public License for more details.
21 
22  You should have received a copy of the GNU General Public License
23  along with Gyoto. If not, see <http://www.gnu.org/licenses/>.
24  */
25 
26 #ifndef __GyotoKappaDistributionSynchrotronSpectrum_H_
27 #define __GyotoKappaDistributionSynchrotronSpectrum_H_
28 #include "GyotoSpectrum.h"
29 #include <GyotoBlackBodySpectrum.h>
30 
31 namespace Gyoto {
32  namespace Spectrum {
33  class KappaDistributionSynchrotron;
34  }
35 }
36 
51  protected:
52  SmartPointer<Spectrum::BlackBody> spectrumBB_;
54  double angle_B_pem_;
55  double cyclotron_freq_;
56  double thetae_;
57  double kappaindex_;
58  double hypergeometric_;
60  double gamma_max_;
61 
62  public:
64 
66  KappaDistributionSynchrotron(const KappaDistributionSynchrotron &);
67  virtual KappaDistributionSynchrotron * clone() const;
68 
69  using Gyoto::Spectrum::Generic::operator();
76  virtual double operator()(double nu) const;
85 #ifndef GYOTO_SWIGIMPORTED
86  virtual double operator()(double nu,double ,double ds) const;
87 #endif
88  // NB: the second argument, opacity in the Spectrum API
89  // is useless here
90 
91  double numberdensityCGS() const;
92  void numberdensityCGS(double rho);
93  double angle_B_pem() const;
94  void angle_B_pem(double aa);
95  double cyclotron_freq() const;
96  void cyclotron_freq(double ff);
97  double thetae() const;
98  void thetae(double th);
99  double kappaindex() const;
100  void kappaindex(double ind);
101  double hypergeometric() const;
102  void hypergeometric(double hh);
103  bool angle_averaged() const;
104  void angle_averaged(bool ang);
105  double gamma_max() const;
106  void gamma_max(double gammaMax);
107 
114  double jnuCGS(double nu) const;
120  double alphanuCGS(double nu) const;
121 
126  void radiativeQ(double jnu[], // output
127  double anu[], // output
128  double const nu_ems[],
129  size_t nbnu
130  ) ;
131 };
132 
133 #endif
#define GYOTO_OBJECT
Declare class::properties and class::getProperties()
Definition: GyotoObject.h:84
double numberdensityCGS_
Number density in CGS UNITS (careful)
Definition: GyotoKappaDistributionSynchrotronSpectrum.h:53
double hypergeometric_
Hypergeometric function evaluation.
Definition: GyotoKappaDistributionSynchrotronSpectrum.h:58
Kappa-distribution synchrotron spectrum from Pandya et al. (2016)
Definition: GyotoKappaDistributionSynchrotronSpectrum.h:49
bool angle_averaged_
Boolean for angle averaging.
Definition: GyotoKappaDistributionSynchrotronSpectrum.h:59
double cyclotron_freq_
Cyclotron frequency (e*B / 2*pi*me*c)
Definition: GyotoKappaDistributionSynchrotronSpectrum.h:55
double gamma_max_
Maximum gamma value to integer the distribution.
Definition: GyotoKappaDistributionSynchrotronSpectrum.h:60
virtual KappaDistributionSynchrotron * clone() const
Cloner.
SmartPointer< Spectrum::BlackBody > spectrumBB_
blackbody emission
Definition: GyotoKappaDistributionSynchrotronSpectrum.h:52
void radiativeQ(double jnu[], double anu[], double const nu_ems[], size_t nbnu)
Spectrum of a simple object (e.g. Star)
Namespace for the Gyoto library.
Definition: GyotoAstrobj.h:43
Pointers performing reference counting.
Definition: GyotoProperty.h:45
double angle_B_pem_
Angle between Bfield and emission direction (rad)
Definition: GyotoKappaDistributionSynchrotronSpectrum.h:54
double kappaindex_
Kappa distribution index.
Definition: GyotoKappaDistributionSynchrotronSpectrum.h:57
Spectrum emitted by an Astrobj.
Definition: GyotoSpectrum.h:134
double thetae_
Dimensionless electron temperature.
Definition: GyotoKappaDistributionSynchrotronSpectrum.h:56
I_nu(nu, T) = cst_*2*h*nu^3/c^2/(exp(h*nu/k*T)-1.);.