void makePics(int iShow){

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

  TFile *hFileMCT = new TFile("mcThrownV3_stage3.root");
  TFile *hFileMCS = new TFile("mcSeenV3_stage3.root");
  TFile *hFileReal = new TFile("realV3_stage3.root");
  //TFile *hFileReal = new TFile("A.root");

  TH1D *hMassMCT; hFileMCT->GetObject("hMassKpKmPi0",hMassMCT);
  TH1D *hMassMCS; hFileMCS->GetObject("hMassKpKmPi0",hMassMCS);
  TH1D *hMassReal; hFileReal->GetObject("hMassKpKmPi0",hMassReal);

  TH1D *hMassRealKpStar; hFileReal->GetObject("hMassKpStar",hMassRealKpStar);
  TH1D *hMassRealKmStar; hFileReal->GetObject("hMassKmStar",hMassRealKmStar);
  TH1D *hMassRealKpKm; hFileReal->GetObject("hMassKpKm",hMassRealKpKm);
 
  
  TH1D *hCosMCT; hFileMCT->GetObject("hCosGJ",hCosMCT);
  TH1D *hCosMCS; hFileMCS->GetObject("hCosGJ",hCosMCS);
  TH1D *hCosReal; hFileReal->GetObject("hCosGJ",hCosReal);

  TH1D *hCosHMCT; hFileMCT->GetObject("hCosH",hCosHMCT);
  TH1D *hCosHMCS; hFileMCS->GetObject("hCosH",hCosHMCS);
  TH1D *hCosHReal; hFileReal->GetObject("hCosH",hCosHReal);

  TH1D *hPhiMCT; hFileMCT->GetObject("hPhiGJ",hPhiMCT);
  TH1D *hPhiMCS; hFileMCS->GetObject("hPhiGJ",hPhiMCS);
  TH1D *hPhiReal; hFileReal->GetObject("hPhiGJ",hPhiReal);

  TH1D *hPhiHMCT; hFileMCT->GetObject("hPhiH",hPhiHMCT);
  TH1D *hPhiHMCS; hFileMCS->GetObject("hPhiH",hPhiHMCS);
  TH1D *hPhiHReal; hFileReal->GetObject("hPhiH",hPhiHReal);

  TH1D *hEff = (TH1D *)hMassReal->Clone("hEff");
  TH1D *hMassRealCorr = (TH1D *)hMassReal->Clone("hMassRealCorr");
  TH1D *hMassBack = (TH1D *)hMassReal->Clone("hMassBack");
  hEff->Sumw2();
  hEff->Divide(hMassMCS,hMassMCT);
  hEff->SetMaximum(0.015);

  if (iShow ==17){
    hMassRealKpStar->GetXaxis()->SetTitle("Mass(K^{+}#pi^{0})");
    hMassRealKpStar->GetYaxis()->SetTitle("Counts");
    hMassRealKpStar->SetLineColor(2);
    hMassRealKpStar->SetLineWidth(2);
    hMassRealKpStar->Draw();
  }
  if (iShow ==18){
    hMassRealKmStar->GetXaxis()->SetTitle("Mass(K^{-}#pi^{0})");
    hMassRealKmStar->GetYaxis()->SetTitle("Counts");
    hMassRealKmStar->SetLineColor(2);
    hMassRealKmStar->SetLineWidth(2);
    hMassRealKmStar->Draw();
  }
  if (iShow ==19){
    hMassRealKpKm->GetXaxis()->SetTitle("Mass(K^{+}K^{-})");
    hMassRealKpKm->GetYaxis()->SetTitle("Counts");
    hMassRealKpKm->SetLineColor(4);
    hMassRealKpKm->SetLineWidth(2);
    hMassRealKpKm->SetMaximum(145);
    hMassRealKpKm->Draw("e");
  }

  if (iShow == 5) {
    hCosMCT->GetXaxis()->SetTitle("cos(#theta)");
    hCosMCT->GetYaxis()->SetTitle("Counts");
    hCosMCT->SetLineColor(4);
    hCosMCT->SetLineWidth(2);
    hCosMCT->GetYaxis()->SetTitleOffset(1.4);
    hCosMCT->Draw();
  }
  if (iShow == 6) {
    hCosMCS->GetXaxis()->SetTitle("cos(#theta)");
    hCosMCS->GetYaxis()->SetTitle("Counts");
    hCosMCS->SetLineColor(2);
    hCosMCS->SetLineWidth(2);
    hCosMCS->GetYaxis()->SetTitleOffset(1.4);
    hCosMCS->Draw();
  }
  if (iShow == 7) {
    hCosReal->GetXaxis()->SetTitle("cos(#theta)");
    hCosReal->GetYaxis()->SetTitle("Counts");
    hCosReal->SetLineColor(1);
    hCosReal->SetLineWidth(2);
    hCosReal->GetYaxis()->SetTitleOffset(1.4);
    hCosReal->Draw();
  }
  /////////
  if (iShow == 8) {
    hCosHMCT->GetXaxis()->SetTitle("cos(#theta_{H})");
    hCosHMCT->GetYaxis()->SetTitle("Counts");
    hCosHMCT->SetLineColor(4);
    hCosHMCT->SetLineWidth(2);
    hCosHMCT->GetYaxis()->SetTitleOffset(1.4);
    hCosHMCT->Draw();
  }
  if (iShow == 9) {
    hCosHMCS->GetXaxis()->SetTitle("cos(#theta_{H})");
    hCosHMCS->GetYaxis()->SetTitle("Counts");
    hCosHMCS->SetLineColor(2);
    hCosHMCS->SetLineWidth(2);
    hCosHMCS->GetYaxis()->SetTitleOffset(1.4);
    hCosHMCS->Draw();
  }
  if (iShow == 10) {
    hCosHReal->GetXaxis()->SetTitle("cos(#theta_{H})");
    hCosHReal->GetYaxis()->SetTitle("Counts");
    hCosHReal->SetLineColor(1);
    hCosHReal->SetLineWidth(2);
    hCosHReal->GetYaxis()->SetTitleOffset(1.4);
    hCosHReal->Draw();
  }
  ///////////////
  if (iShow == 11) {
    hPhiHMCT->GetXaxis()->SetTitle("#phi_{H}");
    hPhiHMCT->GetYaxis()->SetTitle("Counts");
    hPhiHMCT->SetLineColor(4);
    hPhiHMCT->SetLineWidth(2);
    hPhiHMCT->GetYaxis()->SetTitleOffset(1.4);
    hPhiHMCT->Draw();
  }
  if (iShow == 12) {
    hPhiHMCS->GetXaxis()->SetTitle("#phi_{H}");
    hPhiHMCS->GetYaxis()->SetTitle("Counts");
    hPhiHMCS->SetLineColor(2);
    hPhiHMCS->SetLineWidth(2);
    hPhiHMCS->GetYaxis()->SetTitleOffset(1.4);
    hPhiHMCS->Draw();
  }
  if (iShow == 13) {
    hPhiHReal->GetXaxis()->SetTitle("#phi_{H}");
    hPhiHReal->GetYaxis()->SetTitle("Counts");
    hPhiHReal->SetLineColor(1);
    hPhiHReal->SetLineWidth(2);
    hPhiHReal->GetYaxis()->SetTitleOffset(1.4);
    hPhiHReal->Draw();
  }
  /////////
  if (iShow == 14) {
    hPhiMCT->GetXaxis()->SetTitle("#phi");
    hPhiMCT->GetYaxis()->SetTitle("Counts");
    hPhiMCT->SetLineColor(4);
    hPhiMCT->SetLineWidth(2);
    hPhiMCT->GetYaxis()->SetTitleOffset(1.4);
    hPhiMCT->Draw();
  }
  if (iShow == 15) {
    hPhiMCS->GetXaxis()->SetTitle("#phi");
    hPhiMCS->GetYaxis()->SetTitle("Counts");
    hPhiMCS->SetLineColor(2);
    hPhiMCS->SetLineWidth(2);
    hPhiMCS->GetYaxis()->SetTitleOffset(1.4);
    hPhiMCS->Draw();
  }
  if (iShow == 16) {
    hPhiReal->GetXaxis()->SetTitle("#phi");
    hPhiReal->GetYaxis()->SetTitle("Counts");
    hPhiReal->SetLineColor(1);
    hPhiReal->SetLineWidth(2);
    hPhiReal->GetYaxis()->SetTitleOffset(1.4);
    hPhiReal->Draw();
  }
  /////////
  if (iShow == 1) {
    hMassMCS->GetXaxis()->SetTitle("Mass/GeV");
    hMassMCS->GetYaxis()->SetTitle("Counts");
    hMassMCS->GetYaxis()->SetTitleOffset(1.4);
    hMassMCS->GetXaxis()->SetRangeUser(1.21,1.33);
    hMassMCS->SetMinimum(0);
    hMassMCS->Draw();
  }
  if (iShow == 2) {
    hMassMCT->GetXaxis()->SetTitle("Mass/GeV");
    hMassMCT->GetYaxis()->SetTitle("Counts");
    hMassMCT->GetYaxis()->SetTitleOffset(1.4);
    hMassMCT->GetXaxis()->SetRangeUser(1.21,1.33);
    hMassMCT->SetMinimum(0);
    hMassMCT->Draw();
  }
  if (iShow == 3) {
    hEff->GetXaxis()->SetTitle("Mass/GeV");
    hEff->GetYaxis()->SetTitle("Efficiency");
    hEff->GetYaxis()->SetTitleOffset(1.4);
    hEff->GetXaxis()->SetRangeUser(1.21,1.33);
    hEff->Draw();
  }
  hMassRealCorr->Sumw2();
  hMassRealCorr->Divide(hMassReal,hEff);

  hMassRealCorr->GetXaxis()->SetTitle("Mass/GeV");
  hMassRealCorr->GetYaxis()->SetTitle("Y/#epsilon");
  hMassRealCorr->GetYaxis()->SetTitleOffset(1.5);
  hMassRealCorr->GetXaxis()->SetRangeUser(1.21,1.33);

  TF1 *f1 = new TF1("f1","gaus(0)+pol1(3)");
  TF1 *f0 = new TF1("f0","gaus");
  f1->SetLineWidth(2);
  f1->SetLineColor(4);
  if (iShow == 4) {
    f0->SetParameter(0,8000);
    f0->SetParameter(1,1.28);
    f0->SetParameter(2,0.01);
    hMassRealCorr->Fit("f0","B","R",1.24,1.3);
    double par0 = f0->GetParameter(0);
    double par1 = f0->GetParameter(1);
    double par2 = f0->GetParameter(2);
    f1->SetParameter(0,par0);
    f1->SetParameter(1,par1);
    f1->SetParameter(2,par2);
    hMassRealCorr->Fit("f1","B","R",1.24,1.31);

    f1->SetParameter(0,0.0);
    hMassBack->Eval(f1);


    hMassRealCorr->Draw();
    hMassBack->GetXaxis()->SetRangeUser(1.24,1.31);
    hMassBack->SetLineColor(2);
    hMassBack->SetLineWidth(2);
    hMassBack->Draw("csame");
    int binLow = hMassBack->FindBin(1.24);
    int binHigh = hMassBack->FindBin(1.31);
    cout<<"binLow, binHigh = "<<binLow<<" , "<<binHigh<<endl;
    double yAll =  hMassRealCorr->Integral(binLow,binHigh);
    double yTot =  hMassRealCorr->Integral();
    double yBack =  hMassBack->Integral(binLow,binHigh);
    double ySig = yAll-yBack;
    double yBA = yBack/yAll;
    double ySA = ySig/yAll;
    double yST = ySig/yTot;
    cout<<"yBack/yAll = "<<yBA<<endl;
    cout<<"ySig/yAll = "<<ySA<<endl;
    cout<<"ySig/yTot = "<<yST<<endl;
  }
}