Changeset 19982

Timestamp:

08/10/20 17:58:02 (4 years ago)

Author:

tbretz

Message:

Added histogram titles.

File:

• trunk/FACT++/src/spectrum.cc (modified) (20 diffs)

trunk/FACT++/src/spectrum.cc

@@ -1139 +1139 @@
     WriteHistogram(connection, {
              .name     = "OnTime",
-             .title    = "Effective on time",
+             .title    = "Effective on time per Zd bin",
              .dir      = "Zd",
              .binningx = binning_theta,
@@ -1156 +1156 @@
     WriteHistogram(connection, {
              .name     = "CountN",
-             .title    = "Simulated Zenith Distance",
+             .title    = "Simulated Number of Events in each Zd bin",
              .dir      = "Zd",
              .binningx = binning_theta,
@@ -1173 +1173 @@
     WriteHistogram(connection, {
              .name     = "ZdWeight",
-             .title    = "Zenith Distance Weight",
+             .title    = "Calculated Weight for each Zd bin",
              .dir      = "Zd",
              .binningx = binning_theta,
@@ -1531 +1531 @@
         hist_sim.axisy = "Counts";
 
-        hist_sim.title = "";
+        hist_sim.title = "Events simulated with Corsika";
         hist_sim.v     = "SimCountN";
         hist_sim.err   = "ErrSimCountN";
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Events after cleaning";
         hist_sim.v     = "TrigCountN";
         hist_sim.err   = "ErrTrigCountN";
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Events [true positive] after all cuts (quality, background suppression, theta-square)";
         hist_sim.v     = "SignalN";
         hist_sim.err   = "ErrSignalN";
         WriteHistogram(connection, hist_sim);
 
+        hist_sim.title = "Events wrongly classified as background [false positive] after all Cuts";
+        hist_sim.v     = "BackgroundN";
+        hist_sim.err   = "ErrBackgroundN";
+        WriteHistogram(connection, hist_sim);
+
+        hist_sim.title = "Excess events [true positive - false positive]";
+        hist_sim.v     = "ExcessN";
+        hist_sim.err   = "ErrExcessN";
+        WriteHistogram(connection, hist_sim);
+
 
         hist_sim.axisy = *ib=="theta"?"dN/dE [cm^{-2} s^{-1}]":"dN/dE [cm^{-2} s^{-1} TeV^{-1}]";
 
-        hist_sim.title = "";
+        hist_sim.title = "Weighted simulated events (corsika) per simulated area and time";
         hist_sim.v     = "SimFluxW";
         hist_sim.err   = "ErrSimFluxW";
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Weighted events after trigger (cered) per simulated area and time";
         hist_sim.v     = "TrigFluxW";
         hist_sim.err   = "ErrTrigFluxW";
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Weighted excess events [true positive - false positive]";
         hist_sim.v     = "ExcessFluxW";
         hist_sim.err   = "ErrExcessFluxW";
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Weighted events [true positive] after all cuts (quality, background suppression, theta-square)";
         hist_sim.v     = "SignalFluxW";
         hist_sim.err   = "ErrSignalFluxW";
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Weighted events wrongly classified as background [false positive] after all cuts per simulated area and time";
         hist_sim.v     = "BackgroundFluxW";
         hist_sim.err   = "ErrBackgroundFluxW";
@@ -1575 +1585 @@
 
 
-        hist_sim.title = "";
+        hist_sim.title = "Average simulated energy for weighted events after cuts";
         hist_sim.v     = "AvgEnergySimW";
         hist_sim.err   = "";
@@ -1581 +1591 @@
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Average estimated energy for weighted events after cuts";
         hist_sim.v     = "AvgEnergyEstW";
         hist_sim.err   = "";
@@ -1588 +1598 @@
 
 
-        hist_sim.title = "";
+        hist_sim.title = "Cut effiency times simulated area";
         hist_sim.v     = "EffectiveAreaN";
         hist_sim.err   = "ErrEffectiveAreaN";
@@ -1594 +1604 @@
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Cut efficiency for eighted spectrum times simulated area";
         hist_sim.v     = "EffectiveAreaW";
         hist_sim.err   = "ErrEffectiveAreaW";
@@ -1601 +1611 @@
 
 
-        hist_sim.title = "";
+        hist_sim.title = "Bias of energy estimation for weighted events";
         hist_sim.v     = "BiasW";
         hist_sim.err   = "ErrBiasW";
@@ -1607 +1617 @@
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Resolution of energy estimation for weighted events";
         hist_sim.v     = "ResolutionW";
         hist_sim.err   = "";
@@ -1613 +1623 @@
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Ratio of event after trigger (ceres) to simulated events (corsika)";
         hist_sim.v     = "TriggerEfficiencyN";
         hist_sim.err   = "ErrTriggerEfficiencyN";
@@ -1619 +1629 @@
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Ratio of weighted event after trigger (ceres) to weighted simulated events (corsika)";
         hist_sim.v     = "TriggerEfficiencyW";
         hist_sim.err   = "ErrTriggerEfficiencyW";
@@ -1625 +1635 @@
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Ratio of event after cuts to simulated events (corsika)";
         hist_sim.v     = "CutEfficiencyN";
         hist_sim.err   = "ErrCutEfficiencyN";
@@ -1631 +1641 @@
         WriteHistogram(connection, hist_sim);
 
-        hist_sim.title = "";
+        hist_sim.title = "Ratio of weighted event after cuts to the weighted simulated events (corsika)";
         hist_sim.v     = "CutEfficiencyW";
         hist_sim.err   = "ErrCutEfficiencyW";
@@ -1819 +1829 @@
         hist_est.axisy   = "Counts";
 
-        hist_est.title   = "";
+        hist_est.title   = "Events [true positive] after all cuts (quality, background suppression, theta-square)";
         hist_est.v       = "SignalN";
         hist_est.err     = "ErrSignalN";
         WriteHistogram(connection, hist_est);
 
-        hist_est.title   = "";
+        hist_est.title   = "Events wrongly classified as background [false positive] after all Cuts";
         hist_est.v       = "BackgroundN";
         hist_est.err     = "ErrBackgroundN";
         WriteHistogram(connection, hist_est);
 
-        hist_est.title   = "";
+        hist_est.title   = "Excess events [true positive - false positive]";
         hist_est.v       = "ExcessN";
         hist_est.err     = "ErrExcessN";
@@ -1835 +1845 @@
 
 
-        hist_est.title   = "";
+        hist_est.title   = "Average simulated energy for weighted events after cuts";
         hist_est.v       = "AvgEnergySimW";
         hist_est.err     = "";
@@ -1841 +1851 @@
         WriteHistogram(connection, hist_est);
 
-        hist_est.title   = "";
+        hist_est.title   = "Average estimated energy for weighted events after cuts";
         hist_est.v       = "AvgEnergyEstW";
         hist_est.err     = "";
@@ -1850 +1860 @@
         hist_est.axisy   = "dN/dE [cm^{-2} s^{-1} TeV^{-1}]";
 
-        hist_est.title   = "";
+        hist_est.title   = "Weighted events [true positive] after all cuts (quality, background suppression, theta-square)";
         hist_est.v       = "SignalFluxW";
         hist_est.err     = "ErrSignalFluxW";
         WriteHistogram(connection, hist_est);
 
-        hist_est.title   = "";
+        hist_est.title   = "Weighted events wrongly classified as background [false positive] after all cuts per simulated area and time";
         hist_est.v       = "BackgroundFluxW";
         hist_est.err     = "ErrBackgroundFluxW";
         WriteHistogram(connection, hist_est);
 
-        hist_est.title   = "";
+        hist_est.title   = "Weighted excess events [true positive - false positive]";
         hist_est.v       = "ExcessFluxW";
         hist_est.err     = "ErrExcessFluxW";
@@ -1866 +1876 @@
 
 
-        hist_est.title   = "";
+        hist_est.title   = "Bias of energy estimation for weighted events";
         hist_est.v       = "BiasW";
         hist_est.err     = "ErrBiasW";
@@ -1872 +1882 @@
         WriteHistogram(connection, hist_est);
 
-        hist_est.title   = "";
+        hist_est.title   = "Resolution of energy estimation for weighted events";
         hist_est.v       = "ResolutionW";
         hist_est.err     = "";