Changeset 3808 for trunk/MagicSoft/Mars

Timestamp:

04/22/04 21:08:57 (21 years ago)

Author:

jlopez

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• mastro/MAstroCamera.cc (modified) (3 diffs)
• mastro/MAstroCamera.h (modified) (3 diffs)
• mtemp/MFindStars.cc (added)
• mtemp/MFindStars.h (added)
• mtemp/MStarLocalCam.cc (added)
• mtemp/MStarLocalCam.h (added)
• mtemp/MStarLocalPos.cc (added)
• mtemp/MStarLocalPos.h (added)

trunk/MagicSoft/Mars/mastro/MAstroCamera.cc

@@ -89 +89 @@
                                // HOW TO GET RID OF IT? Move MHCamera to mgeom?
 
+#include "MStarLocalPos.h"
+
 ClassImp(MAstroCamera);
 
@@ -116 +118 @@
 
     delete fMirror0;
+
 }
 
@@ -335 +338 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Options:
+//
+//  '*' Draw the mean of the reflections on all mirrors
+//  '.' 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::FillStarList()
+{
+    fStars->SetOwner();
+    fStars->Delete();
+
+    if (!fTime || !fObservatory || !fMirrors || !fStars)
+    {
+        cout << "Missing data..." << endl;
+        return;
+    }
+
+    // Get camera
+    MHCamera *camera=(MHCamera*)FindObjectInPad("MHCamera", gPad);
+    if (camera)
+    {
+        if (!camera->GetGeometry() || camera->GetGeometry()->IsA()!=fGeom->IsA())
+            camera->SetGeometry(*fGeom);
+    }
+    else
+    {
+        camera = new MHCamera(*fGeom);
+        camera->SetName("MHCamera");
+        camera->SetStats(0);
+        camera->SetInvDeepBlueSeaPalette();
+        camera->SetBit(kCanDelete);
+        camera->Draw();
+    }
+
+    const TRotation rot(GetGrid(kTRUE));
+
+    MVector3 *radec;
+    TIter Next(&fList);
+
+    while ((radec=(MVector3*)Next()))
+    {
+        const Double_t mag = radec->Magnitude();
+
+        TVector3 mean;
+        Int_t num = 0;
+
+        TVector3 star(*radec);
+        star *= rot;
+        const TVector3 spot = fMirror0->GetReflection(star, fGeom->GetCameraDist())*1000;
+
+        MStarLocalPos *starpos = new MStarLocalPos;
+        starpos->SetExpValues(mag,mean(0),mean(1));
+        fStars->Add(starpos);
+    }
+}
+
 // ------------------------------------------------------------------------
 //

trunk/MagicSoft/Mars/mastro/MAstroCamera.h

@@ -7 +7 @@
 
 class TClonesArray;
+class TList;
 
 class MTime;
@@ -18 +19 @@
     MGeomCam     *fGeom;
     TClonesArray *fMirrors;
+    TList        *fStars;
 
     MGeomMirror  *fMirror0;     //!
@@ -32 +34 @@
     void SetMirrors(TClonesArray &arr);
     void SetGeom(const MGeomCam &cam);
+    void SetStarList(TList* s) { fStars = s; }
+
+    void FillStarList();
+
+    TList *GetList() const { return fStars; }
 
     ClassDef(MAstroCamera, 1) // Display class to display stars on the camera