Changes between Version 22 and Version 23 of DatabaseBasedAnalysis/Examples

Timestamp:

08/16/18 20:33:20 (7 years ago)

Author:

tbretz

Comment:

--

DatabaseBasedAnalysis/Examples

@@ -574 +574 @@
 Note that this query can easily run 30min or more!
 
-== Optim Disp Plot1 ==
+== Optim Disp (Leakage==0) ==
 
 {{{#!Spoiler
@@ -716 +716 @@
 
 [[Image(Result.png)]]
+
+== Optim Disp (Leakage>0) ==
+
+{{{#!Spoiler
+{{{#!cpp
+#include <iostream>
+#include <iomanip>
+
+#include <TMath.h>
+#include <TH2.h>
+#include <TChain.h>
+#include <TCanvas.h>
+
+void optimdisp3()
+{
+    // 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("simulation.root");
+
+    // Define variables for all leaves to be accessed
+    // By definition rootifysql writes only doubles
+    double FileId, EvtNumber, X, Y, MeanX, MeanY, Width, Length, CosDelta, SinDelta,
+        M3Long, SlopeLong, Leakage1, SlopeSpreadWeighted, Size,
+        ConcCore, ConcCOG, NumIslands, NumUsedPixels, Zd;
+
+    // Connect the variables to the cordesponding leaves
+    c.SetBranchAddress("FileId", &FileId);
+    c.SetBranchAddress("EvtNumber", &EvtNumber);
+    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("Zd", &Zd);
+    //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 -------------------
+
+    TH2F hist("", "",
+              32, 0.800-0.0250, 2.400-0.0250,
+              60, 3.000-0.0250, 6.000-0.0250);
+
+
+    for (float l0=0.8; l0<2.4; l0+=0.05)
+        for (float l1=3; l1<6; l1+=0.05)
+        {
+            int cnt_on = 0;
+
+            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;
+
+                // The abberation correction does increase also Width and Length by 1.02
+
+                bool cutq = NumIslands<3.5 && NumUsedPixels>5.5 && Leakage1>0;
+                if (!cutq)
+                    continue;
+
+                bool cut0 = area < log10(Size)*898-1535;
+                if (!cut0)
+                    continue;
+
+                int angle = 0;
+
+                // -------------------- 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*1.02;
+                double dy = MeanY - py*1.02;
+
+                double norm = sqrt(dx*dx + dy*dy);
+                double dist = norm*mm2deg;
+
+                double lx = min(max((CosDelta*dy - SinDelta*dx)/norm, -1.), 1.);
+                double ly = min(max((CosDelta*dx + SinDelta*dy)/norm, -1.), 1.);
+
+                double alpha = asin(lx);
+                double sgn   = TMath::Sign(1., ly);
+
+                // ------------------------------- Application ----------------------------------
+
+                double m3l   = M3Long*sgn*mm2deg;
+                double slope = SlopeLong*sgn/mm2deg;
+
+                // --------------------------------- Analysis -----------------------------------
+
+                double xi = 1.36 + 0.085*slope + l0 *(1-1/(1+l1*Leakage1));
+
+                // Optim: 1.36 + 0.085
+
+                double sign1 = m3l+0.07;
+                double sign2 = (dist-0.5)*7.2-slope;
+
+                double disp  = (sign1<0 || sign2<0 ? -xi : xi)*(1-Width/Length);
+
+                double thetasq = disp*disp + dist*dist - 2*disp*dist*sqrt(1-lx*lx);
+
+                if (thetasq<0.024)
+                    cnt_on++;
+            }
+
+            hist.Fill(l0, l1, cnt_on);
+
+            cout << l0 << " " << l1 << " " << cnt_on << endl;
+        }
+
+    TCanvas *canv = new TCanvas;
+    canv->SetTopMargin(0.01);
+
+    hist.SetStats(kFALSE);
+    hist.SetContour(100);
+    hist.SetXTitle("L_{0} / deg");
+    hist.SetYTitle("L_{1} / deg");
+    hist.GetYaxis()->SetTitleOffset(1.2);
+    hist.DrawCopy("colz cont2");
+
+}
+}}}
+}}}
+
+Results in
+
+[[Image(Result2.png)]]
 
 == Energy Optimization ==