Changeset 14722

Timestamp:

12/04/12 10:10:03 (12 years ago)

Author:

tbretz

Message:

Added date on x-axis, fixed a problem with summer/winter time (hopefully)

File:

• trunk/FACT++/src/makeplots.cc (modified) (5 diffs)

trunk/FACT++/src/makeplots.cc

@@ -89 +89 @@
         "the plots will refer to the night 1974-09-09/1974-09-10.\n"
         "The advantage is that specification of the date as in\n"
-        "1974-09-09 is enough.\n"
+        "1974-09-09 is enough. Time axis starts and ends at nautical\n"
+        "twilight which is 12deg below horizon.\n"
         "\n"
         "Usage: makeplots sql-datetime [--ra={ra} --dec={dec}]\n";
@@ -129 +130 @@
     const double no_limits   = conf.Get<bool>("no-limits");
 
-    // -12: astronomical twilight
+    // -12: nautical
     ln_rst_time sun_set;   // Sun set with the same date than th provided date
     ln_rst_time sun_rise;  // Sun rise on the following day
@@ -138 +139 @@
     const double mjd = floor(time.Mjd())+49718+0.5;
 
-    cout << "Time: " << time << endl;
+    cout << "Time: " << time     << endl;
+    cout << "Base: " << Time(jd-0.5) << endl;
     cout << "Set:  " << Time(sun_set.set)   << endl;
     cout << "Rise: " << Time(sun_rise.rise) << endl;
@@ -170 +172 @@
     TH1S hframe("", "", 1, (mjd+jd0)*24*3600, (mjd+jd1)*24*3600);
     hframe.SetStats(kFALSE);
-    hframe.GetXaxis()->SetTimeFormat("%Hh%M'");
-    hframe.GetXaxis()->SetTitle("Time [UTC]");
+    hframe.GetXaxis()->SetTimeFormat("%Hh%M%F1995-01-01 00:00:00 GMT");
+    hframe.GetXaxis()->SetTitle((Time(jd).GetAsStr("%d/%m/%Y")+"  -  "+Time(jd+1).GetAsStr("%d/%m/%Y")+"  [UTC]").c_str());
     hframe.GetXaxis()->CenterTitle();
     hframe.GetYaxis()->CenterTitle();
@@ -240 +242 @@
 
         // Create graphs
-        TGraph g1, g2, g3, g4, gm;
+        TGraph g1, g2, g3, g4, gr, gm;
         g1.SetName(name.data());
         g2.SetName(name.data());