Changeset 1609 for trunk/MagicSoft/Mars/macros/multidimdist.C

Timestamp:

11/14/02 12:15:13 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/multidimdist.C (modified) (3 diffs)

trunk/MagicSoft/Mars/macros/multidimdist.C

@@ -40 +40 @@
 
     //
-    // Now setup the tasks and tasklist:
-    // ---------------------------------
+    // This is an example program which reads image parameters for gammas 
+    // and hadrons, builds a 'matrix' of look-alike events,
+    // and then finds quality numbers and acceptances for the same sample
+
+    // Create an empty Parameter List and an empty Task List
+    // The tasklist is identified in the eventloop by its name
     //
-    // The first argument is the tree you want to read.
-    //   Events:     Cosmic ray events
-    //   PedEvents:  Pedestal Events
-    //   CalEvents:  Calibration Events
+    MParList  plist;
+    MTaskList tlistg;
+    plist.AddToList(&tlistg);
+
     //
-    MReadMarsFile  read("Events", "~/Mars/star.root");
+    // ------------- user attention -----------------
+    //
+
+    // cut used both for the matrices and the sample
+    // for more information see the documentation of MF and MDataChain
+    MF filterenergy("MMcEvt.fEnergy > 0");
+
+    // matrix limitation for look-alike events (approximate number)
+    MFEventSelector selector;
+    selector.SetNumSelectEvts(2000);
+
+    // setup an AND-Filterlist from the two filters to be used
+    // in the event selection for the filling of the matrices
+    MFilterList flist;
+    flist.AddToList(&filterenergy);
+    flist.AddToList(&selector);
+
+    //
+    // ---------------------------------------------------------------
+    //  Now set up the tasks and tasklist (first event loop, gammas)
+    // ---------------------------------------------------------------
+    //
+
+    // --------------- user change -----------------
+    //  Give the names of the star-files to be read
+    //   Here you give the trainings sample(s) for
+    //                 the hadrons
+    // ---------------------------------------------
+    MReadMarsFile read("Events");
+    read.AddFile("star_gammas.root");
+    read.AddFile("star_protons.root");
     read.DisableAutoScheme();
-    tlist.AddToList(&read);
+    tlistg.AddToList(&read);
 
     MFParticleId fgamma("MMcEvt", '=', kGAMMA);
@@ -65 +99 @@
     matrix.AddColumn("abs(MHillas.fM3Long)");
     matrix.AddColumn("abs(MHillas.fM3Trans)");
-    matrix.AddColumn("abs(HillasSource.fHeadTail)");
+    matrix.AddColumn("abs(MHillasSrc.fHeadTail)");
     matrix.AddColumn("MHillas.fConc");
     matrix.AddColumn("MHillas.fConc1");
-    matrix.AddColumn("HillasSource.fDist");
+    matrix.AddColumn("MHillasSrc.fDist");
     matrix.AddColumn("log10(MHillas.fSize)");
     plist.AddToList(&matrix);
@@ -107 +141 @@
     plist.Replace(&tlist2);
 
-    MReadMarsFile read2("Events", "~/Mars/star.root");
-//    read2.AddFile("~/Mars/star2.root");
+    MReadMarsFile read2("Events");
+    read2.("gammas2.root");
+    read2.AddFile("hadrons2.root");
     read2.DisableAutoScheme();
     tlist2.AddToList(&read2);
 
     MMultiDimDistCalc calc;
-    calc.SetUseNumRows(15);
-    //calc.SetUseKernelMethod(kTRUE);
+    calc.SetUseNumRows(0);
+    calc.SetUseKernelMethod(kTRUE);
     tlist2.AddToList(&calc);