#ifndef STAR_H
#define STAR_H

#include <stdlib.h>
#include <string.h>
#include <cstdarg>
#include <string>
#include <math.h>

#include <CCfits>
#include "Geocen.h"

using namespace std ;
using namespace CCfits ;

/***************************************************************
* @author Richard J. Mathar
* @since 2006-02-06
*/
class Star {
public:
	/*
	* right ascension [rad]
	*/
	double ra ;
	/*
	* declination [rad]
	*/
	double dec ;

	/*
	* proper motion in right asc [arcsec/yr]
	*/
	double pmra ;

	/*
	* proper motion in declination asc [arcsec/yr]
	*/
	double pmdec ;

	/*
	* effective BB temperature [K]
	*/
	double T ;

	/*
	* K band magnitude
	*/
	double Kmag ;

	Star() ;

	Star(FITS &f, const int indx) ;
	Star(const double alpha, const double delta, const double pmAlpha=0., const double pmDelta=0.) ;

	vector<double> point(const Geocen & pos, const double ha) const ;

	vector<double> pointDt(const Geocen & pos, const double ha) const ;

#if 0
	double flux(double wn) const ;
#endif
	double flux(double wnstart, double wnstop) const;

protected:
private:
	/*
	* The scale factor to convert a black body flux to the Jansky scale:
	* The BB fluxes are to be multiplied with the scale factor to get
	* the flux at the earth.
	*/
	double bbscale ;
} ;

#endif /* STAR_H */