GyotoFlaredDiskSynchrotron.h

Go to the documentation of this file.

/*
    Copyright 2019-2021 Frederic Vincent, Thibaut Paumard, Nicolas Aimar

    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 __GyotoFlaredDiskSynchrotron_H_ 
#define __GyotoFlaredDiskSynchrotron_H_ 

#include <iostream>
#include <fstream>
#include <iomanip>

namespace Gyoto{
  namespace Astrobj { class FlaredDiskSynchrotron; }
  class GridData2D;
}

#include <GyotoStandardAstrobj.h>
#include <GyotoGridData2D.h>
#include <GyotoKappaDistributionSynchrotronSpectrum.h>

class Gyoto::Astrobj::FlaredDiskSynchrotron
: public Astrobj::Standard,
  public GridData2D,
  public Hook::Listener
{
  friend class Gyoto::SmartPointer<Gyoto::Astrobj::FlaredDiskSynchrotron>;
 private:
  SmartPointer<Spectrum::KappaDistributionSynchrotron> spectrumKappaSynch_;
  std::string filename_; 
  double hoverR_; 
  double numberDensityMax_cgs_; 
  double temperatureMax_; 
  double BMax_cgs_; 
  double beta_;
  double * density_; 

  double * velocity_; 
  double * Bvector_; 
  double * time_array_; 
  double magnetizationParameter_; 
  double deltat_;
  double gamm1_; 
  bool flag_; 

 public:
  GYOTO_OBJECT;
  GYOTO_OBJECT_THREAD_SAFETY;
  
  // Constructors - Destructor
  // -------------------------
  FlaredDiskSynchrotron(); 
  
  FlaredDiskSynchrotron(const FlaredDiskSynchrotron& ) ;
  virtual FlaredDiskSynchrotron* clone () const; 
  
  virtual ~FlaredDiskSynchrotron() ;                        
  
  // Accessors
  // ---------
 public:
  void file(std::string const &f) ;
  std::string file() const;
  void hoverR(double const hor) ;
  double hoverR() const;
  /*
    timeTranslation shifts the value of GridData2D::tmin_ and tmax_,
    allowing to scan the full simulation without having to change
    the value of the Screen observation time (which is typically
    not provided in M unit in the XML). 
    Choosing a negative timeTranslation, i.e. performing tmin_,tmax_-=dt, 
    amounts to increasing the Screen observation time by the same value, 
    tobs+=dt.

   */
  void timeTranslation_inMunit(double const dt) ;
  double timeTranslation_inMunit() const ;
  void magnetizationParameter(double rr);
  double magnetizationParameter() const;
  void kappaIndex(double index);
  double kappaIndex()const;
  double numberDensityMax() const;
  double numberDensityMax(std::string const &unit) const;
  void numberDensityMax(double dens) ;
  void numberDensityMax(double dens, std::string const &unit);
  void temperatureMax(double tt);
  double temperatureMax() const;
  void polytropicIndex(double gamma);
  double polytropicIndex() const;
  void betaAtMax(double beta);
  double betaAtMax() const;
  void copyDensity(double const *const density,
                   size_t const naxes[3]);
  double const * getDensity() const;
  void copyVelocity(double const *const velocity,
                    size_t const naxes[3]);
  double const * getVelocity() const;
  void copyBvector(double const *const Bvector,
        size_t const naxes[3]);
  double const * getBvector() const;
  void copyTimeArray(double const *const time_array, size_t const ntimes);
  double const * getTimeArray() const;
 public:
  using Generic::metric;
  std::vector<size_t> fitsRead(std::string filename) ;
  virtual double operator()(double const coord[4]) ;
  virtual 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;
  virtual void getVelocity(double const pos[4], double vel[4]) ;



};

#endif