Changeset 2307 for trunk/MagicSoft

Timestamp:

08/20/03 11:37:37 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (3 diffs)
• mhist/MHHadronness.cc (modified) (5 diffs)
• mhist/MHMatrix.cc (modified) (1 diff)
• mranforest/MRanTree.cc (modified) (4 diffs)
• mranforest/MRanTree.h (modified) (2 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -1 +1 @@
                                                  -*-*- END OF LINE -*-*-
+
+ 2003/08/20: Thomas Bretz
+
+   * mhist/MHCamEvent.[h,cc]:
+     - added PrintOutlayers
+     
+   * mhist/MHHadronness.cc:
+     - changed PrintOutput
+     - set some Grids in plots
+     
+   * mhist/MHMatrix.cc:
+     - removed obsolete call to MEvtLoop::SetName
+     - Changed name to something more intuitive
+     
+   * mranforest/MRanForest.cc, mranforest/MRanTree.h:
+     - removed ^M
+     
+   * mranforest/MRanTree.cc:
+     - removed obsolete check for gRandom
+     - replaced mvar*numdata by a new variable mn
+
 
 
@@ -10 +31 @@
 
 
+
  2003/08/19: Wolfgang Wittek
 
@@ -52 +74 @@
       - MHFindSignificance included
       - MHCT1Supercuts included
+
 
 

trunk/MagicSoft/Mars/mhist/MHHadronness.cc

@@ -285 +285 @@
     }
 
+    //*fLog << dbg << val1x << "/" << val1y << "  " << val2x << "/" << val2y << endl;
+
     return val2x-val1x == 0 ? 0 : val1y - (val2y-val1y)/(val2x-val1x) * (val1x-0.5);
 }
@@ -451 +453 @@
 {
     *fLog << all;
-    *fLog << "Hadronness histograms:" << endl;
-    *fLog << "---------------------" << endl;
+    *fLog << underline << GetDescriptor() << endl;
 
     if (fGraph->GetN()==0)
@@ -506 +507 @@
     pad->cd(2);
     gPad->SetBorderMode(0);
+    gPad->SetGridx();
+    gPad->SetGridy();
     fIntGhness->Draw();
     fIntPhness->Draw("same");
@@ -511 +514 @@
     pad->cd(3);
     gPad->SetBorderMode(0);
+    gPad->SetGridx();
+    gPad->SetGridy();
     fQfac->Draw("A*");
     gPad->Modified();
@@ -526 +531 @@
     pad->cd(4);
     gPad->SetBorderMode(0);
+    gPad->SetGridx();
+    gPad->SetGridy();
     fGraph->Draw("AC");
     gPad->Modified();

trunk/MagicSoft/Mars/mhist/MHMatrix.cc

@@ -633 +633 @@
 
     //MProgressBar bar;
-    MEvtLoop evtloop;
+    MEvtLoop evtloop("MHMatrix::Fill-EvtLoop");
     evtloop.SetParList(plist);
-    evtloop.SetName("EvtLoopMatrix");
     //evtloop.SetProgressBar(&bar);
-
 
     if (!evtloop.Eventloop())

trunk/MagicSoft/Mars/mranforest/MRanTree.cc

@@ -140 +140 @@
                               const TArrayF &winbag)
 {
-    if(!gRandom)
-    {
-        *fLog << err << dbginf << "gRandom not initialized... aborting." << endl;
-        return kFALSE;
-    }
-
     const Int_t nrnodes = fBestSplit.GetSize();
     const Int_t numdata = (nrnodes-1)/2;
@@ -182 +176 @@
     for(Int_t mt=0;mt<fNumTry;mt++)
     {
-        Int_t mvar=Int_t(gRandom->Rndm()*mdim);
+        const Int_t mvar=Int_t(gRandom->Rndm()*mdim);
+        const Int_t mn  = mvar*numdata;
 
         // Gini index = rrn/rrd+rln/rld
@@ -197 +192 @@
         for(Int_t nsp=ndstart;nsp<=ndend-1;nsp++)
         {
-            const Int_t &nc=datasort[mvar*numdata+nsp];
+            const Int_t &nc=datasort[mn+nsp];
             const Int_t &k=hadtrue[nc];
 
@@ -210 +205 @@
             wr[k]-=u;
 
-            if (datarang[mvar*numdata+nc]>=datarang[mvar*numdata+datasort[mvar*numdata+nsp+1]])
+            if (datarang[mn+nc]>=datarang[mn+datasort[mn+nsp+1]])
                 continue;
             if (TMath::Min(rrd,rld)<=1.0e-5)

trunk/MagicSoft/Mars/mranforest/MRanTree.h

@@ -37 +37 @@
     TArrayF fGiniDec;
 
+    Int_t FindBestSplit(const TArrayI &datasort, const TArrayI &datarang,
+                        const TArrayI &hadtrue,
+                        Int_t ndstart, Int_t ndend, TArrayF &tclasspop,
+                        Int_t &msplit, Float_t &decsplit, Int_t &nbest,
+                        const TArrayF &winbag);
+
+    void MoveData(TArrayI &datasort, Int_t ndstart, Int_t ndend,
+                  TArrayI &idmove, TArrayI &ncase, Int_t msplit,
+                  Int_t nbest, Int_t &ndendl);
+
+    void BuildTree(TArrayI &datasort, const TArrayI &datarang,
+                   const TArrayI &hadtrue,
+                   TArrayI &bestsplit,TArrayI &bestsplitnext,
+                   TArrayF &tclasspop,
+                   const TArrayF &winbag,
+                   Int_t ninbag);
+
 public:
     MRanTree(const char *name=NULL, const char *title=NULL);
@@ -66 +83 @@
                   TArrayF &tclasspop, TArrayI &jinbag, const TArrayF &winbag);
 
-    Int_t FindBestSplit(const TArrayI &datasort, const TArrayI &datarang,
-                        const TArrayI &hadtrue,
-                        Int_t ndstart, Int_t ndend, TArrayF &tclasspop,
-                        Int_t &msplit, Float_t &decsplit, Int_t &nbest,
-                        const TArrayF &winbag);
-
-    void MoveData(TArrayI &datasort, Int_t ndstart, Int_t ndend,
-                  TArrayI &idmove, TArrayI &ncase, Int_t msplit,
-                  Int_t nbest, Int_t &ndendl);
-
-    void BuildTree(TArrayI &datasort, const TArrayI &datarang,
-                   const TArrayI &hadtrue,
-                   TArrayI &bestsplit,TArrayI &bestsplitnext,
-                   TArrayF &tclasspop,
-                   const TArrayF &winbag,
-                   Int_t ninbag);
-
     Double_t TreeHad(const TVector &event);
     Double_t TreeHad(const TMatrixRow &event);