Changeset 8657 for trunk/MagicSoft/Mars/mimage

Timestamp:

08/06/07 15:44:16 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mimage

Files:

• MHHillasExt.cc (modified) (15 diffs)
• MHHillasExt.h (modified) (5 diffs)

trunk/MagicSoft/Mars/mimage/MHHillasExt.cc

@@ -18 +18 @@
 !   Author(s): Thomas Bretz, 2001 <mailto:tbretz@astro.uni-wuerzburg.de>
 !
-!   Copyright: MAGIC Software Development, 2000-2004
+!   Copyright: MAGIC Software Development, 2000-2007
 !
 !
@@ -28 +28 @@
 //
 // This class contains histograms for every Hillas parameter
+//
+// Class Version 2:
+// ----------------
+//  - fHMaxDist
+//  + fHSlopeL
 //
 /////////////////////////////////////////////////////////////////////////////
@@ -34 +39 @@
 #include <math.h>
 
-#include <TH1.h>
 #include <TPad.h>
 #include <TLegend.h>
@@ -47 +51 @@
 
 #include "MBinning.h"
+
+#include "MHillas.h"
 #include "MHillasExt.h"
 #include "MHillasSrc.h"
@@ -59 +65 @@
 //
 MHHillasExt::MHHillasExt(const char *name, const char *title)
-    : fMm2Deg(1), fUseMmScale(kTRUE), fHilName("MHillasExt")
+    : fHillas(0), fHillasExt(0), fMm2Deg(1),
+    fUseMmScale(kTRUE), fHilName("MHillasExt")
 {
     //
@@ -75 +82 @@
     fHM3Long.UseCurrentStyle();
     fHM3Trans.UseCurrentStyle();
-    fHMaxDist.UseCurrentStyle();
+    fHSlopeL.UseCurrentStyle();
 
     fHAsym.SetName("Asymmetry");
     fHM3Long.SetName("M3l");
     fHM3Trans.SetName("M3t");
-    fHMaxDist.SetName("MaxDist");
+    fHSlopeL.SetName("SlopeL");
 
     fHAsym.SetTitle("Asymmetry");
     fHM3Long.SetTitle("3^{rd} Moment Longitudinal");
     fHM3Trans.SetTitle("3^{rd} Moment Transverse");
-    fHMaxDist.SetTitle("Distance of max distant pixel");
+    fHSlopeL.SetTitle("Longitudinal time-slope vs. Dist");
 
     fHAsym.SetXTitle("Asym [mm]");
     fHM3Long.SetXTitle("3^{rd} M_{l} [mm]");
     fHM3Trans.SetXTitle("3^{rd} M_{t} [mm]");
-    fHMaxDist.SetXTitle("D_{max} [mm]");
+    fHSlopeL.SetXTitle("D [mm]");
 
     fHAsym.SetYTitle("Counts");
     fHM3Long.SetYTitle("Counts");
     fHM3Trans.SetYTitle("Counts");
-    fHMaxDist.SetYTitle("Counts");
+    fHSlopeL.SetYTitle("S_{l} [ns/mm]");
 
     fHAsym.SetFillStyle(4000);
     fHM3Long.SetFillStyle(4000);
     fHM3Trans.SetFillStyle(4000);
-    fHMaxDist.SetFillStyle(4000);
+    //fHSlopeL.SetFillStyle(4000);
 
     fHAsym.SetDirectory(NULL);
     fHM3Long.SetDirectory(NULL);
     fHM3Trans.SetDirectory(NULL);
-    fHMaxDist.SetDirectory(NULL);
+    fHSlopeL.SetDirectory(NULL);
 
     fHM3Trans.SetLineColor(kBlue);
 
-    MBinning bins;
-
-    bins.SetEdges(51, -326, 326);
-    bins.Apply(fHM3Long);
-    bins.Apply(fHM3Trans);
-
-    bins.SetEdges(51, -593, 593);
-    bins.Apply(fHAsym);
-
-    bins.SetEdges(100, 0, 593);
-    bins.Apply(fHMaxDist);
+    MBinning binsx, binsy;
+
+    binsx.SetEdges(51, -326, 326);
+    binsx.Apply(fHM3Long);
+    binsx.Apply(fHM3Trans);
+
+    binsx.SetEdges(51, -593, 593);
+    binsx.Apply(fHAsym);
+
+    binsx.SetEdges(100,     0,  445);
+    binsy.SetEdges(100, -0.04, 0.04);
+    MH::SetBinning(&fHSlopeL, &binsx, &binsy);
 }
 
@@ -141 +149 @@
     }
 
+    fHillas = (MHillas*)plist->FindObject("MHillas");
+    if (!fHillas)
+    {
+        *fLog << err << "MHillas not found in parameter list... aborting." << endl;
+        return kFALSE;
+    }
+
     const MGeomCam *geom = (MGeomCam*)plist->FindObject("MGeomCam");
     if (!geom)
@@ -150 +165 @@
     }
 
-    ApplyBinning(*plist, "Asym",    &fHAsym);
-    ApplyBinning(*plist, "M3Long",  &fHM3Long);
+    ApplyBinning(*plist, "Asym", &fHAsym);
+    ApplyBinning(*plist, "M3Long", &fHM3Long);
     ApplyBinning(*plist, "M3Trans", &fHM3Trans);
-    ApplyBinning(*plist, "MaxDist", &fHMaxDist);
+    ApplyBinning(*plist, "Dist", "Slope", &fHSlopeL);
 
     return kTRUE;
