void Pulls() {

  TStopwatch timer;
  timer.Start();
  
    TStyle *plain = new TStyle("Plain","Plain Style(no colors/fill areas)");
  plain->SetPadColor(0);
  plain->SetCanvasColor(0);
  plain->SetTitleColor(0);
  plain->SetStatColor(0);
  plain->SetOptStat(1110);
  plain->SetOptFit(0011);
  plain->SetStatW(0.225);
  plain->SetStatH(0.06);
  plain->SetOptTitle(0);
  plain->cd();
  
  TFile *f = new TFile("out.root","read");
/// ///////////////////////////////////////////////////////////////////////////
///                                       ETA                                //
/// ///////////////////////////////////////////////////////////////////////////
  TH1F *resol_eta = (TH1F*) f->Get("resolution_eta");
  TH1F *res_eta = (TH1F*) f->Get("residual_eta");
  TH1F *pull_eta =(TH1F*) f->Get("pull_eta");
  
  TCanvas *canvas_eta;
  canvas_eta= new TCanvas("vot", "gg");
  canvas_eta -> Divide(3,1);
  
  canvas_eta->cd(1);
  TF1 *fRes_eta = new TF1("eta-res","gaus");
  fRes_eta->SetLineColor(2);
  res_eta->Fit("eta-res","","",res_eta->GetXaxis()->GetXmin(),res_eta->GetXaxis()->GetXmax());
  res_eta->GetXaxis()->SetTitle("#eta^{reco} - #eta^{mc}");
  res_eta->GetXaxis()->SetTitleOffset(1);
  res_eta->Draw();
  fRes_eta->Draw("same");

  canvas_eta->cd(2);
  TF1 *fPull_eta = new TF1("eta-pull","gaus");
  fPull_eta->SetLineColor(2);
  pull_eta->Fit("eta-pull","","",pull_eta->GetXaxis()->GetXmin(),pull_eta->GetXaxis()->GetXmax());
  pull_eta->GetXaxis()->SetTitle("[#eta^{reco} - #eta^{mc}]/#sigma_{eta}");
  pull_eta->GetXaxis()->SetTitleOffset(1);
  pull_eta->Draw();
  fPull_eta->Draw("same");
  
  canvas_eta->cd(3);
  TF1 *fResol_eta = new TF1("eta-resol","gaus");
  fResol_eta->SetLineColor(2);
  resol_eta->Fit("eta-resol","","",resol_eta->GetXaxis()->GetXmin(),resol_eta->GetXaxis()->GetXmax());
  resol_eta->GetXaxis()->SetTitle("[#eta^{reco} - #eta^{mc}]/#eta^{mc}*100%");
  resol_eta->GetXaxis()->SetTitleOffset(1);
  resol_eta->Draw();
  fResol_eta->Draw("same");
  
  canvas_eta->SaveAs("eta.pdf");

/// ///////////////////////////////////////////////////////////////////////////
///                                       IP Prarmeter                       //
/// ///////////////////////////////////////////////////////////////////////////

  TH1F *res_ip_abs = (TH1F*) f->Get("res_ip_abs");
  TH1F *res_ip     = (TH1F*) f->Get("res_ip");
  
  TCanvas *canvas_ip;
  canvas_ip= new TCanvas("ip", "ip");
  canvas_ip -> Divide(2,1);
  
  canvas_ip->cd(1);
  TF1 *fRes_ip = new TF1("eta-ip","gaus");
  fRes_ip->SetLineColor(2);
  res_ip->Fit("eta-ip","","",res_ip->GetXaxis()->GetXmin(),res_ip->GetXaxis()->GetXmax());
  res_ip->GetXaxis()->SetTitle("IP^{ESD} - IP^{KFParticle}");
  res_ip->GetXaxis()->SetTitleOffset(1);
  res_ip->Draw();
  fRes_ip->Draw("same");

  canvas_ip->cd(2);
  TF1 *fRes_ip_abs = new TF1("eta-ip_abs","gaus");
  fRes_ip_abs->SetLineColor(2);
  res_ip_abs->Fit("eta-ip_abs","","",res_ip_abs->GetXaxis()->GetXmin(),res_ip_abs->GetXaxis()->GetXmax());
  res_ip_abs->GetXaxis()->SetTitle("|IP^{ESD}| - |IP^{KFParticle}|");
  res_ip_abs->GetXaxis()->SetTitleOffset(1);
  res_ip_abs->Draw();
  fRes_ip_abs->Draw("same");

/*  canvas_ip->cd(3);
  TF1 *fRes_ip_P = new TF1("eta-ip_P","gaus");
  fRes_ip_P->SetLineColor(2);
  res_ip_P->Fit("eta-ip_P","","",res_ip_P->GetXaxis()->GetXmin(),res_ip_P->GetXaxis()->GetXmax());
  res_ip_P->GetXaxis()->SetTitle("P_{diff} [GeV/c]");
  res_ip_P->GetXaxis()->SetTitleOffset(1);
  res_ip_P->Draw();
  fRes_ip_P->Draw("same");

  canvas_ip->cd(4);
  TF1 *fRes_ip_P_all = new TF1("eta-ip_P_all","gaus");
  fRes_ip_P_all->SetLineColor(2);
  res_ip_P_all->Fit("eta-ip_P_all","","",res_ip_P_all->GetXaxis()->GetXmin(),res_ip_P_all->GetXaxis()->GetXmax());
  res_ip_P_all->GetXaxis()->SetTitle("P_{all} [GeV/c]");
  res_ip_P_all->GetXaxis()->SetTitleOffset(1);
  res_ip_P_all->Draw();
  fRes_ip_P_all->Draw("same");*/

  canvas_ip->SaveAs("ip.pdf");

  TCanvas *canvas_ip_err;
  canvas_ip_err= new TCanvas("ip_err", "ip_err");
  canvas_ip_err -> Divide(2,1);
  
  plain->SetOptTitle(1);
  plain->SetStatY(0.9);

  canvas_ip_err->cd(1);
  TF1 *fRes_ip_err = new TF1("eta-ip_err","gaus");
  fRes_ip_err->SetLineColor(2);
  delta_ip->Fit("eta-ip_err","","",delta_ip->GetXaxis()->GetXmin(),delta_ip->GetXaxis()->GetXmax());
  delta_ip->GetXaxis()->SetTitle("#DeltaIP^{KFParticle}/IP^{KFParticle}*100%");
  delta_ip->GetXaxis()->SetTitleOffset(1);
delta_ip->SetTitle("without the vertex cov. matrix (setted to 0)");
  delta_ip->Draw();
  fRes_ip_err->Draw("same");

  canvas_ip_err->cd(2);
  TF1 *fRes_ip_err_abs = new TF1("eta-ip_err_abs","gaus");
  fRes_ip_err_abs->SetLineColor(2);
  delta_ip_cv->Fit("eta-ip_err_abs","","",delta_ip_cv->GetXaxis()->GetXmin(),delta_ip_cv->GetXaxis()->GetXmax());
  delta_ip_cv->GetXaxis()->SetTitle("#DeltaIP^{KFParticle}/IP^{KFParticle}*100%");
delta_ip_cv->SetTitle("with the vertex cov. matrix");
  delta_ip_cv->GetXaxis()->SetTitleOffset(1);
  delta_ip_cv->Draw();
  fRes_ip_err_abs->Draw("same");

  canvas_ip_err->SaveAs("ip_err.pdf");

}