Changeset 7670 for trunk/MagicSoft/Mars/mjtrain

Timestamp:

04/30/06 19:16:08 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mjtrain/MJTrainSeparation.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/mjtrain/MJTrainSeparation.cc

@@ -90 +90 @@
     TH2D res2(g);
 
+    h.SetTitle("Hadronness-Distribution vs. Size");
     res1.SetTitle("Significance Li/Ma");
     res2.SetTitle("Significance Li/Ma * log_{10}(excess)");
@@ -351 +352 @@
 }
 
-Bool_t MJTrainSeparation::AutoTrain()
+Bool_t MJTrainSeparation::AutoTrain(MDataSet &set, UInt_t &seton, UInt_t &setoff)
 {
     Double_t num, min, max;
-    if (!GetEventsProduced(fDataSetTrain, num, min, max))
+    if (!GetEventsProduced(set, num, min, max))
         return kFALSE;
 
@@ -386 +387 @@
     // With an average data rate after star of
     Double_t data=0;
-    const Double_t r = GetDataRate(fDataSetTrain, data); //[Hz]
+    const Double_t r = GetDataRate(set, data); //[Hz]
 
     *fLog << "Events measured per second effective on time:             " << r << "Hz" << endl;
+    *fLog << "Total effective on time:                                  " << data/r  << endl;
 
     // this yields a number of n events to be read for training
@@ -399 +401 @@
         return kFALSE;
 
-    Double_t nummc = GetNumMC(fDataSetTrain);
+    Double_t nummc = GetNumMC(set);
 
     *fLog << "Events available in MC sample:                            " << nummc << endl;
@@ -419 +421 @@
     }
 
-    if (fNumTrainOn>0 && fNumTrainOn<on)
-    {
-        fNumTrainOff = TMath::Nint(off*fNumTrainOn/on);
-        *fLog << "Less MC events requested... scaling by " << fNumTrainOn/on << endl;
+    if (seton>0 && seton<on)
+    {
+        setoff = TMath::Nint(off*seton/on);
+        *fLog << "Less MC events requested... scaling by " << seton/on << endl;
     }
     else
     {
-        fNumTrainOn  = TMath::Nint(on);
-        fNumTrainOff = TMath::Nint(off);
-    }
-
-    *fLog << "Target number of MC events:   " << fNumTrainOn << endl;
-    *fLog << "Target number of data events: " << fNumTrainOff << endl;
-
-    fNumTestOn  = fNumTrainOn;
-    fNumTestOff = fNumTrainOff;
+        seton  = TMath::Nint(on);
+        setoff = TMath::Nint(off);
+    }
+
+    *fLog << "Target number of MC events:   " << seton  << endl;
+    *fLog << "Target number of data events: " << setoff << endl;
 
     /*
@@ -470 +469 @@
 
     if (fAutoTrain)
-        if (!AutoTrain())
+    {
+        fLog->Separator("Auto-Training -- Train-Data");
+        if (!AutoTrain(fDataSetTrain, fNumTrainOn, fNumTrainOff))
             return kFALSE;
+        fLog->Separator("Auto-Training -- Test-Data");
+        if (!AutoTrain(fDataSetTest,  fNumTestOn,  fNumTestOff))
+            return kFALSE;
+    }
 
     // --------------------- Setup files --------------------