//Root Script to create and plot the deviation Map using Runge Kutta 
//Author: Felix Boehmer

{
#include <fstream>
#include <vector>
#include <string>

 
gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
basiclibs();
gSystem->Load("libGeoBase");
gSystem->Load("libParBase");
gSystem->Load("libBase");
gSystem->Load("libMCStack");
gSystem->Load("libField");
gSystem->Load("libGen");
gSystem->Load("libPassive");
gSystem->Load("libtpc");
gSystem->Load("libgenfit");
gROOT->Reset();
gStyle->SetPalette(1);

//define variables
int rBins, zBins;
std::string efield_file, out_file;
out_file = "InhDriftTest4.dat";
efield_file = "Felix4.dat";
rBins = 3;
zBins = 4;

//std::cout<<"Enter name of the E-Field-file (usually DOLFIN output): ";
//getline(std::cin,efield_file);
//std::cout<<std::endl;
//std::ifstream infile(efield_file.c_str(), std::fstream::in);
//if (!infile.good())
//{
//  std::cerr<<"Not Existing"<<std::endl;
//  return(0);
//}
//std::cout<<"Enter the desired number of bins in r-direction: ";
//std::cin>>rBins;
//std::cout<<std::endl;
//std::cout<<"Enter the desired number of bins in z-direction: ";
//std::cin>>zBins;
//std::cout<<std::endl;
//std::cout<<"Enter name for the output file: ";
//getline(std::cin,out_file);
//std::cout<<std::endl;

TpcInhFieldDrifter* drifter = new TpcInhFieldDrifter(efield_file, out_file,
						     rBins, zBins);
std::cout<<"\nStarting the Drifter..."<<std::endl;
drifter->run();

}