mkFit.cc

Go to the documentation of this file.

#include "RecoTracker/MkFitCore/src/Matriplex/MatriplexCommon.h"

// CCCC #include "fittestMPlex.h"
#include "RecoTracker/MkFitCMS/standalone/buildtestMPlex.h"

#include "RecoTracker/MkFitCore/interface/HitStructures.h"
#include "RecoTracker/MkFitCore/interface/MkBuilder.h"
#include "RecoTracker/MkFitCore/src/MkFitter.h"
#include "RecoTracker/MkFitCMS/interface/MkStdSeqs.h"
#include "RecoTracker/MkFitCMS/standalone/MkStandaloneSeqs.h"

#include "RecoTracker/MkFitCore/interface/Config.h"
#include "RecoTracker/MkFitCore/standalone/ConfigStandalone.h"

#include "RecoTracker/MkFitCore/standalone/Event.h"

#ifdef WITH_ROOT
#include "RecoTracker/MkFitCore/standalone/Validation.h"
#include "RecoTracker/MkFitCore/standalone/RntDumper/RntDumper.h"
#endif

//#define DEBUG
#include "RecoTracker/MkFitCore/src/Debug.h"
#include "RecoTracker/MkFitCMS/standalone/Shell.h"

#include "oneapi/tbb/task_arena.h"
#include "oneapi/tbb/parallel_for.h"
#include <oneapi/tbb/global_control.h>

#if defined(USE_VTUNE_PAUSE)
#include "ittnotify.h"
#endif

#include <cstdlib>
#include <limits>
#include <list>
#include <sstream>
#include <memory>

using namespace mkfit;

//==============================================================================

void initGeom() {
  std::cout << "Constructing geometry '" << Config::geomPlugin << "'\n";

  // NB: we currently assume that each node is a layer, and that layers
  // are added starting from the center
  // NB: z is just a dummy variable, VUSolid is actually infinite in size.  *** Therefore, set it to the eta of simulation ***

  /*
  // Test CMSSW loading
  IterationsInfo ii;
  unsigned int algorithms[]={ 4,22,23,5,24,7,8,9,10,6 }; // 10 iterations
  ii.resize(10);
  for (int i = 0; i < 10; ++i) {
      ii[i].m_track_algorithm = algorithms[i];
  }
  auto xx = configJson_Load_File(ii, "mkfit-phase1-initialStep.json");
  printf("%d\n", xx->m_iteration_index);
  */

  mkfit::execTrackerInfoCreatorPlugin(Config::geomPlugin, Config::TrkInfo, Config::ItrInfo);

  ConfigJson cj(Config::json_verbose);

  if (Config::json_dump_before)
    cj.dump(Config::ItrInfo);

  if (!Config::json_load_filenames.empty()) {
    ConfigJsonPatcher::PatchReport report;

    for (auto& fn : Config::json_load_filenames) {
      // This is for testing only ... we drop the loaded IterationConfig
      // as further code will always use IterationsInfo[ iter_index ].
      cj.patchLoad_File(Config::ItrInfo, fn, &report);
    }

    printf(
        "mkFit.cc/%s--JSON-Load read %d JSON entities from %d files, replaced %d parameters.\n"
        "   NOTE that these changes were NOT APPLIED to actual configuration that is going to be used.\n",
        __func__,
        report.n_json_entities,
        report.n_files,
        report.n_replacements);
  }

  if (!Config::json_patch_filenames.empty()) {
    ConfigJsonPatcher::PatchReport report;

    cj.patch_Files(Config::ItrInfo, Config::json_patch_filenames, &report);

    printf("mkFit.cc/%s--JOSN-Patch read %d JSON entities from %d files, replaced %d parameters.\n",
           __func__,
           report.n_json_entities,
           report.n_files,
           report.n_replacements);
  }

  if (Config::json_dump_after)
    cj.dump(Config::ItrInfo);

  if (!Config::json_save_iters_fname_fmt.empty()) {
    cj.save_Iterations(
        Config::ItrInfo, Config::json_save_iters_fname_fmt, Config::json_save_iters_include_iter_info_preamble);
  }

  Config::ItrInfo.setupStandardFunctionsFromNames();

  // Test functions for ConfigJsonPatcher
  // cj.test_Direct (Config::ItrInfo[0]);
  // cj.test_Patcher(Config::ItrInfo[0]);
}

namespace {
  int g_start_event = 1;

  bool g_run_fit_std = false;

  bool g_run_build_default = true;
  bool g_run_build_cmssw = false;
  bool g_run_build_bh = false;
  bool g_run_build_std = false;
  bool g_run_build_ce = false;
  bool g_run_build_mimi = false;

  std::string g_operation = "simulate_and_process";
  ;
  std::string g_input_file = "";
  std::string g_output_file = "";

  seedOptsMap g_seed_opts;
  void init_seed_opts() {
    g_seed_opts["sim"] = {simSeeds, "Use simtracks for seeds"};
    g_seed_opts["cmssw"] = {cmsswSeeds, "Use external CMSSW seeds"};
    g_seed_opts["find"] = {findSeeds, "Use mplex seed finder for seeds"};
  }

  cleanOptsMap g_clean_opts;
  void init_clean_opts() {
    g_clean_opts["none"] = {noCleaning, "No cleaning applied to external CMSSW seeds"};
    g_clean_opts["n2"] = {cleanSeedsN2, "Apply N^2 cleaning by Mario to external CMSSW seeds"};
    g_clean_opts["pure"] = {
        cleanSeedsPure,
        "Only use external CMSSW seeds that have produced a CMSSW track \n    must enable: --read-cmssw-tracks"};
    g_clean_opts["badlabel"] = {cleanSeedsBadLabel, "Remove seeds with label()<0 in external CMSSW seeds"};
  }

