source: trunk/MagicSoft/Mars/mastro/MAstro.h@ 7566

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