{
#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 = 13;
Double_t x[n],  vd1[n],  vd2[n], current[n], power1[n], power2[n]; 
ifstream fin("maintenance_B3.txt");
for (int i=0; i<n; i++) {
fin >> x[i] >> vd1[i]>> vd2[i]>>current[i]>>power1[i]>>power2[i]; 
current[i]=current[i]*100.;}
fin.close();
g1 = new TGraph(n,x, vd2);
g2 = new TGraph(n,x, current);
g3 = new TGraph(n,x, power2);
g1->SetTitle();g1->SetLineColor(4);
g1->SetMarkerColor(4); g1->SetMarkerStyle(8); g1->SetMarkerSize(1.0); 
g2->SetLineColor(2); g2->SetMarkerColor(2); g2->SetMarkerStyle(8);
g2->SetMarkerSize(1.0); 
g3->SetLineColor(1); g3->SetMarkerColor(1); g3->SetMarkerStyle(8);
g3->SetMarkerSize(1.0); 
Double_t ymax= 300., xmax=60., ymin=0;
g1->SetMaximum(ymax);
g1->SetMinimum(ymin);
g1->GetYaxis()->SetTitle("V_{d} [V] and power density [#muW/mm^{2}]");
g1->GetXaxis()->SetTitle("maintenance days per 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=7);
g1->Draw("ALP");
g2->Draw("SAME LP");
g3->Draw("SAME LP");
//draw an axis on the right side
double rightmax=5.;
TGaxis *axis = new TGaxis(60.,0.,60.,ymax,0,rightmax,520,"+L");
axis->SetLineColor(1);axis->SetLabelColor(1);
axis->SetLabelSize(0.04);axis->SetLabelFont(132);
axis->SetTitleSize(0.04);axis->SetTitleFont(132);
axis->SetTitle("leak current / module [mA]");
axis->Draw();
//----------------------------------------------------
TLine *z=new TLine( 0,0,xmax,0);z->SetLineWidth(2); z->Draw();
TLine *z=new TLine(0,0,0,ymax);z->SetLineWidth(2);z->Draw();
TLine *z=new TLine(0,ymax,xmax,ymax);z->SetLineWidth(2);z->Draw();
TLine *z=new TLine(xmax,0,xmax,ymax);z->SetLineWidth(2);z->Draw();
//
leg = new TLegend(0.12 ,0.65 ,0.55,0.88 );
leg->SetFillColor(10);leg->SetBorderSize(2);leg->SetTextFont(132);
leg->SetTextSize(0.04);
leg->SetHeader("Dependence on maintenance length");
leg->AddEntry(v,"values at the end of 10th year",""); 
leg->AddEntry(g1,"full depletion voltage","l"); 
leg->AddEntry(g2,"leak current (right scale)","L"); 
leg->AddEntry(g3,"power density at 0^{o}C","L"); 
leg->Draw();  
}