Changes between Version 18 and Version 19 of DatabaseBasedAnalysis

Timestamp:

08/04/18 16:45:55 (6 years ago)

Author:

tbretz

Comment:

--

DatabaseBasedAnalysis

@@ -448 +448 @@
 }}}
 
-
-
-
-Some statistics:
-
-
-== Simple Analysis ==
-
-
+You can of course include all the calculations into your query already
+
+{{{
+SELECT
+   Events.*,
+   Angle,
+   weight,
+   @PX      := cosa*X - sina*Y,
+   @PY      := cosa*Y + sina*X,
+   @DX      := MeanX-@PX/1.02,
+   @DY      := MeanY-@PY/1.02,
+   @Norm    := SQRT(@DX*@DX + @DY*@DY),
+   @Dist    := @Norm*0.0117193246260285378 AS Dist,
+   PI()*Width*Length*0.0117193246260285378*0.0117193246260285378 AS Area,
+   @LX      := TRUNCATE((CosDelta*@DY - SinDelta*@DX)/@Norm, 6),
+   @LY      := TRUNCATE((CosDelta*@DX + SinDelta*@DY)/@Norm, 6),
+   @Alpha   := ASIN(@LX) AS Alpha,
+   @Sign    := SIGN(@LY) AS Sign,
+   @M3L     := M3Long*@Sign*0.0117193246260285378,
+   @Slope   := SlopeLong*@Sign/0.0117193246260285378 AS Slope,
+   @Xi      := 1.39 + 0.154*@Slope + 1.679*(1-1/(1+4.86*Leakage1)),
+   @Sign1   := @M3L+0.07,
+   @Sign2   := (@Dist-0.5)*7.2-@Slope,
+   @Disp    := IF (SIGN(@Sign1)<0 || SIGN(@Sign2)<0, -@Xi, @Xi) * (1-Width/Length),
+   @ThetaSq := (@Disp*@Disp + @Dist*@Dist - 2*@Disp*@Dist*SQRT(1-@LX*@LX)) AS ThetaSq
+FROM RunInfo
+LEFT JOIN Events   USING (FileId)
+LEFT JOIN Position USING (FileId, EvtNumber)
+CROSS JOIN
+(
+   SELECT   0 AS Angle UNION ALL
+   SELECT  60 AS Angle UNION ALL
+   SELECT 120 AS Angle UNION ALL
+   SELECT 180 AS Angle UNION ALL
+   SELECT 240 AS Angle UNION ALL
+   SELECT 300 AS Angle
+) Wobble
+WHERE 
+   fSourceKey=5 
+AND 
+   fRunTypeKey=1 
+AND 
+   FileId BETWEEN 131101000 AND 131107000
+AND 
+   fZenithDistanceMax<35
+AND 
+   fR750Ref/fR750Cor>0.9
+}}}
+
+Or you use the existing table for the standard 60° Wobble positions and do just {{{CROSS JOIN Wobble}}}.
+
+This will give you all you need in ''crab.root'' ({{{rootifysql --out=crab.root}}}), but significantly increases computing time and the output file will be about six times larger.
+
+A simple macro just applying all the cuts would then be enough to do a theta-square plot
+
+{{{
+void ganymed3()
+{
+    // Create chain for the tree Result
+    // This is just easier than using TFile/TTree
+    TChain c("Result");
+
+    // Add the input file to the
+    c.AddFile("crab.root");
+
+    // Set some constants (they could be included in the database
+    // in the future)
+    c.SetAlias("mm2deg", "+0.0117193246260285378");
+
+    // Define the cuts
+    c.SetAlias("CutQ", "NumIslands<3.5 && NumUsedPixels>5.5 && Leakage1<0.1");
+
+    c.SetAlias("Cut0",
+               "SlopeSpreadWeighted<0.68 &&"
+               "SlopeSpreadWeighted>0.18 &&"
+               "log10(Area)<(log10(Size)-2)*1.1-1.55 &&"
+               "ConcCore>0.13 &&"
+               "ConcCOG >0.15");
+
+    // Do one plot for each wobble position
+    c.Draw("ThetaSq", "(ThetaSq<1 && CutQ && Cut0 && Angle==0)*( weight)");
+    c.Draw("ThetaSq", "(ThetaSq<1 && CutQ && Cut0 && Angle!=0)*(-weight)", "same");
+}
+}}} 
+