@@ -168 +183 @@
 
     const Double_t scale = TMath::Sign(fUseMmScale?1:fMm2Deg, (src ? src->GetCosDeltaAlpha() : 1));
+    const Double_t dist  = src ? src->GetDist() : fHillas->GetDist0();
 
     fHAsym.Fill(scale*fHillasExt->GetAsym(), w);
     fHM3Long.Fill(scale*fHillasExt->GetM3Long(), w);
     fHM3Trans.Fill(scale*fHillasExt->GetM3Trans(), w);
-    fHMaxDist.Fill(TMath::Abs(scale*fHillasExt->GetMaxDist()), w);
+    fHSlopeL.Fill(scale*dist, fHillasExt->GetSlopeLong()/scale, w);
 
     return kTRUE;
@@ -197 +213 @@
     MH::ScaleAxis(&fHM3Long,  scale);
     MH::ScaleAxis(&fHM3Trans, scale);
-    MH::ScaleAxis(&fHMaxDist, scale);
+    MH::ScaleAxis(&fHSlopeL,  scale, 1./scale);
 
     if (mmscale)
@@ -204 +220 @@
         fHM3Long.SetXTitle("3^{rd} M_{l} [mm]");
         fHM3Trans.SetXTitle("3^{rd} M_{t} [mm]");
-        fHMaxDist.SetXTitle("D_{max} [mm]");
+        fHSlopeL.SetXTitle("D [mm]");
+        fHSlopeL.SetYTitle("S_{l} [ns/mm]");
     }
     else
@@ -211 +228 @@
         fHM3Long.SetXTitle("3^{rd} M_{l} [\\circ]");
         fHM3Trans.SetXTitle("3^{rd} M_{t} [\\circ]");
-        fHMaxDist.SetXTitle("D_{max} [\\circ]");
+        fHSlopeL.SetXTitle("D [\\circ]");
+        fHSlopeL.SetYTitle("S_{l} [ns/\\circ]");
     }
 
@@ -262 +280 @@
         fHM3Long.SetName("M3lSame");
         fHM3Trans.SetName("M3tSame");
-        fHMaxDist.SetName("MaxDistSame");
+        fHSlopeL.SetName("SlopeLSame");
 
         fHAsym.SetDirectory(0);
         fHM3Long.SetDirectory(0);
         fHM3Trans.SetDirectory(0);
-        fHMaxDist.SetDirectory(0);
+        fHSlopeL.SetDirectory(0);
 
         fHM3Long.SetLineColor(kMagenta);
         fHM3Trans.SetLineColor(kCyan);
         fHAsym.SetLineColor(kBlue);
-        fHMaxDist.SetLineColor(kBlue);
+        fHSlopeL.SetMarkerColor(kBlue);
     }
 
@@ -288 +306 @@
     pad->cd(2);
     gPad->SetBorderMode(0);
-    RemoveFromPad("MaxDistSame");
-    fHMaxDist.Draw(same?"same":"");
+    //RemoveFromPad("SlopeLSame");
+    //fHSlopeL.Draw(same?"same":"");
+    if (same)
+    {
+        TH2 *h=dynamic_cast<TH2*>(gPad->FindObject("SlopeL"));
+        if (h)
+        {
+            // This causes crashes in THistPainter::PaintTable
+            // if the z-axis is not kept. No idea why...
+            h->SetDrawOption("z");
+            h->SetMarkerColor(kBlack);
+        }
+
+        RemoveFromPad("SlopeLSame");
+        fHSlopeL.SetMarkerColor(kGreen);
+        fHSlopeL.Draw("same");
+    }
+    else
+        fHSlopeL.Draw("colz");
 
     delete pad->GetPad(4);
@@ -302 +337 @@
     if (name.Contains("M3T", TString::kIgnoreCase))
         return const_cast<TH1F*>(&fHM3Trans);
-    if (name.Contains("MaxDist", TString::kIgnoreCase))
-        return const_cast<TH1F*>(&fHMaxDist);
+    if (name.Contains("SlopeL", TString::kIgnoreCase))
+        return const_cast<TH2F*>(&fHSlopeL);
 
     return NULL;

trunk/MagicSoft/Mars/mimage/MHHillasExt.h

@@ -2 +2 @@
 #define MARS_MHHillasExt
 
-#ifndef ROOT_TH1
-#include <TH1.h>
+#ifndef ROOT_TH2
+#include <TH2.h>
 #endif
 #ifndef MARS_MH
@@ -9 +9 @@
 #endif
 
+class MHillas;
 class MHillasExt;
 
@@ -14 +15 @@
 {
 private:
+    MHillas    *fHillas;    //! Pointer to the MHillas container
     MHillasExt *fHillasExt; //! Pointer to the MHillasExt container
 
@@ -19 +21 @@
     TH1F fHM3Long;  // [mm]    3rd moment (e-weighted) along major axis
     TH1F fHM3Trans; // [mm]    3rd moment (e-weighted) along minor axis
-    TH1F fHMaxDist; // [mm]    Distance between shower center maximum distant pixel
+    TH2F fHSlopeL;  //
 
     Float_t fMm2Deg;
@@ -46 +48 @@
     void Draw(Option_t *opt=NULL);
 
-    ClassDef(MHHillasExt, 1) // Container which holds histograms for the extended hillas parameters
+    ClassDef(MHHillasExt, 2) // Container which holds histograms for the extended hillas parameters
 };