Changeset 19908


Ignore:
Timestamp:
Dec 15, 2019, 6:17:23 PM (4 months ago)
Author:
tbretz
Message:
Fixed some axis labels for the spectrum histograms.
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/FACT++/src/spectrum.cc

    r19907 r19908  
    21182118            for (auto it=types.cbegin(); it!=types.cend(); it++)
    21192119            {
     2120                const bool integral = *ib=="Energy" && !it->empty();
     2121
    21202122                hist.dir = *ib=="Theta" ? "Data/Theta" : (it->empty() ? "Data/Energy/Differential" : "Data/Energy/Integral");
    21212123
    21222124                hist.axisy = "Counts";
    2123                 if (*ib=="Energy")
     2125                if (integral)
    21242126                    hist.axisy += " (E>E_{lo})";
    21252127
     
    21432145                hist.err   = "";
    21442146                hist.axisy = "#sigma";
    2145                 if (*ib=="Energy")
     2147                if (integral)
    21462148                    hist.axisy += " (E>E_{lo})";
    21472149                WriteHistogram(connection, hist);
     
    21512153                hist.err   = "";
    21522154                hist.axisy = "<E_{est}>/GeV";
    2153                 if (*ib=="Energy")
     2155                if (integral)
    21542156                    hist.axisy += " (E>E_{lo})";
    21552157                WriteHistogram(connection, hist);
     
    21592161                hist.err   = "ErrExcessRatio"+*it;
    21602162                hist.axisy = "Ratio";
    2161                 if (*ib=="Energy")
     2163                if (integral)
    21622164                    hist.axisy += " (E>E_{lo})";
    21632165                WriteHistogram(connection, hist);
     
    21652167
    21662168            hist.dir = *ib=="Theta" ? "Data/Theta" : "Data/Energy/Differential";
    2167             hist.axisy = "dN/dE ";
    2168             if (*ib=="Energy")
    2169                 hist.axisy += "(E>E_{lo}) ";
    2170             hist.axisy += *ib=="Theta" ? "[cm^{-2} s^{-1}]" : "[cm^{-2} s^{-1} TeV^{-1}]";
     2169            hist.axisy = *ib=="Theta" ? "dN/dE [cm^{-2} s^{-1}]" : "dN/dE [cm^{-2} s^{-1} TeV^{-1}]";
    21712170
    21722171            hist.name = "Spectrum";
     
    21752174            WriteHistogram(connection, hist);
    21762175
    2177             hist.name = "SigmaFlux";
    2178             hist.v    = "SigmaFlux";
    2179             hist.err  = "";
     2176            hist.name  = "SigmaFlux";
     2177            hist.v     = "SigmaFlux";
     2178            hist.err   = "";
     2179            hist.axisy = "Relative standard deviations";
    21802180            WriteHistogram(connection, hist);
    21812181
     
    21902190                WriteHistogram(connection, hist);
    21912191
     2192                hist.dir   = "Data/Energy/Differential";
     2193                hist.name  = "IntegratedSpectrum";
     2194                hist.v     = "IntegratedFlux";
     2195                hist.err   = "ErrIntegratedFlux";
     2196                WriteHistogram(connection, hist);
     2197
    21922198                hist.dir   = "Data/Energy/Integral";
    21932199                hist.name  = "SigmaFlux";
    21942200                hist.v     = "SigmaFluxI";
    21952201                hist.err   = "";
    2196                 WriteHistogram(connection, hist);
    2197 
    2198                 hist.dir   = "Data/Energy/Differential";
    2199                 hist.name  = "IntegratedSpectrum";
    2200                 hist.v     = "IntegratedFlux";
    2201                 hist.err   = "ErrIntegratedFlux";
     2202                hist.axisy = "Relative standard deviations (E>E_{lo})";
    22022203                WriteHistogram(connection, hist);
    22032204            }
Note: See TracChangeset for help on using the changeset viewer.