#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();        
c->SetGrid();


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




g1->SetLineColor(kAzure-3);
g2->SetLineColor(kBlack);


auto mg = new TMultiGraph();
mg->Add(g1.GetPtr());
mg->Add(g2.GetPtr());

//c->SetLogx();


mg->SetTitle("Development of Iota");
mg->GetYaxis()->SetRangeUser(0.86,1.17);
//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.03);
//mg.GetXaxis()->SetTickLength(-0.02);
//mg.GetXaxis()->SetLabelSize(0.02);
//mg.GetXaxis()->SetLabelOffset(0.02);
//mg.GetXaxis()->SetTitleOffset(1.4);

mg->GetYaxis()->SetTitle("Iota");
//mg.GetYaxis()->CenterTitle();
//mg.GetYaxis()->SetTitleSize(0.04);
//mg.GetYaxis()->SetTickLength(-0.013);
//mg.GetYaxis()->SetLabelSize(0.03);
//mg.GetYaxis()->SetLabelOffset(0.013);
//mg.GetYaxis()->SetTitleOffset(.9);

mg->DrawClone("A");

auto legend = new TLegend(0.15,0.72,0.40,0.85);
legend->AddEntry(g2,"Reference Case");
//legend->AddEntry(g1,"New Optimization","l");
//legend->AddEntry(GOrig,"Orig. fixed boundary","l");
legend->SetMargin(0.38);
legend->SetTextSize(0.025);
legend->Draw();



c->Modified();
}