Changeset 1524 for trunk/MagicSoft/Mars/mhist/MH3.cc

Timestamp:

09/16/02 10:10:17 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MH3.cc (modified) (11 diffs)

trunk/MagicSoft/Mars/mhist/MH3.cc

@@ -62 +62 @@
 #include <fstream.h>
 
+#include <TPad.h>
+#include <TStyle.h>
+#include <TCanvas.h>
+
 #include <TH2.h>
 #include <TH3.h>
 #include <TProfile.h>
 #include <TProfile2D.h>
-#include <TPad.h>
-#include <TCanvas.h>
 
 #include "MLog.h"
@@ -73 +75 @@
 
 #include "MParList.h"
-
+#include "MBinning.h"
 #include "MDataChain.h"
 
@@ -202 +204 @@
             return kFALSE;
         }
+        if (binsz->IsLogarithmic())
+            fHist->SetBit(kIsLogz);
         fHist->SetZTitle(fData[2]->GetTitle());
         if (!fData[2]->PreProcess(plist))
@@ -212 +216 @@
             return kFALSE;
         }
+        if (binsy->IsLogarithmic())
+            fHist->SetBit(kIsLogy);
         fHist->SetYTitle(fData[1]->GetTitle());
         if (!fData[1]->PreProcess(plist))
@@ -222 +228 @@
             return kFALSE;
         }
+        if (binsx->IsLogarithmic())
+            fHist->SetBit(kIsLogx);
         fHist->SetXTitle(fData[0]->GetTitle());
         if (!fData[0]->PreProcess(plist))
@@ -306 +314 @@
     return kFALSE;
 }
+/*
+// --------------------------------------------------------------------------
+//
+// Set the palette you wanna use:
+//  - you could set the root "Pretty Palette Violet->Red" by
+//    gStyle->SetPalette(1, 0), but in some cases this may look
+//    confusing
+//  - The maximum colors root allowes us to set by ourself
+//    is 50 (idx: 51-100). This colors are set to a grayscaled
+//    palette
+//  - the number of contours must be two less than the number
+//    of palette entries
+//
+void MHStarMap::PrepareDrawing() const
+{
+    const Int_t numg = 32; // number of gray scaled colors
+    const Int_t numw = 32; // number of white
+
+    Int_t palette[numg+numw];
+
+    //
+    // The first half of the colors are white.
+    // This is some kind of optical background supression
+    //
+    gROOT->GetColor(51)->SetRGB(1, 1, 1);
+
+    Int_t i;
+    for (i=0; i<numw; i++)
+        palette[i] = 51;
+
+    //
+    // now the (gray) scaled part is coming
+    //
+    for (;i<numw+numg; i++)
+    {
+        const Float_t gray = 1.0-(float)(i-numw)/(numg-1.0);
+
+        gROOT->GetColor(52+i)->SetRGB(gray, gray, gray);
+        palette[i] = 52+i;
+    }
+
+    //
+    // Set the palette and the number of contour levels
+    //
+    gStyle->SetPalette(numg+numw, palette);
+    fStarMap->SetContour(numg+numw-2);
+}
+*/
+// --------------------------------------------------------------------------
+//
+// Setup a inversed deep blue sea palette for the fCenter histogram.
+//
+void MH3::SetColors() const
+{
+    // FIXME: This must be redone each time the canvas is repainted....
+    gStyle->SetPalette(51, NULL);
+    Int_t c[50];
+    for (int i=0; i<50; i++)
+        c[49-i] = gStyle->GetColorPalette(i);
+    gStyle->SetPalette(50, c);
+}
 
 // --------------------------------------------------------------------------
 //
 // Draw clone of histogram. So that the object can be deleted
+//
+// Possible options are:
+//   PROFX: Draw a x-profile into the histogram (for 2D histograms only)
+//   PROFY: Draw a y-profile into the histogram (for 2D histograms only)
+//   ONLY:  Draw the profile histogram only (for 2D histograms only)
+//
+// If the kIsLog?-Bit is set the axis is displayed lkogarithmically.
+// eg this is set when applying a logarithmic MBinning
+//
 // and the histogram is still visible in the canvas.
 // The cloned object are deleted together with the canvas if the canvas is
@@ -324 +402 @@
     gROOT->SetSelectedPad(NULL);
 
-    fHist->DrawCopy(opt);
-
     TString str(opt);
+
+    if (str.Contains("COL", TString::kIgnoreCase))
+        SetColors();
+
+    Bool_t only = str.Contains("ONLY", TString::kIgnoreCase) && fDimension==2;
+    if (!only)
+        fHist->DrawCopy(opt);
+
     if (str.Contains("PROFX", TString::kIgnoreCase) && fDimension==2)
     {
         TProfile *p = ((TH2*)fHist)->ProfileX();
-        p->Draw("same");
+        p->Draw(only?"":"same");
         p->SetBit(kCanDelete);
         p->SetDirectory(NULL);
@@ -337 +421 @@
     {
         TProfile *p = ((TH2*)fHist)->ProfileY();
-        p->Draw("same");
+        p->Draw(only?"":"same");
         p->SetBit(kCanDelete);
         p->SetDirectory(NULL);
     }
 
+    if (fHist->TestBit(kIsLogx)) c.SetLogx();
+    if (fHist->TestBit(kIsLogy)) c.SetLogy();
+    if (fHist->TestBit(kIsLogz)) c.SetLogz();
+
     c.Modified();
     c.Update();
@@ -351 +439 @@
 //
 // Creates a new canvas and draws the histogram into it.
+//
+// Possible options are:
+//   PROFX: Draw a x-profile into the histogram (for 2D histograms only)
+//   PROFY: Draw a y-profile into the histogram (for 2D histograms only)
+//   ONLY:  Draw the profile histogram only (for 2D histograms only)
+//
+// If the kIsLog?-Bit is set the axis is displayed lkogarithmically.
+// eg this is set when applying a logarithmic MBinning
+//
 // Be careful: The histogram belongs to this object and won't get deleted
 // together with the canvas.
@@ -359 +456 @@
         MH::MakeDefCanvas(fHist);
 
-    fHist->Draw(opt);
-
     TString str(opt);
+
+    if (str.Contains("COL", TString::kIgnoreCase))
+        SetColors();
+
+    Bool_t only = str.Contains("ONLY", TString::kIgnoreCase) && fDimension==2;
+    if (!only)
+        fHist->Draw(opt);
+
     if (str.Contains("PROFX", TString::kIgnoreCase) && fDimension==2)
     {
         TProfile *p = ((TH2*)fHist)->ProfileX();
-        p->Draw("same");
+        p->Draw(only?"":"same");
         p->SetBit(kCanDelete);
         p->SetDirectory(NULL);
@@ -372 +475 @@
     {
         TProfile *p = ((TH2*)fHist)->ProfileY();
-        p->Draw("same");
+        p->Draw(only?"":"same");
         p->SetBit(kCanDelete);
         p->SetDirectory(NULL);
     }
+
+    if (fHist->TestBit(kIsLogx)) gPad->SetLogx();
+    if (fHist->TestBit(kIsLogy)) gPad->SetLogy();
+    if (fHist->TestBit(kIsLogz)) gPad->SetLogz();
 
     gPad->Modified();