void polPlot(int iShow){ 
  gROOT->Reset();//Reset session memory
  gROOT->SetStyle("Plain");//Make style plain
  gStyle->SetOptStat(0);//Remove the statistics box from histogram

  //Make a canvas
  TCanvas * myCanvas = new TCanvas("myCanvas","Me canvas title");

  //For this example, I will plot 5 points
  const Int_t nPoints = 4;
  //Create vectors that have the x and y coordinates 
  Double_t xVec[nPoints] = {0,45,90,135};
  Double_t yVec[nPoints] = {0.346,0.3476,0.3498,0.3746};

  //Create vectors that have the x and y errors
  Double_t xErrVec[nPoints] = {0,0,0,0}; //Common for x-error = 0
  Double_t yErrVec[nPoints] = {0.0137,0.0108,0.0106,0.0100};
  //0 degrees:     0.3460 +/- 0.0137
  //45 degrees:     0.3476 +/- 0.0108
  //90 degrees:     -0.3498 +/- 0.0106
  //135 degrees:     -0.3746 +/- 0.0100


  //full kill
  //0 degrees:     0.3739 +/- 0.0158
  //45 degrees:    0.3419 +/- 0.0115
  //90 degrees:   -0.3531 +/- 0.0127
  //135 degrees:  -0.3668 +/- 0.0106
  Double_t yVec3[nPoints] = {0.3739,0.3419,0.3531,0.3668};
  Double_t yErrVec3[nPoints] = {0.0158,0.0115,0.0127,0.0106};

  //full kill
  //SigRho0 = 0.882885 +/- 0.0373469
  //SigRho45 = 0.940842 +/- 0.0315752
  //SigRho90 = -0.911846 +/- 0.0328524
  //SigRho135 = -0.917904 +/- 0.0265056
  Double_t yVec31[nPoints] = {0.882885,0.940842,0.911846,0.917904};
  Double_t yErrVec31[nPoints] = {0.0373469,0.0315752,0.0328524,0.0265056};


  //Sim fit full kill
  //pol0   = 0.361605 +/- 0.0146987
  //pol45  = 0.341629 +/- 0.0114203
  //pol90  = 0.349121 +/- 0.0118854
  //pol135 = 0.366529 +/- 0.010547
  Double_t yVec4[nPoints] = {0.361605,0.341629,0.349121,0.366529};
  Double_t yErrVec4[nPoints] = {0.0146987,0.0114203,0.0118854,0.010547};

  //Sim fit full kill
  //SigRho0 = 0.912874 +/- 0.0371069
  //SigRho45 = 0.941665 +/- 0.0314788
  //SigRho90 = 0.922316 +/- 0.0313992
  //SigRho135 = 0.918617 +/- 0.0264336
  Double_t yVec41[nPoints] = {0.912874,0.941665,0.922316,0.918617};
  Double_t yErrVec41[nPoints] = {0.0371069,0.0314788,0.0313992,0.0264336};

  //Create the TGraphErrors object
  TGraphErrors * myGraph = new TGraphErrors(nPoints,xVec,yVec,xErrVec,yErrVec);
  myGraph->SetTitle("");
  myGraph->SetMarkerStyle(21);
  myGraph->SetMarkerSize(1.0);
  myGraph->SetMarkerColor(1);
  myGraph->SetLineColor(1);
  myGraph->SetLineWidth(2);

  TGraphErrors * myGraph3 = new TGraphErrors(nPoints,xVec,yVec3,xErrVec,yErrVec3);
  myGraph3->SetTitle("");
  myGraph3->SetMarkerStyle(21);
  myGraph3->SetMarkerSize(1.0);
  myGraph3->SetMarkerColor(2);
  myGraph3->SetLineColor(2);
  myGraph3->SetLineWidth(2);

  TGraphErrors * myGraph31 = new TGraphErrors(nPoints,xVec,yVec31,xErrVec,yErrVec31);
  myGraph31->SetTitle("");
  myGraph31->SetMarkerStyle(21);
  myGraph31->SetMarkerSize(1.0);
  myGraph31->SetMarkerColor(2);
  myGraph31->SetLineColor(2);
  myGraph31->SetLineWidth(2);

  TGraphErrors * myGraph4 = new TGraphErrors(nPoints,xVec,yVec4,xErrVec,yErrVec4);
  myGraph4->SetTitle("");
  myGraph4->SetMarkerStyle(21);
  myGraph4->SetMarkerSize(1.0);
  myGraph4->SetMarkerColor(4);
  myGraph4->SetLineColor(4);
  myGraph4->SetLineWidth(2);

  TGraphErrors * myGraph41 = new TGraphErrors(nPoints,xVec,yVec41,xErrVec,yErrVec41);
  myGraph41->SetTitle("");
  myGraph41->SetMarkerStyle(21);
  myGraph41->SetMarkerSize(1.0);
  myGraph41->SetMarkerColor(4);
  myGraph41->SetLineColor(4);
  myGraph41->SetLineWidth(2);

  //Draw the graph
  TH1D *h1 = new TH1D("h1","",100,-10,150);
  h1->GetXaxis()->SetTitle("Orientation (degree)");
  h1->GetYaxis()->SetTitle("Polarization");
  h1->GetYaxis()->SetTitleOffset(1.3);


  TH1D *h11 = new TH1D("h11","",100,-10,150);
  h11->GetXaxis()->SetTitle("Orientation (degree)");
  h11->GetYaxis()->SetTitle("#Sigma(#rho)");
  h11->GetYaxis()->SetTitleOffset(1.3);

  if (iShow == 1) {
    h1->SetMaximum(0.39);
    h1->SetMinimum(0.32);
    h1->Draw();
    myGraph->Draw("P");//A->Draw axis, L->Polyline, P->Draw marker
    
    myGraph->GetXaxis()->SetTitle("Orientation (degree)");
    myGraph->GetYaxis()->SetTitle("Polarization");
    myGraph->GetYaxis()->SetTitleOffset(1.3);
    //  myGraph->SetMinimum(0);
  }

  if (iShow == 3) {
    h1->SetMaximum(0.39);
    h1->SetMinimum(0.32);
    h1->Draw();
    myGraph3->Draw("P");//A->Draw axis, L->Polyline, P->Draw marker   
    myGraph3->GetXaxis()->SetTitle("Orientation (degree)");
    myGraph3->GetYaxis()->SetTitle("Polarization");
    myGraph3->GetYaxis()->SetTitleOffset(1.3);

    myGraph3->Fit("pol0");
    //myGraph3->Draw("P");

    //myGraph4->Draw("P");
    //  myGraph->SetMinimum(0);
  }

  if (iShow == 4) {
    h1->SetMaximum(0.39);
    h1->SetMinimum(0.32);
    h1->Draw();
    myGraph4->Draw("P");//A->Draw axis, L->Polyline, P->Draw marker   
    myGraph4->GetXaxis()->SetTitle("Orientation (degree)");
    myGraph4->GetYaxis()->SetTitle("Polarization");
    myGraph4->GetYaxis()->SetTitleOffset(1.3);
    myGraph4->Fit("pol0");
  }

  if (iShow == 41) {
    h11->SetMaximum(1.0);
    h11->SetMinimum(0.8);
    h11->Draw();
    myGraph41->Draw("P");//A->Draw axis, L->Polyline, P->Draw marker   
    myGraph41->GetXaxis()->SetTitle("Orientation (degree)");
    //myGraph41->GetYaxis()->SetTitle("Polarization");
    myGraph41->GetYaxis()->SetTitleOffset(1.3);
    myGraph41->Fit("pol0");
  }

  if (iShow == 31) {
    h11->SetMaximum(1.0);
    h11->SetMinimum(0.8);
    h11->Draw();
    myGraph31->Draw("P");//A->Draw axis, L->Polyline, P->Draw marker   
    myGraph31->GetXaxis()->SetTitle("Orientation (degree)");
    //myGraph41->GetYaxis()->SetTitle("Polarization");
    myGraph31->GetYaxis()->SetTitleOffset(1.3);
    myGraph31->Fit("pol0");
  }
  return;
}