1 | #ifdef HAVE_NOVA
|
---|
2 | #include "Prediction.h"
|
---|
3 | #endif
|
---|
4 |
|
---|
5 | #ifdef HAVE_SQL
|
---|
6 | #include "Database.h"
|
---|
7 | #endif
|
---|
8 |
|
---|
9 | #include <sys/stat.h> //for file stats
|
---|
10 | #include <sys/statvfs.h> //for file statvfs
|
---|
11 |
|
---|
12 | #include <numeric> // std::accumulate
|
---|
13 |
|
---|
14 | #include "Dim.h"
|
---|
15 | #include "Event.h"
|
---|
16 | #include "Shell.h"
|
---|
17 | #include "StateMachineDim.h"
|
---|
18 | #include "Connection.h"
|
---|
19 | #include "Configuration.h"
|
---|
20 | #include "Console.h"
|
---|
21 | #include "DimWriteStatistics.h"
|
---|
22 | #include "PixelMap.h"
|
---|
23 |
|
---|
24 | #include "tools.h"
|
---|
25 |
|
---|
26 | #include "LocalControl.h"
|
---|
27 |
|
---|
28 | #include "HeadersFAD.h"
|
---|
29 | #include "HeadersBIAS.h"
|
---|
30 | #include "HeadersFTM.h"
|
---|
31 | #include "HeadersFSC.h"
|
---|
32 | #include "HeadersGPS.h"
|
---|
33 | #include "HeadersSQM.h"
|
---|
34 | #include "HeadersMCP.h"
|
---|
35 | #include "HeadersLid.h"
|
---|
36 | #include "HeadersDrive.h"
|
---|
37 | #include "HeadersPower.h"
|
---|
38 | #include "HeadersPFmini.h"
|
---|
39 | #include "HeadersBiasTemp.h"
|
---|
40 | #include "HeadersAgilent.h"
|
---|
41 | #include "HeadersFeedback.h"
|
---|
42 | #include "HeadersRateScan.h"
|
---|
43 | #include "HeadersRateControl.h"
|
---|
44 | #include "HeadersTNGWeather.h"
|
---|
45 | #include "HeadersMagicLidar.h"
|
---|
46 | #include "HeadersMagicWeather.h"
|
---|
47 | #include "HeadersTemperature.h"
|
---|
48 | #include "HeadersRainSensor.h"
|
---|
49 |
|
---|
50 | #include <boost/filesystem.hpp>
|
---|
51 |
|
---|
52 | using namespace std;
|
---|
53 |
|
---|
54 | // ------------------------------------------------------------------------
|
---|
55 |
|
---|
56 | #include "DimDescriptionService.h"
|
---|
57 | #include "DimState.h"
|
---|
58 |
|
---|
59 | // ------------------------------------------------------------------------
|
---|
60 | /*
|
---|
61 | template<class T>
|
---|
62 | class buffer : public deque<T>
|
---|
63 | {
|
---|
64 | int32_t max_size;
|
---|
65 |
|
---|
66 | public:
|
---|
67 | buffer(int32_t max=-1) : max_size(max) { }
|
---|
68 | const T &operator=(const T &t) const { push_back(t); if (max_size>0 && deque<T>::size()>max_size) deque<T>::pop_front(); }
|
---|
69 | operator T() const { return deque<T>::size()>0 ? deque<T>::back() : T(); }
|
---|
70 | bool valid() const { return deque<T>::size()>0; }
|
---|
71 | };
|
---|
72 | */
|
---|
73 |
|
---|
74 | // ------------------------------------------------------------------------
|
---|
75 |
|
---|
76 | namespace HTML
|
---|
77 | {
|
---|
78 | const static string kWhite = "#ffffff";
|
---|
79 | const static string kYellow = "#fffff0";
|
---|
80 | const static string kRed = "#fff8f0";
|
---|
81 | const static string kGreen = "#f0fff0";
|
---|
82 | const static string kBlue = "#f0f0ff";
|
---|
83 | };
|
---|
84 |
|
---|
85 | // ========================================================================
|
---|
86 | // ========================================================================
|
---|
87 | // ========================================================================
|
---|
88 |
|
---|
89 | class Sun
|
---|
90 | {
|
---|
91 | public:
|
---|
92 | Time time;
|
---|
93 |
|
---|
94 | // This is always the time of the next...
|
---|
95 | Time fSunRise00;
|
---|
96 | Time fSunRise06;
|
---|
97 | Time fSunRise12;
|
---|
98 | Time fSunRise18;
|
---|
99 |
|
---|
100 | Time fSunSet00;
|
---|
101 | Time fSunSet06;
|
---|
102 | Time fSunSet12;
|
---|
103 | Time fSunSet18;
|
---|
104 |
|
---|
105 | int state;
|
---|
106 | string description;
|
---|
107 | string color;
|
---|
108 |
|
---|
109 | bool isday;
|
---|
110 | bool visible;
|
---|
111 |
|
---|
112 | #ifdef HAVE_LIBNOVA
|
---|
113 | Nova::RstTime Rst(double jd, double hrz=LN_SOLAR_STANDART_HORIZON)
|
---|
114 | {
|
---|
115 | Nova::RstTime rs = Nova::GetSolarRst(jd-0.5, hrz);
|
---|
116 | if (jd>rs.rise || jd>rs.set)
|
---|
117 | {
|
---|
118 | const Nova::RstTime rs2 = Nova::GetSolarRst(jd+0.5, hrz);
|
---|
119 | if (jd>rs.rise)
|
---|
120 | rs.rise = rs2.rise;
|
---|
121 | if (jd>rs.set)
|
---|
122 | rs.set = rs2.set;
|
---|
123 | }
|
---|
124 | return rs;
|
---|
125 | }
|
---|
126 | #endif
|
---|
127 |
|
---|
128 | public:
|
---|
129 | Sun() : time(Time::none)
|
---|
130 | {
|
---|
131 | }
|
---|
132 |
|
---|
133 | // Could be done more efficient: Only recalcuate if
|
---|
134 | // the current time exceeds at least on of the stored times
|
---|
135 | Sun(const Time &t) : time(t)
|
---|
136 | {
|
---|
137 | #ifdef HAVE_LIBNOVA
|
---|
138 | // get Julian day from local time
|
---|
139 | const double JD = time.JD();
|
---|
140 |
|
---|
141 | // >0deg : day
|
---|
142 | // -6deg - 0deg : civil twilight
|
---|
143 | // -12deg - -6deg : nautical twilight
|
---|
144 | // -18deg - -12deg : astronomical twilight
|
---|
145 | // <-18deg : night
|
---|
146 |
|
---|
147 | const Nova::RstTime sun00 = Rst(JD);
|
---|
148 | const Nova::RstTime sun06 = Rst(JD, -6);
|
---|
149 | const Nova::RstTime sun12 = Rst(JD, -12);
|
---|
150 | const Nova::RstTime sun18 = Rst(JD, -18);
|
---|
151 |
|
---|
152 | fSunRise00 = sun00.rise;
|
---|
153 | fSunRise06 = sun06.rise;
|
---|
154 | fSunRise12 = sun12.rise;
|
---|
155 | fSunRise18 = sun18.rise;
|
---|
156 |
|
---|
157 | fSunSet00 = sun00.set;
|
---|
158 | fSunSet06 = sun06.set;
|
---|
159 | fSunSet12 = sun12.set;
|
---|
160 | fSunSet18 = sun18.set;
|
---|
161 |
|
---|
162 | array<double,8> arr =
|
---|
163 | {{
|
---|
164 | sun00.set,
|
---|
165 | sun06.set,
|
---|
166 | sun12.set,
|
---|
167 | sun18.set,
|
---|
168 | sun18.rise,
|
---|
169 | sun12.rise,
|
---|
170 | sun06.rise,
|
---|
171 | sun00.rise,
|
---|
172 | }};
|
---|
173 |
|
---|
174 |
|
---|
175 | state = std::min_element(arr.begin(), arr.end())-arr.begin();
|
---|
176 |
|
---|
177 | string name[] =
|
---|
178 | {
|
---|
179 | "day time",
|
---|
180 | "civil twilight",
|
---|
181 | "nautical twilight",
|
---|
182 | "astron. twilight",
|
---|
183 | "dark time",
|
---|
184 | "astron. twilight",
|
---|
185 | "nautical twilight",
|
---|
186 | "civil twilight"
|
---|
187 | };
|
---|
188 |
|
---|
189 | description = name[state];
|
---|
190 |
|
---|
191 | const string txt = fSunRise18<fSunSet18 ?
|
---|
192 | time.MinutesTo(fSunRise18)+"↑" :
|
---|
193 | time.MinutesTo(fSunSet18)+"↓";
|
---|
194 |
|
---|
195 | description += " ["+txt+"]";
|
---|
196 |
|
---|
197 | isday = state==0;
|
---|
198 |
|
---|
199 | switch (state)
|
---|
200 | {
|
---|
201 | case 0: color = HTML::kRed; break;
|
---|
202 | case 1: case 2: color = HTML::kYellow; break;
|
---|
203 | case 3: case 4: case 5: color = HTML::kGreen; break;
|
---|
204 | case 6: case 7: color = HTML::kYellow; break;
|
---|
205 | }
|
---|
206 |
|
---|
207 | visible = state==0;
|
---|
208 |
|
---|
209 | /*
|
---|
210 | // Warning: return code of 1 means circumpolar and is not checked!
|
---|
211 | Nova::RstTime sun_day = Nova::GetSolarRst(JD-0.5);
|
---|
212 | Nova::RstTime sun_civil = Nova::GetSolarRst(JD-0.5, -6);
|
---|
213 | Nova::RstTime sun_astronomical = Nova::GetSolarRst(JD-0.5, -12);
|
---|
214 | Nova::RstTime sun_dark = Nova::GetSolarRst(JD-0.5, -18);
|
---|
215 |
|
---|
216 | fSetDayTime = Time(sun_day.set);
|
---|
217 | fSetCivil = Time(sun_civil.set);
|
---|
218 | fSetAstronomical = Time(sun_astronomical.set);
|
---|
219 | fSetDarkTime = Time(sun_dark.set);
|
---|
220 |
|
---|
221 | fRiseDayTime = Time(sun_day.rise);
|
---|
222 | fRiseCivil = Time(sun_civil.rise);
|
---|
223 | fRiseAstronomical = Time(sun_astronomical.rise);
|
---|
224 | fRiseDarkTime = Time(sun_dark.rise);
|
---|
225 |
|
---|
226 | const bool is_day = JD>sun_day.rise;
|
---|
227 | const bool is_night = JD>sun_dark.set;
|
---|
228 |
|
---|
229 | sun_day = Nova::GetSolarRst(JD+0.5);
|
---|
230 | sun_civil = Nova::GetSolarRst(JD+0.5, -6);
|
---|
231 | sun_astronomical = Nova::GetSolarRst(JD+0.5, -12);
|
---|
232 | sun_dark = Nova::GetSolarRst(JD+0.5, -18);
|
---|
233 |
|
---|
234 | if (is_day)
|
---|
235 | {
|
---|
236 | fRiseDayTime = Time(sun_day.rise);
|
---|
237 | fRiseCivil = Time(sun_civil.rise);
|
---|
238 | fRiseAstronomical = Time(sun_astronomical.rise);
|
---|
239 | fRiseDarkTime = Time(sun_dark.rise);
|
---|
240 | }
|
---|
241 |
|
---|
242 | if (is_night)
|
---|
243 | {
|
---|
244 | fSetDayTime = Time(sun_day.set);
|
---|
245 | fSetCivil = Time(sun_civil.set);
|
---|
246 | fSetAstronomical = Time(sun_astronomical.set);
|
---|
247 | fSetDarkTime = Time(sun_dark.set);
|
---|
248 | }
|
---|
249 |
|
---|
250 | // case 0: midnight to sun-rise | !is_day && !is_night | rise/set | -> isday=0
|
---|
251 | // case 1: sun-rise to sun-set | is_day && !is_night | set /rise | -> isday=1
|
---|
252 | // case 2: sun-set to midnight | is_day && is_night | rise/set | -> isday=0
|
---|
253 |
|
---|
254 | isday = is_day^is_night;
|
---|
255 |
|
---|
256 | Time fRiseDayTime; // 0: Start of day time (=end of civil twilight)
|
---|
257 | Time fRiseCivil; // -6: End of nautical twilight
|
---|
258 | Time fRiseAstronomical; // -12: End of astron. twilight
|
---|
259 | Time fRiseDarkTime; // -18: End of dark time
|
---|
260 |
|
---|
261 | Time fSetDayTime; // 0: End of day time (=start of civil twilight)
|
---|
262 | Time fSetCivil; // -6: Start of nautical twilight
|
---|
263 | Time fSetAstronomical; // -12: Start of astron. twilight
|
---|
264 | Time fSetDarkTime; // -18: Start of dark time
|
---|
265 |
|
---|
266 | state = isday ? 4 : 0; // 0 [-> Day time ]
|
---|
267 | if (time>fSetDayTime) state++; // 1 [-> Civil twilight]
|
---|
268 | if (time>fSetCivil) state++; // 2 [-> Naut. twilight]
|
---|
269 | if (time>fSetAstronomical) state++; // 3 [-> Astro. twilight]
|
---|
270 | if (time>fSetDarkTime) state++; // 4 [-> Dark time ]
|
---|
271 |
|
---|
272 | if (time>fRiseDarkTime) state++; // 5 [-> Astro. twilight]
|
---|
273 | if (time>fRiseAstronomical) state++; // 6 [-> Naut. twilight]
|
---|
274 | if (time>fRiseCivil) state++; // 7 [-> Civil twilight]
|
---|
275 | if (time>fRiseDayTime) state++; // 8 [-> Day time ]
|
---|
276 |
|
---|
277 | string name[] =
|
---|
278 | {
|
---|
279 | "dark time", // 0
|
---|
280 | "astron. twilight", // 1
|
---|
281 | "civil twilight", // 2
|
---|
282 | "sunrise", // 3
|
---|
283 | "day time", // 4
|
---|
284 | "sunset", // 5
|
---|
285 | "civil twilight", // 6
|
---|
286 | "astron. twilight", // 7
|
---|
287 | "dark time" // 8
|
---|
288 | };
|
---|
289 |
|
---|
290 | description = name[state];
|
---|
291 |
|
---|
292 | const string arr = isday ?
|
---|
293 | fSetDarkTime.MinutesTo(time)+"↓" :
|
---|
294 | fRiseDarkTime.MinutesTo(time)+"↑";
|
---|
295 |
|
---|
296 | description += " ["+arr+"]";
|
---|
297 |
|
---|
298 | switch (state)
|
---|
299 | {
|
---|
300 | case 0: case 1: color = HTML::kGreen; break;
|
---|
301 | case 2: case 3: color = HTML::kYellow; break;
|
---|
302 | case 4: color = HTML::kRed; break;
|
---|
303 | case 5: case 6: color = HTML::kYellow; break;
|
---|
304 | case 7: case 8: color = HTML::kGreen; break;
|
---|
305 | }
|
---|
306 |
|
---|
307 | visible = state>=3 && state<=5;
|
---|
308 | */
|
---|
309 | #endif
|
---|
310 | }
|
---|
311 | };
|
---|
312 |
|
---|
313 | class Moon
|
---|
314 | {
|
---|
315 | public:
|
---|
316 | Time time;
|
---|
317 |
|
---|
318 | double ra;
|
---|
319 | double dec;
|
---|
320 |
|
---|
321 | double zd;
|
---|
322 | double az;
|
---|
323 |
|
---|
324 | double disk;
|
---|
325 |
|
---|
326 | bool visible;
|
---|
327 |
|
---|
328 | Time fRise;
|
---|
329 | Time fTransit;
|
---|
330 | Time fSet;
|
---|
331 |
|
---|
332 | string description;
|
---|
333 | string color;
|
---|
334 |
|
---|
335 | int state;
|
---|
336 |
|
---|
337 | Moon() : time(Time::none)
|
---|
338 | {
|
---|
339 | }
|
---|
340 |
|
---|
341 | // Could be done more efficient: Only recalcuate if
|
---|
342 | // the current time exceeds at least on of the stored times
|
---|
343 | Moon(const Time &t) : time(t)
|
---|
344 | {
|
---|
345 | #ifdef HAVE_LIBNOVA
|
---|
346 | const double JD = time.JD();
|
---|
347 |
|
---|
348 | Nova::RstTime moon = Nova::GetLunarRst(JD-0.5);
|
---|
349 |
|
---|
350 | fRise = Time(moon.rise);
|
---|
351 | fTransit = Time(moon.transit);
|
---|
352 | fSet = Time(moon.set);
|
---|
353 |
|
---|
354 | //visible =
|
---|
355 | // ((JD>moon.rise && JD<moon.set ) && moon.rise<moon.set) ||
|
---|
356 | // ((JD<moon.set || JD>moon.rise) && moon.rise>moon.set);
|
---|
357 |
|
---|
358 | const bool is_up = JD>moon.rise;
|
---|
359 | const bool is_sinking = JD>moon.transit;
|
---|
360 | const bool is_dn = JD>moon.set;
|
---|
361 |
|
---|
362 | moon = Nova::GetLunarRst(JD+0.5);
|
---|
363 | if (is_up)
|
---|
364 | fRise = Time(moon.rise);
|
---|
365 | if (is_sinking)
|
---|
366 | fTransit = Time(moon.transit);
|
---|
367 | if (is_dn)
|
---|
368 | fSet = Time(moon.set);
|
---|
369 |
|
---|
370 | const Nova::EquPosn pos = Nova::GetLunarEquCoords(JD);
|
---|
371 | const Nova::ZdAzPosn hrz = Nova::GetHrzFromEqu(pos, JD);
|
---|
372 |
|
---|
373 | az = hrz.az;
|
---|
374 | zd = hrz.zd;
|
---|
375 |
|
---|
376 | ra = pos.ra/15;
|
---|
377 | dec = pos.dec;
|
---|
378 |
|
---|
379 | disk = Nova::GetLunarDisk(JD)*100;
|
---|
380 | state = 0;
|
---|
381 | if (fRise <fTransit && fRise <fSet) state = 0; // not visible
|
---|
382 | if (fTransit<fSet && fTransit<fRise) state = 1; // before culm
|
---|
383 | if (fSet <fRise && fSet <fTransit) state = 2; // after culm
|
---|
384 |
|
---|
385 | visible = state!=0;
|
---|
386 |
|
---|
387 | // 0: not visible
|
---|
388 | // 1: visible before cul
|
---|
389 | // 2: visible after cul
|
---|
390 |
|
---|
391 | if (!visible || disk<25)
|
---|
392 | color = HTML::kGreen;
|
---|
393 | else
|
---|
394 | color = disk>75 ? HTML::kRed : HTML::kYellow;
|
---|
395 |
|
---|
396 | const string arr = fSet<fRise ?
|
---|
397 | fSet.MinutesTo(time) +"↓" :
|
---|
398 | fRise.MinutesTo(time)+"↑";
|
---|
399 |
|
---|
400 | ostringstream out;
|
---|
401 | out << setprecision(2);
|
---|
402 | out << (visible?"visible ":"") << (disk<0.1?0:disk) << "% [" << arr << "]";
|
---|
403 |
|
---|
404 | description = out.str();
|
---|
405 | #endif
|
---|
406 | }
|
---|
407 |
|
---|
408 | double Angle(double r, double d) const
|
---|
409 | {
|
---|
410 | const double theta0 = M_PI/2-d*M_PI/180;
|
---|
411 | const double phi0 = r*M_PI/12;
|
---|
412 |
|
---|
413 | const double theta1 = M_PI/2-dec*M_PI/180;
|
---|
414 | const double phi1 = ra*M_PI/12;
|
---|
415 |
|
---|
416 | const double x0 = sin(theta0) * cos(phi0);
|
---|
417 | const double y0 = sin(theta0) * sin(phi0);
|
---|
418 | const double z0 = cos(theta0);
|
---|
419 |
|
---|
420 | const double x1 = sin(theta1) * cos(phi1);
|
---|
421 | const double y1 = sin(theta1) * sin(phi1);
|
---|
422 | const double z1 = cos(theta1);
|
---|
423 |
|
---|
424 | double arg = x0*x1 + y0*y1 + z0*z1;
|
---|
425 | if(arg > 1.0) arg = 1.0;
|
---|
426 | if(arg < -1.0) arg = -1.0;
|
---|
427 |
|
---|
428 | return acos(arg) * 180/M_PI;
|
---|
429 | }
|
---|
430 |
|
---|
431 | static string Color(double angle)
|
---|
432 | {
|
---|
433 | if (angle<10 || angle>150)
|
---|
434 | return HTML::kRed;
|
---|
435 | if (angle<20 || angle>140)
|
---|
436 | return HTML::kYellow;
|
---|
437 | return HTML::kGreen;
|
---|
438 | }
|
---|
439 | };
|
---|
440 |
|
---|
441 | // ========================================================================
|
---|
442 | // ========================================================================
|
---|
443 | // ========================================================================
|
---|
444 |
|
---|
445 | class StateMachineSmartFACT : public StateMachineDim
|
---|
446 | {
|
---|
447 | public:
|
---|
448 | static bool fIsServer;
|
---|
449 |
|
---|
450 | private:
|
---|
451 | enum states_t
|
---|
452 | {
|
---|
453 | kStateDimNetworkNA = 1,
|
---|
454 | kStateRunning,
|
---|
455 | };
|
---|
456 |
|
---|
457 | // ------------------------- History classes -----------------------
|
---|
458 |
|
---|
459 | struct EventElement
|
---|
460 | {
|
---|
461 | Time time;
|
---|
462 | string msg;
|
---|
463 |
|
---|
464 | EventElement(const Time &t, const string &s) : time(t), msg(s) { }
|
---|
465 | };
|
---|
466 |
|
---|
467 | class EventHist : public list<EventElement>
|
---|
468 | {
|
---|
469 | const boost::posix_time::time_duration deltat; //boost::posix_time::pos_infin
|
---|
470 | const uint64_t max;
|
---|
471 |
|
---|
472 | public:
|
---|
473 | EventHist(const boost::posix_time::time_duration &dt=boost::posix_time::hours(12), uint64_t mx=UINT64_MAX) : deltat(dt), max(mx) { }
|
---|
474 |
|
---|
475 | void add(const string &s, const Time &t=Time())
|
---|
476 | {
|
---|
477 | while (!empty() && (front().time+deltat<t || size()>max))
|
---|
478 | pop_front();
|
---|
479 |
|
---|
480 | emplace_back(t, s);
|
---|
481 | }
|
---|
482 |
|
---|
483 | void clean()
|
---|
484 | {
|
---|
485 | for (auto it=begin(); it!=end();)
|
---|
486 | if (!it->time)
|
---|
487 | {
|
---|
488 | const auto is = it++;
|
---|
489 | erase(is);
|
---|
490 | }
|
---|
491 | }
|
---|
492 |
|
---|
493 | string get() const
|
---|
494 | {
|
---|
495 | ostringstream out;
|
---|
496 |
|
---|
497 | string last = "";
|
---|
498 | for (auto it=begin(); it!=end(); it++)
|
---|
499 | {
|
---|
500 | const string tm = it->time.GetAsStr("%H:%M:%S ");
|
---|
501 | out << (tm!=last?tm:"--:--:-- ") << it->msg << "<br/>";
|
---|
502 | last = tm;
|
---|
503 | }
|
---|
504 |
|
---|
505 | return out.str();
|
---|
506 | }
|
---|
507 | string rget() const
|
---|
508 | {
|
---|
509 | ostringstream out;
|
---|
510 |
|
---|
511 | for (auto it=rbegin(); it!=rend(); it++)
|
---|
512 | out << it->time.GetAsStr("%H:%M:%S ") << it->msg << "<br/>";
|
---|
513 |
|
---|
514 | return out.str();
|
---|
515 | }
|
---|
516 | };
|
---|
517 |
|
---|
518 | // ------------------------- Internal variables -----------------------
|
---|
519 |
|
---|
520 | const Time fRunTime;
|
---|
521 |
|
---|
522 | PixelMap fPixelMap;
|
---|
523 |
|
---|
524 | string fDatabase;
|
---|
525 |
|
---|
526 | Time fLastUpdate;
|
---|
527 | Time fLastAstroCalc;
|
---|
528 |
|
---|
529 | string fPath;
|
---|
530 |
|
---|
531 | // ----------------------------- Data storage -------------------------
|
---|
532 |
|
---|
533 | EventHist fControlMessageHist;
|
---|
534 | EventHist fControlAlarmHist;
|
---|
535 |
|
---|
536 | int32_t fMcpConfigurationState; // For consistency
|
---|
537 | int64_t fMcpConfigurationMaxTime;
|
---|
538 | int64_t fMcpConfigurationMaxEvents;
|
---|
539 | string fMcpConfigurationName;
|
---|
540 | Time fMcpConfigurationRunStart;
|
---|
541 | EventHist fMcpConfigurationHist;
|
---|
542 | bool fLastRunFinishedWithZeroEvents;
|
---|
543 |
|
---|
544 | enum weather_t { kWeatherBegin=0, kTemp = kWeatherBegin, kDew, kHum, kPress, kWind, kGusts, kDir, kWeatherEnd = kDir+1 };
|
---|
545 | deque<float> fMagicWeatherHist[kWeatherEnd];
|
---|
546 |
|
---|
547 | deque<float> fTngWeatherDustHist;
|
---|
548 | Time fTngWeatherDustTime;
|
---|
549 |
|
---|
550 | deque<float> fRainSensorDataHist;
|
---|
551 | Time fRainSensorDataTime;
|
---|
552 |
|
---|
553 | vector<float> fBiasControlVoltageVec;
|
---|
554 |
|
---|
555 | float fBiasControlPowerTot;
|
---|
556 | float fBiasControlVoltageMed;
|
---|
557 | float fBiasControlCurrentMed;
|
---|
558 | float fBiasControlCurrentMax;
|
---|
559 |
|
---|
560 | deque<float> fBiasControlCurrentHist;
|
---|
561 | deque<float> fFscControlTemperatureHist;
|
---|
562 |
|
---|
563 | float fFscControlHumidityAvg;
|
---|
564 |
|
---|
565 | deque<float> fPfMiniHumidityHist;
|
---|
566 | deque<float> fPfMiniTemperatureHist;
|
---|
567 |
|
---|
568 | deque<float> fTemperatureControlHist;
|
---|
569 |
|
---|
570 | float fDriveControlPointingZd;
|
---|
571 | string fDriveControlPointingAz;
|
---|
572 | string fDriveControlSourceName;
|
---|
573 | float fDriveControlMoonDist;
|
---|
574 |
|
---|
575 | deque<float> fDriveControlTrackingDevHist;
|
---|
576 |
|
---|
577 | int64_t fFadControlNumEvents;
|
---|
578 | int64_t fFadControlStartRun;
|
---|
579 | int32_t fFadControlDrsStep;
|
---|
580 | vector<uint32_t> fFadControlDrsRuns;
|
---|
581 |
|
---|
582 | deque<float> fFtmControlTriggerRateHist;
|
---|
583 | int32_t fFtmControlTriggerRateTooLow;
|
---|
584 | int fFtmControlState;
|
---|
585 |
|
---|
586 | float fFtmPatchThresholdMed;
|
---|
587 | float fFtmBoardThresholdMed;
|
---|
588 |
|
---|
589 | bool fFtmControlFtuOk;
|
---|
590 |
|
---|
591 | deque<float> fRateControlThreshold;
|
---|
592 |
|
---|
593 | uint64_t fRateScanDataId;
|
---|
594 | uint8_t fRateScanBoard;
|
---|
595 | deque<float> fRateScanDataHist[41];
|
---|
596 |
|
---|
597 | set<string> fErrorList;
|
---|
598 | EventHist fErrorHist;
|
---|
599 | EventHist fChatHist;
|
---|
600 |
|
---|
601 | uint64_t fFreeSpace;
|
---|
602 |
|
---|
603 | Sun fSun;
|
---|
604 | Moon fMoon;
|
---|
605 |
|
---|
606 | // --------------------------- File header ----------------------------
|
---|
607 |
|
---|
608 | Time fAudioTime;
|
---|
609 | string fAudioName;
|
---|
610 |
|
---|
611 | string Header(const Time &d)
|
---|
612 | {
|
---|
613 | ostringstream msg;
|
---|
614 | msg << d.JavaDate() << '\t' << fAudioTime.JavaDate() << '\t' << fAudioName;
|
---|
615 | return msg.str();
|
---|
616 | }
|
---|
617 |
|
---|
618 | string Header(const EventImp &d)
|
---|
619 | {
|
---|
620 | return Header(d.GetTime());
|
---|
621 | }
|
---|
622 |
|
---|
623 | void SetAudio(const string &name)
|
---|
624 | {
|
---|
625 | fAudioName = name;
|
---|
626 | fAudioTime = Time();
|
---|
627 | }
|
---|
628 |
|
---|
629 | // ------------- Initialize variables before the Dim stuff ------------
|
---|
630 |
|
---|
631 | DimVersion fDimDNS;
|
---|
632 | DimControl fDimControl;
|
---|
633 | DimDescribedState fDimMcp;
|
---|
634 | DimDescribedState fDimDataLogger;
|
---|
635 | DimDescribedState fDimDriveControl;
|
---|
636 | DimDescribedState fDimTimeCheck;
|
---|
637 | DimDescribedState fDimMagicWeather;
|
---|
638 | DimDescribedState fDimMagicLidar;
|
---|
639 | DimDescribedState fDimTngWeather;
|
---|
640 | DimDescribedState fDimTemperature;
|
---|
641 | DimDescribedState fDimRainSensor;
|
---|
642 | DimDescribedState fDimFeedback;
|
---|
643 | DimDescribedState fDimBiasControl;
|
---|
644 | DimDescribedState fDimFtmControl;
|
---|
645 | DimDescribedState fDimFadControl;
|
---|
646 | DimDescribedState fDimFscControl;
|
---|
647 | DimDescribedState fDimPfMiniControl;
|
---|
648 | DimDescribedState fDimBiasTemp;
|
---|
649 | DimDescribedState fDimGpsControl;
|
---|
650 | DimDescribedState fDimSqmControl;
|
---|
651 | DimDescribedState fDimAgilentControl24;
|
---|
652 | DimDescribedState fDimAgilentControl50;
|
---|
653 | DimDescribedState fDimAgilentControl80;
|
---|
654 | DimDescribedState fDimPwrControl;
|
---|
655 | DimDescribedState fDimLidControl;
|
---|
656 | DimDescribedState fDimRateControl;
|
---|
657 | DimDescribedState fDimRateScan;
|
---|
658 | DimDescribedState fDimChat;
|
---|
659 | DimDescribedState fDimSkypeClient;
|
---|
660 |
|
---|
661 | // -------------------------------------------------------------------
|
---|
662 |
|
---|
663 | string GetDir(const double angle)
|
---|
664 | {
|
---|
665 | static const char *dir[] =
|
---|
666 | {
|
---|
667 | "N", "NNE", "NE", "ENE",
|
---|
668 | "E", "ESE", "SE", "SSE",
|
---|
669 | "S", "SSW", "SW", "WSW",
|
---|
670 | "W", "WNW", "NW", "NNW"
|
---|
671 | };
|
---|
672 |
|
---|
673 | const uint16_t idx = uint16_t(floor(angle/22.5+16.5))%16;
|
---|
674 | return dir[idx];
|
---|
675 | }
|
---|
676 |
|
---|
677 | // -------------------------------------------------------------------
|
---|
678 |
|
---|
679 | bool CheckDataSize(const EventImp &d, const char *name, size_t size, bool min=false)
|
---|
680 | {
|
---|
681 | if (d.GetSize()==0)
|
---|
682 | return false;
|
---|
683 |
|
---|
684 | if ((!min && d.GetSize()==size) || (min && d.GetSize()>size))
|
---|
685 | return true;
|
---|
686 |
|
---|
687 | ostringstream msg;
|
---|
688 | msg << name << " - Received service has " << d.GetSize() << " bytes, but expected ";
|
---|
689 | if (min)
|
---|
690 | msg << "more than ";
|
---|
691 | msg << size << ".";
|
---|
692 | Warn(msg);
|
---|
693 | return false;
|
---|
694 | }
|
---|
695 |
|
---|
696 | // -------------------------------------------------------------------
|
---|
697 |
|
---|
698 | template<class T>
|
---|
699 | void WriteBinaryVec(const Time &tm, const string &fname, const vector<T> &vec, double scale, double offset=0, const string &title="", const string &col="")
|
---|
700 | {
|
---|
701 | if (vec.empty())
|
---|
702 | return;
|
---|
703 |
|
---|
704 | ostringstream out;
|
---|
705 | out << tm.JavaDate() << '\n';
|
---|
706 | out << offset << '\n';
|
---|
707 | out << offset+scale << '\n';
|
---|
708 | out << setprecision(3);
|
---|
709 | if (!title.empty())
|
---|
710 | out << title << '\x7f';
|
---|
711 | else
|
---|
712 | {
|
---|
713 | const Statistics stat(vec[0]);
|
---|
714 | out << stat.min << '\n';
|
---|
715 | out << stat.med << '\n';
|
---|
716 | out << stat.max << '\x7f';
|
---|
717 | }
|
---|
718 | if (!col.empty())
|
---|
719 | out << col;
|
---|
720 | for (auto it=vec.cbegin(); it!=vec.cend(); it++)
|
---|
721 | {
|
---|
722 | // The valid range is from 1 to 127
|
---|
723 | // \0 is used to seperate different curves
|
---|
724 | vector<uint8_t> val(it->size());
|
---|
725 | for (uint64_t i=0; i<it->size(); i++)
|
---|
726 | {
|
---|
727 | float range = nearbyint(126*(double(it->at(i))-offset)/scale); // [-2V; 2V]
|
---|
728 | if (range>126)
|
---|
729 | range=126;
|
---|
730 | if (range<0)
|
---|
731 | range=0;
|
---|
732 | val[i] = (uint8_t)range;
|
---|
733 | }
|
---|
734 |
|
---|
735 | const char *ptr = reinterpret_cast<char*>(val.data());
|
---|
736 | out.write(ptr, val.size()*sizeof(uint8_t));
|
---|
737 | out << '\x7f';
|
---|
738 | }
|
---|
739 |
|
---|
740 | ofstream(fPath+"/"+fname+".bin") << out.str();
|
---|
741 | }
|
---|
742 | /*
|
---|
743 | template<class T>
|
---|
744 | void WriteBinaryVec(const EventImp &d, const string &fname, const vector<T> &vec, double scale, double offset=0, const string &title="")
|
---|
745 | {
|
---|
746 | WriteBinaryVec(d.GetTime(), fname, vec, scale, offset, title);
|
---|
747 | }
|
---|
748 |
|
---|
749 | template<class T>
|
---|
750 | void WriteBinary(const Time &tm, const string &fname, const T &t, double scale, double offset=0)
|
---|
751 | {
|
---|
752 | WriteBinaryVec(tm, fname, vector<T>(&t, &t+1), scale, offset);
|
---|
753 | }
|
---|
754 |
|
---|
755 | template<class T>
|
---|
756 | void WriteBinary(const EventImp &d, const string &fname, const T &t, double scale, double offset=0)
|
---|
757 | {
|
---|
758 | WriteBinaryVec(d.GetTime(), fname, vector<T>(&t, &t+1), scale, offset);
|
---|
759 | }*/
|
---|
760 |
|
---|
761 | template<class T>
|
---|
762 | void WriteHist(const EventImp &d, const string &fname, const T &t, double scale, double offset=0)
|
---|
763 | {
|
---|
764 | WriteBinaryVec(d.GetTime(), fname, vector<T>(&t, &t+1), scale, offset, "", "000");
|
---|
765 | }
|
---|
766 |
|
---|
767 | template<class T>
|
---|
768 | void WriteCam(const EventImp &d, const string &fname, const T &t, double scale, double offset=0)
|
---|
769 | {
|
---|
770 | WriteBinaryVec(d.GetTime(), fname, vector<T>(&t, &t+1), scale, offset, "", "");
|
---|
771 | }
|
---|
772 |
|
---|
773 |
|
---|
774 | // -------------------------------------------------------------------
|
---|
775 |
|
---|
776 | struct Statistics
|
---|
777 | {
|
---|
778 | float min;
|
---|
779 | float max;
|
---|
780 | float med;
|
---|
781 | float avg;
|
---|
782 | //float rms;
|
---|
783 |
|
---|
784 | template<class T>
|
---|
785 | Statistics(const T &t, size_t offset_min=0, size_t offset_max=0)
|
---|
786 | : min(0), max(0), med(0), avg(0)
|
---|
787 | {
|
---|
788 | if (t.empty())
|
---|
789 | return;
|
---|
790 |
|
---|
791 | T copy(t);
|
---|
792 | sort(copy.begin(), copy.end());
|
---|
793 |
|
---|
794 | if (offset_min>t.size())
|
---|
795 | offset_min = 0;
|
---|
796 | if (offset_max>t.size())
|
---|
797 | offset_max = 0;
|
---|
798 |
|
---|
799 | min = copy[offset_min];
|
---|
800 | max = copy[copy.size()-1-offset_max];
|
---|
801 | avg = accumulate (t.begin(), t.end(), 0.)/t.size();
|
---|
802 |
|
---|
803 | const size_t p = copy.size()/2;
|
---|
804 | med = copy.size()%2 ? copy[p] : (copy[p-1]+copy[p])/2.;
|
---|
805 | }
|
---|
806 | };
|
---|
807 |
|
---|
808 | void HandleControlMessageImp(const EventImp &d)
|
---|
809 | {
|
---|
810 | if (d.GetSize()==0)
|
---|
811 | return;
|
---|
812 |
|
---|
813 | fControlMessageHist.add(d.GetText(), d.GetTime());
|
---|
814 |
|
---|
815 | ostringstream out;
|
---|
816 | out << setprecision(3);
|
---|
817 | out << Header(d) << '\n';
|
---|
818 | out << HTML::kWhite << '\t';
|
---|
819 | out << "<->" << fControlMessageHist.get() << "</->";
|
---|
820 | out << '\n';
|
---|
821 |
|
---|
822 | ofstream(fPath+"/scriptlog.data") << out.str();
|
---|
823 | }
|
---|
824 |
|
---|
825 | int HandleDimControlMessage(const EventImp &d)
|
---|
826 | {
|
---|
827 | if (d.GetSize()==0)
|
---|
828 | return GetCurrentState();
|
---|
829 |
|
---|
830 | if (d.GetQoS()==MessageImp::kAlarm)
|
---|
831 | {
|
---|
832 | if (d.GetSize()<2)
|
---|
833 | for (auto it=fControlAlarmHist.begin(); it!=fControlAlarmHist.end(); it++)
|
---|
834 | it->time = Time(Time::none);
|
---|
835 | else
|
---|
836 | fControlAlarmHist.add(d.GetText(), d.GetTime());
|
---|
837 | }
|
---|
838 |
|
---|
839 | if (d.GetQoS()==MessageImp::kComment && d.GetSize()>1)
|
---|
840 | HandleControlMessageImp(d);
|
---|
841 |
|
---|
842 | return GetCurrentState();
|
---|
843 | }
|
---|
844 |
|
---|
845 | int HandleControlStateChange(const EventImp &d)
|
---|
846 | {
|
---|
847 | if (d.GetSize()==0)
|
---|
848 | return StateMachineImp::kSM_KeepState;
|
---|
849 |
|
---|
850 | if (fDimControl.scriptdepth>0)
|
---|
851 | return StateMachineImp::kSM_KeepState;
|
---|
852 |
|
---|
853 | if (d.GetQoS()>=2)
|
---|
854 | return StateMachineImp::kSM_KeepState;
|
---|
855 |
|
---|
856 | #if BOOST_VERSION < 104600
|
---|
857 | const string file = boost::filesystem::path(fDimControl.file).filename();
|
---|
858 | #else
|
---|
859 | const string file = boost::filesystem::path(fDimControl.file).filename().string();
|
---|
860 | #endif
|
---|
861 |
|
---|
862 | // [0] DimControl::kIdle
|
---|
863 | // [1] DimControl::kLoading
|
---|
864 | // [2] DimControl::kCompiling
|
---|
865 | // [3] DimControl::kRunning
|
---|
866 | if (d.GetQoS()==1)
|
---|
867 | {
|
---|
868 | fControlMessageHist.clear();
|
---|
869 | HandleControlMessageImp(Event(d, "========================================", 41));
|
---|
870 | }
|
---|
871 |
|
---|
872 | HandleControlMessageImp(Event(d, ("----- "+fDimControl.shortmsg+" -----").data(), fDimControl.shortmsg.length()+13));
|
---|
873 | if (!file.empty() && d.GetQoS()<2)
|
---|
874 | HandleControlMessageImp(Event(d, file.data(), file.length()+1));
|
---|
875 |
|
---|
876 | // Note that this will also "ding" just after program startup
|
---|
877 | // if the dimctrl is still in state -3
|
---|
878 | if (d.GetQoS()==0)
|
---|
879 | {
|
---|
880 | HandleControlMessageImp(Event(d, "========================================", 41));
|
---|
881 | if (fDimControl.last.second!=DimState::kOffline)
|
---|
882 | SetAudio("ding");
|
---|
883 | }
|
---|
884 |
|
---|
885 | return StateMachineImp::kSM_KeepState;
|
---|
886 | }
|
---|
887 |
|
---|
888 | void AddMcpConfigurationHist(const EventImp &d, const string &msg)
|
---|
889 | {
|
---|
890 | fMcpConfigurationHist.add(msg, d.GetTime());
|
---|
891 |
|
---|
892 | ostringstream out;
|
---|
893 | out << d.GetJavaDate() << '\n';
|
---|
894 | out << HTML::kWhite << '\t';
|
---|
895 | out << "<->" << fMcpConfigurationHist.rget() << "</->";
|
---|
896 | out << '\n';
|
---|
897 |
|
---|
898 | ofstream(fPath+"/observations.data") << out.str();
|
---|
899 | }
|
---|
900 |
|
---|
901 | int HandleFscControlStateChange(const EventImp &d)
|
---|
902 | {
|
---|
903 | const int32_t &last = fDimFscControl.last.second;
|
---|
904 | const int32_t &state = fDimFscControl.state();
|
---|
905 |
|
---|
906 | if (last==DimState::kOffline || state==DimState::kOffline)
|
---|
907 | return StateMachineImp::kSM_KeepState;
|
---|
908 |
|
---|
909 | if (last<FSC::State::kConnected && state==FSC::State::kConnected)
|
---|
910 | {
|
---|
911 | AddMcpConfigurationHist(d, "<B>FSC swiched on</B>");
|
---|
912 | //SetAudio("startup");
|
---|
913 | }
|
---|
914 |
|
---|
915 | if (last==FSC::State::kConnected && state<FSC::State::kConnected)
|
---|
916 | {
|
---|
917 | AddMcpConfigurationHist(d, "<B>FSC swiched off</B>");
|
---|
918 | //SetAudio("shutdown");
|
---|
919 | }
|
---|
920 |
|
---|
921 | return StateMachineImp::kSM_KeepState;
|
---|
922 | }
|
---|
923 |
|
---|
924 | int HandleMcpConfiguration(const EventImp &d)
|
---|
925 | {
|
---|
926 | if (!CheckDataSize(d, "Mcp:Configuration", 16, true))
|
---|
927 | {
|
---|
928 | fMcpConfigurationState = DimState::kOffline;
|
---|
929 | fMcpConfigurationMaxTime = 0;
|
---|
930 | fMcpConfigurationMaxEvents = 0;
|
---|
931 | fMcpConfigurationName = "";
|
---|
932 | fMcpConfigurationRunStart = Time(Time::none);
|
---|
933 | return GetCurrentState();
|
---|
934 | }
|
---|
935 |
|
---|
936 | // If a run ends...
|
---|
937 | if (fMcpConfigurationState==MCP::State::kTakingData && d.GetQoS()==MCP::State::kIdle)
|
---|
938 | {
|
---|
939 | // ...and no script is running just play a simple 'tick'
|
---|
940 | // ...and a script is running just play a simple 'tick'
|
---|
941 | if (/*fDimControl.state()<-2 &&*/ fDimControl.scriptdepth==0)
|
---|
942 | SetAudio("dong");
|
---|
943 | else
|
---|
944 | SetAudio("losticks");
|
---|
945 |
|
---|
946 | fLastRunFinishedWithZeroEvents = fFadControlNumEvents==0;
|
---|
947 |
|
---|
948 | ostringstream out;
|
---|
949 | out << "<#darkred>" << d.Ptr<char>(16);
|
---|
950 | if (!fDriveControlSourceName.empty())
|
---|
951 | out << " [" << fDriveControlSourceName << ']';
|
---|
952 | out << " (N=" << fFadControlNumEvents << ')';
|
---|
953 | out << "</#>";
|
---|
954 |
|
---|
955 | AddMcpConfigurationHist(d, out.str());
|
---|
956 | }
|
---|
957 |
|
---|
958 | if (d.GetQoS()==MCP::State::kTakingData)
|
---|
959 | {
|
---|
960 | fMcpConfigurationRunStart = Time();
|
---|
961 | SetAudio("losticks");
|
---|
962 |
|
---|
963 | ostringstream out;
|
---|
964 | out << "<#darkgreen>" << fMcpConfigurationName;
|
---|
965 | if (!fDriveControlSourceName.empty())
|
---|
966 | out << " [" << fDriveControlSourceName << ']';
|
---|
967 | if (fFadControlStartRun>0)
|
---|
968 | out << " (Run " << fFadControlStartRun << ')';
|
---|
969 | out << "</#>";
|
---|
970 |
|
---|
971 | AddMcpConfigurationHist(d, out.str());
|
---|
972 | }
|
---|
973 |
|
---|
974 | fMcpConfigurationState = d.GetQoS();
|
---|
975 | fMcpConfigurationMaxTime = d.Get<uint64_t>();
|
---|
976 | fMcpConfigurationMaxEvents = d.Get<uint64_t>(8);
|
---|
977 | fMcpConfigurationName = d.Ptr<char>(16);
|
---|
978 |
|
---|
979 | return GetCurrentState();
|
---|
980 | }
|
---|
981 |
|
---|
982 | void WriteWeather(const EventImp &d, const string &name, int i, float min, float max)
|
---|
983 | {
|
---|
984 | const Statistics stat(fMagicWeatherHist[i]);
|
---|
985 |
|
---|
986 | ostringstream out;
|
---|
987 | out << setprecision(3);
|
---|
988 | out << d.GetJavaDate() << '\n';
|
---|
989 |
|
---|
990 | out << HTML::kWhite << '\t' << fMagicWeatherHist[i].back() << '\n';
|
---|
991 | out << HTML::kWhite << '\t' << stat.min << '\n';
|
---|
992 | out << HTML::kWhite << '\t' << stat.avg << '\n';
|
---|
993 | out << HTML::kWhite << '\t' << stat.max << '\n';
|
---|
994 |
|
---|
995 | ofstream(fPath+"/"+name+".data") << out.str();
|
---|
996 |
|
---|
997 | WriteHist(d, "hist-magicweather-"+name, fMagicWeatherHist[i], max-min, min);
|
---|
998 | }
|
---|
999 |
|
---|
1000 | int HandleMagicWeatherData(const EventImp &d)
|
---|
1001 | {
|
---|
1002 | if (!CheckDataSize(d, "MagicWeather:Data", 7*4+2))
|
---|
1003 | return GetCurrentState();
|
---|
1004 |
|
---|
1005 | // Store a history of the last 300 entries
|
---|
1006 | for (int i=kWeatherBegin; i<kWeatherEnd; i++)
|
---|
1007 | {
|
---|
1008 | fMagicWeatherHist[i].push_back(d.Ptr<float>(2)[i]);
|
---|
1009 | if (fMagicWeatherHist[i].size()>300)
|
---|
1010 | fMagicWeatherHist[i].pop_front();
|
---|
1011 | }
|
---|
1012 |
|
---|
1013 | ostringstream out;
|
---|
1014 | out << d.GetJavaDate() << '\n';
|
---|
1015 | if (fSun.time.IsValid() && fMoon.time.IsValid())
|
---|
1016 | {
|
---|
1017 | out << fSun.color << '\t' << fSun.description << '\n';
|
---|
1018 | out << setprecision(2);
|
---|
1019 | out << (fSun.isday?HTML::kWhite:fMoon.color) << '\t' << fMoon.description << '\n';
|
---|
1020 | }
|
---|
1021 | else
|
---|
1022 | out << "\n\n";
|
---|
1023 | out << setprecision(3);
|
---|
1024 | for (int i=0; i<6; i++)
|
---|
1025 | out << HTML::kWhite << '\t' << fMagicWeatherHist[i].back() << '\n';
|
---|
1026 | out << HTML::kWhite << '\t' << GetDir(fMagicWeatherHist[kDir].back()) << '\n';
|
---|
1027 | out << HTML::kWhite << '\t';
|
---|
1028 | if (!fTngWeatherDustHist.empty())
|
---|
1029 | out << fTngWeatherDustHist.back() << '\t' << fTngWeatherDustTime.GetAsStr("%H:%M") << '\n';
|
---|
1030 | else
|
---|
1031 | out << "\t\n";
|
---|
1032 |
|
---|
1033 | ofstream(fPath+"/weather.data") << out.str();
|
---|
1034 |
|
---|
1035 | WriteWeather(d, "temp", kTemp, -5, 35);
|
---|
1036 | WriteWeather(d, "dew", kDew, -5, 35);
|
---|
1037 | WriteWeather(d, "hum", kHum, 0, 100);
|
---|
1038 | WriteWeather(d, "wind", kWind, 0, 100);
|
---|
1039 | WriteWeather(d, "gusts", kGusts, 0, 100);
|
---|
1040 | WriteWeather(d, "press", kPress, 700, 1000);
|
---|
1041 |
|
---|
1042 | return GetCurrentState();
|
---|
1043 | }
|
---|
1044 |
|
---|
1045 | int HandleTngWeatherDust(const EventImp &d)
|
---|
1046 | {
|
---|
1047 | if (!CheckDataSize(d, "TngWeather:Dust", 4))
|
---|
1048 | return GetCurrentState();
|
---|
1049 |
|
---|
1050 | fTngWeatherDustTime = d.GetTime();
|
---|
1051 |
|
---|
1052 | fTngWeatherDustHist.push_back(d.GetFloat());
|
---|
1053 | if (fTngWeatherDustHist.size()>300)
|
---|
1054 | fTngWeatherDustHist.pop_front();
|
---|
1055 |
|
---|
1056 | const Statistics stat(fTngWeatherDustHist);
|
---|
1057 |
|
---|
1058 | const double scale = stat.max>0 ? pow(10, ceil(log10(stat.max))) : 0;
|
---|
1059 |
|
---|
1060 | WriteHist(d, "hist-tng-dust", fTngWeatherDustHist, scale);
|
---|
1061 |
|
---|
1062 | ostringstream out;
|
---|
1063 | out << d.GetJavaDate() << '\n';
|
---|
1064 |
|
---|
1065 | ofstream(fPath+"/tngdust.data") << out.str();
|
---|
1066 |
|
---|
1067 | return GetCurrentState();
|
---|
1068 | }
|
---|
1069 |
|
---|
1070 | int HandleRainSensorData(const EventImp &d)
|
---|
1071 | {
|
---|
1072 | if (!CheckDataSize(d, "RainSensor:Data", 12)) // F:1;X:1
|
---|
1073 | return GetCurrentState();
|
---|
1074 |
|
---|
1075 | fRainSensorDataTime = d.GetTime();
|
---|
1076 |
|
---|
1077 | fRainSensorDataHist.push_back(d.GetFloat());
|
---|
1078 | if (fRainSensorDataHist.size()>300)
|
---|
1079 | fRainSensorDataHist.pop_front();
|
---|
1080 |
|
---|
1081 | const Statistics stat(fRainSensorDataHist);
|
---|
1082 |
|
---|
1083 | const double scale = stat.max>0 ? pow(10, ceil(log10(stat.max))) : 0;
|
---|
1084 |
|
---|
1085 | WriteHist(d, "hist-rain-sensor-data", fRainSensorDataHist, scale);
|
---|
1086 |
|
---|
1087 | ostringstream out;
|
---|
1088 | out << d.GetJavaDate() << '\n';
|
---|
1089 |
|
---|
1090 | ofstream(fPath+"/rainsensor.data") << out.str();
|
---|
1091 |
|
---|
1092 | return GetCurrentState();
|
---|
1093 | }
|
---|
1094 |
|
---|
1095 | int HandleDriveControlStateChange(const EventImp &d)
|
---|
1096 | {
|
---|
1097 | const int32_t &last = fDimFscControl.last.second;
|
---|
1098 | const int32_t &state = fDimFscControl.state();
|
---|
1099 |
|
---|
1100 | if (last==DimState::kOffline || state==DimState::kOffline)
|
---|
1101 | return StateMachineImp::kSM_KeepState;
|
---|
1102 |
|
---|
1103 | if (last<Drive::State::kInitialized && state>=Drive::State::kInitialized)
|
---|
1104 | AddMcpConfigurationHist(d, "Drive ready");
|
---|
1105 |
|
---|
1106 | if (last>=Drive::State::kInitialized && state<Drive::State::kInitialized)
|
---|
1107 | AddMcpConfigurationHist(d, "Drive not ready");
|
---|
1108 |
|
---|
1109 | return StateMachineImp::kSM_KeepState;
|
---|
1110 | }
|
---|
1111 |
|
---|
1112 | int HandleDrivePointing(const EventImp &d)
|
---|
1113 | {
|
---|
1114 | if (!CheckDataSize(d, "DriveControl:Pointing", 16))
|
---|
1115 | return GetCurrentState();
|
---|
1116 |
|
---|
1117 | fDriveControlPointingZd = d.Get<double>();
|
---|
1118 |
|
---|
1119 | const double az = d.Get<double>(8);
|
---|
1120 |
|
---|
1121 | fDriveControlPointingAz = GetDir(az);
|
---|
1122 |
|
---|
1123 | ostringstream out;
|
---|
1124 | out << d.GetJavaDate() << '\n';
|
---|
1125 |
|
---|
1126 | out << setprecision(0) << fixed;
|
---|
1127 | out << HTML::kWhite << '\t' << az << '\t' << fDriveControlPointingAz << '\n';
|
---|
1128 | out << HTML::kWhite << '\t' << fDriveControlPointingZd << '\n';
|
---|
1129 |
|
---|
1130 | ofstream(fPath+"/pointing.data") << out.str();
|
---|
1131 |
|
---|
1132 | return GetCurrentState();
|
---|
1133 | }
|
---|
1134 |
|
---|
1135 | int HandleDriveTracking(const EventImp &d)
|
---|
1136 | {
|
---|
1137 | if (!CheckDataSize(d, "DriveControl:Tracking", 96))
|
---|
1138 | return GetCurrentState();
|
---|
1139 |
|
---|
1140 |
|
---|
1141 |
|
---|
1142 | const double Ra = d.Get<double>(0*8);
|
---|
1143 | const double Dec = d.Get<double>(1*8);
|
---|
1144 | const double Zd = d.Get<double>(6*8);
|
---|
1145 | const double Az = d.Get<double>(7*8);
|
---|
1146 |
|
---|
1147 | const double dev = d.Get<double>(11*8);
|
---|
1148 |
|
---|
1149 | fDriveControlTrackingDevHist.push_back(dev);
|
---|
1150 | if (fDriveControlTrackingDevHist.size()>300)
|
---|
1151 | fDriveControlTrackingDevHist.pop_front();
|
---|
1152 |
|
---|
1153 | WriteHist(d, "hist-control-deviation", fDriveControlTrackingDevHist, 120);
|
---|
1154 |
|
---|
1155 | ostringstream out;
|
---|
1156 | out << d.GetJavaDate() << '\n';
|
---|
1157 |
|
---|
1158 | out << HTML::kWhite << '\t' << fDriveControlSourceName << '\n';
|
---|
1159 | out << setprecision(5);
|
---|
1160 | out << HTML::kWhite << '\t' << Ra << '\n';
|
---|
1161 | out << HTML::kWhite << '\t' << Dec << '\n';
|
---|
1162 | out << setprecision(3);
|
---|
1163 | out << HTML::kWhite << '\t' << Zd << '\n';
|
---|
1164 | out << HTML::kWhite << '\t' << Az << '\n';
|
---|
1165 | out << HTML::kWhite << '\t' << dev << '\n';
|
---|
1166 |
|
---|
1167 | fDriveControlMoonDist = -1;
|
---|
1168 |
|
---|
1169 | if (fMoon.visible)
|
---|
1170 | {
|
---|
1171 | const double angle = fMoon.Angle(Ra, Dec);
|
---|
1172 | out << Moon::Color(angle) << '\t' << setprecision(3) << angle << '\n';
|
---|
1173 |
|
---|
1174 | fDriveControlMoonDist = angle;
|
---|
1175 | }
|
---|
1176 | else
|
---|
1177 | out << HTML::kWhite << "\t— \n";
|
---|
1178 |
|
---|
1179 | ofstream(fPath+"/tracking.data") << out.str();
|
---|
1180 |
|
---|
1181 | return GetCurrentState();
|
---|
1182 | }
|
---|
1183 |
|
---|
1184 | int HandleDriveSource(const EventImp &d)
|
---|
1185 | {
|
---|
1186 | if (!CheckDataSize(d, "DriveControl:Source", 5*8+31))
|
---|
1187 | return GetCurrentState();
|
---|
1188 |
|
---|
1189 | const double *ptr = d.Ptr<double>();
|
---|
1190 |
|
---|
1191 | const double ra = ptr[0]; // Ra[h]
|
---|
1192 | const double dec = ptr[1]; // Dec[deg]
|
---|
1193 | const double woff = ptr[2]; // Wobble offset [deg]
|
---|
1194 | const double wang = ptr[3]; // Wobble angle [deg]
|
---|
1195 | const double period = ptr[4]; // Wobble angle [deg]
|
---|
1196 |
|
---|
1197 | fDriveControlSourceName = d.Ptr<char>(5*8);
|
---|
1198 |
|
---|
1199 | ostringstream out;
|
---|
1200 | out << d.GetJavaDate() << '\n';
|
---|
1201 |
|
---|
1202 | out << HTML::kWhite << '\t' << fDriveControlSourceName << '\n';
|
---|
1203 | out << setprecision(5);
|
---|
1204 | out << HTML::kWhite << '\t' << ra << '\n';
|
---|
1205 | out << HTML::kWhite << '\t' << dec << '\n';
|
---|
1206 | out << setprecision(3);
|
---|
1207 | out << HTML::kWhite << '\t' << woff << '\n';
|
---|
1208 | out << HTML::kWhite << '\t' << wang << '\n';
|
---|
1209 | out << HTML::kWhite << '\t' << period << '\n';
|
---|
1210 |
|
---|
1211 | ofstream(fPath+"/source.data") << out.str();
|
---|
1212 |
|
---|
1213 | return GetCurrentState();
|
---|
1214 | }
|
---|
1215 |
|
---|
1216 | int HandleFeedbackCalibratedCurrents(const EventImp &d)
|
---|
1217 | {
|
---|
1218 | if (!CheckDataSize(d, "Feedback:CalibratedCurrents", (416+1+1+1+1+1+416+1+1)*sizeof(float)+sizeof(uint32_t)))
|
---|
1219 | return GetCurrentState();
|
---|
1220 |
|
---|
1221 | const float *ptr = d.Ptr<float>();
|
---|
1222 |
|
---|
1223 | double power_tot = 0;
|
---|
1224 | double power_apd = 0;
|
---|
1225 |
|
---|
1226 | if (fBiasControlVoltageVec.size()>0)
|
---|
1227 | {
|
---|
1228 | // Calibrate the data (subtract offset)
|
---|
1229 | for (int i=0; i<320; i++)
|
---|
1230 | {
|
---|
1231 | // Exclude crazy pixels
|
---|
1232 | if (i==66 || i==191 || i==193)
|
---|
1233 | continue;
|
---|
1234 |
|
---|
1235 | // Group index (0 or 1) of the of the pixel (4 or 5 pixel patch)
|
---|
1236 | const int N = fPixelMap.hv(i).count();
|
---|
1237 |
|
---|
1238 | // Serial resistor of the individual G-APDs
|
---|
1239 | double R5 = 3900/N;
|
---|
1240 |
|
---|
1241 | // This is also valid for the patches with wrong resistors,
|
---|
1242 | // because Iapd is a factor f larger but R a factor f smaller
|
---|
1243 | double Iapd = ptr[i] * 1e-6; // [A]
|
---|
1244 | double Iout = Iapd*N; // [A]
|
---|
1245 |
|
---|
1246 | double UdrpCam = 1000 *Iout; // Voltage seen by everything in Camera
|
---|
1247 | double UdrpApd = (R5+2000)*Iout; // Voltage seen by G-APD
|
---|
1248 |
|
---|
1249 | const double pwrCam = Iout * (fBiasControlVoltageVec[i]-UdrpCam);
|
---|
1250 | const double pwrApd = Iout * (fBiasControlVoltageVec[i]-UdrpApd);
|
---|
1251 |
|
---|
1252 | // Total power participated in the camera at the G-APD
|
---|
1253 | // and the serial resistors (total voltage minus voltage
|
---|
1254 | // drop at resistors in bias crate)
|
---|
1255 | power_tot += pwrCam;
|
---|
1256 |
|
---|
1257 | // Power consumption per G-APD
|
---|
1258 | power_apd += pwrApd;
|
---|
1259 | }
|
---|
1260 | }
|
---|
1261 |
|
---|
1262 | // Divide by number of summed channels, convert to mW
|
---|
1263 | power_apd /= (320-3)*1e-3; // [mW]
|
---|
1264 |
|
---|
1265 | if (power_tot<1e-3)
|
---|
1266 | power_tot = 0;
|
---|
1267 | if (power_apd<1e-3)
|
---|
1268 | power_apd = 0;
|
---|
1269 |
|
---|
1270 | fBiasControlPowerTot = power_tot;
|
---|
1271 |
|
---|
1272 | // --------------------------------------------------------
|
---|
1273 |
|
---|
1274 | // Get the maximum of each patch
|
---|
1275 | vector<float> val(320, 0);
|
---|
1276 | for (int i=0; i<320; i++)
|
---|
1277 | {
|
---|
1278 | const int idx = (fPixelMap.hv(i).hw()/9)*2+fPixelMap.hv(i).group();
|
---|
1279 | val[idx] = ptr[i];
|
---|
1280 | }
|
---|
1281 |
|
---|
1282 | // Write the 160 patch values to a file
|
---|
1283 | WriteCam(d, "cam-biascontrol-current", val, 100);
|
---|
1284 |
|
---|
1285 | // --------------------------------------------------------
|
---|
1286 |
|
---|
1287 | // After being displayed, exclude the patches with
|
---|
1288 | // the crazy pixels from the statsitics
|
---|
1289 |
|
---|
1290 | vector<float> cpy(ptr, ptr+320);
|
---|
1291 | cpy[66] = 0;
|
---|
1292 | cpy[191] = 0;
|
---|
1293 | cpy[193] = 0;
|
---|
1294 | const Statistics stat(cpy);
|
---|
1295 |
|
---|
1296 | // Exclude the three crazy channels
|
---|
1297 | fBiasControlCurrentMed = stat.med;
|
---|
1298 | fBiasControlCurrentMax = stat.max;
|
---|
1299 |
|
---|
1300 | // Store a history of the last 60 entries
|
---|
1301 | fBiasControlCurrentHist.push_back(fBiasControlCurrentMed);
|
---|
1302 | if (fBiasControlCurrentHist.size()>360)
|
---|
1303 | fBiasControlCurrentHist.pop_front();
|
---|
1304 |
|
---|
1305 | // write the history to a file
|
---|
1306 | WriteHist(d, "hist-biascontrol-current", fBiasControlCurrentHist, 125);
|
---|
1307 |
|
---|
1308 | // --------------------------------------------------------
|
---|
1309 |
|
---|
1310 | string col1 = HTML::kGreen;
|
---|
1311 | string col2 = HTML::kGreen;
|
---|
1312 | string col3 = HTML::kGreen;
|
---|
1313 | string col4 = HTML::kGreen;
|
---|
1314 |
|
---|
1315 | if (stat.min>90)
|
---|
1316 | col1 = HTML::kYellow;
|
---|
1317 | if (stat.min>110)
|
---|
1318 | col1 = HTML::kRed;
|
---|
1319 |
|
---|
1320 | if (stat.med>90)
|
---|
1321 | col2 = HTML::kYellow;
|
---|
1322 | if (stat.med>110)
|
---|
1323 | col2 = HTML::kRed;
|
---|
1324 |
|
---|
1325 | if (stat.avg>90)
|
---|
1326 | col3 = HTML::kYellow;
|
---|
1327 | if (stat.avg>110)
|
---|
1328 | col3 = HTML::kRed;
|
---|
1329 |
|
---|
1330 | if (stat.max>90)
|
---|
1331 | col4 = HTML::kYellow;
|
---|
1332 | if (stat.max>110)
|
---|
1333 | col4 = HTML::kRed;
|
---|
1334 |
|
---|
1335 | ostringstream out;
|
---|
1336 | out << setprecision(3);
|
---|
1337 | out << d.GetJavaDate() << '\n';
|
---|
1338 | out << HTML::kGreen << '\t' << "yes" << '\n';
|
---|
1339 | out << col1 << '\t' << stat.min << '\n';
|
---|
1340 | out << col2 << '\t' << stat.med << '\n';
|
---|
1341 | out << col3 << '\t' << stat.avg << '\n';
|
---|
1342 | out << col4 << '\t' << stat.max << '\n';
|
---|
1343 | out << HTML::kWhite << '\t' << power_tot << "W [" << power_apd << "mW]\n";
|
---|
1344 | ofstream(fPath+"/current.data") << out.str();
|
---|
1345 |
|
---|
1346 | // --------------------------------------------------------
|
---|
1347 |
|
---|
1348 | const float Unom = ptr[2*416+6];
|
---|
1349 | const float Utmp = ptr[2*416+7];
|
---|
1350 |
|
---|
1351 | vector<float> Uov(ptr+416+6, ptr+416+6+320);
|
---|
1352 |
|
---|
1353 | WriteCam(d, "cam-feedback-overvoltage", Uov, 0.2, -0.1);
|
---|
1354 |
|
---|
1355 | const Statistics stat2(Uov);
|
---|
1356 |
|
---|
1357 | out.str("");
|
---|
1358 | out << d.GetJavaDate() << '\n';
|
---|
1359 | out << setprecision(3);
|
---|
1360 | out << HTML::kWhite << '\t' << Utmp << '\n';
|
---|
1361 | out << HTML::kWhite << '\t' << Unom << '\n';
|
---|
1362 | out << HTML::kWhite << '\t' << stat2.min << '\n';
|
---|
1363 | out << HTML::kWhite << '\t' << stat2.med << '\n';
|
---|
1364 | out << HTML::kWhite << '\t' << stat2.avg << '\n';
|
---|
1365 | out << HTML::kWhite << '\t' << stat2.max << '\n';
|
---|
1366 | ofstream(fPath+"/feedback.data") << out.str();
|
---|
1367 |
|
---|
1368 | return GetCurrentState();
|
---|
1369 | }
|
---|
1370 |
|
---|
1371 | int HandleBiasCurrent(const EventImp &d)
|
---|
1372 | {
|
---|
1373 | if (fDimFeedback.state()>=Feedback::State::kCalibrated)
|
---|
1374 | return GetCurrentState();
|
---|
1375 |
|
---|
1376 | if (!CheckDataSize(d, "BiasControl:Current", 832))
|
---|
1377 | return GetCurrentState();
|
---|
1378 |
|
---|
1379 | // Convert dac counts to uA
|
---|
1380 | vector<float> v(320);
|
---|
1381 | for (int i=0; i<320; i++)
|
---|
1382 | v[i] = d.Ptr<uint16_t>()[i] * 5000./4096;
|
---|
1383 |
|
---|
1384 | fBiasControlPowerTot = 0;
|
---|
1385 |
|
---|
1386 | // Get the maximum of each patch
|
---|
1387 | vector<float> val(320, 0);
|
---|
1388 | for (int i=0; i<320; i++)
|
---|
1389 | {
|
---|
1390 | const PixelMapEntry &hv = fPixelMap.hv(i);
|
---|
1391 | if (!hv)
|
---|
1392 | continue;
|
---|
1393 |
|
---|
1394 | const int idx = (hv.hw()/9)*2+hv.group();
|
---|
1395 | val[idx] = v[i];
|
---|
1396 | }
|
---|
1397 |
|
---|
1398 | // Write the 160 patch values to a file
|
---|
1399 | WriteCam(d, "cam-biascontrol-current", val, 1000);
|
---|
1400 |
|
---|
1401 | const Statistics stat(v, 0, 3);
|
---|
1402 |
|
---|
1403 | // Exclude the three crazy channels
|
---|
1404 | fBiasControlCurrentMed = stat.med;
|
---|
1405 | fBiasControlCurrentMax = stat.max;
|
---|
1406 |
|
---|
1407 | // Store a history of the last 60 entries
|
---|
1408 | fBiasControlCurrentHist.push_back(fBiasControlCurrentMed);
|
---|
1409 | if (fBiasControlCurrentHist.size()>360)
|
---|
1410 | fBiasControlCurrentHist.pop_front();
|
---|
1411 |
|
---|
1412 | // write the history to a file
|
---|
1413 | WriteHist(d, "hist-biascontrol-current", fBiasControlCurrentHist, 1000);
|
---|
1414 |
|
---|
1415 | ostringstream out;
|
---|
1416 | out << setprecision(3);
|
---|
1417 | out << d.GetJavaDate() << '\n';
|
---|
1418 | out << HTML::kWhite<< '\t' << "no" << '\n';
|
---|
1419 | out << HTML::kWhite << '\t' << stat.min << '\n';
|
---|
1420 | out << HTML::kWhite << '\t' << stat.med << '\n';
|
---|
1421 | out << HTML::kWhite << '\t' << stat.avg << '\n';
|
---|
1422 | out << HTML::kWhite << '\t' << stat.max << '\n';
|
---|
1423 | out << HTML::kWhite << '\t' << "---\n";
|
---|
1424 | ofstream(fPath+"/current.data") << out.str();
|
---|
1425 |
|
---|
1426 | return GetCurrentState();
|
---|
1427 | }
|
---|
1428 |
|
---|
1429 | int HandleBiasVoltage(const EventImp &d)
|
---|
1430 | {
|
---|
1431 | if (!CheckDataSize(d, "BiasControl:Voltage", 1664))
|
---|
1432 | {
|
---|
1433 | fBiasControlVoltageVec.clear();
|
---|
1434 | return GetCurrentState();
|
---|
1435 | }
|
---|
1436 |
|
---|
1437 | fBiasControlVoltageVec.assign(d.Ptr<float>(), d.Ptr<float>()+320);
|
---|
1438 |
|
---|
1439 | const Statistics stat(fBiasControlVoltageVec);
|
---|
1440 |
|
---|
1441 | fBiasControlVoltageMed = stat.med;
|
---|
1442 |
|
---|
1443 | vector<float> val(320, 0);
|
---|
1444 | for (int i=0; i<320; i++)
|
---|
1445 | {
|
---|
1446 | const int idx = (fPixelMap.hv(i).hw()/9)*2+fPixelMap.hv(i).group();
|
---|
1447 | val[idx] = fBiasControlVoltageVec[i];
|
---|
1448 | }
|
---|
1449 |
|
---|
1450 | if (fDimBiasControl.state()==BIAS::State::kVoltageOn || fDimBiasControl.state()==BIAS::State::kRamping)
|
---|
1451 | WriteCam(d, "cam-biascontrol-voltage", val, 10, 65);
|
---|
1452 | else
|
---|
1453 | WriteCam(d, "cam-biascontrol-voltage", val, 75);
|
---|
1454 |
|
---|
1455 | ostringstream out;
|
---|
1456 | out << setprecision(3);
|
---|
1457 | out << d.GetJavaDate() << '\n';
|
---|
1458 | out << HTML::kWhite << '\t' << stat.min << '\n';
|
---|
1459 | out << HTML::kWhite << '\t' << stat.med << '\n';
|
---|
1460 | out << HTML::kWhite << '\t' << stat.avg << '\n';
|
---|
1461 | out << HTML::kWhite << '\t' << stat.max << '\n';
|
---|
1462 | ofstream(fPath+"/voltage.data") << out.str();
|
---|
1463 |
|
---|
1464 | return GetCurrentState();
|
---|
1465 | }
|
---|
1466 |
|
---|
1467 | int HandleFadEvents(const EventImp &d)
|
---|
1468 | {
|
---|
1469 | if (!CheckDataSize(d, "FadControl:Events", 4*4))
|
---|
1470 | {
|
---|
1471 | fFadControlNumEvents = -1;
|
---|
1472 | return GetCurrentState();
|
---|
1473 | }
|
---|
1474 |
|
---|
1475 | fFadControlNumEvents = d.Get<uint32_t>();
|
---|
1476 |
|
---|
1477 | return GetCurrentState();
|
---|
1478 | }
|
---|
1479 |
|
---|
1480 | int HandleFadStartRun(const EventImp &d)
|
---|
1481 | {
|
---|
1482 | if (!CheckDataSize(d, "FadControl:StartRun", 16))
|
---|
1483 | {
|
---|
1484 | fFadControlStartRun = -1;
|
---|
1485 | return GetCurrentState();
|
---|
1486 | }
|
---|
1487 |
|
---|
1488 | fFadControlStartRun = d.Get<int64_t>();
|
---|
1489 |
|
---|
1490 | return GetCurrentState();
|
---|
1491 | }
|
---|
1492 |
|
---|
1493 | int HandleFadDrsRuns(const EventImp &d)
|
---|
1494 | {
|
---|
1495 | if (!CheckDataSize(d, "FadControl:DrsRuns", 5*4))
|
---|
1496 | {
|
---|
1497 | fFadControlDrsStep = -1;
|
---|
1498 | return GetCurrentState();
|
---|
1499 | }
|
---|
1500 |
|
---|
1501 | const uint32_t *ptr = d.Ptr<uint32_t>();
|
---|
1502 | fFadControlDrsStep = ptr[0];
|
---|
1503 | fFadControlDrsRuns[0] = ptr[1];
|
---|
1504 | fFadControlDrsRuns[1] = ptr[2];
|
---|
1505 | fFadControlDrsRuns[2] = ptr[3];
|
---|
1506 |
|
---|
1507 | return GetCurrentState();
|
---|
1508 | }
|
---|
1509 |
|
---|
1510 | int HandleFadConnections(const EventImp &d)
|
---|
1511 | {
|
---|
1512 | if (!CheckDataSize(d, "FadControl:Connections", 41))
|
---|
1513 | {
|
---|
1514 | //fStatusEventBuilderLabel->setText("Offline");
|
---|
1515 | return GetCurrentState();
|
---|
1516 | }
|
---|
1517 |
|
---|
1518 | string rc(40, '-'); // orange/red [45]
|
---|
1519 |
|
---|
1520 | const uint8_t *ptr = d.Ptr<uint8_t>();
|
---|
1521 |
|
---|
1522 | int c[4] = { '.', '.', '.', '.' };
|
---|
1523 |
|
---|
1524 | for (int i=0; i<40; i++)
|
---|
1525 | {
|
---|
1526 | const uint8_t stat1 = ptr[i]&3;
|
---|
1527 | const uint8_t stat2 = ptr[i]>>3;
|
---|
1528 |
|
---|
1529 | if (stat1==0 && stat2==0)
|
---|
1530 | rc[i] = '.'; // gray [46]
|
---|
1531 | else
|
---|
1532 | if (stat1>=2 && stat2==8)
|
---|
1533 | rc[i] = stat1==2?'+':'*'; // green [43] : check [42]
|
---|
1534 |
|
---|
1535 | if (rc[i]<c[i/10])
|
---|
1536 | c[i/10] = rc[i];
|
---|
1537 | }
|
---|
1538 |
|
---|
1539 | string col[4];
|
---|
1540 | for (int i=0; i<4; i++)
|
---|
1541 | switch (c[i])
|
---|
1542 | {
|
---|
1543 | case '.': col[i]=HTML::kWhite; break;
|
---|
1544 | case '-': col[i]=HTML::kRed; break;
|
---|
1545 | case '+': col[i]=HTML::kYellow; break;
|
---|
1546 | case '*': col[i]=HTML::kGreen; break;
|
---|
1547 | }
|
---|
1548 |
|
---|
1549 | ostringstream out;
|
---|
1550 | out << setprecision(3);
|
---|
1551 | out << d.GetJavaDate() << '\n';
|
---|
1552 | out << col[0] << '\t' << rc.substr( 0, 10) << '\n';
|
---|
1553 | out << col[1] << '\t' << rc.substr(10, 10) << '\n';
|
---|
1554 | out << col[2] << '\t' << rc.substr(20, 10) << '\n';
|
---|
1555 | out << col[3] << '\t' << rc.substr(30, 10) << '\n';
|
---|
1556 | ofstream(fPath+"/fad.data") << out.str();
|
---|
1557 |
|
---|
1558 | return GetCurrentState();
|
---|
1559 | }
|
---|
1560 |
|
---|
1561 | /*
|
---|
1562 | int HandleFtmControlStateChange()
|
---|
1563 | {
|
---|
1564 | const int32_t &last = fDimFtmControl.last.second;
|
---|
1565 | const int32_t &state = fDimFtmControl.state();
|
---|
1566 |
|
---|
1567 | // If a new run has been started ensure that the counter
|
---|
1568 | // is reset. The reset in HandleFtmTriggerRates might
|
---|
1569 | // arrive only after the run was started.
|
---|
1570 | if (last!=FTM::State::kTriggerOn && state==MCP::State::kTriggerOn)
|
---|
1571 | fFtmControlTriggerRateTooLow = -1;
|
---|
1572 |
|
---|
1573 | return StateMachineImp::kSM_KeepState;
|
---|
1574 | }*/
|
---|
1575 |
|
---|
1576 |
|
---|
1577 | int HandleFtmTriggerRates(const EventImp &d)
|
---|
1578 | {
|
---|
1579 | if (!CheckDataSize(d, "FtmControl:TriggerRates", 24+160+640+8))
|
---|
1580 | {
|
---|
1581 | fFtmControlTriggerRateTooLow = 0;
|
---|
1582 | return GetCurrentState();
|
---|
1583 | }
|
---|
1584 |
|
---|
1585 | const FTM::DimTriggerRates &dim = d.Ref<FTM::DimTriggerRates>();
|
---|
1586 |
|
---|
1587 | // If the trigger rate is too low...
|
---|
1588 | // ... and the run was not just started (can lead to very small elapsed times)
|
---|
1589 | // ... and the trigger is switched on
|
---|
1590 | // ... and there was no state change (then the trigger was started or stopped)
|
---|
1591 | fFtmControlTriggerRateTooLow =
|
---|
1592 | dim.fTriggerRate<1 && dim.fElapsedTime>0.45 &&
|
---|
1593 | (fFtmControlState&FTM::kFtmStates)==FTM::kFtmRunning &&
|
---|
1594 | (fFtmControlState&FTM::kFtmStates)==(d.GetQoS()&FTM::kFtmStates);
|
---|
1595 |
|
---|
1596 | fFtmControlState = d.GetQoS();
|
---|
1597 |
|
---|
1598 | const float *brates = dim.fBoardRate; // Board rate
|
---|
1599 | const float *prates = dim.fPatchRate; // Patch rate
|
---|
1600 |
|
---|
1601 | // Store a history of the last 60 entries
|
---|
1602 | fFtmControlTriggerRateHist.push_back(dim.fTriggerRate);
|
---|
1603 | if (fFtmControlTriggerRateHist.size()>300)
|
---|
1604 | fFtmControlTriggerRateHist.pop_front();
|
---|
1605 |
|
---|
1606 | // FIXME: Add statistics for all kind of rates
|
---|
1607 |
|
---|
1608 | WriteHist(d, "hist-ftmcontrol-triggerrate",
|
---|
1609 | fFtmControlTriggerRateHist, 100);
|
---|
1610 | WriteCam(d, "cam-ftmcontrol-boardrates",
|
---|
1611 | vector<float>(brates, brates+40), 10);
|
---|
1612 | WriteCam(d, "cam-ftmcontrol-patchrates",
|
---|
1613 | vector<float>(prates, prates+160), 10);
|
---|
1614 |
|
---|
1615 | ostringstream out;
|
---|
1616 | out << setprecision(3);
|
---|
1617 | out << d.GetJavaDate() << '\n';
|
---|
1618 | out << HTML::kWhite << '\t' << dim.fTriggerRate << '\n';
|
---|
1619 |
|
---|
1620 | ofstream(fPath+"/trigger.data") << out.str();
|
---|
1621 |
|
---|
1622 | const Statistics bstat(vector<float>(brates, brates+ 40));
|
---|
1623 | const Statistics pstat(vector<float>(prates, prates+160));
|
---|
1624 |
|
---|
1625 | out.str("");
|
---|
1626 | out << d.GetJavaDate() << '\n';
|
---|
1627 | out << HTML::kWhite << '\t' << bstat.min << '\n';
|
---|
1628 | out << HTML::kWhite << '\t' << bstat.med << '\n';
|
---|
1629 | out << HTML::kWhite << '\t' << bstat.avg << '\n';
|
---|
1630 | out << HTML::kWhite << '\t' << bstat.max << '\n';
|
---|
1631 | ofstream(fPath+"/boardrates.data") << out.str();
|
---|
1632 |
|
---|
1633 | out.str("");
|
---|
1634 | out << d.GetJavaDate() << '\n';
|
---|
1635 | out << HTML::kWhite << '\t' << pstat.min << '\n';
|
---|
1636 | out << HTML::kWhite << '\t' << pstat.med << '\n';
|
---|
1637 | out << HTML::kWhite << '\t' << pstat.avg << '\n';
|
---|
1638 | out << HTML::kWhite << '\t' << pstat.max << '\n';
|
---|
1639 | ofstream(fPath+"/patchrates.data") << out.str();
|
---|
1640 |
|
---|
1641 | return GetCurrentState();
|
---|
1642 | }
|
---|
1643 |
|
---|
1644 | int HandleFtmStaticData(const EventImp &d)
|
---|
1645 | {
|
---|
1646 | if (!CheckDataSize(d, "FtmControl:StaticData", sizeof(FTM::DimStaticData)))
|
---|
1647 | return GetCurrentState();
|
---|
1648 |
|
---|
1649 | // If the FTM is in state Configuring, the clock conditioner
|
---|
1650 | // is always reported to be unlocked
|
---|
1651 | fFtmControlState = d.GetQoS();
|
---|
1652 |
|
---|
1653 | const FTM::DimStaticData &dat = d.Ref<FTM::DimStaticData>();
|
---|
1654 |
|
---|
1655 | vector<uint16_t> vecp(dat.fThreshold, dat.fThreshold+160);
|
---|
1656 | vector<uint16_t> vecb(dat.fMultiplicity, dat.fMultiplicity+40);
|
---|
1657 |
|
---|
1658 | WriteCam(d, "cam-ftmcontrol-thresholds-patch", vecp, 1000);
|
---|
1659 | WriteCam(d, "cam-ftmcontrol-thresholds-board", vecb, 100);
|
---|
1660 |
|
---|
1661 | const Statistics statp(vecp);
|
---|
1662 | const Statistics statb(vecb);
|
---|
1663 |
|
---|
1664 | fFtmPatchThresholdMed = statp.med;
|
---|
1665 | fFtmBoardThresholdMed = statb.med;
|
---|
1666 |
|
---|
1667 | ostringstream out;
|
---|
1668 | out << d.GetJavaDate() << '\n';
|
---|
1669 | out << HTML::kWhite << '\t' << statb.min << '\n';
|
---|
1670 | out << HTML::kWhite << '\t' << statb.med << '\n';
|
---|
1671 | out << HTML::kWhite << '\t' << statb.max << '\n';
|
---|
1672 | ofstream(fPath+"/thresholds-board.data") << out.str();
|
---|
1673 |
|
---|
1674 | out.str("");
|
---|
1675 | out << d.GetJavaDate() << '\n';
|
---|
1676 | out << HTML::kWhite << '\t' << statp.min << '\n';
|
---|
1677 | out << HTML::kWhite << '\t' << statp.med << '\n';
|
---|
1678 | out << HTML::kWhite << '\t' << statp.max << '\n';
|
---|
1679 | ofstream(fPath+"/thresholds-patch.data") << out.str();
|
---|
1680 |
|
---|
1681 | out.str("");
|
---|
1682 | out << d.GetJavaDate() << '\n';
|
---|
1683 | out << HTML::kWhite << '\t' << statb.med << '\n';
|
---|
1684 | out << HTML::kWhite << '\t' << statp.med << '\n';
|
---|
1685 | ofstream(fPath+"/thresholds.data") << out.str();
|
---|
1686 |
|
---|
1687 | out.str("");
|
---|
1688 | out << d.GetJavaDate() << '\n';
|
---|
1689 | out << HTML::kWhite << '\t' << dat.fTriggerInterval << '\n';
|
---|
1690 | out << HTML::kWhite << '\t';
|
---|
1691 | if (dat.HasPedestal())
|
---|
1692 | out << dat.fTriggerSeqPed;
|
---|
1693 | else
|
---|
1694 | out << "–";
|
---|
1695 | out << ':';
|
---|
1696 | if (dat.HasLPext())
|
---|
1697 | out << dat.fTriggerSeqLPext;
|
---|
1698 | else
|
---|
1699 | out << "–";
|
---|
1700 | out << ':';
|
---|
1701 | if (dat.HasLPint())
|
---|
1702 | out << dat.fTriggerSeqLPint;
|
---|
1703 | else
|
---|
1704 | out << "–";
|
---|
1705 | out << '\n';
|
---|
1706 |
|
---|
1707 | out << HTML::kWhite << '\t' << (dat.HasTrigger()?"on":"off") << " / " << (dat.HasExt1()?"on":"off") << " / " << (dat.HasExt2()?"on":"off") << '\n';
|
---|
1708 | out << HTML::kWhite << '\t' << (dat.HasVeto()?"on":"off") << " / " << (dat.HasClockConditioner()?"time cal":"marker") << '\n';
|
---|
1709 | out << HTML::kWhite << '\t' << dat.fMultiplicityPhysics << " / " << dat.fMultiplicityCalib << '\n';
|
---|
1710 | out << HTML::kWhite << '\t' << dat.fWindowPhysics << '\t' << dat.fWindowCalib << '\n';
|
---|
1711 | out << HTML::kWhite << '\t' << dat.fDelayTrigger << '\t' << dat.fDelayTimeMarker << '\n';
|
---|
1712 | out << HTML::kWhite << '\t' << dat.fDeadTime << '\n';
|
---|
1713 |
|
---|
1714 | int64_t vp = dat.fPrescaling[0];
|
---|
1715 | for (int i=1; i<40; i++)
|
---|
1716 | if (vp!=dat.fPrescaling[i])
|
---|
1717 | vp = -1;
|
---|
1718 |
|
---|
1719 | if (vp<0)
|
---|
1720 | out << HTML::kYellow << "\tdifferent\n";
|
---|
1721 | else
|
---|
1722 | out << HTML::kWhite << '\t' << 0.5*vp << "\n";
|
---|
1723 |
|
---|
1724 | ofstream(fPath+"/ftm.data") << out.str();
|
---|
1725 |
|
---|
1726 | // Active FTUs: IsActive(i)
|
---|
1727 | // Enabled Pix: IsEnabled(i)
|
---|
1728 |
|
---|
1729 | return GetCurrentState();
|
---|
1730 | }
|
---|
1731 |
|
---|
1732 | int HandleFtmFtuList(const EventImp &d)
|
---|
1733 | {
|
---|
1734 | if (!CheckDataSize(d, "FtmControl:FtuList", sizeof(FTM::DimFtuList)))
|
---|
1735 | return GetCurrentState();
|
---|
1736 |
|
---|
1737 | const FTM::DimFtuList &sdata = d.Ref<FTM::DimFtuList>();
|
---|
1738 |
|
---|
1739 | ostringstream out;
|
---|
1740 | out << d.GetJavaDate() << '\n';
|
---|
1741 |
|
---|
1742 | int cnt = 0;
|
---|
1743 | for (int i=0; i<4; i++)
|
---|
1744 | {
|
---|
1745 | out << HTML::kWhite << '\t';
|
---|
1746 | for (int j=0; j<10; j++)
|
---|
1747 | if (sdata.IsActive(i*10+j))
|
---|
1748 | {
|
---|
1749 | if (sdata.fPing[i*10+j]==1)
|
---|
1750 | {
|
---|
1751 | out << '*';
|
---|
1752 | cnt++;
|
---|
1753 | }
|
---|
1754 | else
|
---|
1755 | out << sdata.fPing[i*10+j];
|
---|
1756 | }
|
---|
1757 | else
|
---|
1758 | out << '-';
|
---|
1759 | out << '\n';
|
---|
1760 | }
|
---|
1761 |
|
---|
1762 | fFtmControlFtuOk = cnt==40;
|
---|
1763 |
|
---|
1764 | ofstream(fPath+"/ftu.data") << out.str();
|
---|
1765 |
|
---|
1766 | return GetCurrentState();
|
---|
1767 | }
|
---|
1768 |
|
---|
1769 | int HandleFadEventData(const EventImp &d)
|
---|
1770 | {
|
---|
1771 | if (!CheckDataSize(d, "FadControl:EventData", 23048))
|
---|
1772 | return GetCurrentState();
|
---|
1773 |
|
---|
1774 | //const uint32_t run = d.GetUint();
|
---|
1775 | //const uint32_t evt = d.GetUint(4);
|
---|
1776 |
|
---|
1777 | const float *dat = d.Ptr<float>(8+1440*sizeof(float)*2);
|
---|
1778 |
|
---|
1779 | /*
|
---|
1780 | vector<float> max(320, 0);
|
---|
1781 | for (int i=0; i<1440; i++)
|
---|
1782 | {
|
---|
1783 | if (i%9==8)
|
---|
1784 | continue;
|
---|
1785 |
|
---|
1786 | const int idx = (fPixelMap.hw(i).hw()/9)*2+fPixelMap.hw(i).group();
|
---|
1787 | const double v = dat[i]/1000;
|
---|
1788 | //if (v>max[idx])
|
---|
1789 | // max[idx]=v;
|
---|
1790 |
|
---|
1791 | max[idx] += v/4;
|
---|
1792 | } */
|
---|
1793 |
|
---|
1794 | vector<float> val(1440);
|
---|
1795 | for (int i=0; i<1440; i++)
|
---|
1796 | val[i] = dat[i%9==8 ? i-2 : i]/1000;
|
---|
1797 |
|
---|
1798 | vector<float> sorted(val);
|
---|
1799 | nth_element(sorted.begin(), sorted.begin()+3, sorted.end(),
|
---|
1800 | std::greater<float>());
|
---|
1801 |
|
---|
1802 | const uint32_t trig = d.GetQoS() & FAD::EventHeader::kLPext;
|
---|
1803 |
|
---|
1804 | const float min = fFadControlDrsRuns[0]==0 ? -1 : 0;
|
---|
1805 |
|
---|
1806 | float scale = 2;
|
---|
1807 | if (trig&FAD::EventHeader::kLPext)
|
---|
1808 | scale = 1;
|
---|
1809 | if (trig&FAD::EventHeader::kPedestal)
|
---|
1810 | scale = 0.25;
|
---|
1811 | if (trig==0)
|
---|
1812 | scale = max(0.25f, sorted[3]);
|
---|
1813 |
|
---|
1814 | // assume it is drs-gain
|
---|
1815 | //if ((trig&FAD::EventHeader::kPedestal) && fFadControlDrsRuns[0]>0 && fFadControlDrsRuns[1]==0)
|
---|
1816 | // min = 0.75;
|
---|
1817 |
|
---|
1818 | WriteCam(d, "cam-fadcontrol-eventdata", val, scale, min);
|
---|
1819 |
|
---|
1820 | return GetCurrentState();
|
---|
1821 | }
|
---|
1822 |
|
---|
1823 | int HandleStats(const EventImp &d)
|
---|
1824 | {
|
---|
1825 | if (!CheckDataSize(d, "Stats", 4*8))
|
---|
1826 | {
|
---|
1827 | fFreeSpace = UINT64_MAX;
|
---|
1828 | return GetCurrentState();
|
---|
1829 | }
|
---|
1830 |
|
---|
1831 | const DimWriteStatistics::Stats &s = d.Ref<DimWriteStatistics::Stats>();
|
---|
1832 | fFreeSpace = s.freeSpace;
|
---|
1833 |
|
---|
1834 | return GetCurrentState();
|
---|
1835 | }
|
---|
1836 |
|
---|
1837 | int HandleFscTemperature(const EventImp &d)
|
---|
1838 | {
|
---|
1839 | if (!CheckDataSize(d, "FscControl:Temperature", 240))
|
---|
1840 | return GetCurrentState();
|
---|
1841 |
|
---|
1842 | const float *ptr = d.Ptr<float>(4);
|
---|
1843 |
|
---|
1844 | double avg = 0;
|
---|
1845 | double rms = 0;
|
---|
1846 | double min = 99;
|
---|
1847 | double max = -99;
|
---|
1848 |
|
---|
1849 | int num = 0;
|
---|
1850 | for (const float *t=ptr; t<ptr+31; t++)
|
---|
1851 | {
|
---|
1852 | if (*t==0)
|
---|
1853 | continue;
|
---|
1854 |
|
---|
1855 | if (*t>max)
|
---|
1856 | max = *t;
|
---|
1857 |
|
---|
1858 | if (*t<min)
|
---|
1859 | min = *t;
|
---|
1860 |
|
---|
1861 | avg += *t;
|
---|
1862 | rms += *t * *t;
|
---|
1863 |
|
---|
1864 | num++;
|
---|
1865 | }
|
---|
1866 |
|
---|
1867 | avg /= num;
|
---|
1868 | rms /= num;
|
---|
1869 | rms += avg*avg;
|
---|
1870 | rms = rms<0 ? 0 : sqrt(rms);
|
---|
1871 |
|
---|
1872 | // Clean broken reports
|
---|
1873 | static double pre_rms1 = 1.5;
|
---|
1874 | static double pre_rms2 = 0;
|
---|
1875 |
|
---|
1876 | const double cut = pre_rms1 + 0.1;
|
---|
1877 |
|
---|
1878 | const bool reject = rms>cut && pre_rms2<cut;
|
---|
1879 |
|
---|
1880 | pre_rms2 = pre_rms1;
|
---|
1881 | pre_rms1 = rms;
|
---|
1882 |
|
---|
1883 | if (reject)
|
---|
1884 | return GetCurrentState();
|
---|
1885 |
|
---|
1886 |
|
---|
1887 | if (!fMagicWeatherHist[kTemp].empty())
|
---|
1888 | {
|
---|
1889 | fFscControlTemperatureHist.push_back(avg-fMagicWeatherHist[kTemp].back());
|
---|
1890 | if (fFscControlTemperatureHist.size()>300)
|
---|
1891 | fFscControlTemperatureHist.pop_front();
|
---|
1892 | }
|
---|
1893 |
|
---|
1894 | const Statistics stat(fFscControlTemperatureHist);
|
---|
1895 |
|
---|
1896 | ostringstream out;
|
---|
1897 | out << setprecision(3);
|
---|
1898 | out << d.GetJavaDate() << '\n';
|
---|
1899 | out << HTML::kWhite << '\t' << fFscControlHumidityAvg << '\n';
|
---|
1900 | out << HTML::kWhite << '\t' << stat.min << '\n';
|
---|
1901 | out << HTML::kWhite << '\t' << stat.avg << '\n';
|
---|
1902 | out << HTML::kWhite << '\t' << stat.max << '\n';
|
---|
1903 |
|
---|
1904 | ofstream(fPath+"/fsc.data") << out.str();
|
---|
1905 |
|
---|
1906 | WriteHist(d, "hist-fsccontrol-temperature",
|
---|
1907 | fFscControlTemperatureHist, 10);
|
---|
1908 |
|
---|
1909 | out.str("");
|
---|
1910 | out << setprecision(3);
|
---|
1911 | out << d.GetJavaDate() << '\n';
|
---|
1912 | out << HTML::kWhite << '\t' << max << '\n';
|
---|
1913 | out << HTML::kWhite << '\t' << avg << '\n';
|
---|
1914 | out << HTML::kWhite << '\t' << min << '\n';
|
---|
1915 |
|
---|
1916 | ofstream(fPath+"/camtemp.data") << out.str();
|
---|
1917 |
|
---|
1918 | return GetCurrentState();
|
---|
1919 | }
|
---|
1920 |
|
---|
1921 | int HandleFscBiasTemp(const EventImp &d)
|
---|
1922 | {
|
---|
1923 | if (!CheckDataSize(d, "FscControl:BiasTemp", 323*4))
|
---|
1924 | return GetCurrentState();
|
---|
1925 |
|
---|
1926 | const float *ptr = d.Ptr<float>(4);
|
---|
1927 | const float avg = d.Get<float>(321*4);
|
---|
1928 | //const float rms = d.Get<float>(322*4);
|
---|
1929 |
|
---|
1930 | vector<double> tout(320);
|
---|
1931 | for (int i=0; i<320; i++)
|
---|
1932 | {
|
---|
1933 | const int idx = (fPixelMap.hv(i).hw()/9)*2+fPixelMap.hv(i).group();
|
---|
1934 | tout[idx] = ptr[i];
|
---|
1935 | }
|
---|
1936 |
|
---|
1937 | WriteCam(d, "cam-fsccontrol-temperature", tout, 3, avg-1.75);
|
---|
1938 |
|
---|
1939 | return GetCurrentState();
|
---|
1940 | }
|
---|
1941 |
|
---|
1942 | int HandleFscHumidity(const EventImp &d)
|
---|
1943 | {
|
---|
1944 | if (!CheckDataSize(d, "FscControl:Humidity", 5*4))
|
---|
1945 | return GetCurrentState();
|
---|
1946 |
|
---|
1947 | const float *ptr = d.Ptr<float>(4);
|
---|
1948 |
|
---|
1949 | double avg =0;
|
---|
1950 | int num = 0;
|
---|
1951 |
|
---|
1952 | for (const float *t=ptr; t<ptr+4; t++)
|
---|
1953 | if (*t>0 && *t<=100 && t!=ptr+2 /*excl broken sensor*/)
|
---|
1954 | {
|
---|
1955 | avg += *t;
|
---|
1956 | num++;
|
---|
1957 | }
|
---|
1958 |
|
---|
1959 | fFscControlHumidityAvg = num>0 ? avg/num : 0;
|
---|
1960 |
|
---|
1961 | return GetCurrentState();
|
---|
1962 | }
|
---|
1963 |
|
---|
1964 | int HandlePfMiniData(const EventImp &d)
|
---|
1965 | {
|
---|
1966 | if (!CheckDataSize(d, "PfMini:Data", sizeof(PFmini::Data)))
|
---|
1967 | return GetCurrentState();
|
---|
1968 |
|
---|
1969 | const PFmini::Data &data = d.Ref<PFmini::Data>();
|
---|
1970 |
|
---|
1971 | ostringstream out;
|
---|
1972 |
|
---|
1973 | out << fixed << setprecision(1);
|
---|
1974 | out << d.GetJavaDate() << '\n';
|
---|
1975 |
|
---|
1976 | out << HTML::kGreen << '\t' << data.temp << '\n';
|
---|
1977 | out << HTML::kGreen << '\t' << data.hum << '\n';
|
---|
1978 |
|
---|
1979 | ofstream(fPath+"/pfmini.data") << out.str();
|
---|
1980 |
|
---|
1981 | fPfMiniTemperatureHist.push_back(data.temp);
|
---|
1982 | if (fPfMiniTemperatureHist.size()>60*4) // 1h
|
---|
1983 | fPfMiniTemperatureHist.pop_front();
|
---|
1984 |
|
---|
1985 | fPfMiniHumidityHist.push_back(data.hum);
|
---|
1986 | if (fPfMiniHumidityHist.size()>60*4) // 1h
|
---|
1987 | fPfMiniHumidityHist.pop_front();
|
---|
1988 |
|
---|
1989 | WriteHist(d, "hist-pfmini-temp",
|
---|
1990 | fPfMiniTemperatureHist, 45, 0);
|
---|
1991 |
|
---|
1992 | WriteHist(d, "hist-pfmini-hum",
|
---|
1993 | fPfMiniHumidityHist, 100, 0);
|
---|
1994 |
|
---|
1995 | return GetCurrentState();
|
---|
1996 | }
|
---|
1997 |
|
---|
1998 | int HandleBiasTemp(const EventImp &d)
|
---|
1999 | {
|
---|
2000 | if (!CheckDataSize(d, "BiasTemp:Data", sizeof(BiasTemp::Data)))
|
---|
2001 | return GetCurrentState();
|
---|
2002 |
|
---|
2003 | const BiasTemp::Data &data = d.Ref<BiasTemp::Data>();
|
---|
2004 |
|
---|
2005 | ostringstream out;
|
---|
2006 |
|
---|
2007 | out << fixed << setprecision(1);
|
---|
2008 | out << d.GetJavaDate() << '\n';
|
---|
2009 |
|
---|
2010 | out << HTML::kGreen << '\t' << data.time << '\n';
|
---|
2011 | out << HTML::kGreen << '\t' << data.avg << '\n';
|
---|
2012 | out << HTML::kGreen << '\t' << data.rms << '\n';
|
---|
2013 |
|
---|
2014 | ofstream(fPath+"/biastemp.data") << out.str();
|
---|
2015 |
|
---|
2016 | return GetCurrentState();
|
---|
2017 | }
|
---|
2018 |
|
---|
2019 | int HandleGpsNema(const EventImp &d)
|
---|
2020 | {
|
---|
2021 | if (!CheckDataSize(d, "GpsControl:Nema", sizeof(GPS::NEMA)))
|
---|
2022 | return GetCurrentState();
|
---|
2023 |
|
---|
2024 | const GPS::NEMA &nema = d.Ref<GPS::NEMA>();
|
---|
2025 |
|
---|
2026 | ostringstream out;
|
---|
2027 |
|
---|
2028 | out << fixed;
|
---|
2029 | out << d.GetJavaDate() << '\n';
|
---|
2030 |
|
---|
2031 | switch (nema.qos)
|
---|
2032 | {
|
---|
2033 | case 1: out << HTML::kGreen << "\tGPS fix [1]\n"; break;
|
---|
2034 | case 2: out << HTML::kGreen << "\tDifferential fix [2]\n"; break;
|
---|
2035 | default: out << HTML::kRed << "\tinvalid [" << nema.qos << "]\n"; break;
|
---|
2036 | }
|
---|
2037 |
|
---|
2038 | out << HTML::kWhite << '\t' << nema.count << '\n';
|
---|
2039 | out << HTML::kWhite << '\t' << Time(floor(Time().Mjd())+nema.time).GetAsStr("%H:%M:%S") << '\n';
|
---|
2040 | out << HTML::kWhite << '\t' << setprecision(4) << nema.lat << '\n';
|
---|
2041 | out << HTML::kWhite << '\t' << setprecision(4) << nema.lng << '\n';
|
---|
2042 | out << HTML::kWhite << '\t' << setprecision(1) << nema.height << "\n";
|
---|
2043 | out << HTML::kWhite << '\t' << setprecision(1) << nema.hdop << "\n";
|
---|
2044 | out << HTML::kWhite << '\t' << setprecision(1) << nema.geosep << "\n";
|
---|
2045 |
|
---|
2046 | ofstream(fPath+"/gps.data") << out.str();
|
---|
2047 |
|
---|
2048 | return GetCurrentState();
|
---|
2049 | }
|
---|
2050 |
|
---|
2051 | int HandleSqmData(const EventImp &d)
|
---|
2052 | {
|
---|
2053 | if (!CheckDataSize(d, "SqmControl:Data", sizeof(SQM::Data)))
|
---|
2054 | return GetCurrentState();
|
---|
2055 |
|
---|
2056 | const SQM::Data &data = d.Ref<SQM::Data>();
|
---|
2057 |
|
---|
2058 | ostringstream out;
|
---|
2059 |
|
---|
2060 | out << fixed;
|
---|
2061 | out << d.GetJavaDate() << '\n';
|
---|
2062 | out << HTML::kWhite << '\t' << setprecision(2) << data.mag << '\n';
|
---|
2063 | out << HTML::kWhite << '\t' << data.freq << '\n';
|
---|
2064 | out << HTML::kWhite << '\t' << data.counts << '\n';
|
---|
2065 | out << HTML::kWhite << '\t' << setprecision(3) << data.period << '\n';
|
---|
2066 | out << HTML::kWhite << '\t' << setprecision(1) << data.temp << "\n";
|
---|
2067 |
|
---|
2068 | ofstream(fPath+"/sqm.data") << out.str();
|
---|
2069 |
|
---|
2070 | return GetCurrentState();
|
---|
2071 | }
|
---|
2072 |
|
---|
2073 | string GetTempColor(float t)
|
---|
2074 | {
|
---|
2075 | if (t>25 && t<30)
|
---|
2076 | return HTML::kGreen;
|
---|
2077 |
|
---|
2078 | if (t<20 || t>35)
|
---|
2079 | return HTML::kRed;
|
---|
2080 |
|
---|
2081 | return HTML::kYellow;
|
---|
2082 | }
|
---|
2083 |
|
---|
2084 | int HandleTemperatureData(const EventImp &d)
|
---|
2085 | {
|
---|
2086 | if (!CheckDataSize(d, "Temperature:Data", 3*sizeof(float)))
|
---|
2087 | return GetCurrentState();
|
---|
2088 |
|
---|
2089 | const float *temp = d.Ptr<float>();
|
---|
2090 |
|
---|
2091 | ostringstream out;
|
---|
2092 |
|
---|
2093 | out << fixed << setprecision(1);
|
---|
2094 | out << d.GetJavaDate() << '\n';
|
---|
2095 |
|
---|
2096 | out << GetTempColor(temp[1]) << '\t' << temp[1] << '\n';
|
---|
2097 | out << GetTempColor(temp[0]) << '\t' << temp[0] << '\n';
|
---|
2098 | out << GetTempColor(temp[2]) << '\t' << temp[2] << '\n';
|
---|
2099 |
|
---|
2100 | ofstream(fPath+"/temperature.data") << out.str();
|
---|
2101 |
|
---|
2102 | fTemperatureControlHist.push_back(temp[0]);
|
---|
2103 | if (fTemperatureControlHist.size()>60) // 1h
|
---|
2104 | fTemperatureControlHist.pop_front();
|
---|
2105 |
|
---|
2106 | WriteHist(d, "hist-temperaturecontrol",
|
---|
2107 | fTemperatureControlHist, 45, 0);
|
---|
2108 |
|
---|
2109 | return GetCurrentState();
|
---|
2110 | }
|
---|
2111 |
|
---|
2112 | int HandleAgilentData(const EventImp &d, const string &ext)
|
---|
2113 | {
|
---|
2114 | if (!CheckDataSize(d, ("Agilent"+ext+":Data").c_str(), 4*sizeof(float)))
|
---|
2115 | return GetCurrentState();
|
---|
2116 |
|
---|
2117 | const float *data = d.Ptr<float>();
|
---|
2118 |
|
---|
2119 | ostringstream out;
|
---|
2120 |
|
---|
2121 | out << fixed << setprecision(1);
|
---|
2122 | out << d.GetJavaDate() << '\n';
|
---|
2123 |
|
---|
2124 | out << HTML::kWhite << '\t' << data[0] << '\n';
|
---|
2125 | out << HTML::kWhite << '\t' << data[1] << '\n';
|
---|
2126 | out << HTML::kWhite << '\t' << data[2] << '\n';
|
---|
2127 | out << HTML::kWhite << '\t' << data[3] << '\n';
|
---|
2128 |
|
---|
2129 | ofstream(fPath+"/agilent"+ext+".data") << out.str();
|
---|
2130 |
|
---|
2131 | return GetCurrentState();
|
---|
2132 | }
|
---|
2133 |
|
---|
2134 | int HandleRateScanData(const EventImp &d)
|
---|
2135 | {
|
---|
2136 | if (!CheckDataSize(d, "RateScan:Data", 824))
|
---|
2137 | return GetCurrentState();
|
---|
2138 |
|
---|
2139 | const uint64_t id = d.Get<uint64_t>();
|
---|
2140 | const float *rate = d.Ptr<float>(20);
|
---|
2141 |
|
---|
2142 | if (fRateScanDataId!=id)
|
---|
2143 | {
|
---|
2144 | for (int i=0; i<41; i++)
|
---|
2145 | fRateScanDataHist[i].clear();
|
---|
2146 | fRateScanDataId = id;
|
---|
2147 | }
|
---|
2148 | fRateScanDataHist[0].push_back(log10(rate[0]));
|
---|
2149 |
|
---|
2150 | double max = 0;
|
---|
2151 | for (int i=1; i<41; i++)
|
---|
2152 | {
|
---|
2153 | fRateScanDataHist[i].push_back(log10(rate[i]));
|
---|
2154 | if (rate[i]>max)
|
---|
2155 | max = rate[i];
|
---|
2156 | }
|
---|
2157 |
|
---|
2158 | // Cycle by time!
|
---|
2159 | fRateScanBoard ++;
|
---|
2160 | fRateScanBoard %= 40;
|
---|
2161 |
|
---|
2162 | WriteHist(d, "hist-ratescan", fRateScanDataHist[0], 10, -2);
|
---|
2163 | WriteCam(d, "cam-ratescan-board", fRateScanDataHist[fRateScanBoard+1], 10, -4);
|
---|
2164 |
|
---|
2165 | ostringstream out;
|
---|
2166 | out << setprecision(3);
|
---|
2167 | out << d.GetJavaDate() << '\n';
|
---|
2168 | out << HTML::kWhite << '\t' << fFtmBoardThresholdMed << '\n';
|
---|
2169 | out << HTML::kWhite << '\t' << fFtmPatchThresholdMed << '\n';
|
---|
2170 | out << HTML::kWhite << '\t' << floor(pow(10, fRateScanDataHist[0].back())+.5) << '\n';
|
---|
2171 | out << HTML::kWhite << '\t' << floor(max+.5) << '\n';
|
---|
2172 |
|
---|
2173 | ofstream(fPath+"/ratescan.data") << out.str();
|
---|
2174 |
|
---|
2175 | out.str("");
|
---|
2176 | out << d.GetJavaDate() << '\n';
|
---|
2177 | out << HTML::kWhite << '\t' << int(fRateScanBoard) << '\n';
|
---|
2178 | out << HTML::kWhite << '\t' << pow(10, fRateScanDataHist[fRateScanBoard+1].back()) << '\n';
|
---|
2179 |
|
---|
2180 | ofstream(fPath+"/ratescan_board.data") << out.str();
|
---|
2181 |
|
---|
2182 | return GetCurrentState();
|
---|
2183 | }
|
---|
2184 |
|
---|
2185 | int HandleRateControlThreshold(const EventImp &d)
|
---|
2186 | {
|
---|
2187 | if (!CheckDataSize(d, "RateControl:Threshold", 18))
|
---|
2188 | return GetCurrentState();
|
---|
2189 |
|
---|
2190 | const uint16_t th = d.Get<uint16_t>();
|
---|
2191 |
|
---|
2192 | fRateControlThreshold.push_back(th);
|
---|
2193 | if (fRateControlThreshold.size()>300)
|
---|
2194 | fRateControlThreshold.pop_front();
|
---|
2195 |
|
---|
2196 | WriteHist(d, "hist-ratecontrol-threshold", fRateControlThreshold, 1000);
|
---|
2197 |
|
---|
2198 | return GetCurrentState();
|
---|
2199 | }
|
---|
2200 |
|
---|
2201 | int HandleChatMsg(const EventImp &d)
|
---|
2202 | {
|
---|
2203 | if (d.GetSize()==0 || d.GetQoS()!=MessageImp::kComment)
|
---|
2204 | return GetCurrentState();
|
---|
2205 |
|
---|
2206 | if (Time()<d.GetTime()+boost::posix_time::minutes(1))
|
---|
2207 | SetAudio("message");
|
---|
2208 |
|
---|
2209 | fChatHist.add(d.GetText(), d.GetTime());
|
---|
2210 |
|
---|
2211 | ostringstream out;
|
---|
2212 | out << setprecision(3);
|
---|
2213 | out << Header(d) << '\n';
|
---|
2214 | out << HTML::kWhite << '\t';
|
---|
2215 | out << "<->" << fChatHist.rget() << "</->";
|
---|
2216 | out << '\n';
|
---|
2217 |
|
---|
2218 | ofstream(fPath+"/chat.data") << out.str();
|
---|
2219 |
|
---|
2220 | return GetCurrentState();
|
---|
2221 | }
|
---|
2222 |
|
---|
2223 | // -------------------------------------------------------------------
|
---|
2224 |
|
---|
2225 | int HandleDoTest(const EventImp &d)
|
---|
2226 | {
|
---|
2227 | ostringstream out;
|
---|
2228 | out << d.GetJavaDate() << '\n';
|
---|
2229 |
|
---|
2230 | switch (d.GetQoS())
|
---|
2231 | {
|
---|
2232 | case -3: out << HTML::kWhite << "\tNot running\n"; break;
|
---|
2233 | case -2: out << HTML::kBlue << "\tLoading\n"; break;
|
---|
2234 | case -1: out << HTML::kBlue << "\tStarted\n"; break;
|
---|
2235 | default: out << HTML::kGreen << "\tRunning [" << d.GetQoS() << "]\n"; break;
|
---|
2236 | }
|
---|
2237 |
|
---|
2238 | ofstream(fPath+"/dotest.data") << out.str();
|
---|
2239 |
|
---|
2240 | return StateMachineImp::kSM_KeepState;
|
---|
2241 | }
|
---|
2242 |
|
---|
2243 | // -------------------------------------------------------------------
|
---|
2244 |
|
---|
2245 | /*
|
---|
2246 | bool CheckEventSize(size_t has, const char *name, size_t size)
|
---|
2247 | {
|
---|
2248 | if (has==size)
|
---|
2249 | return true;
|
---|
2250 |
|
---|
2251 | ostringstream msg;
|
---|
2252 | msg << name << " - Received event has " << has << " bytes, but expected " << size << ".";
|
---|
2253 | Fatal(msg);
|
---|
2254 | return false;
|
---|
2255 | }*/
|
---|
2256 |
|
---|
2257 | int Print() const
|
---|
2258 | {
|
---|
2259 | Out() << fDimDNS << endl;
|
---|
2260 | Out() << fDimMcp << endl;
|
---|
2261 | Out() << fDimControl << endl;
|
---|
2262 | Out() << fDimDataLogger << endl;
|
---|
2263 | Out() << fDimDriveControl << endl;
|
---|
2264 | Out() << fDimTimeCheck << endl;
|
---|
2265 | Out() << fDimFadControl << endl;
|
---|
2266 | Out() << fDimFtmControl << endl;
|
---|
2267 | Out() << fDimBiasControl << endl;
|
---|
2268 | Out() << fDimFeedback << endl;
|
---|
2269 | Out() << fDimRateControl << endl;
|
---|
2270 | Out() << fDimFscControl << endl;
|
---|
2271 | Out() << fDimAgilentControl24 << endl;
|
---|
2272 | Out() << fDimAgilentControl50 << endl;
|
---|
2273 | Out() << fDimAgilentControl80 << endl;
|
---|
2274 | Out() << fDimPwrControl << endl;
|
---|
2275 | Out() << fDimLidControl << endl;
|
---|
2276 | Out() << fDimMagicWeather << endl;
|
---|
2277 | Out() << fDimTngWeather << endl;
|
---|
2278 | Out() << fDimMagicLidar << endl;
|
---|
2279 | Out() << fDimTemperature << endl;
|
---|
2280 | Out() << fDimRainSensor << endl;
|
---|
2281 | Out() << fDimRateScan << endl;
|
---|
2282 | Out() << fDimChat << endl;
|
---|
2283 | Out() << fDimSkypeClient << endl;
|
---|
2284 |
|
---|
2285 | return GetCurrentState();
|
---|
2286 | }
|
---|
2287 |
|
---|
2288 | string GetStateHtml(const DimState &state, int green) const
|
---|
2289 | {
|
---|
2290 | if (!state.online())
|
---|
2291 | return HTML::kWhite+"\t—\n";
|
---|
2292 |
|
---|
2293 | if (&state==&fDimControl)
|
---|
2294 | return HTML::kGreen +'\t'+(state.state()==0?"Idle":fDimControl.shortmsg)+'\n';
|
---|
2295 |
|
---|
2296 | const State rc = state.description();
|
---|
2297 |
|
---|
2298 | // Sate not found in list, server online (-3: offline; -2: not found)
|
---|
2299 | if (rc.index==-2)
|
---|
2300 | {
|
---|
2301 | ostringstream out;
|
---|
2302 | out << HTML::kWhite << '\t' << state.state() << '\n';
|
---|
2303 | return out.str();
|
---|
2304 | }
|
---|
2305 |
|
---|
2306 | //ostringstream msg;
|
---|
2307 | //msg << HTML::kWhite << '\t' << rc.name << " [" << rc.index << "]\n";
|
---|
2308 | //return msg.str();
|
---|
2309 |
|
---|
2310 | if (rc.index<0)
|
---|
2311 | return HTML::kWhite + "\t—\n";
|
---|
2312 |
|
---|
2313 | string col = HTML::kGreen;
|
---|
2314 | if (rc.index<green)
|
---|
2315 | col = HTML::kYellow;
|
---|
2316 | if (rc.index>0xff)
|
---|
2317 | col = HTML::kRed;
|
---|
2318 |
|
---|
2319 | return col + '\t' + rc.name + '\n';
|
---|
2320 | }
|
---|
2321 |
|
---|
2322 | bool SetError(bool b, const string &err)
|
---|
2323 | {
|
---|
2324 | if (!b)
|
---|
2325 | {
|
---|
2326 | fErrorList.erase(err);
|
---|
2327 | return 0;
|
---|
2328 | }
|
---|
2329 |
|
---|
2330 | const bool isnew = fErrorList.insert(err).second;
|
---|
2331 | if (isnew)
|
---|
2332 | fErrorHist.add(err);
|
---|
2333 |
|
---|
2334 | return isnew;
|
---|
2335 | }
|
---|
2336 |
|
---|
2337 | #ifdef HAVE_NOVA
|
---|
2338 |
|
---|
2339 | //vector<pair<Nova::EquPosn, double>> fMoonCoords;
|
---|
2340 |
|
---|
2341 | vector<Nova::SolarObjects> fCoordinates;
|
---|
2342 |
|
---|
2343 | void CalcCoordinates(double jd)
|
---|
2344 | {
|
---|
2345 | jd = floor(jd);
|
---|
2346 |
|
---|
2347 | fCoordinates.clear();
|
---|
2348 | for (double h=0; h<1; h+=1./(24*12))
|
---|
2349 | fCoordinates.emplace_back(jd+h);
|
---|
2350 | }
|
---|
2351 |
|
---|
2352 | pair<vector<float>, pair<Time, float>> GetVisibility(Nova::EquPosn *src=0)
|
---|
2353 | {
|
---|
2354 | const double sunset = fSun.fSunSet12.JD()-1;
|
---|
2355 | const double sunrise = fSun.fSunRise12.JD();
|
---|
2356 |
|
---|
2357 | Nova::EquPosn moon;
|
---|
2358 | Nova::EquPosn *pos = src ? src : &moon;
|
---|
2359 |
|
---|
2360 | double max = 0;
|
---|
2361 | double maxjd = 0;
|
---|
2362 |
|
---|
2363 | int cnt = 0;
|
---|
2364 |
|
---|
2365 | vector<float> alt;
|
---|
2366 | for (auto it=fCoordinates.begin(); it!=fCoordinates.end(); it++)
|
---|
2367 | {
|
---|
2368 | if (src==0)
|
---|
2369 | moon = it->fMoonEqu;
|
---|
2370 |
|
---|
2371 | const Nova::HrzPosn hrz = Nova::GetHrzFromEqu(*pos, it->fJD);
|
---|
2372 |
|
---|
2373 | if (it->fJD>sunset && it->fJD<sunrise)
|
---|
2374 | alt.push_back(hrz.alt);
|
---|
2375 |
|
---|
2376 | if (hrz.alt>max)
|
---|
2377 | {
|
---|
2378 | max = hrz.alt;
|
---|
2379 | maxjd = it->fJD;
|
---|
2380 | }
|
---|
2381 |
|
---|
2382 | if (it->fJD>sunset && it->fJD<sunrise && hrz.alt>15)
|
---|
2383 | cnt++;
|
---|
2384 | }
|
---|
2385 |
|
---|
2386 | if (max<=15 || cnt==0)
|
---|
2387 | return make_pair(vector<float>(), make_pair(Time(), 0));
|
---|
2388 |
|
---|
2389 | return make_pair(alt, make_pair(maxjd, maxjd>sunset&&maxjd<sunrise?max:0));
|
---|
2390 | }
|
---|
2391 |
|
---|
2392 | pair<vector<float>, pair<Time, float>> GetLightCondition(const Nova::EquPosn &src_pos)
|
---|
2393 | {
|
---|
2394 | const double sunset = fSun.fSunSet12.JD()-1;
|
---|
2395 | const double sunrise = fSun.fSunRise12.JD();
|
---|
2396 |
|
---|
2397 | double max = -1;
|
---|
2398 | double maxjd = 0;
|
---|
2399 |
|
---|
2400 | int cnt = 0;
|
---|
2401 |
|
---|
2402 | vector<float> vec;
|
---|
2403 | for (auto it=fCoordinates.begin(); it!=fCoordinates.end(); it++)
|
---|
2404 | {
|
---|
2405 | double cur = -1;
|
---|
2406 |
|
---|
2407 | if (it->fJD>sunset && it->fJD<sunrise)
|
---|
2408 | {
|
---|
2409 | cur = FACT::PredictI(*it, src_pos);
|
---|
2410 | vec.push_back(cur);
|
---|
2411 | }
|
---|
2412 |
|
---|
2413 | if (cur>max)
|
---|
2414 | {
|
---|
2415 | max = cur;
|
---|
2416 | maxjd = it->fJD;
|
---|
2417 | }
|
---|
2418 |
|
---|
2419 | if (it->fJD>sunset && it->fJD<sunrise && cur>0)
|
---|
2420 | cnt++;
|
---|
2421 | }
|
---|
2422 |
|
---|
2423 | if (max<=0 || cnt==0)
|
---|
2424 | return make_pair(vector<float>(), make_pair(Time(), 0));
|
---|
2425 |
|
---|
2426 | return make_pair(vec, make_pair(maxjd, maxjd>sunset&&maxjd<sunrise?max:-1));
|
---|
2427 | }
|
---|
2428 | #endif
|
---|
2429 |
|
---|
2430 | void UpdateAstronomy()
|
---|
2431 | {
|
---|
2432 | Time now;
|
---|
2433 |
|
---|
2434 | #ifdef HAVE_NOVA
|
---|
2435 | CalcCoordinates(now.JD());
|
---|
2436 | #endif
|
---|
2437 |
|
---|
2438 | fSun = Sun (now);
|
---|
2439 | fMoon = Moon(now);
|
---|
2440 |
|
---|
2441 | vector<string> color(8, HTML::kWhite);
|
---|
2442 | color[fSun.state] = HTML::kBlue;
|
---|
2443 |
|
---|
2444 | ostringstream out;
|
---|
2445 | out << setprecision(3);
|
---|
2446 | out << now.JavaDate() << '\n';
|
---|
2447 | out << color[4] << '\t' << fSun.fSunRise18.GetAsStr("%H:%M") << '\n';
|
---|
2448 | out << color[5] << '\t' << fSun.fSunRise12.GetAsStr("%H:%M") << '\n';
|
---|
2449 | out << color[6] << '\t' << fSun.fSunRise06.GetAsStr("%H:%M") << '\n';
|
---|
2450 | out << color[7] << '\t' << fSun.fSunRise00.GetAsStr("%H:%M") << '\n';
|
---|
2451 |
|
---|
2452 | out << color[0] << '\t' << fSun.fSunSet00.GetAsStr("%H:%M") << '\n';
|
---|
2453 | out << color[1] << '\t' << fSun.fSunSet06.GetAsStr("%H:%M") << '\n';
|
---|
2454 | out << color[2] << '\t' << fSun.fSunSet12.GetAsStr("%H:%M") << '\n';
|
---|
2455 | out << color[3] << '\t' << fSun.fSunSet18.GetAsStr("%H:%M") << '\n';
|
---|
2456 |
|
---|
2457 | ofstream(fPath+"/sun.data") << out.str();
|
---|
2458 |
|
---|
2459 | color.assign(3, HTML::kWhite);
|
---|
2460 | color[fMoon.state%3] = HTML::kBlue;
|
---|
2461 |
|
---|
2462 | out.str("");
|
---|
2463 | out << now.JavaDate() << '\n';
|
---|
2464 |
|
---|
2465 | out << color[0] << '\t' << fMoon.fRise.GetAsStr("%H:%M") << '\n';
|
---|
2466 | out << color[1] << '\t' << fMoon.fTransit.GetAsStr("%H:%M") << '\n';
|
---|
2467 | out << color[2] << '\t' << fMoon.fSet.GetAsStr("%H:%M") << '\n';
|
---|
2468 |
|
---|
2469 | out << (fSun.isday?HTML::kWhite:fMoon.color) << '\t' << fMoon.description << '\n';
|
---|
2470 |
|
---|
2471 | if (!fMoon.visible)
|
---|
2472 | out << HTML::kWhite << "\t—\t\n";
|
---|
2473 | else
|
---|
2474 | {
|
---|
2475 | string col = HTML::kWhite;
|
---|
2476 | if (!fSun.isday)
|
---|
2477 | {
|
---|
2478 | col = HTML::kGreen;
|
---|
2479 | if (fMoon.zd>25)
|
---|
2480 | col = HTML::kYellow;
|
---|
2481 | if (fMoon.zd>45 && fMoon.zd<80)
|
---|
2482 | col = HTML::kRed;
|
---|
2483 | if (fMoon.zd>=80)
|
---|
2484 | col = HTML::kRed;
|
---|
2485 | }
|
---|
2486 | out << col << '\t' << fMoon.zd << '\t' << GetDir(fMoon.az) << '\n';
|
---|
2487 | }
|
---|
2488 |
|
---|
2489 | ostringstream out2, out3, out4;
|
---|
2490 | out2 << setprecision(3);
|
---|
2491 | out2 << now.JavaDate() << '\n';
|
---|
2492 | out3 << now.JavaDate() << '\n';
|
---|
2493 | out4 << now.JavaDate() << '\n';
|
---|
2494 |
|
---|
2495 | struct Entry
|
---|
2496 | {
|
---|
2497 | string name;
|
---|
2498 | float value;
|
---|
2499 | int color;
|
---|
2500 | Entry(const string &n, float v, int c) : name(n), value(v), color(c%8) { }
|
---|
2501 |
|
---|
2502 | const string &Col() const
|
---|
2503 | {
|
---|
2504 | // If this list is updatd the number count in the constructor needs
|
---|
2505 | // to be updated, too
|
---|
2506 | static const string hcol[] = { "888", "8cf", "c8f", "bbb", "8fc", "cf8", "f8c", "fc8" };
|
---|
2507 | return hcol[color];
|
---|
2508 | }
|
---|
2509 |
|
---|
2510 | vector<float> GetColor(double scale, double offset=0) const
|
---|
2511 | {
|
---|
2512 | vector<float> rc(3);
|
---|
2513 | rc[0] = double(Col()[0])*scale/126+offset;
|
---|
2514 | rc[1] = double(Col()[1])*scale/126+offset;
|
---|
2515 | rc[2] = double(Col()[2])*scale/126+offset;
|
---|
2516 | return rc;
|
---|
2517 | }
|
---|
2518 | };
|
---|
2519 |
|
---|
2520 | multimap<Time, Entry> culmination;
|
---|
2521 | multimap<Time, Entry> lightcond;
|
---|
2522 | vector<vector<float>> alt;
|
---|
2523 | vector<vector<float>> cur;
|
---|
2524 |
|
---|
2525 | #ifdef HAVE_NOVA
|
---|
2526 | int ccol = 0;
|
---|
2527 | int lcol = 0;
|
---|
2528 |
|
---|
2529 | /*const*/ pair<vector<float>, pair<Time, float>> vism = GetVisibility();
|
---|
2530 | if (!vism.first.empty())
|
---|
2531 | {
|
---|
2532 | const Entry entry("Moon", vism.second.second, ccol);
|
---|
2533 | culmination.insert(make_pair(vism.second.first, entry));
|
---|
2534 | const vector<float> col = entry.GetColor(75, 15);
|
---|
2535 | vism.first.insert(vism.first.begin(), col.begin(), col.end());
|
---|
2536 | alt.push_back(vism.first);
|
---|
2537 |
|
---|
2538 | ccol++;
|
---|
2539 | }
|
---|
2540 | #endif
|
---|
2541 |
|
---|
2542 | #ifdef HAVE_SQL
|
---|
2543 | try
|
---|
2544 | {
|
---|
2545 | const mysqlpp::StoreQueryResult res =
|
---|
2546 | Database(fDatabase).query("SELECT fSourceName, fRightAscension, fDeclination FROM Source WHERE fSourceTypeKEY=1").store();
|
---|
2547 |
|
---|
2548 | out << HTML::kWhite << '\t';
|
---|
2549 | out2 << HTML::kWhite << '\t';
|
---|
2550 | out3 << HTML::kWhite << '\t';
|
---|
2551 | out4 << HTML::kWhite << '\t';
|
---|
2552 |
|
---|
2553 | for (vector<mysqlpp::Row>::const_iterator v=res.begin(); v<res.end(); v++)
|
---|
2554 | {
|
---|
2555 | const string name = (*v)[0].c_str();
|
---|
2556 | const double ra = (*v)[1];
|
---|
2557 | const double dec = (*v)[2];
|
---|
2558 | #ifdef HAVE_NOVA
|
---|
2559 | Nova::EquPosn pos;
|
---|
2560 | pos.ra = ra*15;
|
---|
2561 | pos.dec = dec;
|
---|
2562 |
|
---|
2563 | const Nova::ZdAzPosn hrz = Nova::GetHrzFromEqu(pos, now.JD());
|
---|
2564 |
|
---|
2565 | /*const*/ pair<vector<float>, pair<Time, float>> vis = GetVisibility(&pos);
|
---|
2566 | if (!vis.first.empty())
|
---|
2567 | {
|
---|
2568 | const Entry entry(name, vis.second.second, ccol);
|
---|
2569 | culmination.insert(make_pair(vis.second.first, entry));
|
---|
2570 | const vector<float> col = entry.GetColor(75, 15);
|
---|
2571 | vis.first.insert(vis.first.begin(), col.begin(), col.end());
|
---|
2572 | alt.push_back(vis.first);
|
---|
2573 |
|
---|
2574 | ccol++;
|
---|
2575 |
|
---|
2576 | /*const*/ pair<vector<float>, pair<Time, float>> lc = GetLightCondition(pos);
|
---|
2577 | if (!lc.first.empty())
|
---|
2578 | {
|
---|
2579 | const Entry entry2(name, lc.second.second, lcol);
|
---|
2580 | lightcond.insert(make_pair(lc.second.first, entry2));
|
---|
2581 | const vector<float> col2 = entry2.GetColor(100);
|
---|
2582 | lc.first.insert(lc.first.begin(), col2.begin(), col2.end());
|
---|
2583 | cur.push_back(lc.first);
|
---|
2584 |
|
---|
2585 | lcol++;
|
---|
2586 | }
|
---|
2587 | }
|
---|
2588 |
|
---|
2589 | string col = HTML::kWhite;
|
---|
2590 | if (hrz.zd<85)
|
---|
2591 | col = HTML::kRed;
|
---|
2592 | if (hrz.zd<65)
|
---|
2593 | col = HTML::kYellow;
|
---|
2594 | if (hrz.zd<30)
|
---|
2595 | col = HTML::kGreen;
|
---|
2596 |
|
---|
2597 | out2 << "<tr bgcolor='" << col << "'>";
|
---|
2598 | out2 << "<td>" << name << "</td>";
|
---|
2599 | if (hrz.zd<85)
|
---|
2600 | {
|
---|
2601 | out2 << "<td>" << hrz.zd << "°</td>";
|
---|
2602 | out2 << "<td>" << GetDir(hrz.az) << "</td>";
|
---|
2603 | }
|
---|
2604 | else
|
---|
2605 | out2 << "<td/><td/>";
|
---|
2606 | out2 << "</tr>";
|
---|
2607 | #endif
|
---|
2608 | const int32_t angle = fMoon.Angle(ra, dec);
|
---|
2609 |
|
---|
2610 | out << "<tr bgcolor='" << Moon::Color(angle) << "'>";
|
---|
2611 | out << "<td>" << name << "</td>";
|
---|
2612 | out << "<td>" << round(angle) << "°</td>";
|
---|
2613 | out << "</tr>";
|
---|
2614 | }
|
---|
2615 |
|
---|
2616 | for (auto it=culmination.begin(); it!=culmination.end(); it++)
|
---|
2617 | {
|
---|
2618 | const Entry &e = it->second;
|
---|
2619 | if (it!=culmination.begin())
|
---|
2620 | out3 << ", ";
|
---|
2621 | out3 << "<B#" << e.Col() << ">" << e.name << "</B>";
|
---|
2622 | if (e.value>0)
|
---|
2623 | out3 << " [" << nearbyint(90-e.value) << "°]";
|
---|
2624 | }
|
---|
2625 |
|
---|
2626 | out4 << setprecision(3);
|
---|
2627 |
|
---|
2628 | for (auto it=lightcond.begin(); it!=lightcond.end(); it++)
|
---|
2629 | {
|
---|
2630 | const Entry &e = it->second;
|
---|
2631 | if (it!=lightcond.begin())
|
---|
2632 | out4 << ", ";
|
---|
2633 | out4 << "<B#" << e.Col() << ">" << e.name << "</B>";
|
---|
2634 | if (e.value>0)
|
---|
2635 | out4 << " [" << nearbyint(e.value) << "]";
|
---|
2636 | }
|
---|
2637 |
|
---|
2638 | const Time st = fSun.fSunSet12;;
|
---|
2639 | const Time rs = fSun.fSunRise12;
|
---|
2640 |
|
---|
2641 | ostringstream title;
|
---|
2642 | title << st.GetAsStr("%H:%M");
|
---|
2643 | title << " / ";
|
---|
2644 | title << ((rs>st?rs-st:st-rs)/20).minutes();
|
---|
2645 | title << "' / ";
|
---|
2646 | title << rs.GetAsStr("%H:%M");
|
---|
2647 |
|
---|
2648 | out << '\n';
|
---|
2649 | out2 << '\n';
|
---|
2650 | out3 << '\n';
|
---|
2651 | out4 << '\n';
|
---|
2652 | out << HTML::kWhite << '\t' << Time()-now << '\n';
|
---|
2653 | out2 << HTML::kWhite << '\t' << Time()-now << '\n';
|
---|
2654 |
|
---|
2655 | WriteBinaryVec(now, "hist-visibility", alt, 75, 15, "Alt "+title.str());
|
---|
2656 | WriteBinaryVec(now, "hist-current-prediction", cur, 100, 0, "I " +title.str());
|
---|
2657 | }
|
---|
2658 | catch (const exception &e)
|
---|
2659 | {
|
---|
2660 | out << '\n';
|
---|
2661 | out2 << '\n';
|
---|
2662 | out << HTML::kWhite << '\t' << "ERROR - "+string(e.what()) << '\n';
|
---|
2663 | out2 << HTML::kWhite << '\t' << "ERROR - "+string(e.what()) << '\n';
|
---|
2664 | out3 << HTML::kWhite << '\t' << "ERROR - "+string(e.what()) << '\n';
|
---|
2665 | out4 << HTML::kWhite << '\t' << "ERROR - "+string(e.what()) << '\n';
|
---|
2666 | }
|
---|
2667 | #endif
|
---|
2668 |
|
---|
2669 | ofstream(fPath+"/moon.data") << out.str();
|
---|
2670 | ofstream(fPath+"/source-list.data") << out2.str();
|
---|
2671 | ofstream(fPath+"/visibility.data") << out3.str();
|
---|
2672 | ofstream(fPath+"/current-prediction.data") << out4.str();
|
---|
2673 | }
|
---|
2674 |
|
---|
2675 | int Execute()
|
---|
2676 | {
|
---|
2677 | Time now;
|
---|
2678 | if (now-fLastUpdate<boost::posix_time::seconds(1))
|
---|
2679 | return fDimDNS.online() ? kStateRunning : kStateDimNetworkNA;
|
---|
2680 | fLastUpdate=now;
|
---|
2681 |
|
---|
2682 | // ==============================================================
|
---|
2683 |
|
---|
2684 | bool reqscript = false;
|
---|
2685 |
|
---|
2686 | #ifdef HAVE_SQL
|
---|
2687 | try
|
---|
2688 | {
|
---|
2689 | const string query = Tools::Form("SELECT COUNT(*) FROM calendar.Data WHERE NOT u LIKE 'moon' AND y=%d AND m=%d AND d=%d",
|
---|
2690 | now.NightAsInt()/10000, (now.NightAsInt()/100)%100-1, now.NightAsInt()%100);
|
---|
2691 |
|
---|
2692 | const mysqlpp::StoreQueryResult res = Database(fDatabase).query(query).store();
|
---|
2693 |
|
---|
2694 | const uint32_t cnt = res[0][0];
|
---|
2695 |
|
---|
2696 | reqscript = cnt>0 && (fSun.state==3 || fSun.state==4);
|
---|
2697 | }
|
---|
2698 | catch (const exception &e)
|
---|
2699 | {
|
---|
2700 | Out() << e.what() << endl;
|
---|
2701 | }
|
---|
2702 | #endif
|
---|
2703 | // ==============================================================
|
---|
2704 |
|
---|
2705 | struct statvfs vfs;
|
---|
2706 | statvfs("/newdata", &vfs);
|
---|
2707 |
|
---|
2708 | const uint64_t free_newdata = vfs.f_bsize*vfs.f_bavail;
|
---|
2709 |
|
---|
2710 | // ==============================================================
|
---|
2711 |
|
---|
2712 | const bool data_taking =
|
---|
2713 | fDimMcp.state()==MCP::State::kTriggerOn ||
|
---|
2714 | fDimMcp.state()==MCP::State::kTakingData;
|
---|
2715 |
|
---|
2716 | const bool data_run =
|
---|
2717 | fMcpConfigurationName=="data" ||
|
---|
2718 | fMcpConfigurationName=="data-rt";
|
---|
2719 |
|
---|
2720 | const bool bias_on =
|
---|
2721 | fDimBiasControl.state()==BIAS::State::kRamping ||
|
---|
2722 | fDimBiasControl.state()==BIAS::State::kOverCurrent ||
|
---|
2723 | fDimBiasControl.state()==BIAS::State::kVoltageOn;
|
---|
2724 |
|
---|
2725 | const bool calibrated =
|
---|
2726 | fDimFeedback.state()>=Feedback::State::kCalibrated;
|
---|
2727 |
|
---|
2728 | const bool haderr = !fErrorList.empty();
|
---|
2729 |
|
---|
2730 | bool newerr = false;
|
---|
2731 |
|
---|
2732 | newerr |= SetError(!fDimDNS.online(),
|
---|
2733 | "<b><#darkred>DIM network not available</#></b>");
|
---|
2734 | newerr |= SetError(!fDimControl.online(),
|
---|
2735 | "<b>no dimctrl server available</b>");
|
---|
2736 | newerr |= SetError(fDimDataLogger.state()<20 || fDimDataLogger.state()>40,
|
---|
2737 | "<b>datalogger not ready</b>");
|
---|
2738 |
|
---|
2739 | newerr |= SetError(fDimControl.state()!=3 && reqscript,
|
---|
2740 | "<b>No script running during datataking time.</b>");
|
---|
2741 |
|
---|
2742 | //newerr |= SetError(fDimDriveControl.state()==Drive::State::kLocked,
|
---|
2743 | // "<b><#darkred>Drive in LOCKED state, drive was automatically parked</#></b>");
|
---|
2744 |
|
---|
2745 | newerr |= SetError(fDimDriveControl.state()>0xff && data_taking && data_run,
|
---|
2746 | "Drive in ERROR state during data-run");
|
---|
2747 | newerr |= SetError(fDriveControlMoonDist>155,
|
---|
2748 | "Moon within the field-of-view of the cones");
|
---|
2749 | newerr |= SetError(fDriveControlMoonDist>=0 && fDriveControlMoonDist<3,
|
---|
2750 | "Moon within the field-of-view of the camera");
|
---|
2751 |
|
---|
2752 | newerr |= SetError(fDimBiasControl.state()<BIAS::State::kRamping && data_taking && data_run,
|
---|
2753 | "BIAS not operating during data-run");
|
---|
2754 | newerr |= SetError(fDimBiasControl.state()==BIAS::State::kOverCurrent,
|
---|
2755 | "BIAS channels in OverCurrent");
|
---|
2756 | newerr |= SetError(fDimBiasControl.state()==BIAS::State::kNotReferenced,
|
---|
2757 | "BIAS voltage not at reference");
|
---|
2758 |
|
---|
2759 | newerr |= SetError(fDimFeedback.state()==Feedback::State::kOnStandby,
|
---|
2760 | "Feedback in standby due to high currents");
|
---|
2761 |
|
---|
2762 |
|
---|
2763 | newerr |= SetError(bias_on && calibrated && fBiasControlCurrentMed>115,
|
---|
2764 | "Median current (excl. crazy) exceeds 115µA/pix");
|
---|
2765 | newerr |= SetError(bias_on && calibrated && fBiasControlCurrentMax>160,
|
---|
2766 | "Maximum current (excl. crazy) exceeds 160µA/pix");
|
---|
2767 |
|
---|
2768 | newerr |= SetError(fFscControlHumidityAvg>60,
|
---|
2769 | "Average camera humidity exceed 60%");
|
---|
2770 |
|
---|
2771 | newerr |= SetError(!fPfMiniHumidityHist.empty() && fPfMiniHumidityHist.back()>50,
|
---|
2772 | "Camera humidity inside camera exceeds 50% (PFmini)");
|
---|
2773 | newerr |= SetError(!fTemperatureControlHist.empty() && (fTemperatureControlHist.back()<18.0 || fTemperatureControlHist.back()>22.0),
|
---|
2774 | "Container temperature outside [18.0;22.0]°C");
|
---|
2775 |
|
---|
2776 | newerr |= SetError(!fMagicWeatherHist[kHum].empty() && fMagicWeatherHist[kHum].back()>98 && fDimLidControl.state()==Lid::State::kOpen,
|
---|
2777 | "Outside humidity exceeds 98% while lid is open");
|
---|
2778 | newerr |= SetError(!fMagicWeatherHist[kGusts].empty() && fMagicWeatherHist[kGusts].back()>50 && (fDimDriveControl.state()==Drive::State::kTracking||fDimDriveControl.state()==Drive::State::kOnTrack),
|
---|
2779 | "Wind gusts exceed 50km/h during tracking");
|
---|
2780 |
|
---|
2781 | newerr |= SetError(fDimFscControl.state()>=FSC::State::kConnected && !fFscControlTemperatureHist.empty() && fFscControlTemperatureHist.back()>15,
|
---|
2782 | "Sensor temperature exceeds outside temperature by more than 15°C");
|
---|
2783 |
|
---|
2784 | newerr |= SetError(fFtmControlTriggerRateTooLow>0,
|
---|
2785 | "Trigger rate below 1Hz while trigger switched on");
|
---|
2786 |
|
---|
2787 | newerr |= SetError(fFtmControlState!=FTM::kFtmConfig && (fFtmControlState&FTM::kFtmLocked)==0,
|
---|
2788 | "FTM - clock conditioner not locked!");
|
---|
2789 |
|
---|
2790 | newerr |= SetError(fDimTimeCheck.state()==1,
|
---|
2791 | "Warning NTP time difference of drive PC exceeds 1s");
|
---|
2792 | newerr |= SetError(fDimTimeCheck.state()<1,
|
---|
2793 | "Warning timecheck not running");
|
---|
2794 |
|
---|
2795 | newerr |= SetError(fDimBiasControl.state()==BIAS::State::kVoltageOn &&
|
---|
2796 | fDimFeedback.state()<Feedback::State::kCalibrating &&
|
---|
2797 | fBiasControlVoltageMed>3,
|
---|
2798 | "Bias voltage switched on, but bias crate not calibrated");
|
---|
2799 |
|
---|
2800 | newerr |= SetError(fLastRunFinishedWithZeroEvents,
|
---|
2801 | "Last run finshed, but contained zero events.");
|
---|
2802 |
|
---|
2803 | newerr |= SetError(fFreeSpace<uint64_t(50000000000),
|
---|
2804 | "Less than 50GB disk space left on newdaq.");
|
---|
2805 |
|
---|
2806 | newerr |= SetError(free_newdata<uint64_t(800000000000),
|
---|
2807 | "Less than 800GB disk space left on /newdata.");
|
---|
2808 |
|
---|
2809 | newerr |= SetError(fDimPwrControl.state()==Power::State::kCoolingFailure,
|
---|
2810 | "Cooling unit reports failure!");
|
---|
2811 |
|
---|
2812 | for (auto it=fControlAlarmHist.begin(); it!=fControlAlarmHist.end(); it++)
|
---|
2813 | newerr |= SetError(it->time.IsValid(), it->msg);
|
---|
2814 | fControlAlarmHist.clean();;
|
---|
2815 |
|
---|
2816 | fLastRunFinishedWithZeroEvents = false;
|
---|
2817 |
|
---|
2818 | // FTM in Connected instead of Idle --> power cyclen
|
---|
2819 |
|
---|
2820 | /* // Check offline and disconnected status?
|
---|
2821 | Out() << fDimMcp << endl;
|
---|
2822 | Out() << fDimControl << endl;
|
---|
2823 | Out() << fDimDataLogger << endl;
|
---|
2824 | Out() << fDimDriveControl << endl;
|
---|
2825 | Out() << fDimFadControl << endl;
|
---|
2826 | Out() << fDimFtmControl << endl;
|
---|
2827 | Out() << fDimBiasControl << endl;
|
---|
2828 | Out() << fDimFeedback << endl;
|
---|
2829 | Out() << fDimRateControl << endl;
|
---|
2830 | Out() << fDimFscControl << endl;
|
---|
2831 | Out() << fDimMagicWeather << endl;
|
---|
2832 | Out() << fDimRateScan << endl;
|
---|
2833 | Out() << fDimChat << endl;
|
---|
2834 | */
|
---|
2835 |
|
---|
2836 | // FTU in error
|
---|
2837 | // FAD lost
|
---|
2838 |
|
---|
2839 | // --------------------------------------------------------------
|
---|
2840 | ostringstream out;
|
---|
2841 |
|
---|
2842 | if (newerr)
|
---|
2843 | {
|
---|
2844 | SetAudio("error");
|
---|
2845 |
|
---|
2846 | out << now.JavaDate() << '\n';
|
---|
2847 | out << HTML::kWhite << '\t';
|
---|
2848 | out << "<->" << fErrorHist.rget() << "<->";
|
---|
2849 | out << '\n';
|
---|
2850 |
|
---|
2851 | ofstream(fPath+"/errorhist.data") << out.str();
|
---|
2852 | }
|
---|
2853 |
|
---|
2854 | out.str("");
|
---|
2855 | out << Header(now) << '\t' << (!fErrorList.empty()) << '\t' << (fDimControl.state()>0) << '\n';
|
---|
2856 | out << setprecision(3);
|
---|
2857 | out << HTML::kWhite << '\t';
|
---|
2858 | for (auto it=fErrorList.begin(); it!=fErrorList.end(); it++)
|
---|
2859 | out << *it << "<br/>";
|
---|
2860 | out << '\n';
|
---|
2861 |
|
---|
2862 | if (haderr || !fErrorList.empty())
|
---|
2863 | ofstream(fPath+"/error.data") << out.str();
|
---|
2864 |
|
---|
2865 | // ==============================================================
|
---|
2866 |
|
---|
2867 | out.str("");
|
---|
2868 | out << Header(now) << '\t' << (!fErrorList.empty()) << '\t' << (fDimControl.state()>0) << '\n';
|
---|
2869 | out << setprecision(3);
|
---|
2870 |
|
---|
2871 | // -------------- System status --------------
|
---|
2872 | if (fDimDNS.online() && fDimMcp.state()>=MCP::State::kIdle) // Idle
|
---|
2873 | {
|
---|
2874 | string col = HTML::kBlue;
|
---|
2875 | switch (fMcpConfigurationState)
|
---|
2876 | {
|
---|
2877 | case MCP::State::kIdle:
|
---|
2878 | case DimState::kOffline:
|
---|
2879 | col = HTML::kWhite;
|
---|
2880 | break;
|
---|
2881 | case MCP::State::kConfiguring1:
|
---|
2882 | case MCP::State::kConfiguring2:
|
---|
2883 | case MCP::State::kConfiguring3:
|
---|
2884 | case MCP::State::kConfigured:
|
---|
2885 | case MCP::State::kTriggerOn:
|
---|
2886 | col = HTML::kBlue;
|
---|
2887 | break;
|
---|
2888 | case MCP::State::kTakingData:
|
---|
2889 | col = HTML::kBlue;
|
---|
2890 | if (fDimFadControl.state()==FAD::State::kRunInProgress)
|
---|
2891 | col = HTML::kGreen;
|
---|
2892 | break;
|
---|
2893 | }
|
---|
2894 |
|
---|
2895 | const bool other =
|
---|
2896 | fDimRateControl.state()==RateControl::State::kSettingGlobalThreshold ||
|
---|
2897 | fDimLidControl.state()==Lid::State::kMoving ||
|
---|
2898 | fDimRateScan.state()==RateScan::State::kInProgress;
|
---|
2899 |
|
---|
2900 | if (other)
|
---|
2901 | col = HTML::kBlue;
|
---|
2902 |
|
---|
2903 | out << col << '\t';
|
---|
2904 |
|
---|
2905 | if (!other)
|
---|
2906 | {
|
---|
2907 | const string conf = fMcpConfigurationName.length()>0?" ["+fMcpConfigurationName+"]":"";
|
---|
2908 | switch (fMcpConfigurationState)
|
---|
2909 | {
|
---|
2910 | case MCP::State::kIdle:
|
---|
2911 | out << "Idle" << conf;
|
---|
2912 | break;
|
---|
2913 | case MCP::State::kConfiguring1:
|
---|
2914 | case MCP::State::kConfiguring2:
|
---|
2915 | case MCP::State::kConfiguring3:
|
---|
2916 | out << "Configuring" << conf;
|
---|
2917 | break;
|
---|
2918 | case MCP::State::kConfigured:
|
---|
2919 | out << "Configured" << conf;
|
---|
2920 | break;
|
---|
2921 | case MCP::State::kTriggerOn:
|
---|
2922 | case MCP::State::kTakingData:
|
---|
2923 | out << fMcpConfigurationName;
|
---|
2924 | if (fFadControlDrsRuns[2]>0)
|
---|
2925 | out << "(" << fFadControlDrsRuns[2] << ")";
|
---|
2926 | break;
|
---|
2927 |
|
---|
2928 | case MCP::State::kCrateReset0:
|
---|
2929 | out << "Crate reset phase 0";
|
---|
2930 | break;
|
---|
2931 | case MCP::State::kCrateReset1:
|
---|
2932 | out << "Crate reset phase 1";
|
---|
2933 | break;
|
---|
2934 | case MCP::State::kCrateReset2:
|
---|
2935 | out << "Crate reset phase 2";
|
---|
2936 | break;
|
---|
2937 | case MCP::State::kCrateReset3:
|
---|
2938 | out << "Crate reset phase 3";
|
---|
2939 | break;
|
---|
2940 | }
|
---|
2941 | }
|
---|
2942 | else
|
---|
2943 | if (fDimRateControl.state()==RateControl::State::kSettingGlobalThreshold)
|
---|
2944 | out << "Calibrating threshold";
|
---|
2945 | else
|
---|
2946 | if (fDimRateScan.state()==RateScan::State::kInProgress)
|
---|
2947 | out << "Rate scan in progress";
|
---|
2948 | else
|
---|
2949 | if (fDimLidControl.state()==Lid::State::kMoving)
|
---|
2950 | out << "Lid moving";
|
---|
2951 |
|
---|
2952 |
|
---|
2953 | if (fMcpConfigurationState>MCP::State::kConfigured &&
|
---|
2954 | fDimRateControl.state()!=RateControl::State::kSettingGlobalThreshold)
|
---|
2955 | {
|
---|
2956 | ostringstream evt;
|
---|
2957 | if (fMcpConfigurationMaxEvents>0)
|
---|
2958 | {
|
---|
2959 | const int64_t de = int64_t(fMcpConfigurationMaxEvents) - int64_t(fFadControlNumEvents);
|
---|
2960 | if (de>=0 && fMcpConfigurationState==MCP::State::kTakingData)
|
---|
2961 | evt << de;
|
---|
2962 | else
|
---|
2963 | evt << fMcpConfigurationMaxEvents;
|
---|
2964 | }
|
---|
2965 | else
|
---|
2966 | {
|
---|
2967 | if (fMcpConfigurationState==MCP::State::kTakingData)
|
---|
2968 | {
|
---|
2969 | if (fFadControlNumEvents>2999)
|
---|
2970 | evt << floor(fFadControlNumEvents/1000) << 'k';
|
---|
2971 | else
|
---|
2972 | evt << fFadControlNumEvents;
|
---|
2973 | }
|
---|
2974 | }
|
---|
2975 |
|
---|
2976 | ostringstream tim;
|
---|
2977 | if (fMcpConfigurationMaxTime>0)
|
---|
2978 | {
|
---|
2979 | const uint32_t dt = (Time()-fMcpConfigurationRunStart).total_seconds();
|
---|
2980 | if (dt<=fMcpConfigurationMaxTime && fMcpConfigurationState==MCP::State::kTakingData)
|
---|
2981 | tim << fMcpConfigurationMaxTime-dt << 's';
|
---|
2982 | else
|
---|
2983 | tim << fMcpConfigurationMaxTime << 's';
|
---|
2984 | }
|
---|
2985 | else
|
---|
2986 | {
|
---|
2987 | if (fMcpConfigurationState==MCP::State::kTakingData)
|
---|
2988 | tim << fMcpConfigurationRunStart.SecondsTo();
|
---|
2989 | }
|
---|
2990 |
|
---|
2991 | const bool has_evt = !evt.str().empty();
|
---|
2992 | const bool has_tim = !tim.str().empty();
|
---|
2993 |
|
---|
2994 | if (has_evt || has_tim)
|
---|
2995 | out << " [";
|
---|
2996 | out << evt.str();
|
---|
2997 | if (has_evt && has_tim)
|
---|
2998 | out << '/';
|
---|
2999 | out << tim.str();
|
---|
3000 | if (has_evt || has_tim)
|
---|
3001 | out << ']';
|
---|
3002 | }
|
---|
3003 | }
|
---|
3004 | else
|
---|
3005 | out << HTML::kWhite;
|
---|
3006 | out << '\n';
|
---|
3007 |
|
---|
3008 | // ------------------ Drive -----------------
|
---|
3009 | if (fDimDNS.online() && fDimDriveControl.state()>=Drive::State::kInitialized) // Armed, Moving, Tracking, OnTrack, Error
|
---|
3010 | {
|
---|
3011 | const uint32_t dev = !fDriveControlTrackingDevHist.empty() ? round(fDriveControlTrackingDevHist.back()) : 0;
|
---|
3012 | const State rc = fDimDriveControl.description();
|
---|
3013 | string col = HTML::kGreen;
|
---|
3014 | if (fDimDriveControl.state()==Drive::State::kInitialized) // Armed
|
---|
3015 | col = HTML::kWhite;
|
---|
3016 | if (fDimDriveControl.state()>Drive::State::kInitialized && // Moving
|
---|
3017 | fDimDriveControl.state()<Drive::State::kTracking)
|
---|
3018 | col = HTML::kBlue;
|
---|
3019 | if (fDimDriveControl.state()==Drive::State::kTracking || // Tracking
|
---|
3020 | fDimDriveControl.state()==Drive::State::kOnTrack)
|
---|
3021 | {
|
---|
3022 | if (dev>60) // ~1.5mm
|
---|
3023 | col = HTML::kYellow;
|
---|
3024 | if (dev>120) // ~1/4 of a pixel ~ 2.5mm
|
---|
3025 | col = HTML::kRed;
|
---|
3026 | }
|
---|
3027 | if (fDimDriveControl.state()>0xff)
|
---|
3028 | col = HTML::kRed;
|
---|
3029 | out << col << '\t';
|
---|
3030 |
|
---|
3031 | //out << rc.name << '\t';
|
---|
3032 | out << fDriveControlPointingAz << ' ';
|
---|
3033 | out << fDriveControlPointingZd << "°";
|
---|
3034 | out << setprecision(2);
|
---|
3035 | if (fDimDriveControl.state()==Drive::State::kTracking ||
|
---|
3036 | fDimDriveControl.state()==Drive::State::kOnTrack) // Tracking
|
---|
3037 | {
|
---|
3038 | out << " ± " << dev << '"';
|
---|
3039 | if (!fDriveControlSourceName.empty())
|
---|
3040 | out << " [" << fDriveControlSourceName << ']';
|
---|
3041 | }
|
---|
3042 | if (fDimDriveControl.state()>Drive::State::kInitialized && // Moving
|
---|
3043 | fDimDriveControl.state()<Drive::State::kTracking)
|
---|
3044 | out << " ⟳";
|
---|
3045 | out << setprecision(3);
|
---|
3046 | }
|
---|
3047 | else
|
---|
3048 | out << HTML::kWhite << '\t';
|
---|
3049 |
|
---|
3050 | if (fSun.time.IsValid() && fMoon.time.IsValid())
|
---|
3051 | {
|
---|
3052 | if (fSun.visible)
|
---|
3053 | {
|
---|
3054 | out << " ☼";
|
---|
3055 | if (fDimDriveControl.state()<Drive::State::kInitialized)
|
---|
3056 | out << " [" << fSun.fSunSet12.MinutesTo() << "↓]";
|
---|
3057 | }
|
---|
3058 | else
|
---|
3059 | if (!fSun.visible && fMoon.visible)
|
---|
3060 | {
|
---|
3061 | out << " ☾";
|
---|
3062 | if (fDimDriveControl.state()<Drive::State::kInitialized)
|
---|
3063 | out << " [" << fMoon.disk << "%]";
|
---|
3064 | }
|
---|
3065 | }
|
---|
3066 | if (fDimDNS.online() && fDimDriveControl.state()>0xff)
|
---|
3067 | out << " <ERR>";
|
---|
3068 | if (fDimDNS.online() && fDimDriveControl.state()==Drive::State::kLocked)
|
---|
3069 | out << " ⊗";
|
---|
3070 | out << '\n';
|
---|
3071 |
|
---|
3072 | // ------------------- FSC ------------------
|
---|
3073 | if (fDimDNS.online() && fDimFscControl.state()>FSC::State::kDisconnected && !fFscControlTemperatureHist.empty())
|
---|
3074 | {
|
---|
3075 | string col = HTML::kGreen;
|
---|
3076 | if (fFscControlTemperatureHist.back()>9)
|
---|
3077 | col = HTML::kYellow;
|
---|
3078 | if (fFscControlTemperatureHist.back()>15)
|
---|
3079 | col = HTML::kRed;
|
---|
3080 |
|
---|
3081 | out << col << '\t' << fFscControlTemperatureHist.back() << '\n';
|
---|
3082 | }
|
---|
3083 | else
|
---|
3084 | out << HTML::kWhite << '\n';
|
---|
3085 |
|
---|
3086 | // --------------- MagicWeather -------------
|
---|
3087 | if (fDimDNS.online() && fDimMagicWeather.state()==MagicWeather::State::kReceiving && !fMagicWeatherHist[kWeatherBegin].empty())
|
---|
3088 | {
|
---|
3089 | /*
|
---|
3090 | const float diff = fMagicWeatherHist[kTemp].back()-fMagicWeatherHist[kDew].back();
|
---|
3091 | string col1 = HTML::kRed;
|
---|
3092 | if (diff>0.3)
|
---|
3093 | col1 = HTML::kYellow;
|
---|
3094 | if (diff>0.7)
|
---|
3095 | col1 = HTML::kGreen;
|
---|
3096 | */
|
---|
3097 |
|
---|
3098 | const float wind = fMagicWeatherHist[kGusts].back();
|
---|
3099 | const float hum = fMagicWeatherHist[kHum].back();
|
---|
3100 | string col = HTML::kGreen;
|
---|
3101 | if (wind>35 || hum>95)
|
---|
3102 | col = HTML::kYellow;
|
---|
3103 | if (wind>45 || hum>98)
|
---|
3104 | col = HTML::kRed;
|
---|
3105 |
|
---|
3106 | if (fDimRainSensor.state()==RainSensor::State::kValid && !fRainSensorDataHist.empty() && fRainSensorDataHist.back()>0)
|
---|
3107 | {
|
---|
3108 | out << HTML::kRed << '\t';
|
---|
3109 | out << "RAIN" << '\t';
|
---|
3110 | }
|
---|
3111 | else
|
---|
3112 | {
|
---|
3113 | out << col << '\t';
|
---|
3114 | out << fMagicWeatherHist[kHum].back() << '\t';
|
---|
3115 | }
|
---|
3116 | out << setprecision(2);
|
---|
3117 | out << fMagicWeatherHist[kGusts].back() << '\n';
|
---|
3118 | out << setprecision(3);
|
---|
3119 | }
|
---|
3120 | else
|
---|
3121 | out << HTML::kWhite << "\n";
|
---|
3122 |
|
---|
3123 | // --------------- FtmControl -------------
|
---|
3124 | if (fDimDNS.online() && fDimFtmControl.state()==FTM::State::kTriggerOn)
|
---|
3125 | {
|
---|
3126 | string col = HTML::kGreen;
|
---|
3127 | if (!fFtmControlTriggerRateHist.empty())
|
---|
3128 | {
|
---|
3129 | if (fFtmControlTriggerRateHist.back()<15)
|
---|
3130 | col = HTML::kYellow;
|
---|
3131 | if (fFtmControlTriggerRateHist.back()>100)
|
---|
3132 | col = HTML::kRed;
|
---|
3133 |
|
---|
3134 | out << col << '\t' << fFtmControlTriggerRateHist.back() << " Hz";
|
---|
3135 | }
|
---|
3136 |
|
---|
3137 | if (bias_on)
|
---|
3138 | out << " (" << setprecision(4) << fFtmPatchThresholdMed << ')';
|
---|
3139 | out << '\n';
|
---|
3140 | }
|
---|
3141 | else
|
---|
3142 | out << HTML::kWhite << '\n';
|
---|
3143 |
|
---|
3144 | // --------------- BiasControl -------------
|
---|
3145 | const bool bias_off = fDimBiasControl.state()==BIAS::State::kVoltageOff;
|
---|
3146 | const bool bias_oc = fDimBiasControl.state()==BIAS::State::kOverCurrent;
|
---|
3147 |
|
---|
3148 | if (fDimDNS.online() && (bias_on || bias_off))
|
---|
3149 | {
|
---|
3150 |
|
---|
3151 | string col = fBiasControlVoltageMed>3?HTML::kGreen:HTML::kWhite;
|
---|
3152 | if (bias_on)
|
---|
3153 | {
|
---|
3154 | if (fBiasControlCurrentMed>95 || fBiasControlCurrentMax>135)
|
---|
3155 | col = HTML::kYellow;
|
---|
3156 | if (fBiasControlCurrentMed>100 || fBiasControlCurrentMax>140)
|
---|
3157 | col = HTML::kRed;
|
---|
3158 | }
|
---|
3159 |
|
---|
3160 | // Bias in overcurrent => Red
|
---|
3161 | if (bias_oc)
|
---|
3162 | col = HTML::kRed;
|
---|
3163 |
|
---|
3164 | // MCP in ReadyForDatataking/Configuring/Configured/TriggerOn/TakingData
|
---|
3165 | // and Bias not in "data-taking state' => Red
|
---|
3166 | if (fMcpConfigurationState>MCP::State::kIdle && !bias_on)
|
---|
3167 | col = HTML::kWhite;
|
---|
3168 |
|
---|
3169 | const bool cal = fDimFeedback.state()>=Feedback::State::kCalibrated;
|
---|
3170 |
|
---|
3171 | // Feedback is currently calibrating => Blue
|
---|
3172 | if (fDimFeedback.state()==Feedback::State::kCalibrating)
|
---|
3173 | {
|
---|
3174 | out << HTML::kBlue << '\t';
|
---|
3175 | out << "***\t";
|
---|
3176 | out << "***\t";
|
---|
3177 | }
|
---|
3178 | else
|
---|
3179 | {
|
---|
3180 | out << col << '\t';
|
---|
3181 | out << setprecision(fBiasControlCurrentMed<100?2:3);
|
---|
3182 | out << (bias_off ? 0 : (fBiasControlCurrentMed<10?fBiasControlCurrentMed:floor(fBiasControlCurrentMed))) << '\t';
|
---|
3183 | if (bias_oc)
|
---|
3184 | out << "(OC) ";
|
---|
3185 | else
|
---|
3186 | {
|
---|
3187 | if (cal)
|
---|
3188 | {
|
---|
3189 | out << setprecision(fBiasControlCurrentMax<100?2:3);
|
---|
3190 | out << (bias_off ? 0 : (fBiasControlCurrentMax<10?fBiasControlCurrentMax:floor(fBiasControlCurrentMax)));
|
---|
3191 | }
|
---|
3192 | else
|
---|
3193 | out << "— ";
|
---|
3194 | }
|
---|
3195 | out << '\t';
|
---|
3196 | }
|
---|
3197 | if (cal && fDimFeedback.state()!=Feedback::State::kCalibrating)
|
---|
3198 | out << setprecision(2) << fBiasControlPowerTot << " W";
|
---|
3199 | else
|
---|
3200 | out << setprecision(3) << (bias_off ? 0 : fBiasControlVoltageMed) << " V";
|
---|
3201 | out << '\n';
|
---|
3202 | }
|
---|
3203 | else
|
---|
3204 | out << HTML::kWhite << '\n';
|
---|
3205 |
|
---|
3206 | ofstream(fPath+"/fact.data") << out.str();
|
---|
3207 |
|
---|
3208 | // ==============================================================
|
---|
3209 |
|
---|
3210 | out.str("");
|
---|
3211 | out << Header(now) << '\t' << (!fErrorList.empty()) << '\t' << (fDimControl.state()>0) << '\n';
|
---|
3212 |
|
---|
3213 | if (!fDimDNS.online())
|
---|
3214 | out << HTML::kWhite << "\tOffline\n\n\n\n\n\n\n\n\n\n\n\n\n";
|
---|
3215 | else
|
---|
3216 | {
|
---|
3217 | ostringstream dt;
|
---|
3218 | dt << (Time()-fRunTime);
|
---|
3219 |
|
---|
3220 | out << HTML::kGreen << '\t' << fDimDNS.version() << '\n';
|
---|
3221 |
|
---|
3222 | out << GetStateHtml(fDimControl, 0);
|
---|
3223 | out << GetStateHtml(fDimMcp, MCP::State::kConnected);
|
---|
3224 | out << GetStateHtml(fDimDataLogger, 1);
|
---|
3225 | out << GetStateHtml(fDimDriveControl, Drive::State::kConnected);
|
---|
3226 | out << GetStateHtml(fDimTimeCheck, 1);
|
---|
3227 | out << GetStateHtml(fDimFadControl, FAD::State::kConnected);
|
---|
3228 | out << GetStateHtml(fDimFtmControl, FTM::State::kConnected);
|
---|
3229 | out << GetStateHtml(fDimBiasControl, BIAS::State::kConnected);
|
---|
3230 | out << GetStateHtml(fDimFeedback, Feedback::State::kConnected);
|
---|
3231 | out << GetStateHtml(fDimRateControl, RateControl::State::kConnected);
|
---|
3232 | out << GetStateHtml(fDimFscControl, FSC::State::kConnected);
|
---|
3233 | out << GetStateHtml(fDimPfMiniControl, PFmini::State::kConnected);
|
---|
3234 | out << GetStateHtml(fDimBiasTemp, BiasTemp::State::kConnected);
|
---|
3235 | out << GetStateHtml(fDimGpsControl, GPS::State::kConnected);
|
---|
3236 | out << GetStateHtml(fDimSqmControl, SQM::State::kConnected);
|
---|
3237 | out << GetStateHtml(fDimAgilentControl24, Agilent::State::kVoltageOff);
|
---|
3238 | out << GetStateHtml(fDimAgilentControl50, Agilent::State::kVoltageOff);
|
---|
3239 | out << GetStateHtml(fDimAgilentControl80, Agilent::State::kVoltageOff);
|
---|
3240 | out << GetStateHtml(fDimPwrControl, Power::State::kSystemOff);
|
---|
3241 | out << GetStateHtml(fDimLidControl, Lid::State::kConnected);
|
---|
3242 | out << GetStateHtml(fDimRateScan, RateScan::State::kConnected);
|
---|
3243 | out << GetStateHtml(fDimMagicWeather, MagicWeather::State::kConnected);
|
---|
3244 | out << GetStateHtml(fDimTngWeather, TNGWeather::State::kConnected);
|
---|
3245 | out << GetStateHtml(fDimMagicLidar, MagicLidar::State::kConnected);
|
---|
3246 | out << GetStateHtml(fDimTemperature, Temperature::State::kValid);
|
---|
3247 | out << GetStateHtml(fDimRainSensor, RainSensor::State::kValid);
|
---|
3248 | out << GetStateHtml(fDimChat, 0);
|
---|
3249 | out << GetStateHtml(fDimSkypeClient, 1);
|
---|
3250 |
|
---|
3251 | string col = HTML::kRed;
|
---|
3252 | if (fFreeSpace>uint64_t(199999999999))
|
---|
3253 | col = HTML::kYellow;
|
---|
3254 | if (fFreeSpace>uint64_t(999999999999))
|
---|
3255 | col = HTML::kGreen;
|
---|
3256 | if (fFreeSpace==UINT64_MAX)
|
---|
3257 | col = HTML::kWhite;
|
---|
3258 |
|
---|
3259 | out << col << '\t' << Tools::Scientific(fFreeSpace) << "B\n";
|
---|
3260 |
|
---|
3261 | col = HTML::kRed;
|
---|
3262 | if (free_newdata > uint64_t(999999999999))
|
---|
3263 | col = HTML::kYellow;
|
---|
3264 | if (free_newdata > uint64_t(149999999999))
|
---|
3265 | col = HTML::kGreen;
|
---|
3266 | if (free_newdata == UINT64_MAX)
|
---|
3267 | col = HTML::kWhite;
|
---|
3268 |
|
---|
3269 | out << col << '\t' << Tools::Scientific(free_newdata) << "B\n";
|
---|
3270 |
|
---|
3271 | out << HTML::kGreen << '\t' << dt.str().substr(0, dt.str().length()-7) << '\n';
|
---|
3272 | }
|
---|
3273 |
|
---|
3274 | ofstream(fPath+"/status.data") << out.str();
|
---|
3275 |
|
---|
3276 | if (now-fLastAstroCalc>boost::posix_time::seconds(15))
|
---|
3277 | {
|
---|
3278 | UpdateAstronomy();
|
---|
3279 | fLastAstroCalc = now;
|
---|
3280 | }
|
---|
3281 |
|
---|
3282 | return fDimDNS.online() ? kStateRunning : kStateDimNetworkNA;
|
---|
3283 | }
|
---|
3284 |
|
---|
3285 |
|
---|
3286 | public:
|
---|
3287 | StateMachineSmartFACT(ostream &out=cout) : StateMachineDim(out, fIsServer?"SMART_FACT":""),
|
---|
3288 | fLastAstroCalc(boost::date_time::neg_infin),
|
---|
3289 | fPath("www/smartfact/data"),
|
---|
3290 | fMcpConfigurationState(DimState::kOffline),
|
---|
3291 | fMcpConfigurationMaxTime(0),
|
---|
3292 | fMcpConfigurationMaxEvents(0),
|
---|
3293 | fLastRunFinishedWithZeroEvents(false),
|
---|
3294 | fTngWeatherDustTime(Time::none),
|
---|
3295 | fBiasControlVoltageMed(0),
|
---|
3296 | fBiasControlCurrentMed(0),
|
---|
3297 | fBiasControlCurrentMax(0),
|
---|
3298 | fFscControlHumidityAvg(0),
|
---|
3299 | fDriveControlMoonDist(-1),
|
---|
3300 | fFadControlNumEvents(0),
|
---|
3301 | fFadControlDrsRuns(3),
|
---|
3302 | fFtmControlState(FTM::kFtmLocked),
|
---|
3303 | fRateScanDataId(0),
|
---|
3304 | fRateScanBoard(0),
|
---|
3305 | fFreeSpace(UINT64_MAX),
|
---|
3306 | // ---
|
---|
3307 | fDimMcp ("MCP"),
|
---|
3308 | fDimDataLogger ("DATA_LOGGER"),
|
---|
3309 | fDimDriveControl ("DRIVE_CONTROL"),
|
---|
3310 | fDimTimeCheck ("TIME_CHECK"),
|
---|
3311 | fDimMagicWeather ("MAGIC_WEATHER"),
|
---|
3312 | fDimMagicLidar ("MAGIC_LIDAR"),
|
---|
3313 | fDimTngWeather ("TNG_WEATHER"),
|
---|
3314 | fDimTemperature ("TEMPERATURE"),
|
---|
3315 | fDimRainSensor ("RAIN_SENSOR"),
|
---|
3316 | fDimFeedback ("FEEDBACK"),
|
---|
3317 | fDimBiasControl ("BIAS_CONTROL"),
|
---|
3318 | fDimFtmControl ("FTM_CONTROL"),
|
---|
3319 | fDimFadControl ("FAD_CONTROL"),
|
---|
3320 | fDimFscControl ("FSC_CONTROL"),
|
---|
3321 | fDimPfMiniControl ("PFMINI_CONTROL"),
|
---|
3322 | fDimBiasTemp ("BIAS_TEMP"),
|
---|
3323 | fDimGpsControl ("GPS_CONTROL"),
|
---|
3324 | fDimSqmControl ("SQM_CONTROL"),
|
---|
3325 | fDimAgilentControl24("AGILENT_CONTROL_24V"),
|
---|
3326 | fDimAgilentControl50("AGILENT_CONTROL_50V"),
|
---|
3327 | fDimAgilentControl80("AGILENT_CONTROL_80V"),
|
---|
3328 | fDimPwrControl ("PWR_CONTROL"),
|
---|
3329 | fDimLidControl ("LID_CONTROL"),
|
---|
3330 | fDimRateControl ("RATE_CONTROL"),
|
---|
3331 | fDimRateScan ("RATE_SCAN"),
|
---|
3332 | fDimChat ("CHAT"),
|
---|
3333 | fDimSkypeClient ("SKYPE_CLIENT")
|
---|
3334 | {
|
---|
3335 | fDimDNS.Subscribe(*this);
|
---|
3336 | fDimControl.Subscribe(*this);
|
---|
3337 | fDimMcp.Subscribe(*this);
|
---|
3338 | fDimDataLogger.Subscribe(*this);
|
---|
3339 | fDimDriveControl.Subscribe(*this);
|
---|
3340 | fDimTimeCheck.Subscribe(*this);
|
---|
3341 | fDimMagicWeather.Subscribe(*this);
|
---|
3342 | fDimMagicLidar.Subscribe(*this);
|
---|
3343 | fDimTngWeather.Subscribe(*this);
|
---|
3344 | fDimTemperature.Subscribe(*this);
|
---|
3345 | fDimRainSensor.Subscribe(*this);
|
---|
3346 | fDimFeedback.Subscribe(*this);
|
---|
3347 | fDimBiasControl.Subscribe(*this);
|
---|
3348 | fDimFtmControl.Subscribe(*this);
|
---|
3349 | fDimFadControl.Subscribe(*this);
|
---|
3350 | fDimFscControl.Subscribe(*this);
|
---|
3351 | fDimPfMiniControl.Subscribe(*this);
|
---|
3352 | fDimBiasTemp.Subscribe(*this);
|
---|
3353 | fDimGpsControl.Subscribe(*this);
|
---|
3354 | fDimSqmControl.Subscribe(*this);
|
---|
3355 | fDimAgilentControl24.Subscribe(*this);
|
---|
3356 | fDimAgilentControl50.Subscribe(*this);
|
---|
3357 | fDimAgilentControl80.Subscribe(*this);
|
---|
3358 | fDimPwrControl.Subscribe(*this);
|
---|
3359 | fDimLidControl.Subscribe(*this);
|
---|
3360 | fDimRateControl.Subscribe(*this);
|
---|
3361 | fDimRateScan.Subscribe(*this);
|
---|
3362 | fDimChat.Subscribe(*this);
|
---|
3363 | fDimSkypeClient.Subscribe(*this);
|
---|
3364 |
|
---|
3365 | fDimFscControl.SetCallback(bind(&StateMachineSmartFACT::HandleFscControlStateChange, this, placeholders::_1));
|
---|
3366 | //fDimFtmControl.SetCallback(bind(&StateMachineSmartFACT::HandleFtmControlStateChange, this));
|
---|
3367 | fDimDriveControl.SetCallback(bind(&StateMachineSmartFACT::HandleDriveControlStateChange, this, placeholders::_1));
|
---|
3368 | fDimControl.SetCallback(bind(&StateMachineSmartFACT::HandleControlStateChange, this, placeholders::_1));
|
---|
3369 | fDimControl.AddCallback("dotest.dim", bind(&StateMachineSmartFACT::HandleDoTest, this, placeholders::_1));
|
---|
3370 |
|
---|
3371 | Subscribe("DIM_CONTROL/MESSAGE")
|
---|
3372 | (bind(&StateMachineSmartFACT::HandleDimControlMessage, this, placeholders::_1));
|
---|
3373 |
|
---|
3374 | Subscribe("MCP/CONFIGURATION")
|
---|
3375 | (bind(&StateMachineSmartFACT::HandleMcpConfiguration, this, placeholders::_1));
|
---|
3376 |
|
---|
3377 | Subscribe("DRIVE_CONTROL/POINTING_POSITION")
|
---|
3378 | (bind(&StateMachineSmartFACT::HandleDrivePointing, this, placeholders::_1));
|
---|
3379 | Subscribe("DRIVE_CONTROL/TRACKING_POSITION")
|
---|
3380 | (bind(&StateMachineSmartFACT::HandleDriveTracking, this, placeholders::_1));
|
---|
3381 | Subscribe("DRIVE_CONTROL/SOURCE_POSITION")
|
---|
3382 | (bind(&StateMachineSmartFACT::HandleDriveSource, this, placeholders::_1));
|
---|
3383 |
|
---|
3384 | Subscribe("FSC_CONTROL/TEMPERATURE")
|
---|
3385 | (bind(&StateMachineSmartFACT::HandleFscTemperature, this, placeholders::_1));
|
---|
3386 | Subscribe("FSC_CONTROL/HUMIDITY")
|
---|
3387 | (bind(&StateMachineSmartFACT::HandleFscHumidity, this, placeholders::_1));
|
---|
3388 | Subscribe("FSC_CONTROL/BIAS_TEMP")
|
---|
3389 | (bind(&StateMachineSmartFACT::HandleFscBiasTemp, this, placeholders::_1));
|
---|
3390 |
|
---|
3391 | Subscribe("PFMINI_CONTROL/DATA")
|
---|
3392 | (bind(&StateMachineSmartFACT::HandlePfMiniData, this, placeholders::_1));
|
---|
3393 |
|
---|
3394 | Subscribe("BIAS_TEMP/DATA")
|
---|
3395 | (bind(&StateMachineSmartFACT::HandleBiasTemp, this, placeholders::_1));
|
---|
3396 |
|
---|
3397 | Subscribe("GPS_CONTROL/NEMA")
|
---|
3398 | (bind(&StateMachineSmartFACT::HandleGpsNema, this, placeholders::_1));
|
---|
3399 |
|
---|
3400 | Subscribe("SQM_CONTROL/DATA")
|
---|
3401 | (bind(&StateMachineSmartFACT::HandleSqmData, this, placeholders::_1));
|
---|
3402 |
|
---|
3403 | Subscribe("TEMPERATURE/DATA")
|
---|
3404 | (bind(&StateMachineSmartFACT::HandleTemperatureData, this, placeholders::_1));
|
---|
3405 |
|
---|
3406 | Subscribe("AGILENT_CONTROL_24V/DATA")
|
---|
3407 | (bind(&StateMachineSmartFACT::HandleAgilentData, this, placeholders::_1, "24"));
|
---|
3408 | Subscribe("AGILENT_CONTROL_50V/DATA")
|
---|
3409 | (bind(&StateMachineSmartFACT::HandleAgilentData, this, placeholders::_1, "50"));
|
---|
3410 | Subscribe("AGILENT_CONTROL_80V/DATA")
|
---|
3411 | (bind(&StateMachineSmartFACT::HandleAgilentData, this, placeholders::_1, "80"));
|
---|
3412 |
|
---|
3413 | Subscribe("MAGIC_WEATHER/DATA")
|
---|
3414 | (bind(&StateMachineSmartFACT::HandleMagicWeatherData, this, placeholders::_1));
|
---|
3415 | Subscribe("TNG_WEATHER/DUST")
|
---|
3416 | (bind(&StateMachineSmartFACT::HandleTngWeatherDust, this, placeholders::_1));
|
---|
3417 | Subscribe("RAIN_SENSOR/DATA")
|
---|
3418 | (bind(&StateMachineSmartFACT::HandleRainSensorData, this, placeholders::_1));
|
---|
3419 |
|
---|
3420 | Subscribe("FEEDBACK/CALIBRATED_CURRENTS")
|
---|
3421 | (bind(&StateMachineSmartFACT::HandleFeedbackCalibratedCurrents, this, placeholders::_1));
|
---|
3422 |
|
---|
3423 | Subscribe("BIAS_CONTROL/VOLTAGE")
|
---|
3424 | (bind(&StateMachineSmartFACT::HandleBiasVoltage, this, placeholders::_1));
|
---|
3425 | Subscribe("BIAS_CONTROL/CURRENT")
|
---|
3426 | (bind(&StateMachineSmartFACT::HandleBiasCurrent, this, placeholders::_1));
|
---|
3427 |
|
---|
3428 | Subscribe("FAD_CONTROL/CONNECTIONS")
|
---|
3429 | (bind(&StateMachineSmartFACT::HandleFadConnections, this, placeholders::_1));
|
---|
3430 | Subscribe("FAD_CONTROL/EVENTS")
|
---|
3431 | (bind(&StateMachineSmartFACT::HandleFadEvents, this, placeholders::_1));
|
---|
3432 | Subscribe("FAD_CONTROL/START_RUN")
|
---|
3433 | (bind(&StateMachineSmartFACT::HandleFadStartRun, this, placeholders::_1));
|
---|
3434 | Subscribe("FAD_CONTROL/DRS_RUNS")
|
---|
3435 | (bind(&StateMachineSmartFACT::HandleFadDrsRuns, this, placeholders::_1));
|
---|
3436 | Subscribe("FAD_CONTROL/EVENT_DATA")
|
---|
3437 | (bind(&StateMachineSmartFACT::HandleFadEventData, this, placeholders::_1));
|
---|
3438 | Subscribe("FAD_CONTROL/STATS")
|
---|
3439 | (bind(&StateMachineSmartFACT::HandleStats, this, placeholders::_1));
|
---|
3440 |
|
---|
3441 | Subscribe("DATA_LOGGER/STATS")
|
---|
3442 | (bind(&StateMachineSmartFACT::HandleStats, this, placeholders::_1));
|
---|
3443 |
|
---|
3444 | Subscribe("FTM_CONTROL/TRIGGER_RATES")
|
---|
3445 | (bind(&StateMachineSmartFACT::HandleFtmTriggerRates, this, placeholders::_1));
|
---|
3446 | Subscribe("FTM_CONTROL/STATIC_DATA")
|
---|
3447 | (bind(&StateMachineSmartFACT::HandleFtmStaticData, this, placeholders::_1));
|
---|
3448 | Subscribe("FTM_CONTROL/FTU_LIST")
|
---|
3449 | (bind(&StateMachineSmartFACT::HandleFtmFtuList, this, placeholders::_1));
|
---|
3450 |
|
---|
3451 | Subscribe("RATE_CONTROL/THRESHOLD")
|
---|
3452 | (bind(&StateMachineSmartFACT::HandleRateControlThreshold,this, placeholders::_1));
|
---|
3453 |
|
---|
3454 | Subscribe("RATE_SCAN/DATA")
|
---|
3455 | (bind(&StateMachineSmartFACT::HandleRateScanData, this, placeholders::_1));
|
---|
3456 |
|
---|
3457 | Subscribe("CHAT/MESSAGE")
|
---|
3458 | (bind(&StateMachineSmartFACT::HandleChatMsg, this, placeholders::_1));
|
---|
3459 |
|
---|
3460 |
|
---|
3461 | // =================================================================
|
---|
3462 |
|
---|
3463 | // State names
|
---|
3464 | AddStateName(kStateDimNetworkNA, "DimNetworkNotAvailable",
|
---|
3465 | "The Dim DNS is not reachable.");
|
---|
3466 |
|
---|
3467 | AddStateName(kStateRunning, "Running", "");
|
---|
3468 |
|
---|
3469 | // =================================================================
|
---|
3470 |
|
---|
3471 | AddEvent("PRINT")
|
---|
3472 | (bind(&StateMachineSmartFACT::Print, this))
|
---|
3473 | ("Print a list of the states of all connected servers.");
|
---|
3474 |
|
---|
3475 | }
|
---|
3476 | int EvalOptions(Configuration &conf)
|
---|
3477 | {
|
---|
3478 | if (!fPixelMap.Read(conf.Get<string>("pixel-map-file")))
|
---|
3479 | {
|
---|
3480 | Error("Reading mapping table from "+conf.Get<string>("pixel-map-file")+" failed.");
|
---|
3481 | return 1;
|
---|
3482 | }
|
---|
3483 |
|
---|
3484 | fPath = conf.Get<string>("path");
|
---|
3485 | fDatabase = conf.Get<string>("source-database");
|
---|
3486 |
|
---|
3487 | struct stat st;
|
---|
3488 | if (stat(fPath.c_str(), &st))
|
---|
3489 | {
|
---|
3490 | Error(fPath+" does not exist!");
|
---|
3491 | return 2;
|
---|
3492 | }
|
---|
3493 |
|
---|
3494 | if ((st.st_mode&S_IFDIR)==0)
|
---|
3495 | {
|
---|
3496 | Error(fPath+" not a directory!");
|
---|
3497 | return 3;
|
---|
3498 | }
|
---|
3499 |
|
---|
3500 | if ((st.st_mode&S_IWUSR)==0)
|
---|
3501 | {
|
---|
3502 | Error(fPath+" has no write permission!");
|
---|
3503 | return 4;
|
---|
3504 | }
|
---|
3505 |
|
---|
3506 | if ((st.st_mode&S_IXUSR)==0)
|
---|
3507 | {
|
---|
3508 | Error(fPath+" has no execute permission!");
|
---|
3509 | return 5;
|
---|
3510 | }
|
---|
3511 |
|
---|
3512 | ostringstream out;
|
---|
3513 | out << Time().JavaDate() << '\n';
|
---|
3514 |
|
---|
3515 | ofstream(fPath+"/error.data") << out.str();
|
---|
3516 |
|
---|
3517 | return -1;
|
---|
3518 | }
|
---|
3519 | };
|
---|
3520 |
|
---|
3521 | bool StateMachineSmartFACT::fIsServer = false;
|
---|
3522 |
|
---|
3523 | // ------------------------------------------------------------------------
|
---|
3524 |
|
---|
3525 | #include "Main.h"
|
---|
3526 |
|
---|
3527 | template<class T>
|
---|
3528 | int RunShell(Configuration &conf)
|
---|
3529 | {
|
---|
3530 | StateMachineSmartFACT::fIsServer = !conf.Get<bool>("client");
|
---|
3531 | return Main::execute<T, StateMachineSmartFACT>(conf);
|
---|
3532 | }
|
---|
3533 |
|
---|
3534 | void SetupConfiguration(Configuration &conf)
|
---|
3535 | {
|
---|
3536 | po::options_description control("Smart FACT");
|
---|
3537 | control.add_options()
|
---|
3538 | ("pixel-map-file", var<string>()->required(), "Pixel mapping file. Used here to get the default reference voltage")
|
---|
3539 | ("path", var<string>("www/smartfact/data"), "Output path for the data-files")
|
---|
3540 | ("source-database", var<string>(""), "Database link as in\n\tuser:password@server[:port]/database[/comp].")
|
---|
3541 | ("client", po_bool(false), "For a standalone client choose this option.")
|
---|
3542 | ;
|
---|
3543 |
|
---|
3544 | conf.AddOptions(control);
|
---|
3545 | }
|
---|
3546 |
|
---|
3547 | /*
|
---|
3548 | Extract usage clause(s) [if any] for SYNOPSIS.
|
---|
3549 | Translators: "Usage" and "or" here are patterns (regular expressions) which
|
---|
3550 | are used to match the usage synopsis in program output. An example from cp
|
---|
3551 | (GNU coreutils) which contains both strings:
|
---|
3552 | Usage: cp [OPTION]... [-T] SOURCE DEST
|
---|
3553 | or: cp [OPTION]... SOURCE... DIRECTORY
|
---|
3554 | or: cp [OPTION]... -t DIRECTORY SOURCE...
|
---|
3555 | */
|
---|
3556 | void PrintUsage()
|
---|
3557 | {
|
---|
3558 | cout <<
|
---|
3559 | "SmartFACT is a tool writing the files needed for the SmartFACT web interface.\n"
|
---|
3560 | "\n"
|
---|
3561 | "The default is that the program is started without user intercation. "
|
---|
3562 | "All actions are supposed to arrive as DimCommands. Using the -c "
|
---|
3563 | "option, a local shell can be initialized. With h or help a short "
|
---|
3564 | "help message about the usuage can be brought to the screen.\n"
|
---|
3565 | "\n"
|
---|
3566 | "Usage: smartfact [-c type] [OPTIONS]\n"
|
---|
3567 | " or: smartfact [OPTIONS]\n";
|
---|
3568 | cout << endl;
|
---|
3569 | }
|
---|
3570 |
|
---|
3571 | void PrintHelp()
|
---|
3572 | {
|
---|
3573 | Main::PrintHelp<StateMachineSmartFACT>();
|
---|
3574 |
|
---|
3575 | /* Additional help text which is printed after the configuration
|
---|
3576 | options goes here */
|
---|
3577 |
|
---|
3578 | /*
|
---|
3579 | cout << "bla bla bla" << endl << endl;
|
---|
3580 | cout << endl;
|
---|
3581 | cout << "Environment:" << endl;
|
---|
3582 | cout << "environment" << endl;
|
---|
3583 | cout << endl;
|
---|
3584 | cout << "Examples:" << endl;
|
---|
3585 | cout << "test exam" << endl;
|
---|
3586 | cout << endl;
|
---|
3587 | cout << "Files:" << endl;
|
---|
3588 | cout << "files" << endl;
|
---|
3589 | cout << endl;
|
---|
3590 | */
|
---|
3591 | }
|
---|
3592 |
|
---|
3593 | int main(int argc, const char* argv[])
|
---|
3594 | {
|
---|
3595 | Configuration conf(argv[0]);
|
---|
3596 | conf.SetPrintUsage(PrintUsage);
|
---|
3597 | Main::SetupConfiguration(conf);
|
---|
3598 | SetupConfiguration(conf);
|
---|
3599 |
|
---|
3600 | if (!conf.DoParse(argc, argv, PrintHelp))
|
---|
3601 | return 127;
|
---|
3602 |
|
---|
3603 | if (!conf.Has("console"))
|
---|
3604 | return RunShell<LocalStream>(conf);
|
---|
3605 |
|
---|
3606 | if (conf.Get<int>("console")==0)
|
---|
3607 | return RunShell<LocalShell>(conf);
|
---|
3608 | else
|
---|
3609 | return RunShell<LocalConsole>(conf);
|
---|
3610 |
|
---|
3611 | return 0;
|
---|
3612 | }
|
---|