#include <algorithm>
#include <omp.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <vector>
#include <iostream>

#include "util.h"
#include "riemann-track-finder.h"

using namespace std;

int main(int argc, char **argv) {

	if(argc == 7){
		RiemannTrackFinderOptions opts(argv[1], argv[2]);
		opts.max_sz_chi_square = atof(argv[3]);
		opts.max_dist_from_plane = atof(argv[4]);
		opts.max_dist_from_sz_line = atof(argv[5]);
		opts.verbosity = atof(argv[6]);

		RiemannTrackFinder finder(opts);
		finder.find_tracks();
	}
	else{
		RiemannTrackFinderOptions opts("data/hit_points_with_layer_information_3.txt", "output.txt");
		//RiemannTrackFinderOptions opts("data/hit_points_with_layer_information.txt", "output.txt");
		//	RiemannTrackFinderOptions opts("data/mvd-hits.csv", "output.txt");

		//opts.verbosity = 90;
		opts.verbosity = 9999;
		RiemannTrackFinder finder(opts);
		finder.find_tracks();
	}

	return 0;
}  // main