| 1 | #ifndef MARS_MAstro
|
|---|
| 2 | #define MARS_MAstro
|
|---|
| 3 |
|
|---|
| 4 | #ifndef ROOT_TROOT
|
|---|
| 5 | #include <TROOT.h>
|
|---|
| 6 | #endif
|
|---|
| 7 | #ifndef ROOT_TVector2
|
|---|
| 8 | #include <TVector2.h>
|
|---|
| 9 | #endif
|
|---|
| 10 |
|
|---|
| 11 | class TVector3;
|
|---|
| 12 |
|
|---|
| 13 | class MAstro
|
|---|
| 14 | {
|
|---|
| 15 | private:
|
|---|
| 16 | static Double_t Round(Double_t val);
|
|---|
| 17 | static Double_t Trunc(Double_t val);
|
|---|
| 18 |
|
|---|
| 19 | static Double_t Kepler(Double_t m, Double_t ecc);
|
|---|
| 20 |
|
|---|
| 21 | public:
|
|---|
| 22 | // Angle/Time conversion functions
|
|---|
| 23 | static Double_t Hms2Sec(Int_t deg, UInt_t min, Double_t sec, char sgn='+');
|
|---|
| 24 | static Double_t Dms2Rad(Int_t deg, UInt_t min, Double_t sec, Char_t sgn='+');
|
|---|
| 25 | static Double_t Hms2Rad(Int_t hor, UInt_t min, Double_t sec, Char_t sgn='+');
|
|---|
| 26 | static Double_t Dms2Deg(Int_t deg, UInt_t min, Double_t sec, Char_t sgn='+');
|
|---|
| 27 | static Double_t Hms2Deg(Int_t hor, UInt_t min, Double_t sec, Char_t sgn='+');
|
|---|
| 28 | static Double_t Dms2Hor(Int_t deg, UInt_t min, Double_t sec, Char_t sgn='+');
|
|---|
| 29 | static Double_t Hms2Hor(Int_t hor, UInt_t min, Double_t sec, Char_t sgn='+');
|
|---|
| 30 |
|
|---|
| 31 | static void Day2Hms(Double_t rad, Char_t &sgn, UShort_t °, UShort_t &min, UShort_t &sec);
|
|---|
| 32 | static void Rad2Dms(Double_t rad, Char_t &sgn, UShort_t °, UShort_t &min, UShort_t &sec);
|
|---|
| 33 | static void Rad2Hms(Double_t rad, Char_t &sgn, UShort_t °, UShort_t &min, UShort_t &sec);
|
|---|
| 34 | static void Deg2Dms(Double_t rad, Char_t &sgn, UShort_t °, UShort_t &min, UShort_t &sec);
|
|---|
| 35 | static void Deg2Hms(Double_t rad, Char_t &sgn, UShort_t °, UShort_t &min, UShort_t &sec);
|
|---|
| 36 | static void Hor2Dms(Double_t rad, Char_t &sgn, UShort_t °, UShort_t &min, UShort_t &sec);
|
|---|
| 37 | static void Hor2Hms(Double_t rad, Char_t &sgn, UShort_t °, UShort_t &min, UShort_t &sec);
|
|---|
| 38 |
|
|---|
| 39 | static void Day2Hm(Double_t rad, Char_t &sgn, UShort_t °, Double_t &min);
|
|---|
| 40 | static void Rad2Dm(Double_t rad, Char_t &sgn, UShort_t °, Double_t &min);
|
|---|
| 41 | static void Rad2Hm(Double_t rad, Char_t &sgn, UShort_t °, Double_t &min);
|
|---|
| 42 | static void Deg2Dm(Double_t rad, Char_t &sgn, UShort_t °, Double_t &min);
|
|---|
| 43 | static void Deg2Hm(Double_t rad, Char_t &sgn, UShort_t °, Double_t &min);
|
|---|
| 44 | static void Hor2Dm(Double_t rad, Char_t &sgn, UShort_t °, Double_t &min);
|
|---|
| 45 | static void Hor2Hm(Double_t rad, Char_t &sgn, UShort_t °, Double_t &min);
|
|---|
| 46 |
|
|---|
| 47 | // Angle treatment functions
|
|---|
| 48 | static Bool_t String2Angle(TString &str, Double_t &ret);
|
|---|
| 49 | static Bool_t Coordinate2Angle(const TString &str, Double_t &ret);
|
|---|
| 50 |
|
|---|
| 51 | static Double_t AngularDistance(Double_t theta0, Double_t phi0, Double_t theta1, Double_t phi1);
|
|---|
| 52 |
|
|---|
| 53 | // Time treatment functions
|
|---|
| 54 | static void Mjd2Ymd(UInt_t mjd, UShort_t &y, Byte_t &m, Byte_t &d);
|
|---|
| 55 | static Int_t Ymd2Mjd(UShort_t y, Byte_t m, Byte_t d);
|
|---|
| 56 |
|
|---|
| 57 | static Double_t UT2GMST(Double_t ut1);
|
|---|
| 58 |
|
|---|
| 59 | // Rotation angle between local and sky coordinate system
|
|---|
| 60 | static void RotationAngle(Double_t sinl, Double_t cosl, Double_t theta, Double_t phi, Double_t &sin, Double_t &cos);
|
|---|
| 61 | static Double_t RotationAngle(Double_t sinl, Double_t cosl, Double_t theta, Double_t phi);
|
|---|
| 62 |
|
|---|
| 63 | static Double_t GetMoonPhase(Double_t mjd);
|
|---|
| 64 | static Double_t GetMoonPeriod(Double_t mjd);
|
|---|
| 65 | static Int_t GetMagicPeriod(Double_t mjd);
|
|---|
| 66 |
|
|---|
| 67 | static TVector2 GetDistOnPlain(const TVector3 &v0, TVector3 v1, Double_t dist=1);
|
|---|
| 68 |
|
|---|
| 69 | ClassDef(MAstro, 0)
|
|---|
| 70 | };
|
|---|
| 71 |
|
|---|
| 72 | #endif
|
|---|
