GyotoPageThorneDisk.h

Go to the documentation of this file.

 
/*
  Copyright 2011-2016, 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 __GyotoPageThorneDisk_H_ 
#define __GyotoPageThorneDisk_H_ 
 
#include <iostream>
#include <fstream>
#include <iomanip>
 
namespace Gyoto{
  namespace Astrobj { class PageThorneDisk; }
}
 
//#include <GyotoMetric.h>
#include <GyotoThinDisk.h>
#include <GyotoBlackBodySpectrum.h>
 
class Gyoto::Astrobj::PageThorneDisk
: public Astrobj::ThinDisk,
  public Hook::Listener
{
  friend class Gyoto::SmartPointer<Gyoto::Astrobj::PageThorneDisk>;
 private:
  double aa_; 
  double aa2_; 
  double x0_; 
  double x1_; 
  double x2_; 
  double x3_; 
  double mdot_; 
  bool uniflux_; 
  SmartPointer<Spectrum::BlackBody> spectrumBB_; 
  
  // Constructors - Destructor
  // -------------------------
 public:
  GYOTO_OBJECT;
  GYOTO_OBJECT_THREAD_SAFETY;
 
  PageThorneDisk(); 
  
  PageThorneDisk(const PageThorneDisk& ) ;
  virtual PageThorneDisk* clone () const; 
  
  virtual ~PageThorneDisk() ;                        
  
  // Accessors
  // ---------
 public:
  using ThinDisk::metric;
  virtual void metric(SmartPointer<Metric::Generic>);
 
  void mdot(double v);
  double mdot() const;
  void uniFlux(bool t) ;
  bool uniFlux() const ;
 
 private:
  virtual void updateSpin() ;
 
 public:
  using ThinDisk::emission;
  virtual double emission(double nu_em, double dsem,
                          state_t const &c_ph, double const c_obj[8]=NULL) const;
 
  virtual double bolometricEmission(double nuem, double dsem,
                                    double const c_obj[8]) const;
 
  virtual void processHitQuantities(Photon* ph, state_t const &coord_ph_hit,
                                   double const *coord_obj_hit, double dt,
                                   Astrobj::Properties* data) const;
 
  Gyoto::Quantity_t getDefaultQuantities();
 
  // Hook::Listener API //
 public:
  virtual void tell(Gyoto::Hook::Teller *msg);
 
  virtual 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;
 
};
 
#endif