GyotoThermalBremsstrahlungSpectrum.h

Go to the documentation of this file.

/*
    Copyright 2014 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 __GyotoThermalBremsstrahlungSpectrum_H_ 
#define __GyotoThermalBremsstrahlungSpectrum_H_ 
#include "GyotoSpectrum.h"
#include <GyotoBlackBodySpectrum.h>

namespace Gyoto {
  namespace Spectrum {
    class ThermalBremsstrahlung;
  }
}

class Gyoto::Spectrum::ThermalBremsstrahlung : public Gyoto::Spectrum::Generic {
  friend class Gyoto::SmartPointer<Gyoto::Spectrum::ThermalBremsstrahlung>;
 protected:
  SmartPointer<Spectrum::BlackBody> spectrumBB_; 
  double T_; 
  double Tm1_; 
  double Tm05_; 
  double numberdensityCGS_; 

 public:
  GYOTO_OBJECT;

  ThermalBremsstrahlung();
  ThermalBremsstrahlung(const ThermalBremsstrahlung &);
  virtual ThermalBremsstrahlung * clone() const; 

  using Gyoto::Spectrum::Generic::operator();
  virtual double operator()(double nu) const;
#ifndef GYOTO_SWIGIMPORTED
  virtual double operator()(double nu,double ,double ds) const;
#endif
  // NB: the second argument, opacity in the Spectrum API
  // is useless here

  double temperature() const;
  void temperature(double tt);
  double numberdensityCGS() const;
  void numberdensityCGS(double rho);

  double jnuCGS(double nu) const;
  double alphanuCGS(double nu) const;

  void radiativeQ(double jnu[], // output
                  double alphanu[], // output
                  double const nu_ems[],
                  size_t nbnu
                  );

};

#endif