Changeset 7675

Timestamp:

04/30/06 21:48:49 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mjtrain/MJTrainSeparation.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mjtrain/MJTrainSeparation.cc

@@ -92 +92 @@
     h.SetTitle("Hadronness-Distribution vs. Size");
     res1.SetTitle("Significance Li/Ma");
-    res2.SetTitle("Significance Li/Ma * log_{10}(excess)");
+    res1.SetXTitle("Size [phe]");
+    res1.SetYTitle("Hadronness");
+    res2.SetTitle("Significance-Distribution");
+    res2.SetXTitle("Size-Cut [phe]");
+    res2.SetYTitle("Hadronness-Cut");
+    res1.SetContour(50);
+    res2.SetContour(50);
 
     const Int_t nx = h.GetNbinsX();
@@ -128 +134 @@
 
             res1.SetBinContent(x+1, y+1, sig1);
-            res2.SetBinContent(x+1, y+1, sig2);
         }
 
@@ -134 +139 @@
         gr2.SetPoint(x, h.GetXaxis()->GetBinCenter(x+1), h.GetYaxis()->GetBinCenter(maxy2+1));
 
+        delete hx;
+        delete gx;
+    }
+
+    for (int x=0; x<nx; x++)
+    {
+        TH1 *hx = h.ProjectionY("H_py", x+1);
+        TH1 *gx = g.ProjectionY("G_py", x+1);
+        for (int y=0; y<ny; y++)
+        {
+            const Float_t s = gx->Integral(1, y+1);
+            const Float_t b = hx->Integral(1, y+1);
+            const Float_t sig = MMath::SignificanceLiMa(s+b, b);
+            res2.SetBinContent(x+1, y+1, sig);
+        }
         delete hx;
         delete gx;
@@ -203 +223 @@
     gPad->SetGridy();
     gPad->AddExec("color", "gStyle->SetPalette(1, 0);");
+    res2.SetMaximum(res2.GetMaximum()*1.075);
     res2.DrawCopy("colz");
 }