{

Double_t TMar = 0.05;
Double_t BMar = 0.2;
Double_t RMar = 0.03;
Double_t LMar = 0.01;
Double_t Height = 1.0;
Double_t Total;
Total = 5*Height+TMar+BMar;

Int_t n = 1000;
Double_t xdmu[n], yd[n];
Double_t xpmu[n], yp[n];
Double_t xd1dk[n], xp1dk[n];
Double_t xd2dk[n], xp2dk[n];
Double_t xd3dk[n], xp3dk[n];
Double_t xd4dk[n], xp4dk[n];
Double_t xd4be[n], xp4be[n];
Double_t xd5dk[n], xp5dk[n];

Int_t m = 2;
Double_t xwm[m];
Double_t yw[m];
Double_t xmb1[m];
Double_t xmb2[m];
Double_t x1dk[m];
Double_t x1b1[m];
Double_t x1b2[m];
Double_t x2dk[m];
Double_t x2b1[m];
Double_t x2b2[m];
Double_t x3dk[m];
Double_t x3b1[m];
Double_t x3b2[m];
Double_t x4dk[m];
Double_t x4b1[m];
Double_t x4b2[m];
Double_t x4be[m];
Double_t x4bb1[m];
Double_t x4bb2[m];
Double_t x5dk[m];
Double_t x5b1[m];
Double_t x5b2[m];
Double_t yb[m];
xwm[0] = -0.050;
xwm[1] = 1.0;
yw[0] = 0.0;
yw[1] = 0.0;
xmb1[0] = -0.050;
xmb1[1] = -0.050;
xmb2[0] = 1.0;
xmb2[1] = 1.0;

x1dk[0] = 4.7;
x1dk[1] = 5.7;
x1b1[0] = 4.7;
x1b1[1] = 4.7;
x1b2[0] = 5.7;
x1b2[1] = 5.7;

x2dk[0] = 0.33;
x2dk[1] = 1.33;
x2b1[0] = 0.33;
x2b1[1] = 0.33;
x2b2[0] = 1.33;
x2b2[1] = 1.33;

x3dk[0] = 0.03;
x3dk[1] = 1.03;
x3b1[0] = 0.03;
x3b1[1] = 0.03;
x3b2[0] = 1.03;
x3b2[1] = 1.03;

x4dk[0] = 7.35;
x4dk[1] = 8.35;
x4b1[0] = 7.35;
x4b1[1] = 7.35;
x4b2[0] = 8.35;
x4b2[1] = 8.35;

x4be[0] = 2.55;
x4be[1] = 3.55;
x4bb1[0] = 2.55;
x4bb1[1] = 2.55;
x4bb2[0] = 3.55;
x4bb2[1] = 3.55;

x5dk[0] = 3.63;
x5dk[1] = 4.63;
x5b1[0] = 3.63;
x5b1[1] = 3.63;
x5b2[0] = 4.63;
x5b2[1] = 4.63;

yb[0] = 0.0;
yb[1] = 200.;

TGraph *gwm = new TGraph(m,xwm,yw);
TGraph *gwmb1 = new TGraph(m,xmb1,yb);
TGraph *gwmb2 = new TGraph(m,xmb2,yb);
gwm->SetLineWidth(6);
gwm->SetLineColor(28);
gwmb1->SetLineWidth(4);
gwmb1->SetLineColor(28);
gwmb2->SetLineWidth(4);
gwmb2->SetLineColor(28);

TGraph *gw1dk = new TGraph(m,x1dk,yw);
TGraph *gw1b1 = new TGraph(m,x1b1,yb);
TGraph *gw1b2 = new TGraph(m,x1b2,yb);
gw1dk->SetLineWidth(6);
gw1dk->SetLineColor(6);
gw1b1->SetLineWidth(4);
gw1b1->SetLineColor(6);
gw1b2->SetLineWidth(4);
gw1b2->SetLineColor(6);

TGraph *gw2dk = new TGraph(m,x2dk,yw);
TGraph *gw2b1 = new TGraph(m,x2b1,yb);
TGraph *gw2b2 = new TGraph(m,x2b2,yb);
gw2dk->SetLineWidth(6);
gw2dk->SetLineColor(6);
gw2b1->SetLineWidth(4);
gw2b1->SetLineColor(6);
gw2b2->SetLineWidth(4);
gw2b2->SetLineColor(6);

TGraph *gw3dk = new TGraph(m,x3dk,yw);
TGraph *gw3b1 = new TGraph(m,x3b1,yb);
TGraph *gw3b2 = new TGraph(m,x3b2,yb);
gw3dk->SetLineWidth(6);
gw3dk->SetLineColor(6);
gw3b1->SetLineWidth(4);
gw3b1->SetLineColor(6);
gw3b2->SetLineWidth(4);
gw3b2->SetLineColor(6);

TGraph *gw4dk = new TGraph(m,x4dk,yw);
TGraph *gw4b1 = new TGraph(m,x4b1,yb);
TGraph *gw4b2 = new TGraph(m,x4b2,yb);
gw4dk->SetLineWidth(6);
gw4dk->SetLineColor(6);
gw4b1->SetLineWidth(4);
gw4b1->SetLineColor(6);
gw4b2->SetLineWidth(4);
gw4b2->SetLineColor(6);

TGraph *gw4be = new TGraph(m,x4be,yw);
TGraph *gw4bb1 = new TGraph(m,x4bb1,yb);
TGraph *gw4bb2 = new TGraph(m,x4bb2,yb);
gw4be->SetLineWidth(6);
gw4be->SetLineColor(3);
gw4bb1->SetLineWidth(4);
gw4bb1->SetLineColor(3);
gw4bb2->SetLineWidth(4);
gw4bb2->SetLineColor(3);

TGraph *gwm5 = new TGraph(m,xwm,yw);
TGraph *gw5dk = new TGraph(m,x5dk,yw);
TGraph *gw5b1 = new TGraph(m,x5b1,yb);
TGraph *gw5b2 = new TGraph(m,x5b2,yb);
TGraph *gw5bdk = new TGraph(m,x5dk,yw);
TGraph *gw5bb1 = new TGraph(m,x5b1,yb);
TGraph *gw5bb2 = new TGraph(m,x5b2,yb);
gwm5->SetLineWidth(3);
gwm5->SetLineColor(28);
gw5dk->SetLineWidth(3);
gw5dk->SetLineColor(6);
gw5dk->SetLineStyle(2);
gw5b1->SetLineWidth(4);
gw5b1->SetLineColor(6);
gw5b1->SetLineStyle(2);
gw5b2->SetLineWidth(4);
gw5b2->SetLineColor(6);
gw5b2->SetLineStyle(2);
gw5bdk->SetLineWidth(3);
gw5bdk->SetLineColor(3);
gw5bb1->SetLineWidth(4);
gw5bb1->SetLineColor(3);
gw5bb2->SetLineWidth(4);
gw5bb2->SetLineColor(3);


for(Int_t i=0; i<n; i++){
  xdmu[i] = (i-25.)/1000.+0.020;
  xd1dk[i] = (i-25.)/1000.+4.75;
  xd2dk[i] = (i-25.)/1000.+0.38;
  xd3dk[i] = (i-25.)/1000.+0.080;
  xd4dk[i] = (i-25.)/1000.+7.4;
  xd4be[i] = (i-25.)/1000.+2.6;
  xd5dk[i] = (i-25.)/1000.+3.68;
  if(i < 50){
    yd[i] = (1000./pow(50.,2))*pow(i,2);
  }
  else if(i < 1000){
    yd[i] = (1000./(pow(950.,2)))*pow(i-1000,2)*exp((50.-i)/950.);
  }
  else{
    yd[i] = 0.;
  }
}

for(Int_t i=0; i<n; i++){
  xpmu[i] = (i-2.)/1000.;
  xp1dk[i] = (i-2.)/1000.+4.75;
  xp2dk[i] = (i-2.)/1000.+0.38;
  xp3dk[i] = (i-2.)/1000.+0.080;
  xp4dk[i] = (i-2.)/1000.+7.4;
  xp4be[i] = (i-2.)/1000.+2.6;
  xp5dk[i] = (i-2.)/1000.+3.68;
  if(i < 5){
    yp[i] = (500./pow(5.,2))*pow(i,2);
  }
  else if(i < 100){
    yp[i] = (500./(pow(95.,2)))*pow(i-100,2)*exp((5.-i)/95.);
  }
  else{
    yp[i] = 0.;
  }
}

TH2D* hwin;
hwin = new TH2D("hwin","",1100,-1.,10.,110,0.,1100.);
hwin->SetStats(kFALSE);
hwin->GetYaxis()->SetLabelSize(0.0);
hwin->GetYaxis()->SetNdivisions(000,kTRUE);
hwin->GetXaxis()->SetLabelSize(0.2);
hwin->GetXaxis()->SetNdivisions(209,kTRUE);
TGraph *gdmu = new TGraph(n,xdmu,yd);
gdmu->SetLineWidth(3);
gdmu->SetLineColor(4);
TGraph *gd1dk = new TGraph(n,xd1dk,yd);
gd1dk->SetLineWidth(3);
gd1dk->SetLineColor(4);
TGraph *gd2dk = new TGraph(n,xd2dk,yd);
gd2dk->SetLineWidth(3);
gd2dk->SetLineColor(4);
TGraph *gd3dk = new TGraph(n,xd3dk,yd);
gd3dk->SetLineWidth(3);
gd3dk->SetLineColor(4);
TGraph *gd4dk = new TGraph(n,xd4dk,yd);
gd4dk->SetLineWidth(3);
gd4dk->SetLineColor(4);
TGraph *gd4be = new TGraph(n,xd4be,yd);
gd4be->SetLineWidth(3);
gd4be->SetLineColor(4);
TGraph *gd5dk = new TGraph(n,xd5dk,yd);
gd5dk->SetLineWidth(3);
gd5dk->SetLineColor(4);

TGraph *gpmu = new TGraph(n,xpmu,yp);
gpmu->SetLineWidth(2);
gpmu->SetLineColor(2);
TGraph *gp1dk = new TGraph(n,xp1dk,yp);
gp1dk->SetLineWidth(2);
gp1dk->SetLineColor(2);
TGraph *gp2dk = new TGraph(n,xp2dk,yp);
gp2dk->SetLineWidth(2);
gp2dk->SetLineColor(2);
TGraph *gp3dk = new TGraph(n,xp3dk,yp);
gp3dk->SetLineWidth(2);
gp3dk->SetLineColor(2);
TGraph *gp4dk = new TGraph(n,xp4dk,yp);
gp4dk->SetLineWidth(2);
gp4dk->SetLineColor(2);
TGraph *gp4be = new TGraph(n,xp4be,yp);
gp4be->SetLineWidth(2);
gp4be->SetLineColor(2);
TGraph *gp5dk = new TGraph(n,xp5dk,yp);
gp5dk->SetLineWidth(2);
gp5dk->SetLineColor(2);

TCanvas* cwin;
cwin = new TCanvas("cwin","",800,800);
cwin->SetFillColor(0);
cwin->SetBorderSize(0.0);
cwin->SetBorderMode(0);
cwin->SetFrameBorderMode(0);
cwin->Divide(1,5);

cwin->cd(1);
TPad *p1;
Double_t p1x1 = 0.0;
Double_t p1y1 = (4*Height+BMar)/Total;
Double_t p1x2 = 1.0;
Double_t p1y2 = 1.0;
p1 = gPad;
p1->SetPad("p1","",p1x1,p1y1,p1x2,p1y2);
p1->SetFillColor(0);
p1->SetBorderMode(0);
p1->SetBorderSize(0.0);
p1->SetTopMargin(TMar/(TMar + Height));
p1->SetBottomMargin(0.0);
p1->SetLeftMargin(LMar);
p1->SetRightMargin(RMar);
p1->SetFrameBorderMode(0);

cwin->cd(2);
TPad *p2;
Double_t p2x1 = 0.0;
Double_t p2y1 = (3*Height+BMar)/Total;
Double_t p2x2 = 1.0;
Double_t p2y2 = (4*Height+BMar)/Total;
p2 = gPad;
p2->SetPad("p2","",p2x1,p2y1,p2x2,p2y2);
p2->SetFillColor(0);
p2->SetBorderMode(0);
p2->SetBorderSize(0.0);
p2->SetTopMargin(0.0);
p2->SetBottomMargin(0.0);
p2->SetLeftMargin(LMar);
p2->SetRightMargin(RMar);
p2->SetFrameBorderMode(0);

cwin->cd(3);
TPad *p3;
Double_t p3x1 = 0.0;
Double_t p3y1 = (2*Height+BMar)/Total;
Double_t p3x2 = 1.0;
Double_t p3y2 = (3*Height+BMar)/Total;
p3 = gPad;
p3->SetPad("p3","",p3x1,p3y1,p3x2,p3y2);
p3->SetFillColor(0);
p3->SetBorderMode(0);
p3->SetBorderSize(0.0);
p3->SetTopMargin(0.0);
p3->SetBottomMargin(0.0);
p3->SetLeftMargin(LMar);
p3->SetRightMargin(RMar);
p3->SetFrameBorderMode(0);

cwin->cd(4);
TPad *p4;
Double_t p4x1 = 0.0;
Double_t p4y1 = (1*Height+BMar)/Total;
Double_t p4x2 = 1.0;
Double_t p4y2 = (2*Height+BMar)/Total;
p4 = gPad;
p4->SetPad("p4","",p4x1,p4y1,p4x2,p4y2);
p4->SetFillColor(0);
p4->SetBorderMode(0);
p4->SetBorderSize(0.0);
p4->SetTopMargin(0.0);
p4->SetBottomMargin(0.0);
p4->SetLeftMargin(LMar);
p4->SetRightMargin(RMar);
p4->SetFrameBorderMode(0);

cwin->cd(5);
TPad *p5;
Double_t p5x1 = 0.0;
Double_t p5y1 = 0.0;
Double_t p5x2 = 1.0;
Double_t p5y2 = (Height+BMar)/Total;
p5 = gPad;
p5->SetPad("p5","",p5x1,p5y1,p5x2,p5y2);
p5->SetFillColor(0);
p5->SetBorderMode(0);
p5->SetBorderSize(0.0);
p5->SetTopMargin(0.0);
p5->SetBottomMargin(BMar/(BMar + Height));
p5->SetLeftMargin(LMar);
p5->SetRightMargin(RMar);
p5->SetFrameBorderMode(0);

cwin->cd(1);
hwin->Draw();
gdmu->Draw("C");
gpmu->Draw("C");
gd1dk->Draw("C");
gp1dk->Draw("C");
gwm->Draw("C");
gwmb1->Draw("C");
gwmb2->Draw("C");
gw1dk->Draw("C");
gw1b1->Draw("C");
gw1b2->Draw("C");

cwin->cd(2);
hwin->Draw();
gdmu->Draw("C");
gpmu->Draw("C");
gd2dk->Draw("C");
gp2dk->Draw("C");
gwm->Draw("C");
gwmb1->Draw("C");
gwmb2->Draw("C");
gw2dk->Draw("C");
gw2b1->Draw("C");
gw2b2->Draw("C");

cwin->cd(3);
hwin->Draw();
gdmu->Draw("C");
gpmu->Draw("C");
gd3dk->Draw("C");
gp3dk->Draw("C");
gwm->Draw("C");
gwmb1->Draw("C");
gwmb2->Draw("C");
gw3dk->Draw("C");
gw3b1->Draw("C");
gw3b2->Draw("C");

cwin->cd(4);
hwin->Draw();
gdmu->Draw("C");
gpmu->Draw("C");
gd4dk->Draw("C");
gp4dk->Draw("C");
gd4be->Draw("C");
gp4be->Draw("C");
gwm->Draw("C");
gwmb1->Draw("C");
gwmb2->Draw("C");
gw4dk->Draw("C");
gw4b1->Draw("C");
gw4b2->Draw("C");
gwm->Draw("C");
gwmb1->Draw("C");
gwmb2->Draw("C");
gw4be->Draw("C");
gw4bb1->Draw("C");
gw4bb2->Draw("C");

cwin->cd(5);
hwin->Draw();
gdmu->Draw("C");
gpmu->Draw("C");
gd5dk->Draw("C");
gp5dk->Draw("C");
gwm5->Draw("C");
gwmb1->Draw("C");
gwmb2->Draw("C");
gw5bdk->Draw("C");
gw5bb1->Draw("C");
gw5bb2->Draw("C");
gw5dk->Draw("C");
gw5b1->Draw("C");
gw5b2->Draw("C");



}