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

Timestamp:

07/25/02 13:58:19 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/mhist/MHHillas.cc

@@ -17 +17 @@
 !
 !   Author(s): Thomas Bretz  2001 <mailto:tbretz@uni-sw.gwdg.de>
+!              Wolfgang Wittek  2002 <mailto:wittek@mppmu.mpg.de>
 !
 !   Copyright: MAGIC Software Development, 2000-2002
@@ -27 +28 @@
 // MHHillas
 //
-// This class contains histograms for every Hillas parameter
+// This class contains histograms for the source independent image parameters
 //
 /////////////////////////////////////////////////////////////////////////////
@@ -36 +37 @@
 
 #include <TH1.h>
+#include <TH2.h>
 #include <TPad.h>
+#include <TStyle.h>
 #include <TCanvas.h>
 
@@ -42 +45 @@
 #include "MLogManip.h"
 
+#include "MParList.h"
+
+#include "MHillas.h"
 #include "MGeomCam.h"
-#include "MHillas.h"
-#include "MParList.h"
+#include "MBinning.h"
+
 
 ClassImp(MHHillas);
@@ -53 +59 @@
 //
 MHHillas::MHHillas(const char *name, const char *title)
-    : fMm2Deg(-1), fUseMmScale(kTRUE)
+    : fMm2Deg(1), fUseMmScale(kFALSE)
 {
     //
@@ -59 +65 @@
     //
     fName  = name  ? name  : "MHHillas";
-    fTitle = title ? title : "Container for Hillas histograms";
-
-    //
-    // loop over all Pixels and create two histograms
-    // one for the Low and one for the High gain
-    // connect all the histogram with the container fHist
-    //
+    fTitle = title ? title : "Source independent image parameters";
+
+    fLength = new TH1F("Length", "Length of Ellipse", 100, 0, 300);
+    fLength->SetDirectory(NULL);
+    fLength->SetXTitle("Length [mm]");
+    fLength->SetYTitle("Counts");
+
     fWidth  = new TH1F("Width",  "Width of Ellipse",  100, 0, 300);
-    fLength = new TH1F("Length", "Length of Ellipse", 100, 0, 300);
-
-    fLength->SetDirectory(NULL);
     fWidth->SetDirectory(NULL);
-
-    fLength->SetXTitle("Length [mm]");
     fWidth->SetXTitle("Width [mm]");
-
-    fLength->SetYTitle("Counts");
     fWidth->SetYTitle("Counts");
-}
-
-// --------------------------------------------------------------------------
-//
-// Delete the four histograms
+
+    fDistC  = new TH1F("DistC",  "Distance from center of camera",  100, 0, 600);
+    fDistC->SetDirectory(NULL);
+    fDistC->SetXTitle("Distance [mm]");
+    fDistC->SetYTitle("Counts");
+
+    fDelta  = new TH1F("Delta",  "Angle (Main axis - x-axis)", 100, -90, 90);
+    fDelta->SetDirectory(NULL);
+    fDelta->SetXTitle("Delta [\\circ]");
+    fDelta->SetYTitle("Counts");
+
+    MBinning bins;
+    bins.SetEdgesLog(50, 1, 1e7);
+
+    fSize  = new TH1F;
+    fSize->SetName("Size");
+    fSize->SetTitle("Number of Photons");
+    fSize->SetDirectory(NULL);
+    fSize->SetXTitle("Size");
+    fSize->SetYTitle("Counts");
+    fSize->GetXaxis()->SetTitleOffset(1.2);
+    fSize->GetXaxis()->SetLabelOffset(-0.015);
+
+    MH::SetBinning(fSize, &bins);
+
+    fCenter = new TH2F("Center", "Center of Ellipse", 60, -600, 600, 60, -600, 600);
+    fCenter->SetDirectory(NULL);
+    fCenter->SetXTitle("x [mm]");
+    fCenter->SetYTitle("y [mm]");
+    fCenter->SetZTitle("Counts");
+}
+
+// --------------------------------------------------------------------------
+//
+// Delete the histograms
 //
 MHHillas::~MHHillas()
 {
+    delete fLength;
     delete fWidth;
-    delete fLength;
+
+    delete fDistC;
+    delete fDelta;
+
+    delete fSize;
+    delete fCenter;
 }
 
