Changeset 2133 for trunk/MagicSoft/Mars

Timestamp:

05/22/03 18:27:22 (22 years ago)

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mhist/MHMatrix.cc (modified) (1 diff)
• mhist/MHMatrix.h (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -1 +1 @@
                                                  -*-*- END OF LINE -*-*-
 
-
+ 2003/05/22: Abelardo Moralejo
+
+   * mhist/MHMatrix.[h,cc]
+     - add member function ShuffleRows() to randomize the order of the
+       matrix rows. This is useful for instance for the random forest, 
+       (See RanForest.C) when we feed a hadron training sample with 
+       both protons and helium nuclei: if they are ordered (first all
+       events of one type, then those of the other) the method does not
+       seem to work well. Any other kind of ordering might be harmful 
+       as well (in theta, phi or whatever).
 
  2003/05/22: Wolfgang Wittek

trunk/MagicSoft/Mars/mhist/MHMatrix.cc

@@ -1083 +1083 @@
     return i!=0;
 }
+
+// --------------------------------------------------------------------------
+//
+// ShuffleEvents. Shuffles the order of the matrix rows.
+// 
+//
+void MHMatrix::ShuffleRows(UInt_t seed)
+{
+  TRandom rnd(seed);
+
+  for (Int_t irow = 0; irow < fNumRow; irow++)
+    {
+      Int_t jrow = rnd.Integer(fNumRow);
+
+      for (Int_t icol = 0; icol < fM.GetNcols(); icol++)
+        {
+          Real_t tmp = fM(irow,icol);
+          fM(irow,icol) = fM(jrow,icol);
+          fM(jrow,icol) = tmp;
+        }
+    }
+
+  *fLog << warn << this->GetName() << " : Attention! Matrix rows have been shuffled." << endl;
+
+}

trunk/MagicSoft/Mars/mhist/MHMatrix.h

@@ -112 +112 @@
     Bool_t ReadEnv(const TEnv &env, TString prefix, Bool_t print=kFALSE);
 
+    void ShuffleRows(UInt_t seed);
+
     ClassDef(MHMatrix, 1) // Multidimensional Matrix to store events
 };