void multigraph()
{
   gStyle->SetOptFit();
   auto c1 = new TCanvas("c1","multigraph",1600,1000);
   c1->SetGrid();
   auto *mg = new TMultiGraph();

   const Int_t n = 40;
   Double_t x[]  = {0,0,1307,2338,2514,2441,2297,2139,1987,1849,1724,1612,1512,1423,1343,1271,1206,1147,1093,1044,999,957,919,884,851,821,
                    792,766,741,718,696,675,656,637,620,603,588,573,559};
   Double_t x1[] = {19073,10564,7271,5524,4447,3719,3194,2799,2490,2243,2040,1871,1727,1604,1497,1404,1322,1248,1183,1124,1070,1022,977,937,
                   899,864,832,803,775,749,725,702,681,661,642,624,607,591,576,562};
   Double_t y[]  = {1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40};

   auto g2 = new TGraph(n, x1, y);
   g2->SetPoint(n,x1[n-1],y[0]);

   mg->Add(g2);
   g2->SetFillColor(kGreen-2);

   auto g1 = new TGraph(n, x, y);
  // g1->SetFillColor(kRed-2);
   g1->SetPoint(n,x[0],y[n-1]);
   mg->Add(g1);

   mg->SetTitle("; m^{2}_{12};tan#beta");
  // mg->GetXaxis()->CenterTitle(true);
  // mg->GetYaxis()->CenterTitle(true);
   mg->Draw("ACF");

   c1->SetCanvasSize(1200,1000);
   //c1->RedrawAxis("G");
   mg->GetHistogram()->DrawCopy("sameaxig");
   //c1->SaveAs("multigraph.png");
}