Index: /trunk/FACT++/src/spectrum.cc
===================================================================
--- /trunk/FACT++/src/spectrum.cc	(revision 19907)
+++ /trunk/FACT++/src/spectrum.cc	(revision 19908)
@@ -2118,8 +2118,10 @@
             for (auto it=types.cbegin(); it!=types.cend(); it++)
             {
+                const bool integral = *ib=="Energy" && !it->empty();
+
                 hist.dir = *ib=="Theta" ? "Data/Theta" : (it->empty() ? "Data/Energy/Differential" : "Data/Energy/Integral");
 
                 hist.axisy = "Counts";
-                if (*ib=="Energy")
+                if (integral)
                     hist.axisy += " (E>E_{lo})";
 
@@ -2143,5 +2145,5 @@
                 hist.err   = "";
                 hist.axisy = "#sigma";
-                if (*ib=="Energy")
+                if (integral)
                     hist.axisy += " (E>E_{lo})";
                 WriteHistogram(connection, hist);
@@ -2151,5 +2153,5 @@
                 hist.err   = "";
                 hist.axisy = "<E_{est}>/GeV";
-                if (*ib=="Energy")
+                if (integral)
                     hist.axisy += " (E>E_{lo})";
                 WriteHistogram(connection, hist);
@@ -2159,5 +2161,5 @@
                 hist.err   = "ErrExcessRatio"+*it;
                 hist.axisy = "Ratio";
-                if (*ib=="Energy")
+                if (integral)
                     hist.axisy += " (E>E_{lo})";
                 WriteHistogram(connection, hist);
@@ -2165,8 +2167,5 @@
 
             hist.dir = *ib=="Theta" ? "Data/Theta" : "Data/Energy/Differential";
-            hist.axisy = "dN/dE ";
-            if (*ib=="Energy")
-                hist.axisy += "(E>E_{lo}) ";
-            hist.axisy += *ib=="Theta" ? "[cm^{-2} s^{-1}]" : "[cm^{-2} s^{-1} TeV^{-1}]";
+            hist.axisy = *ib=="Theta" ? "dN/dE [cm^{-2} s^{-1}]" : "dN/dE [cm^{-2} s^{-1} TeV^{-1}]";
 
             hist.name = "Spectrum";
@@ -2175,7 +2174,8 @@
             WriteHistogram(connection, hist);
 
-            hist.name = "SigmaFlux";
-            hist.v    = "SigmaFlux";
-            hist.err  = "";
+            hist.name  = "SigmaFlux";
+            hist.v     = "SigmaFlux";
+            hist.err   = "";
+            hist.axisy = "Relative standard deviations";
             WriteHistogram(connection, hist);
 
@@ -2190,14 +2190,15 @@
                 WriteHistogram(connection, hist);
 
+                hist.dir   = "Data/Energy/Differential";
+                hist.name  = "IntegratedSpectrum";
+                hist.v     = "IntegratedFlux";
+                hist.err   = "ErrIntegratedFlux";
+                WriteHistogram(connection, hist);
+
                 hist.dir   = "Data/Energy/Integral";
                 hist.name  = "SigmaFlux";
                 hist.v     = "SigmaFluxI";
                 hist.err   = "";
-                WriteHistogram(connection, hist);
-
-                hist.dir   = "Data/Energy/Differential";
-                hist.name  = "IntegratedSpectrum";
-                hist.v     = "IntegratedFlux";
-                hist.err   = "ErrIntegratedFlux";
+                hist.axisy = "Relative standard deviations (E>E_{lo})";
                 WriteHistogram(connection, hist);
             }