Changeset 717 for trunk/MagicSoft/Mars/mhist

Timestamp:

04/03/01 13:28:36 (24 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHStarMap.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mhist/MHStarMap.cc

@@ -11 +11 @@
 #include "MHStarMap.h"
 
-#include <TH2.h>
-#include <TStyle.h>
-#include <TCanvas.h>
+#include <TH2.h>      // TH2F
+#include <TStyle.h>   // gStyle
+#include <TColor.h>   // SetRGB
+#include <TCanvas.h>  // TCanvas
 
 #include "MHillas.h"
@@ -39 +40 @@
     //
     fStarMap = new TH2F("Star Map", "Counts",
-                        50, -300, 300,
-                        50, -300, 300);
+                        150, -300, 300,
+                        150, -300, 300);
 }
 
@@ -50 +51 @@
 void MHStarMap::Draw(Option_t *)
 {
+    //
+    // Creates a new canvas, creates a useful palette and
+    // draws the histogram in the new created canvas
+    //
+
     TCanvas *c = new TCanvas("Star Map", "Star Map created from Hillas Parameters", 500, 500);
 
-    gStyle->SetPalette(1, 0);
+    //
+    // Set the palette you wanna use:
+    //  - you could set the root "Pretty Palette Violet->Red" by
+    //    gStyle->SetPalette(1, 0), but in some cases this may look
+    //    confusing
+    //  - The maximum colors root allowes us to set by ourself
+    //    is 50 (idx: 51-100). This colors are set to a grayscaled
+    //    palette
+    //  - the number of contours must be two less than the number
+    //    of palette entries
+    //
+
+    const Int_t numg = 32; // number of gray scaled colors
+    const Int_t numw = 32; // number of white
+
+    Int_t palette[numg+numw];
+
+    //
+    // The first half of the colors are white.
+    // This is some kind of optical background supression
+    //
+    gROOT->GetColor(51)->SetRGB(1, 1, 1);
+
+    Int_t i;
+    for (i=0; i<numw; i++)
+        palette[i] = 51;
+
+    //
+    // now the (gray) scaled part is coming
+    //
+    for (;i<numw+numg; i++)
+    {
+        const Float_t gray = 1.0-(float)(i-numw)/(numg-1.0);
+
+        gROOT->GetColor(52+i)->SetRGB(gray, gray, gray);
+        palette[i] = 52+i;
+    }
+
+    //
+    // Set the palette and the number of contour levels
+    //
+    gStyle->SetPalette(numg+numw, palette);
+    fStarMap->SetContour(numg+numw-2);
+
+    // gStyle->SetPalette(1, 0);
     fStarMap->Draw("colz");
 
@@ -68 +118 @@
     const float t = dist*(sin(theta)-cos(theta)*m);
 
-    float y0 = -m*300+t;
-    for (int i=-297; i<300; i+=3)
-    {
-        const float y1 = m*i+t;
+    if (m>-1 && m<1)
+        for (int x=-298; x<298; x+=4)
+        {
+            const float y = m*x+t;
 
-//        if (y1<=y0-3 || y1>=y0+3)
-//            continue;
+            fStarMap->Fill(x, y);
+        }
+    else
+        for (int y=-298; y<298; y+=4)
+        {
+            const float x = (y-t)/m;
 
-        fStarMap->Fill(i, y1);
-        y0 = y1;
-    }
+            fStarMap->Fill(x, y);
+        }
 }