Changes between Version 254 and Version 255 of DatabaseBasedAnalysis/Spectrum

Timestamp:

12/10/19 14:47:05 (5 years ago)

Author:

tbretz

Comment:

--

DatabaseBasedAnalysis/Spectrum

@@ -272 +272 @@
 `%0:bins` is a placeholder for the bin boundaries, e.g. `5, 10, 15, 20, 25, 30` (five bins between 5° and 30° plus underflow and overflow). It should be identical to the binning used for data files in !DataFiles.
 
+=== Analysis Query ===
+
+The analysis of Data and !MontaCarlo files must be done totally identical to produce reasonable results. Therefore, the exact same query (or code) should be used for the analysis of both. In this example, the query has to provide two columns, `Weight` and `LogEnergyEst` for all gamma-line events. The weight must be +1 for events in the on-region and -0.2 for an event in the off-region (corresponding to the number of five wobble positions). !LogEnergyEst must contain \(\log_{10}\) of the estimated energy for each event. Only events surviving background suppression and spatial (theta) cuts should be considered.
+
+{{{#!Spoiler
+{{{#!sql
+WITH Table0 AS
+(
+   SELECT
+      %0:columns  -- this could be removed if we can join events via the same columns (without CorsikaNumResuse)
+      Weight,
+      Size,
+      NumUsedPixels,
+      NumIslands,
+      Leakage1,
+      MeanX,
+      MeanY,
+      CosDelta,
+      SinDelta,
+      M3Long,
+      SlopeLong,
+      Width/Length      AS WdivL,
+      PI()*Width*Length AS Area,
+      cosa*X - sina*Y   AS PX,
+      cosa*Y + sina*X   AS PY
+   FROM
+      %1:files
+   LEFT JOIN
+      %2:runinfo USING (FileId)
+   LEFT JOIN 
+      %3:events USING (FileId)  -- This could be replaced by a user uploaded temporary table
+   LEFT JOIN 
+      %4:positions USING (FileId, EvtNumber)
+   CROSS JOIN 
+      Wobble
+   WHERE
+      NumUsedPixels>5.5
+   AND
+      NumIslands<3.5
+   AND
+      Leakage1<0.1
+),
+
+Table1 AS
+(
+   SELECT
+      %0:columns
+      Weight,
+      Size, CosDelta, SinDelta, M3Long, SlopeLong, Leakage1, WdivL,
+      MeanX - PX/1.02e0 AS DX,
+      MeanY - PY/1.02e0 AS DY
+   FROM
+      Table0
+   WHERE
+      Area < LOG10(Size)*898e0 - 1535e0
+),
+
+Table2 AS
+(
+   SELECT
+      %0:columns
+      Weight,
+      Size, CosDelta, SinDelta, DX, DY, M3Long, SlopeLong, Leakage1, WdivL,
+      SQRT(DX*DX + DY*DY) AS Norm
+   FROM
+      Table1
+),
+
+Table3 AS
+(
+   SELECT
+      %0:columns
+      Weight,
+      Size, M3Long, SlopeLong, Leakage1, WdivL, Norm,
+      LEAST(GREATEST((CosDelta*DY - SinDelta*DX)/Norm, -1), 1) AS LX,
+      SIGN(CosDelta*DX + SinDelta*DY) AS Sign
+   FROM
+      Table2
+),
+
+Table5 AS
+(
+   SELECT
+      %0:columns
+      Weight,
+      Size, Leakage1, WdivL, LX,
+      Norm          *0.0117193246260285378e0 AS Dist,
+      M3Long   *Sign*0.0117193246260285378e0 AS M3L,
+      SlopeLong*Sign/0.0117193246260285378e0 AS Slope
+   FROM
+      Table3
+),
+
+Table6 AS
+(
+   SELECT
+      %0:columns
+      Weight,
+      Size, WdivL, Dist, LX, M3L, Slope,
+      1.39252e0 + 0.154247e0*Slope + 1.67972e0*(1-1/(1+4.86232e0*Leakage1)) AS Xi
+   FROM
+      Table5
+),
+
+Table7 AS
+(
+   SELECT
+      %0:columns
+      Weight,
+      Size, Dist, LX,
+      IF (M3L<-0.07 OR (Dist-0.5e0)*7.2e0-Slope<0, -Xi, Xi) * (1-WdivL) AS Disp
+   FROM
+      Table6
+)
+
+SELECT
+   %0:columns
+   Weight,
+   (Disp*Disp + Dist*Dist - 2*Disp*Dist*SQRT(1-LX*LX)) AS ThetaSq,
+   %5:estimator AS LogEnergyEst
+FROM
+   Table7
+HAVING
+   ThetaSq<0.024
+}}}
+}}}
+
+The placeholder `%0:column` is currently used to request !MonteCarlo true values (such as Energy) in addition to the anaylsis values only for the analysis of simulated data. The must be no comma behind! In additon,  `%1:files` is a placeholder for the table containing the `FileId`s to analyze, `%2:runinfo` for the table with the run info data, `%3:events` for the table with the image parameters and `%4:positions` for the table with the source positions in the camera.
 
 === Analyze Data ===