  matchOptsMap g_match_opts;
  void init_match_opts() {
    g_match_opts["trkparam"] = {trkParamBased,
                                "Use track parameter-based matching for validating against CMSSW tracks"};
    g_match_opts["hits"] = {hitBased, "Use hit-based matching for validating against CMSSW tracks"};
    g_match_opts["label"] = {labelBased, "Only allowed with pure seeds: stricter hit-based matching"};
  }

  const char* b2a(bool b) { return b ? "true" : "false"; }

}  // namespace

//==============================================================================

// Getters and setters of enum configs (from command line using anon. namespace above)

template <typename T, typename U>
std::string getOpt(const T& c_opt, const U& g_opt_map) {
  static const std::string empty("");

  for (const auto& g_opt_pair : g_opt_map) {
    if (g_opt_pair.second.first == c_opt)
      return g_opt_pair.first;
  }
  std::cerr << "No match for option " << c_opt << std::endl;
  return empty;
}

template <typename T, typename U>
void setOpt(const std::string& cmd_ln_str, T& c_opt, const U& g_opt_map, const std::string& ex_txt) {
  if (g_opt_map.count(cmd_ln_str))
    c_opt = g_opt_map.at(cmd_ln_str).first;
  else {
    std::cerr << cmd_ln_str << " is not a valid " << ex_txt << " option!! Exiting..." << std::endl;
    exit(1);
  }
}

template <typename U>
void listOpts(const U& g_opt_map) {
  for (const auto& g_opt_pair : g_opt_map) {
    std::cout << "  " << g_opt_pair.first.c_str() << " : " << g_opt_pair.second.second.c_str() << std::endl;
  }
}

//==============================================================================

