Changeset 1216 for trunk/MagicSoft/Mars/mhist/MHHillas.cc

Timestamp:

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

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHHillas.cc (modified) (6 diffs)

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;