Changes between Version 2 and Version 3 of DatabaseBasedAnalysis/RandomForest

Timestamp:

08/15/18 17:35:15 (7 years ago)

Author:

tbretz

Comment:

--

DatabaseBasedAnalysis/RandomForest

@@ -240 +240 @@
 }}}
 
-
-
+== Plotting Result == 
+
+{{{#!Spoiler
+{{{#!cpp
+#include <fstream>
+#include <TH2.h>
+#include <TF1.h>
+#include <TGraph.h>
+#include <TCanvas.h>
+
+void readsim()
+{
+    TH2F h2svse ("H_SvsE",      "", 26*3, 2.2, 4.8, 3*30, 1.3, 4.3);
+    TH2F h2est  ("H_EstVsMC",   "", 26*3, 2.2, 4.8, 26*3, 2.2, 4.8);
+    TH2F h2bias ("H_BiasLog",   "", 26*3, 2.2, 4.8, 100, -1, 1);
+    TH2F h2lin  ("H_BiasLin",   "", 26*3, 2.2, 4.8, 100, -1, 2);
+
+    TH2F h2size ("H_ResSize",   "", 26*3, 2.2, 4.8, 100, -1, 1);
+    TH2F h2dist ("H_ResDist",   "", 50,  0, 2.5,    100, -1, 1);
+    TH2F h2slope("H_ResSlope",  "", 50, -10, 10,    100, -1, 1);
+    TH2F h2zd   ("H_ResZd",     "", 90, 0,   90,    100, -1, 1);
+
+    ifstream fin0("sim-test-cuts.csv");
+    ifstream fin1("sim-test-results-cuts.csv");
+
+    TString str;
+    str.ReadLine(fin0);
+    str.ReadLine(fin1);
+
+    while (1)
+    {
+        if (!fin0 || !fin1)
+            break;
+
+        double Emc, Size, Zd, Dist, Disp, Slope, M3L, Leakage, Width, Length, Eest;
+        fin0 >> Emc >> Size >> Zd >> Dist >> Disp >> Slope >> M3L >> Leakage >> Width >> Length;
+        fin1 >> Eest;
+
+        h2svse.Fill(Emc,  Size);
+        h2est.Fill( Emc,  Eest);
+        h2bias.Fill(Eest, Eest-Emc);
+        h2lin.Fill( Eest, (pow(10, Eest)-pow(10, Emc))/pow(10, Emc));
+
+        h2size.Fill( Size,  Eest-Emc);
+        h2dist.Fill( Dist,  Eest-Emc);
+        h2slope.Fill(Slope, Eest-Emc);
+        h2zd.Fill(   Zd,    Eest-Emc);
+    }
+
+
+    TF1 fx("f", "x", -100, 100);
+    TF1 f0("f", "0", -100, 100);
+
+    TCanvas *canv = new TCanvas("Canvas", "Energy estimation", 800, 900);
+    canv->Divide(2,4);
+
+    canv->cd(1);
+    gPad->SetTopMargin(0.005);
+    gPad->SetGridy();
+    h2size.SetStats(kFALSE);
+    h2size.SetXTitle("lg(E_{mc}/GeV)");
+    h2size.SetYTitle("lg(E_{est}/GeV)-lg(E_{mc}/GeV)");
+    h2size.DrawCopy("colz");
+    f0.DrawCopy("same");
+
+    canv->cd(3);
+    gPad->SetTopMargin(0.005);
+    gPad->SetGridy();
+    h2dist.SetStats(kFALSE);
+    h2dist.SetXTitle("Dist / degree");
+    h2dist.SetYTitle("lg(E_{est}/GeV)-lg(E_{mc}/GeV)");
+    h2dist.DrawCopy("colz");
+    f0.DrawCopy("same");
+
+    canv->cd(5);
+    gPad->SetTopMargin(0.005);
+    gPad->SetGridy();
+    h2slope.SetStats(kFALSE);
+    h2slope.SetXTitle("Slope / ns/degree");
+    h2slope.SetYTitle("lg(E_{est}/GeV)-lg(E_{mc}/GeV)");
+    h2slope.DrawCopy("colz");
+    f0.DrawCopy("same");
+
+    canv->cd(7);
+    gPad->SetTopMargin(0.005);
+    gPad->SetGridy();
+    h2zd.SetStats(kFALSE);
+    h2zd.SetXTitle("Zenith Distance / degree");
+    h2zd.SetYTitle("lg(E_{est}/GeV)-lg(E_{mc}/GeV)");
+    h2zd.DrawCopy("colz");
+    f0.DrawCopy("same");
+
+    canv->cd(2);
+    gPad->SetTopMargin(0.005);
+    gPad->SetGridy();
+    h2svse.SetStats(kFALSE);
+    h2svse.SetXTitle("lg(Size)");
+    h2svse.SetYTitle("lg(E_{mc}/GeV)");
+    h2svse.DrawCopy("colz");
+    fx.DrawCopy("same");
+
+    canv->cd(4);
+    gPad->SetTopMargin(0.005);
+    gPad->SetGridy();
+    h2est.SetStats(kFALSE);
+    h2est.SetXTitle("lg(E_{mc})");
+    h2est.SetYTitle("lg(E_{est}/GeV)");
+    h2est.DrawCopy("colz");
+    fx.DrawCopy("same");
+
+    canv->cd(6);
+    gPad->SetTopMargin(0.005);
+    gPad->SetGridy();
+    h2bias.SetStats(kFALSE);
+    h2bias.SetXTitle("lg(E_{est}/GeV)");
+    h2bias.SetYTitle("lg(E_{est}/GeV)-lg(E_{mc}/GeV)");
+    h2bias.DrawCopy("colz");
+    f0.DrawCopy("same");
+
+
+    canv->cd(8);
+    gPad->SetTopMargin(0.005);
+    gPad->SetGridy();
+
+    TGraph grlin;
+    TGraph grlog;
+
+    for (int x=1; x<h2bias.GetNbinsX(); x++)
+    {
+        TH1D *p = h2bias.ProjectionY("_py", x, x);
+        if (p->GetRMS()>0)
+            grlog.SetPoint(grlog.GetN(), h2bias.GetXaxis()->GetBinCenter(x), pow(10, p->GetRMS())-1);
+        delete p;
+
+        p = h2lin.ProjectionY("_py", x, x);
+        if (p->GetRMS()>0)
+            grlin.SetPoint(grlin.GetN(), h2lin.GetXaxis()->GetBinCenter(x), p->GetRMS());
+        delete p;
+    }
+
+    grlog.SetMarkerColor(kBlue);
+    grlin.SetMinimum(0);
+    grlin.GetHistogram()->SetXTitle("log(E_{est}/GeV)");
+    grlin.GetHistogram()->SetYTitle("RMS");
+    grlin.DrawClone("A*");
+    grlog.DrawClone("*");
+}
+}}}
+}}}
+
+