ROOT logo
//
// Jet analysis task.
//
// Author: S.Aiola

#include <TClonesArray.h>
#include <TH2F.h>
#include <TH3F.h>
#include <THnSparse.h>
#include <TList.h>
#include <TLorentzVector.h>

#include "AliVCluster.h"
#include "AliVParticle.h"
#include "AliEmcalJet.h"
#include "AliRhoParameter.h"
#include "AliLog.h"
#include "AliJetContainer.h"

#include "AliAnalysisTaskSAJF.h"

ClassImp(AliAnalysisTaskSAJF)

//________________________________________________________________________
AliAnalysisTaskSAJF::AliAnalysisTaskSAJF() : 
  AliAnalysisTaskEmcalJet("AliAnalysisTaskSAJF", kTRUE),
  fHistoType(1),
  fHistRejectionReason(0),
  fHistTracksJetPt(0),
  fHistClustersJetPt(0),
  fHistTracksPtDist(0),
  fHistClustersPtDist(0),
  fHistJetObservables(0),
  fHistJetPtEtaPhi(0),
  fHistJetPtArea(0),
  fHistJetPtEP(0),
  fHistJetPtNEF(0),
  fHistJetPtZ(0),
  fHistJetPtLeadingPartPt(0),
  fHistJetCorrPtEtaPhi(0),
  fHistJetCorrPtArea(0),
  fHistJetCorrPtEP(0),
  fHistJetCorrPtNEF(0),
  fHistJetCorrPtZ(0),
  fHistJetCorrPtLeadingPartPt(0),
  fHistJetPtCorrPt(0),
  fHistJetPtMCPt(0),
  fHistJetMCPtCorrPt(0)
{
  // Default constructor.

  SetMakeGeneralHistograms(kTRUE);
}

//________________________________________________________________________
AliAnalysisTaskSAJF::AliAnalysisTaskSAJF(const char *name) : 
  AliAnalysisTaskEmcalJet(name, kTRUE),
  fHistoType(1),
  fHistRejectionReason(0),
  fHistTracksJetPt(0),
  fHistClustersJetPt(0),
  fHistTracksPtDist(0),
  fHistClustersPtDist(0),
  fHistJetObservables(0),
  fHistJetPtEtaPhi(0),
  fHistJetPtArea(0),
  fHistJetPtEP(0),
  fHistJetPtNEF(0),
  fHistJetPtZ(0),
  fHistJetPtLeadingPartPt(0),
  fHistJetCorrPtEtaPhi(0),
  fHistJetCorrPtArea(0),
  fHistJetCorrPtEP(0),
  fHistJetCorrPtNEF(0),
  fHistJetCorrPtZ(0),
  fHistJetCorrPtLeadingPartPt(0),
  fHistJetPtCorrPt(0),
  fHistJetPtMCPt(0),
  fHistJetMCPtCorrPt(0)
{
  // Standard constructor.

  SetMakeGeneralHistograms(kTRUE);
}

//________________________________________________________________________
void AliAnalysisTaskSAJF::AllocateTHnSparse()
{
  TString title[20]= {""};
  Int_t nbins[20]  = {0};
  Double_t min[20] = {0.};
  Double_t max[20] = {0.};
  Int_t dim = 0;

  if (fForceBeamType != kpp) {
    title[dim] = "Centrality (%)";
    nbins[dim] = 22;
    min[dim] = -5;
    max[dim] = 105;
    dim++;
    
    title[dim] = "#phi_{jet} - #psi_{RP}";
    nbins[dim] = 100;
    min[dim] = 0;
    max[dim] = TMath::Pi();
    dim++;
  }

  title[dim] = "#eta_{jet}";
  nbins[dim] = 100;
  min[dim] = -1;
  max[dim] = 1;
  dim++;

  title[dim] = "#phi_{jet} (rad)";
  nbins[dim] = 201;
  min[dim] = 0;
  max[dim] = 2*TMath::Pi()*nbins[dim]/(nbins[dim]-1);
  dim++;

  title[dim] = "p_{T} (GeV/c)";
  nbins[dim] = fNbins;
  min[dim] = fMinBinPt;
  max[dim] = fMaxBinPt;
  dim++;

  if (fIsEmbedded) {
    title[dim] = "p_{T}^{MC} (GeV/c)";
    nbins[dim] = fNbins;
    min[dim] = fMinBinPt;
    max[dim] = fMaxBinPt;
    dim++;
  }

  if (!GetRhoName().IsNull()) {
    title[dim] = "p_{T}^{corr} (GeV/c)";
    nbins[dim] = fNbins*2;
    min[dim] = -fMaxBinPt;
    max[dim] = fMaxBinPt;
    dim++;
  }

  title[dim] = "A_{jet}";
  nbins[dim] = 150;
  min[dim] = 0;
  max[dim] = 1.5;
  dim++;

  title[dim] = "NEF";
  nbins[dim] = 102;
  min[dim] = 0;
  max[dim] = 1.02;
  dim++;

  title[dim] = "Z";
  nbins[dim] = 102;
  min[dim] = 0;
  max[dim] = 1.02;
  dim++;

  title[dim] = "No. of constituents";
  nbins[dim] = 250;
  min[dim] = -0.5;
  max[dim] = 249.5;
  dim++;

  title[dim] = "p_{T,particle}^{leading} (GeV/c)";
  nbins[dim] = 120;
  min[dim] = 0;
  max[dim] = 120;
  dim++;

  fHistJetObservables = new THnSparseD("fHistJetObservables","fHistJetObservables",dim,nbins,min,max);
  fOutput->Add(fHistJetObservables);
  for (Int_t i = 0; i < dim; i++)
    fHistJetObservables->GetAxis(i)->SetTitle(title[i]);
}

