Changeset 19802 for trunk/FACT++/src/rootifysql.cc

Timestamp:

10/27/19 11:11:01 (5 years ago)

Author:

tbretz

Message:

Moved the common Splitting algorithm to a common class 'Splitting' (added) and added the possibility to rename columns (csv2root)

File:

• trunk/FACT++/src/rootifysql.cc (modified) (11 diffs)

trunk/FACT++/src/rootifysql.cc

@@ -1 +1 @@
 #include "Database.h"
-
-#include <random>
 
 #include <boost/regex.hpp>
@@ -10 +8 @@
 #include "tools.h"
 #include "Time.h"
-#include "Configuration.h"
+#include "Splitting.h"
 
 #include <TROOT.h>
@@ -62 +60 @@
         ;
 
-    po::options_description split("Splitting options");
-    split.add_options()
-        ("split-sequence,S", vars<uint16_t>(),            "Split data sequentially into several trees/files (e.g. 1, 1, 2)")
-        ("split-quantile,Q", vars<double>(),              "Split data randomly into several trees/files (e.g. 0.5, 1)")
-        ("seed", var<uint64_t>(mt19937_64::default_seed), "Seed value in case of random split")
-        ;
-
     po::positional_options_description p;
     p.add("file", 1); // The 1st positional options (n=1)
@@ -76 +67 @@
     conf.AddOptions(ascii);
     conf.AddOptions(root);
-    conf.AddOptions(split);
+    conf.AddOptions(Tools::Splitting::options());
     conf.SetArgumentPositions(p);
 }
@@ -125 +116 @@
         "/*comment*/ or introduced with # (shell script style) or -- (SQL style).\n"
         "\n"
-        "For several purposes, it might be convenient to split the output to several "
-        "different root-trees or ascii files. This can be done using the --split-sequence (-S) "
-        "and the --split-quantile (-Q) options. If a split sequence is defined as "
-        "-S 1 -S 2 -S 1 the events are split by 1:2:1 in this sequence order. If "
-        "quantiled are given as -Q 0.5 -Q 0.6, the first tree will contain 50% of "
-        "the second one 10% and the third one 40%. The corresponding seed value can "
-        "be set with --seed.\n"
+        << Tools::Splitting::usage() <<
         "\n"
         "In case of success, 0 is returned, a value>0 otherwise.\n"
@@ -522 +507 @@
     // ----------------------- Setup splitting ---------------------------------
 
-    vector<uint16_t> split_seq   = conf.Vec<uint16_t>("split-sequence");
-    vector<double>   split_quant = conf.Vec<double>("split-quantile");
-
-    if (!split_seq.empty() && !split_quant.empty())
-        throw runtime_error("Only splitting by --split-sequence or --split-quantile is allowed.");
-
-    const size_t num_split = split_seq.size()+split_quant.size()==0 ? 0 :
-        ::max(split_seq.size(), split_quant.size()+1);
-
-    map<size_t, size_t> split_lut;
-    for (size_t i=0; i<split_seq.size(); i++)
-    {
-        const size_t sz = split_lut.size();
-        for (size_t j=0; j<split_seq[i]; j++)
-            split_lut.emplace(j+sz, i);
-    }
-
-    for (size_t i=0; i<split_quant.size(); i++)
-        if (split_quant[i]<0 || split_quant[i]>=1)
-            throw runtime_error("Splitting quantiles must be in the range [0;1)");
-
-    for (size_t i=1; i<split_quant.size(); i++)
-    {
-        if (split_quant[i]<=split_quant[i-1])
-            throw runtime_error("Splitting quantiles must be in increasing order.");
-    }
-
-    // -------------------------------------------------------------------------
-
     const auto vars = conf.GetWildcardOptions("var.*");
 
@@ -558 +514 @@
 
     // -------------------------------------------------------------------------
+
+    /*const*/ Tools::Splitting split(conf);
 
     if (verbose>0)
@@ -826 +784 @@
         cout << "Opening file '" << path << "' [compression=" << compression << "]...\n";
         cout << "Writing data to tree '" << tree << "'" << (nofill?" (--skipped--)":"") << endl;
-        if (num_split)
-        {
-            cout << "Splitting configured " << (split_seq.empty()?"randomly":"in sequence") << " into " << num_split << " branches." << endl;
-            if (!split_quant.empty())
-                cout << "Seed value configured as " << conf.Get<uint64_t>("seed") << "." << endl;
-        }
+        split.print();
     }
 
@@ -875 +828 @@
     vector<TTree*> ttree;
 
-    if (num_split==0)
+    if (split.empty())
         ttree.emplace_back(new TTree(tree.c_str(), query.c_str()));
     else
-        for (size_t i=0; i<num_split; i++)
+        for (size_t i=0; i<split.size(); i++)
             ttree.emplace_back(new TTree((tree+"["+to_string(i)+"]").c_str(), query.c_str()));
 
@@ -955 +908 @@
     {
         vector<string> names;
-        if (num_split==0)
+        if (split.empty())
             names.emplace_back(write);
         else
-            for (size_t i=0; i<num_split; i++)
+            for (size_t i=0; i<split.size(); i++)
                 names.emplace_back(write+"-"+to_string(i));
 
@@ -1017 +970 @@
     // ---------------------- Fill TTree with DB data --------------------------
 
-    const uniform_real_distribution<double> distribution(0,1);
-    mt19937_64 generator;
-    generator.seed(conf.Get<uint64_t>("seed"));
-    auto rndm = bind(distribution, generator);
-
     size_t count = 0;
     size_t skip  = 0;
     do
     {
-        size_t index = 0;
-        if (!split_lut.empty())
-            index = split_lut[count % split_lut.size()];
-        if (!split_quant.empty())
-        {
-            const float r = rndm();
-            for (; r>=split_quant[index]; index++)
-                if (index==split_quant.size())
-                    break;
-        }
-
-        count++;
+        size_t index = split.index(count++);
 
         ostringstream rtxt;