void pl_pla_resolution(UInt_t uTdcNumber1, UInt_t uChNumber1, UInt_t uTdcNumber2, UInt_t uChNumber2,
                       UInt_t uTdcNumber3, UInt_t uChNumber3, UInt_t uTdcNumber4, UInt_t uChNumber4)
{

  delete gROOT->FindObjectAny("tHistogram");

  FairRootManager* tManager = FairRootManager::Instance();
  if(!tManager)
  {
    cout<<"FairRootManager could not be retrieved. Abort macro execution."<<endl;
    return;
  }

  TTree* tOutputTree = tManager->GetOutTree();
  if(!tOutputTree)
  {
    cout<<"Output tree could not be retrieved from FairRootManager. Abort macro execution."<<endl;
    return;
  }

  TBranch* tBranch = tOutputTree->GetBranch("TofCalibData");
  if(!tBranch)
  {
    cout<<"Branch 'TofCalibData' not found in output tree. Abort macro execution."<<endl;
    return;
  }

  void* vOldAddress = tBranch->GetAddress();

  TDirectory* tOldDirectory = gDirectory;
  gROOT->cd();

  TClonesArray* tArray = new TClonesArray("TTofCalibData");
  tBranch->SetAddress(&tArray);

  TPaveStats* tStats = 0;
  TPaletteAxis* tPalette = 0;

  TCanvas* tCanvas1 = new TCanvas("tCanvas1","",0,0,500,500);
  tCanvas1->SetFillColor(0);
  tCanvas1->SetGridx(0);
  tCanvas1->SetGridy(0);
  tCanvas1->SetTopMargin(0.05);
  tCanvas1->SetRightMargin(0.05);
  tCanvas1->SetBottomMargin(0.15);
  tCanvas1->SetLeftMargin(0.15);

  TH1D* tHistogram = new TH1D("tHistogram","",2000,-10000,10000);
  tHistogram->GetXaxis()->SetTitle(Form("dt PLA (%u,%u;%u,%u) - (%u,%u;%u,%u) [ps]",uTdcNumber1,uChNumber1,uTdcNumber2,uChNumber2,uTdcNumber3,uChNumber3,uTdcNumber4,uChNumber4));
  tHistogram->GetXaxis()->SetTitleSize(0.06);
  tHistogram->GetXaxis()->SetTitleOffset(1.20);
  tHistogram->GetXaxis()->SetLabelSize(0.055);
  tHistogram->GetXaxis()->SetLabelOffset(0.018);
  tHistogram->GetYaxis()->SetTitleSize(0.06);
  tHistogram->GetYaxis()->SetTitleOffset(1.27);
  tHistogram->GetYaxis()->SetLabelSize(0.055);
  tHistogram->GetYaxis()->SetLabelOffset(0.01);
  tHistogram->SetLineColor(kBlue);
  tHistogram->SetLineWidth(2);
  tHistogram->GetYaxis()->SetTitle("");
  tHistogram->GetXaxis()->SetNdivisions(509);

  Long64_t lBranchEntries = tBranch->GetEntries();

  for(Long64_t lBranchEntry = 0; lBranchEntry < lBranchEntries; lBranchEntry++)
  {
    tArray->Clear("C");

    tBranch->GetEntry(lBranchEntry);

    Int_t iDataIndex1 = -1;
    Int_t iDataIndex2 = -1;
    Int_t iDataIndex3 = -1;
    Int_t iDataIndex4 = -1;
    
    Int_t iArrayEntries = tArray->GetEntriesFast();

    for(Int_t iArrayEntry = 0; iArrayEntry < iArrayEntries; iArrayEntry++)
    {
      TTofCalibData* tCalibTdcData = (TTofCalibData*)tArray->At(iArrayEntry);

      if(tCalibTdcData->GetBoard() == uTdcNumber1 && tCalibTdcData->GetChannel() == uChNumber1)
      {
        iDataIndex1 = iArrayEntry;
      }
      else if(tCalibTdcData->GetBoard() == uTdcNumber2 && tCalibTdcData->GetChannel() == uChNumber2)
      {
        iDataIndex2 = iArrayEntry;
      }

      else if(tCalibTdcData->GetBoard() == uTdcNumber3 && tCalibTdcData->GetChannel() == uChNumber3)
      {
        iDataIndex3 = iArrayEntry;
      }

      else if(tCalibTdcData->GetBoard() == uTdcNumber4 && tCalibTdcData->GetChannel() == uChNumber4)
      {
        iDataIndex4 = iArrayEntry;
      }


    }

    if(iDataIndex1 != -1 && iDataIndex2 != -1 && iDataIndex3 != -1 && iDataIndex4 != -1)
    {
      tHistogram->Fill( 0.5*(((TTofCalibData*)tArray->At(iDataIndex1))->GetTime()+((TTofCalibData*)tArray->At(iDataIndex2))->GetTime())
                       -0.5*(((TTofCalibData*)tArray->At(iDataIndex3))->GetTime()+((TTofCalibData*)tArray->At(iDataIndex4))->GetTime()) );
    }

  }

  tBranch->SetAddress(vOldAddress);

  tArray->Clear("C");
  tArray = NULL;

  (tHistogram->GetYaxis())->SetRangeUser(0.2,tHistogram->GetMaximum()+1000.);
  tHistogram->Draw();
  gPad->Update();
  tHistogram->Fit("gaus","Q");
  (tHistogram->GetFunction("gaus"))->FixParameter(0,(tHistogram->GetFunction("gaus"))->GetParameter(0));
  (tHistogram->GetFunction("gaus"))->FixParameter(1,(tHistogram->GetFunction("gaus"))->GetParameter(1));
  (tHistogram->GetFunction("gaus"))->SetLineWidth(3);
  gPad->SetLogy();
  tStats = (TPaveStats*)tHistogram->FindObject("stats");
  tStats->SetOptStat(1010);
  tStats->SetOptFit(11);
  tStats->SetX1NDC(0.50);
  tStats->SetY1NDC(0.70);
  tStats->SetX2NDC(0.95);
  tStats->SetY2NDC(0.95);
  tStats->SetTextSizePixels(25);
  tStats = 0;
  tPalette = 0;
  gPad->Modified();
  gPad->Update();

  tOldDirectory->cd();
}