Changeset 19915 for trunk/FACT++/spectrum

Timestamp:

12/16/19 20:03:18 (5 years ago)

Author:

tbretz

Message:

Improved comments

Location:

trunk/FACT++/spectrum

Files:

• data.sql (modified) (2 diffs)
• simulation.sql (modified) (3 diffs)
• summary-est.sql (modified) (1 diff)
• summary-sim.sql (modified) (1 diff)

trunk/FACT++/spectrum/data.sql

@@ -1 +1 @@
 SELECT
+   -- Convert variable to bin index
    INTERVAL(Theta, %107:theta)  AS `.theta`,
    INTERVAL(LogEnergyEst, %108:sparse)  AS  `.sparse_est`,
@@ -7 +8 @@
    COUNT(IF(Weight<0, 1, NULL))  AS  `Background`,
 
-   -- Average Energy
+   -- Average Energy: SumEnergyEst/SumW
    SUM(Weight*POW(10, LogEnergyEst))  AS  SumEnergyEst,
    SUM(Weight)                        AS  SumW

trunk/FACT++/spectrum/simulation.sql

@@ -3 +3 @@
    SELECT
       Weight, Energy, LogEnergyEst,
+
+      -- Convert variable to bin index
       INTERVAL(Theta, %107:theta)  AS `.theta`,
       INTERVAL(LogEnergyEst, %108:sparse)  AS `.sparse_est`,
@@ -38 +40 @@
 */
 
-   -- Binning in estimated energy: Signal, Background, Excess
+   -- Binning in estimated energy: Signal, Background
    SUM(  IF(Weight>0,        ZdWeight*SpectralWeight,    0))  AS  SignalW,
    SUM(  IF(Weight<0,        ZdWeight*SpectralWeight,    0))  AS  BackgroundW,
@@ -44 +46 @@
    SUM(  IF(Weight<0, POW(ErrZdWeight*SpectralWeight,2), 0))  AS  BackgroundW2,
 
-   -- Energy Estimation
+   -- Energy Estimation: Bias=ResidualW/SignalW; Resolution=sigma(Bias)
    SUM(IF(Weight>0,     Residual*   ZdWeight*SpectralWeight, 0))  AS  ResidualW,
    SUM(IF(Weight>0, POW(Residual,2)*ZdWeight*SpectralWeight, 0))  AS  ResidualW2,
 
+   -- Average Energy: SumEnergyX/SignalW
    SUM(IF(Weight>0, Energy               *ZdWeight*SpectralWeight, 0))  AS  SumEnergySimW,
    SUM(IF(Weight>0, POW(10, LogEnergyEst)*ZdWeight*SpectralWeight, 0))  AS  SumEnergyEstW

trunk/FACT++/spectrum/summary-est.sql

@@ -11 +11 @@
    SELECT
       *,
-      (POW(10,hi)-POW(10,lo))/1000  AS  Width,
-      Area*ObsTime                  AS  AreaTime
+      (POW(10,hi)-POW(10,lo))/1000  AS  Width,    -- Bin Windth in TeV
+      Area*ObsTime                  AS  AreaTime  -- Total scale dA*dT
    FROM
       BinningEnergy_%100:binning

trunk/FACT++/spectrum/summary-sim.sql

@@ -11 +11 @@
    SELECT
       *,
-      (%103:binwidth)  AS  Width,
-      Area*ObsTime  AS  AreaTime
+      (%103:binwidth)  AS  Width,   -- Bin Windth in TeV or just 1
+      Area*ObsTime  AS  AreaTime    -- Total scale dA*dT
    FROM
       %101:binning