GyotoOscilTorus.h

Go to the documentation of this file.

/*
 *   Copyright (c) 2013, 2018 Frederic Vincent & Thibaut Paumard
 *
 *
 */

#ifndef __OscilTorus_h
#define __OscilTorus_h

#include <GyotoAstrobj.h>
#include <GyotoStandardAstrobj.h>
#include <GyotoKerrBL.h>


namespace Gyoto {
  namespace Astrobj {
    class OscilTorus;
  };
  class FactoryMessenger;
  namespace Spectrum {
    class Generic;
  }
};

class Gyoto::Astrobj::OscilTorus
: public Gyoto::Astrobj::Standard,
  public Hook::Listener {
  friend class Gyoto::SmartPointer<Gyoto::Astrobj::OscilTorus>;

 private:
  // Members corresponding to properties:  
  double c_;
  unsigned long mode_;
  double polycst_; 
  double polyindex_; 
  double central_density_; 
  enum perturb_t {Radial=1, Vertical=2, X=3, Plus=4, Breathing=5};
  perturb_t perturb_kind_;
  std::string emitting_area_; 
  double perturb_intens_; 

  // Cached values:
  SmartPointer<Gyoto::Metric::KerrBL> kerrbl_;
  std::vector<double> tt_;
  std::vector<double> area_; // tt_ and area_ contain area of cross section at time tt
  size_t nbt_; 
  int with_cross_; 

  double sigma_; 
  double alpha_; 
  double w1_; 
  double w2_; 
  double omr2_; 
  double omth2_;
  double Omegac_; 
  double lc_;
  double g_rr_; 
  double g_thth_;
  bool hold_;
  
 public:
  GYOTO_OBJECT;
  OscilTorus();
  OscilTorus(const OscilTorus &o);
  virtual ~OscilTorus();
  virtual OscilTorus * clone() const ;

  GYOTO_OBJECT_ACCESSORS_UNIT(largeRadius);
  GYOTO_OBJECT_ACCESSORS(unsigned long, mode);
  GYOTO_OBJECT_ACCESSORS(double, polyCst);
  GYOTO_OBJECT_ACCESSORS(double, polyIndex);
  GYOTO_OBJECT_ACCESSORS(double, centralDensity);
  GYOTO_OBJECT_ACCESSORS_STRING(perturbKind);
  GYOTO_OBJECT_ACCESSORS(double, perturbIntens);
  GYOTO_OBJECT_ACCESSORS_STRING(emittingArea);
  using Generic::metric;
  virtual void metric(Gyoto::SmartPointer<Gyoto::Metric::Generic>);

  virtual double operator()(double const coord[4]) ;
  virtual void getVelocity(double const pos[4], double vel[4]) ;
  double emission(double nuem,double,state_t const &,double const *) const;

#ifdef GYOTO_USE_XERCES
  virtual void setParameters(Gyoto::FactoryMessenger *fmp) ;
#endif

  virtual void updateCachedValues();
  void computeXbYb(const double * pos, double & xb, double & yb) ;

  // Hook::Listener API //
 public:
  virtual void tell(Gyoto::Hook::Teller *msg);

};

#endif