Changeset 19971 for trunk/FACT++/spectrum

Timestamp:

07/20/20 16:20:59 (4 years ago)

Author:

tbretz

Message:

Added a binning in Impact Parameter for crosschecks

Location:

trunk/FACT++/spectrum

Files:

• display.C (modified) (1 diff)
• simulation.sql (modified) (3 diffs)

trunk/FACT++/spectrum/display.C

@@ -91 +91 @@
     h1->SetMarkerStyle(23);
     h1->DrawCopy("P same");
+
+    // --------------------------------
+
+    c = new TCanvas("Impact", "Impact");
+    c->Divide(2,2);
+    c->cd(1);
+    file.GetObject("MC/dense/TrueEnergy/Impact", h2);
+    h2->DrawCopy("colz");
+    h2->DrawCopy("same");
+    c->cd(2);
+    file.GetObject("MC/sparse/TrueEnergy/Impact", h2);
+    h2->DrawCopy("colz");
+    h2->DrawCopy("same");
+    c->cd(3);
+    file.GetObject("MC/theta/Impact", h2);
+    h2->DrawCopy("colz");
+    h2->DrawCopy("same");
 
     // --------------------------------

trunk/FACT++/spectrum/simulation.sql

@@ -10 +10 @@
       INTERVAL(LOG10(Energy), %108:sparse)  AS `.sparse_sim`,
       INTERVAL(LOG10(Energy), %109:dense)  AS `.dense_sim`,
+      INTERVAL(Impact/100, %110:impact)  AS `.impact`,
 
-      (%110:spectrum)/POW(Energy, SpectralIndex) AS SpectralWeight,  -- FIXME: Is this correct for files with different Slopes?
+      (%111:spectrum)/POW(Energy, SpectralIndex) AS SpectralWeight,  -- FIXME: Is this correct for files with different Slopes?
       LogEnergyEst - log10(Energy) AS Residual
    FROM
@@ -27 +28 @@
    `.dense_est`,
    `.dense_sim`,
+   `.impact`,
 
    -- Without any weight applied
@@ -59 +61 @@
    ThetaDist USING(`.theta`)
 GROUP BY
-   `.theta`, `.sparse_est`, `.sparse_sim`, `.dense_est`, `.dense_sim`
+   `.theta`, `.sparse_est`, `.sparse_sim`, `.dense_est`, `.dense_sim`, `.impact`