@@ -101 +136 @@
 Bool_t MHHillas::SetupFill(const MParList *plist)
 {
-    const MBinning* binsw = (MBinning*)plist->FindObject("BinningWidth");
-    const MBinning* binsl = (MBinning*)plist->FindObject("BinningLength");
-    if (!binsw || !binsl)
+    const MBinning *binsw = (MBinning*)plist->FindObject("BinningWidth");
+    const MBinning *binsl = (MBinning*)plist->FindObject("BinningLength");
+    const MBinning *binsd = (MBinning*)plist->FindObject("BinningDist");
+    const MBinning *binsc = (MBinning*)plist->FindObject("BinningCamera");
+
+    if (!binsw || !binsl || !binsd || !binsc)
     {
         *fLog << err << dbginf << "At least one MBinning not found... aborting." << endl;
@@ -111 +149 @@
     SetBinning(fWidth,  binsw);
     SetBinning(fLength, binsl);
+    SetBinning(fDistC,  binsd);
+    SetBinning(fCenter, binsc, binsc);
 
     const MGeomCam *geom = (MGeomCam*)plist->FindObject("MGeomCam");
@@ -119 +159 @@
     }
 
+    fMm2Deg     = geom->GetConvMm2Deg();
+    fUseMmScale = kFALSE;
+
     fLength->SetXTitle("Length [\\circ]");
     fWidth->SetXTitle("Width [\\circ]");
-
-    fMm2Deg     = geom->GetConvMm2Deg();
-    fUseMmScale = kFALSE;
-
-    return kTRUE;
-}
-
-// --------------------------------------------------------------------------
-//
-// Fill the four histograms with data from a MHillas-Container.
-// Be careful: Only call this with an object of type MHillas
-//
-Bool_t MHHillas::Fill(const MParContainer *par)
-{
-    const MHillas &h = *(MHillas*)par;
-
-    if (fUseMmScale)
-    {
-        fWidth ->Fill(h.GetWidth());
-        fLength->Fill(h.GetLength());
-    }
-    else
-    {
-        fWidth ->Fill(fMm2Deg*h.GetWidth());
-        fLength->Fill(fMm2Deg*h.GetLength());
-    }
+    fDistC->SetXTitle("Distance [\\circ]");
+    fCenter->SetXTitle("x [\\circ]");
+    fCenter->SetYTitle("y [\\circ]");
 
     return kTRUE;
@@ -190 +210 @@
     if (fUseMmScale)
     {
+        MH::ScaleAxis(fLength, 1./fMm2Deg);
+        MH::ScaleAxis(fWidth,  1./fMm2Deg);
+        MH::ScaleAxis(fDistC,  1./fMm2Deg);
+        MH::ScaleAxis(fCenter, 1./fMm2Deg, 1./fMm2Deg);
+
         fLength->SetXTitle("Length [mm]");
         fWidth->SetXTitle("Width [mm]");
-
-        fLength->Scale(1./fMm2Deg);
-        fWidth->Scale(1./fMm2Deg);
+        fCenter->SetXTitle("x [mm]");
+        fCenter->SetYTitle("y [mm]");
     }
     else
     {
+        MH::ScaleAxis(fLength, fMm2Deg);
+        MH::ScaleAxis(fWidth,  fMm2Deg);
+        MH::ScaleAxis(fDistC,  fMm2Deg);
+        MH::ScaleAxis(fCenter, fMm2Deg, fMm2Deg);
+
         fLength->SetXTitle("Length [\\circ]");
         fWidth->SetXTitle("Width [\\circ]");
-
-        fLength->Scale(fMm2Deg);
-        fWidth->Scale(fMm2Deg);
+        fDistC->SetXTitle("Distance [\\circ]");
+        fCenter->SetXTitle("x [\\circ]");
+        fCenter->SetYTitle("y [\\circ]");
     }
 
@@ -210 +239 @@
 // --------------------------------------------------------------------------
 //
-// Draw clones of all four histograms. So that the object can be deleted
+// Fill the histograms with data from a MHillas-Container.
+// Be careful: Only call this with an object of type MHillas
+//
+Bool_t MHHillas::Fill(const MParContainer *par)
+{
+    const MHillas &h = *(MHillas*)par;
+
+    const Double_t d = sqrt(h.GetMeanX()*h.GetMeanX() + h.GetMeanY()*h.GetMeanY());
+
+    if (fUseMmScale)
+    {
+        fLength->Fill(h.GetLength());
+        fWidth ->Fill(h.GetWidth());
+        fDistC ->Fill(d);
+        fCenter->Fill(h.GetMeanX(), h.GetMeanY());
+    }
+    else
+    {
+        fLength->Fill(fMm2Deg*h.GetLength());
+        fWidth ->Fill(fMm2Deg*h.GetWidth());
+        fDistC ->Fill(fMm2Deg*d);
+        fCenter->Fill(fMm2Deg*h.GetMeanX(), fMm2Deg*h.GetMeanY());
+    }
+
+    fDelta->Fill(kRad2Deg*h.GetDelta());
+    fSize->Fill(h.GetSize());
+
+    return kTRUE;
+}
+
+// --------------------------------------------------------------------------
+//
+// Setup a inversed deep blue sea palette for the fCenter histogram.
+//
+void MHHillas::SetColors() const
+{
+    gStyle->SetPalette(51, NULL);
+    Int_t c[50];
+    for (int i=0; i<50; i++)
+        c[49-i] = gStyle->GetColorPalette(i);
+    gStyle->SetPalette(50, c);
+}
+
+// --------------------------------------------------------------------------
+//
+// Draw clones of four histograms. So that the object can be deleted
 // and the histograms are still visible in the canvas.
 // The cloned object are deleted together with the canvas if the canvas is
@@ -218 +292 @@
 TObject *MHHillas::DrawClone(Option_t *opt) const
 {
-    TCanvas *c = MakeDefCanvas("Hillas", "Histograms of Hillas Parameters",
-                               350, 500);
-    c->Divide(1, 2);
+    TCanvas *c = MakeDefCanvas("Hillas", fTitle, 700, 750);
+    c->Divide(2,3);
 
     gROOT->SetSelectedPad(NULL);
 
-    //
-    // This is necessary to get the expected bahviour of DrawClone
-    //
     c->cd(1);
     fLength->DrawCopy();
@@ -232 +302 @@
     c->cd(2);
     fWidth->DrawCopy();
+
+    c->cd(3);
+    gPad->SetLogx();
+    fSize->DrawCopy();
+
+    c->cd(4);
+    fDelta->DrawCopy();
+
+    c->cd(5);
+    fDistC->DrawCopy();
+
+    c->cd(6);
+    SetColors();
+    fCenter->DrawCopy("colz");
 
     c->Modified();
@@ -248 +332 @@
 {
     if (!gPad)
-        MakeDefCanvas("Hillas", "Histograms of Hillas Parameters", 350, 500);
-
-    gPad->Divide(1, 2);
+        MakeDefCanvas("Hillas", fTitle, 700, 750);
+
+    gPad->Divide(2,3);
 
     gPad->cd(1);
@@ -258 +342 @@
     fWidth->Draw();
 
+    gPad->cd(3);
+    gPad->SetLogx();
+    fSize->Draw();
+
+    gPad->cd(4);
+    fDelta->Draw();
+
+    gPad->cd(5);
+    fDistC->Draw();
+
+    gPad->cd(6);
+    SetColors();
+    fCenter->Draw("colz");
+
     gPad->Modified();
     gPad->Update();