Changeset 2265 for trunk/MagicSoft/Mars/mhist/MHCamera.cc

Timestamp:

07/06/03 21:02:56 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHCamera.cc (modified) (15 diffs)

trunk/MagicSoft/Mars/mhist/MHCamera.cc

@@ -47 +47 @@
 #include <TRandom.h>
 #include <TPaveText.h>
+#include <TPaveStats.h>
 #include <TClonesArray.h>
+#include <THistPainter.h>
+
+#include "MLog.h"
+#include "MLogManip.h"
 
 #include "MH.h"
@@ -76 +81 @@
 //  Default Constructor. To be used by the root system ONLY.
 //
-MHCamera::MHCamera() : TH1D(), fGeomCam(NULL), fColors(kItemsLegend)
+MHCamera::MHCamera() : TH1D(), fGeomCam(NULL), fColors(kItemsLegend), fOptStat(-1)
 {
     SetDirectory(NULL);
@@ -215 +220 @@
    fTsumwx2 += z*x*x;
    return bin;
+}
+
+Stat_t MHCamera::GetMean(Int_t axis) const
+{
+    Stat_t mean = 0;
+    for (int i=1; i<fNcells-1; i++)
+        mean += fArray[i];
+
+    return mean/(fNcells-2);
+}
+
+Stat_t MHCamera::GetRMS(Int_t axis) const
+{
+    const Int_t n = fNcells-2;
+
+    Stat_t sum = 0;
+    Stat_t sq  = 0;
+    for (int i=1; i<n+1; i++)
+    {
+        sum += fArray[i];
+        sq  += fArray[i]*fArray[i];
+    }
+
+    sum /= n;
+    sq  /= n;
+
+    return sqrt(sq-sum*sum);
 }
 
@@ -314 +346 @@
 void MHCamera::SetRange()
 {
-    const Float_t range = fGeomCam->GetMaxRadius();
+    const Float_t range = fGeomCam->GetMaxRadius()*1.05;
 
     //
@@ -373 +405 @@
     {
         if (IsUsed(i) && iscol)
+        {
+            if (TMath::IsNaN(fArray[i+1]))
+                gLog << warn << "MHCamera::Update: " << GetName() << " <" << GetTitle() << "> - Pixel Index #" << i << " contents is NaN (Not a Number)..." << endl;
+
             hex.SetFillColor(GetColor(fArray[i+1], min, max, islog));
+        }
         else
             hex.SetFillColor(10);
@@ -381 +418 @@
             hex.PaintHexagon(pix.GetX(), pix.GetY(), pix.GetD());
         else
-            if (IsUsed(i))
+            if (IsUsed(i) && !TMath::IsNaN(fArray[i+1]))
             {
                 Float_t size = pix.GetD()*(fArray[i+1]-min)/(max-min);
@@ -413 +450 @@
 void MHCamera::PaintAxisTitle()
 {
-    Float_t fRange = fGeomCam->GetMaxRadius();
-
-    TLatex *ptitle = new TLatex(1.2*fRange, .97*fRange, GetYaxis()->GetTitle());
-
-    ptitle->SetTextSize(0.03);
-    ptitle->SetTextAlign(33);
+    const Float_t range = fGeomCam->GetMaxRadius()*1.05;
+    const Float_t w = (1 + 1.5/sqrt((float)(fNcells-2)))*range;
+
+    TLatex *ptitle = new TLatex(w, -.90*range, GetYaxis()->GetTitle());
+
+    ptitle->SetTextSize(0.05);
+    ptitle->SetTextAlign(21);
 
     // box with the histogram title
@@ -430 +468 @@
 // ------------------------------------------------------------------------
 //
-// Paint the histogram title
-//
-void MHCamera::PaintTitle()
-{
-//    *-*-*-*-*-*-*-*-*-*Draw the histogram title*-*-*-*-*-*-*-*-*-*-*-*-*
-//                       ========================
-   //if (Hoption.Same) return;
-#if ROOT_VERSION_CODE > ROOT_VERSION(3,05,00)
-    if (TestBit(kNoTitle))
-        return;
-#endif
-
-    const Int_t nt = strlen(GetTitle());
-
-    TPaveText *title = (TPaveText*)gPad->FindObject("title");
-    if (nt == 0 || gStyle->GetOptTitle() <= 0)
-    {
-        if (title)
-            delete title;
-        return;
-    }
-
-    Double_t ht = gStyle->GetTitleH();
-    Double_t wt = gStyle->GetTitleW();
-
-    if (ht <= 0)
-        ht = 0.05;
-    if (wt <= 0)
-    {
-        TLatex l;
-        l.SetTextSize(ht);
-        l.SetTitle(GetTitle());
-        Double_t wndc = l.GetXsize()/(gPad->GetX2() - gPad->GetX1());
-        wt = TMath::Min(0.7, 0.02+wndc);
-    }
-    if (title)
-    {
-        TText *t0 = (TText*)title->GetLine(0);
-        if (t0)
-        {
-            if (!strcmp(t0->GetTitle(), GetTitle()))
-                return;
-
-            t0->SetTitle(GetTitle());
-            if (wt > 0)
-                title->SetX2NDC(title->GetX1NDC()+wt);
-        }
-        return;
-    }
-
-    TPaveText *ptitle = new TPaveText(
-                                      gStyle->GetTitleX(),
-                                      gStyle->GetTitleY()-ht,
-                                      gStyle->GetTitleX()+wt,
-                                      gStyle->GetTitleY(),"blNDC");
-
-    // box with the histogram title
-#if ROOT_VERSION_CODE > ROOT_VERSION(3,05,01)
-    ptitle->SetFillColor(gStyle->GetTitleFillColor());
-    ptitle->SetTextFont(gStyle->GetTitleFont(""));
-    if (gStyle->GetTitleFont("")%10 > 2)
-        ptitle->SetTextSize(gStyle->GetTitleFontSize());
-#endif
-    ptitle->SetFillStyle(gStyle->GetTitleStyle());
-    ptitle->SetName("title");
-    ptitle->SetBorderSize(gStyle->GetTitleBorderSize());
-    ptitle->SetTextColor(gStyle->GetTitleTextColor());
-    ptitle->AddText(GetTitle());
-    ptitle->SetBit(kCanDelete);
-    ptitle->Draw();
-    ptitle->Paint();
-}
-
-// ------------------------------------------------------------------------
-//
 // Paints the camera.
 //
@@ -515 +478 @@
     {
         opt.ReplaceAll("hist", "");
-
         Int_t mode = gStyle->GetOptStat();
         TVirtualPad *save = gPad;
         gPad=NULL;
-        gStyle->SetOptStat(1000011);
+        gStyle->SetOptStat(fOptStat<0 ? mode : fOptStat);
         gPad=save;
         TH1D::Paint(o);
@@ -534 +496 @@
     Bool_t iscol = isbox ? !opt.Contains("nocol") : 1;
 
+    THistPainter paint;
+    if (!TestBit(TH1::kNoStats))
+    {
+        paint.SetHistogram(this);
+        paint.PaintStat(fOptStat<0?gStyle->GetOptStat():fOptStat, NULL);
+    }
+
     // Update Contents of the pixels and paint legend
     Update(gPad->GetLogy(), isbox, iscol);
 
     // Paint primitives (pixels, color legend, photons, ...)
-    PaintTitle();
+    paint.PaintTitle();
     PaintAxisTitle();
 }
@@ -619 +588 @@
         const MGeomPix &h = (*fGeomCam)[i];
         TText *nt = txt.DrawText(h.GetX(), h.GetY(), num);
-        nt->SetTextSize(0.3*h.GetD()/fGeomCam->GetMaxRadius());
+        nt->SetTextSize(0.3*h.GetD()/fGeomCam->GetMaxRadius()/1.05);
     }
 }
@@ -642 +611 @@
         const MGeomPix &h = (*fGeomCam)[i];
         TText *nt = txt.DrawText(h.GetX(), h.GetY(), num);
-        nt->SetTextSize(0.3*h.GetD()/fGeomCam->GetMaxRadius());
+        nt->SetTextSize(0.3*h.GetD()/fGeomCam->GetMaxRadius()/1.05);
     }
 }
