{
#include <fstream.h>
#include <math.h>
c1 = new TCanvas("c1","c1",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 x[n], dummy[n]={n*0.};
Double_t y1[n],  y2[n], y3[n], y4[n], y5[n], z1[n],  z2[n], z3[n], z4[n], z5[n],dmy; 
ifstream fin("Vd_scenarios.txt");
for (int i=0; i<n; i++)fin >> x[i] >> y1[i]>> y2[i]>>y3[i]>>y4[i]>>y5[i]>>z1[i]>> z2[i]>>z3[i]>>z4[i]>>z5[i]; fin.close();
gd = new TGraph(n,x,dummy);
g1 = new TGraph(n,x,y1);
g2 = new TGraph(n,x,y2);
g3 = new TGraph(n,x,y3);
g4 = new TGraph(n,x,y4);
g5 = new TGraph(n,x,y5);
h1 = new TGraph(n,x,z1);
h2 = new TGraph(n,x,z2);
h3 = new TGraph(n,x,z3);
h4 = new TGraph(n,x,z4);
h5 = new TGraph(n,x,z5);
//
ifstream fin("newest-intL.txt");
for (int i=0; i<n; i++)fin >> x[i] >> y1[i]; fin.close();
g0 = new TGraph(n,x,y1);
//
Double_t ymax= 350., xmax=12.0;
gd->SetMaximum(ymax);
gd->SetTitle(0);
gd->GetYaxis()->SetTitle("V_{D} (Volt) and  Int Lum (fb^{-1} )");
gd->GetYaxis()->SetTitleSize(0.04); gd->GetYaxis()->SetTitleFont(132);
gd->GetYaxis()->SetLabelSize(0.04); gd->GetYaxis()->SetLabelFont(132);
gd->GetXaxis()->SetTitleSize(0.04); gd->GetXaxis()->SetTitleFont(132);
gd->GetXaxis()->SetLabelSize(0.04); gd->GetXaxis()->SetLabelFont(132);
gd->GetYaxis()->SetTitleOffset(1.2);gd->GetXaxis()->SetTitleOffset(1.0);
gd->GetXaxis()->SetTitle("LHC Year");
gd->GetXaxis()->SetRangeUser(0.,xmax);
gd->GetXaxis()->SetNdivisions(Int_n=13);
gd->Draw("ALP");
//
for (i=0;i<12; i++) {
box = new TBox(i+0.14,0.,i+0.45, ymax); 
box->SetFillColor(13); box->SetFillStyle(3003); box->Draw(); 
}
//
g1->SetTitle();
g1->SetMarkerColor(2); g1->SetMarkerStyle(1); g1->SetLineColor(2);
g2->SetMarkerColor(3); g2->SetMarkerStyle(1); g2->SetLineColor(3);
g3->SetMarkerColor(4); g3->SetMarkerStyle(1); g3->SetLineColor(4);
g4->SetMarkerColor(6); g4->SetMarkerStyle(1); g4->SetLineColor(6);
g5->SetMarkerColor(1); g5->SetMarkerStyle(1); g5->SetLineColor(1);
h1->SetMarkerColor(2); h1->SetMarkerStyle(1); h1->SetLineColor(2);
h2->SetMarkerColor(3); h2->SetMarkerStyle(1); h2->SetLineColor(3);
h3->SetMarkerColor(4); h3->SetMarkerStyle(1); h3->SetLineColor(4);
h4->SetMarkerColor(6); h4->SetMarkerStyle(1); h4->SetLineColor(6);
h5->SetMarkerColor(1); h5->SetMarkerStyle(1); h5->SetLineColor(1);
g0->SetMarkerColor(1); g0->SetMarkerStyle(1); g0->SetLineColor(1);
g0->SetLineStyle(5);
h2->Draw("SAME LP");
h3->Draw("SAME LP");
h4->Draw("SAME LP");
h5->Draw("SAME LP");
h1->Draw("SAME LP");
g0->Draw("SAME LP");
//
b= new TLine(0,0, 0,ymax); b->SetLineWidth(2); b->Draw();
b= new TLine(0,0, xmax,0); b->SetLineWidth(2); b->Draw();
b= new TLine(0,ymax, xmax,ymax); b->SetLineWidth(2); b->Draw();
b= new TLine(xmax,0, xmax,ymax); b->SetLineWidth(2); b->Draw();
b= new TLine(10.,0, 10.,ymax); b->SetLineWidth(1); b->Draw();
t=new TLatex(10.08,  20,"optional years");t->SetTextSize(0.035);
t->SetTextFont(132); t->Draw();   
//
p = new TPave(0.2,225,7.2,345); p->SetFillColor(10); 
p->SetLineColor(1);p->SetBorderSize(2); p->Draw(); 
//
Double_t h=330; 
Double_t q1=0.5, q2=0.8, q3=1.1, q4=1.5 , q5=4.2, q6=2.4, q7=4.5, q8=6.3; 
t=new TLatex(q1-0.2, h,"cooling");
t->SetTextFont(132); t->SetTextSize(0.035); t->Draw();  
t=new TLatex(q4, h,"----first 3 yrs----"); 
t->SetTextFont(132); t->SetTextSize(0.035); t->Draw();   
t=new TLatex(q5, h,"----next 9 yrs----"); 
t->SetTextFont(132); t->SetTextSize(0.035); t->Draw();   
//
h = 310;
t=new TLatex(q1-0.2, h,"scenario");
t->SetTextFont(132); t->SetTextSize(0.035); t->Draw();   
t=new TLatex(q4, h,"50 116 50  23 126  50 116 50  23 126 days");
t->SetTextFont(132); t->SetTextSize(0.035); t->Draw();   
//
h =  290; 
u=new TLatex(q3,h,"A   -7  -7  -7  20 -22   -7  -7  -7  20 -22"); 
u->SetTextFont(132); u->SetTextSize(0.035); u->Draw();   
u=new TLatex(q8,h,"^{o}C"); 
u->SetTextFont(132); u->SetTextSize(0.035); u->Draw();   
xx = new TLine(q1,h+5, q2,h+5); xx->SetLineColor(2);  xx->Draw();
//
h =  275; 
u=new TLatex(q3,h,"B    0   0   0  20 -15    -7  -7  -7  20 -22"); 
u->SetTextFont(132); u->SetTextSize(0.035); u->Draw();   
u=new TLatex(q8,h,"^{o}C"); 
u->SetTextFont(132); u->SetTextSize(0.035); u->Draw();   
xx = new TLine(q1,h+5, q2,h+5); xx->SetLineColor(3);  xx->Draw();
//
h =  260; 
u=new TLatex(q3,h,"D  15 15 15  20    0    -7  -7   -7  20 -22");
u->SetTextFont(132); u->SetTextSize(0.035); u->Draw();   
u=new TLatex(q8,h,"^{o}C"); 
u->SetTextFont(132); u->SetTextSize(0.035); u->Draw();   
xx = new TLine(q1,h+5, q2,h+5); xx->SetLineColor(4);  xx->Draw();
//
h =  245; 
u=new TLatex(q3,h,"G    0   0   0  20 -15  -10 -10 -10 20 -25");
u->SetTextFont(132); u->SetTextSize(0.035); u->Draw();   
u=new TLatex(q8,h,"^{o}C"); 
u->SetTextFont(132); u->SetTextSize(0.035); u->Draw();   
xx = new TLine(q1,h+5, q2,h+5); xx->SetLineColor(6);  xx->Draw();
//
h =  230; 
u=new TLatex(q3,h,"I     0   0   0  20  -15     0    0   0  20 -15");
u->SetTextFont(132); u->SetTextSize(0.035); u->Draw();   
u=new TLatex(q8,h,"^{o}C"); 
u->SetTextFont(132); u->SetTextSize(0.035); u->Draw();   
xx = new TLine(q1,h+5, q2,h+5); xx->SetLineColor(1);  xx->Draw();
//
t=new TLatex(9.0,  85,"Hamburg model");t->SetTextSize(0.035);
t->SetTextFont(132); t->SetTextSize(0.045); t->Draw();     
t=new TLatex(4.1, 165,"integrated");t->SetTextSize(0.03);
t->SetTextFont(132); t->SetTextSize(0.045); t->Draw();   
t=new TLatex(4.3, 150,"luminosity");t->SetTextSize(0.03);
t->SetTextFont(132); t->SetTextSize(0.045); t->Draw();   
}