Changeset 19908 for trunk/FACT++
- Timestamp:
- 12/15/19 18:17:23 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/FACT++/src/spectrum.cc
r19907 r19908 2118 2118 for (auto it=types.cbegin(); it!=types.cend(); it++) 2119 2119 { 2120 const bool integral = *ib=="Energy" && !it->empty(); 2121 2120 2122 hist.dir = *ib=="Theta" ? "Data/Theta" : (it->empty() ? "Data/Energy/Differential" : "Data/Energy/Integral"); 2121 2123 2122 2124 hist.axisy = "Counts"; 2123 if ( *ib=="Energy")2125 if (integral) 2124 2126 hist.axisy += " (E>E_{lo})"; 2125 2127 … … 2143 2145 hist.err = ""; 2144 2146 hist.axisy = "#sigma"; 2145 if ( *ib=="Energy")2147 if (integral) 2146 2148 hist.axisy += " (E>E_{lo})"; 2147 2149 WriteHistogram(connection, hist); … … 2151 2153 hist.err = ""; 2152 2154 hist.axisy = "<E_{est}>/GeV"; 2153 if ( *ib=="Energy")2155 if (integral) 2154 2156 hist.axisy += " (E>E_{lo})"; 2155 2157 WriteHistogram(connection, hist); … … 2159 2161 hist.err = "ErrExcessRatio"+*it; 2160 2162 hist.axisy = "Ratio"; 2161 if ( *ib=="Energy")2163 if (integral) 2162 2164 hist.axisy += " (E>E_{lo})"; 2163 2165 WriteHistogram(connection, hist); … … 2165 2167 2166 2168 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}]"; 2171 2170 2172 2171 hist.name = "Spectrum"; … … 2175 2174 WriteHistogram(connection, hist); 2176 2175 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"; 2180 2180 WriteHistogram(connection, hist); 2181 2181 … … 2190 2190 WriteHistogram(connection, hist); 2191 2191 2192 hist.dir = "Data/Energy/Differential"; 2193 hist.name = "IntegratedSpectrum"; 2194 hist.v = "IntegratedFlux"; 2195 hist.err = "ErrIntegratedFlux"; 2196 WriteHistogram(connection, hist); 2197 2192 2198 hist.dir = "Data/Energy/Integral"; 2193 2199 hist.name = "SigmaFlux"; 2194 2200 hist.v = "SigmaFluxI"; 2195 2201 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})"; 2202 2203 WriteHistogram(connection, hist); 2203 2204 }
Note:
See TracChangeset
for help on using the changeset viewer.