#include "TCanvas.h"
#include "TROOT.h"
#include "TLegend.h"
#include "AtlasStyle.C"

void band() {

  gStyle->SetOptStat(0);  
 
  TCanvas* canvas = new TCanvas; 
  canvas->cd(); canvas->SetTopMargin(0.05); canvas->SetBottomMargin(0.14); canvas->SetLeftMargin(0.13); canvas->SetRightMargin(0.035); 

  Double_t mass[17]; Double_t median[17]; Double_t minus[17]; Double_t plus[17]; 
   
  //====| mass points |=====================
  mass[0] = 3.0;  mass[1] = 3.5;  mass[2] = 4.0;  mass[3] = 4.5;  mass[4] = 5.0;  mass[5] = 5.5;  mass[6] = 6.0;  mass[7] = 6.5; mass[8] = 7.0; mass[9] = 6.5;
  mass[10] = 6.0; mass[11] = 5.5; mass[12] = 5.0; mass[13] = 4.5; mass[14] = 4.0; mass[15] = 3.5; mass[16] = 3.0; 

  //====| expected limit |=====================
  median[0] = 0.0280;  median[1] = 0.0289;  median[2] = 0.0385;  median[3] = 0.0602;  median[4] = 0.0772;  median[5] = 0.0972;  median[6] = 0.1461;
  median[7] = 0.2502;  median[8] = 0.5000;  median[9] = 0.7499;  median[10] = 0.8540; median[11] = 0.9029; median[12] = 0.9229; median[13] = 0.9399;
  median[14] = 0.9615; median[15] = 0.9711; median[16] = 0.9721;

  TGraph* gr_med = new TGraph(17,mass,median); gr_med->SetLineWidth(2); gr_med->SetLineStyle(7); 
 
  //====| -/+ error band |=====================
  minus[0] = 0.0140;  minus[1] = 0.0149;  minus[2] = 0.0195;  minus[3] = 0.0302;  minus[4] = 0.0402;  minus[5] = 0.0552;  minus[6] = 0.0801; 
  minus[7] = 0.1472;  minus[8] = 0.1470;  minus[9] = 0.1471;  minus[10] = 0.0800; minus[11] = 0.0521; minus[12] = 0.0402; minus[13] = 0.0301; 
  minus[14] = 0.0195; minus[15] = 0.0149; minus[16] = 0.0139; 

  plus[0] = 0.0310;  plus[1] = 0.0321; plus[2] = 0.0475;  plus[3] = 0.0768;  plus[4] = 0.1068;  plus[5] = 0.1558;  plus[6] = 0.2079;  plus[7] = 0.1532; 
  plus[8] = 0.1532;  plus[9] = 0.1531; plus[10] = 0.2080; plus[11] = 0.1559; plus[12] = 0.1068; plus[13] = 0.0769; plus[14] = 0.0435; plus[15] = 0.0321;
  plus[16] = 0.0311;

  TGraphAsymmErrors* gr_band = new TGraphAsymmErrors(17,mass,median,0,0,minus,plus);  
  gr_band->GetYaxis()->SetRangeUser(0.0,1.0); gr_band->GetXaxis()->SetRangeUser(300.,1100.); 
  gr_band->GetXaxis()->SetTitle("var_{1}");   gr_band->GetYaxis()->SetTitle("var_{2}");
  gr_band->GetYaxis()->SetTitleOffset(1.2);   gr_band->GetXaxis()->SetTitleOffset(1.2); 
  gr_band->SetLineColor(kGreen+3);            gr_band->SetFillColor(kGreen+3);

  gr_band->Draw("AC 3"); 
  gr_med->Draw("Same C"); 
  
  //====| Legend |=====================
  TLegend* leg1 = new TLegend(0.470,0.525,0.075,0.810); 
  TLegend* leg2 = new TLegend(0.180,0.400,0.520,0.580); 
  
  leg1->AddEntry((TObject*)0, "graph with assymetric error band", "");
  leg1->SetTextSize(0.05); leg1->SetFillColor(0); leg1->SetFillStyle(0); leg1->SetShadowColor(2); leg1->SetBorderSize(0); leg1->Draw(); 

  leg2->AddEntry(gr_med,"median","l"); leg2->AddEntry(gr_band,"error band","f");
  leg2->SetTextSize(0.038); leg2->SetFillColor(0); leg2->SetFillStyle(0); leg2->SetShadowColor(2); leg2->SetBorderSize(0); leg2->Draw();

  gPad->Update(); canvas->RedrawAxis(); canvas->SaveAs("band.png"); canvas->Update(); cout << " " << endl;
}