Changeset 2114 for trunk/MagicSoft/Mars/mranforest

Timestamp:

05/13/03 15:16:54 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mranforest

Files:

• MRanForest.cc (modified) (2 diffs)
• MRanForest.h (modified) (1 diff)
• MRanForestCalc.cc (modified) (1 diff)
• MRanTree.cc (modified) (5 diffs)
• MRanTree.h (modified) (1 diff)

trunk/MagicSoft/Mars/mranforest/MRanForest.cc

@@ -85 +85 @@
 void MRanForest::SetPriors(Float_t prior_had, Float_t prior_gam)
 {
-    Float_t sum;
-
-    sum=prior_gam+prior_had;
+    const Float_t sum=prior_gam+prior_had;
 
     prior_gam/=sum;
@@ -98 +96 @@
 }
 
-Double_t MRanForest::CalcHadroness(TVector &event)
+Double_t MRanForest::CalcHadroness(const TVector &event)
 {
     Double_t hadroness=0;
     Int_t ntree=0;
+
+    TIter forest(fForest);
+
     MRanTree *tree;
-
-    TIter forest(fForest);
-    forest.Reset();
-
     while ((tree=(MRanTree*)forest.Next()))
     {

trunk/MagicSoft/Mars/mranforest/MRanForest.h

@@ -96 +96 @@
  
     // use forest to calculate hadronness of event
-    Double_t CalcHadroness(TVector &event);
+    Double_t CalcHadroness(const TVector &event);
 
     Bool_t AsciiWrite(ostream &out) const;

trunk/MagicSoft/Mars/mranforest/MRanForestCalc.cc

@@ -137 +137 @@
 Bool_t MRanForestCalc::Process()
 {
-    const Double_t ncols = fData->GetNumEntries();
-    TVector event(ncols);
-
-    for (int i=0; i<fData->GetNumEntries(); i++)
-        event(i) = (*fData)(i);
+    // first copy the data from the data array to a vector event
+    TVector event;
+    *fData >> event;
 
     Double_t hadroness=fRanForest->CalcHadroness(event);

trunk/MagicSoft/Mars/mranforest/MRanTree.cc

@@ -178 +178 @@
     // (Gini-index as criterium crit)
 
-    critmax=-1.0e20;
+    critmax=-1.0e20;  // FIXME: Replace by a constant from limits.h
 
     // random split selection, number of trials = fNumTry
@@ -301 +301 @@
     for(Int_t n=ndstart;n<=ndend;n++)
         ncase[n]=datasort[msplit*numdata+n];
-
-    return;
 }
 
@@ -459 +457 @@
             fBestSplit[kn]/=(classpop[0*nrnodes+kn]+classpop[1*nrnodes+kn]);
         }
-
-    return;
 }
 
@@ -468 +464 @@
 }
 
-Double_t MRanTree::TreeHad(TVector &event)
-{
-    Int_t kt=0;
-    // to optimize on storage space node status and node class
-    // are coded into fBestVar:
-    // status of node kt = TMath::Sign(1,fBestVar[kt])
-    // hadronness assigned to node kt = fBestSplit[kt]
-
-    for (Int_t k=0;k<fNumNodes;k++)
-    {
-        if (fBestVar[kt]<0)
-            break;
-
-        Int_t m=fBestVar[kt];
-
-        if (event(m)<=fBestSplit[kt])
-            kt=fTreeMap1[kt];
-        else
-            kt=fTreeMap2[kt];
-    }
-
-    return fBestSplit[kt];
-}
-
-Double_t MRanTree::TreeHad()
-{
-    const Double_t ncols = fData->GetNumEntries();
-    TVector event(ncols);
-
-    for (int i=0; i<fData->GetNumEntries(); i++)
-        event(i) = (*fData)(i);
-
+Double_t MRanTree::TreeHad(const TVector &event)
+{
     Int_t kt=0;
     // to optimize on storage space node status and node class
@@ -513 +479 @@
             break;
 
-        Int_t m=fBestVar[kt];
-
-        if (event(m)<=fBestSplit[kt])
-            kt=fTreeMap1[kt];
-        else
-            kt=fTreeMap2[kt];
-
+        const Int_t m=fBestVar[kt];
+
+        kt = event(m)<=fBestSplit[kt] ? fTreeMap1[kt] : fTreeMap2[kt];
     }
 
     return fBestSplit[kt];
+}
+
+Double_t MRanTree::TreeHad()
+{
+    TVector event;
+    *fData >> event;
+
+    return TreeHad(event);
 }
 

trunk/MagicSoft/Mars/mranforest/MRanTree.h

@@ -81 +81 @@
                    TArrayF &wl,Int_t ninbag);
 
-    Double_t TreeHad(TVector &event);
+    Double_t TreeHad(const TVector &event);
     Double_t TreeHad();