void sys(int iShow) {

  gROOT->Reset();
  gROOT->SetStyle("Plain");
  gStyle->SetOptStat(0);

  TFile *hFile = new TFile("myOutFile.root");

  TH1D *hYStar1;  hFile->GetObject("hYStar1",hYStar1);
  TH1D *hYStar2;  hFile->GetObject("hYStar2",hYStar2);

  TH1D *hYStarLow;  hFile->GetObject("hYStarLow",hYStarLow);
  TH1D *hYStarHigh;  hFile->GetObject("hYStarHigh",hYStarHigh);

  TH1D *hYStarSlow;  hFile->GetObject("hYStarSlow",hYStarSlow);
  TH1D *hYStarFast;  hFile->GetObject("hYStarFast",hYStarFast);

  TH1D *hYStarLowT;  hFile->GetObject("hYStarLowT",hYStarLowT);
  TH1D *hYStarHighT;  hFile->GetObject("hYStarHighT",hYStarHighT);

  TH1D *hYStarSlowT;  hFile->GetObject("hYStarSlowT",hYStarSlowT);
  TH1D *hYStarFastT;  hFile->GetObject("hYStarFastT",hYStarFastT);

  TH1D* hDiv11B = (TH1D*)hYStarLow->Clone("hDiv11B");
  TH1D* hDiv21B = (TH1D*)hYStarLow->Clone("hDiv21B");
  TH1D* hDiv31B = (TH1D*)hYStarLow->Clone("hDiv31B");
  TH1D* hDiv41B = (TH1D*)hYStarLow->Clone("hDiv41B");

  hDiv11B->Reset();
  hDiv21B->Reset();
  hDiv31B->Reset();
  hDiv41B->Reset();

  TH1D* hDiv = (TH1D*)hYStarLow->Clone("hDiv");
  TH1D* hDiff = (TH1D*)hYStarLow->Clone("hDiff");
  hDiff->Sumw2();
  hDiff->Add(hYStarLow,hYStar2,1.0,-1.0);
  hDiff->Scale(100);
  hDiv->Divide(hDiff,hYStar2);

  TH1D* hDiv11 = (TH1D*)hYStarLowT->Clone("hDiv11");
  hDiv11->Divide(hYStarLowT,hYStar2);

  TH1D* hDiv2 = (TH1D*)hYStarHigh->Clone("hDiv2");
  TH1D* hDiff2 = (TH1D*)hYStarHigh->Clone("hDiff2");
  hDiff2->Sumw2();
  hDiff2->Add(hYStarHigh,hYStar2,1.0,-1.0);
  hDiff2->Scale(100);
  hDiv2->Divide(hDiff2,hYStar2);

  TH1D* hDiv21 = (TH1D*)hYStarHighT->Clone("hDiv21");
  hDiv21->Divide(hYStarHighT,hYStar2);

  TH1D* hDiv3 = (TH1D*)hYStarSlow->Clone("hDiv3");
  TH1D* hDiff3 = (TH1D*)hYStarSlow->Clone("hDiff3");
  hDiff3->Sumw2();
  hDiff3->Add(hYStarSlow,hYStar2,1.0,-1.0);
  hDiff3->Scale(100);
  hDiv3->Divide(hDiff3,hYStar2);

  TH1D* hDiv31 = (TH1D*)hYStarSlowT->Clone("hDiv31");
  hDiv31->Divide(hYStarSlowT,hYStar2);

  TH1D* hDiv4 = (TH1D*)hYStarFast->Clone("hDiv4");
  TH1D* hDiff4 = (TH1D*)hYStarFast->Clone("hDiff4");
  hDiff4->Sumw2();
  hDiff4->Add(hYStarFast,hYStar2,1.0,-1.0);
  hDiff4->Scale(100);
  hDiv4->Divide(hDiff4,hYStar2);

  TH1D* hDiv41 = (TH1D*)hYStarFastT->Clone("hDiv41");
  hDiv41->Divide(hYStarFastT,hYStar2);

  //hDiv11->SetMaximum(1.0);
  //hDiv21->SetMaximum(1.0);
  //hDiv31->SetMaximum(1.0);
  //hDiv41->SetMaximum(1.0);

  if (iShow == 0){
    hYStar2->Draw();
  }
  hDiv->GetXaxis()->SetTitle("Invariant mass Y^{*}_{Low} GeV");
  hDiv->GetYaxis()->SetTitle("Percent Error Y^{*}_{Low}");
  hDiv->GetXaxis()->SetRangeUser(1.8,3.2);

  hDiv2->GetXaxis()->SetTitle("Invariant mass Y^{*}_{High} GeV");
  hDiv2->GetYaxis()->SetTitle("Percent Error Y^{*}_{High}");
  hDiv2->GetXaxis()->SetRangeUser(1.8,3.2);
  
  hDiv3->GetXaxis()->SetTitle("Invariant mass Y^{*}_{Slow} GeV");
  hDiv3->GetYaxis()->SetTitle("Percent Error Y^{*}_{Slow}");
  hDiv3->GetXaxis()->SetRangeUser(1.8,3.2);  

  hDiv4->GetXaxis()->SetTitle("Invariant mass Y^{*}_{Fast} GeV");
  hDiv4->GetYaxis()->SetTitle("Percent Error Y^{*}_{Fast}");
  hDiv4->GetXaxis()->SetRangeUser(1.8,3.2);

  //hDiv->Sumw2();

  //hDiv->Scale(100);
  //hDiv2->Scale(100);
  //hDiv3->Scale(100);
  //hDiv4->Scale(100);

  if (iShow == 1){
    hDiv->Draw();
  }
  if (iShow == 2){

    hDiv2->Draw();
  }
  if (iShow == 3){

    hDiv3->Draw();
  }
  if (iShow == 4){

    hDiv4->Draw();
  }

  /////////
  if (iShow == 11){
    hDiv11->Draw();
  }
  if (iShow == 21){
    hDiv21->Draw();
  }
  if (iShow == 31){
    hDiv31->Draw();
  }
  if (iShow == 41){
    hDiv41->Draw();
  }

  if (iShow == -1){
    c1->Clear();
    c1->Divide(2,2);
    c1->cd(1); hDiv->Draw("hist");
    c1->cd(2); hDiv2->Draw("hist");
    c1->cd(3); hDiv3->Draw("hist");
    c1->cd(4); hDiv4->Draw("hist");
    
  }

  int nBins = hDiv11->GetNbinsX();
  for (int iBin =1; iBin <= nBins; iBin++) {
    double f11 = hDiv11->GetBinContent(iBin);
    double f21 = hDiv21->GetBinContent(iBin);
    double f31 = hDiv31->GetBinContent(iBin);
    double f41 = hDiv41->GetBinContent(iBin);
    
    if (f11 > f21 && f11 > f31 && f11 > f41) {
      hDiv11B->SetBinContent(iBin,f11);
      
    }
    if (f21 > f11 && f21 > f31 && f21 > f41) {
      hDiv21B->SetBinContent(iBin,f21);
    }
    if (f31 > f11 && f31 > f21 && f31 > f41) {
      hDiv31B->SetBinContent(iBin,f31);
    }
    if (f41 > f11 && f41 > f21 && f41 > f31) {
      hDiv41B->SetBinContent(iBin,f41);
    }

  }

  if (iShow == -2){
    c1->Clear();
    c1->Divide(2,2);
    c1->cd(1); hDiv11B->Draw();
    c1->cd(2); hDiv21B->Draw();
    c1->cd(3); hDiv31B->Draw();
    c1->cd(4); hDiv41B->Draw();
    
  }

}