Changeset 14851 for trunk/Mars

Timestamp:

02/04/13 09:53:50 (12 years ago)

Author:

tbretz

Message:

Made GetDevAbs more accurate.

File:

• trunk/Mars/mastro/MAstro.cc (modified) (1 diff)

trunk/Mars/mastro/MAstro.cc

@@ -919 +919 @@
 // All values given in deg, the return value, too.
 //
-Double_t MAstro::GetDevAbs(Double_t nomzd, Double_t devzd, Double_t devaz)
-{
-    const Double_t pzd = nomzd * TMath::DegToRad();
-    const Double_t azd = devzd * TMath::DegToRad();
-    const Double_t aaz = devaz * TMath::DegToRad();
-
-    const double el = TMath::Pi()/2-pzd;
-
-    const double dphi2 = aaz/2.;
-    const double cos2  = TMath::Cos(dphi2)*TMath::Cos(dphi2);
-    const double sin2  = TMath::Sin(dphi2)*TMath::Sin(dphi2);
-    const double d     = TMath::Cos(azd)*cos2 - TMath::Cos(2*el)*sin2;
-
-    return TMath::ACos(d)*TMath::RadToDeg();
+Double_t MAstro::GetDevAbs(Double_t nomzd, Double_t meszd, Double_t devaz)
+{
+    nomzd *= TMath::DegToRad();
+    meszd *= TMath::DegToRad();
+    devaz *= TMath::DegToRad();
+
+    const double x = sin(meszd) * sin(nomzd) * cos(dvaz);
+    const double y = cos(meszd) * cos(nomzd);
+
+    return acos(x + y) * TMath::RadToDeg();
 }