Gyoto
GyotoEquatorialHotSpot.h
Go to the documentation of this file.
1 
7 /*
8  Copyright 2013, 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 __GyotoEquatorialHotSpot_h
27 #define __GyotoEquatorialHotSpot_h
28 
29 #include <GyotoAstrobj.h>
30 #include <GyotoThinDisk.h>
31 #include <GyotoWorldline.h>
33 #include <string>
34 
35 
36 namespace Gyoto {
37  namespace Astrobj {
38  class EquatorialHotSpot;
39  };
40 };
41 
44  public Gyoto::Worldline {
46 
47  private:
48  double sizespot_;
49  enum beaming_t {IsotropicBeaming=0, NormalBeaming=1,
50  RadialBeaming=2, IsotropicConstant=3};
51  beaming_t beaming_;
52  double beamangle_;
53  SmartPointer<Spectrum::ThermalSynchrotron> spectrumThermalSynch_; // Thermal distribution synchrotron spectrum
54  std::string magneticConfig_;
55 
56  public:
58  GYOTO_WORLDLINE;
63  virtual ~EquatorialHotSpot();
64  virtual EquatorialHotSpot * clone() const ;
65 
66  // Accessors for the Property list
67  void spotRadSize(double t);
68  double spotRadSize() const;
69 
70  void beaming(std::string const &b);
71  std::string beaming() const;
72 
73  void beamAngle(double t);
74  double beamAngle() const;
75 
76  void magneticConfiguration(std::string config);
77  std::string magneticConfiguration() const;
78 
79  //
80 
81  double getMass() const;
82  using Generic::metric;
84  void setInitialCondition(double coord[8]);
85 
86  void getVelocity(double const pos[4], double vel[4]);
87 
88  double emission(double nu_em, double dsem,
89  state_t const &,
90  double const coord_obj[8]) const;
91 
92  void radiativeQ(double *Inu, double *Qnu, double *Unu,
93  double *Vnu,
94  Eigen::Matrix4d *Onu,
95  double const *nuem , size_t nbnu,
96  double dsem,
97  state_t const &cph,
98  double const *co) const;
99 
100  void radiativeQ(double Inu[], // output
101  double Taunu[], // output
102  double const nu_ems[], size_t nbnu, // input
103  double dsem,
104  state_t const &coord_ph,
105  double const coord_obj[8]) const;
106 
107  // needed for legacy XML files
108  virtual int setParameter(std::string name,
109  std::string content,
110  std::string unit);
111 #ifdef GYOTO_USE_XERCES
112  // needed for wait_pos_
114  virtual void fillProperty(Gyoto::FactoryMessenger *fmp, Property const &p) const ;
115 #endif
116 
117 #endif
118 };
Astronomical objects (light emitters)
#define GYOTO_OBJECT
Declare class::properties and class::getProperties()
Definition: GyotoObject.h:84
Thermal synchrotron spectrum.
Geometrically thin disks and rings.
Timelike or null geodesics.
Definition: GyotoEquatorialHotSpot.h:44
virtual int setParameter(std::string name, std::string content, std::string unit)
Set parameter by name.
virtual EquatorialHotSpot * clone() const
Cloner.
double emission(double nu_em, double dsem, state_t const &, double const coord_obj[8]) const
Specific intensity Iν
virtual void fillProperty(Gyoto::FactoryMessenger *fmp, Property const &p) const
Output a single Property to XML.
double getMass() const
Get mass of particule.
void setParameters(FactoryMessenger *fmp)
Main loop in Subcontractor_t function.
void radiativeQ(double *Inu, double *Qnu, double *Unu, double *Vnu, Eigen::Matrix4d *Onu, double const *nuem, size_t nbnu, double dsem, state_t const &cph, double const *co) const
Compute the increment of Stokes parameters and transmission matrix. Polarised version of RadiaveQ.
void metric(SmartPointer< Metric::Generic > gg)
Set the Metric gg_.
void getVelocity(double const pos[4], double vel[4])
Get fluid 4-velocity at point.
std::string magneticConfig_
Specify the magnetic field configuration for polarisation.
Definition: GyotoEquatorialHotSpot.h:54
void radiativeQ(double Inu[], double Taunu[], double const nu_ems[], size_t nbnu, double dsem, state_t const &coord_ph, double const coord_obj[8]) const
emission and transmission together
virtual double deltaMax(double coord[8])
Get max step constraint for adaptive integration.
virtual SmartPointer< Metric::Generic > metric() const
Get the Metric gg_.
Geometrically thin disks and rings.
Definition: GyotoThinDisk.h:71
Factory / SmartPointee::Subcontractor_t interface.
Definition: GyotoFactoryMessenger.h:92
Property that can be set and got using standard methods.
Definition: GyotoProperty.h:608
Pointers performing reference counting.
Definition: GyotoSmartPointer.h:135
Timelike or null geodesics.
Definition: GyotoWorldline.h:240
double deltaMax() const
Get delta_max_.
Namespace for the Gyoto library.
Definition: GyotoAstrobj.h:44