void qa(Int_t runNumber) {
  TStopwatch timer;
  timer.Start();

  gEnv->SetValue("XSec.GSI.DelegProxy","2");
  // Select ROOT version
  TProof::Mgr("aliprod@alicecaf")->SetROOTVersion("v5-24-00b-caf");
  // Login to CAF
  TProof::Open("aliprod@alicecaf");

  // Enable AliRoot
  gProof->UploadPackage("/afs/cern.ch/alice/caf/sw/ALICE/PARs/v4-17-Release.rec/AF-v4-17-rec.par");
  gProof->EnablePackage("AF-v4-17-rec.par");

  gProof->Load(Form("%s/AliAnalysisTaskCosmic.cxx++g",
		    gSystem->ExpandPathName("$ALICE_ROOT/PWG1/cosmic"));

  //____________________________________________//
  // Make the analysis manager
  AliAnalysisManager *mgr = new AliAnalysisManager("AnalysisQAManager");
  AliESDInputHandler* esdH = new AliESDInputHandler();
  mgr->SetInputEventHandler(esdH);  
  mgr->SetDebugLevel(10);

  //____________________________________________//
  // 1st Cosmic task
  AliAnalysisTaskCosmic *task1 = new AliAnalysisTaskCosmic("TaskCosmic");
  mgr->AddTask(task1);

  // Create containers for input/output
  AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("chist1",TList::Class(),AliAnalysisManager::kOutputContainer,
							    Form("run%d.root",runNumber));

  //____________________________________________//
  mgr->ConnectInput  (task1,  0, mgr->GetCommonInputContainer());
  mgr->ConnectOutput (task1,  1, coutput1);
  
  if (!mgr->InitAnalysis()) return;
  mgr->PrintStatus();
  mgr->StartAnalysis("proof",
		     Form("ALIREC/aliprod/run%d",runNumber));

  timer.Stop();
  timer.Print();

  plot(runNumber);
}

void plot(Int_t runNumber)
{
  TFile* f = new TFile(Form("run%d.root",runNumber), "read");

  // pt, phi ....
  TCanvas* c1 = new TCanvas("c1", "pt,eta,phi", 10, 10, 1100, 800);
  c1->Divide(2,2);
  c1->cd(1);
  c1->GetPad(1)->SetLogy();
    
  fhPtP = (TH1F*) chist1->FindObject("fhPtPC");
  fhPtP->SetLineColor(2);
  fhPtP->Draw();
  fhPtN = (TH1F*) chist1->FindObject("fhPtNC");
  fhPtN->SetLineColor(4);
  fhPtN->Draw("same");

  c1->cd(2);
  fhPhiP = (TH1F*) chist1->FindObject("fhPhiPC");
  fhPhiP->SetLineColor(2);
  fhPhiP->Draw();
  fhPhiN = (TH1F*) chist1->FindObject("fhPhiNC");
  fhPhiN->SetLineColor(4);
  fhPhiN->Draw("same");

  c1->cd(3);
  fhThetaP = (TH1F*) chist1->FindObject("fhThetaPC");
  fhThetaP->SetLineColor(2);
  fhThetaP->Draw();
  fhThetaN = (TH1F*) chist1->FindObject("fhThetaNC");
  fhThetaN->SetLineColor(4);
  fhThetaN->Draw("same");

    
  // delta phi
  TCanvas* c2 = new TCanvas("c2", "#Delta#phi", 10, 10, 1100, 800);
  c2->Divide(2,2);
  // pos/neg charge
  c2->cd(1);
  c2->GetPad(1)->SetLogy();
  fhDPhiP = (TH1F*) chist1->FindObject("fhDPhiPC");
  fhDPhiN = (TH1F*) chist1->FindObject("fhDPhiNC");
  fhDPhiP->SetLineColor(2);
  fhDPhiN->SetLineColor(4);
  fhDPhiP->SetXTitle("#Delta#phi [rad]");
  fhDPhiP->Draw();
  fhDPhiN->Draw("same");

  // pos/neg z
  c2->cd(2);
  c2->GetPad(2)->SetLogy();

  fhDPhiZP = (TH1F*) chist1->FindObject("fhDPhiPZ");
  fhDPhiZN = (TH1F*) chist1->FindObject("fhDPhiNZ");
  fhDPhiBA = (TH1F*) chist1->FindObject("fhDPhi_Bad");
  fhDPhiZP->SetLineColor(2);
  fhDPhiZN->SetLineColor(4);
  fhDPhiBA->SetLineColor(6);

  fhDPhiZP->SetXTitle("#Delta#phi [rad]");
  fhDPhiZP->Draw();
  fhDPhiZN->Draw("same");
  //  fhDPhiBA->Draw("same");

    
  // delta theta
  // pos/neg charge
  c2->cd(3);
  c2->GetPad(3)->SetLogy();
  fhDThetaP = (TH1F*) chist1->FindObject("fhDThetaPC");
  fhDThetaN = (TH1F*) chist1->FindObject("fhDThetaNC");
  fhDThetaP->SetLineColor(2);
  fhDThetaN->SetLineColor(4);
 

  fhDThetaP->SetXTitle("#Delta#theta [rad]");
  fhDThetaP->Draw();
  fhDThetaN->Draw("same");

  // pos/neg z
  c2->cd(4);
  c2->GetPad(4)->SetLogy();

  fhDThetaZP = (TH1F*) chist1->FindObject("fhDThetaPZ");
  fhDThetaZN = (TH1F*) chist1->FindObject("fhDThetaNZ");
  fhDThetaBA = (TH1F*) chist1->FindObject("fhDTheta_Bad");
  fhDThetaZP->SetLineColor(2);
  fhDThetaZN->SetLineColor(4);
  fhDThetaBA->SetLineColor(6);
  fhDThetaZP->SetXTitle("#Delta#theta [rad]");
  fhDThetaZP->Draw();

  fhDThetaZN->Draw("same");
  //    fhDThetaBA->Draw("same");

    
  // delta Pt
  TCanvas* c3 = new TCanvas("c3", "#Delta p_{T}", 10, 10, 1100, 800);
  c3->Divide(2,2);
  // pos/neg charge
  c3->cd(1);
  c3->GetPad(1)->SetLogy();
  fhDPtP = (TH1F*) chist1->FindObject("fhDPtPC");
  fhDPtN = (TH1F*) chist1->FindObject("fhDPtNC");
  fhDPtP->SetLineColor(2);
  fhDPtN->SetLineColor(4);
    
  fhDPtP->Draw();
  fhDPtN->Draw("same");

  // pos/neg z
  c3->cd(2);
  c3->GetPad(2)->SetLogy();

  fhDPtZP = (TH1F*) chist1->FindObject("fhDPtPZ");
  fhDPtZN = (TH1F*) chist1->FindObject("fhDPtNZ");
  fhDPtZP->SetLineColor(2);
  fhDPtZN->SetLineColor(4);
  fhDPtZP->Draw();
  fhDPtZN->Draw("same");

  c3->cd(3);

  fpDPtP = (TH1F*) chist1->FindObject("fpDPtPC");
  fpDPtN = (TH1F*) chist1->FindObject("fpDPtNC");
  fpDPtP->SetLineColor(2);
  fpDPtN->SetLineColor(4);
  fpDPtP->Draw();
  fpDPtN->Draw("same");

  c3->cd(4);

  fpDPtZP = (TH1F*) chist1->FindObject("fpDPtPZ");
  fpDPtZN = (TH1F*) chist1->FindObject("fpDPtNZ");
  fpDPtZP->SetLineColor(2);
  fpDPtZN->SetLineColor(4);
  fpDPtZP->Draw();
  fpDPtZN->Draw("same");
    

  // dZ
  TCanvas* c4 = new TCanvas("c4", "#Delta Z", 10, 10, 1100, 800);
  c4->Divide(2,2);
  // pos/neg charge
  c4->cd(1);
  c4->GetPad(1)->SetLogy();
  fhDZP = (TH1F*) chist1->FindObject("fhDZPC");
  fhDZN = (TH1F*) chist1->FindObject("fhDZNC");
  fhDZP->SetLineColor(2);
  fhDZN->SetLineColor(4);
    
  fhDZP->Draw();
  fhDZN->Draw("same");
  
  // pos/neg z
  c4->cd(2);
  c4->GetPad(2)->SetLogy();

  fhDZZP = (TH1F*) chist1->FindObject("fhDZPZ");
  fhDZZN = (TH1F*) chist1->FindObject("fhDZNZ");
  fhDZBA = (TH1F*) chist1->FindObject("fhDZ_Bad");

  fhDZZP->SetLineColor(2);
  fhDZZN->SetLineColor(4);
  fhDZBA->SetLineColor(6);
  
  fhDZZP->Draw();
  fhDZZN->Draw("same");
  fhDZBA->Draw("same");

  // dX
  // pos/neg charge
  c4->cd(3);
  c4->GetPad(3)->SetLogy();
  fhDXP = (TH1F*) chist1->FindObject("fhDXPC");
  fhDXN = (TH1F*) chist1->FindObject("fhDXNC");
  fhDXP->SetLineColor(2);
  fhDXN->SetLineColor(4);
    
  fhDXP->Draw();
  fhDXN->Draw("same");

  // pos/neg z
  c4->cd(4);
  c4->GetPad(4)->SetLogy();

  fhDXZP = (TH1F*) chist1->FindObject("fhDXPZ");
  fhDXZN = (TH1F*) chist1->FindObject("fhDXNZ");
  fhDXBA = (TH1F*) chist1->FindObject("fhDX_Bad");

  fhDXZP->SetLineColor(2);
  fhDXZN->SetLineColor(4);
  fhDXBA->SetLineColor(6);

  fhDXZP->Draw();
  fhDXZN->Draw("same");
  fhDXBA->Draw("same");

  // dY
  TCanvas* c4a = new TCanvas("c4a", "#Delta Y", 10, 10, 1100, 800);
  c4a->Divide(2,2);
  // pos/neg charge
  c4a->cd(1);
  c4a->GetPad(1)->SetLogy();
  fhDYP = (TH1F*) chist1->FindObject("fhDYPC");
  fhDYN = (TH1F*) chist1->FindObject("fhDYNC");
  fhDYP->SetLineColor(2);
  fhDYN->SetLineColor(4);
  
  fhDYP->Draw();
  fhDYN->Draw("same");

  // pos/neg z
  c4a->cd(2);
  c4a->GetPad(2)->SetLogy();

  fhDYZP = (TH1F*) chist1->FindObject("fhDYPZ");
  fhDYZN = (TH1F*) chist1->FindObject("fhDYNZ");
  fhDYBA = (TH1F*) chist1->FindObject("fhDY_Bad");
  
  fhDYZP->SetLineColor(2);
  fhDYZN->SetLineColor(4);
  fhDYBA->SetLineColor(6);

  fhDYZP->Draw();
  fhDYZN->Draw("same");
  fhDYBA->Draw("same");

  // delta Pt
  TCanvas* c5 = new TCanvas("c5", "#Delta p_{T} (n-sigma)", 10, 10, 1100, 800);
  c5->Divide(2,2);
  c5->cd(1);

  fpDPtSP = (TH1F*) chist1->FindObject("fpDPtSPC");
  fpDPtSN = (TH1F*) chist1->FindObject("fpDPtSNC");
  fpDPtSP->SetLineColor(2);
  fpDPtSN->SetLineColor(4);
  fpDPtSP->Draw();
  fpDPtSN->Draw("same");


  c5->cd(2);

  fpDPtSPZ = (TH1F*) chist1->FindObject("fpDPtSPZ");
  fpDPtSNZ = (TH1F*) chist1->FindObject("fpDPtSNZ");
  fpDPtSPZ->SetLineColor(2);
  fpDPtSNZ->SetLineColor(4);
  fpDPtSPZ->Draw();
  fpDPtSNZ->Draw("same");

  TCanvas* c6 = new TCanvas("c6", "Dz vs z", 10, 10, 1100, 800);
  fhDZvsZ = (TH2F*) chist1->FindObject("fhDZvsZ");
  fhDZvsZ->SetXTitle("z_{in} * sign(z_{in}) * sign(z_{out}) [cm]");
  fhDZvsZ->SetYTitle("#DeltaZ [cm]");

  gStyle->SetPalette(1);
    
  fhDZvsZ->Draw("colz");
}