Changeset 14589 for trunk

Timestamp:

11/08/12 15:51:36 (12 years ago)

Author:

tbretz

Message:

The new algorithm didn't corrctly tak into account that it was el not zd

File:

• trunk/Cosy/tpoint/TPointStar.cc (modified) (1 diff)

trunk/Cosy/tpoint/TPointStar.cc

@@ -62 +62 @@
 Double_t TPointStar::GetResidual(Double_t *err) const
 {
-    /*
-     TVector3 v1, v2;
-     v1.SetMagThetaPhi(1, TMath::Pi()/2-fRawEl, fRawAz);
-     v2.SetMagThetaPhi(1, TMath::Pi()/2-fStarEl, fStarAz);
-
-     return v1.Angle(v2)*TMath::RadToDeg();
-     */
-
     const Double_t del = fRawEl-fStarEl;
     const Double_t daz = fRawAz-fStarAz;
 
-    const double x = sin(fRawEl) * sin(fStarEl) * cos(fStarAz-fRawAz);
-    const double y = cos(fRawEl) * cos(fStarEl);
+    const double x = cos(fRawEl) * cos(fStarEl) * cos(fStarAz-fRawAz);
+    const double y = sin(fRawEl) * sin(fStarEl);
 
-    const Double_t d = acos(x + y);
+    const Double_t d = x + y;
 
     if (err)