Changeset 4189 for trunk/MagicSoft/Mars/mhist

Timestamp:

05/26/04 11:18:32 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhist

Files:

• MHCamera.cc (modified) (2 diffs)
• MHCamera.h (modified) (2 diffs)

trunk/MagicSoft/Mars/mhist/MHCamera.cc

@@ -68 +68 @@
 #include <TH1.h>
 #include <TF1.h>
+#include <TCanvas.h>
 
 #include "MLog.h"
@@ -1850 +1851 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Draw the MHCamera into the MStatusDisplay: 
+// 
+// 1) Draw it as histogram (MHCamera::DrawCopy("hist")
+// 2) Draw it as a camera, with MHCamera::SetPrettyPalette() set. 
+// 3) If "rad" is not zero, draw its values vs. the radius from the camera center. 
+//    (DrawRadialProfile())
+// 4) Depending on the variable "fit", draw the values projection on the y-axis
+//    (DrawProjection()):
+//    0: don't draw
+//    1: Draw fit to Single Gauss (for distributions flat-fielded over the whole camera)
+//    2: Draw and fit to Double Gauss (for distributions different for inner and outer pixels)
+//    3: Draw and fit to Triple Gauss (for distributions with inner, outer pixels and outliers)
+//    4: Draw and fit to Polynomial grade 0: (for the probability distributions)
+//    >4: Draw and don;t fit.
+//
+void MHCamera::CamDraw(TCanvas &c, const Int_t x, const Int_t y, 
+                       const Int_t fit, const Int_t rad, TObject *notify)
+{
+
+  c.cd(x);
+  gPad->SetBorderMode(0);
+  gPad->SetTicks();
+  MHCamera *obj1=(MHCamera*)DrawCopy("hist");
+  obj1->SetDirectory(NULL);
+  
+  if (notify)
+    obj1->AddNotify(notify);
+
+  c.cd(x+y);
+  gPad->SetBorderMode(0);
+  obj1->SetPrettyPalette();
+  obj1->Draw();
+
+  if (rad)
+    {
+      c.cd(x+2*y);
+      gPad->SetBorderMode(0);
+      gPad->SetTicks();
+      DrawRadialProfile();
+    }
+  
+  if (!fit)
+    return;
+  
+  c.cd(rad ? x+3*y : x+2*y);
+  gPad->SetBorderMode(0);
+  gPad->SetTicks();
+  DrawProjection(fit);
+}
+
+

trunk/MagicSoft/Mars/mhist/MHCamera.h

@@ -16 +16 @@
 #ifndef ROOT_TH1
 #include <TH1.h>
+#endif
+#ifndef ROOT_TCanvas
+#include <TCanvas.h>
 #endif
 
@@ -249 +252 @@
     }
     TProfile *RadialProfileS(const TArrayI &sector, const TArrayI &aidx, const char *name="_rad", const Int_t nbins=25) const;
+
+    void CamDraw(TCanvas &c, const Int_t x, const Int_t y, 
+                 const Int_t fit, const Int_t rad=0,
+                 TObject *notify=NULL);             
     
     const MGeomCam &GetGeomCam() const { return *fGeomCam; }