Changes between Version 17 and Version 18 of DatabaseBasedAnalysis

Timestamp:

08/04/18 16:36:00 (6 years ago)

Author:

tbretz

Comment:

--

DatabaseBasedAnalysis

@@ -318 +318 @@
 using {{{--list.MyList=Crab.txt}}} as command-line option to rootifysql. Both queries are similar in execution time.
 
+Let's assume the output file is ''crab-data-only.root'' ({{{rootifysql --out=crab-data-only.root}}}). Requesting the data and writing the file took me about 60s.
+
+To run an analysis on the data you can use the following root macro "''ganymed.C''". It produces a theta-square plot. Its execution took about 5s ({{{root ganymed.C++}}})
+
+{{{
+#include <iostream>
+
+#include <TMath.h>
+#include <TH1.h>
+#include <TChain.h>
+#include <TStopwatch.h>
+
+void ganymed()
+{
+    // 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-data-only.root");
+
+    // Define variables for all leaves to be accessed
+    // By definition rootifysql writes only doubles
+    double X, Y, MeanX, MeanY, Width, Length, CosDelta, SinDelta,
+        M3Long, SlopeLong, Leakage1, SlopeSpreadWeighted, Size,
+        ConcCore, ConcCOG, NumIslands, NumUsedPixels;
+
+    // Connect the variables to the cordesponding leaves
+    c.SetBranchAddress("X", &X);
+    c.SetBranchAddress("Y", &Y);
+    c.SetBranchAddress("MeanX", &MeanX);
+    c.SetBranchAddress("MeanY", &MeanY);
+    c.SetBranchAddress("Width", &Width);
+    c.SetBranchAddress("Length", &Length);
+    c.SetBranchAddress("CosDelta", &CosDelta);
+    c.SetBranchAddress("SinDelta", &SinDelta);
+    c.SetBranchAddress("M3Long", &M3Long);
+    c.SetBranchAddress("SlopeLong", &SlopeLong);
+    c.SetBranchAddress("Leakage1", &Leakage1);
+    c.SetBranchAddress("NumIslands", &NumIslands);
+    c.SetBranchAddress("NumUsedPixels", &NumUsedPixels);
+    c.SetBranchAddress("SlopeSpreadWeighted", &SlopeSpreadWeighted);
+    c.SetBranchAddress("Size", &Size);
+    c.SetBranchAddress("ConcCOG", &ConcCOG);
+    c.SetBranchAddress("ConcCore", &ConcCore);
+
+    // Set some constants (they could be included in the database
+    // in the future)
+    double mm2deg = +0.0117193246260285378;
+    double abberation = 1.02;
+
+    // -------------------- Source dependent parameter calculation -------------------
+
+    // Create a histogram for on- and off-data
+    TH1F hon("on",   "", 55, 0, 1);
+    TH1F hoff("off", "", 55, 0, 1);
+
+    // Loop over all events
+    TStopwatch clock;
+    for (int i=0; i<c.GetEntries(); i++)
+    {
+        // read the i-th event from the file
+        c.GetEntry(i);
+
+        // First calculate all cuts to speedup the analysis
+        double area = TMath::Pi()*Width*Length;
+
+        bool cutq = NumIslands<3.5 && NumUsedPixels>5.5 && Leakage1<0.1;
+
+        bool cut0 =
+            SlopeSpreadWeighted<0.68 &&
+            SlopeSpreadWeighted>0.18 &&
+            log10(area*mm2deg*mm2deg)<(log10(Size)-2)*1.1-1.55 &&
+            ConcCore>0.13 &&
+            ConcCOG >0.15;
+
+        if (!cutq || !cut0)
+            continue;
+
+        // Loop over all wobble positions in the camera
+        for (int angle=0; angle<360; angle+=60)
+        {
+            // -------------------- Source dependent parameter calculation -------------------
+
+            double cr = cos(angle*TMath::DegToRad());
+            double sr = sin(angle*TMath::DegToRad());
+
+            double px = cr*X-sr*Y;
+            double py = cr*Y+sr*X;
+
+            double dx = MeanX - px/abberation;
+            double dy = MeanY - py/abberation;
+
+            double dist = sqrt(dx*dx + dy*dy);
+
+            double alpha = asin((CosDelta*dy - SinDelta*dx)/dist);
+            double sgn   = TMath::Sign(1., (CosDelta*dx + SinDelta*dy)/dist);
+
+            // ------------------------------- Application ----------------------------------
+
+            double m3l   = M3Long*sgn*mm2deg;
+            double slope = SlopeLong*sgn/mm2deg;
+
+            // --------------------------------- Aanlysis -----------------------------------
+
+            double xi = 1.39 + 0.154*slope + 1.679*(1-1/(1+4.86*Leakage1));
+
+            double sign1 = m3l+0.07;
+            double sign2 = (dist*mm2deg-0.5)*7.2-slope;
+
+            double disp  = (sign1<0 || sign2<0 ? -xi : xi)*(1-Width/Length)/mm2deg;
+
+            double thetasq = (disp*disp + dist*dist - 2*disp*dist*cos(alpha))*mm2deg*mm2deg;
+
+            // Fill the on- and off-histograms
+            if (angle==0)
+                hon.Fill(thetasq);
+            else
+                hoff.Fill(thetasq, 1./5);
+        }
+    }
+    clock.Print();
+
+    // Plot the result
+    hon.SetMinimum(0);
+    hon.DrawCopy();
+    hoff.DrawCopy("same");
+}
+}}}
+