Changeset 3497 for trunk/MagicSoft/Mars/mastro

Timestamp:

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

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mastro

Files:

• MObservatory.cc (modified) (2 diffs)
• MObservatory.h (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());
+}

trunk/MagicSoft/Mars/mastro/MObservatory.h

@@ -5 +5 @@
 #include "MParContainer.h"
 #endif
+
+class MTime;
 
 class MObservatory : public MParContainer
@@ -56 +58 @@
     void RotationAngle(Double_t theta, Double_t phi, Double_t &sin, Double_t &cos) const;
     Double_t RotationAngle(Double_t theta, Double_t phi) const;
+    Double_t RotationAngle(Double_t ra, Double_t dec, const MTime &t) const;
 
     LocationName_t GetObservatoryKey() const { return fObservatoryKey; }