Changeset 6869 for trunk/MagicSoft/Mars/mimage/MHNewImagePar.cc

Timestamp:

03/21/05 10:39:58 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mimage/MHNewImagePar.cc (modified) (7 diffs)

trunk/MagicSoft/Mars/mimage/MHNewImagePar.cc

@@ -56 +56 @@
 //
 MHNewImagePar::MHNewImagePar(const char *name, const char *title)
+    : fMm2Deg(1), fUseMmScale(kTRUE)
 {
     fName  = name  ? name  : "MHNewImagePar";
@@ -142 +143 @@
     bins.Apply(fHistCorePix);
 
-    bins.SetEdges(75, 0, 0.249);
-    bins.Apply(fHistUsedArea);
-    bins.Apply(fHistCoreArea);
+    //bins.SetEdges(75, 0, 0.249);
+    //bins.Apply(fHistUsedArea);
+    //bins.Apply(fHistCoreArea);
 }
 
@@ -154 +155 @@
 Bool_t MHNewImagePar::SetupFill(const MParList *plist)
 {
+    MGeomCam *geom = (MGeomCam*)plist->FindObject("MGeomCam");
+    if (!geom)
+        *fLog << warn << GetDescriptor() << ": No Camera Geometry available. Using mm-scale for histograms." << endl;
+    else
+    {
+        fMm2Deg = geom->GetConvMm2Deg();
+        SetMmScale(kFALSE);
+    }
+
+    const MBinning *bins = (MBinning*)plist->FindObject("BinningArea");
+    if (!bins)
+    {
+        float r = geom ? 1.5 : 0.133;
+
+        MBinning b;
+        b.SetEdges(50, 0, r);
+        b.Apply(fHistUsedArea);
+        b.Apply(fHistCoreArea);
+    }
+    else
+    {
+        bins->Apply(fHistUsedArea);
+        bins->Apply(fHistCoreArea);
+    }
+
     ApplyBinning(*plist, "Leakage", &fHistLeakage1);
     ApplyBinning(*plist, "Leakage", &fHistLeakage2);
@@ -160 +186 @@
     ApplyBinning(*plist, "Pixels",  &fHistCorePix);
 
-    ApplyBinning(*plist, "Area",    &fHistUsedArea);
-    ApplyBinning(*plist, "Area",    &fHistCoreArea);
+    //ApplyBinning(*plist, "Area",    &fHistUsedArea);
+    //ApplyBinning(*plist, "Area",    &fHistCoreArea);
 
     ApplyBinning(*plist, "Conc",    &fHistConc);
@@ -182 +208 @@
     }
 
+    const Double_t scale = fUseMmScale ? 1e-6 : fMm2Deg*fMm2Deg;
+
     const MNewImagePar &h = *(MNewImagePar*)par;
 
@@ -190 +218 @@
     fHistCorePix.Fill(h.GetNumCorePixels(), w);
 
-    fHistUsedArea.Fill(h.GetUsedArea()/1000000, w);
-    fHistCoreArea.Fill(h.GetCoreArea()/1000000, w);
+    fHistUsedArea.Fill(h.GetUsedArea()*scale, w);
+    fHistCoreArea.Fill(h.GetCoreArea()*scale, w);
 
     fHistConc.Fill(h.GetConc(), w);
@@ -197 +225 @@
 
     return kTRUE;
+}
+
+// --------------------------------------------------------------------------
+//
+// With this function you can convert the histogram ('on the fly') between
+// degrees and millimeters.
+//
+void MHNewImagePar::SetMmScale(Bool_t mmscale)
+{
+    if (fUseMmScale == mmscale)
+        return;
+
+    if (fMm2Deg<0)
+    {
+        *fLog << warn << dbginf << "Warning - Sorry, no conversion factor for conversion available." << endl;
+        return;
+    }
+
+    const Double_t scale = mmscale ? 1./(fMm2Deg*fMm2Deg) : (fMm2Deg*fMm2Deg);
+    MH::ScaleAxis(&fHistUsedArea,  scale);
+    MH::ScaleAxis(&fHistCoreArea,  scale);
+
+    if (mmscale)
+    {
+        fHistUsedArea.SetXTitle("A [m^{2}]");
+        fHistCoreArea.SetXTitle("A [m^{2}]");
+    }
+    else
+    {
+        fHistUsedArea.SetXTitle("A [deg^{2}]");
+        fHistCoreArea.SetXTitle("A [deg^{2}]");
+    }
+
+    fUseMmScale = mmscale;
+}
+
+// --------------------------------------------------------------------------
+//
+// Use this function to setup your own conversion factor between degrees
+// and millimeters. The conversion factor should be the one calculated in
+// MGeomCam. Use this function with Caution: You could create wrong values
+// by setting up your own scale factor.
+//
+void MHNewImagePar::SetMm2Deg(Float_t mmdeg)
+{
+    if (mmdeg<0)
+    {
+        *fLog << warn << dbginf << "Warning - Conversion factor < 0 - nonsense. Ignored." << endl;
+        return;
+    }
+
+    if (fMm2Deg>=0)
+        *fLog << warn << dbginf << "Warning - Conversion factor already set. Overwriting" << endl;
+
+    fMm2Deg = mmdeg;
 }