void makePlot(int iShow){

  gROOT->Reset();
  gROOT->SetStyle("Plain");
  gStyle->SetOptStat(0);
  char  hId[120];
  char  fString[120];
  char  tString[120];


  //TFile *hFileData = new TFile("./outS3Data_All.root");
  //TFile *hFileThrown = new TFile("./outS3Thrown_All.root");
  //TFile *hFileSeen = new TFile("./outS3Seen_All.root");

  TFile *hFileData = new TFile("./outS3Data_m1285.root");
  TFile *hFileThrown = new TFile("./outS3Thrown_m1285.root");
  TFile *hFileSeen = new TFile("./outS3Seen_m1285.root");

  TH2D *hThetaPhiGJD; hFileData->GetObject("hThetaPhiGJ",hThetaPhiGJD);
  TH2D *hThetaPhiGJT; hFileThrown->GetObject("hThetaPhiGJ",hThetaPhiGJT);
  TH2D *hThetaPhiGJS; hFileSeen->GetObject("hThetaPhiGJ",hThetaPhiGJS);

  TH2D *hThetaPhiGJ2D; hFileData->GetObject("hThetaPhiGJ2",hThetaPhiGJ2D);
  TH2D *hThetaPhiGJ2T; hFileThrown->GetObject("hThetaPhiGJ2",hThetaPhiGJ2T);
  TH2D *hThetaPhiGJ2S; hFileSeen->GetObject("hThetaPhiGJ2",hThetaPhiGJ2S);

  TH2D *hCosThetaPhiGJD; hFileData->GetObject("hCosThetaPhiGJ",hCosThetaPhiGJD);
  TH2D *hCosThetaPhiGJT; hFileThrown->GetObject("hCosThetaPhiGJ",hCosThetaPhiGJT);
  TH2D *hCosThetaPhiGJS; hFileSeen->GetObject("hCosThetaPhiGJ",hCosThetaPhiGJS);

  int rX = 8;
  int rY = 50;
  //hThetaPhiGJD->Rebin2D(rX,rY);
  //hThetaPhiGJT->Rebin2D(rX,rY);
  //hThetaPhiGJS->Rebin2D(rX,rY);

  hThetaPhiGJ2D->Rebin2D(rX,rY);
  hThetaPhiGJ2T->Rebin2D(rX,rY);
  hThetaPhiGJ2S->Rebin2D(rX,rY);

  hCosThetaPhiGJD->Rebin2D(rX,rY);
  hCosThetaPhiGJT->Rebin2D(rX,rY);
  hCosThetaPhiGJS->Rebin2D(rX,rY);

  TH2D * hEffThetaGJ = (TH2D *)hThetaPhiGJS->Clone("hEffThetaGJ");
  TH2D * hEffTheta2GJ = (TH2D *)hThetaPhiGJ2S->Clone("hEffTheta2GJ");
  TH2D * hEffCosThetaGJ = (TH2D *)hCosThetaPhiGJS->Clone("hEffCosThetaGJ");

  TH2D * hCorThetaGJ = (TH2D *)hThetaPhiGJD->Clone("hCorThetaGJ");
  TH2D * hCorTheta2GJ = (TH2D *)hThetaPhiGJ2D->Clone("hCorTheta2GJ");
  TH2D * hCorCosThetaGJ = (TH2D *)hCosThetaPhiGJD->Clone("hCorCosThetaGJ");

  hEffThetaGJ->Divide(hThetaPhiGJT);
  hEffTheta2GJ->Divide(hThetaPhiGJ2T);
  hEffCosThetaGJ->Divide(hCosThetaPhiGJT);

  hCorThetaGJ->Divide(hEffThetaGJ);
  hCorTheta2GJ->Divide(hEffTheta2GJ);
  hCorCosThetaGJ->Divide(hEffCosThetaGJ);

  ///////////////

  TH2D *hThetaPhiD; hFileData->GetObject("hThetaPhi",hThetaPhiD);
  TH2D *hThetaPhiT; hFileThrown->GetObject("hThetaPhi",hThetaPhiT);
  TH2D *hThetaPhiS; hFileSeen->GetObject("hThetaPhi",hThetaPhiS);

  TH2D *hThetaPhi2D; hFileData->GetObject("hThetaPhi2",hThetaPhi2D);
  TH2D *hThetaPhi2T; hFileThrown->GetObject("hThetaPhi2",hThetaPhi2T);
  TH2D *hThetaPhi2S; hFileSeen->GetObject("hThetaPhi2",hThetaPhi2S);

  TH2D *hCosThetaPhiD; hFileData->GetObject("hCosThetaPhi",hCosThetaPhiD);
  TH2D *hCosThetaPhiT; hFileThrown->GetObject("hCosThetaPhi",hCosThetaPhiT);
  TH2D *hCosThetaPhiS; hFileSeen->GetObject("hCosThetaPhi",hCosThetaPhiS);

  hThetaPhiD->Rebin2D(rX,rY);
  hThetaPhiT->Rebin2D(rX,rY);
  hThetaPhiS->Rebin2D(rX,rY);

  hThetaPhi2D->Rebin2D(rX,rY);
  hThetaPhi2T->Rebin2D(rX,rY);
  hThetaPhi2S->Rebin2D(rX,rY);

  hCosThetaPhiD->Rebin2D(rX,rY);
  hCosThetaPhiT->Rebin2D(rX,rY);
  hCosThetaPhiS->Rebin2D(rX,rY);

  TH2D * hEffTheta = (TH2D *)hThetaPhiS->Clone("hEffTheta");
  TH2D * hEffTheta2 = (TH2D *)hThetaPhi2S->Clone("hEffTheta2");
  TH2D * hEffCosTheta = (TH2D *)hCosThetaPhiS->Clone("hEffCosTheta");

  TH2D * hCorTheta = (TH2D *)hThetaPhiD->Clone("hCorTheta");
  TH2D * hCorTheta2 = (TH2D *)hThetaPhi2D->Clone("hCorTheta2");
  TH2D * hCorCosTheta = (TH2D *)hCosThetaPhiD->Clone("hCorCosTheta");

  hEffTheta->Sumw2();
  hEffTheta2->Sumw2();
  hEffCosTheta->Sumw2();

  hCorTheta->Sumw2();
  hCorTheta2->Sumw2();
  hCorCosTheta->Sumw2();

  hEffTheta->Divide(hThetaPhiT);
  hEffTheta2->Divide(hThetaPhi2T);
  hEffCosTheta->Divide(hCosThetaPhiT);

  hCorTheta->Divide(hEffTheta);
  hCorTheta2->Divide(hEffTheta2);
  hCorCosTheta->Divide(hEffCosTheta);

  ///////////////

  if (iShow == 1){
    hCorThetaGJ->Draw("pol colz");
  } 
  if (iShow == 2){
    //hCorTheta2GJ->SetMaximum(700);
    hCorTheta2GJ->Draw("colz");
  } 
  if (iShow == 3){
    hCorCosThetaGJ->Draw("colz");
  } 

  /////////////

  if (iShow == 4){
    hCorTheta->Draw("pol colz");
  } 
  if (iShow == 5){
    //hCorTheta2->Rebin2D(2,2);
    //hCorTheta2->SetMaximum(700);
    hCorTheta2->Draw("colz");
  } 
  if (iShow == 6){
    //hCorCosTheta->Rebin2D(2,2);
    //hCorCosTheta->SetMaximum(700);
    hCorCosTheta->Draw("colz");
  } 


}