Changeset 1504 for trunk/MagicSoft

Timestamp:

08/16/02 11:47:30 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (2 diffs)
• mhist/MH.cc (modified) (19 diffs)
• mhist/MH.h (modified) (1 diff)
• mhist/MHHillas.cc (modified) (5 diffs)
• mhist/MHHillasExt.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -7 +7 @@
 
    * mhist/MHHillas.cc, mhist/MHHillasExt.cc:
-     - added legends. (still work to be done with the stat boxes)
+     - use the new MH::Draw[Copy] function
+     - added names to the histograms in MHHillasExt
 
    * manalysis/MCerPhotPix.h:
@@ -15 +16 @@
      - added comment
 
+   * mhist/MH.[h,cc]:
+     - added function to draw two histograms in one pad with
+       a layout of two stat boxes and a legend: Draw[Copy]
+     - added comments
+     - added another sanity check in SetBinning(TH1*,TH1*)
 
 

trunk/MagicSoft/Mars/mhist/MH.cc

@@ -54 +54 @@
 #include <TGaxis.h>
 #include <TCanvas.h>
+#include <TLegend.h>
+#include <TPaveStats.h>
 
 #include "MLog.h"
@@ -133 +135 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Applies a given binning to a 1D-histogram
+//
 void MH::SetBinning(TH1 *h, const MBinning *binsx)
 {
@@ -161 +167 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Applies given binnings to the two axis of a 2D-histogram
+//
 void MH::SetBinning(TH2 *h, const MBinning *binsx, const MBinning *binsy)
 {
@@ -194 +204 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Applies given binnings to the three axis of a 3D-histogram
+//
 void MH::SetBinning(TH3 *h, const MBinning *binsx, const MBinning *binsy, const MBinning *binsz)
 {
@@ -232 +246 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Applies given binning (the n+1 edges)  to the axis of a 1D-histogram
+//
 void MH::SetBinning(TH1 *h, const TArrayD &binsx)
 {
@@ -239 +257 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Applies given binning (the n+1 edges) to the two axis of a
+// 2D-histogram
+//
 void MH::SetBinning(TH2 *h, const TArrayD &binsx, const TArrayD &binsy)
 {
@@ -248 +271 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Applies given binning (the n+1 edges) to the three axis of a
+// 3D-histogram
+//
 void MH::SetBinning(TH3 *h, const TArrayD &binsx, const TArrayD &binsy, const TArrayD &binsz)
 {
@@ -259 +287 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Applies the binning of a TAxis (eg from a root histogram) to the axis
+// of a 1D-histogram
+//
 void MH::SetBinning(TH1 *h, const TAxis *binsx)
 {
@@ -270 +303 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Applies the binnings of the TAxis' (eg from a root histogram) to the
+// two axis' of a 2D-histogram
+//
 void MH::SetBinning(TH2 *h, const TAxis *binsx, const TAxis *binsy)
 {
@@ -285 +323 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Applies the binnings of the TAxis' (eg from a root histogram) to the
+// three axis' of a 3D-histogram
+//
 void MH::SetBinning(TH3 *h, const TAxis *binsx, const TAxis *binsy, const TAxis *binsz)
 {
@@ -304 +347 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Applies the binnings of one root-histogram x to another one h
+// Both histograms must be of the same type: TH1, TH2 or TH3
+//
 void MH::SetBinning(TH1 *h, TH1 *x)
 {
@@ -311 +359 @@
         return;
     }
+    if (h->InheritsFrom(TH3::Class()) || x->InheritsFrom(TH3::Class()))
+        return;
     if (h->InheritsFrom(TH2::Class()) && x->InheritsFrom(TH2::Class()))
     {
@@ -316 +366 @@
         return;
     }
+    if (h->InheritsFrom(TH2::Class()) || x->InheritsFrom(TH2::Class()))
+        return;
     if (h->InheritsFrom(TH1::Class()) && x->InheritsFrom(TH1::Class()))
     {
@@ -323 +375 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Multiplies all entries in a TArrayD by a float f
+//
 void MH::ScaleArray(TArrayD &bins, Double_t f)
 {
@@ -329 +385 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Scales the binning of a TAxis by a float f
+//
 TArrayD MH::ScaleAxis(TAxis &axe, Double_t f)
 {
@@ -341 +401 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Scales the binning of one, two or three axis of a histogram by a float f
+//
 void MH::ScaleAxis(TH1 *h, Double_t fx, Double_t fy, Double_t fz)
 {
@@ -364 +428 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Tries to find a MBinning container with the name "Binning"+name
+// in the given parameter list. If it was found it is applied to the
+// given histogram. This is only valid for 1D-histograms
+//
 Bool_t MH::ApplyBinning(const MParList &plist, TString name, TH1 *h)
 {
@@ -442 +512 @@
 }
 
+// --------------------------------------------------------------------------
+//
+//  Returns the lowest entry in a histogram which is greater than gt (eg >0)
+//
 Double_t MH::GetMinimumGT(const TH1 &h, Double_t gt)
 {
@@ -456 +530 @@
     return min;
 }
+
+// --------------------------------------------------------------------------
+//
+// Draws a copy of the two given histograms. Uses title as the pad title.
+// Also layout the two statistic boxes and a legend.
+//
+void MH::DrawCopy(const TH1 &hist1, const TH1 &hist2, const TString title)
+{
+    //
+    // Draw first histogram
+    //
+    TH1 *h1 = (TH1*)((TH1&)hist1).DrawCopy();
+    gPad->Update();
+
+    //
+    // Rename first statistics box
+    //
+    TPaveStats &s1 = *(TPaveStats*)gPad->FindObject("stats");
+    s1.SetX1NDC(s1.GetX1NDC()-0.01);
+    s1.SetName((TString)"Stat"+hist1.GetTitle());
+
+    //
+    // Draw second histogram
+    //
+    TH1 *h2 = (TH1*)((TH1&)hist2).DrawCopy("sames");
+
+    gPad->Update();
+
+    //
+    // Set new position of second statistics box
+    //
+    TPaveStats &s2 = *(TPaveStats*)gPad->FindObject("stats");
+    s2.SetX1NDC(s1.GetX1NDC()-(s2.GetX2NDC()-s2.GetX1NDC())-0.01);
+    s2.SetX2NDC(s1.GetX1NDC()-0.01);
+
+    //
+    // Draw Legend
+    //
+    TLegend &l = *new TLegend(s1.GetX1NDC(),
+                              s1.GetY1NDC()-0.015-(s1.GetY2NDC()-s1.GetY1NDC())/2,
+                              s1.GetX2NDC(),
+                              s1.GetY1NDC()-0.01
+                            );
+    l.AddEntry(h1, h1->GetTitle());
+    l.AddEntry(h2, h2->GetTitle());
+    l.SetTextSize(s2.GetTextSize());
+    l.SetTextFont(s2.GetTextFont());
+    l.SetBorderSize(s2.GetBorderSize());
+    l.SetBit(kCanDelete);
+
+    // FIXME: Don't change the hist title. Change TPaveText::title
+    h1->SetTitle(" "+title+" ");
+
+    l.Draw();
+}
+
+// --------------------------------------------------------------------------
+//
+// Draws the two given histograms. Uses title as the pad title.
+// Also layout the two statistic boxes and a legend.
+//
+void MH::Draw(TH1 &hist1, TH1 &hist2, const TString title)
+{
+    //
+    // Draw first histogram
+    //
+    hist1.Draw();
+    // FIXME: Don't change the hist title. Change TPaveText::title
+    //h1->SetTitle(" "+title+" ");
+
+    gPad->Update();
+
+    //
+    // Rename first statistics box
+    //
+    TPaveStats &s1 = *(TPaveStats*)gPad->FindObject("stats");
+    s1.SetX1NDC(s1.GetX1NDC()-0.01);
+    s1.SetName((TString)"Stat"+hist1.GetTitle());
+
+    //
+    // Draw second histogram
+    //
+    hist2.Draw("sames");
+
+    gPad->Update();
+
+    //
+    // Set new position of second statistics box
+    //
+    TPaveStats &s2 = *(TPaveStats*)gPad->FindObject("stats");
+    s2.SetX1NDC(s1.GetX1NDC()-(s2.GetX2NDC()-s2.GetX1NDC())-0.01);
+    s2.SetX2NDC(s1.GetX1NDC()-0.01);
+
+    //
+    // Draw Legend
+    //
+    TLegend &l = *new TLegend(s1.GetX1NDC(),
+                              s1.GetY1NDC()-0.015-(s1.GetY2NDC()-s1.GetY1NDC())/2,
+                              s1.GetX2NDC(),
+                              s1.GetY1NDC()-0.01
+                             );
+    l.AddEntry(&hist1, hist1.GetTitle());
+    l.AddEntry(&hist2, hist2.GetTitle());
+    l.SetTextSize(s2.GetTextSize());
+    l.SetTextFont(s2.GetTextFont());
+    l.SetBorderSize(s2.GetBorderSize());
+    l.SetBit(kCanDelete);
+    l.Draw();
+}

trunk/MagicSoft/Mars/mhist/MH.h

@@ -50 +50 @@
     static void    ScaleAxis(TH1 *bins, Double_t fx=1, Double_t fy=1, Double_t fz=1);
 
+    static void DrawCopy(const TH1 &hist1, const TH1 &hist2, const TString title);
+    static void Draw(TH1 &hist1, TH1 &hist2, const TString title);
+
     static void FindGoodLimits(Int_t nbins, Int_t &newbins, Double_t &xmin, Double_t &xmax, Bool_t isInteger);
     static Double_t GetMinimumGT(const TH1 &h, Double_t gt=0);

trunk/MagicSoft/Mars/mhist/MHHillas.cc

@@ -41 +41 @@
 #include <TStyle.h>
 #include <TCanvas.h>
-#include <TLegend.h>
 
 #include "MLog.h"
@@ -295 +294 @@
 
     c->cd(1);
-
-    TObject *obj;
-    TLegend *l = new TLegend(0.75, 0.8, 0.9, 0.9);
-    obj = fWidth->DrawCopy();
-    ((TH1*)obj)->SetTitle(" Width / Length ");
-    l->AddEntry(obj, "Width");
-    obj = fLength->DrawCopy("same");
-    l->AddEntry(obj,  "Length");
-    l->SetBit(kCanDelete);
-    l->Draw();
+    DrawCopy(*fWidth, *fLength, "Width / Length");
 
     c->cd(2);
@@ -311 +301 @@
 
     c->cd(3);
-    l = new TLegend(0.70, 0.8, 0.9, 0.9);
-    obj = fCorePix->DrawCopy();
-    ((TH1*)obj)->SetTitle(" Number of core/used Pixels ");
-    l->AddEntry(obj, "# Core Pixels");
-    obj = fUsedPix->DrawCopy("same");
-    l->AddEntry(obj,  "# Used Pixels");
-    l->SetBit(kCanDelete);
-    l->Draw();
+    DrawCopy(*fCorePix, *fUsedPix, "Number of core/used Pixels");
 
     c->cd(4);
@@ -350 +333 @@
 
     gPad->cd(1);
-    fWidth->Draw();
-    fLength->Draw("same");
+    MH::Draw(*fWidth, *fLength, "Width / Length");
 
     gPad->cd(2);
@@ -358 +340 @@
 
     gPad->cd(3);
-    fCorePix->Draw();
-    fUsedPix->Draw("same");
+    MH::Draw(*fCorePix, *fUsedPix, "Number of core/used Pixels");
 
     gPad->cd(4);

trunk/MagicSoft/Mars/mhist/MHHillasExt.cc

@@ -80 +80 @@
     fHM3Trans.SetDirectory(NULL);
 
+    fHConc.SetName("Conc2");
+    fHConc1.SetName("Conc1");
+    fHAsym.SetName("Asymmetry");
+    fHM3Long.SetName("3rd Mom Long");
+    fHM3Trans.SetName("3rd Mom Trans");
+
     fHConc.SetTitle("Ratio: Conc");
     fHConc1.SetTitle("Ratio: Conc1");
@@ -261 +267 @@
 
     c.cd(1);
-    TObject *obj;
-    TLegend *l = new TLegend(0.75, 0.8, 0.9, 0.9);
-    obj = ((TH1F)fHConc1).DrawCopy();
-    ((TH1*)obj)->SetTitle(" Concentrations ");
-    l->AddEntry(obj, "Conc1");
-    obj = ((TH1F)fHConc).DrawCopy("same");
-    l->AddEntry(obj,  "Conc2");
-    l->SetBit(kCanDelete);
-    l->Draw();
+    DrawCopy(fHConc1, fHConc, "Concentrations");
 
     c.cd(2);
-    ((TH1F)fHAsym).DrawCopy();
+    ((TH1&)fHAsym).DrawCopy();
 
     c.cd(3);
-    ((TH1F)fHM3Long).DrawCopy();
+    ((TH1&)fHM3Long).DrawCopy();
 
     c.cd(4);
-    ((TH1F)fHM3Trans).DrawCopy();
+    ((TH1&)fHM3Trans).DrawCopy();
 
     c.Modified();
@@ -300 +298 @@
 
     gPad->cd(1);
-    fHConc1.Draw();
-    fHConc.Draw("same");
-
-    TLegend *l = new TLegend(0.75, 0.8, 0.9, 0.9);
-    l->AddEntry((TH1*)&fHConc1, "Conc 1");
-    l->AddEntry((TH1*)&fHConc,  "Conc");
-    l->SetBit(kCanDelete);
-    l->Draw();
+    MH::Draw(fHConc1, fHConc, "Concentrations");
 
     gPad->cd(2);