{
#include <fstream.h>
#include <math.h>
c1 = new TCanvas("c1","A Simple Graph with error bars",0,0,700,500);
c1->SetFillStyle(4100);
c1->SetGrid();
c1->GetFrame()->SetBorderSize(12);
TStyle *default = new TStyle("Default","Default Style");
//default->SetPadBorderMode(0);
default->SetCanvasColor(0);
default->cd(); // this becomes now the current style gStyle
// create the arrays for the points
Int_t n = 409;
Double_t year[n], neff1[n], neff2[n], neff3[n], neff4[n], nefftotal[n];
////////////
ifstream fin("Neff_Hamburg.txt");
//ifstream fin("Neff_TDR.txt");
for (int i=0; i<n; i++){
fin >> year[i] >> neff1[i]>>neff2[i]>>neff3[i]>>neff4[i]>>nefftotal[i]; 
neff1[i]=neff1[i]/1.E12; nefftotal[i]=nefftotal[i]/1.E12;
neff2[i]=-neff2[i]/1.E12; neff3[i]=-neff3[i]/1.E12; neff4[i]=-neff4[i]/1.E12;
}
fin.close();
g1 = new TGraph(n,year,nefftotal);
g2 = new TGraph(n,year,neff1);
g3 = new TGraph(n,year,neff2);
g4 = new TGraph(n,year,neff3);
g5 = new TGraph(n,year,neff4);
//
g1->SetTitle();
g1->SetLineWidth(2);
g1->SetMarkerColor(1); g1->SetMarkerStyle(1); g1->SetLineColor(1);
g2->SetLineColor(2);g3->SetLineColor(3);g4->SetLineColor(4);g5->SetLineColor(6);
g2->SetMarkerColor(2);g3->SetMarkerColor(3);g4->SetMarkerColor(4);g5->SetMarkerColor(6);
Double_t ymax= 2., xmax=12.;
Double_t ymin= -8.;
g1->SetMaximum(ymax);
g1->SetMinimum(ymin);
g1->GetYaxis()->SetTitle("N_{eff} / 10^{12} [cm^{-3}]");
g1->GetYaxis()->SetTitleSize(0.04); g1->GetYaxis()->SetTitleFont(132);
g1->GetYaxis()->SetLabelSize(0.04); g1->GetYaxis()->SetLabelFont(132);
g1->GetXaxis()->SetTitleSize(0.04); g1->GetXaxis()->SetTitleFont(132);
g1->GetXaxis()->SetLabelSize(0.04); g1->GetXaxis()->SetLabelFont(132);
g1->GetYaxis()->SetTitleOffset(1.0);
g1->GetXaxis()->SetTitle("LHC Years");
g1->GetXaxis()->SetTitleOffset(1.0);
g1->GetXaxis()->SetRangeUser(0.,xmax);
g1->GetXaxis()->SetNdivisions(Int_n=13);
g1->Draw("ALP");
g2->Draw("SAME LP");
g3->Draw("SAME LP");
g4->Draw("SAME LP");
g5->Draw("SAME LP");
//
TLine *b= new TLine(0,ymin, 0,ymax); b->SetLineWidth(2); b->Draw();
TLine *b= new TLine(0,ymin, xmax,ymin); b->SetLineWidth(2); b->Draw();
TLine *b= new TLine(0,ymax, xmax,ymax); b->SetLineWidth(2); b->Draw();
TLine *b= new TLine(xmax,ymin, xmax,ymax); b->SetLineWidth(2); b->Draw();
TLine *b= new TLine(10.,ymin, 10.,ymax); b->SetLineWidth(1); b->Draw();
TText *t=new TText(10.08,  -7.6,"optional years");t->SetTextSize(0.035);
t->SetTextFont(132); t->Draw();   
//
leg = new TLegend(0.15 ,0.15 ,0.4,0.4);
leg->SetFillColor(10);leg->SetBorderSize(2);leg->SetTextFont(132);
leg->SetHeader("Hamburg Model");
//leg->SetHeader("TDR Model");
leg->AddEntry(g2,"Intrinsic donor","L");
leg->AddEntry(g3,"Stable acceptor","L");
leg->AddEntry(g4,"Unstable acceptor","l");
leg->AddEntry(g5,"Reverse annealing","l");
leg->AddEntry(g1,"Total","l");
leg->Draw(); 
}