#include <cassert>
#include <cmath>
#include <vector>

#include "TH1.h"
#include "TTree.h"
#include "TNtuple.h"
#include "TH1D.h"
#include "TCanvas.h"
#include "TRandom.h"

#include <TGraph.h>
#include <TF1.h>
#include <TMath.h>
#include <TStyle.h>

#include <iostream>
#include <fstream>
#include <string>


void ANTS()
{  

TCanvas* c = new TCanvas();        


//------------------------------------------ Read -------------------------------------------------- 
auto g1 = ROOT::RDF::MakeCsvDataFrame("8000_values.vacuum.txt", false, ' ').Graph("Col0", "Col4");
//auto g2 = ROOT::RDF::MakeCsvDataFrame("Ref_values.vacuum", false, ' ').Graph("Col1", "Col4");




g1->SetLineColor(kAzure);
//g2->SetLineColor(kRed);



TMultiGraph mg;
mg.Add(g1.GetPtr());
//mg.Add(g2.GetPtr());

//c->SetLogx();

mg.DrawClone("AL");

mg.SetTitle("Development of Iota");
//mg.GetYaxis()->SetRangeUser(Min-margin,);
mg.GetXaxis()->SetRangeUser(5.935982,6.21);
//mg.SetAxisRange(5.935982, 6.21,"X");

/*mg.GetXaxis()->SetTitle("Major radius");
mg.GetXaxis()->CenterTitle();
mg.GetXaxis()->SetTitleSize(0.05);
mg.GetXaxis()->SetTickLength(-0.02);
mg.GetXaxis()->SetLabelSize(0.05);
mg.GetXaxis()->SetLabelOffset(0.02);
mg.GetXaxis()->SetTitleOffset(1.4);

mg.GetYaxis()->SetTitle("Iota");
mg.GetYaxis()->CenterTitle();
mg.GetYaxis()->SetTitleSize(0.05);
mg.GetYaxis()->SetTickLength(-0.013);
mg.GetYaxis()->SetLabelSize(0.05);
mg.GetYaxis()->SetLabelOffset(0.013);
mg.GetYaxis()->SetTitleOffset(.9);
*/
c->Modified();
}