void pl_over_CS(Double_t dChi2Lim, Double_t dSigmaDX04, Double_t dSigmaDY04, Double_t dSigmaDT04)
{
//  gStyle->SetOptStat(0);
  gStyle->SetOptStat("emr");
//  gStyle->SetOptFit(0);
  gStyle->SetLabelSize(0.07, "x");
  gStyle->SetLabelSize(0.07, "y");
  gStyle->SetLabelSize(0.06, "z");

  Bool_t gbChi2LimFitRange = kFALSE;


  TCanvas* tCanvas = new TCanvas("tCanvas", "", 0., 0., 900., 1200.);
  tCanvas->Divide(3, 4);

  TH1* tHisto;
  TH1* tHisto1;
  TH2* tHisto2;
  TProfile* tProfile;
  TLine* tLine;
  TPaveStats* tStatsBox;

  if(gbChi2LimFitRange)
  {
    TFile::Open("./calib_ana.cor_out.root", "READ");
    if(gFile)
    {
      tHisto = dynamic_cast<TH1*>(gFile->FindObjectAny("hSelTypeNNResidualX_Width"));
      dSigmaDX04 = tHisto->GetBinContent(1);
      tHisto = dynamic_cast<TH1*>(gFile->FindObjectAny("hSelTypeNNResidualY_Width"));
      dSigmaDY04 = tHisto->GetBinContent(1);
      tHisto = dynamic_cast<TH1*>(gFile->FindObjectAny("hSelTypeNNResidualT_Width"));
      dSigmaDT04 = tHisto->GetBinContent(1);

      gFile->Close();
      delete gFile;
    }
    else
    {
      cout<<"Could not open output calibration file. Using default 04 sigmas."<<endl;
    }
  }

  Double_t dFitRangeDX04 = (gbChi2LimFitRange ? TMath::Sqrt(dChi2Lim) : 1.)*dSigmaDX04;
  Double_t dFitRangeDY04 = (gbChi2LimFitRange ? TMath::Sqrt(dChi2Lim) : 1.)*dSigmaDY04;
  Double_t dFitRangeDT04 = (gbChi2LimFitRange ? TMath::Sqrt(dChi2Lim) : 1.)*dSigmaDT04;


  TH1F* tAllDX04SigmaCluSize = new TH1F("tAllDX04SigmaCluSize",
                                        "; cluster size [strips]; global #sigma(#Deltax) [cm]",
                                        20, 0.5, 20.5);

  TH1F* tAllDY04SigmaCluSize = new TH1F("tAllDY04SigmaCluSize",
                                        "; cluster size [strips]; global #sigma(#Deltay) [cm]",
                                        20, 0.5, 20.5);

  TH1F* tAllDT04SigmaCluSize = new TH1F("tAllDT04SigmaCluSize",
                                        "; cluster size [strips]; global #sigma(#Deltat) [cm]",
                                        20, 0.5, 20.5);


  for(Int_t iClusterSize = 1; iClusterSize <= 20; iClusterSize++)
  {
    tHisto2 = dynamic_cast<TH2*>(gROOT->FindObjectAny("hAllDX04CluSize"));
    tHisto1 = tHisto2->ProjectionY("_py", iClusterSize, iClusterSize);

    if(tHisto1->GetEntries() >= 1000.)
    {
      if(!static_cast<Int_t>(tHisto1->Fit("gaus", "Q0", "", -dFitRangeDX04, dFitRangeDX04)))
      {
        tAllDX04SigmaCluSize->SetBinContent(iClusterSize, (tHisto1->GetFunction("gaus"))->GetParameter("Sigma"));
        tAllDX04SigmaCluSize->SetBinError(iClusterSize, (tHisto1->GetFunction("gaus"))->GetParError(2));
      }
    }

    delete tHisto1;


    tHisto2 = dynamic_cast<TH2*>(gROOT->FindObjectAny("hAllDY04CluSize"));
    tHisto1 = tHisto2->ProjectionY("_py", iClusterSize, iClusterSize);

    if(tHisto1->GetEntries() >= 1000.)
    {
      if(!static_cast<Int_t>(tHisto1->Fit("gaus", "Q0", "", -dFitRangeDY04, dFitRangeDY04)))
      {
        tAllDY04SigmaCluSize->SetBinContent(iClusterSize, (tHisto1->GetFunction("gaus"))->GetParameter("Sigma"));
        tAllDY04SigmaCluSize->SetBinError(iClusterSize, (tHisto1->GetFunction("gaus"))->GetParError(2));
      }
    }

    delete tHisto1;


    tHisto2 = dynamic_cast<TH2*>(gROOT->FindObjectAny("hAllDT04CluSize"));
    tHisto1 = tHisto2->ProjectionY("_py", iClusterSize, iClusterSize);

    if(tHisto1->GetEntries() >= 1000.)
    {
      if(!static_cast<Int_t>(tHisto1->Fit("gaus", "Q0", "", -dFitRangeDT04, dFitRangeDT04)))
      {
        tAllDT04SigmaCluSize->SetBinContent(iClusterSize, (tHisto1->GetFunction("gaus"))->GetParameter("Sigma"));
        tAllDT04SigmaCluSize->SetBinError(iClusterSize, (tHisto1->GetFunction("gaus"))->GetParError(2));
      }
    }

    delete tHisto1;
  }


  tCanvas->cd(1);
  gPad->SetLogz();
  tHisto2 = dynamic_cast<TH2*>(gROOT->FindObjectAny("hAllDX04CluSize"));
  tHisto2->GetYaxis()->SetRangeUser(-20., 20.);
  tHisto2->Draw("colz");
  tProfile = tHisto2->ProfileX("_pfx", 1, -1, "s");
  tProfile->Draw("SAME");
  gPad->Modified();
  gPad->Update();
  tStatsBox = dynamic_cast<TPaveStats*>(tHisto2->FindObject("stats"));
  tStatsBox->SetX1NDC(0.500);
  tStatsBox->SetY1NDC(0.62);
  tStatsBox->SetX2NDC(0.843);
  tStatsBox->SetY2NDC(0.90);
  gPad->Modified();
  gPad->Update();

  tCanvas->cd(2);
  gPad->SetLogz();
  tHisto2 = dynamic_cast<TH2*>(gROOT->FindObjectAny("hAllDY04CluSize"));
  tHisto2->GetYaxis()->SetRangeUser(-20., 20.);
  tHisto2->Draw("colz");
  tProfile = tHisto2->ProfileX("_pfx", 1, -1, "s");
  tProfile->Draw("SAME");
  gPad->Modified();
  gPad->Update();
  tStatsBox = dynamic_cast<TPaveStats*>(tHisto2->FindObject("stats"));
  tStatsBox->SetX1NDC(0.500);
  tStatsBox->SetY1NDC(0.62);
  tStatsBox->SetX2NDC(0.843);
  tStatsBox->SetY2NDC(0.90);
  gPad->Modified();
  gPad->Update();

  tCanvas->cd(3);
  gPad->SetLogz();
  tHisto2 = dynamic_cast<TH2*>(gROOT->FindObjectAny("hAllDT04CluSize"));
  tHisto2->GetYaxis()->SetRangeUser(-3., 3.);
  tHisto2->Draw("colz");
  tProfile = tHisto2->ProfileX("_pfx", 1, -1, "s");
  tProfile->Draw("SAME");
  gPad->Modified();
  gPad->Update();
  tStatsBox = dynamic_cast<TPaveStats*>(tHisto2->FindObject("stats"));
  tStatsBox->SetX1NDC(0.500);
  tStatsBox->SetY1NDC(0.62);
  tStatsBox->SetX2NDC(0.843);
  tStatsBox->SetY2NDC(0.90);
  gPad->Modified();
  gPad->Update();

  tCanvas->cd(4);
  tAllDX04SigmaCluSize->GetYaxis()->SetRangeUser(0., 3.0);
  tAllDX04SigmaCluSize->SetStats(0);
  tAllDX04SigmaCluSize->Draw();
  gPad->Modified();
  gPad->Update();

  tCanvas->cd(5);
  tAllDY04SigmaCluSize->GetYaxis()->SetRangeUser(0., 3.0);
  tAllDY04SigmaCluSize->SetStats(0);
  tAllDY04SigmaCluSize->Draw();
  gPad->Modified();
  gPad->Update();

  tCanvas->cd(6);
  tAllDT04SigmaCluSize->GetYaxis()->SetRangeUser(0., 0.15);
  tAllDT04SigmaCluSize->SetStats(0);
  tAllDT04SigmaCluSize->Draw();
  gPad->Modified();
  gPad->Update();

  tCanvas->cd(7);
  gPad->SetLogy();
  tHisto = dynamic_cast<TH1*>(gROOT->FindObjectAny("hSelHitTupleSel2SelMul"));
  tHisto->SetLineColor(kRed);
  tHisto->Draw();
  gPad->Modified();
  gPad->Update();
  tStatsBox = dynamic_cast<TPaveStats*>(tHisto->FindObject("stats"));
  tStatsBox->SetX1NDC(0.50);
  tStatsBox->SetY1NDC(0.71);
  tStatsBox->SetX2NDC(0.98);
  tStatsBox->SetY2NDC(0.90);
  gPad->Modified();
  gPad->Update();

  tCanvas->cd(8);
  gPad->SetLogy();
  tHisto = dynamic_cast<TH1*>(gROOT->FindObjectAny("hSelHitTupleDutSelMul"));
  tHisto->SetLineColor(kRed);
  tHisto->Draw();
  gPad->Modified();
  gPad->Update();
  tStatsBox = dynamic_cast<TPaveStats*>(tHisto->FindObject("stats"));
  tStatsBox->SetX1NDC(0.50);
  tStatsBox->SetY1NDC(0.71);
  tStatsBox->SetX2NDC(0.98);
  tStatsBox->SetY2NDC(0.90);
  gPad->Modified();
  gPad->Update();

  tCanvas->cd(9);
  gPad->SetLogy();
  tHisto = dynamic_cast<TH1*>(gROOT->FindObjectAny("hSelHitTupleMRefSelMul"));
  tHisto->SetLineColor(kRed);
  tHisto->Draw();
  gPad->Modified();
  gPad->Update();
  tStatsBox = dynamic_cast<TPaveStats*>(tHisto->FindObject("stats"));
  tStatsBox->SetX1NDC(0.50);
  tStatsBox->SetY1NDC(0.71);
  tStatsBox->SetX2NDC(0.98);
  tStatsBox->SetY2NDC(0.90);
  gPad->Modified();
  gPad->Update();

  tCanvas->cd(10);
  gPad->SetLogy();
  tHisto = dynamic_cast<TH1*>(gROOT->FindObjectAny("hSel2CluSize"));
  tHisto->Draw();
  gPad->Modified();
  gPad->Update();
  tStatsBox = dynamic_cast<TPaveStats*>(tHisto->FindObject("stats"));
  tStatsBox->SetX1NDC(0.50);
  tStatsBox->SetY1NDC(0.75);
  tStatsBox->SetX2NDC(0.98);
  tStatsBox->SetY2NDC(0.94);
  tHisto1 = dynamic_cast<TH1*>(gROOT->FindObjectAny("hSelHitTupleSel2SelCluSize"));
  tHisto1->SetLineColor(kRed);
  tHisto1->Draw("SAMES");
  gPad->Modified();
  gPad->Update();
  tStatsBox = dynamic_cast<TPaveStats*>(tHisto1->FindObject("stats"));
  tStatsBox->SetX1NDC(0.50);
  tStatsBox->SetY1NDC(0.51);
  tStatsBox->SetX2NDC(0.98);
  tStatsBox->SetY2NDC(0.71);
  gPad->Modified();
  gPad->Update();

  tCanvas->cd(11);
  gPad->SetLogy();
  tHisto = dynamic_cast<TH1*>(gROOT->FindObjectAny("hDUTCluSize"));
  tHisto->Draw();
  gPad->Modified();
  gPad->Update();
  tStatsBox = dynamic_cast<TPaveStats*>(tHisto->FindObject("stats"));
  tStatsBox->SetX1NDC(0.50);
  tStatsBox->SetY1NDC(0.75);
  tStatsBox->SetX2NDC(0.98);
  tStatsBox->SetY2NDC(0.94);
  tHisto2 = dynamic_cast<TH2*>(gROOT->FindObjectAny("hAllSelHitTupleDutCluSizeTIS"));
  tHisto1 = tHisto2->ProjectionY();
  tHisto1->SetLineColor(kRed);
  tHisto1->Draw("SAMES");
  gPad->Modified();
  gPad->Update();
  tStatsBox = dynamic_cast<TPaveStats*>(tHisto1->FindObject("stats"));
  tStatsBox->SetX1NDC(0.50);
  tStatsBox->SetY1NDC(0.51);
  tStatsBox->SetX2NDC(0.98);
  tStatsBox->SetY2NDC(0.71);
  gPad->Modified();
  gPad->Update();

  tCanvas->cd(12);
  gPad->SetLogy();
  tHisto = dynamic_cast<TH1*>(gROOT->FindObjectAny("hMRefCluSize"));
  tHisto->Draw();
  gPad->Modified();
  gPad->Update();
  tStatsBox = dynamic_cast<TPaveStats*>(tHisto->FindObject("stats"));
  tStatsBox->SetX1NDC(0.50);
  tStatsBox->SetY1NDC(0.75);
  tStatsBox->SetX2NDC(0.98);
  tStatsBox->SetY2NDC(0.94);
  tHisto1 = dynamic_cast<TH1*>(gROOT->FindObjectAny("hSelHitTupleMRefSelCluSize"));
  tHisto1->SetLineColor(kRed);
  tHisto1->Draw("SAMES");
  gPad->Modified();
  gPad->Update();
  tStatsBox = dynamic_cast<TPaveStats*>(tHisto1->FindObject("stats"));
  tStatsBox->SetX1NDC(0.50);
  tStatsBox->SetY1NDC(0.51);
  tStatsBox->SetX2NDC(0.98);
  tStatsBox->SetY2NDC(0.71);
  gPad->Modified();
  gPad->Update();


  tCanvas->SaveAs("pl_over_CS.pdf");
}