Changeset 9555 for trunk/MagicSoft/Cosy/main/MMoonPointing.cc

Timestamp:

02/25/10 11:30:27 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Cosy/main/MMoonPointing.cc (modified) (1 diff)

trunk/MagicSoft/Cosy/main/MMoonPointing.cc

@@ -128 +128 @@
     const Double_t fcos = const_cast<TGraph2D&>(fCos).Interpolate(moon.Az(), moon.Zd());
     const Double_t fsin = const_cast<TGraph2D&>(fSin).Interpolate(moon.Az(), moon.Zd());
-    if (fcos*fcos + fsin*fsin<0.98)
+
+    // This is a sanity check for the case the Interpolation failed
+    const Double_t norm = TMath::Hypot(fcos, fsin);
+    if (TMath::Abs(norm-1)<0.001)
         return kFALSE;
+
+    // Just a small numerical improvement
+    fsin /= norm;
+    fcos /= norm;
 
     // Calculate Moon Shadow Position: