{
#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 x[n]; Double_t y1[n],  y2[n], y3[n], y4[n], y5[n], z1[n],  z2[n], z3[n], z4[n], z5[n], z2B[n]; 
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();
//
g1 = new TGraph(n,x,z2);  // pick up only scenario-B, Hamburg model
g1->SetTitle();
g1->SetLineWidth(2);
g1->SetMarkerColor(2); g1->SetMarkerStyle(1); g1->SetLineColor(2);
Double_t ymax= 600., xmax = 12.;
g1->SetMaximum(ymax);
g1->GetYaxis()->SetTitle("V_{d} and Bias Voltages V_{B} (Volt)");
g1->GetXaxis()->SetTitle("LHC year");
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.2);  g1->GetXaxis()->SetTitleOffset(1.0);
g1->GetXaxis()->SetRangeUser(0.,xmax);
g1->GetXaxis()->SetNdivisions(Int_n=13);
g1->Draw("ALP");
//
ifstream fin("BiasVoltage.txt");
for (int i=0; i<n; i++) fin >> x[i] >> y2[i]>> z2[i] >> z2B[i] ;fin.close();
//
int k = 0;
for (int i=0; i<n; i++){
if(x[i]-int(x[i]) < 0.06) continue;
if(x[i]-int(x[i]) > 0.60) continue;
k++;
if(k<26)continue;
x[i-25]=x[i-26];
h0 = new TGraph(26,&x[i-25],&z2[i-25]);
h1 = new TGraph(26,&x[i-25],&z2B[i-25]);
h0->SetLineWidth(2);h0->SetMarkerColor(4); h0->SetMarkerStyle(1); 
h0->SetLineColor(4);h0->Draw("SAME LP");
h1->SetLineWidth(2);h1->SetMarkerColor(8); h1->SetMarkerStyle(1); 
h1->SetLineColor(8);h1->Draw("SAME LP");
k=0;
}
fin.close();
//
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(); 
}
//
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();
TLatex *t=new TLatex(10.08,  20,"optional years");t->SetTextSize(0.035);
t->SetTextFont(132); t->Draw();   
//
leg = new TLegend(0.12 ,0.70 ,0.60,0.88 );
leg->SetFillColor(10);leg->SetBorderSize(2);leg->SetTextFont(132);
leg->SetTextSize(0.04);
leg->SetHeader("   Setting the bias voltage (Hamburg Model)");
leg->AddEntry(g1,"Full depletion voltage","l"); 
leg->AddEntry(h0,"Vd + 100V (case A)","L"); 
leg->AddEntry(h1,"fixed 430V (case B)","L"); 
leg->Draw();   
}