void test_standard() {
  printf("Running test_standard(), operation=\"%s\", best_out_of=%d\n",
         g_operation.c_str(),
         Config::finderReportBestOutOfN);

  if (Config::seedInput == cmsswSeeds)
    printf("- reading seeds from file\n");

  initGeom();

  if (Config::nEvents <= 0) {
    return;
  }

  DataFile data_file;
  if (g_operation == "read") {
    int evs_in_file = data_file.openRead(g_input_file, Config::TrkInfo.n_layers());
    int evs_available = evs_in_file - g_start_event + 1;
    if (Config::nEvents == -1) {
      Config::nEvents = evs_available;
    } else if (Config::nEvents > evs_available and not Config::loopOverFile) {
      printf("Requested number of events %d, only %d available.\n", Config::nEvents, evs_available);
      Config::nEvents = evs_available;
    }

    if (g_start_event > 1) {
      data_file.skipNEvents(g_start_event - 1);
    }
  }

  constexpr int NT = 5;
  double t_sum[NT] = {0};
  double t_skip[NT] = {0};
  std::vector<double> t_sum_iter(Config::nItersCMSSW, 0.0);
  std::vector<double> t_skip_iter(Config::nItersCMSSW, 0.0);
  double time = dtime();

  std::atomic<int> nevt{g_start_event};
  std::atomic<int> seedstot{0}, simtrackstot{0}, candstot{0};
  std::atomic<int> maxHits_all{0}, maxLayer_all{0};

  MkBuilder::populate();

  std::vector<std::unique_ptr<Event>> evs(Config::numThreadsEvents);
  std::vector<std::unique_ptr<Validation>> vals(Config::numThreadsEvents);
  std::vector<std::unique_ptr<MkBuilder>> mkbs(Config::numThreadsEvents);
  std::vector<std::shared_ptr<EventOfHits>> eohs(Config::numThreadsEvents);
  std::vector<std::shared_ptr<FILE>> fps;
  fps.reserve(Config::numThreadsEvents);

  const std::string valfile("valtree");

  for (int i = 0; i < Config::numThreadsEvents; ++i) {
    std::ostringstream serial;
    if (Config::numThreadsEvents > 1) {
      serial << "_" << i;
    }
    vals[i].reset(Validation::make_validation(valfile + serial.str() + ".root", &Config::TrkInfo));
    mkbs[i] = MkBuilder::make_builder(Config::silent);
    eohs[i].reset(new EventOfHits(Config::TrkInfo));
    evs[i].reset(new Event(*vals[i], 0, Config::TrkInfo.n_layers()));
    if (g_operation == "read") {
      fps.emplace_back(fopen(g_input_file.c_str(), "r"), [](FILE* fp) {
        if (fp)
          fclose(fp);
      });
    }
  }

  tbb::task_arena arena(Config::numThreadsFinder);

  dprint("parallel_for step size " << (Config::nEvents + Config::numThreadsEvents - 1) / Config::numThreadsEvents);

  time = dtime();

  int events_per_thread = (Config::nEvents + Config::numThreadsEvents - 1) / Config::numThreadsEvents;

  arena.execute([&]() {
    TBB_PARALLEL_FOR(
        tbb::blocked_range<int>(0, Config::numThreadsEvents, 1),
        [&](const tbb::blocked_range<int>& threads) {
          int thisthread = threads.begin();

          assert(threads.begin() == threads.end() - 1 && thisthread < Config::numThreadsEvents);

          // std::vector<Track> plex_tracks;
          auto& ev = *evs[thisthread].get();
          auto& mkb = *mkbs[thisthread].get();
          auto& eoh = *eohs[thisthread].get();
          auto fp = fps[thisthread].get();

          int evstart = thisthread * events_per_thread;
          int evend = std::min(Config::nEvents, evstart + events_per_thread);

          dprint("thisthread " << thisthread << " events " << Config::nEvents << " events/thread " << events_per_thread
                               << " range " << evstart << ":" << evend);

          for (int evt = evstart; evt < evend; ++evt) {
            ev.reset(nevt++);

            if (!Config::silent) {
              std::lock_guard<std::mutex> printlock(Event::printmutex);
              printf("\n");
              printf("Processing event %d\n", ev.evtID());
            }

            ev.read_in(data_file, fp);

            // skip events with zero seed tracks!
            if (ev.seedTracks_.empty())
              continue;

            // plex_tracks.resize(ev.simTracks_.size());

            StdSeq::loadHitsAndBeamSpot(ev, eoh);

            if (Config::useDeadModules) {
              StdSeq::loadDeads(eoh, mkfit::internal::deadvectors);
            }

            double t_best[NT] = {0}, t_cur[NT] = {0};
            std::vector<double> t_cur_iter;
            simtrackstot += ev.simTracks_.size();
            seedstot += ev.seedTracks_.size();

            int ncands_thisthread = 0;
            int maxHits_thisthread = 0;
            int maxLayer_thisthread = 0;
            for (int b = 0; b < Config::finderReportBestOutOfN; ++b) {
              t_cur[0] = 0;  // t_cur[0] = (g_run_fit_std) ? runFittingTestPlex(ev, plex_tracks) : 0;
              t_cur[1] = (g_run_build_bh) ? runBuildingTestPlexBestHit(ev, eoh, mkb) : 0;
              t_cur[2] = (g_run_build_default || g_run_build_std) ? runBuildingTestPlexStandard(ev, eoh, mkb) : 0;
              t_cur[3] = (g_run_build_default || g_run_build_ce) ? runBuildingTestPlexCloneEngine(ev, eoh, mkb) : 0;
              if (g_run_build_mimi)
                t_cur_iter = runBtpCe_MultiIter(ev, eoh, mkb, Config::nItersCMSSW);
              t_cur[4] = (g_run_build_mimi) ? t_cur_iter[Config::nItersCMSSW] : 0;
              if (g_run_build_cmssw)
                runBuildingTestPlexDumbCMSSW(ev, eoh, mkb);
              if (g_run_build_ce || g_run_build_mimi) {
                ncands_thisthread = mkb.total_cands();
                auto const& ln = mkb.max_hits_layer(eoh);
                maxHits_thisthread = ln.first;
                maxLayer_thisthread = ln.second;
              }
              for (int i = 0; i < NT; ++i)
                t_best[i] = (b == 0) ? t_cur[i] : std::min(t_cur[i], t_best[i]);

              if (!Config::silent) {
                std::lock_guard<std::mutex> printlock(Event::printmutex);
                if (Config::finderReportBestOutOfN > 1) {
                  printf("----------------------------------------------------------------\n");
                  printf("Best-of-times:");
                  for (int i = 0; i < NT; ++i)
                    printf("  %.5f/%.5f", t_cur[i], t_best[i]);
                  printf("\n");
                }
                printf("----------------------------------------------------------------\n");
              }
            }

            candstot += ncands_thisthread;
            if (maxHits_thisthread > maxHits_all) {
              maxHits_all = maxHits_thisthread;
              maxLayer_all = maxLayer_thisthread;
            }
            if (!Config::silent) {
              std::lock_guard<std::mutex> printlock(Event::printmutex);
              printf("Matriplex fit = %.5f  --- Build  BHMX = %.5f  STDMX = %.5f  CEMX = %.5f  MIMI = %.5f\n",
                     t_best[0],
                     t_best[1],
                     t_best[2],
                     t_best[3],
                     t_best[4]);
            }

            {
              static std::mutex sum_up_lock;
              std::lock_guard<std::mutex> locker(sum_up_lock);

              for (int i = 0; i < NT; ++i)
                t_sum[i] += t_best[i];
              if (evt > 0)
                for (int i = 0; i < NT; ++i)
                  t_skip[i] += t_best[i];
              if (g_run_build_mimi) {
                for (int i = 0; i < Config::nItersCMSSW; ++i)
                  t_sum_iter[i] += t_cur_iter[i];
                if (evt > 0)
                  for (int i = 0; i < Config::nItersCMSSW; ++i)
                    t_skip_iter[i] += t_cur_iter[i];
              }
            }
          }
        },
        tbb::simple_partitioner());
  });

  time = dtime() - time;

  printf("\n");
  printf("================================================================\n");
  printf("=== TOTAL for %d events\n", Config::nEvents);
  printf("================================================================\n");

  printf("Total Matriplex fit = %.5f  --- Build  BHMX = %.5f  STDMX = %.5f  CEMX = %.5f  MIMI = %.5f\n",
         t_sum[0],
         t_sum[1],
         t_sum[2],
         t_sum[3],
         t_sum[4]);
  printf("Total event > 1 fit = %.5f  --- Build  BHMX = %.5f  STDMX = %.5f  CEMX = %.5f  MIMI = %.5f\n",
         t_skip[0],
         t_skip[1],
         t_skip[2],
         t_skip[3],
         t_skip[4]);
  printf("Total event loop time %.5f simtracks %d seedtracks %d builtcands %d maxhits %d on lay %d\n",
         time,
         simtrackstot.load(),
         seedstot.load(),
         candstot.load(),
         maxHits_all.load(),
         maxLayer_all.load());
  //fflush(stdout);
  if (g_run_build_mimi) {
    printf("================================================================\n");
    for (int i = 0; i < Config::nItersCMSSW; ++i)
      std::cout << " Iteration " << i << " build time = " << t_sum_iter[i] << " \n";
    printf("================================================================\n");
    for (int i = 0; i < Config::nItersCMSSW; ++i)
      std::cout << " Iteration " << i << " build time (event > 1) = " << t_skip_iter[i] << " \n";
    printf("================================================================\n");
  }
  if (Config::quality_val) {
    printf("Sum up of quality-val:\n");
    StdSeq::Quality::s_quality_sum.quality_print();
  }
  if (g_operation == "read") {
    data_file.close();
  }

  for (auto& val : vals) {
    val->fillConfigTree();
    val->saveTTrees();
  }
}

//==============================================================================
// Command line argument parsing
//==============================================================================

typedef std::list<std::string> lStr_t;
typedef lStr_t::iterator lStr_i;

