Changeset 4966 for trunk/MagicSoft/Mars/mastro

Timestamp:

09/13/04 08:57:58 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mastro

Files:

• MAstro.cc (modified) (2 diffs)
• MAstro.h (modified) (2 diffs)

trunk/MagicSoft/Mars/mastro/MAstro.cc

@@ -36 +36 @@
 
 #include "MTime.h"    // MTime::GetGmst
+
+#include "MAstroCatalog.h" // FIXME: replace by MVector3!
 
 using namespace std;
@@ -540 +542 @@
     return (Int_t)TMath::Floor(period)-284;
 }
+
+// --------------------------------------------------------------------------
+//
+// Returns the distance in x,y between two polar-vectors (eg. Alt/Az, Ra/Dec)
+// projected on aplain in a distance dist. For Magic this this the distance
+// of the camera plain (1700mm) dist also determins the unit in which
+// the TVector2 is returned.
+//
+// v0 is the reference vector (eg. the vector to the center of the camera)
+// v1 is the vector to which we determin the distance on the plain
+//
+//  (see also MStarCamTrans::Loc0LocToCam())
+//
+TVector2 MAstro::GetDistOnPlain(const TVector3 &v0, TVector3 v1, Double_t dist)
+{
+    v1.RotateZ(-v0.Phi());
+    v1.RotateY(-v0.Theta());
+    v1.RotateZ(-TMath::Pi()/2); // exchange x and y
+    v1 *= dist/v1.Z();
+
+    return v1.XYvector(); //TVector2(v1.Y(), -v1.X());//v1.XYvector();
+}

trunk/MagicSoft/Mars/mastro/MAstro.h

@@ -5 +5 @@
 #include <TROOT.h>
 #endif
+#ifndef ROOT_TVector2
+#include <TVector2.h>
+#endif
+
+class TVector3;
 
 class MAstro
@@ -60 +65 @@
     static Int_t    GetMagicPeriod(Double_t mjd);
 
+    static TVector2 GetDistOnPlain(const TVector3 &v0, TVector3 v1, Double_t dist=1);
+
     ClassDef(MAstro, 0)
 };