Changeset 3704 for trunk/MagicSoft/Mars/mastro/MAstroCamera.cc

Timestamp:

04/10/04 16:02:20 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mastro/MAstroCamera.cc (modified) (10 diffs)

trunk/MagicSoft/Mars/mastro/MAstroCamera.cc

@@ -33 +33 @@
 // PRELIMINARY!!
 //
+// The caluclation of the position of the reflection in the camera is
+// done by:
+//   - Rotation of the star-field such that the camera is looking into
+//     the pointing direction
+//   - Calculation of the reflected star-light vector by calling
+//     MGeomMirror::GetReflection (which returns the point at which
+//     the vector hits the camera plain)
+//   - Depending on the draw-option you get each reflected point, the
+//     reflection on a virtual ideal mirror or the reflection on each
+//     individual mirror
+//
+// GUI: You can use the the cursor keys to change the pointing position
+//      and plus/minus to change the time by a quarter of an hour.
+//
 /////////////////////////////////////////////////////////////////////////////
 #include "MAstroCamera.h"
@@ -58 +72 @@
 
 // --------------------------------------------------------------------------
+//
+// Create a virtual MGeomMirror which is in the center of the coordinate
+// system and has a normal vector in z-direction.
+//
 MAstroCamera::MAstroCamera() : fGeom(0), fMirrors(0)
 {
@@ -65 +83 @@
 
 // --------------------------------------------------------------------------
+//
+// Delete fGeom, fMirrors and the virtual 0-Mirror fMirror0
+//
 MAstroCamera::~MAstroCamera()
 {
@@ -76 +97 @@
 
 // --------------------------------------------------------------------------
+//
+// Set a list of mirrors. The Mirrors must be of type MGeomMirror and
+// stored in a TClonesArray
+//
 void MAstroCamera::SetMirrors(TClonesArray &arr)
 {
     if (arr.GetClass()!=MGeomMirror::Class())
+    {
+        cout << "ERROR - TClonesArray doesn't contain objects of type MGeomMirror... ignored." << endl;
         return;
+    }
 
     const Int_t n = arr.GetSize();
@@ -94 +122 @@
 
 // --------------------------------------------------------------------------
+//
+// Set the camera geometry. The MGeomCam object is cloned.
+//
 void MAstroCamera::SetGeom(const MGeomCam &cam)
 {
@@ -103 +134 @@
 
 // --------------------------------------------------------------------------
+//
+// Convert To Pad coordinates (see MAstroCatalog)
+//
 Int_t MAstroCamera::ConvertToPad(const TVector3 &w, TVector2 &v) const
 {
@@ -124 +158 @@
 
 // --------------------------------------------------------------------------
+//
+// Find an object with a given name in the list of primitives of this pad.
+//
 TObject *FindObjectInPad(const char *name, TVirtualPad *pad)
 {
@@ -152 +189 @@
 //
 //  '*' Draw the mean of the reflections on all mirrors
-//  '.' Draw a dot for the reflection on each mirror
+//  '.' Draw a dot for the reflection on each individual mirror
+//  'h' To create a TH2D of the star-light which is displayed
+//  'c' Use the underlaying MHCamera as histogram
+//  '0' Draw the reflection on a virtual perfect mirror
+//
+// If the Pad contains an object MHCamera of type MHCamera it is used.
+// Otherwise a new object is created.
 //
 void MAstroCamera::AddPrimitives(TString o)
@@ -165 +208 @@
         o = "*.";
 
+    const Bool_t hasnull = o.Contains("0", TString::kIgnoreCase);
     const Bool_t hashist = o.Contains("h", TString::kIgnoreCase);
     const Bool_t hasmean = o.Contains("*", TString::kIgnoreCase);
@@ -257 +301 @@
         // transform meters into millimeters (camera display works with mm)
         mean *= 1./num;
-
         DrawStar(mean(0), mean(1), *radec, !hasmean, Form("x=%.1fmm y=%.1fmm", mean(0), mean(1)));
+
+        if (hasnull)
+        {
+            TVector3 star(*radec);
+            star *= rot;
+            const TVector3 spot = fMirror0->GetReflection(star, fGeom->GetCameraDist())*1000;
+            DrawStar(spot(0), spot(1), *radec, !hasmean, Form("x=%.1fmm y=%.1fmm", mean(0), mean(1)));
+        }
     }
 }