Changeset 1273 for trunk/MagicSoft/Cosy/gui/MGSkyPosition.cc

Timestamp:

04/11/02 13:20:10 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Cosy/gui/MGSkyPosition.cc (modified) (13 diffs)

trunk/MagicSoft/Cosy/gui/MGSkyPosition.cc

@@ -16 +16 @@
 
 #include "coord.h"
+#include "SlaStars.h"
 #include "SlaPlanets.h"
 
 ClassImp(MGSkyPosition);
+
+void MGSkyPosition::InitArc(TArc *arc, Int_t fillstyle, Int_t fillcolor, Int_t linecolor)
+{
+    arc->SetFillStyle(fillstyle);      // (s. TAttFill)
+    arc->SetFillColor(fillcolor);      // (s. TAttFill)
+    arc->SetLineColor(linecolor);
+    arc->Draw();
+}
 
 void MGSkyPosition::InitPlanets()
@@ -34 +43 @@
     //   Neptun     7.8
     //
-    fSlalib = new SlaPlanets;
+    fSlaPlanet = new SlaPlanets;
+    fSlaStar   = new SlaStars;
 
     //
@@ -41 +51 @@
     //   1, 12, 13, 14, 15, 16, 17, 18, 19, 10
     //
-    fMars = new TArc(0, 0, 1);
-    fMars->SetFillStyle(1001);     // solid      (s. TAttFill)
-    fMars->SetFillColor(13);       // light grey (s. TAttFill)
-    fMars->SetLineColor(12);
-    fMars->Draw();
-
-    fVenus = new TArc(0, 0, 2);
-    fVenus->SetFillStyle(1001);    // solid      (s. TAttFill)
-    fVenus->SetFillColor(15);      // light grey (s. TAttFill)
-    fVenus->SetLineColor(14);
-    fVenus->Draw();
-
+    fMars    = new TArc(0, 0, 1);
+    fVenus   = new TArc(0, 0, 2);
     fJupiter = new TArc(0, 0, 2);
-    fJupiter->SetFillStyle(1001);  // solid      (s. TAttFill)
-    fJupiter->SetFillColor(17);    // light grey (s. TAttFill)
-    fJupiter->SetLineColor(16);
-    fJupiter->Draw();
-
-    fMoon = new TArc(0, 0, 3, 290, 70);
-    fMoon->SetFillStyle(1001);     // solid      (s. TAttFill)
-    fMoon->SetFillColor(19);       // light grey (s. TAttFill)
-    fMoon->SetLineColor(15);
-    fMoon->Draw();
+    fMoon    = new TArc(0, 0, 3, 290, 70);
+    fSun     = new TArc(0, 0, 2);
+
+    InitArc(fMars,    1001, 13, 12);
+    InitArc(fVenus,   1001, 15, 14);
+    InitArc(fJupiter, 1001, 17, 16);
+    InitArc(fMoon,    1001, 19, 15);
+    InitArc(fSun,     1001, 10, 19);
 
     for (int i=0; i<4; i++)
@@ -74 +72 @@
         fList->Add(fSunL[i]);
     }
-
-    fSun = new TArc(0, 0, 2);
-    fSun->SetFillStyle(1001);     // solid      (s. TAttFill)
-    fSun->SetFillColor(10);       // white      (s. TAttFill)
-    fSun->SetLineColor(19);
-    fSun->Draw();
 
     fList->Add(fMars);
@@ -168 +160 @@
     fList->Add(fLin1b);
     fList->Add(fLin2);
+
+    for (int i=0; i<6; i++)
+    {
+        fDot[i] = new TArc(0, 0, 1);
+        InitArc(fDot[i], 1001, 10, 19);
+        fList->Add(fDot[i]);
+    }
 }
 
@@ -173 +172 @@
     : MGEmbeddedCanvas("SkyPosition", p, w, 110)
 {
+    cout << "MGSkyPosition: Initilizing." << endl;
     DrawCoordinateSystem();
     InitPlanets();
@@ -178 +178 @@
     InitText();
     InitCanvas();
+    cout << "MGSkyPosition: Initilizing done." << endl;
 }
 
 MGSkyPosition::~MGSkyPosition()
 {
-    delete fSlalib;
+    delete fSlaPlanet;
+    delete fSlaStar;
 
     cout << "MGSkyPosition destroyed." << endl;
@@ -203 +205 @@
 }
 
-void MGSkyPosition::UpdatePosition(Float_t zd, Float_t az)
-{
+void MGSkyPosition::SetDot(TArc *arc, RaDec &radec, Int_t off)
+{
+    RaDec rd(radec.Ra()+off*360/24*kDeg2Rad, radec.Dec());
+    ZdAz zdaz = fSlaStar->CalcZdAz(rd);
+
+    const float s = sin(zdaz.Az());
+    const float c = cos(zdaz.Az());
+
+    const float x = s*zdaz.Zd()*kRad2Deg;             // [deg]
+    const float y = c*zdaz.Zd()*kRad2Deg;             // [deg]
+
+    arc->SetX1(x);
+    arc->SetY1(y);
+
+    SetModified();
+}
+
+void MGSkyPosition::UpdatePosition(RaDec &radec, Float_t zd, Float_t az)
+{
+    //
+    // zd [deg]
+    // az [deg]
+    //
     static int X =~0;
     static int Y =~0;
@@ -210 +233 @@
     static int Ry=~0;
 
-    const float rad = az*kDeg2Rad;
+    const float rad = az*kDeg2Rad;    // [rad]
 
     const float s = sin(rad);
     const float c = cos(rad);
 
-    const float x = s*zd;
-    const float y = c*zd;
+    const float x = s*zd;             // [deg]
+    const float y = c*zd;             // [deg]
 
     const int pixx = (int)(x/fPix);
@@ -238 +261 @@
     SetLin1(s*95., c*95.);
     SetLin2(x+dy, y-dx, x-dy, y+dx);
+
+    if (zd<80)
+    {
+        for(int i=0; i<3; i++)
+            SetDot(fDot[i], radec, i-3);
+        for(int i=3; i<6; i++)
+            SetDot(fDot[i], radec, i-2);
+    }
 
     SetModified();
@@ -279 +310 @@
     // calculate actual time for planet positions
     //
-    fSlalib->SetMjd2Now();
+    fSlaPlanet->SetMjd2Now();
+    fSlaStar->SetMjd2Now();
 
     UpdatePlanet(kESun,     fSun);
@@ -287 +319 @@
     UpdatePlanet(kEMars,    fMars);
 
-    UpdatePosition(pos.Zd(), pos.Az());
+    RaDec radec = fSlaStar->CalcRaDec(pos*kDeg2Rad);
+
+    UpdatePosition(radec, pos.Zd(), pos.Az());
     UpdateText(pos.Zd(), pos.Az());
 
@@ -298 +332 @@
     // calculate actual position of planet
     //
-    fSlalib->UpdatePlanetPos((ePlanets_t)planet);
-
-    ZdAz zdaz = fSlalib->GetPlanetPos((ePlanets_t)planet);
+    fSlaPlanet->UpdatePlanetPos((ePlanets_t)planet);
+
+    ZdAz zdaz = fSlaPlanet->GetPlanetPos((ePlanets_t)planet);
 
     //