File tree Expand file tree Collapse file tree 1 file changed +6
-6
lines changed
Expand file tree Collapse file tree 1 file changed +6
-6
lines changed Original file line number Diff line number Diff line change @@ -38,7 +38,6 @@ NOTE: the use of time_t (and its default implementation as a 32 int)
3838*/
3939
4040#include " TTimeStamp.h"
41- #include " TMath.h"
4241#include " TString.h"
4342#include " TError.h"
4443#include " Riostream.h"
@@ -191,12 +190,13 @@ Double_t TTimeStamp::AsGMST(Double_t UT1Offset) const
191190
192191Double_t TTimeStamp::AsGAST (Double_t UT1Offset) const
193192{
193+ Double_t Pi = 3.14159265358979323846 ;
194194 Double_t D = (AsJulianDate () + UT1Offset/86400000.0 ) - 2451545.0 ;
195- Double_t epsilon = (23.4393 - 0.0000004 *D)* TMath::Pi ()/ 180.0 ;
196- Double_t L = (280.47 + 0.98565 *D)* TMath::Pi ()/ 180.0 ;
197- Double_t Omega = (125.04 - 0.052954 *D)* TMath::Pi ()/ 180.0 ;
198- Double_t Deltapsi = -0.000319 * TMath::Sin (Omega) - 0.000024 * TMath::Sin (2.0 * L);
199- Double_t eqeq = Deltapsi* TMath::Cos (epsilon);
195+ Double_t epsilon = (23.4393 - 0.0000004 * D) * Pi / 180.0 ;
196+ Double_t L = (280.47 + 0.98565 * D) * Pi / 180.0 ;
197+ Double_t Omega = (125.04 - 0.052954 * D) * Pi / 180.0 ;
198+ Double_t Deltapsi = -0.000319 * std::sin (Omega) - 0.000024 * std::sin (2.0 * L);
199+ Double_t eqeq = Deltapsi * std::cos (epsilon);
200200 Double_t rval = fmod (AsGMST (UT1Offset) + eqeq, 24.0 );
201201 return rval < 0 ? rval + 24.0 : rval;
202202}
You can’t perform that action at this time.
0 commit comments