source: tags/Mars-V0.10.2/mastro/MAstro.h

Last change on this file was 8066, checked in by tbretz, 18 years ago
*** empty log message ***
File size: 4.3 KB
Line 
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
11class TArrayD;
12class TVector3;
13class MTime;
14
15class MAstro
16{
17private:
18 static Double_t Round(Double_t val);
19 static Double_t Trunc(Double_t val);
20
21 static Double_t Kepler(Double_t m, Double_t ecc);
22
23public:
24 virtual ~MAstro() { }
25
26 // Unit conversion
27 static Double_t Parsec2Meter() { return 3.0856776e+16; }
28 static Double_t Parsec2Lighyear() { return 3.2615644; }
29
30
31 // Angle/Time conversion functions
32 static Double_t Hms2Sec(Int_t deg, UInt_t min, Double_t sec, char sgn='+');
33 static Double_t Dms2Rad(Int_t deg, UInt_t min, Double_t sec, Char_t sgn='+');
34 static Double_t Hms2Rad(Int_t hor, UInt_t min, Double_t sec, Char_t sgn='+');
35 static Double_t Dms2Deg(Int_t deg, UInt_t min, Double_t sec, Char_t sgn='+');
36 static Double_t Hms2Deg(Int_t hor, UInt_t min, Double_t sec, Char_t sgn='+');
37 static Double_t Dms2Hor(Int_t deg, UInt_t min, Double_t sec, Char_t sgn='+');
38 static Double_t Hms2Hor(Int_t hor, UInt_t min, Double_t sec, Char_t sgn='+');
39
40 static void Day2Hms(Double_t day, Char_t &sgn, UShort_t &hor, UShort_t &min, UShort_t &sec);
41 static void Rad2Dms(Double_t rad, Char_t &sgn, UShort_t &deg, UShort_t &min, UShort_t &sec);
42 static void Rad2Hms(Double_t rad, Char_t &sgn, UShort_t &hor, UShort_t &min, UShort_t &sec);
43 static void Deg2Dms(Double_t deg, Char_t &sgn, UShort_t &deg, UShort_t &min, UShort_t &sec);
44 static void Deg2Hms(Double_t deg, Char_t &sgn, UShort_t &hor, UShort_t &min, UShort_t &sec);
45 static void Hor2Dms(Double_t hor, Char_t &sgn, UShort_t &deg, UShort_t &min, UShort_t &sec);
46 static void Hor2Hms(Double_t hor, Char_t &sgn, UShort_t &hor, UShort_t &min, UShort_t &sec);
47
48 static void Day2Hm(Double_t day, Char_t &sgn, UShort_t &hor, Double_t &min);
49 static void Rad2Dm(Double_t rad, Char_t &sgn, UShort_t &deg, Double_t &min);
50 static void Rad2Hm(Double_t rad, Char_t &sgn, UShort_t &hor, Double_t &min);
51 static void Deg2Dm(Double_t deg, Char_t &sgn, UShort_t &deg, Double_t &min);
52 static void Deg2Hm(Double_t deg, Char_t &sgn, UShort_t &hor, Double_t &min);
53 static void Hor2Dm(Double_t rad, Char_t &sgn, UShort_t &deg, Double_t &min);
54 static void Hor2Hm(Double_t rad, Char_t &sgn, UShort_t &hor, Double_t &min);
55
56 // Print funtions
57 static TString GetStringDeg(Double_t deg, const char *fmt="%c%d:%02d:%02d");
58 static TString GetStringHor(Double_t deg, const char *fmt="%c%d:%02d:%02d");
59
60 // Angle treatment functions
61 static Bool_t String2Angle(TString &str, Double_t &ret);
62 static Bool_t Coordinate2Angle(const TString &str, Double_t &ret);
63 static TString Angle2Coordinate(Double_t val);
64
65 static Double_t AngularDistance(Double_t theta0, Double_t phi0, Double_t theta1, Double_t phi1);
66
67 // Time treatment functions
68 static void Mjd2Ymd(UInt_t mjd, UShort_t &y, Byte_t &m, Byte_t &d);
69 static Int_t Ymd2Mjd(UShort_t y, Byte_t m, Byte_t d);
70
71 static Double_t UT2GMST(Double_t ut1);
72
73 // Rotation angle between local and sky coordinate system
74 static void RotationAngle(Double_t sinl, Double_t cosl, Double_t theta, Double_t phi, Double_t &sin, Double_t &cos);
75 static Double_t RotationAngle(Double_t sinl, Double_t cosl, Double_t theta, Double_t phi);
76
77 static void GetMoonRaDec(Double_t mjd, Double_t &ra, Double_t &dec);
78
79 static Double_t GetMoonPhase(Double_t mjd);
80 static Double_t GetMoonPeriod(Double_t mjd);
81 static Int_t GetMagicPeriod(Double_t mjd);
82
83 // Estimate some parameters around the sun
84 static Double_t GetSunRaDec(Double_t mjd, Double_t &ra, Double_t &dec);
85 static Double_t GetEquationOfTime(Double_t mjd);
86 static Double_t GetNoonTime(Double_t mjd, Double_t longit);
87 static Double_t GetTimeFromNoonToAlt(Double_t mjd, Double_t latit, Double_t alt);
88 static TArrayD GetSunRiseSet(Double_t mjd, Double_t longit, Double_t latit, Double_t alt=0);
89
90 // Get distance between v1 and v0 in a plain perpendicular to v0 in distance dist
91 static TVector2 GetDistOnPlain(const TVector3 &v0, TVector3 v1, Double_t dist=1);
92
93 // Calculate absolute misspointing from error in zd and error in az
94 static Double_t GetDevAbs(Double_t nomzd, Double_t devzd, Double_t devaz);
95
96 static Int_t GetEasterOffset(UShort_t year);
97
98 ClassDef(MAstro, 0)
99};
100
101#endif
Note: See TracBrowser for help on using the repository browser.