Changeset 6489 for trunk/MagicSoft/Mars/mimage/MHImagePar.cc

Timestamp:

02/15/05 16:18:28 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mimage/MHImagePar.cc (modified) (8 diffs)

trunk/MagicSoft/Mars/mimage/MHImagePar.cc

@@ -66 +66 @@
     fHistSatHi.SetDirectory(NULL);
     fHistSatHi.UseCurrentStyle();
+    fHistSatHi.SetLineColor(kBlue);
     fHistSatHi.SetFillStyle(4000);
 
@@ -72 +73 @@
     fHistSatLo.SetXTitle("Pixels");
     fHistSatLo.SetYTitle("Counts");
-    fHistSatLo.SetLineColor(kBlue);
     fHistSatLo.SetDirectory(NULL);
     fHistSatLo.UseCurrentStyle();
@@ -83 +83 @@
     fHistIslands.SetDirectory(NULL);
     fHistIslands.UseCurrentStyle();
-    fHistIslands.SetLineColor(kBlue);
+    //fHistIslands.SetLineColor(kBlue);
     fHistIslands.SetFillStyle(4000);
 
+    fHistSizeSubIslands.SetName("SizeSub");
+    fHistSizeSubIslands.SetTitle("Size of Sub Islands");
+    fHistSizeSubIslands.SetXTitle("S [phe]");
+    fHistSizeSubIslands.SetYTitle("Counts");
+    fHistSizeSubIslands.SetDirectory(NULL);
+    fHistSizeSubIslands.UseCurrentStyle();
+    fHistSizeSubIslands.SetLineColor(kBlue);
+    fHistSizeSubIslands.SetFillStyle(4000);
+
+    fHistSizeMainIsland.SetName("SizeMain");
+    fHistSizeMainIsland.SetTitle("Size of Main Island");
+    fHistSizeMainIsland.SetXTitle("S [phe]");
+    fHistSizeMainIsland.SetYTitle("Counts");
+    fHistSizeMainIsland.SetDirectory(NULL);
+    fHistSizeMainIsland.UseCurrentStyle();
+    fHistSizeMainIsland.SetFillStyle(4000);
+
+    fHistNumSP.SetName("NumSP");
+    fHistNumSP.SetTitle("Number of single core pixels");
+    fHistNumSP.SetXTitle("N");
+    fHistNumSP.SetYTitle("Counts");
+    fHistNumSP.SetDirectory(NULL);
+    fHistNumSP.UseCurrentStyle();
+    //fHistNumSP.SetLineColor(kBlue);
+    fHistNumSP.SetFillStyle(4000);
+
+    fHistSizeSP.SetName("SizeSP");
+    fHistSizeSP.SetTitle("Size of single core pixels");
+    fHistSizeSP.SetXTitle("S [phe]");
+    fHistSizeSP.SetYTitle("Counts");
+    fHistSizeSP.SetDirectory(NULL);
+    fHistSizeSP.UseCurrentStyle();
+    //fHistSizeSP.SetLineColor(kBlue);
+    fHistSizeSP.SetFillStyle(4000);
+
+
     MBinning bins;
 
     bins.SetEdges(60, -0.5, 59.5);
+    bins.Apply(fHistSatLo);
     bins.Apply(fHistSatHi);
-    bins.Apply(fHistSatHi);
+    bins.Apply(fHistNumSP);
+
+    bins.SetEdgesLog(50, 1, 1e7);
+    bins.Apply(fHistSizeSubIslands);
+    bins.Apply(fHistSizeMainIsland);
+    bins.Apply(fHistSizeSP);
 
     bins.SetEdges(15, 0.5, 15.5);
