Changeset 9341 for trunk/MagicSoft/Mars/mimage

Timestamp:

02/15/09 15:44:45 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mimage

Files:

• MHNewImagePar2.h (modified) (1 diff)
• MHVsSize.cc (modified) (5 diffs)
• MHVsSize.h (modified) (3 diffs)

trunk/MagicSoft/Mars/mimage/MHNewImagePar2.h

@@ -20 +20 @@
     TH1F fHistBorder2;  //
 
-    Float_t fMm2Deg;
-    Bool_t  fUseMmScale;
-
 public:
     MHNewImagePar2(const char *name=NULL, const char *title=NULL);

trunk/MagicSoft/Mars/mimage/MHVsSize.cc

@@ -59 +59 @@
 //
 MHVsSize::MHVsSize(const char *name, const char *title)
-    : fHillas(NULL), fHillasExt(NULL), fNewImagePar(NULL), fMm2Deg(1), fUseMmScale(kTRUE)
+    : fGeom(0), fHillas(0), fHillasExt(0), fNewImagePar(0)
 {
     //
@@ -88 +88 @@
     fM3Long.SetXTitle("Size [phe]");
 
-    fLength.SetYTitle("Length [mm]");
-    fWidth.SetYTitle("Width [mm]");
-    fDist.SetYTitle("Distance [mm]");
+    fLength.SetYTitle("Length [\\circ]");
+    fWidth.SetYTitle("Width [\\circ]");
+    fDist.SetYTitle("Distance [\\circ]");
     fConc1.SetYTitle("Conc1 [ratio]");
-    fArea.SetYTitle("Conc1 [mm^{2}]");
-    fM3Long.SetYTitle("M3Long [mm]");
+    fArea.SetYTitle("Area [\\circ^{2}]");
+    fM3Long.SetYTitle("M3Long [\\circ]");
 
     MBinning binse, binsl, binsd, binsc, binsa, binsm;
     binse.SetEdgesLog( 50,   10, 1e5);
-    binsl.SetEdges(   100,    0, 296.7/2);
-    binsd.SetEdges(   100,    0, 600);
+    binsl.SetEdges(   100,    0, 0.5);
+    binsd.SetEdges(   100,    0, 2.0);
     binsc.SetEdgesLog(100, 3e-3, 1);
-    binsa.SetEdges(   100,    0, 445*45);
-    binsm.SetEdges(   100, -445, 445);
+    binsa.SetEdges(   100,    0, 0.25);
+    binsm.SetEdges(   100, -1.5, 1.5);
 
     MH::SetBinning(&fLength, &binse, &binsl);
@@ -130 +130 @@
 Bool_t MHVsSize::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);
+    fGeom = (MGeomCam*)plist->FindObject("MGeomCam");
+    if (!fGeom)
+    {
+        *fLog << err << "MGeomCam not found... abort." << endl;
+        return kFALSE;
     }
 
@@ -173 +171 @@
 // --------------------------------------------------------------------------
 //
-// 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 MHVsSize::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;
-}
-
-// --------------------------------------------------------------------------
-//
-// With this function you can convert the histogram ('on the fly') between
-// degrees and millimeters.
-//
-void MHVsSize::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;
-    MH::ScaleAxis(&fLength, 1, scale);
-    MH::ScaleAxis(&fWidth,  1, scale);
-    MH::ScaleAxis(&fDist,   1, scale);
-    MH::ScaleAxis(&fM3Long, 1, scale);
-    MH::ScaleAxis(&fArea,   1, scale*scale);
-
-    if (mmscale)
-    {
-        fLength.SetYTitle("Length [mm]");
-        fWidth.SetYTitle("Width [mm]");
-        fDist.SetYTitle("Distance [mm]");
-        fArea.SetYTitle("Area [mm^{2}]");
-        fM3Long.SetYTitle("M3Long [mm]");
-    }
-    else
-    {
-        fLength.SetYTitle("Length [\\circ]");
-        fWidth.SetYTitle("Width [\\circ]");
-        fDist.SetYTitle("Distance [\\circ]");
-        fArea.SetYTitle("Area [\\circ^{2}]");
-        fM3Long.SetYTitle("M3Long [\\circ]");
-    }
-
-    fUseMmScale = mmscale;
-}
-
-// --------------------------------------------------------------------------
-//
 // Fill the histograms with data from a MHillas-Container.
 // Be careful: Only call this with an object of type MHillas
@@ -249 +183 @@
     }
 
-    const Double_t scale = fUseMmScale ? 1 : fMm2Deg;
+    const Double_t scale = fGeom->GetConvMm2Deg();
 
     fLength.Fill(fHillas->GetSize(), scale*fHillas->GetLength(),     w);

trunk/MagicSoft/Mars/mimage/MHVsSize.h

@@ -10 +10 @@
 #endif
 
+class MGeomCam;
 class MHillas;
 class MHillasExt;
@@ -17 +18 @@
 {
 private:
+    MGeomCam     *fGeom;        //! Conversion from mm to deg
     MHillas      *fHillas;      //!
     MHillasExt   *fHillasExt;   //!
@@ -30 +32 @@
     TH2F fArea;    // 
 
-    Float_t fMm2Deg;
-    Bool_t  fUseMmScale;
-
 public:
     MHVsSize(const char *name=NULL, const char *title=NULL);
-
-    void SetMmScale(Bool_t mmscale=kTRUE);
-    virtual void SetMm2Deg(Float_t mmdeg);
 
     Bool_t SetupFill(const MParList *pList);