1 | #include "MGAccuracy.h"
|
---|
2 |
|
---|
3 | #include <iostream.h> // cout
|
---|
4 |
|
---|
5 | #include <TArc.h>
|
---|
6 | #include <TLine.h>
|
---|
7 | #include <TText.h>
|
---|
8 | #include <TWbox.h>
|
---|
9 | #include <TGaxis.h>
|
---|
10 | #include <TGraph.h>
|
---|
11 | #include <TCanvas.h>
|
---|
12 |
|
---|
13 | #include "coord.h"
|
---|
14 | #include "MTime.h"
|
---|
15 | #include "MAstro.h"
|
---|
16 |
|
---|
17 | ClassImp(MGAccuracy);
|
---|
18 |
|
---|
19 | void MGAccuracy::DrawText(const char *c1, const char *c2, const char *c3, const char*txt)
|
---|
20 | {
|
---|
21 | fCanvas->SetEditable(kTRUE);
|
---|
22 | fCanvas->cd();
|
---|
23 |
|
---|
24 | TText text;
|
---|
25 | text.SetTextAlign(22); // centered, centered (s.TAttText)
|
---|
26 | text.DrawText(-150, 265, txt);
|
---|
27 |
|
---|
28 | text.SetTextAlign(11); // left, bottom (s.TAttText)
|
---|
29 | text.SetTextColor(3);
|
---|
30 | text.DrawText(220, -220, c1);
|
---|
31 |
|
---|
32 | text.SetTextColor(5);
|
---|
33 | text.DrawText(220, -250, c2);
|
---|
34 |
|
---|
35 | text.SetTextColor(2);
|
---|
36 | text.DrawText(220, -280, c3);
|
---|
37 |
|
---|
38 | fCanvas->SetEditable(kFALSE);
|
---|
39 | }
|
---|
40 |
|
---|
41 | void MGAccuracy::DrawCircles(Double_t rad)
|
---|
42 | {
|
---|
43 | fCanvas->SetEditable(kTRUE);
|
---|
44 | fCanvas->cd();
|
---|
45 |
|
---|
46 | TArc arc;
|
---|
47 | arc.SetFillStyle(4000); // transparent
|
---|
48 | arc.SetFillColor(39);
|
---|
49 | arc.SetLineColor(3); // green
|
---|
50 | arc.SetLineStyle(2); // dashed (s. TAttLine)
|
---|
51 | arc.DrawArc(0, 0, rad/2.); // 0.5se
|
---|
52 |
|
---|
53 | arc.SetLineColor(5); // yellow
|
---|
54 | arc.DrawArc(0, 0, rad); // 1.0se
|
---|
55 |
|
---|
56 | arc.SetLineColor(2); // red
|
---|
57 | arc.DrawArc(0, 0, rad*2.); // 2.0se
|
---|
58 |
|
---|
59 | fCanvas->SetEditable(kFALSE);
|
---|
60 | }
|
---|
61 |
|
---|
62 | void MGAccuracy::DrawCoordinateSystem()
|
---|
63 | {
|
---|
64 | TWbox box;
|
---|
65 | box.DrawWbox(-290, 290, -10, 240, 18, 2, 1);
|
---|
66 |
|
---|
67 | // TLine line;
|
---|
68 | // line.DrawLine(-65*4, 0, 65*4, 0);
|
---|
69 | // line.DrawLine( 0, -65*4, 0, 65*4);
|
---|
70 |
|
---|
71 | //
|
---|
72 | // Can be replaced by TGaxis axe; in a later root version
|
---|
73 | // than 3.01/06. I talked to Rene
|
---|
74 | //
|
---|
75 | TGaxis *axe;
|
---|
76 | axe = new TGaxis(-60*4, 0, 60*4, 0, -4, 4, 30204, "+-N");
|
---|
77 | axe->SetTitle("Az"); // \xb0
|
---|
78 | axe->SetBit(kCanDelete);
|
---|
79 | axe->Draw();
|
---|
80 |
|
---|
81 | axe = new TGaxis( 0, -60*4, 0, 60*4, -4, 4, 304, "+-N");
|
---|
82 | axe->SetTitle("Zd"); // \xb0
|
---|
83 | axe->SetLabelOffset(-0.02);
|
---|
84 | axe->SetBit(kCanDelete);
|
---|
85 | axe->Draw();
|
---|
86 | }
|
---|
87 |
|
---|
88 | void MGAccuracy::InitText()
|
---|
89 | {
|
---|
90 | fTxt = new TText(280, 280, "0' / 0'");
|
---|
91 | fTxt->SetTextAlign(33); // right, top
|
---|
92 | fTxt->SetTextColor(10);
|
---|
93 | fTxt->Draw();
|
---|
94 |
|
---|
95 | fList->Add(fTxt);
|
---|
96 | }
|
---|
97 |
|
---|
98 | void MGAccuracy::InitBar()
|
---|
99 | {
|
---|
100 | fBar = new TLine(0, 0, 0, 0);
|
---|
101 | fBar->SetLineColor(kBlack);
|
---|
102 | fBar->SetLineStyle(1);
|
---|
103 | fBar->SetLineWidth(5);
|
---|
104 | fBar->Draw();
|
---|
105 |
|
---|
106 | fList->Add(fBar);
|
---|
107 | }
|
---|
108 |
|
---|
109 | void MGAccuracy::InitCross()
|
---|
110 | {
|
---|
111 | fLin1 = new TLine(0, 0, 0, 0);
|
---|
112 | fLin2 = new TLine(0, 0, 0, 0);
|
---|
113 |
|
---|
114 | fLin1->SetLineColor(10); // white (s. TAttFill)
|
---|
115 | fLin2->SetLineColor(10); // white
|
---|
116 | fLin1->SetLineStyle(1); // solid (s. TAttLine)
|
---|
117 | fLin2->SetLineStyle(1);
|
---|
118 |
|
---|
119 | fLin1->SetLineWidth(2);
|
---|
120 | fLin2->SetLineWidth(2);
|
---|
121 |
|
---|
122 | fLin1->Draw();
|
---|
123 | fLin2->Draw();
|
---|
124 |
|
---|
125 | fList->Add(fLin1);
|
---|
126 | fList->Add(fLin2);
|
---|
127 | }
|
---|
128 |
|
---|
129 | MGAccuracy::MGAccuracy(const TGWindow* p, const UInt_t w)
|
---|
130 | : MGEmbeddedCanvas("Accuracy", p, w, 300), fColorScale(360*60/16384)
|
---|
131 | {
|
---|
132 | // FIXME: Overload MapWindow in newer Root versions to remove
|
---|
133 | // the contents of the graph!
|
---|
134 | fGraph = new TGraph;
|
---|
135 | fGraph->SetPoint(0, 0, 0);
|
---|
136 | fGraph->SetLineColor(kBlue);
|
---|
137 | fGraph->SetMarkerColor(kBlue);
|
---|
138 | fGraph->SetMarkerStyle(kFullDotMedium);
|
---|
139 | fGraph->Draw("LP");
|
---|
140 | fList->Add(fGraph);
|
---|
141 |
|
---|
142 | DrawCoordinateSystem();
|
---|
143 |
|
---|
144 | InitText();
|
---|
145 | InitCross();
|
---|
146 | InitBar();
|
---|
147 |
|
---|
148 | InitCanvas();
|
---|
149 |
|
---|
150 | SetNoContextMenu();
|
---|
151 |
|
---|
152 | MTime t(-1);
|
---|
153 | fTime = t.GetAxisTime();
|
---|
154 | }
|
---|
155 |
|
---|
156 | MGAccuracy::~MGAccuracy()
|
---|
157 | {
|
---|
158 | }
|
---|
159 |
|
---|
160 | // dist [deg]
|
---|
161 | void MGAccuracy::UpdateText(Float_t dist)
|
---|
162 | {
|
---|
163 | dist *= 3600.; // [sec]
|
---|
164 |
|
---|
165 | int rs = (int)floor(fmod(dist, 60.));
|
---|
166 |
|
---|
167 | dist /= 60.; // [min]
|
---|
168 | int rm = (int)dist;
|
---|
169 |
|
---|
170 | char txt[100];
|
---|
171 | rm ? sprintf(txt, "%d'%02d\"", rm, rs) : sprintf(txt, "%d\"", rs);
|
---|
172 |
|
---|
173 | fTxt->SetText(fTxt->GetX(), fTxt->GetY(), txt);
|
---|
174 |
|
---|
175 | fBar->SetY2(dist*60); // [sec]
|
---|
176 | if (dist<fColorScale) //1*360*60/16384
|
---|
177 | fBar->SetLineColor(kGreen);
|
---|
178 | else
|
---|
179 | if (dist<2*fColorScale) //2*360*60/16384
|
---|
180 | fBar->SetLineColor(kYellow);
|
---|
181 | else
|
---|
182 | fBar->SetLineColor(kRed);
|
---|
183 |
|
---|
184 | SetModified();
|
---|
185 | }
|
---|
186 |
|
---|
187 | // dist [deg]
|
---|
188 | void MGAccuracy::UpdateGraph(Float_t dist)
|
---|
189 | {
|
---|
190 | MTime t(-1);
|
---|
191 | const Double_t dtime = t.GetAxisTime()-fTime; // range [-0.5h, 0h]
|
---|
192 |
|
---|
193 | dist *= 60; // min
|
---|
194 |
|
---|
195 | static int odist = -1;
|
---|
196 | if (odist==(int)(dist*10*60) && dtime<10)
|
---|
197 | return;
|
---|
198 |
|
---|
199 | odist = (int)(dist*10*60);
|
---|
200 |
|
---|
201 | fGraph->SetPoint(fGraph->GetN(), dtime, dist*60);
|
---|
202 |
|
---|
203 | const Double_t ntime = dtime;
|
---|
204 | for (int i=0; i<fGraph->GetN(); i++)
|
---|
205 | {
|
---|
206 | Double_t x, y;
|
---|
207 | fGraph->GetPoint(i, x, y);
|
---|
208 | fGraph->SetPoint(i, x-ntime, y);
|
---|
209 | //cout << i << ": " << x-ntime << " / " << y << endl;
|
---|
210 | }
|
---|
211 | while (fGraph->GetN()>0)
|
---|
212 | {
|
---|
213 | Double_t x, y;
|
---|
214 | fGraph->GetPoint(0, x, y);
|
---|
215 |
|
---|
216 | if (x==-ntime && y==0)
|
---|
217 | {
|
---|
218 | fGraph->RemovePoint(0);
|
---|
219 | continue;
|
---|
220 | }
|
---|
221 |
|
---|
222 | if (x>-4.75*60)
|
---|
223 | break;
|
---|
224 |
|
---|
225 | fGraph->RemovePoint(0);
|
---|
226 | }
|
---|
227 |
|
---|
228 | fTime = t.GetAxisTime();
|
---|
229 |
|
---|
230 | SetModified();
|
---|
231 | }
|
---|
232 |
|
---|
233 | void MGAccuracy::Update(Float_t pzd, Float_t azd, Float_t aaz)
|
---|
234 | {
|
---|
235 | const double dist = MAstro::GetDevAbs(pzd, azd, aaz);
|
---|
236 |
|
---|
237 | UpdateText(dist);
|
---|
238 | UpdateGraph(dist);
|
---|
239 | }
|
---|
240 |
|
---|
241 | void MGAccuracy::UpdateCross(Float_t x, Float_t y)
|
---|
242 | {
|
---|
243 | //
|
---|
244 | // x["], y["]
|
---|
245 | //
|
---|
246 | static int X = ~0;
|
---|
247 | static int Y = ~0;
|
---|
248 |
|
---|
249 | int pixx = (int)(x/fPix); // [pix]
|
---|
250 | int pixy = (int)(y/fPix); // [pix]
|
---|
251 |
|
---|
252 | if (X==pixx && Y==pixy)
|
---|
253 | return;
|
---|
254 |
|
---|
255 | X = pixx;
|
---|
256 | Y = pixy;
|
---|
257 |
|
---|
258 | fLin1->SetX1(x-5.);
|
---|
259 | fLin1->SetX2(x+5.);
|
---|
260 |
|
---|
261 | fLin2->SetX1(x-5.);
|
---|
262 | fLin2->SetX2(x+5.);
|
---|
263 |
|
---|
264 | fLin1->SetY1(y-5.);
|
---|
265 | fLin1->SetY2(y+5.);
|
---|
266 |
|
---|
267 | fLin2->SetY1(y+5.);
|
---|
268 | fLin2->SetY2(y-5.);
|
---|
269 |
|
---|
270 | SetModified();
|
---|
271 | }
|
---|
272 |
|
---|
273 | void MGAccuracy::Update(ZdAz &pos, ZdAz &dev)
|
---|
274 | {
|
---|
275 | // Estimate the distance in az direction on the camera plane
|
---|
276 | const double daz = MAstro::GetDevAbs(pos.Zd(), 0, dev.Az())*3600.;
|
---|
277 |
|
---|
278 | //cout << "--> T: " << dev.Az()*60 << " " << dev.Zd()*60. << endl;
|
---|
279 |
|
---|
280 | UpdateCross(TMath::Sign(daz, dev.Az()), dev.Zd()*3600.);
|
---|
281 | Update(pos.Zd(), dev.Zd(), dev.Az());
|
---|
282 |
|
---|
283 | UpdateCanvas();
|
---|
284 | }
|
---|