//________________________________________________________________________
void AliAnalysisTaskSAJF::AllocateTHX()
{
  fHistJetPtEtaPhi = new TH3*[fNcentBins];
  fHistJetPtArea = new TH2*[fNcentBins];
  fHistJetPtEP = new TH2*[fNcentBins];
  fHistJetPtNEF = new TH2*[fNcentBins];
  fHistJetPtZ = new TH2*[fNcentBins];
  fHistJetPtLeadingPartPt = new TH2*[fNcentBins];
  fHistJetCorrPtEtaPhi = new TH3*[fNcentBins];
  fHistJetCorrPtArea = new TH2*[fNcentBins];
  fHistJetCorrPtEP = new TH2*[fNcentBins];
  fHistJetCorrPtNEF = new TH2*[fNcentBins];
  fHistJetCorrPtZ = new TH2*[fNcentBins];
  fHistJetCorrPtLeadingPartPt = new TH2*[fNcentBins];
  fHistJetPtCorrPt = new TH2*[fNcentBins];
  fHistJetPtMCPt = new TH2*[fNcentBins];
  fHistJetMCPtCorrPt = new TH2*[fNcentBins];

  for (Int_t i = 0; i < fNcentBins; i++) {
    TString histname;

    histname = "fHistJetPtEtaPhi_";
    histname += i;
    fHistJetPtEtaPhi[i] = new TH3F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, 20, -1, 1, 41, 0, 2*TMath::Pi()*41/40);
    fHistJetPtEtaPhi[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
    fHistJetPtEtaPhi[i]->GetYaxis()->SetTitle("#eta");
    fHistJetPtEtaPhi[i]->GetZaxis()->SetTitle("#phi_{jet} (rad)");
    fOutput->Add(fHistJetPtEtaPhi[i]);
      
    histname = "fHistJetPtArea_";
    histname += i;
    fHistJetPtArea[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, 150, 0, 1.5);
    fHistJetPtArea[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
    fHistJetPtArea[i]->GetYaxis()->SetTitle("A_{jet}");
    fHistJetPtArea[i]->GetZaxis()->SetTitle("counts");
    fOutput->Add(fHistJetPtArea[i]);

    histname = "fHistJetPtEP_";
    histname += i;
    fHistJetPtEP[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, 100, 0, TMath::Pi());
    fHistJetPtEP[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
    fHistJetPtEP[i]->GetYaxis()->SetTitle("#phi_{jet} - #psi_{RP}");
    fHistJetPtEP[i]->GetZaxis()->SetTitle("counts");
    fOutput->Add(fHistJetPtEP[i]);

    histname = "fHistJetPtNEF_";
    histname += i;
    fHistJetPtNEF[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, 102, 0, 1.02);
    fHistJetPtNEF[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
    fHistJetPtNEF[i]->GetYaxis()->SetTitle("NEF");
    fHistJetPtNEF[i]->GetZaxis()->SetTitle("counts");
    fOutput->Add(fHistJetPtNEF[i]);

    histname = "fHistJetPtZ_";
    histname += i;
    fHistJetPtZ[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, 102, 0, 1.02);
    fHistJetPtZ[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
    fHistJetPtZ[i]->GetYaxis()->SetTitle("z");
    fHistJetPtZ[i]->GetZaxis()->SetTitle("counts");
    fOutput->Add(fHistJetPtZ[i]);

    histname = "fHistJetPtLeadingPartPt_";
    histname += i;
    fHistJetPtLeadingPartPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, 120, 0, 120);
    fHistJetPtLeadingPartPt[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
    fHistJetPtLeadingPartPt[i]->GetYaxis()->SetTitle("p_{T,particle}^{leading} (GeV/c)");
    fHistJetPtLeadingPartPt[i]->GetZaxis()->SetTitle("counts");
    fOutput->Add(fHistJetPtLeadingPartPt[i]);

    if (!GetRhoName().IsNull()) {
      histname = "fHistJetCorrPtEtaPhi_";
      histname += i;
      fHistJetCorrPtEtaPhi[i] = new TH3F(histname.Data(), histname.Data(), fNbins*2, -fMaxBinPt, fMaxBinPt, 20, -1, 1, 41, 0, 2*TMath::Pi()*201/200);
      fHistJetCorrPtEtaPhi[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
      fHistJetCorrPtEtaPhi[i]->GetYaxis()->SetTitle("#eta");
      fHistJetCorrPtEtaPhi[i]->GetZaxis()->SetTitle("#phi_{jet} (rad)");
      fOutput->Add(fHistJetCorrPtEtaPhi[i]);
      
      histname = "fHistJetCorrPtArea_";
      histname += i;
      fHistJetCorrPtArea[i] = new TH2F(histname.Data(), histname.Data(), fNbins*2, -fMaxBinPt, fMaxBinPt, 150, 0, 1.5);
      fHistJetCorrPtArea[i]->GetXaxis()->SetTitle("p_{T}^{corr} (GeV/c)");
      fHistJetCorrPtArea[i]->GetYaxis()->SetTitle("A_{jet}");
      fHistJetCorrPtArea[i]->GetZaxis()->SetTitle("counts");
      fOutput->Add(fHistJetCorrPtArea[i]);

      histname = "fHistJetCorrPtEP_";
      histname += i;
      fHistJetCorrPtEP[i] = new TH2F(histname.Data(), histname.Data(), fNbins*2, -fMaxBinPt, fMaxBinPt, 100, 0, TMath::Pi());
      fHistJetCorrPtEP[i]->GetXaxis()->SetTitle("p_{T}^{corr} (GeV/c)");
      fHistJetCorrPtEP[i]->GetYaxis()->SetTitle("#phi_{jet} - #psi_{RP}");
      fHistJetCorrPtEP[i]->GetZaxis()->SetTitle("counts");
      fOutput->Add(fHistJetCorrPtEP[i]);

      histname = "fHistJetCorrPtNEF_";
      histname += i;
      fHistJetCorrPtNEF[i] = new TH2F(histname.Data(), histname.Data(), fNbins*2, -fMaxBinPt, fMaxBinPt, 102, 0, 1.02);
      fHistJetCorrPtNEF[i]->GetXaxis()->SetTitle("p_{T}^{corr} (GeV/c)");
      fHistJetCorrPtNEF[i]->GetYaxis()->SetTitle("NEF");
      fHistJetCorrPtNEF[i]->GetZaxis()->SetTitle("counts");
      fOutput->Add(fHistJetCorrPtNEF[i]);

      histname = "fHistJetCorrPtZ_";
      histname += i;
      fHistJetCorrPtZ[i] = new TH2F(histname.Data(), histname.Data(), fNbins*2, -fMaxBinPt, fMaxBinPt, 102, 0, 1.02);
      fHistJetCorrPtZ[i]->GetXaxis()->SetTitle("p_{T}^{corr} (GeV/c)");
      fHistJetCorrPtZ[i]->GetYaxis()->SetTitle("z");
      fHistJetCorrPtZ[i]->GetZaxis()->SetTitle("counts");
      fOutput->Add(fHistJetCorrPtZ[i]);

      histname = "fHistJetCorrPtLeadingPartPt_";
      histname += i;
      fHistJetCorrPtLeadingPartPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins*2, -fMaxBinPt, fMaxBinPt, 120, 0, 120);
      fHistJetCorrPtLeadingPartPt[i]->GetXaxis()->SetTitle("p_{T}^{corr} (GeV/c)");
      fHistJetCorrPtLeadingPartPt[i]->GetYaxis()->SetTitle("p_{T,particle}^{leading} (GeV/c)");
      fHistJetCorrPtLeadingPartPt[i]->GetZaxis()->SetTitle("counts");
      fOutput->Add(fHistJetCorrPtLeadingPartPt[i]);

      histname = "fHistJetPtCorrPt_";
      histname += i;
      fHistJetPtCorrPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, fNbins*2, -fMaxBinPt, fMaxBinPt);
      fHistJetPtCorrPt[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
      fHistJetPtCorrPt[i]->GetYaxis()->SetTitle("p_{T}^{corr} (GeV/c)");
      fHistJetPtCorrPt[i]->GetZaxis()->SetTitle("counts");
      fOutput->Add(fHistJetPtCorrPt[i]);

      if (fIsEmbedded) {
	histname = "fHistJetMCPtCorrPt_";
	histname += i;
	fHistJetMCPtCorrPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, fNbins*2, -fMaxBinPt, fMaxBinPt);
	fHistJetMCPtCorrPt[i]->GetXaxis()->SetTitle("p_{T}^{MC} (GeV/c)");
	fHistJetMCPtCorrPt[i]->GetYaxis()->SetTitle("p_{T}^{corr} (GeV/c)");
	fHistJetMCPtCorrPt[i]->GetZaxis()->SetTitle("counts");
	fOutput->Add(fHistJetMCPtCorrPt[i]);
      }
    }

    if (fIsEmbedded) {
      histname = "fHistJetPtMCPt_";
      histname += i;
      fHistJetPtMCPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, fNbins, fMinBinPt, fMaxBinPt);
      fHistJetPtMCPt[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
      fHistJetPtMCPt[i]->GetYaxis()->SetTitle("p_{T}^{MC} (GeV/c)");
      fHistJetPtMCPt[i]->GetZaxis()->SetTitle("counts");
      fOutput->Add(fHistJetPtMCPt[i]);
    }
  }
}

//________________________________________________________________________
void AliAnalysisTaskSAJF::UserCreateOutputObjects()
{
  // Create user output.

  AliAnalysisTaskEmcalJet::UserCreateOutputObjects();

  if (fHistoType == 0) 
    AllocateTHX();
  else
    AllocateTHnSparse();

  fHistTracksJetPt = new TH2*[fNcentBins];
  fHistClustersJetPt = new TH2*[fNcentBins];
  fHistTracksPtDist = new TH2*[fNcentBins];
  fHistClustersPtDist = new TH2*[fNcentBins];
  fHistRejectionReason = new TH2*[fNcentBins];

  for (Int_t i = 0; i < fNcentBins; i++) {
    TString histname;

    if (fParticleCollArray.GetEntriesFast() > 0) {
      histname = "fHistTracksJetPt_";
      histname += i;
      fHistTracksJetPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins / 2, fMinBinPt, fMaxBinPt / 2, fNbins, fMinBinPt, fMaxBinPt);
      fHistTracksJetPt[i]->GetXaxis()->SetTitle("p_{T,track} (GeV/c)");
      fHistTracksJetPt[i]->GetYaxis()->SetTitle("p_{T,jet} (GeV/c)");
      fHistTracksJetPt[i]->GetZaxis()->SetTitle("counts");
      fOutput->Add(fHistTracksJetPt[i]);
      
      histname = "fHistTracksPtDist_";
      histname += i;
      fHistTracksPtDist[i] = new TH2F(histname.Data(), histname.Data(), fNbins / 2, fMinBinPt, fMaxBinPt / 2, 100, 0, 5);
      fHistTracksPtDist[i]->GetXaxis()->SetTitle("p_{T,track} (GeV/c)");
      fHistTracksPtDist[i]->GetYaxis()->SetTitle("d");
      fHistTracksPtDist[i]->GetZaxis()->SetTitle("counts");
      fOutput->Add(fHistTracksPtDist[i]);
    }

    if (fClusterCollArray.GetEntriesFast() > 0) {
      histname = "fHistClustersJetPt_";
      histname += i;
      fHistClustersJetPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins / 2, fMinBinPt, fMaxBinPt / 2, fNbins, fMinBinPt, fMaxBinPt);
      fHistClustersJetPt[i]->GetXaxis()->SetTitle("p_{T,clus} (GeV/c)");
      fHistClustersJetPt[i]->GetYaxis()->SetTitle("p_{T,jet} (GeV/c)");
      fHistClustersJetPt[i]->GetZaxis()->SetTitle("counts");
      fOutput->Add(fHistClustersJetPt[i]);

      histname = "fHistClustersPtDist_";
      histname += i;
      fHistClustersPtDist[i] = new TH2F(histname.Data(), histname.Data(), fNbins / 2, fMinBinPt, fMaxBinPt / 2, 100, 0, 5);
      fHistClustersPtDist[i]->GetXaxis()->SetTitle("p_{T,clus} (GeV/c)");
      fHistClustersPtDist[i]->GetYaxis()->SetTitle("d");
      fHistClustersPtDist[i]->GetZaxis()->SetTitle("counts");
      fOutput->Add(fHistClustersPtDist[i]);
    }

    histname = "fHistRejectionReason_";
    histname += i;
    fHistRejectionReason[i] = new TH2F(histname, histname, 32, 0, 32, 100, 0, 250);
    fHistRejectionReason[i]->GetXaxis()->SetTitle("Rejection reason");
    fHistRejectionReason[i]->GetYaxis()->SetTitle("p_{T,jet} (GeV/c)");
    fHistRejectionReason[i]->GetZaxis()->SetTitle("counts");
    SetRejectionReasonLabels(fHistRejectionReason[i]->GetXaxis());
    fOutput->Add(fHistRejectionReason[i]);
  }

  PostData(1, fOutput); // Post data for ALL output slots >0 here, to get at least an empty histogram

}

//________________________________________________________________________
Bool_t AliAnalysisTaskSAJF::FillHistograms()
{
  // Fill histograms.

  AliJetContainer *jets = static_cast<AliJetContainer*>(fJetCollArray.At(0));

  if (!jets) return kFALSE;
  
  AliEmcalJet* jet = 0;

  jets->ResetCurrentID();
  while ((jet = jets->GetNextJet())) {
    if (!jet) {
      AliError("Could not receive jet!");
      continue;
    }

    if (!jets->AcceptJet(jet)) {
      fHistRejectionReason[fCentBin]->Fill(jets->GetRejectionReasonBitPosition(), jet->Pt());
      continue;
    }

    Float_t ptLeading = GetLeadingHadronPt(jet);
    Float_t corrPt = jet->Pt() - fRhoVal * jet->Area();

    // Fill THnSparse
    Double_t ep = jet->Phi() - fEPV0;
    while (ep < 0) ep += TMath::Pi();
    while (ep >= TMath::Pi()) ep -= TMath::Pi();

    FillJetHisto(fCent, ep, jet->Eta(), jet->Phi(), jet->Pt(), jet->MCPt(), corrPt, jet->Area(), 
		 jet->NEF(), ptLeading/jet->Pt(), jet->GetNumberOfConstituents(), ptLeading);

    if (fTracks) {
      for (Int_t it = 0; it < jet->GetNumberOfTracks(); it++) {
	AliVParticle *track = jet->TrackAt(it, fTracks);
	if (track) {
	  fHistTracksJetPt[fCentBin]->Fill(track->Pt(), jet->Pt());
	  Double_t dist = TMath::Sqrt((track->Eta() - jet->Eta()) * (track->Eta() - jet->Eta()) + (track->Phi() - jet->Phi()) * (track->Phi() - jet->Phi()));
	  fHistTracksPtDist[fCentBin]->Fill(track->Pt(), dist);
	}
      }
    }

    if (fCaloClusters) {
      for (Int_t ic = 0; ic < jet->GetNumberOfClusters(); ic++) {
	AliVCluster *cluster = jet->ClusterAt(ic, fCaloClusters);
	
	if (cluster) {
	  TLorentzVector nPart;
	  cluster->GetMomentum(nPart, fVertex);

	  fHistClustersJetPt[fCentBin]->Fill(nPart.Pt(), jet->Pt());
	  Double_t dist = TMath::Sqrt((nPart.Eta() - jet->Eta()) * (nPart.Eta() - jet->Eta()) + (nPart.Phi() - jet->Phi()) * (nPart.Phi() - jet->Phi()));
	  fHistClustersPtDist[fCentBin]->Fill(nPart.Pt(), dist);
	}
      }
    }
  } //jet loop 

  return kTRUE;
}

//________________________________________________________________________
void AliAnalysisTaskSAJF::FillJetHisto(Double_t cent, Double_t ep, Double_t eta, Double_t phi, Double_t pt, Double_t MCpt, Double_t corrpt, Double_t area, 
				       Double_t NEF, Double_t z, Int_t n, Double_t leadingpt)
{
  if (fHistoType == 0) {
    fHistJetPtEtaPhi[fCentBin]->Fill(pt,eta,phi);
    fHistJetPtArea[fCentBin]->Fill(pt,area);
    fHistJetPtEP[fCentBin]->Fill(pt,ep);
    fHistJetPtNEF[fCentBin]->Fill(pt,NEF);
    fHistJetPtZ[fCentBin]->Fill(pt,z);
    fHistJetPtLeadingPartPt[fCentBin]->Fill(pt,leadingpt);
    if (fHistJetCorrPtEtaPhi[fCentBin]) {
      fHistJetCorrPtEtaPhi[fCentBin]->Fill(corrpt,eta,phi);
      fHistJetCorrPtArea[fCentBin]->Fill(corrpt,area);
      fHistJetCorrPtEP[fCentBin]->Fill(corrpt,ep);
      fHistJetCorrPtNEF[fCentBin]->Fill(corrpt,NEF);
      fHistJetCorrPtZ[fCentBin]->Fill(corrpt,z);
      fHistJetCorrPtLeadingPartPt[fCentBin]->Fill(corrpt,leadingpt);
      fHistJetPtCorrPt[fCentBin]->Fill(pt,corrpt);
      if (fIsEmbedded)
	fHistJetMCPtCorrPt[fCentBin]->Fill(MCpt,corrpt);
    }
    if (fIsEmbedded)
      fHistJetPtMCPt[fCentBin]->Fill(pt,MCpt);
  }
  else {
 
    Double_t contents[20]={0};

    for (Int_t i = 0; i < fHistJetObservables->GetNdimensions(); i++) {
      TString title(fHistJetObservables->GetAxis(i)->GetTitle());
      if (title=="Centrality (%)")
	contents[i] = cent;
      else if (title=="#phi_{jet} - #psi_{RP}")
	contents[i] = ep;
      else if (title=="#eta_{jet}")
	contents[i] = eta;
      else if (title=="#phi_{jet} (rad)")
	contents[i] = phi;
      else if (title=="p_{T} (GeV/c)")
	contents[i] = pt;
      else if (title=="p_{T}^{MC} (GeV/c)")
	contents[i] = MCpt;
      else if (title=="p_{T}^{corr} (GeV/c)")
	contents[i] = corrpt;
      else if (title=="A_{jet}")
	contents[i] = area;
      else if (title=="NEF")
	contents[i] = NEF;
      else if (title=="Z")
	contents[i] = z;
      else if (title=="No. of constituents")
	contents[i] = n;
      else if (title=="p_{T,particle}^{leading} (GeV/c)")
	contents[i] = leadingpt;
      else 
	AliWarning(Form("Unable to fill dimension %s!",title.Data()));
    }

    fHistJetObservables->Fill(contents);
  }
}
 AliAnalysisTaskSAJF.cxx:1
 AliAnalysisTaskSAJF.cxx:2
 AliAnalysisTaskSAJF.cxx:3
 AliAnalysisTaskSAJF.cxx:4
 AliAnalysisTaskSAJF.cxx:5
 AliAnalysisTaskSAJF.cxx:6
 AliAnalysisTaskSAJF.cxx:7
 AliAnalysisTaskSAJF.cxx:8
 AliAnalysisTaskSAJF.cxx:9
 AliAnalysisTaskSAJF.cxx:10
 AliAnalysisTaskSAJF.cxx:11
 AliAnalysisTaskSAJF.cxx:12
 AliAnalysisTaskSAJF.cxx:13
 AliAnalysisTaskSAJF.cxx:14
 AliAnalysisTaskSAJF.cxx:15
 AliAnalysisTaskSAJF.cxx:16
 AliAnalysisTaskSAJF.cxx:17
 AliAnalysisTaskSAJF.cxx:18
 AliAnalysisTaskSAJF.cxx:19
 AliAnalysisTaskSAJF.cxx:20
 AliAnalysisTaskSAJF.cxx:21
 AliAnalysisTaskSAJF.cxx:22
 AliAnalysisTaskSAJF.cxx:23
 AliAnalysisTaskSAJF.cxx:24
 AliAnalysisTaskSAJF.cxx:25
 AliAnalysisTaskSAJF.cxx:26
 AliAnalysisTaskSAJF.cxx:27
 AliAnalysisTaskSAJF.cxx:28
 AliAnalysisTaskSAJF.cxx:29
 AliAnalysisTaskSAJF.cxx:30
 AliAnalysisTaskSAJF.cxx:31
 AliAnalysisTaskSAJF.cxx:32
 AliAnalysisTaskSAJF.cxx:33
 AliAnalysisTaskSAJF.cxx:34
 AliAnalysisTaskSAJF.cxx:35
 AliAnalysisTaskSAJF.cxx:36
 AliAnalysisTaskSAJF.cxx:37
 AliAnalysisTaskSAJF.cxx:38
 AliAnalysisTaskSAJF.cxx:39
 AliAnalysisTaskSAJF.cxx:40
 AliAnalysisTaskSAJF.cxx:41
 AliAnalysisTaskSAJF.cxx:42
 AliAnalysisTaskSAJF.cxx:43
 AliAnalysisTaskSAJF.cxx:44
 AliAnalysisTaskSAJF.cxx:45
 AliAnalysisTaskSAJF.cxx:46
 AliAnalysisTaskSAJF.cxx:47
 AliAnalysisTaskSAJF.cxx:48
 AliAnalysisTaskSAJF.cxx:49
 AliAnalysisTaskSAJF.cxx:50
 AliAnalysisTaskSAJF.cxx:51
 AliAnalysisTaskSAJF.cxx:52
 AliAnalysisTaskSAJF.cxx:53
 AliAnalysisTaskSAJF.cxx:54
 AliAnalysisTaskSAJF.cxx:55
 AliAnalysisTaskSAJF.cxx:56
 AliAnalysisTaskSAJF.cxx:57
 AliAnalysisTaskSAJF.cxx:58
 AliAnalysisTaskSAJF.cxx:59
 AliAnalysisTaskSAJF.cxx:60
 AliAnalysisTaskSAJF.cxx:61
 AliAnalysisTaskSAJF.cxx:62
 AliAnalysisTaskSAJF.cxx:63
 AliAnalysisTaskSAJF.cxx:64
 AliAnalysisTaskSAJF.cxx:65
 AliAnalysisTaskSAJF.cxx:66
 AliAnalysisTaskSAJF.cxx:67
 AliAnalysisTaskSAJF.cxx:68
 AliAnalysisTaskSAJF.cxx:69
 AliAnalysisTaskSAJF.cxx:70
 AliAnalysisTaskSAJF.cxx:71
 AliAnalysisTaskSAJF.cxx:72
 AliAnalysisTaskSAJF.cxx:73
 AliAnalysisTaskSAJF.cxx:74
 AliAnalysisTaskSAJF.cxx:75
 AliAnalysisTaskSAJF.cxx:76
 AliAnalysisTaskSAJF.cxx:77
 AliAnalysisTaskSAJF.cxx:78
 AliAnalysisTaskSAJF.cxx:79
 AliAnalysisTaskSAJF.cxx:80
 AliAnalysisTaskSAJF.cxx:81
 AliAnalysisTaskSAJF.cxx:82
 AliAnalysisTaskSAJF.cxx:83
 AliAnalysisTaskSAJF.cxx:84
 AliAnalysisTaskSAJF.cxx:85
 AliAnalysisTaskSAJF.cxx:86
 AliAnalysisTaskSAJF.cxx:87
 AliAnalysisTaskSAJF.cxx:88
 AliAnalysisTaskSAJF.cxx:89
 AliAnalysisTaskSAJF.cxx:90
 AliAnalysisTaskSAJF.cxx:91
 AliAnalysisTaskSAJF.cxx:92
 AliAnalysisTaskSAJF.cxx:93
 AliAnalysisTaskSAJF.cxx:94
 AliAnalysisTaskSAJF.cxx:95
 AliAnalysisTaskSAJF.cxx:96
 AliAnalysisTaskSAJF.cxx:97
 AliAnalysisTaskSAJF.cxx:98
 AliAnalysisTaskSAJF.cxx:99
 AliAnalysisTaskSAJF.cxx:100
 AliAnalysisTaskSAJF.cxx:101
 AliAnalysisTaskSAJF.cxx:102
 AliAnalysisTaskSAJF.cxx:103
 AliAnalysisTaskSAJF.cxx:104
 AliAnalysisTaskSAJF.cxx:105
 AliAnalysisTaskSAJF.cxx:106
 AliAnalysisTaskSAJF.cxx:107
 AliAnalysisTaskSAJF.cxx:108
 AliAnalysisTaskSAJF.cxx:109
 AliAnalysisTaskSAJF.cxx:110
 AliAnalysisTaskSAJF.cxx:111
 AliAnalysisTaskSAJF.cxx:112
 AliAnalysisTaskSAJF.cxx:113
 AliAnalysisTaskSAJF.cxx:114
 AliAnalysisTaskSAJF.cxx:115
 AliAnalysisTaskSAJF.cxx:116
 AliAnalysisTaskSAJF.cxx:117
 AliAnalysisTaskSAJF.cxx:118
 AliAnalysisTaskSAJF.cxx:119
 AliAnalysisTaskSAJF.cxx:120
 AliAnalysisTaskSAJF.cxx:121
 AliAnalysisTaskSAJF.cxx:122
 AliAnalysisTaskSAJF.cxx:123
 AliAnalysisTaskSAJF.cxx:124
 AliAnalysisTaskSAJF.cxx:125
 AliAnalysisTaskSAJF.cxx:126
 AliAnalysisTaskSAJF.cxx:127
 AliAnalysisTaskSAJF.cxx:128
 AliAnalysisTaskSAJF.cxx:129
 AliAnalysisTaskSAJF.cxx:130
 AliAnalysisTaskSAJF.cxx:131
 AliAnalysisTaskSAJF.cxx:132
 AliAnalysisTaskSAJF.cxx:133
 AliAnalysisTaskSAJF.cxx:134
 AliAnalysisTaskSAJF.cxx:135
 AliAnalysisTaskSAJF.cxx:136
 AliAnalysisTaskSAJF.cxx:137
 AliAnalysisTaskSAJF.cxx:138
 AliAnalysisTaskSAJF.cxx:139
 AliAnalysisTaskSAJF.cxx:140
 AliAnalysisTaskSAJF.cxx:141
 AliAnalysisTaskSAJF.cxx:142
 AliAnalysisTaskSAJF.cxx:143
 AliAnalysisTaskSAJF.cxx:144
 AliAnalysisTaskSAJF.cxx:145
 AliAnalysisTaskSAJF.cxx:146
 AliAnalysisTaskSAJF.cxx:147
 AliAnalysisTaskSAJF.cxx:148
 AliAnalysisTaskSAJF.cxx:149
 AliAnalysisTaskSAJF.cxx:150
 AliAnalysisTaskSAJF.cxx:151
 AliAnalysisTaskSAJF.cxx:152
 AliAnalysisTaskSAJF.cxx:153
 AliAnalysisTaskSAJF.cxx:154
 AliAnalysisTaskSAJF.cxx:155
 AliAnalysisTaskSAJF.cxx:156
 AliAnalysisTaskSAJF.cxx:157
 AliAnalysisTaskSAJF.cxx:158
 AliAnalysisTaskSAJF.cxx:159
 AliAnalysisTaskSAJF.cxx:160
 AliAnalysisTaskSAJF.cxx:161
 AliAnalysisTaskSAJF.cxx:162
 AliAnalysisTaskSAJF.cxx:163
 AliAnalysisTaskSAJF.cxx:164
 AliAnalysisTaskSAJF.cxx:165
 AliAnalysisTaskSAJF.cxx:166
 AliAnalysisTaskSAJF.cxx:167
 AliAnalysisTaskSAJF.cxx:168
 AliAnalysisTaskSAJF.cxx:169
 AliAnalysisTaskSAJF.cxx:170
 AliAnalysisTaskSAJF.cxx:171
 AliAnalysisTaskSAJF.cxx:172
 AliAnalysisTaskSAJF.cxx:173
 AliAnalysisTaskSAJF.cxx:174
 AliAnalysisTaskSAJF.cxx:175
 AliAnalysisTaskSAJF.cxx:176
 AliAnalysisTaskSAJF.cxx:177
 AliAnalysisTaskSAJF.cxx:178
 AliAnalysisTaskSAJF.cxx:179
 AliAnalysisTaskSAJF.cxx:180
 AliAnalysisTaskSAJF.cxx:181
 AliAnalysisTaskSAJF.cxx:182
 AliAnalysisTaskSAJF.cxx:183
 AliAnalysisTaskSAJF.cxx:184
 AliAnalysisTaskSAJF.cxx:185
 AliAnalysisTaskSAJF.cxx:186
 AliAnalysisTaskSAJF.cxx:187
 AliAnalysisTaskSAJF.cxx:188
 AliAnalysisTaskSAJF.cxx:189
 AliAnalysisTaskSAJF.cxx:190
 AliAnalysisTaskSAJF.cxx:191
 AliAnalysisTaskSAJF.cxx:192
 AliAnalysisTaskSAJF.cxx:193
 AliAnalysisTaskSAJF.cxx:194
 AliAnalysisTaskSAJF.cxx:195
 AliAnalysisTaskSAJF.cxx:196
 AliAnalysisTaskSAJF.cxx:197
 AliAnalysisTaskSAJF.cxx:198
 AliAnalysisTaskSAJF.cxx:199
 AliAnalysisTaskSAJF.cxx:200
 AliAnalysisTaskSAJF.cxx:201
 AliAnalysisTaskSAJF.cxx:202
 AliAnalysisTaskSAJF.cxx:203
 AliAnalysisTaskSAJF.cxx:204
 AliAnalysisTaskSAJF.cxx:205
 AliAnalysisTaskSAJF.cxx:206
 AliAnalysisTaskSAJF.cxx:207
 AliAnalysisTaskSAJF.cxx:208
 AliAnalysisTaskSAJF.cxx:209
 AliAnalysisTaskSAJF.cxx:210
 AliAnalysisTaskSAJF.cxx:211
 AliAnalysisTaskSAJF.cxx:212
 AliAnalysisTaskSAJF.cxx:213
 AliAnalysisTaskSAJF.cxx:214
 AliAnalysisTaskSAJF.cxx:215
 AliAnalysisTaskSAJF.cxx:216
 AliAnalysisTaskSAJF.cxx:217
 AliAnalysisTaskSAJF.cxx:218
 AliAnalysisTaskSAJF.cxx:219
 AliAnalysisTaskSAJF.cxx:220
 AliAnalysisTaskSAJF.cxx:221
 AliAnalysisTaskSAJF.cxx:222
 AliAnalysisTaskSAJF.cxx:223
 AliAnalysisTaskSAJF.cxx:224
 AliAnalysisTaskSAJF.cxx:225
 AliAnalysisTaskSAJF.cxx:226
 AliAnalysisTaskSAJF.cxx:227
 AliAnalysisTaskSAJF.cxx:228
 AliAnalysisTaskSAJF.cxx:229
 AliAnalysisTaskSAJF.cxx:230
 AliAnalysisTaskSAJF.cxx:231
 AliAnalysisTaskSAJF.cxx:232
 AliAnalysisTaskSAJF.cxx:233
 AliAnalysisTaskSAJF.cxx:234
 AliAnalysisTaskSAJF.cxx:235
 AliAnalysisTaskSAJF.cxx:236
 AliAnalysisTaskSAJF.cxx:237
 AliAnalysisTaskSAJF.cxx:238
 AliAnalysisTaskSAJF.cxx:239
 AliAnalysisTaskSAJF.cxx:240
 AliAnalysisTaskSAJF.cxx:241
 AliAnalysisTaskSAJF.cxx:242
 AliAnalysisTaskSAJF.cxx:243
 AliAnalysisTaskSAJF.cxx:244
 AliAnalysisTaskSAJF.cxx:245
 AliAnalysisTaskSAJF.cxx:246
 AliAnalysisTaskSAJF.cxx:247
 AliAnalysisTaskSAJF.cxx:248
 AliAnalysisTaskSAJF.cxx:249
 AliAnalysisTaskSAJF.cxx:250
 AliAnalysisTaskSAJF.cxx:251
 AliAnalysisTaskSAJF.cxx:252
 AliAnalysisTaskSAJF.cxx:253
 AliAnalysisTaskSAJF.cxx:254
 AliAnalysisTaskSAJF.cxx:255
 AliAnalysisTaskSAJF.cxx:256
 AliAnalysisTaskSAJF.cxx:257
 AliAnalysisTaskSAJF.cxx:258
 AliAnalysisTaskSAJF.cxx:259
 AliAnalysisTaskSAJF.cxx:260
 AliAnalysisTaskSAJF.cxx:261
 AliAnalysisTaskSAJF.cxx:262
 AliAnalysisTaskSAJF.cxx:263
 AliAnalysisTaskSAJF.cxx:264
 AliAnalysisTaskSAJF.cxx:265
 AliAnalysisTaskSAJF.cxx:266
 AliAnalysisTaskSAJF.cxx:267
 AliAnalysisTaskSAJF.cxx:268
 AliAnalysisTaskSAJF.cxx:269
 AliAnalysisTaskSAJF.cxx:270
 AliAnalysisTaskSAJF.cxx:271
 AliAnalysisTaskSAJF.cxx:272
 AliAnalysisTaskSAJF.cxx:273
 AliAnalysisTaskSAJF.cxx:274
 AliAnalysisTaskSAJF.cxx:275
 AliAnalysisTaskSAJF.cxx:276
 AliAnalysisTaskSAJF.cxx:277
 AliAnalysisTaskSAJF.cxx:278
 AliAnalysisTaskSAJF.cxx:279
 AliAnalysisTaskSAJF.cxx:280
 AliAnalysisTaskSAJF.cxx:281
 AliAnalysisTaskSAJF.cxx:282
 AliAnalysisTaskSAJF.cxx:283
 AliAnalysisTaskSAJF.cxx:284
 AliAnalysisTaskSAJF.cxx:285
 AliAnalysisTaskSAJF.cxx:286
 AliAnalysisTaskSAJF.cxx:287
 AliAnalysisTaskSAJF.cxx:288
 AliAnalysisTaskSAJF.cxx:289
 AliAnalysisTaskSAJF.cxx:290
 AliAnalysisTaskSAJF.cxx:291
 AliAnalysisTaskSAJF.cxx:292
 AliAnalysisTaskSAJF.cxx:293
 AliAnalysisTaskSAJF.cxx:294
 AliAnalysisTaskSAJF.cxx:295
 AliAnalysisTaskSAJF.cxx:296
 AliAnalysisTaskSAJF.cxx:297
 AliAnalysisTaskSAJF.cxx:298
 AliAnalysisTaskSAJF.cxx:299
 AliAnalysisTaskSAJF.cxx:300
 AliAnalysisTaskSAJF.cxx:301
 AliAnalysisTaskSAJF.cxx:302
 AliAnalysisTaskSAJF.cxx:303
 AliAnalysisTaskSAJF.cxx:304
 AliAnalysisTaskSAJF.cxx:305
 AliAnalysisTaskSAJF.cxx:306
 AliAnalysisTaskSAJF.cxx:307
 AliAnalysisTaskSAJF.cxx:308
 AliAnalysisTaskSAJF.cxx:309
 AliAnalysisTaskSAJF.cxx:310
 AliAnalysisTaskSAJF.cxx:311
 AliAnalysisTaskSAJF.cxx:312
 AliAnalysisTaskSAJF.cxx:313
 AliAnalysisTaskSAJF.cxx:314
 AliAnalysisTaskSAJF.cxx:315
 AliAnalysisTaskSAJF.cxx:316
 AliAnalysisTaskSAJF.cxx:317
 AliAnalysisTaskSAJF.cxx:318
 AliAnalysisTaskSAJF.cxx:319
 AliAnalysisTaskSAJF.cxx:320
 AliAnalysisTaskSAJF.cxx:321
 AliAnalysisTaskSAJF.cxx:322
 AliAnalysisTaskSAJF.cxx:323
 AliAnalysisTaskSAJF.cxx:324
 AliAnalysisTaskSAJF.cxx:325
 AliAnalysisTaskSAJF.cxx:326
 AliAnalysisTaskSAJF.cxx:327
 AliAnalysisTaskSAJF.cxx:328
 AliAnalysisTaskSAJF.cxx:329
 AliAnalysisTaskSAJF.cxx:330
 AliAnalysisTaskSAJF.cxx:331
 AliAnalysisTaskSAJF.cxx:332
 AliAnalysisTaskSAJF.cxx:333
 AliAnalysisTaskSAJF.cxx:334
 AliAnalysisTaskSAJF.cxx:335
 AliAnalysisTaskSAJF.cxx:336
 AliAnalysisTaskSAJF.cxx:337
 AliAnalysisTaskSAJF.cxx:338
 AliAnalysisTaskSAJF.cxx:339
 AliAnalysisTaskSAJF.cxx:340
 AliAnalysisTaskSAJF.cxx:341
 AliAnalysisTaskSAJF.cxx:342
 AliAnalysisTaskSAJF.cxx:343
 AliAnalysisTaskSAJF.cxx:344
 AliAnalysisTaskSAJF.cxx:345
 AliAnalysisTaskSAJF.cxx:346
 AliAnalysisTaskSAJF.cxx:347
 AliAnalysisTaskSAJF.cxx:348
 AliAnalysisTaskSAJF.cxx:349
 AliAnalysisTaskSAJF.cxx:350
 AliAnalysisTaskSAJF.cxx:351
 AliAnalysisTaskSAJF.cxx:352
 AliAnalysisTaskSAJF.cxx:353
 AliAnalysisTaskSAJF.cxx:354
 AliAnalysisTaskSAJF.cxx:355
 AliAnalysisTaskSAJF.cxx:356
 AliAnalysisTaskSAJF.cxx:357
 AliAnalysisTaskSAJF.cxx:358
 AliAnalysisTaskSAJF.cxx:359
 AliAnalysisTaskSAJF.cxx:360
 AliAnalysisTaskSAJF.cxx:361
 AliAnalysisTaskSAJF.cxx:362
 AliAnalysisTaskSAJF.cxx:363
 AliAnalysisTaskSAJF.cxx:364
 AliAnalysisTaskSAJF.cxx:365
 AliAnalysisTaskSAJF.cxx:366
 AliAnalysisTaskSAJF.cxx:367
 AliAnalysisTaskSAJF.cxx:368
 AliAnalysisTaskSAJF.cxx:369
 AliAnalysisTaskSAJF.cxx:370
 AliAnalysisTaskSAJF.cxx:371
 AliAnalysisTaskSAJF.cxx:372
 AliAnalysisTaskSAJF.cxx:373
 AliAnalysisTaskSAJF.cxx:374
 AliAnalysisTaskSAJF.cxx:375
 AliAnalysisTaskSAJF.cxx:376
 AliAnalysisTaskSAJF.cxx:377
 AliAnalysisTaskSAJF.cxx:378
 AliAnalysisTaskSAJF.cxx:379
 AliAnalysisTaskSAJF.cxx:380
 AliAnalysisTaskSAJF.cxx:381
 AliAnalysisTaskSAJF.cxx:382
 AliAnalysisTaskSAJF.cxx:383
 AliAnalysisTaskSAJF.cxx:384
 AliAnalysisTaskSAJF.cxx:385
 AliAnalysisTaskSAJF.cxx:386
 AliAnalysisTaskSAJF.cxx:387
 AliAnalysisTaskSAJF.cxx:388
 AliAnalysisTaskSAJF.cxx:389
 AliAnalysisTaskSAJF.cxx:390
 AliAnalysisTaskSAJF.cxx:391
 AliAnalysisTaskSAJF.cxx:392
 AliAnalysisTaskSAJF.cxx:393
 AliAnalysisTaskSAJF.cxx:394
 AliAnalysisTaskSAJF.cxx:395
 AliAnalysisTaskSAJF.cxx:396
 AliAnalysisTaskSAJF.cxx:397
 AliAnalysisTaskSAJF.cxx:398
 AliAnalysisTaskSAJF.cxx:399
 AliAnalysisTaskSAJF.cxx:400
 AliAnalysisTaskSAJF.cxx:401
 AliAnalysisTaskSAJF.cxx:402
 AliAnalysisTaskSAJF.cxx:403
 AliAnalysisTaskSAJF.cxx:404
 AliAnalysisTaskSAJF.cxx:405
 AliAnalysisTaskSAJF.cxx:406
 AliAnalysisTaskSAJF.cxx:407
 AliAnalysisTaskSAJF.cxx:408
 AliAnalysisTaskSAJF.cxx:409
 AliAnalysisTaskSAJF.cxx:410
 AliAnalysisTaskSAJF.cxx:411
 AliAnalysisTaskSAJF.cxx:412
 AliAnalysisTaskSAJF.cxx:413
 AliAnalysisTaskSAJF.cxx:414
 AliAnalysisTaskSAJF.cxx:415
 AliAnalysisTaskSAJF.cxx:416
 AliAnalysisTaskSAJF.cxx:417
 AliAnalysisTaskSAJF.cxx:418
 AliAnalysisTaskSAJF.cxx:419
 AliAnalysisTaskSAJF.cxx:420
 AliAnalysisTaskSAJF.cxx:421
 AliAnalysisTaskSAJF.cxx:422
 AliAnalysisTaskSAJF.cxx:423
 AliAnalysisTaskSAJF.cxx:424
 AliAnalysisTaskSAJF.cxx:425
 AliAnalysisTaskSAJF.cxx:426
 AliAnalysisTaskSAJF.cxx:427
 AliAnalysisTaskSAJF.cxx:428
 AliAnalysisTaskSAJF.cxx:429
 AliAnalysisTaskSAJF.cxx:430
 AliAnalysisTaskSAJF.cxx:431
 AliAnalysisTaskSAJF.cxx:432
 AliAnalysisTaskSAJF.cxx:433
 AliAnalysisTaskSAJF.cxx:434
 AliAnalysisTaskSAJF.cxx:435
 AliAnalysisTaskSAJF.cxx:436
 AliAnalysisTaskSAJF.cxx:437
 AliAnalysisTaskSAJF.cxx:438
 AliAnalysisTaskSAJF.cxx:439
 AliAnalysisTaskSAJF.cxx:440
 AliAnalysisTaskSAJF.cxx:441
 AliAnalysisTaskSAJF.cxx:442
 AliAnalysisTaskSAJF.cxx:443
 AliAnalysisTaskSAJF.cxx:444
 AliAnalysisTaskSAJF.cxx:445
 AliAnalysisTaskSAJF.cxx:446
 AliAnalysisTaskSAJF.cxx:447
 AliAnalysisTaskSAJF.cxx:448
 AliAnalysisTaskSAJF.cxx:449
 AliAnalysisTaskSAJF.cxx:450
 AliAnalysisTaskSAJF.cxx:451
 AliAnalysisTaskSAJF.cxx:452
 AliAnalysisTaskSAJF.cxx:453
 AliAnalysisTaskSAJF.cxx:454
 AliAnalysisTaskSAJF.cxx:455
 AliAnalysisTaskSAJF.cxx:456
 AliAnalysisTaskSAJF.cxx:457
 AliAnalysisTaskSAJF.cxx:458
 AliAnalysisTaskSAJF.cxx:459
 AliAnalysisTaskSAJF.cxx:460
 AliAnalysisTaskSAJF.cxx:461
 AliAnalysisTaskSAJF.cxx:462
 AliAnalysisTaskSAJF.cxx:463
 AliAnalysisTaskSAJF.cxx:464
 AliAnalysisTaskSAJF.cxx:465
 AliAnalysisTaskSAJF.cxx:466
 AliAnalysisTaskSAJF.cxx:467
 AliAnalysisTaskSAJF.cxx:468
 AliAnalysisTaskSAJF.cxx:469
 AliAnalysisTaskSAJF.cxx:470
 AliAnalysisTaskSAJF.cxx:471
 AliAnalysisTaskSAJF.cxx:472
 AliAnalysisTaskSAJF.cxx:473
 AliAnalysisTaskSAJF.cxx:474
 AliAnalysisTaskSAJF.cxx:475
 AliAnalysisTaskSAJF.cxx:476
 AliAnalysisTaskSAJF.cxx:477
 AliAnalysisTaskSAJF.cxx:478
 AliAnalysisTaskSAJF.cxx:479
 AliAnalysisTaskSAJF.cxx:480
 AliAnalysisTaskSAJF.cxx:481
 AliAnalysisTaskSAJF.cxx:482
 AliAnalysisTaskSAJF.cxx:483
 AliAnalysisTaskSAJF.cxx:484
 AliAnalysisTaskSAJF.cxx:485
 AliAnalysisTaskSAJF.cxx:486
 AliAnalysisTaskSAJF.cxx:487
 AliAnalysisTaskSAJF.cxx:488
 AliAnalysisTaskSAJF.cxx:489
 AliAnalysisTaskSAJF.cxx:490
 AliAnalysisTaskSAJF.cxx:491
 AliAnalysisTaskSAJF.cxx:492
 AliAnalysisTaskSAJF.cxx:493
 AliAnalysisTaskSAJF.cxx:494
 AliAnalysisTaskSAJF.cxx:495
 AliAnalysisTaskSAJF.cxx:496
 AliAnalysisTaskSAJF.cxx:497
 AliAnalysisTaskSAJF.cxx:498
 AliAnalysisTaskSAJF.cxx:499
 AliAnalysisTaskSAJF.cxx:500
 AliAnalysisTaskSAJF.cxx:501
 AliAnalysisTaskSAJF.cxx:502
 AliAnalysisTaskSAJF.cxx:503
 AliAnalysisTaskSAJF.cxx:504
 AliAnalysisTaskSAJF.cxx:505
 AliAnalysisTaskSAJF.cxx:506
 AliAnalysisTaskSAJF.cxx:507
 AliAnalysisTaskSAJF.cxx:508
 AliAnalysisTaskSAJF.cxx:509
 AliAnalysisTaskSAJF.cxx:510
 AliAnalysisTaskSAJF.cxx:511
 AliAnalysisTaskSAJF.cxx:512
 AliAnalysisTaskSAJF.cxx:513
 AliAnalysisTaskSAJF.cxx:514
 AliAnalysisTaskSAJF.cxx:515
 AliAnalysisTaskSAJF.cxx:516
 AliAnalysisTaskSAJF.cxx:517
 AliAnalysisTaskSAJF.cxx:518
 AliAnalysisTaskSAJF.cxx:519
 AliAnalysisTaskSAJF.cxx:520
 AliAnalysisTaskSAJF.cxx:521
 AliAnalysisTaskSAJF.cxx:522
 AliAnalysisTaskSAJF.cxx:523
 AliAnalysisTaskSAJF.cxx:524
 AliAnalysisTaskSAJF.cxx:525