Changeset 9555

Timestamp:

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

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Cosy/main

Files:

• MMoonPointing.cc (modified) (1 diff)
• MTracking.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:

trunk/MagicSoft/Cosy/main/MTracking.cc

@@ -239 +239 @@
         ZdAz srcpos, pointpos;
 
-        const ZdAz za     = sla.GetZdAzRad();
-        const ZdAz offset = moon.CalcPosition(za, srcpos, pointpos);
-        sla.ApplyOffsetZdAz(pointpos-za);
+        const ZdAz za = sla.GetZdAzRad();
+        if (moon.CalcPosition(za, srcpos, pointpos))
+            sla.Set(pointpos);
 
         // Ra/Dec, Zd/Az from pointpos