Changeset 3497 for trunk/MagicSoft/Mars/mastro/MObservatory.cc

Timestamp:

03/15/04 12:04:09 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mastro/MObservatory.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/mastro/MObservatory.cc

@@ -34 +34 @@
 #include "MObservatory.h"
 
+#include <TVector3.h>
+
+#include "MTime.h"
 #include "MAstro.h"
 
@@ -169 +172 @@
     return TMath::ASin((fCosLatitude*sinp) / denom);
 }
+
+// --------------------------------------------------------------------------
+//
+// RotationAngle 
+// 
+// calculates the angle for the rotation of the sky image in the camera;
+// this angle is a function of the sky coordinates, the observatory
+// location and the time
+//
+//  ra  [rad]: Right ascension
+//  dec [rad]: Declination
+//
+// Return RotationAngle in rad
+//
+Double_t MObservatory::RotationAngle(Double_t ra, Double_t dec, const MTime &t) const
+{
+    const Double_t alpha = t.GetGmst() + GetElong();
+
+    TVector3 v;
+    v.SetMagThetaPhi(1, TMath::Pi()/2-dec, alpha-ra);
+    v.RotateY(GetPhi()-TMath::Pi()/2);
+
+    return RotationAngle(v.Theta(), v.Phi());
+}