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 day, Char_t &sgn, UShort_t &hor, 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 &hor, UShort_t &min, UShort_t &sec);
|
---|
34 | static void Deg2Dms(Double_t deg, Char_t &sgn, UShort_t °, UShort_t &min, UShort_t &sec);
|
---|
35 | static void Deg2Hms(Double_t deg, Char_t &sgn, UShort_t &hor, UShort_t &min, UShort_t &sec);
|
---|
36 | static void Hor2Dms(Double_t hor, Char_t &sgn, UShort_t °, UShort_t &min, UShort_t &sec);
|
---|
37 | static void Hor2Hms(Double_t hor, Char_t &sgn, UShort_t &hor, UShort_t &min, UShort_t &sec);
|
---|
38 |
|
---|
39 | static void Day2Hm(Double_t day, Char_t &sgn, UShort_t &hor, 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 &hor, Double_t &min);
|
---|
42 | static void Deg2Dm(Double_t deg, Char_t &sgn, UShort_t °, Double_t &min);
|
---|
43 | static void Deg2Hm(Double_t deg, Char_t &sgn, UShort_t &hor, 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 &hor, Double_t &min);
|
---|
46 |
|
---|
47 | // Print funtions
|
---|
48 | static TString GetStringDeg(Double_t deg, const char *fmt="%c%d:%02d:%02d");
|
---|
49 | static TString GetStringHor(Double_t deg, const char *fmt="%c%d:%02d:%02d");
|
---|
50 |
|
---|
51 | // Angle treatment functions
|
---|
52 | static Bool_t String2Angle(TString &str, Double_t &ret);
|
---|
53 | static Bool_t Coordinate2Angle(const TString &str, Double_t &ret);
|
---|
54 | static TString Angle2Coordinate(Double_t val);
|
---|
55 |
|
---|
56 | static Double_t AngularDistance(Double_t theta0, Double_t phi0, Double_t theta1, Double_t phi1);
|
---|
57 |
|
---|
58 | // Time treatment functions
|
---|
59 | static void Mjd2Ymd(UInt_t mjd, UShort_t &y, Byte_t &m, Byte_t &d);
|
---|
60 | static Int_t Ymd2Mjd(UShort_t y, Byte_t m, Byte_t d);
|
---|
61 |
|
---|
62 | static Double_t UT2GMST(Double_t ut1);
|
---|
63 |
|
---|
64 | // Rotation angle between local and sky coordinate system
|
---|
65 | static void RotationAngle(Double_t sinl, Double_t cosl, Double_t theta, Double_t phi, Double_t &sin, Double_t &cos);
|
---|
66 | static Double_t RotationAngle(Double_t sinl, Double_t cosl, Double_t theta, Double_t phi);
|
---|
67 |
|
---|
68 | static Double_t GetMoonPhase(Double_t mjd);
|
---|
69 | static Double_t GetMoonPeriod(Double_t mjd);
|
---|
70 | static Int_t GetMagicPeriod(Double_t mjd);
|
---|
71 |
|
---|
72 | // Get distance between v1 and v0 in a plain perpendicular to v0 in distance dist
|
---|
73 | static TVector2 GetDistOnPlain(const TVector3 &v0, TVector3 v1, Double_t dist=1);
|
---|
74 |
|
---|
75 | // Calculate absolute misspointing from error in zd and error in az
|
---|
76 | static Double_t GetDevAbs(Double_t nomzd, Double_t devzd, Double_t devaz);
|
---|
77 |
|
---|
78 | ClassDef(MAstro, 0)
|
---|
79 | };
|
---|
80 |
|
---|
81 | #endif
|
---|