Changes between Version 237 and Version 238 of DatabaseBasedAnalysis/Spectrum

Timestamp:

12/09/19 18:28:20 (5 years ago)

Author:

tbretz

Comment:

--

DatabaseBasedAnalysis/Spectrum

@@ -194 +194 @@
 == Get Zenith Angle Histogram ==
 
+{{{!#Spoiler
+{{{!#sql
+CREATE TEMPORARY TABLE ThetaHist
+(
+   `.theta`    SMALLINT UNSIGNED NOT NULL,
+   lo          DOUBLE            NOT NULL COMMENT 'Lower edge of zenith distance bin in degree',
+   hi          DOUBLE            NOT NULL COMMENT 'Upper edge of zenith distance bin in degree',
+   CountN      INT UNSIGNED      NOT NULL,
+   OnTime      FLOAT             NOT NULL,
+   ZdWeight    DOUBLE            NOT NULL COMMENT 'tau(delta theta)',
+   ErrZdWeight DOUBLE            NOT NULL COMMENT 'sigma(tau)',
+   PRIMARY KEY (`.theta`) USING HASH
+) ENGINE=Memory
+AS
+(
+   WITH EventCount AS
+   (
+      SELECT
+         INTERVAL(DEGREES(Theta), %0:bins) AS `.theta`,
+         COUNT(*) AS CountN
+      FROM
+         MonteCarloFiles
+      LEFT JOIN
+         factmc.OriginalMC USING(FileId)
+      GROUP BY
+         `.theta`
+   )
+   SELECT
+      `.theta`, lo, hi,
+      CountN,
+      OnTime,
+      OnTime/CountN AS ZdWeight,
+      (OnTime/CountN)*SQRT(POW(1/300, 2) + 1/CountN) AS ErrZdWeight\
+   FROM
+      ObservationTime
+   LEFT JOIN
+      EventCount USING(`.theta`)
+   LEFT JOIN
+      BinningTheta ON `.theta`=bin
+   ORDER BY
+      `.theta`
+)
+}}}
+}}}
+
 == Analyze Data ==