@@ -103 +145 @@
 Bool_t MHImagePar::SetupFill(const MParList *plist)
 {
+    ApplyBinning(*plist, "Pixels",  &fHistSatLo);
     ApplyBinning(*plist, "Pixels",  &fHistSatHi);
-    ApplyBinning(*plist, "Pixels",  &fHistSatHi);
+    ApplyBinning(*plist, "Pixels",  &fHistNumSP);
 
     ApplyBinning(*plist, "Islands", &fHistIslands);
+
+    ApplyBinning(*plist, "Size",    &fHistSizeSubIslands);
+    ApplyBinning(*plist, "Size",    &fHistSizeMainIsland);
+    ApplyBinning(*plist, "Size",    &fHistSizeSP);
 
     return kTRUE;
@@ -126 +173 @@
     const MImagePar &h = *(MImagePar*)par;
 
-    fHistSatHi.Fill(h.GetNumSatPixelsHG(), w);
-    fHistSatLo.Fill(h.GetNumSatPixelsLG(), w);
-
-    fHistIslands.Fill(h.GetNumIslands(),   w);
+    fHistSatHi.Fill(h.GetNumSatPixelsHG(),    w);
+    fHistSatLo.Fill(h.GetNumSatPixelsLG(),    w);
+    fHistNumSP.Fill(h.GetNumSinglePixels(),   w);
+    fHistSizeSP.Fill(h.GetSizeSinglePixels(), w);
+    fHistSizeSubIslands.Fill(h.GetSizeSubIslands(), w);
+    fHistSizeMainIsland.Fill(h.GetSizeMainIsland(), w);
+    fHistIslands.Fill(h.GetNumIslands(),      w);
 
     return kTRUE;
@@ -136 +186 @@
 void MHImagePar::Paint(Option_t *o)
 {
-    /*
-     if (TString(o)==(TString)"log1" && fHistSatHi.GetMaximum()>0)
-     gPad->SetLogy();
-     if (TString(o)==(TString)"log2" && fHistIslands.GetMaximum()>0)
-     gPad->SetLogy();
-     */
      if (fHistSatHi.GetMaximum()>0 && gPad->GetPad(1))
          gPad->GetPad(1)->SetLogy();
-     if (fHistIslands.GetMaximum()>0 && gPad->GetPad(2))
-         gPad->GetPad(2)->SetLogy();
+     if (fHistIslands.GetMaximum()>0 && gPad->GetPad(3))
+         gPad->GetPad(3)->SetLogy();
+
+     TVirtualPad *pad = gPad->GetPad(2);
+     if (pad)
+     {
+         if (fHistNumSP.GetMaximum()>0 && pad->GetPad(1))
+             pad->GetPad(1)->SetLogy();
+         if (fHistSizeSP.GetMaximum()>0 && pad->GetPad(2))
+         {
+             pad->GetPad(2)->SetLogx();
+             pad->GetPad(2)->SetLogy();
+         }
+     }
+     if (fHistSizeMainIsland.GetMaximum()>0 && gPad->GetPad(4))
+     {
+         gPad->GetPad(4)->SetLogx();
+         gPad->GetPad(4)->SetLogy();
+     }
 }
 
@@ -168 +229 @@
 
     if (!same)
-        pad->Divide(1,2);
+        pad->Divide(2,2);
     else
         fHistIslands.SetLineColor(kGreen);
@@ -177 +238 @@
         gPad->SetBorderMode(0);
         MH::DrawSame(fHistSatHi, fHistSatLo, "Saturating Pixels");
-        fHistSatHi.SetMinimum(); // switch off to allow log-scale
-        fHistSatLo.SetMinimum(); // switch off to allow log-scale
-        fHistSatLo.SetMaximum(0.1);   // dummy value to allow log-scale
-        //AppendPad("log1");
+        fHistSatHi.SetMinimum();    // switch off to allow log-scale
+        fHistSatLo.SetMinimum();    // switch off to allow log-scale
+        fHistSatLo.SetMaximum(0.1); // dummy value to allow log-scale
+
+        pad->cd(4);
+        gPad->SetBorderMode(0);
+        MH::DrawSame(fHistSizeMainIsland, fHistSizeSubIslands, "Sizes...");
+        fHistSizeMainIsland.SetMinimum();    // switch off to allow log-scale
+        fHistSizeSubIslands.SetMinimum();    // switch off to allow log-scale
+        fHistSizeSubIslands.SetMaximum(0.1); // dummy value to allow log-scale
     }
 
     pad->cd(2);
     gPad->SetBorderMode(0);
+    pad->GetPad(2)->Divide(1,2,0,0);
+    pad->GetPad(2)->cd(1);
+    gPad->SetBorderMode(0);
+    fHistNumSP.Draw(same?"same":"");
+    pad->GetPad(2)->cd(2);
+    gPad->SetBorderMode(0);
+    fHistSizeSP.Draw(same?"same":"");
+
+    pad->cd(3);
+    gPad->SetBorderMode(0);
     fHistIslands.Draw(same?"same":"");
-    //AppendPad("log2");
 }