bool has_suffix(const std::string& str, const std::string& suffix) {
  return str.size() >= suffix.size() && str.compare(str.size() - suffix.size(), suffix.size(), suffix) == 0;
}

void next_arg_or_die(lStr_t& args, lStr_i& i, bool allow_single_minus = false) {
  lStr_i j = i;
  if (++j == args.end() || has_suffix(*j, ".C") || ((*j)[0] == '-' && !(*j == "-" && allow_single_minus))) {
    std::cerr << "Error: option " << *i << " requires an argument.\n";
    exit(1);
  }
  i = j;
}

//==============================================================================
// main
//==============================================================================

#include <fenv.h>

int main(int argc, const char* argv[]) {
#ifdef _GNU_SOURCE
  if (Const::nan_etc_sigs_enable) {
    feenableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW);  //FE_ALL_EXCEPT);
  }
#endif

#ifdef USE_VTUNE_PAUSE
  __itt_pause();
#endif

  assert(sizeof(Track::Status) == 4 &&
         "To make sure this is true for icc and gcc<6 when mixing bools/ints in bitfields.");

  // init enum maps
  init_seed_opts();
  init_clean_opts();
  init_match_opts();

  lStr_t mArgs;
  for (int i = 1; i < argc; ++i) {
    mArgs.push_back(argv[i]);
  }
  bool run_shell = false;

  lStr_i i = mArgs.begin();
  while (i != mArgs.end()) {
    lStr_i start = i;

    if (*i == "-h" || *i == "-help" || *i == "--help") {
      printf(
          "\n"
          "Usage: %s [options]\n"
          "Options: defaults defined as (def: DEFAULT VALUE)\n"
          "\n----------------------------------------------------------------------------------------------------------"
          "\n\n"
          "Generic options\n\n"
          "  --geom           <str>   geometry plugin to use (def: %s)\n"
          "  --silent                 suppress printouts inside event loop (def: %s)\n"
          "  --best-out-of    <int>   run test num times, report best time (def: %d)\n"
          "  --input-file             file name for reading (def: %s)\n"
          "  --output-file            file name for writitng (def: %s)\n"
          "  --read-cmssw-tracks      read external cmssw reco tracks if available (def: %s)\n"
          "  --read-simtrack-states   read in simTrackStates for pulls in validation (def: %s)\n"
          "  --num-events     <int>   number of events to run over or simulate (def: %d)\n"
          "                           if using --input-file, must be enabled AFTER on command line\n"
          "  --start-event    <int>   event number to start at when reading from a file (def: %d)\n"
          "  --loop-over-file         after reaching the end of the file, start over from the beginning until "
          "                           <num-events> events have been processed\n"
          "  --shell                  start interactive shell instead of running test_standard()\n"
          "\n"
          "If no --input-file is specified, will trigger simulation\n"
          "  --num-tracks     <int>   number of tracks to generate for each event (def: %d)\n"
          "\n----------------------------------------------------------------------------------------------------------"
          "\n\n"
          "Threading related options\n\n"
          "  --num-thr-sim    <int>   number of threads for simulation (def: %d)\n"
          "  --num-thr        <int>   number of threads for track finding (def: %d)\n"
          "  --num-thr-ev     <int>   number of threads to run the event loop (def: %d)\n"
          "  --seeds-per-task <int>   number of seeds to process in a tbb task (def: %d)\n"
          "  --hits-per-task  <int>   number of layer1 hits per task when using find seeds (def: %d)\n"
          "\n----------------------------------------------------------------------------------------------------------"
          "\n\n"
          "FittingTestMPlex options\n\n"
          "  --fit-std                run standard fitting test (def: %s)\n"
          "  --cf-fitting             enable conformal fit before fitting tracks to get initial estimate of track "
          "parameters and errors (def: %s)\n"
          "  --fit-val                enable ROOT based validation for fittingMPlex  (def: %s)\n"
          "\n----------------------------------------------------------------------------------------------------------"
          "\n\n"
          "BuildingTestMPlex options\n\n"
          " **Specify which building routine you would like to run\n"
          "  --build-cmssw            run dummy validation of CMSSW tracks with MkBuilder stuff (def: %s)\n"
          "  --build-bh               run best-hit building test (def: %s)\n"
          "  --build-std              run standard combinatorial building test (def: %s)\n"
          "  --build-ce               run clone engine combinatorial building test (def: %s)\n"
          "  --build-mimi             run clone engine on multiple-iteration test (def: %s)\n"
          "  --num-iters-cmssw <int>  number of mimi iterations to run (def: set to 3 when --build-mimi is in effect, "
          "0 otherwise)\n"
          "\n"
          " **Seeding options\n"
          "  --seed-input     <str>   which seed collecion used for building (def: %s)\n"
          "  --seed-cleaning  <str>   which seed cleaning to apply if using cmssw seeds (def: %s)\n"
          "  --cf-seeding             enable conformal fit over seeds (def: %s)\n"
          "\n"
          " **Duplicate removal options\n"
          "  --remove-dup            run duplicate removal after building, using both hit and kinematic criteria (def: "
          "%s)\n"
          "  --remove-dup-no-hit     run duplicate removal after building, using kinematic criteria only (def: %s)\n"
          "\n"
          " **Dead module (strip) option\n"
          "  --use-dead-modules          run duplicate removal after building, using both hit and kinematic criteria "
          "(def: %s)\n"
          "\n"
          " **Additional options for building\n"
          "  --use-phiq-arr           use phi-Q arrays in select hit indices (def: %s)\n"
          "  --kludge-cms-hit-errors  make sure err(xy) > 15 mum, err(z) > 30 mum (def: %s)\n"
          "  --backward-fit           perform backward fit during building (def: %s)\n"
          "  --no-backward-search     do not do backward search after backward fit\n"
          "                           (def: do search if backward-fit is enabled and available in given iteration)\n"
          "  --include-pca            do the backward fit to point of closest approach, does not imply "
          "'--backward-fit' (def: %s)\n"
          "  --use-p2p <0|1>          use prop-to-plane (def: %d)\n"
          "  --use-ptms <0|1>         use pT multiple scattering (def: %d)\n"
          "\n----------------------------------------------------------------------------------------------------------"
          "\n\n"
          "Validation options\n\n"
          " **Text file based options\n"
          "  --quality-val            enable printout validation for MkBuilder (def: %s)\n"
          "                             must enable: --dump-for-plots\n"
          "  --dump-for-plots         make shell printouts for plots (def: %s)\n"
          "  --mtv-like-val           configure validation to emulate CMSSW MultiTrackValidator (MTV) (def: %s)\n"
          "  --mtv-require-seeds      configure validation to emulate MTV but require sim tracks to be matched to "
          "seeds (def: %s)\n"
          "\n"
          " **ROOT based options\n"
          "  --sim-val-for-cmssw      enable ROOT based validation for CMSSW tracks with simtracks as reference [eff, "
          "FR, DR] (def: %s)\n"
          "  --sim-val                enable ROOT based validation for seeding, building, and fitting with simtracks "
          "as reference [eff, FR, DR] (def: %s)\n"
          "  --cmssw-val              enable ROOT based validation for building and fitting with CMSSW tracks as "
          "reference [eff, FR, DR] (def: %s)\n"
          "                             must enable: --geom CMS-phase1 --read-cmssw-tracks\n"
          "  --cmssw-match-fw  <str>  which cmssw track matching routine to use if validating against CMSSW tracks, "
          "forward built tracks only (def: %s)\n"
          "                             must enable: --geom CMS-phase1 --cmssw-val --read-cmssw-tracks\n"
          "  --cmssw-match-bk  <str>  which cmssw track matching routine to use if validating against CMSSW tracks, "
          "backward fit tracks only (def: %s)\n"
          "                             must enable: --geom CMS-phase1 --cmssw-val --read-cmssw-tracks --backward-fit "
          "--backward-fit-pca\n"
          "  --inc-shorts             include short reco tracks into FR (def: %s)\n"
          "  --keep-hit-info          keep vectors of hit idxs and branches in trees (def: %s)\n"
          "  --try-to-save-sim-info   two options for this flag [related to validation with simtracks as reference "
          "collection] (def: %s)\n"
          "                              a) IF (--read-simtrack-states is enabled)\n"
          "                                    If a sim track is associated to a reco track, but it does not contain "
          "the last found on the reco track\n"
          "                                    still save sim track kinematic info from generator position\n"
          "                              b) ELSE (--read-simtrack-states is NOT enabled)\n"
          "                                    Save sim track kinematic info from generator position if matched to "
          "reco track\n"
          "\n----------------------------------------------------------------------------------------------------------"
          "\n\n"
          "Combo spaghetti, that's with cole slaw:\n\n"
          " **Building and fitting combo options\n"
          "  --backward-fit-pca       perform backward fit to point of closest approach during building\n"
          "                             == --backward-fit --include-pca\n"
          " **Seed combo options\n"
          "  --cmssw-simseeds         use CMS geom with simtracks for seeds\n"
          "                             == --geom CMS-phase1 --seed-input %s\n"
          "  --cmssw-stdseeds         use CMS geom with CMSSW seeds uncleaned\n"
          "                             == --geom CMS-phase1 --seed-input %s --seed-cleaning %s\n"
          "  --cmssw-n2seeds          use CMS geom with CMSSW seeds cleaned with N^2 routine\n"
          "                             == --geom CMS-phase1 --seed-input %s --seed-cleaning %s\n"
          "  --cmssw-pureseeds        use CMS geom with pure CMSSW seeds (seeds which produced CMSSW reco tracks), "
          "enable read of CMSSW tracks\n"
          "                             == --geom CMS-phase1 --seed-input %s --seed-cleaning %s --read-cmssw-tracks\n"
          "  --cmssw-goodlabelseeds   use CMS geom with CMSSW seeds with label() >= 0\n"
          "                             == --geom CMS-phase1 --seed-input %s --seed-cleaning %s\n"
          "\n"
          " **CMSSW validation combo options\n"
          "  --cmssw-val-fhit-bhit    use CMSSW validation with hit based matching (50 percent after seed) for forward "
          "built tracks\n"
          "                           use CMSSW validation with hit based matching (50 percent after seed) for "
          "backward fit tracks\n"
          "                             == --cmssw-val --read-cmssw-tracks --cmssw-match-fw %s --cmssw-match-bk %s\n"
          "                             must enable: --backward-fit-pca\n"
          "  --cmssw-val-fhit-bprm    use CMSSW validation with hit based matching (50 percent after seed) for forward "
          "built tracks\n"
          "                           use CMSSW validation with track parameter based matching for backward fit "
          "tracks\n"
          "                             == --cmssw-val --read-cmssw-tracks --cmssw-match-fw %s --cmssw-match-bk %s\n"
          "                             must enable: --backward-fit-pca\n"
          "  --cmssw-val-fprm-bhit    use CMSSW validation with track parameter based matching for forward built "
          "tracks\n"
          "                           use CMSSW validation with hit based matching (50 percent after seed) for "
          "backward fit tracks\n"
          "                             == --cmssw-val --read-cmssw-tracks --cmssw-match-fw %s --cmssw-match-bk %s\n"
          "                             must enable: --backward-fit-pca\n"
          "  --cmssw-val-fprm-bprm    use CMSSW validation with track parameter based matching for forward built "
          "tracks\n"
          "                           use CMSSW validation with track parameter based matching for backward fit "
          "tracks\n"
          "                             == --cmssw-val --read-cmssw-tracks --cmssw-match-fw %s --cmssw-match-bk %s\n"
          "                             must enable: --backward-fit-pca\n"
          "  --cmssw-val-label        use CMSSW validation with stricter hit based matching for both forward built and "
          "backward fit tracks, enable read of CMSSW tracks\n"
          "                             == --cmssw-val --read-cmssw-tracks --cmssw-match-fw %s --cmssw-match-bk %s\n"
          "                             must enable: --cmssw-pureseeds --backward-fit-pca\n"
          "\n----------------------------------------------------------------------------------------------------------"
          "\n\n"
          "JSON config patcher options:\n\n"
          "  --json-load  <filename>  load single IterationConfig from given JSON file (def: do not load)\n"
          "                           can be specified multiple times for several files\n"
          "  --json-patch <filename>  patch IterationsInfo from given JSON file (def: do not patch)\n"
          "                           can be specified multiple times for several files\n"
          "  --json-save-iterations <fname-fmt> save per iteration json files\n"
          "                           %%d in fname-fmt gets replaced with iteration index\n"
          "                           %%s in fname-fmt gets replaced with iteration algorithm name\n"
          "                           exactly one of %%d and %%s must be specified\n"
          "  --json-save-iterations-include-iter-info-preamble (def: %s)\n"
          "  --json-verbose     print each patch assignment as it is being made (def: %s)\n"
          "  --json-dump-before print iteration config before patching (def: %s)\n"
          "  --json-dump-after  print iteration config after  patching (def: %s)\n"
          "\n----------------------------------------------------------------------------------------------------------"
          "\n\n",
          argv[0],

          Config::geomPlugin.c_str(),
          b2a(Config::silent),
          Config::finderReportBestOutOfN,
          g_input_file.c_str(),
          g_output_file.c_str(),
          b2a(Config::readCmsswTracks),
          b2a(Config::readSimTrackStates),
          Config::nEvents,
          g_start_event,
          Config::nTracks,

          Config::numThreadsSimulation,
          Config::numThreadsFinder,
          Config::numThreadsEvents,
          Config::numSeedsPerTask,
          Config::numHitsPerTask,

          b2a(g_run_fit_std),
          b2a(Config::cf_fitting),
          b2a(Config::fit_val),

          b2a(g_run_build_cmssw),
          b2a(g_run_build_bh),
          b2a(g_run_build_default || g_run_build_std),
          b2a(g_run_build_default || g_run_build_ce),
          b2a(g_run_build_mimi),

          getOpt(Config::seedInput, g_seed_opts).c_str(),
          getOpt(Config::seedCleaning, g_clean_opts).c_str(),
          b2a(Config::cf_seeding),

          b2a(Config::removeDuplicates && Config::useHitsForDuplicates),
          b2a(Config::removeDuplicates && !Config::useHitsForDuplicates),

          b2a(Config::useDeadModules),

          b2a(Config::usePhiQArrays),
          b2a(Config::kludgeCmsHitErrors),
          b2a(Config::backwardFit),
          b2a(Config::includePCA),
          int(Config::usePropToPlane),
          int(Config::usePtMultScat),

          b2a(Config::quality_val),
          b2a(Config::dumpForPlots),
          b2a(Config::mtvLikeValidation),
          b2a(Config::mtvRequireSeeds),

          b2a(Config::sim_val_for_cmssw),
          b2a(Config::sim_val),
          b2a(Config::cmssw_val),
          getOpt(Config::cmsswMatchingFW, g_match_opts).c_str(),
          getOpt(Config::cmsswMatchingBK, g_match_opts).c_str(),
          b2a(Config::inclusiveShorts),
          b2a(Config::keepHitInfo),
          b2a(Config::tryToSaveSimInfo),

          getOpt(simSeeds, g_seed_opts).c_str(),
          getOpt(cmsswSeeds, g_seed_opts).c_str(),
          getOpt(noCleaning, g_clean_opts).c_str(),
          getOpt(cmsswSeeds, g_seed_opts).c_str(),
          getOpt(cleanSeedsN2, g_clean_opts).c_str(),
          getOpt(cmsswSeeds, g_seed_opts).c_str(),
          getOpt(cleanSeedsPure, g_clean_opts).c_str(),
          getOpt(cmsswSeeds, g_seed_opts).c_str(),
          getOpt(cleanSeedsBadLabel, g_clean_opts).c_str(),

          getOpt(hitBased, g_match_opts).c_str(),
          getOpt(hitBased, g_match_opts).c_str(),
          getOpt(hitBased, g_match_opts).c_str(),
          getOpt(trkParamBased, g_match_opts).c_str(),
          getOpt(trkParamBased, g_match_opts).c_str(),
          getOpt(hitBased, g_match_opts).c_str(),
          getOpt(trkParamBased, g_match_opts).c_str(),
          getOpt(trkParamBased, g_match_opts).c_str(),
          getOpt(labelBased, g_match_opts).c_str(),
          getOpt(labelBased, g_match_opts).c_str(),

          b2a(Config::json_verbose),
          b2a(Config::json_save_iters_include_iter_info_preamble),
          b2a(Config::json_dump_before),
          b2a(Config::json_dump_after));

      printf("List of options for string based inputs \n");
      printf(
          "--geom \n"
          "  CMS-phase1 \n"
          "  CMS-phase2 \n"
          "  CylCowWLids \n"
          "\n");

      printf("--seed-input \n");
      listOpts(g_seed_opts);
      printf("\n");

      printf("--seed-cleaning \n");
      listOpts(g_clean_opts);
      printf("\n");

      printf("--cmssw-matching \n");
      listOpts(g_match_opts);
      printf("\n");

      exit(0);
    }  // end of "help" block

    else if (*i == "--geom") {
      next_arg_or_die(mArgs, i);
      Config::geomPlugin = *i;
    } else if (*i == "--silent") {
      Config::silent = true;
    } else if (*i == "--best-out-of") {
      next_arg_or_die(mArgs, i);
      Config::finderReportBestOutOfN = atoi(i->c_str());
    } else if (*i == "--input-file") {
      next_arg_or_die(mArgs, i);
      g_input_file = *i;
      g_operation = "read";
      Config::nEvents = -1;
    } else if (*i == "--output-file") {
      next_arg_or_die(mArgs, i);
      g_output_file = *i;
      g_operation = "write";
    } else if (*i == "--read-cmssw-tracks") {
      Config::readCmsswTracks = true;
    } else if (*i == "--read-simtrack-states") {
      Config::readSimTrackStates = true;
    } else if (*i == "--num-events") {
      next_arg_or_die(mArgs, i);
      Config::nEvents = atoi(i->c_str());
    } else if (*i == "--start-event") {
      next_arg_or_die(mArgs, i);
      g_start_event = atoi(i->c_str());
    } else if (*i == "--loop-over-file") {
      Config::loopOverFile = true;
    } else if (*i == "--shell") {
#ifndef WITH_ROOT
      std::cerr << "--shell option is only supported when compiled with ROOT.\n";
      exit(1);
#endif
      run_shell = true;
    } else if (*i == "--num-tracks") {
      next_arg_or_die(mArgs, i);
      Config::nTracks = atoi(i->c_str());
    } else if (*i == "--num-thr-sim") {
      next_arg_or_die(mArgs, i);
      Config::numThreadsSimulation = atoi(i->c_str());
    } else if (*i == "--num-thr") {
      next_arg_or_die(mArgs, i);
      Config::numThreadsFinder = atoi(i->c_str());
    } else if (*i == "--num-thr-ev") {
      next_arg_or_die(mArgs, i);
      Config::numThreadsEvents = atoi(i->c_str());
    } else if (*i == "--seeds-per-task") {
      next_arg_or_die(mArgs, i);
      Config::numSeedsPerTask = atoi(i->c_str());
    } else if (*i == "--hits-per-task") {
      next_arg_or_die(mArgs, i);
      Config::numHitsPerTask = atoi(i->c_str());
    } else if (*i == "--fit-std") {
      g_run_build_default = false;
      g_run_fit_std = true;
    } else if (*i == "--cf-fitting") {
      Config::cf_fitting = true;
    } else if (*i == "--fit-val") {
      Config::fit_val = true;
    } else if (*i == "--build-cmssw") {
      g_run_build_default = false;
      g_run_build_cmssw = true;
    } else if (*i == "--build-bh") {
      g_run_build_default = false;
      g_run_build_bh = true;
    } else if (*i == "--build-std") {
      g_run_build_default = false;
      g_run_build_std = true;
    } else if (*i == "--build-ce") {
      g_run_build_default = false;
      g_run_build_ce = true;
    } else if (*i == "--build-mimi") {
      g_run_build_default = false;
      g_run_build_mimi = true;
      if (Config::nItersCMSSW == 0)
        Config::nItersCMSSW = 3;
    } else if (*i == "--num-iters-cmssw") {
      next_arg_or_die(mArgs, i);
      Config::nItersCMSSW = atoi(i->c_str());
    } else if (*i == "--seed-input") {
      next_arg_or_die(mArgs, i);
      setOpt(*i, Config::seedInput, g_seed_opts, "seed input collection");
    } else if (*i == "--seed-cleaning") {
      next_arg_or_die(mArgs, i);
      setOpt(*i, Config::seedCleaning, g_clean_opts, "seed cleaning");
    } else if (*i == "--cf-seeding") {
      Config::cf_seeding = true;
    } else if (*i == "--use-phiq-arr") {
#ifdef CONFIG_PhiQArrays
      Config::usePhiQArrays = true;
#else
      printf("--use-phiq-arr has no effect: recompile with CONFIG_PhiQArrays\n");
#endif
    } else if (*i == "--remove-dup") {
      Config::removeDuplicates = true;
      Config::useHitsForDuplicates = true;
    } else if (*i == "--remove-dup-no-hit") {
      Config::removeDuplicates = true;
      Config::useHitsForDuplicates = false;
    } else if (*i == "--use-dead-modules") {
      Config::useDeadModules = true;
    } else if (*i == "--kludge-cms-hit-errors") {
      Config::kludgeCmsHitErrors = true;
    } else if (*i == "--backward-fit") {
      Config::backwardFit = true;
    } else if (*i == "--no-backward-search") {
      Config::backwardSearch = false;
    } else if (*i == "--include-pca") {
      Config::includePCA = true;
    } else if (*i == "--use-p2p") {
      next_arg_or_die(mArgs, i);
      Config::usePropToPlane = (bool)atoi(i->c_str());
    } else if (*i == "--use-ptms") {
      next_arg_or_die(mArgs, i);
      Config::usePtMultScat = (bool)atoi(i->c_str());
    } else if (*i == "--quality-val") {
      Config::quality_val = true;
    } else if (*i == "--dump-for-plots") {
      Config::dumpForPlots = true;
    } else if (*i == "--mtv-like-val") {
      Config::mtvLikeValidation = true;
      Config::cmsSelMinLayers = 0;
      Config::nMinFoundHits = 0;
    } else if (*i == "--mtv-require-seeds") {
      Config::mtvLikeValidation = true;
      Config::cmsSelMinLayers = 0;
      Config::nMinFoundHits = 0;
      Config::mtvRequireSeeds = true;
    } else if (*i == "--sim-val-for-cmssw") {
      Config::sim_val_for_cmssw = true;
    } else if (*i == "--sim-val") {
      Config::sim_val = true;
    } else if (*i == "--cmssw-val") {
      Config::cmssw_val = true;
    } else if (*i == "--cmssw-match-fw") {
      next_arg_or_die(mArgs, i);
      setOpt(*i, Config::cmsswMatchingFW, g_match_opts, "CMSSW validation track matching for forward built tracks");
    } else if (*i == "--cmssw-match-bk") {
      next_arg_or_die(mArgs, i);
      setOpt(*i, Config::cmsswMatchingBK, g_match_opts, "CMSSW validation track matching for backward fit tracks");
    } else if (*i == "--inc-shorts") {
      Config::inclusiveShorts = true;
    } else if (*i == "--keep-hit-info") {
      Config::keepHitInfo = true;
    } else if (*i == "--try-to-save-sim-info") {
      Config::tryToSaveSimInfo = true;
    } else if (*i == "--backward-fit-pca") {
      Config::backwardFit = true;
      Config::includePCA = true;
    } else if (*i == "--cmssw-simseeds") {
      Config::geomPlugin = "CMS-phase1";
      Config::seedInput = simSeeds;
    } else if (*i == "--cmssw-stdseeds") {
      Config::geomPlugin = "CMS-phase1";
      Config::seedInput = cmsswSeeds;
      Config::seedCleaning = noCleaning;
    } else if (*i == "--cmssw-n2seeds") {
      Config::geomPlugin = "CMS-phase1";
      Config::seedInput = cmsswSeeds;
      Config::seedCleaning = cleanSeedsN2;
    } else if (*i == "--cmssw-pureseeds") {
      Config::geomPlugin = "CMS-phase1";
      Config::seedInput = cmsswSeeds;
      Config::seedCleaning = cleanSeedsPure;
      Config::readCmsswTracks = true;
    } else if (*i == "--cmssw-goodlabelseeds") {
      Config::geomPlugin = "CMS-phase1";
      Config::seedInput = cmsswSeeds;
      Config::seedCleaning = cleanSeedsBadLabel;
    } else if (*i == "--cmssw-val-fhit-bhit") {
      Config::cmssw_val = true;
      Config::readCmsswTracks = true;
      Config::cmsswMatchingFW = hitBased;
      Config::cmsswMatchingBK = hitBased;
    } else if (*i == "--cmssw-val-fhit-bprm") {
      Config::cmssw_val = true;
      Config::readCmsswTracks = true;
      Config::cmsswMatchingFW = hitBased;
      Config::cmsswMatchingBK = trkParamBased;
    } else if (*i == "--cmssw-val-fprm-bhit") {
      Config::cmssw_val = true;
      Config::readCmsswTracks = true;
      Config::cmsswMatchingFW = trkParamBased;
      Config::cmsswMatchingBK = hitBased;
    } else if (*i == "--cmssw-val-fprm-bprm") {
      Config::cmssw_val = true;
      Config::readCmsswTracks = true;
      Config::cmsswMatchingFW = trkParamBased;
      Config::cmsswMatchingBK = trkParamBased;
    } else if (*i == "--cmssw-val-label") {
      Config::cmssw_val = true;
      Config::readCmsswTracks = true;
      Config::cmsswMatchingFW = labelBased;
      Config::cmsswMatchingBK = labelBased;
    } else if (*i == "--json-load") {
      next_arg_or_die(mArgs, i);
      Config::json_load_filenames.push_back(*i);
    } else if (*i == "--json-patch") {
      next_arg_or_die(mArgs, i);
      Config::json_patch_filenames.push_back(*i);
    } else if (*i == "--json-save-iterations") {
      next_arg_or_die(mArgs, i);
      Config::json_save_iters_fname_fmt = *i;
    } else if (*i == "--json-save-iterations-include-iter-info-preamble") {
      Config::json_save_iters_include_iter_info_preamble = true;
    } else if (*i == "--json-verbose") {
      Config::json_verbose = true;
    } else if (*i == "--json-dump-before") {
      Config::json_dump_before = true;
    } else if (*i == "--json-dump-after") {
      Config::json_dump_after = true;
    } else {
      fprintf(stderr, "Error: Unknown option/argument '%s'.\n", i->c_str());
      exit(1);
    }

    mArgs.erase(start, ++i);
  }

  // clang-format off

  // Do some checking of options before going...
  if (Config::seedCleaning != cleanSeedsPure &&
      (Config::cmsswMatchingFW == labelBased || Config::cmsswMatchingBK == labelBased)) {
    std::cerr << "What have you done?!? Can't mix cmssw label matching without pure seeds! Exiting...\n";
    exit(1);
  } else if (Config::mtvLikeValidation && Config::inclusiveShorts) {
    std::cerr << "What have you done?!? Short reco tracks are already accounted for in the MTV-Like Validation! Inclusive "
                 "shorts is only an option for the standard simval, and will break the MTV-Like simval! Exiting...\n";
    exit(1);
  }

  Config::recalculateDependentConstants();

  printf("mkFit configuration complete.\n"
         "  vusize=%d, num_thr_events=%d, num_thr_finder=%d\n"
         "  sizeof(Track)=%zu, sizeof(Hit)=%zu, sizeof(SVector3)=%zu, sizeof(SMatrixSym33)=%zu, sizeof(MCHitInfo)=%zu\n",
         MPT_SIZE, Config::numThreadsEvents, Config::numThreadsFinder,
         sizeof(Track), sizeof(Hit), sizeof(SVector3), sizeof(SMatrixSym33), sizeof(MCHitInfo));

#ifdef WITH_ROOT
  Shell *shell = nullptr;
#endif
  if (run_shell) {
#ifdef WITH_ROOT
    tbb::global_control global_limit(tbb::global_control::max_allowed_parallelism, Config::numThreadsFinder);

    initGeom();
    shell = new Shell(mkfit::internal::deadvectors, g_input_file, g_start_event);
    shell->Run();
#else
    std::cerr << "shell selected on a non-ROOT build.\n";
#endif
  } else {
    test_standard();
  }

#ifdef WITH_ROOT
  RntDumper::FinalizeAll();
  if (run_shell) {
    delete shell;
  }
#endif

  // clang-format on

  return 0;
}