Changeset 1410 for trunk/MagicSoft/Mars/macros/plot.C

Timestamp:

07/16/02 13:24:59 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/plot.C (modified) (6 diffs)

trunk/MagicSoft/Mars/macros/plot.C

@@ -24 +24 @@
 \* ======================================================================== */
 
+// -------------------------------------------------------------------------
+//
+//  plot.C
+//
+//  This macro shows how to fill and display a histogram using Mars
+//
 void plot()
 {
-    //
-    // This is a demonstration program which calculates the Hillas
-    // parameter out of a Magic root file.
-
     //
     // Create a empty Parameter List and an empty Task List
@@ -43 +45 @@
     // ---------------------------------
     //
-    // The first argument is the tree you want to read.
-    //   Events:     Cosmic ray events
-    //   PedEvents:  Pedestal Events
-    //   CalEvents:  Calibration Events
-    //
+
+    // First Task: Read file with image parameters
+    // (created with the star.C macro)
     MReadMarsFile  read("Events", "star.root");
     read.DisableAutoScheme();
     tlist.AddToList(&read);
 
+    // Create a filter for Gammas
     MFParticleId fgamma("MMcEvt", '=', kGAMMA);
     tlist.AddToList(&fgamma);
 
+    // Create a filter for Non-Gammas
     MFParticleId fhadrons("MMcEvt", '!', kGAMMA);
     tlist.AddToList(&fhadrons);
@@ -66 +68 @@
     MBinning bins("BinningMH3X");
     bins.SetEdgesLog(50, 100, 20000);
+    plist.AddToList(&bins);
     //
-    // -------------------------------------------------
+    // -------------------------------------------------------
 
+    // Create a histogram for the data from gammas and from non-gammas
     MH3 h3g(var);
     MH3 h3h(var);
 
+    // Add the histograms to the parameter container list
     plist.AddToList(&h3g);
     plist.AddToList(&h3h);
 
-    plist.AddToList(&bins);
-
+    // Create a task which fills one histogram with the gamma-data
     MFillH fillg(&h3g);
     fillg.SetFilter(&fgamma);
     tlist.AddToList(&fillg);
 
+    // Create a task which fills the other histogram with the non-gamma-data
     MFillH fillh(&h3h);
     fillh.SetFilter(&fhadrons);
@@ -99 +104 @@
     tlist.PrintStatistics();
 
+    // Create a default canvas called Plot and set the
+    // x-axis to logarithmic scale
     MH::MakeDefCanvas("Plot");
+    gPad->SetLogx();
+
+    // Setup some style options of the two histograms
+    // and draw a copy of both
     h3h.GetHist().SetLineColor(kRed);
     h3h.GetHist().SetFillStyle(4000);
@@ -105 +116 @@
     h3h.GetHist().DrawCopy("same");
 
+    // Now create a new histogram, fill it with the division of the
+    // two histograms and draw also a copy of it
     TH1D h;
     MH::SetBinning(&h, &bins);
@@ -111 +124 @@
     h.SetFillStyle(4000);
     h.DrawCopy("same");
-
-    gPad->SetLogx();
 }