@@ -848 +817 @@
 {
     if (TMath::IsNaN(val)) // FIXME: gLog!
-    {
-        cout << "MHCamera::GetColor: " << GetName() << " <" << GetTitle() << "> - Color for val=NaN (Not a Number) requested... set val=min." << endl;
-        val = min;
-    }
+        return 10;
 
     //
@@ -877 +843 @@
 }
 
+TPaveStats *MHCamera::GetStats()
+{
+    TObject *obj = 0;
+
+    TIter Next(fFunctions);
+    while ((obj = Next()))
+        if (obj->InheritsFrom(TPaveStats::Class()))
+            return static_cast<TPaveStats*>(obj);
+
+    return NULL;
+}
+
 // ------------------------------------------------------------------------
 //
@@ -883 +861 @@
 void MHCamera::UpdateLegend(Float_t min, Float_t max, Bool_t islog)
 {
-    const Float_t range = fGeomCam->GetMaxRadius();
-
-    const Float_t H = 0.9*range;
+    TPaveStats *stats = GetStats();
+
+    const Float_t hndc   = 0.92 - (stats ? stats->GetY1NDC() : 1);
+    const Float_t range  = fGeomCam->GetMaxRadius()*1.05;
+    const Float_t H      = (0.75-hndc)*range;
+    const Float_t offset = hndc*range;
+
     const Float_t h = 2./kItemsLegend;
-    const Float_t offset = 0.04*range;
-
     const Float_t w = range/sqrt((float)(fNcells-2));
 
     TBox newbox;
     TText newtxt;
-    newtxt.SetTextSize(0.025);
+    newtxt.SetTextSize(0.03);
     newtxt.SetTextAlign(12);
 #if ROOT_VERSION_CODE > ROOT_VERSION(3,01,06)