void makePics(int iShow){

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

  TFile *fFF = new TFile("./outFit/ff.root");
  TFile *fA  = new TFile("./A.root");

  TH1D *hReal1; fA->GetObject("hMassKpKmPi0Course1",hReal1);
  TH1D *hReal3; fA->GetObject("hMassKpKmPi0Course3",hReal3);

  TH2D *hFitFrac; fFF->GetObject("hFitFrac",hFitFrac);

  TH1D * hJ0 = (TH1D*)hFitFrac->ProjectionX("hJ0",1,1);
  //TH1D * hJ1L1S0a0 = (TH1D*)hFitFrac->ProjectionX("hJ1L1S0a0",2,10);
  TH1D * hJ1L1S0a0 = (TH1D*)hFitFrac->ProjectionX("hJ1L1S0a0",2,4);
  //TH1D * hJ1L1S0a0 = (TH1D*)hFitFrac->ProjectionX("hJ1L1S0a0",5,10);

  //TH1D * hJ1L1S0 = (TH1D*)hFitFrac->ProjectionX("hJ1L1S0",5,7);
  //TH1D * hJ1L0S1 = (TH1D*)hFitFrac->ProjectionX("hJ1L0S1",8,10);
  //TH1D * hJ1L1S1 = (TH1D*)hFitFrac->ProjectionX("hJ1L1S1",11,13);

  //TH1D * hJ1tmp = (TH1D*)hFitFrac->ProjectionX("hJ1tmp",6,11);

  TH1D * hRealJ0 = (TH1D *)hReal3->Clone("hRealJ0");
  TH1D * hRealJ1L1S0a0 = (TH1D *) hReal3->Clone("hRealJ1L1S0a0");
  //TH1D * hRealJ1L1S0 = (TH1D *) hReal3->Clone("hRealJ1L1S0");
  //TH1D * hRealJ1L0S1 = (TH1D *) hReal3->Clone("hRealJ1L0S1");
  //TH1D * hRealJ1L1S1 = (TH1D *) hReal3->Clone("hRealJ1L1S1");
  //TH1D * hRealJ1tmp = (TH1D *) hReal3->Clone("hRealJ1tmp");

  //Kill
  //hRealJ0->SetBinContent(8,0);
  //hRealJ0->SetBinError(8,0);

  //hRealJ1L1S0->SetBinContent(8,0);
  //hRealJ1L1S0->SetBinError(8,0);

  //hRealJ1L0S1->SetBinContent(8,0);
  //hRealJ1L0S1->SetBinError(8,0);

  //hRealJ1L1S1->SetBinContent(8,0);
  //hRealJ1L1S1->SetBinError(8,0);


  
  hRealJ0->Multiply(hRealJ0,hJ0);
  hRealJ0->SetLineWidth(2);
  hRealJ0->SetLineColor(2);
  hRealJ0->SetFillStyle(1001);
  hRealJ0->SetFillColor(2);

  //hRealJ1L1S0->Multiply(hRealJ1L1S0,hJ1L1S0);
  //hRealJ1L1S0->SetLineWidth(2);
  //hRealJ1L1S0->SetLineColor(2);
  //hRealJ1L1S0->SetFillStyle(1001);
  //hRealJ1L1S0->SetFillColor(2);

  hRealJ1L1S0a0->Multiply(hRealJ1L1S0a0,hJ1L1S0a0);
  hRealJ1L1S0a0->SetLineWidth(2);
  hRealJ1L1S0a0->SetLineColor(3);
  hRealJ1L1S0a0->SetFillStyle(1001);
  hRealJ1L1S0a0->SetFillColor(3);


  //hRealJ1L0S1->Multiply(hRealJ1L1S0,hJ1L1S0);
  //hRealJ1L0S1->SetLineWidth(2);
  //hRealJ1L0S1->SetLineColor(3);
  //hRealJ1L0S1->SetFillStyle(1001);
  //hRealJ1L0S1->SetFillColor(3);

  //hRealJ1tmp->Multiply(hRealJ1tmp,hJ1tmp);
  //hRealJ1tmp->SetLineWidth(2);
  //hRealJ1tmp->SetLineColor(7);
  //hRealJ1tmp->SetFillStyle(1001);
  //hRealJ1tmp->SetFillColor(7);

  //hRealJ1L1S1->Multiply(hRealJ1L1S1,hJ1L1S1);
  //hRealJ1L1S1->SetLineWidth(2);
  //hRealJ1L1S1->SetLineColor(4);
  //hRealJ1L1S1->SetFillStyle(1001);
  //hRealJ1L1S1->SetFillColor(4);

  TF1 *fg1 = new TF1("fg1","gaus",0.0,3.0);
  //TF1 *fg2 = new TF1("fg2","gaus",0.0,3.0);
  TF1 *fg2a0 = new TF1("fg2a0","gaus",0.0,3.0);
  //TF1 *fg3 = new TF1("fg3","gaus",0.0,3.0);
  //TF1 *fg4 = new TF1("fg4","gaus",0.0,3.0);
  //TF1 *fg5 = new TF1("fg5","gaus",0.0,3.0);

  fg1->SetLineColor(2);
  //fg2->SetLineColor(2);
  fg2a0->SetLineColor(3);
  //fg3->SetLineColor(3);
  //fg4->SetLineColor(4);
  //fg5->SetLineColor(7);

  double fLow = 1.24;
  double fHigh = 1.31;

  hRealJ0->Fit("fg1","R","",fLow,fHigh);
  //hRealJ1L1S0->Fit("fg2");
  hRealJ1L1S0a0->Fit("fg2a0","R","",fLow,fHigh);
  //hRealJ1L0S1->Fit("fg3");
  //hRealJ1L1S1->Fit("fg4");
  //hRealJ1tmp->Fit("fg5");

  cout<<"----------"<<endl;
  double cJ0 = 1000*fg1->GetParameter(1);
  double cJ0Err = 1000*fg1->GetParError(1);
  cout<<"cJ0 = "<<cJ0<<" +/- "<<cJ0Err<<endl;

  double wJ0 = 2.35*1000*fg1->GetParameter(2);
  double wJ0Err = 2.35*1000*fg1->GetParError(2);
  cout<<"wJ0 = "<<wJ0<<" +/- "<<wJ0Err<<endl;

  //cout<<"----------"<<endl;
  double cJ1L1S0a0 = 1000*fg2a0->GetParameter(1);
  double cJ1L1S0a0Err = 1000*fg2a0->GetParError(1);
  cout<<"cJ1L1S0a0 = "<<cJ1L1S0a0<<" +/- "<<cJ1L1S0a0Err<<endl;

  double wJ1L1S0a0 = 2.35*1000*fg2a0->GetParameter(2);
  double wJ1L1S0a0Err = 2.35*1000*fg2a0->GetParError(2);
  cout<<"wJ1L1S0a0 = "<<wJ1L1S0a0<<" +/- "<<wJ1L1S0a0Err<<endl;

  //cout<<"----------"<<endl;
  //double cJ1L1S0 = 1000*fg2->GetParameter(1);
  //double cJ1L1S0Err = 1000*fg2->GetParError(1);
  //cout<<"cJ1L1S0 = "<<cJ1L1S0<<" +/- "<<cJ1L1S0Err<<endl;

  //double wJ1L1S0 = 2.35*1000*fg2->GetParameter(2);
  //double wJ1L1S0Err = 2.35*1000*fg2->GetParError(2);
  //cout<<"wJ1L1S0 = "<<wJ1L1S0<<" +/- "<<wJ1L1S0Err<<endl;

  //cout<<"----------"<<endl;
  //double cJ1L0S1 = 1000*fg3->GetParameter(1);
  //double cJ1L0S1Err = 1000*fg3->GetParError(1);
  //cout<<"cJ1L0S1 = "<<cJ1L0S1<<" +/- "<<cJ1L0S1Err<<endl;

  //double wJ1L0S1 = 2.35*1000*fg3->GetParameter(2);
  //double wJ1L0S1Err = 2.35*1000*fg3->GetParError(2);
  //cout<<"wJ1L0S1 = "<<wJ1L0S1<<" +/- "<<wJ1L0S1Err<<endl;

  //cout<<"----------"<<endl;
  //double cJ1L1S1 = 1000*fg4->GetParameter(1);
  //double cJ1L1S1Err = 1000*fg4->GetParError(1);
  //cout<<"cJ1L1S1 = "<<cJ1L1S1<<" +/- "<<cJ1L1S1Err<<endl;

  //double wJ1L1S1 = 2.35*1000*fg4->GetParameter(2);
  //double wJ1L1S1Err = 2.35*1000*fg4->GetParError(2);
  //cout<<"wJ1L1S1 = "<<wJ1L1S1<<" +/- "<<wJ1L1S1Err<<endl;

  //cout<<"----------"<<endl;
  //double cJ1tmp = 1000*fg5->GetParameter(1);
  //double cJ1tmpErr = 1000*fg5->GetParError(1);
  //cout<<"cJ1tmp = "<<cJ1tmp<<" +/- "<<cJ1tmpErr<<endl;

  //double wJ1tmp = 2.35*1000*fg5->GetParameter(2);
  //double wJ1tmpErr = 2.35*1000*fg5->GetParError(2);
  //cout<<"wJ1tmp = "<<wJ1tmp<<" +/- "<<wJ1tmpErr<<endl;
  

  hReal3->SetMarkerStyle(20);
  hReal3->GetXaxis()->SetTitle("Mass(K^{+}K^{-}#pi^{0})/GeV");
  hReal3->GetYaxis()->SetTitle("Counts");

  hReal1->SetMarkerStyle(20);
  hReal1->GetXaxis()->SetTitle("Mass(K^{+}K^{-}#pi^{0})/GeV");
  hReal1->GetYaxis()->SetTitle("Counts");

  hRealJ0->GetXaxis()->SetTitle("Mass(K^{+}K^{-}#pi^{0})/GeV");
  hRealJ0->GetYaxis()->SetTitle("Counts");

  hRealJ1L1S0->GetXaxis()->SetTitle("Mass(K^{+}K^{-}#pi^{0})/GeV");
  hRealJ1L1S0->GetYaxis()->SetTitle("Counts");

  hRealJ1L0S1->GetXaxis()->SetTitle("Mass(K^{+}K^{-}#pi^{0})/GeV");
  hRealJ1L0S1->GetYaxis()->SetTitle("Counts");

  hRealJ1L1S1->GetXaxis()->SetTitle("Mass(K^{+}K^{-}#pi^{0})/GeV");
  hRealJ1L1S1->GetYaxis()->SetTitle("Counts");

  if (iShow == 0) {
    hReal3->Draw("e1");

  }

  hReal3->GetXaxis()->SetRangeUser(1.24,1.31);
  if (iShow == 1) {
    hReal3->Draw("e1");
    hRealJ0->Draw("same");
    //hRealJ1L1S0->Draw("same");
    hRealJ1L1S0a0->Draw("same");
    //hRealJ1L0S1->Draw("same");
    //hRealJ1L1S1->Draw("same");

    Double_t xl1=0.12, yl1 = 0.67, xl2=xl1+0.23, yl2=yl1+0.2;
    TLegend *legend = new TLegend(xl1,yl1,xl2,yl2);
    legend->SetFillColor(0);
    legend->SetLineColor(0);
    legend->SetTextSize(0.04);
    legend->AddEntry(hRealJ0,"J=0","f");
    //legend->AddEntry(hRealJ1L1S0,"J=1,L=1,S=0","f");
    //legend->AddEntry(hRealJ1L1S0a0,"J=1,L=1,S=0 a0","f");
    legend->AddEntry(hRealJ1L1S0a0,"J=1","f");
    //legend->AddEntry(hRealJ1L0S1,"J=1,L=0,S=1","f");
    //legend->AddEntry(hRealJ1L1S1,"J=1,L=1,S=1","f");

    legend->Draw();


  }
  Double_t xl1=0.12, yl1 = 0.8, xl2=xl1+0.23, yl2=yl1+0.06;
  TLegend *legend = new TLegend(xl1,yl1,xl2,yl2);
  legend->SetFillColor(0);
  legend->SetLineColor(0);
  legend->SetTextSize(0.04);
 
  if (iShow == 2) {
    hRealJ0->Draw("e");
    legend->AddEntry(hRealJ0,"J=0","f");
    legend->Draw();
  }
  if (iShow == 3) {
    hRealJ1L1S0->Draw("e");
    legend->AddEntry(hRealJ1L1S0,"J=1,L=1,S=0","f");
    legend->Draw();
  }
  if (iShow == 4) {
    hRealJ1L1S0a0->Draw("e");
    legend->AddEntry(hRealJ1L1S0a0,"J=1,L=1,S=0 a0","f");
    legend->Draw();
  }
  if (iShow == 5) {
    hRealJ1L0S1->Draw("e");
    legend->AddEntry(hRealJ1L0S1,"J=1,L=0,S=1","f");
    legend->Draw();
  }
  if (iShow == 6) {
    hRealJ1L1S1->Draw("e");
    legend->AddEntry(hRealJ1L1S1,"J=1,L=1,S=1","f");
    legend->Draw();
  }
}