Changeset 19901 for trunk/FACT++


Ignore:
Timestamp:
12/15/19 11:57:40 (5 years ago)
Author:
tbretz
Message:
Fixed a bug in the calculation of the number of bg events for the integral spectrum, unified the meaning of Background[NW2]
Location:
trunk/FACT++/spectrum
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • trunk/FACT++/spectrum/spectrum.sql

    r19900 r19901  
    1717   SELECT
    1818      *,
    19       SUM(`Signal`)     OVER Integral     AS  SignalI,
    20       SUM(`Background`) OVER Integral /5 AS  BackgroundI,
    21       SUM(SumEnergyEst) OVER Integral     AS  EnergyEstI,
    22       SUM(SumW)         OVER Integral     AS  SumWI
     19      SUM(`Signal`)     OVER Integral  AS  SignalI,
     20      SUM(`Background`) OVER Integral AS  BackgroundI,
     21      SUM(SumEnergyEst) OVER Integral  AS  EnergyEstI,
     22      SUM(SumW)         OVER Integral  AS  SumWI
    2323   FROM
    2424      BinnedData
    … …  
    4444      `Signal` - `Background`         AS  Excess,
    4545      SQRT(`Signal`)                  AS  ErrSignal,
    46       SQRT(`Background`)              AS  ErrBackground,
     46      SQRT(`Background`/5)            AS  ErrBackground,
    4747      ExcErr(`Signal`, `Background`)  AS  ErrExcess,
    4848      LiMa(  `Signal`, `Background`)  AS  Significance,
    … …  
    5050      SignalI - BackgroundI           AS  ExcessI,
    5151      SQRT(SignalI)                   AS  ErrSignalI,
    52       SQRT(BackgroundI)               AS  ErrBackgroundI,
     52      SQRT(BackgroundI/5)             AS  ErrBackgroundI,
    5353      ExcErr(SignalI, BackgroundI)    AS  ErrExcessI,
    5454      LiMa(  SignalI, BackgroundI)    AS  SignificanceI,
    … …  
    123123WHERE
    124124   Excess>0
    125 ORDER BY
    126    %101:bin
    127125)
    128126SELECT
    … …  
    136134FROM
    137135   Spectrum
    138 
    139136WINDOW
    140137   Integral AS (ORDER BY %101:bin DESC)
     138ORDER BY
     139   %101:bin ASC

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

    r19897 r19901  
    2323      `.%100:binning:_est`,
    2424
    25       SUM(SignalN)        AS  SignalN,
    26       SUM(SignalW)        AS  SignalW,
    27       SUM(SignalW2)       AS  SignalW2,
     25      SUM(SignalN)          AS  SignalN,
     26      SUM(SignalW)          AS  SignalW,
     27      SUM(SignalW2)         AS  SignalW2,
    2828
    29       SUM(BackgroundN)    AS  BackgroundN,
    30       SUM(BackgroundW)    AS  BackgroundW,
    31       SUM(BackgroundW2)   AS  BackgroundW2,
     29      SUM(BackgroundN)/5    AS  BackgroundN,
     30      SUM(BackgroundW)/5    AS  BackgroundW,
     31      SUM(BackgroundW2)/25  AS  BackgroundW2,
    3232
    33       SUM(ResidualW)      AS  ResidualW,
    34       SUM(ResidualW2)     AS  ResidualW2,
     33      SUM(ResidualW)        AS  ResidualW,
     34      SUM(ResidualW2)       AS  ResidualW2,
    3535
    36       SUM(SumEnergySimW)  AS  SumEnergySimW,
    37       SUM(SumEnergyEstW)  AS  SumEnergyEstW
     36      SUM(SumEnergySimW)    AS  SumEnergySimW,
     37      SUM(SumEnergyEstW)    AS  SumEnergyEstW
    3838   FROM
    3939      AnalysisMC
    … …  
    4646      *,
    4747
    48       SignalN - BackgroundN/5                           AS  ExcessN,
     48      SignalN - BackgroundN                          AS  ExcessN,
    4949
    50       SQRT(SignalN + BackgroundN/25)                    AS  ErrExcessN,
    51       SQRT(SignalN)                                     AS  ErrSignalN,
    52       SQRT(BackgroundN)                                 AS  ErrBackgroundN,
     50      SQRT(SignalN + BackgroundN)                    AS  ErrExcessN,
     51      SQRT(SignalN)                                  AS  ErrSignalN,
     52      SQRT(BackgroundN)                              AS  ErrBackgroundN,
    5353
    54       SignalW                          /Width/AreaTime  AS  SignalFluxW,
    55       BackgroundW/5                    /Width/AreaTime  AS  BackgroundFluxW,
    56       (SignalW - BackgroundW/5)        /Width/AreaTime  AS  ExcessFluxW,
    57       SQRT(SignalW2)                   /Width/AreaTime  AS  ErrSignalFluxW,
    58       SQRT(BackgroundW2)/5             /Width/AreaTime  AS  ErrBackgroundFluxW,
    59       SQRT(SignalW2 + BackgroundW2/25) /Width/AreaTime  AS  ErrExcessFluxW,
     54      SignalW                       /Width/AreaTime  AS  SignalFluxW,
     55      BackgroundW                   /Width/AreaTime  AS  BackgroundFluxW,
     56      (SignalW - BackgroundW)       /Width/AreaTime  AS  ExcessFluxW,
     57      SQRT(SignalW2)                /Width/AreaTime  AS  ErrSignalFluxW,
     58      SQRT(BackgroundW2)            /Width/AreaTime  AS  ErrBackgroundFluxW,
     59      SQRT(SignalW2 + BackgroundW2) /Width/AreaTime  AS  ErrExcessFluxW,
    6060
    61       SumEnergyEstW/SignalW                             AS  AvgEnergyEstW,
    62       SumEnergySimW/SignalW                             AS  AvgEnergySimW,
     61      SumEnergyEstW/SignalW                          AS  AvgEnergyEstW,
     62      SumEnergySimW/SignalW                          AS  AvgEnergySimW,
    6363
    6464      IF(SignalW>0, ResidualW/SignalW,                                                             NULL)  AS  BiasW,
    … …  
    8484   *,
    8585
    86    IntegralSignalFluxW-IntegralBackgroundFluxW/5           AS  IntegralExcessFluxW,
    87    SQRT(IntegralSignalFluxW2+IntegralBackgroundFluxW2/25)  AS  ErrIntegralExcessFluxW,
    88    SQRT(IntegralSignalFluxW2)                              AS  ErrIntegralSignalFluxW,
    89    SQRT(IntegralBackgroundFluxW2)                          AS  ErrIntegralBackgroundFluxW,
     86   IntegralSignalFluxW-IntegralBackgroundFluxW          AS  IntegralExcessFluxW,
     87   SQRT(IntegralSignalFluxW2+IntegralBackgroundFluxW2)  AS  ErrIntegralExcessFluxW,
     88   SQRT(IntegralSignalFluxW2)                           AS  ErrIntegralSignalFluxW,
     89   SQRT(IntegralBackgroundFluxW2)                       AS  ErrIntegralBackgroundFluxW,
    9090
    91    IntegralEnergyEstW/IntegralSignalW                      AS  AvgIntegralEnergyEstW,
    92    IntegralEnergySimW/IntegralSignalW                      AS  AvgIntegralEnergySimW
     91   IntegralEnergyEstW/IntegralSignalW                   AS  AvgIntegralEnergyEstW,
     92   IntegralEnergySimW/IntegralSignalW                   AS  AvgIntegralEnergySimW
    9393
    9494FROM

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

    r19897 r19901  
    4747      %102:bin,
    4848
    49       SUM(SignalN)        AS  SignalN,
    50       SUM(SignalW)        AS  SignalW,
    51       SUM(SignalW2)       AS  SignalW2,
     49      SUM(SignalN)          AS  SignalN,
     50      SUM(SignalW)          AS  SignalW,
     51      SUM(SignalW2)         AS  SignalW2,
    5252
    53       SUM(BackgroundN)    AS  BackgroundN,
    54       SUM(BackgroundW)    AS  BackgroundW,
    55       SUM(BackgroundW2)   AS  BackgroundW2,
     53      SUM(BackgroundN)/5    AS  BackgroundN,
     54      SUM(BackgroundW)/5    AS  BackgroundW,
     55      SUM(BackgroundW2)/25  AS  BackgroundW2,
    5656
    57       SUM(ResidualW)      AS  ResidualW,
    58       SUM(ResidualW2)     AS  ResidualW2,
     57      SUM(ResidualW)        AS  ResidualW,
     58      SUM(ResidualW2)       AS  ResidualW2,
    5959
    60       SUM(SumEnergySimW)  AS  SumEnergyEstW,
    61       SUM(SumEnergyEstW)  AS  SumEnergySimW
     60      SUM(SumEnergySimW)    AS  SumEnergyEstW,
     61      SUM(SumEnergyEstW)    AS  SumEnergySimW
    6262   FROM
    6363      AnalysisMC
    … …  
    6969   SELECT
    7070      *,
    71       TrigCountN/SimCountN        AS  TriggerEfficiencyN,
    72       TrigSumW/SimSumW            AS  TriggerEfficiencyW,
     71      TrigCountN/SimCountN     AS  TriggerEfficiencyN,
     72      TrigSumW/SimSumW         AS  TriggerEfficiencyW,
    7373
    74       SignalW  - BackgroundW/5    AS  ExcessW,
    75       SignalW2 + BackgroundW2/25  AS  ExcessW2
     74      SignalW  - BackgroundW   AS  ExcessW,
     75      SignalW2 + BackgroundW2  AS  ExcessW2
    7676   FROM
    7777      Simulated
    … …  
    8888      *,
    8989
    90       SQRT(SimCountN)                       AS  ErrSimCountN,
    91       SQRT(TrigCountN)                      AS  ErrTrigCountN,
     90      SQRT(SimCountN)                        AS  ErrSimCountN,
     91      SQRT(TrigCountN)                       AS  ErrTrigCountN,
    9292      SQRT(SignalN)                          AS  ErrSignalN,
    9393
    … …  
    9595      TrigSumW             /Width/AreaTime   AS  TrigFluxW,
    9696      SignalW              /Width/AreaTime   AS  SignalFluxW,
    97       BackgroundW/5        /Width/AreaTime   AS  BackgroundFluxW,
     97      BackgroundW          /Width/AreaTime   AS  BackgroundFluxW,
    9898      ExcessW              /Width/AreaTime   AS  ExcessFluxW,
     99
    99100      SQRT(SimSumW2)       /Width/AreaTime   AS  ErrSimFluxW,
    100101      SQRT(TrigSumW2)      /Width/AreaTime   AS  ErrTrigFluxW,
    101102      SQRT(SignalW2)       /Width/AreaTime   AS  ErrSignalFluxW,
    102       SQRT(BackgroundW2)/5 /Width/AreaTime   AS  ErrBackgroundFluxW,
     103      SQRT(BackgroundW2)  /Width/AreaTime   AS  ErrBackgroundFluxW,
    103104      SQRT(ExcessW2)       /Width/AreaTime   AS  ErrExcessFluxW,
    104105
    … …  
    144145   CutEfficiencyW*SQRT(ExcessW2/POW(ExcessW,2) + TrigSumW2/POW(TrigSumW,2))  AS  ErrCutEfficiencyW,
    145146
    146    IntegralSignalFluxW - IntegralBackgroundFluxW/5         AS  IntegralExcessFluxW,
    147    SQRT(IntegralSignalFluxW2+IntegralBackgroundFluxW2/25)  AS  ErrIntegralExcessFluxW,
    148    SQRT(IntegralSignalFluxW2)                              AS  ErrIntegralSignalFluxW,
    149    SQRT(IntegralBackgroundFluxW2)                          AS  ErrIntegralBackgroundFluxW,
    150    SQRT(IntegralSimFluxW2)                                 AS  ErrIntegralSimFluxW,
     147   IntegralSignalFluxW - IntegralBackgroundFluxW        AS  IntegralExcessFluxW,
     148   SQRT(IntegralSignalFluxW2+IntegralBackgroundFluxW2)  AS  ErrIntegralExcessFluxW,
     149   SQRT(IntegralSignalFluxW2)                           AS  ErrIntegralSignalFluxW,
     150   SQRT(IntegralBackgroundFluxW2)                       AS  ErrIntegralBackgroundFluxW,
     151   SQRT(IntegralSimFluxW2)                              AS  ErrIntegralSimFluxW,
    151152
    152    IntegralEnergyEstW/IntegralSignalW                      AS  AvgIntegralEnergyEstW,
    153    IntegralEnergySimW/IntegralSignalW                      AS  AvgIntegralEnergySimW
     153   IntegralEnergyEstW/IntegralSignalW                   AS  AvgIntegralEnergyEstW,
     154   IntegralEnergySimW/IntegralSignalW                   AS  AvgIntegralEnergySimW
    154155
    155156FROM
Note: See TracChangeset for help on using the changeset viewer.