Changeset 4076 for trunk/MagicSoft/Cosy/gui/MGAccuracy.cc

Timestamp:

05/15/04 16:46:27 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Cosy/gui/MGAccuracy.cc (modified) (9 diffs)

trunk/MagicSoft/Cosy/gui/MGAccuracy.cc

@@ -15 +15 @@
 #include <TList.h>
 #include <TGaxis.h>
+#include <TGraph.h>
 #include <TCanvas.h>
 
 #include "coord.h"
+#include "MTime.h"
 
 ClassImp(MGAccuracy);
@@ -28 +30 @@
     TText text;
     text.SetTextAlign(22);  // centered, centered (s.TAttText)
-    text.DrawText(-80*2, 132.5*2, "Tracking Error [']");
+    text.DrawText(-80*2, 132.5*2, "Ctrl Deviation [min]");
 
 
@@ -75 +77 @@
     axe->SetBit(kCanDelete);
     axe->Draw();
-
-    //
-    // FIXME? Use TAxis?
-    //
-/*
-    TLine line;
-    line.SetLineColor(13);
-    line.SetLineStyle(3);  // dotted  (s. TAttLine)
-    line.DrawLine(-30, -65, -30,  65);
-    line.DrawLine(-65, -30,  65, -30);
-    line.DrawLine( 30, -65,  30,  65);
-    line.DrawLine(-65,  30,  65,  30);
-
-    line.DrawLine(-15, -65, -15,  65);
-    line.DrawLine(-65, -15,  65, -15);
-    line.DrawLine( 15, -65,  15,  65);
-    line.DrawLine(-65,  15,  65,  15);
-
-    line.DrawLine(-45, -65, -45,  65);
-    line.DrawLine(-65, -45,  65, -45);
-    line.DrawLine( 45, -65,  45,  65);
-    line.DrawLine(-65,  45,  65,  45);
-
-    line.SetLineColor(12);
-    line.SetLineStyle(2);  // dashed  (s. TAttLine)
-    line.DrawLine(-60, -65, -60,  65);
-    line.DrawLine(-65, -60,  65, -60);
-    line.DrawLine( 60, -65,  60,  65);
-    line.DrawLine(-65,  60,  65,  60);
-
-    line.SetLineColor(1);  // black
-    line.SetLineStyle(1);  // solid  (s. TAttLine)
-    line.DrawLine(-65,   0, 65,  0);
-    line.DrawLine(  0, -65,  0, 65);
-
-    line.DrawLine(-1,  60, 1,  60);
-    line.DrawLine(-1, -60, 1, -60);
-    line.DrawLine(-1,  30, 1,  30);
-    line.DrawLine(-1, -30, 1, -30);
-
-    line.DrawLine( 60, -1,  60, 1);
-    line.DrawLine(-60, -1, -60, 1);
-    line.DrawLine( 30, -1,  30, 1);
-    line.DrawLine(-30, -1, -30, 1);
-
-    TText text;
-    text.SetTextAlign(22);  // centered, centered (s.TAttText)
-    text.DrawText(60, 5, "dAz[\xb0]");
-    text.DrawText(0, 70, "dZd[\xb0]");
-
-    text.SetTextAlign(23);  // centered, centered (s.TAttText)
-    text.DrawText(-60, -2, "-1'");
-    text.DrawText( 60, -2, "1'");
-    text.DrawText(-30, -2, "-30\"");
-    text.DrawText( 30, -2, "30\"");
-
-    text.SetTextAlign(32);  // centered, centered (s.TAttText)
-    text.DrawText(-2, -60, "-1'");
-    text.DrawText(-2,  60, "1'");
-    text.DrawText(-2, -30, "-30\"");
-    text.DrawText(-2,  30, "30\"");
-    */
 }
 
@@ -183 +123 @@
 : MGEmbeddedCanvas("Accuracy", p, w, 300)
 {
+    // FIXME: Overload MapWindow in newer Root versions to remove
+    //        the contents of the graph!
+    fGraph = new TGraph;
+    fGraph->SetPoint(0, 0, 0);
+    fGraph->SetLineColor(kBlue);
+    fGraph->SetMarkerColor(kBlue);
+    fGraph->SetMarkerStyle(kFullDotMedium);
+    fGraph->Draw("LP");
+    fList->Add(fGraph);
+    //fGraph->SetNameTitle("AccVsT", "Accuracy vs Min of Time");
+    //fGraph->Draw("APL");
+    //fGraph->SetMarkerSize(2);
+
     DrawCoordinateSystem();
 
@@ -192 +145 @@
 
     SetNoContextMenu();
+
+    MTime t(-1);
+    fTime = t.GetAxisTime();
 }
 
@@ -199 +155 @@
 }
 
