Changeset 1216 for trunk/MagicSoft/Mars/mhist

Timestamp:

02/13/02 16:40:37 (23 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhist

Files:

• MHHillas.cc (modified) (6 diffs)
• MHHillas.h (modified) (1 diff)
• MHHillasSrc.cc (modified) (5 diffs)
• MHHillasSrc.h (modified) (1 diff)
• Makefile (modified) (1 diff)

trunk/MagicSoft/Mars/mhist/MHHillas.cc

@@ -39 +39 @@
 #include <TCanvas.h>
 
+#include "MLog.h"
+#include "MLogManip.h"
+
+#include "MGeomCam.h"
 #include "MHillas.h"
 #include "MParList.h"
@@ -49 +53 @@
 //
 MHHillas::MHHillas(const char *name, const char *title)
+    : fMm2Deg(1), fUseMmScale(kFALSE)
 {
     //
@@ -61 +66 @@
     // connect all the histogram with the container fHist
     //
-    fWidth  = new TH1F("Width [mm]",  "Width of Hillas",  100, 0, 300);
-    fLength = new TH1F("Length [mm]", "Length of Hillas", 100, 0, 300);
+    fWidth  = new TH1F("Width",  "Width of Ellipse",  100, 0, 300);
+    fLength = new TH1F("Length", "Length of Ellipse", 100, 0, 300);
 
     fLength->SetDirectory(NULL);
@@ -89 +94 @@
 // instances of MBinning (with the names 'BinningWidth' and 'BinningLength')
 // are found in the parameter list
+// Use this function if you want to set the conversion factor which
+// is used to convert the mm-scale in the camera plain into the deg-scale
+// used for histogram presentations. The conversion factor is part of
+// the camera geometry. Please create a corresponding MGeomCam container.
 //
 Bool_t MHHillas::SetupFill(const MParList *plist)
@@ -94 +103 @@
     const MBinning* binsw = (MBinning*)plist->FindObject("BinningWidth");
     const MBinning* binsl = (MBinning*)plist->FindObject("BinningLength");
-
-    if (binsw)
-        SetBinning(fWidth, binsw);
-
-    if (binsl)
-        SetBinning(fLength, binsl);
+    if (!binsw || !binsl)
+    {
+        *fLog << err << dbginf << "At least one MBinning not found... aborting." << endl;
+        return kFALSE;
+    }
+
+    SetBinning(fWidth,  binsw);
+    SetBinning(fLength, binsl);
+
+    const MGeomCam *geom = (MGeomCam*)plist->FindObject("MGeomCam");
+    if (!geom)
+    {
+        *fLog << warn << dbginf << "No Camera Geometry available. Using mm-scale for histograms." << endl;
+        return kTRUE;
+    }
+
+    fLength->GetXaxis()->SetTitle("Length [\\circ]");
+    fWidth->GetXaxis()->SetTitle("Width [\\circ]");
+
+    fMm2Deg = geom->GetConvMm2Deg();
 
     return kTRUE;
@@ -113 +136 @@
     const MHillas &h = *(MHillas*)par;
 
-    fWidth ->Fill(h.GetWidth());
-    fLength->Fill(h.GetLength());
+    if (fUseMmScale)
+    {
+        fWidth ->Fill(h.GetWidth());
+        fLength->Fill(h.GetLength());
+    }
+    else
+    {
+        fWidth ->Fill(fMm2Deg*h.GetWidth());
+        fLength->Fill(fMm2Deg*h.GetLength());
+    }
 
     return kTRUE;

trunk/MagicSoft/Mars/mhist/MHHillas.h

@@ -15 +15 @@
     TH1F *fLength;
 
+    Float_t fMm2Deg;
+
+    Bool_t fUseMmScale;
+
 public:
     MHHillas(const char *name=NULL, const char *title=NULL);
     ~MHHillas();
+
+    void SetMmScale(Bool_t mmscale=kTRUE) { fUseMmScale = mmscale; }
 
     Bool_t SetupFill(const MParList *pList);

trunk/MagicSoft/Mars/mhist/MHHillasSrc.cc

@@ -38 +38 @@
 #include <TCanvas.h>
 
+#include "MLog.h"
+#include "MLogManip.h"
+
+#include "MGeomCam.h"
+
 #include "MParList.h"
 
@@ -62 +67 @@
     // connect all the histogram with the container fHist
     //
-    fAlpha = new TH1F("Alpha [deg]", "Alpha of Hillas",   90, 0,  90);
-    fDist  = new TH1F("Dist [mm]",   "Dist of Hillas",   100, 0, 600);
+    fAlpha = new TH1F("Alpha", "Alpha of Ellipse",   90, 0,  90);
+    fDist  = new TH1F("Dist",  "Dist of Ellipse",   100, 0, 600);
 
     fAlpha->SetDirectory(NULL);
@@ -90 +95 @@
 // instances of MBinning (with the names 'BinningAlpha' and 'BinningDist')
 // are found in the parameter list
+// Use this function if you want to set the conversion factor which
+// is used to convert the mm-scale in the camera plain into the deg-scale
+// used for histogram presentations. The conversion factor is part of
+// the camera geometry. Please create a corresponding MGeomCam container.
 //
 Bool_t MHHillasSrc::SetupFill(const MParList *plist)
@@ -95 +104 @@
     const MBinning* binsa = (MBinning*)plist->FindObject("BinningAlpha");
     const MBinning* binsd = (MBinning*)plist->FindObject("BinningDist");
-
-    if (binsa)
-        SetBinning(fAlpha, binsa);
-
-    if (binsd)
-        SetBinning(fDist, binsd);
+    if (!binsa || !binsd)
+    {
+        *fLog << err << dbginf << "At least one MBinning not found... aborting." << endl;
+        return kFALSE;
+    }
+
+    SetBinning(fAlpha, binsa);
+    SetBinning(fDist,  binsd);
+
+    const MGeomCam *geom = (MGeomCam*)plist->FindObject("MGeomCam");
+    if (!geom)
+    {
+        *fLog << warn << dbginf << "No Camera Geometry available. Using mm-scale for histograms." << endl;
+        return kTRUE;
+    }
+
+    fDist->GetXaxis()->SetTitle("Dist [\\circ]");
+
+    fMm2Deg = geom->GetConvMm2Deg();
 
     return kTRUE;
@@ -115 +137 @@
 
     fAlpha->Fill(fabs(h.GetAlpha()));
-    fDist ->Fill(h.GetDist());
+    fDist ->Fill(fUseMmScale ? h.GetDist() : fMm2Deg*h.GetDist());
 
     return kTRUE;

trunk/MagicSoft/Mars/mhist/MHHillasSrc.h

@@ -15 +15 @@
     TH1F *fDist;
 
+    Float_t fMm2Deg;
+    Bool_t  fUseMmScale;
+
 public:
     MHHillasSrc(const char *name=NULL, const char *title=NULL);
     ~MHHillasSrc();
+
+    void SetUseMmScale(Bool_t mmscale=kTRUE) { fUseMmScale = mmscale; }
 
     Bool_t SetupFill(const MParList *pList);

trunk/MagicSoft/Mars/mhist/Makefile

@@ -22 +22 @@
 #  connect the include files defined in the config.mk file
 #
-INCLUDES = -I. -I../mbase -I../mraw -I../manalysis -I../mmc
+INCLUDES = -I. -I../mbase -I../mraw -I../manalysis -I../mmc -I../mgui
 
 #------------------------------------------------------------------------------