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