| [601] | 1 | // This may look like C code, but it is really -*- C++ -*- | 
|---|
|  | 2 | //-------------------------------------------------------------------------- | 
|---|
|  | 3 | // File and Version Information: | 
|---|
| [902] | 4 | //      $Id: nupower.h,v 1.7 2000-04-13 08:20:00 ansari Exp $ | 
|---|
| [601] | 5 | // | 
|---|
|  | 6 | // Description: | 
|---|
|  | 7 | // | 
|---|
|  | 8 | // History (add to end): | 
|---|
|  | 9 | //      Sophie   Oct, 1999  - creation | 
|---|
|  | 10 | // | 
|---|
|  | 11 | //------------------------------------------------------------------------ | 
|---|
|  | 12 |  | 
|---|
|  | 13 | #ifndef NUPOWER_H_SEEN | 
|---|
|  | 14 | #define NUPOWER_H_SEEN | 
|---|
|  | 15 |  | 
|---|
|  | 16 | #include "radspec.h" | 
|---|
|  | 17 | #include "convtools.h" | 
|---|
| [901] | 18 | namespace SOPHYA { | 
|---|
| [610] | 19 | // Power law spectra  f = a ((nu-nu0)/dnu)^b | 
|---|
| [601] | 20 | class PowerLawSpectra : public RadSpectra | 
|---|
|  | 21 | { | 
|---|
|  | 22 | public:   //Constructor | 
|---|
|  | 23 |  | 
|---|
| [668] | 24 | PowerLawSpectra(); | 
|---|
| [610] | 25 | PowerLawSpectra(double a, double b, double nu0, double dnu, double numin=0., double numax=9.e49); | 
|---|
| [601] | 26 |  | 
|---|
|  | 27 | virtual ~PowerLawSpectra(); | 
|---|
|  | 28 |  | 
|---|
|  | 29 | virtual double flux(double nu) const; | 
|---|
|  | 30 |  | 
|---|
|  | 31 | // Redefintion of the Print function | 
|---|
|  | 32 | virtual void      Print(ostream& os)  const; | 
|---|
| [902] | 33 |  | 
|---|
| [668] | 34 | virtual inline double giveNorm() const {return _a;}; | 
|---|
|  | 35 | virtual inline double giveNu0() const {return _nu0;}; | 
|---|
|  | 36 | virtual inline double giveDNu() const {return _dnu;}; | 
|---|
|  | 37 | virtual inline double giveExp() const {return _b;}; | 
|---|
|  | 38 |  | 
|---|
| [902] | 39 | virtual inline void setNorm(double a) {_a=a;}; | 
|---|
|  | 40 | virtual inline void setNu0(double nu0) {_nu0=nu0;}; | 
|---|
|  | 41 | virtual inline void setDNu(double dnu) {_dnu=dnu;}; | 
|---|
|  | 42 | virtual inline void setExp(double b) {_b=b;}; | 
|---|
|  | 43 |  | 
|---|
|  | 44 |  | 
|---|
| [601] | 45 |  | 
|---|
|  | 46 | protected: | 
|---|
| [610] | 47 | double _a, _b, _nu0, _dnu; | 
|---|
| [601] | 48 | }; | 
|---|
| [901] | 49 | } | 
|---|
| [601] | 50 | #endif | 
|---|