-void MGAccuracy::UpdateText(Float_t pzd, Float_t azd, Float_t aaz)
-{
-    const Float_t d2r = TMath::Pi()/180.;
-
-    pzd *= d2r;
-    azd *= d2r;
-    aaz *= d2r;
-
-    const double el = TMath::Pi()/2-pzd;
-
-    const double dphi2 = aaz/2.;
-    const double cos2  = cos(dphi2)*cos(dphi2);
-    const double sin2  = sin(dphi2)*sin(dphi2);
-    const double d     = cos(azd)*cos2 - cos(2*el)*sin2;
-
-    //
-    // Original:
-    //   cos(Zd1)*cos(Zd2)+sin(Zd1)*sin(Zd2)*cos(dAz)
-    //
-    // Correct:
-    //   const double d = cos(azd)*cos2 - cos(el1+el2)*sin2;
-    //
-    // Estimated:
-    //   const double d = cos(azd)*cos2 - cos(2*el)*sin2;
-    //
-
-    double dist = acos(d);
-
-    dist *= 3600./d2r; // [min]
+// dist [deg]
+void MGAccuracy::UpdateText(Float_t dist)
+{
+    dist *= 3600.; // [sec]
 
     int rs = (int)floor(fmod(dist, 60.));
 
-    dist /= 60.;
-    int rm = (int)dist;//floor(fmod(dist, 60.));
+    dist /= 60.;   // [min]
+    int rm = (int)dist;
 
     char txt[100];
@@ -239 +170 @@
     fTxt->SetText(fTxt->GetX(), fTxt->GetY(), txt);
 
-    fBar->SetX2(dist*60);
+    fBar->SetY2(dist*60); // [sec]
     if (dist*16384<1*360*60)
         fBar->SetLineColor(kGreen);
@@ -247 +178 @@
         else
             fBar->SetLineColor(kRed);
+
+    SetModified();
+}
+
+// dist [deg]
+void MGAccuracy::UpdateGraph(Float_t dist)
+{
+    MTime t(-1);
+    const Double_t dtime = t.GetAxisTime()-fTime; // range [-0.5h, 0h]
+
+    dist *= 60; // min
+
+    static int odist = -1;
+    if (odist==(int)(dist*10*60) && dtime<10)
+        return;
+
+    odist = (int)(dist*10*60);
+
+    fGraph->SetPoint(fGraph->GetN(), dtime, dist*60);
+
+    const Double_t ntime = dtime;
+    for (int i=0; i<fGraph->GetN(); i++)
+    {
+        Double_t x, y;
+        fGraph->GetPoint(i, x,       y);
+        fGraph->SetPoint(i, x-ntime, y);
+        //cout << i << ":  " << x-ntime << " / " << y << endl;
+    }
+    while (fGraph->GetN()>0)
+    {
+        Double_t x, y;
+        fGraph->GetPoint(0, x, y);
+
+        if (x==-ntime && y==0)
+        {
+            fGraph->RemovePoint(0);
+            continue;
+        }
+
+        if (x>-4.75*60)
+            break;
+
+        fGraph->RemovePoint(0);
+    }
+
+    fTime = t.GetAxisTime();
+
+    SetModified();
+
+    //cout << "N1 == " << fGraph->GetN() << endl;
+
+    //fGraph->GetHistogram()->SetXTitle("Time");
+    //fGraph->GetHistogram()->SetYTitle("\\Delta [arcmin]");
+    //fGraph->GetHistogram()->GetXaxis()->SetTimeFormat("%M %F1995-01-01 00:00:00");
+    //fGraph->GetHistogram()->GetXaxis()->SetTimeDisplay(1);
+    //fGraph->GetHistogram()->GetXaxis()->SetLabelSize(0.033);
+}
+
+void MGAccuracy::Update(Float_t pzd, Float_t azd, Float_t aaz)
+{
+    const Float_t d2r = TMath::Pi()/180.;
+
+    pzd *= d2r;
+    azd *= d2r;
+    aaz *= d2r;
+
+    const double el = TMath::Pi()/2-pzd;
+
+    const double dphi2 = aaz/2.;
+    const double cos2  = cos(dphi2)*cos(dphi2);
+    const double sin2  = sin(dphi2)*sin(dphi2);
+    const double d     = cos(azd)*cos2 - cos(2*el)*sin2;
+
+    //
+    // Original:
+    //   cos(Zd1)*cos(Zd2)+sin(Zd1)*sin(Zd2)*cos(dAz)
+    //
+    // Correct:
+    //   const double d = cos(azd)*cos2 - cos(el1+el2)*sin2;
+    //
+    // Estimated:
+    //   const double d = cos(azd)*cos2 - cos(2*el)*sin2;
+    //
+
+    double dist = acos(d)*TMath::RadToDeg();
+
+    UpdateText(dist);
+    UpdateGraph(dist);
 }
 
 void MGAccuracy::UpdateCross(Float_t x, Float_t y)
 {
+    //
+    // calculate actual time for planet positions
+    // acc [deg]
+    //
+    // x["], y["]
+    //
+    static int X = ~0;
+    static int Y = ~0;
+
+    int pixx = (int)(x/fPix);  // [pix]
+    int pixy = (int)(y/fPix);  // [pix]
+
+    if (X==pixx && Y==pixy)
+        return;
+
+    X = pixx;
+    Y = pixy;
+
     fLin1->SetX1(x-5.);
     fLin1->SetX2(x+5.);
@@ -262 +299 @@
     fLin2->SetY1(y+5.);
     fLin2->SetY2(y-5.);
+
+    SetModified();
 }
 
 void MGAccuracy::Update(ZdAz &pos, ZdAz &acc)
 {
-    //
-    // calculate actual time for planet positions
-    // acc [deg]
-    //
-    static int X = ~0;
-    static int Y = ~0;
-
-    float x = acc.Az()*3600.; // ["]
-    float y = acc.Zd()*3600.; // ["]
-
-    int pixx = (int)(x/fPix);  // [pix]
-    int pixy = (int)(y/fPix);  // [pix]
-
-    if (X==pixx && Y==pixy)
-        return;
-
-    X = pixx;
-    Y = pixy;
-
-    UpdateCross(x, y);
-    UpdateText(pos.Zd(), acc.Zd(), acc.Az());
-
-    SetModified();
+    UpdateCross(acc.Az()*3600., acc.Zd()*3600.);
+    Update(pos.Zd(), acc.Zd(), acc.Az());
+
     UpdateCanvas();
 }