[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / UserTasks / AliAnalysisTaskSAQA.cxx

// $Id$
//
// General QA task.
//
// Author: S.Aiola

#include <TChain.h>
#include <TClonesArray.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TH3F.h>
#include <TList.h>
#include <TLorentzVector.h>

#include "AliAnalysisManager.h"
#include "AliCentrality.h"
#include "AliVCluster.h"
#include "AliVParticle.h"
#include "AliVTrack.h"
#include "AliEmcalJet.h"
#include "AliVEventHandler.h"
#include "AliAODEvent.h"
#include "AliExternalTrackParam.h"
#include "AliTrackerBase.h"
#include "AliLog.h"

#include "AliAnalysisTaskSAQA.h"

ClassImp(AliAnalysisTaskSAQA)

//________________________________________________________________________
AliAnalysisTaskSAQA::AliAnalysisTaskSAQA() : 
  AliAnalysisTaskEmcalJet("AliAnalysisTaskSAQA", kTRUE),
  fCellEnergyCut(0.1),
  fNclusters(0),
  fNtracks(0),
  fNjets(0),
  fHistTracksCent(0),
  fHistClusCent(0),
  fHistJetsCent(0),
  fHistClusTracks(0),
  fHistJetsParts(0),
  fHistCellsCent(0),
  fHistCellsTracks(0),
  fHistTrEmcPhiEta(0),
  fHistTrPhiEtaNonProp(0),
  fHistDeltaEtaPt(0),
  fHistDeltaPhiPt(0),
  fHistNCellsEnergy(0),
  fHistClusTimeEnergy(0),
  fHistCellsAbsIdEnergy(0),
  fHistChVSneCells(0),
  fHistChVSneClus(0),
  fHistChVSneCorrCells(0)
{
  // Default constructor.

  for (Int_t i = 0; i < 5; i++) {
    fHistTrackPhi[i] = 0;
    fHistTrackEta[i] = 0;
  }

  for (Int_t i = 0; i < 4; i++) {
    fHistTrPhiEtaPt[i] = 0;
    fHistClusPhiEtaEnergy[i] = 0;
    fHistJetsPhiEtaPt[i] = 0;
    fHistJetsPtArea[i] = 0;
  }

  SetMakeGeneralHistograms(kTRUE);
}

//________________________________________________________________________
AliAnalysisTaskSAQA::AliAnalysisTaskSAQA(const char *name) : 
  AliAnalysisTaskEmcalJet(name, kTRUE),
  fCellEnergyCut(0.1),
  fNclusters(0),
  fNtracks(0),
  fNjets(0),
  fHistTracksCent(0),
  fHistClusCent(0),
  fHistJetsCent(0),
  fHistClusTracks(0),
  fHistJetsParts(0),
  fHistCellsCent(0),
  fHistCellsTracks(0),
  fHistTrEmcPhiEta(0),
  fHistTrPhiEtaNonProp(0),
  fHistDeltaEtaPt(0),
  fHistDeltaPhiPt(0),
  fHistNCellsEnergy(0),
  fHistClusTimeEnergy(0),
  fHistCellsAbsIdEnergy(0),
  fHistChVSneCells(0),
  fHistChVSneClus(0),
  fHistChVSneCorrCells(0)
{
  // Standard constructor.

  for (Int_t i = 0; i < 5; i++) {
    fHistTrackPhi[i] = 0;
    fHistTrackEta[i] = 0;
  }

  for (Int_t i = 0; i < 4; i++) {
    fHistTrPhiEtaPt[i] = 0;
    fHistClusPhiEtaEnergy[i] = 0;
    fHistJetsPhiEtaPt[i] = 0;
    fHistJetsPtArea[i] = 0;
  }

  SetMakeGeneralHistograms(kTRUE);
}

//________________________________________________________________________
AliAnalysisTaskSAQA::~AliAnalysisTaskSAQA()
{
  // Destructor
}

//________________________________________________________________________
void AliAnalysisTaskSAQA::UserCreateOutputObjects()
{
  // Create histograms

  AliAnalysisTaskEmcalJet::UserCreateOutputObjects();

  if (!fTracksName.IsNull()) {
    fHistTracksCent = new TH2F("fHistTracksCent","Tracks vs. centrality", 100, 0, 100, fNbins, 0, 4000);
    fHistTracksCent->GetXaxis()->SetTitle("Centrality (%)");
    fHistTracksCent->GetYaxis()->SetTitle("No. of tracks");
    fOutput->Add(fHistTracksCent);

    TString histname;

    for (Int_t i = 0; i < 4; i++) {
      histname = "fHistTrPhiEtaPt_";
      histname += i;
      fHistTrPhiEtaPt[i] = new TH3F(histname,histname, 100, -1, 1, 201, 0, TMath::Pi() * 2.01, fNbins, fMinBinPt, fMaxBinPt);
      fHistTrPhiEtaPt[i] ->GetXaxis()->SetTitle("#eta");
      fHistTrPhiEtaPt[i] ->GetYaxis()->SetTitle("#phi");
      fHistTrPhiEtaPt[i] ->GetZaxis()->SetTitle("p_{T} (GeV/c)");
      fOutput->Add(fHistTrPhiEtaPt[i]);
    }

    fHistTrEmcPhiEta = new TH2F("fHistTrEmcPhiEta","Phi-Eta emcal distribution of tracks", 100, -1, 1, 201, 0, TMath::Pi() * 2.01);
    fHistTrEmcPhiEta->GetXaxis()->SetTitle("#eta");
    fHistTrEmcPhiEta->GetYaxis()->SetTitle("#phi");
    fOutput->Add(fHistTrEmcPhiEta);

    fHistTrPhiEtaNonProp = new TH2F("fHistTrPhiEtaNonProp","fHistTrPhiEtaNonProp", 100, -1, 1, 201, 0, TMath::Pi() * 2.01);
    fHistTrPhiEtaNonProp->GetXaxis()->SetTitle("#eta");
    fHistTrPhiEtaNonProp->GetYaxis()->SetTitle("#phi");
    fOutput->Add(fHistTrPhiEtaNonProp);
    
    fHistDeltaEtaPt = new TH2F("fHistDeltaEtaPt","fHistDeltaEtaPt", fNbins, fMinBinPt, fMaxBinPt, 80, -0.5, 0.5);
    fHistDeltaEtaPt->GetXaxis()->SetTitle("p_{T} [GeV/c]");
    fHistDeltaEtaPt->GetYaxis()->SetTitle("#delta#eta");
    fOutput->Add(fHistDeltaEtaPt);
    
    fHistDeltaPhiPt = new TH2F("fHistDeltaPhiPt","fHistDeltaPhiPt", fNbins, fMinBinPt, fMaxBinPt, 256, -1.6, 4.8);
    fHistDeltaPhiPt->GetXaxis()->SetTitle("p_{T} [GeV/c]");
    fHistDeltaPhiPt->GetYaxis()->SetTitle("#delta#phi");
    fOutput->Add(fHistDeltaPhiPt);

    for (Int_t i = 0; i < 5; i++) {
      TString histnamephi("fHistTrackPhi_");
      histnamephi += i;
      fHistTrackPhi[i] = new TH1F(histnamephi.Data(),histnamephi.Data(), 201, 0, TMath::Pi() * 2.01);
      fHistTrackPhi[i]->GetXaxis()->SetTitle("Phi");
      fOutput->Add(fHistTrackPhi[i]);
      
      TString histnameeta("fHistTrackEta_");
      histnameeta += i;
      fHistTrackEta[i] = new TH1F(histnameeta.Data(),histnameeta.Data(), 100, -1, 1);
      fHistTrackEta[i]->GetXaxis()->SetTitle("Eta");
      fOutput->Add(fHistTrackEta[i]);
    }

    fHistTrackPhi[0]->SetLineColor(kRed);
    fHistTrackEta[0]->SetLineColor(kRed);
    fHistTrackPhi[1]->SetLineColor(kBlue);
    fHistTrackEta[1]->SetLineColor(kBlue);
    fHistTrackPhi[2]->SetLineColor(kGreen);
    fHistTrackEta[2]->SetLineColor(kGreen);
    fHistTrackPhi[3]->SetLineColor(kOrange);
    fHistTrackEta[3]->SetLineColor(kOrange);
    fHistTrackPhi[4]->SetLineColor(kBlack);
    fHistTrackEta[4]->SetLineColor(kBlack);
  }

  if (!fCaloName.IsNull()) {
    fHistClusCent = new TH2F("fHistClusCent","Clusters vs. centrality", 100, 0, 100, fNbins, 0, 2000);
    fHistClusCent->GetXaxis()->SetTitle("Centrality (%)");
    fHistClusCent->GetYaxis()->SetTitle("No. of clusters");
    fOutput->Add(fHistClusCent);

    fHistClusTracks = new TH2F("fHistClusTracks","Clusters vs. tracks", fNbins, 0, 4000, fNbins, 0, 2000);
    fHistClusTracks->GetXaxis()->SetTitle("No. of tracks");
    fHistClusTracks->GetYaxis()->SetTitle("No. of clusters");
    fOutput->Add(fHistClusTracks);

    fHistCellsCent = new TH2F("fHistCellsCent","Cells vs. centrality", 100, 0, 100, fNbins, 0, 6000);
    fHistCellsCent->GetXaxis()->SetTitle("Centrality (%)");
    fHistCellsCent->GetYaxis()->SetTitle("No. of EMCal cells");
    fOutput->Add(fHistCellsCent);

    fHistCellsTracks = new TH2F("fHistCellsTracks","Cells vs. tracks", fNbins, 0, 4000, fNbins, 0, 6000);
    fHistCellsTracks->GetXaxis()->SetTitle("No. of tracks");
    fHistCellsTracks->GetYaxis()->SetTitle("No. of EMCal cells");
    fOutput->Add(fHistCellsTracks);

    TString histname;

    for (Int_t i = 0; i < 4; i++) {
      histname = "fHistClusPhiEtaEnergy_";
      histname += i;
      fHistClusPhiEtaEnergy[i] = new TH3F(histname, histname, 100, -1.2, 1.2, 201, 0, TMath::Pi() * 2.01, fNbins, fMinBinPt, fMaxBinPt);
      fHistClusPhiEtaEnergy[i]->GetXaxis()->SetTitle("#eta");
      fHistClusPhiEtaEnergy[i]->GetYaxis()->SetTitle("#phi");
      fHistClusPhiEtaEnergy[i]->GetZaxis()->SetTitle("Energy (GeV)");
      fOutput->Add(fHistClusPhiEtaEnergy[i]);
    }

    fHistClusTimeEnergy = new TH2F("fHistClusTimeEnergy","Time vs. energy of clusters", fNbins, fMinBinPt, fMaxBinPt, fNbins,  -1e-6, 1e-6);
    fHistClusTimeEnergy->GetXaxis()->SetTitle("Energy (GeV)");
    fHistClusTimeEnergy->GetYaxis()->SetTitle("Time");
    fOutput->Add(fHistClusTimeEnergy);

    fHistNCellsEnergy = new TH2F("fHistNCellsEnergy","Number of cells vs. energy of clusters", fNbins, fMinBinPt, fMaxBinPt, 30, 0, 30);
    fHistNCellsEnergy->GetXaxis()->SetTitle("Energy (GeV)");
    fHistNCellsEnergy->GetYaxis()->SetTitle("N_{cells}");
    fOutput->Add(fHistNCellsEnergy); 
     
    fHistCellsAbsIdEnergy = new TH2F("fHistCellsAbsIdEnergy","fHistCellsAbsIdEnergy", 11600,0,11599,(Int_t)(fNbins / 2), fMinBinPt, fMaxBinPt / 2);
    fHistCellsAbsIdEnergy->GetXaxis()->SetTitle("cell abs. Id");
    fHistCellsAbsIdEnergy->GetYaxis()->SetTitle("Energy (GeV)");
    fHistCellsAbsIdEnergy->GetZaxis()->SetTitle("counts");    
    fOutput->Add(fHistCellsAbsIdEnergy);
    
    fHistChVSneCells = new TH2F("fHistChVSneCells","Charged energy vs. neutral (cells) energy", 
				(Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5, (Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5);
    fHistChVSneCells->GetXaxis()->SetTitle("Energy (GeV)");
    fHistChVSneCells->GetYaxis()->SetTitle("Momentum (GeV/c)");
    fOutput->Add(fHistChVSneCells);
    
    fHistChVSneClus = new TH2F("fHistChVSneClus","Charged energy vs. neutral (clusters) energy", 
			       (Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5, (Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5);
    fHistChVSneClus->GetXaxis()->SetTitle("Energy (GeV)");
    fHistChVSneClus->GetYaxis()->SetTitle("Momentum (GeV/c)");
    fOutput->Add(fHistChVSneClus);
    
    fHistChVSneCorrCells = new TH2F("fHistChVSneCorrCells","Charged energy vs. neutral (corrected cells) energy", 
				    (Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5, (Int_t)(fNbins * 2.5), fMinBinPt , fMaxBinPt * 2.5);
    fHistChVSneCorrCells->GetXaxis()->SetTitle("Energy (GeV)");
    fHistChVSneCorrCells->GetYaxis()->SetTitle("Momentum (GeV/c)");
    fOutput->Add(fHistChVSneCorrCells);
  }
       
  if (!fJetsName.IsNull()) {
    fHistJetsCent = new TH2F("fHistJetsCent","Jets vs. centrality", 100, 0, 100, 150, -0.5, 149.5);
    fHistJetsCent->GetXaxis()->SetTitle("Centrality (%)");
    fHistJetsCent->GetYaxis()->SetTitle("No. of jets");
    fOutput->Add(fHistJetsCent);
    
    fHistJetsParts = new TH2F("fHistJetsParts","Jets vs. centrality", fNbins, 0, 6000, 150, -0.5, 149.5);
    fHistJetsParts->GetXaxis()->SetTitle("No. of particles");
    fHistJetsParts->GetYaxis()->SetTitle("No. of jets");
    fOutput->Add(fHistJetsParts);

    TString histname;

    for (Int_t i = 0; i < 4; i++) {
      histname = "fHistJetsPhiEtaPt_";
      histname += i;
      fHistJetsPhiEtaPt[i] = new TH3F(histname.Data(), histname.Data(), 100, -1.2, 1.2, 201, 0, TMath::Pi() * 2.01,  (Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5);
      fHistJetsPhiEtaPt[i]->GetXaxis()->SetTitle("#eta");
      fHistJetsPhiEtaPt[i]->GetYaxis()->SetTitle("#phi");
      fHistJetsPhiEtaPt[i]->GetZaxis()->SetTitle("p_{T} (GeV/c)");
      fOutput->Add(fHistJetsPhiEtaPt[i]);

      histname = "fHistJetsPtArea_";
      histname += i;
      fHistJetsPtArea[i] = new TH2F(histname.Data(), histname.Data(), (Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5, 30, 0, fJetRadius * fJetRadius * TMath::Pi() * 3);
      fHistJetsPtArea[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
      fHistJetsPtArea[i]->GetYaxis()->SetTitle("area");
      fOutput->Add(fHistJetsPtArea[i]);
    }
  }

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

//________________________________________________________________________
Bool_t AliAnalysisTaskSAQA::RetrieveEventObjects()
{
  // Retrieve event objects.

  if (!AliAnalysisTaskEmcalJet::RetrieveEventObjects())
    return kFALSE;

  fNclusters = 0;
  fNtracks = 0;
  fNjets = 0;

  return kTRUE;
}


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

  Float_t clusSum = 0;
  Float_t trackSum = 0;

  if (fTracks) {
    trackSum = DoTrackLoop();

    fHistTracksCent->Fill(fCent, fNtracks);
  } 

  if (fCaloClusters) {
    clusSum = DoClusterLoop();

    fHistClusCent->Fill(fCent, fNclusters);
    fHistClusTracks->Fill(fNtracks, fNclusters);

    Float_t cellSum = 0, cellCutSum = 0;
    
    Int_t ncells = DoCellLoop(cellSum, cellCutSum);

    if (fTracks)
      fHistCellsTracks->Fill(fNtracks, ncells);

    fHistCellsCent->Fill(fCent, ncells);
    
    fHistChVSneCells->Fill(cellSum, trackSum);
    fHistChVSneClus->Fill(clusSum, trackSum);
    fHistChVSneCorrCells->Fill(cellCutSum, trackSum);
  }

  if (fJets) {
    DoJetLoop();
    
    fHistJetsCent->Fill(fCent, fNjets);
    fHistJetsParts->Fill(fNtracks + fNclusters, fNjets);
  }

  return kTRUE;
}

//________________________________________________________________________
Int_t AliAnalysisTaskSAQA::DoCellLoop(Float_t &sum, Float_t &sum_cut)
{
  // Do cell loop.

  AliVCaloCells *cells = InputEvent()->GetEMCALCells();

  if (!cells)
    return 0;

  const Int_t ncells = cells->GetNumberOfCells();

  for (Int_t pos = 0; pos < ncells; pos++) {
    Float_t amp   = cells->GetAmplitude(pos);
    Int_t   absId = cells->GetCellNumber(pos);
    fHistCellsAbsIdEnergy->Fill(absId,amp);
    sum += amp;
    if (amp < fCellEnergyCut)
      continue;
    sum_cut += amp;
  } 

  return ncells;
}

//________________________________________________________________________
Float_t AliAnalysisTaskSAQA::DoClusterLoop()
{
  // Do cluster loop.

  if (!fCaloClusters)
    return 0;

  Float_t sum = 0;

  // Cluster loop
  Int_t nclusters =  fCaloClusters->GetEntriesFast();

  for (Int_t iClusters = 0; iClusters < nclusters; iClusters++) {
    AliVCluster* cluster = static_cast<AliVCluster*>(fCaloClusters->At(iClusters));
    if (!cluster) {
      AliError(Form("Could not receive cluster %d", iClusters));
      continue;
    }  

    if (!AcceptCluster(cluster, kTRUE))
      continue;

    sum += cluster->E();

    TLorentzVector nPart;
    cluster->GetMomentum(nPart, fVertex);

    fHistClusPhiEtaEnergy[fCentBin]->Fill(nPart.Eta(), nPart.Phi(), cluster->E());
    fHistNCellsEnergy->Fill(cluster->E(), cluster->GetNCells());

    fHistClusTimeEnergy->Fill(cluster->E(), cluster->GetTOF());

    fNclusters++;
  }

  return sum;
}

//________________________________________________________________________
Float_t AliAnalysisTaskSAQA::DoTrackLoop()
{
  // Do track loop.

  if (!fTracks)
    return 0;

  Float_t sum = 0;

  Int_t ntracks = fTracks->GetEntriesFast();
  Int_t nclusters = 0;
  if (fCaloClusters)
    nclusters = fCaloClusters->GetEntriesFast();

  for (Int_t i = 0; i < ntracks; i++) {

    AliVParticle* track = static_cast<AliVParticle*>(fTracks->At(i)); // pointer to reconstructed to track  

    if (!track) {
      AliError(Form("Could not retrieve track %d",i)); 
      continue; 
    }

    AliVTrack* vtrack = dynamic_cast<AliVTrack*>(track); 
    
    if (vtrack && !AcceptTrack(vtrack, kTRUE)) 
      continue;

    fNtracks++;

    sum += track->P();
      
    fHistTrPhiEtaPt[fCentBin]->Fill(track->Eta(), track->Phi(), track->Pt());

    Int_t label = track->GetLabel();

    if (label >= 0 && label < 4) {
      fHistTrackEta[label]->Fill(track->Eta());
      fHistTrackPhi[label]->Fill(track->Phi());
    }

    fHistTrackEta[4]->Fill(track->Eta());
    fHistTrackPhi[4]->Fill(track->Phi());

    if (!vtrack)
      continue;

    if (vtrack->GetTrackEtaOnEMCal() == -999 || vtrack->GetTrackPhiOnEMCal() == -999)
      fHistTrPhiEtaNonProp->Fill(vtrack->Eta(), vtrack->Phi());

    fHistTrEmcPhiEta->Fill(vtrack->GetTrackEtaOnEMCal(), vtrack->GetTrackPhiOnEMCal());
    fHistDeltaEtaPt->Fill(vtrack->Pt(), vtrack->Eta() - vtrack->GetTrackEtaOnEMCal());
    fHistDeltaPhiPt->Fill(vtrack->Pt(), vtrack->Phi() - vtrack->GetTrackPhiOnEMCal());
  }
  
  return sum;
}

//________________________________________________________________________
void AliAnalysisTaskSAQA::DoJetLoop()
{
  // Do jet loop.

  if (!fJets)
    return;

  Int_t njets = fJets->GetEntriesFast();

  for (Int_t ij = 0; ij < njets; ij++) {

    AliEmcalJet* jet = static_cast<AliEmcalJet*>(fJets->At(ij));

    if (!jet) {
      AliError(Form("Could not receive jet %d", ij));
      continue;
    }  

    if (!AcceptJet(jet))
      continue;

    fNjets++;

    fHistJetsPhiEtaPt[fCentBin]->Fill(jet->Eta(), jet->Phi(), jet->Pt());
    fHistJetsPtArea[fCentBin]->Fill(jet->Pt(), jet->Area());
  }
}


//________________________________________________________________________
void AliAnalysisTaskSAQA::Terminate(Option_t *) 
{
  // Called once at the end of the analysis.
}
Commit	Line	Data
6e8d91c9	1	// $Id$
c3ba2d3d	2	//
297edd60	3	// General QA task.
c3ba2d3d	4	//
297edd60	5	// Author: S.Aiola
c3ba2d3d	6
	7	#include <TChain.h>
	8	#include <TClonesArray.h>
	9	#include <TH1F.h>
	10	#include <TH2F.h>
3000c095	11	#include <TH3F.h>
c3ba2d3d	12	#include <TList.h>
	13	#include <TLorentzVector.h>
	14
	15	#include "AliAnalysisManager.h"
	16	#include "AliCentrality.h"
	17	#include "AliVCluster.h"
	18	#include "AliVParticle.h"
	19	#include "AliVTrack.h"
	20	#include "AliEmcalJet.h"
	21	#include "AliVEventHandler.h"
079b4732	22	#include "AliAODEvent.h"
	23	#include "AliExternalTrackParam.h"
	24	#include "AliTrackerBase.h"
c3ba2d3d	25	#include "AliLog.h"
	26
	27	#include "AliAnalysisTaskSAQA.h"
	28
	29	ClassImp(AliAnalysisTaskSAQA)
	30
	31	//________________________________________________________________________
	32	AliAnalysisTaskSAQA::AliAnalysisTaskSAQA() :
99cfd012	33	AliAnalysisTaskEmcalJet("AliAnalysisTaskSAQA", kTRUE),
c3ba2d3d	34	fCellEnergyCut(0.1),
7c1d624c	35	fNclusters(0),
	36	fNtracks(0),
	37	fNjets(0),
c3ba2d3d	38	fHistTracksCent(0),
c3ba2d3d	39	fHistClusCent(0),
7c1d624c	40	fHistJetsCent(0),
090a0c3e	41	fHistClusTracks(0),
7c1d624c	42	fHistJetsParts(0),
be7b6e63	43	fHistCellsCent(0),
be7b6e63	44	fHistCellsTracks(0),
6fd5039f	45	fHistTrEmcPhiEta(0),
079b4732	46	fHistTrPhiEtaNonProp(0),
3e951e32	47	fHistDeltaEtaPt(0),
3e951e32	48	fHistDeltaPhiPt(0),
7cf4626b	49	fHistNCellsEnergy(0),
090a0c3e	50	fHistClusTimeEnergy(0),
f483218e	51	fHistCellsAbsIdEnergy(0),
c3ba2d3d	52	fHistChVSneCells(0),
	53	fHistChVSneClus(0),
	54	fHistChVSneCorrCells(0)
	55	{
	56	// Default constructor.
	57
	58	for (Int_t i = 0; i < 5; i++) {
	59	fHistTrackPhi[i] = 0;
	60	fHistTrackEta[i] = 0;
	61	}
6fd5039f	62
6fd5039f	63	for (Int_t i = 0; i < 4; i++) {
a487deae	64	fHistTrPhiEtaPt[i] = 0;
	65	fHistClusPhiEtaEnergy[i] = 0;
	66	fHistJetsPhiEtaPt[i] = 0;
98750b70	67	fHistJetsPtArea[i] = 0;
6fd5039f	68	}
a487deae	69
a487deae	70	SetMakeGeneralHistograms(kTRUE);
c3ba2d3d	71	}
	72
	73	//________________________________________________________________________
	74	AliAnalysisTaskSAQA::AliAnalysisTaskSAQA(const char *name) :
99cfd012	75	AliAnalysisTaskEmcalJet(name, kTRUE),
c3ba2d3d	76	fCellEnergyCut(0.1),
7c1d624c	77	fNclusters(0),
	78	fNtracks(0),
	79	fNjets(0),
c3ba2d3d	80	fHistTracksCent(0),
c3ba2d3d	81	fHistClusCent(0),
7c1d624c	82	fHistJetsCent(0),
090a0c3e	83	fHistClusTracks(0),
7c1d624c	84	fHistJetsParts(0),
be7b6e63	85	fHistCellsCent(0),
be7b6e63	86	fHistCellsTracks(0),
6fd5039f	87	fHistTrEmcPhiEta(0),
079b4732	88	fHistTrPhiEtaNonProp(0),
3e951e32	89	fHistDeltaEtaPt(0),
3e951e32	90	fHistDeltaPhiPt(0),
7cf4626b	91	fHistNCellsEnergy(0),
090a0c3e	92	fHistClusTimeEnergy(0),
f483218e	93	fHistCellsAbsIdEnergy(0),
c3ba2d3d	94	fHistChVSneCells(0),
	95	fHistChVSneClus(0),
	96	fHistChVSneCorrCells(0)
	97	{
	98	// Standard constructor.
	99
	100	for (Int_t i = 0; i < 5; i++) {
	101	fHistTrackPhi[i] = 0;
	102	fHistTrackEta[i] = 0;
	103	}
6fd5039f	104
6fd5039f	105	for (Int_t i = 0; i < 4; i++) {
a487deae	106	fHistTrPhiEtaPt[i] = 0;
	107	fHistClusPhiEtaEnergy[i] = 0;
	108	fHistJetsPhiEtaPt[i] = 0;
98750b70	109	fHistJetsPtArea[i] = 0;
6fd5039f	110	}
a487deae	111
a487deae	112	SetMakeGeneralHistograms(kTRUE);
c3ba2d3d	113	}
	114
	115	//________________________________________________________________________
	116	AliAnalysisTaskSAQA::~AliAnalysisTaskSAQA()
	117	{
	118	// Destructor
	119	}
	120
	121	//________________________________________________________________________
	122	void AliAnalysisTaskSAQA::UserCreateOutputObjects()
	123	{
	124	// Create histograms
	125
a487deae	126	AliAnalysisTaskEmcalJet::UserCreateOutputObjects();
c3ba2d3d	127
a487deae	128	if (!fTracksName.IsNull()) {
	129	fHistTracksCent = new TH2F("fHistTracksCent","Tracks vs. centrality", 100, 0, 100, fNbins, 0, 4000);
	130	fHistTracksCent->GetXaxis()->SetTitle("Centrality (%)");
	131	fHistTracksCent->GetYaxis()->SetTitle("No. of tracks");
	132	fOutput->Add(fHistTracksCent);
c3ba2d3d	133
a487deae	134	TString histname;
be7b6e63	135
a487deae	136	for (Int_t i = 0; i < 4; i++) {
	137	histname = "fHistTrPhiEtaPt_";
	138	histname += i;
	139	fHistTrPhiEtaPt[i] = new TH3F(histname,histname, 100, -1, 1, 201, 0, TMath::Pi() * 2.01, fNbins, fMinBinPt, fMaxBinPt);
	140	fHistTrPhiEtaPt[i] ->GetXaxis()->SetTitle("#eta");
	141	fHistTrPhiEtaPt[i] ->GetYaxis()->SetTitle("#phi");
	142	fHistTrPhiEtaPt[i] ->GetZaxis()->SetTitle("p_{T} (GeV/c)");
	143	fOutput->Add(fHistTrPhiEtaPt[i]);
	144	}
c3ba2d3d	145
a487deae	146	fHistTrEmcPhiEta = new TH2F("fHistTrEmcPhiEta","Phi-Eta emcal distribution of tracks", 100, -1, 1, 201, 0, TMath::Pi() * 2.01);
	147	fHistTrEmcPhiEta->GetXaxis()->SetTitle("#eta");
	148	fHistTrEmcPhiEta->GetYaxis()->SetTitle("#phi");
	149	fOutput->Add(fHistTrEmcPhiEta);
	150
	151	fHistTrPhiEtaNonProp = new TH2F("fHistTrPhiEtaNonProp","fHistTrPhiEtaNonProp", 100, -1, 1, 201, 0, TMath::Pi() * 2.01);
	152	fHistTrPhiEtaNonProp->GetXaxis()->SetTitle("#eta");
	153	fHistTrPhiEtaNonProp->GetYaxis()->SetTitle("#phi");
	154	fOutput->Add(fHistTrPhiEtaNonProp);
	155
	156	fHistDeltaEtaPt = new TH2F("fHistDeltaEtaPt","fHistDeltaEtaPt", fNbins, fMinBinPt, fMaxBinPt, 80, -0.5, 0.5);
	157	fHistDeltaEtaPt->GetXaxis()->SetTitle("p_{T} [GeV/c]");
	158	fHistDeltaEtaPt->GetYaxis()->SetTitle("#delta#eta");
	159	fOutput->Add(fHistDeltaEtaPt);
	160
	161	fHistDeltaPhiPt = new TH2F("fHistDeltaPhiPt","fHistDeltaPhiPt", fNbins, fMinBinPt, fMaxBinPt, 256, -1.6, 4.8);
	162	fHistDeltaPhiPt->GetXaxis()->SetTitle("p_{T} [GeV/c]");
	163	fHistDeltaPhiPt->GetYaxis()->SetTitle("#delta#phi");
	164	fOutput->Add(fHistDeltaPhiPt);
	165
	166	for (Int_t i = 0; i < 5; i++) {
	167	TString histnamephi("fHistTrackPhi_");
	168	histnamephi += i;
	169	fHistTrackPhi[i] = new TH1F(histnamephi.Data(),histnamephi.Data(), 201, 0, TMath::Pi() * 2.01);
	170	fHistTrackPhi[i]->GetXaxis()->SetTitle("Phi");
	171	fOutput->Add(fHistTrackPhi[i]);
	172
	173	TString histnameeta("fHistTrackEta_");
	174	histnameeta += i;
	175	fHistTrackEta[i] = new TH1F(histnameeta.Data(),histnameeta.Data(), 100, -1, 1);
	176	fHistTrackEta[i]->GetXaxis()->SetTitle("Eta");
	177	fOutput->Add(fHistTrackEta[i]);
	178	}
7c1d624c	179
a487deae	180	fHistTrackPhi[0]->SetLineColor(kRed);
	181	fHistTrackEta[0]->SetLineColor(kRed);
	182	fHistTrackPhi[1]->SetLineColor(kBlue);
	183	fHistTrackEta[1]->SetLineColor(kBlue);
	184	fHistTrackPhi[2]->SetLineColor(kGreen);
	185	fHistTrackEta[2]->SetLineColor(kGreen);
	186	fHistTrackPhi[3]->SetLineColor(kOrange);
	187	fHistTrackEta[3]->SetLineColor(kOrange);
	188	fHistTrackPhi[4]->SetLineColor(kBlack);
	189	fHistTrackEta[4]->SetLineColor(kBlack);
	190	}
7c1d624c	191
a487deae	192	if (!fCaloName.IsNull()) {
a5621834	193	fHistClusCent = new TH2F("fHistClusCent","Clusters vs. centrality", 100, 0, 100, fNbins, 0, 2000);
6fd5039f	194	fHistClusCent->GetXaxis()->SetTitle("Centrality (%)");
	195	fHistClusCent->GetYaxis()->SetTitle("No. of clusters");
	196	fOutput->Add(fHistClusCent);
090a0c3e	197
	198	fHistClusTracks = new TH2F("fHistClusTracks","Clusters vs. tracks", fNbins, 0, 4000, fNbins, 0, 2000);
	199	fHistClusTracks->GetXaxis()->SetTitle("No. of tracks");
	200	fHistClusTracks->GetYaxis()->SetTitle("No. of clusters");
	201	fOutput->Add(fHistClusTracks);
	202
be7b6e63	203	fHistCellsCent = new TH2F("fHistCellsCent","Cells vs. centrality", 100, 0, 100, fNbins, 0, 6000);
	204	fHistCellsCent->GetXaxis()->SetTitle("Centrality (%)");
	205	fHistCellsCent->GetYaxis()->SetTitle("No. of EMCal cells");
	206	fOutput->Add(fHistCellsCent);
	207
	208	fHistCellsTracks = new TH2F("fHistCellsTracks","Cells vs. tracks", fNbins, 0, 4000, fNbins, 0, 6000);
	209	fHistCellsTracks->GetXaxis()->SetTitle("No. of tracks");
	210	fHistCellsTracks->GetYaxis()->SetTitle("No. of EMCal cells");
	211	fOutput->Add(fHistCellsTracks);
	212
a487deae	213	TString histname;
	214
	215	for (Int_t i = 0; i < 4; i++) {
	216	histname = "fHistClusPhiEtaEnergy_";
	217	histname += i;
	218	fHistClusPhiEtaEnergy[i] = new TH3F(histname, histname, 100, -1.2, 1.2, 201, 0, TMath::Pi() * 2.01, fNbins, fMinBinPt, fMaxBinPt);
	219	fHistClusPhiEtaEnergy[i]->GetXaxis()->SetTitle("#eta");
	220	fHistClusPhiEtaEnergy[i]->GetYaxis()->SetTitle("#phi");
	221	fHistClusPhiEtaEnergy[i]->GetZaxis()->SetTitle("Energy (GeV)");
	222	fOutput->Add(fHistClusPhiEtaEnergy[i]);
	223	}
	224
781da0a3	225	fHistClusTimeEnergy = new TH2F("fHistClusTimeEnergy","Time vs. energy of clusters", fNbins, fMinBinPt, fMaxBinPt, fNbins, -1e-6, 1e-6);
090a0c3e	226	fHistClusTimeEnergy->GetXaxis()->SetTitle("Energy (GeV)");
	227	fHistClusTimeEnergy->GetYaxis()->SetTitle("Time");
	228	fOutput->Add(fHistClusTimeEnergy);
	229
7cf4626b	230	fHistNCellsEnergy = new TH2F("fHistNCellsEnergy","Number of cells vs. energy of clusters", fNbins, fMinBinPt, fMaxBinPt, 30, 0, 30);
a487deae	231	fHistNCellsEnergy->GetXaxis()->SetTitle("Energy (GeV)");
7cf4626b	232	fHistNCellsEnergy->GetYaxis()->SetTitle("N_{cells}");
a487deae	233	fOutput->Add(fHistNCellsEnergy);
a487deae	234
f483218e	235	fHistCellsAbsIdEnergy = new TH2F("fHistCellsAbsIdEnergy","fHistCellsAbsIdEnergy", 11600,0,11599,(Int_t)(fNbins / 2), fMinBinPt, fMaxBinPt / 2);
	236	fHistCellsAbsIdEnergy->GetXaxis()->SetTitle("cell abs. Id");
	237	fHistCellsAbsIdEnergy->GetYaxis()->SetTitle("Energy (GeV)");
	238	fHistCellsAbsIdEnergy->GetZaxis()->SetTitle("counts");
	239	fOutput->Add(fHistCellsAbsIdEnergy);
6fd5039f	240
98750b70	241	fHistChVSneCells = new TH2F("fHistChVSneCells","Charged energy vs. neutral (cells) energy",
be7b6e63	242	(Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5, (Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5);
a487deae	243	fHistChVSneCells->GetXaxis()->SetTitle("Energy (GeV)");
a487deae	244	fHistChVSneCells->GetYaxis()->SetTitle("Momentum (GeV/c)");
6fd5039f	245	fOutput->Add(fHistChVSneCells);
6fd5039f	246
98750b70	247	fHistChVSneClus = new TH2F("fHistChVSneClus","Charged energy vs. neutral (clusters) energy",
be7b6e63	248	(Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5, (Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5);
a487deae	249	fHistChVSneClus->GetXaxis()->SetTitle("Energy (GeV)");
a487deae	250	fHistChVSneClus->GetYaxis()->SetTitle("Momentum (GeV/c)");
6fd5039f	251	fOutput->Add(fHistChVSneClus);
6fd5039f	252
98750b70	253	fHistChVSneCorrCells = new TH2F("fHistChVSneCorrCells","Charged energy vs. neutral (corrected cells) energy",
be7b6e63	254	(Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5, (Int_t)(fNbins * 2.5), fMinBinPt , fMaxBinPt * 2.5);
a487deae	255	fHistChVSneCorrCells->GetXaxis()->SetTitle("Energy (GeV)");
a487deae	256	fHistChVSneCorrCells->GetYaxis()->SetTitle("Momentum (GeV/c)");
6fd5039f	257	fOutput->Add(fHistChVSneCorrCells);
c3ba2d3d	258	}
a487deae	259
e44e8726	260	if (!fJetsName.IsNull()) {
a487deae	261	fHistJetsCent = new TH2F("fHistJetsCent","Jets vs. centrality", 100, 0, 100, 150, -0.5, 149.5);
	262	fHistJetsCent->GetXaxis()->SetTitle("Centrality (%)");
	263	fHistJetsCent->GetYaxis()->SetTitle("No. of jets");
	264	fOutput->Add(fHistJetsCent);
	265
	266	fHistJetsParts = new TH2F("fHistJetsParts","Jets vs. centrality", fNbins, 0, 6000, 150, -0.5, 149.5);
	267	fHistJetsParts->GetXaxis()->SetTitle("No. of particles");
	268	fHistJetsParts->GetYaxis()->SetTitle("No. of jets");
	269	fOutput->Add(fHistJetsParts);
6fd5039f	270
e44e8726	271	TString histname;
	272
	273	for (Int_t i = 0; i < 4; i++) {
a487deae	274	histname = "fHistJetsPhiEtaPt_";
e44e8726	275	histname += i;
a487deae	276	fHistJetsPhiEtaPt[i] = new TH3F(histname.Data(), histname.Data(), 100, -1.2, 1.2, 201, 0, TMath::Pi() * 2.01, (Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5);
	277	fHistJetsPhiEtaPt[i]->GetXaxis()->SetTitle("#eta");
	278	fHistJetsPhiEtaPt[i]->GetYaxis()->SetTitle("#phi");
	279	fHistJetsPhiEtaPt[i]->GetZaxis()->SetTitle("p_{T} (GeV/c)");
	280	fOutput->Add(fHistJetsPhiEtaPt[i]);
e44e8726	281
	282	histname = "fHistJetsPtArea_";
	283	histname += i;
a487deae	284	fHistJetsPtArea[i] = new TH2F(histname.Data(), histname.Data(), (Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5, 30, 0, fJetRadius * fJetRadius * TMath::Pi() * 3);
a487deae	285	fHistJetsPtArea[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
e44e8726	286	fHistJetsPtArea[i]->GetYaxis()->SetTitle("area");
	287	fOutput->Add(fHistJetsPtArea[i]);
	288	}
6fd5039f	289	}
6fd5039f	290
c3ba2d3d	291	PostData(1, fOutput); // Post data for ALL output slots >0 here, to get at least an empty histogram
	292	}
	293
	294	//________________________________________________________________________
6fd5039f	295	Bool_t AliAnalysisTaskSAQA::RetrieveEventObjects()
c3ba2d3d	296	{
c3604199	297	// Retrieve event objects.
	298
	299	if (!AliAnalysisTaskEmcalJet::RetrieveEventObjects())
6fd5039f	300	return kFALSE;
c3ba2d3d	301
7c1d624c	302	fNclusters = 0;
	303	fNtracks = 0;
	304	fNjets = 0;
	305
6fd5039f	306	return kTRUE;
c3ba2d3d	307	}
c3ba2d3d	308
7c1d624c	309
c3ba2d3d	310	//________________________________________________________________________
6fd5039f	311	Bool_t AliAnalysisTaskSAQA::FillHistograms()
c3ba2d3d	312	{
c3604199	313	// Fill histograms.
c3604199	314
a487deae	315	Float_t clusSum = 0;
a487deae	316	Float_t trackSum = 0;
090a0c3e	317
a487deae	318	if (fTracks) {
a487deae	319	trackSum = DoTrackLoop();
c3ba2d3d	320
a487deae	321	fHistTracksCent->Fill(fCent, fNtracks);
a487deae	322	}
be7b6e63	323
a487deae	324	if (fCaloClusters) {
a487deae	325	clusSum = DoClusterLoop();
c3ba2d3d	326
7c1d624c	327	fHistClusCent->Fill(fCent, fNclusters);
	328	fHistClusTracks->Fill(fNtracks, fNclusters);
	329
6fd5039f	330	Float_t cellSum = 0, cellCutSum = 0;
6fd5039f	331
be7b6e63	332	Int_t ncells = DoCellLoop(cellSum, cellCutSum);
	333
	334	if (fTracks)
a487deae	335	fHistCellsTracks->Fill(fNtracks, ncells);
be7b6e63	336
be7b6e63	337	fHistCellsCent->Fill(fCent, ncells);
6fd5039f	338
	339	fHistChVSneCells->Fill(cellSum, trackSum);
	340	fHistChVSneClus->Fill(clusSum, trackSum);
	341	fHistChVSneCorrCells->Fill(cellCutSum, trackSum);
a487deae	342	}
c3ba2d3d	343
a487deae	344	if (fJets) {
a487deae	345	DoJetLoop();
6fd5039f	346
a487deae	347	fHistJetsCent->Fill(fCent, fNjets);
a487deae	348	fHistJetsParts->Fill(fNtracks + fNclusters, fNjets);
6fd5039f	349	}
c3ba2d3d	350
6fd5039f	351	return kTRUE;
c3ba2d3d	352	}
	353
	354	//________________________________________________________________________
be7b6e63	355	Int_t AliAnalysisTaskSAQA::DoCellLoop(Float_t &sum, Float_t &sum_cut)
c3ba2d3d	356	{
c3604199	357	// Do cell loop.
c3604199	358
c3ba2d3d	359	AliVCaloCells *cells = InputEvent()->GetEMCALCells();
	360
	361	if (!cells)
be7b6e63	362	return 0;
c3ba2d3d	363
c3604199	364	const Int_t ncells = cells->GetNumberOfCells();
c3ba2d3d	365
c3ba2d3d	366	for (Int_t pos = 0; pos < ncells; pos++) {
f483218e	367	Float_t amp = cells->GetAmplitude(pos);
	368	Int_t absId = cells->GetCellNumber(pos);
	369	fHistCellsAbsIdEnergy->Fill(absId,amp);
c3ba2d3d	370	sum += amp;
c3ba2d3d	371	if (amp < fCellEnergyCut)
c3ba2d3d	372	continue;
c3ba2d3d	373	sum_cut += amp;
c3ba2d3d	374	}
be7b6e63	375
be7b6e63	376	return ncells;
c3ba2d3d	377	}
	378
	379	//________________________________________________________________________
	380	Float_t AliAnalysisTaskSAQA::DoClusterLoop()
	381	{
c3604199	382	// Do cluster loop.
c3604199	383
c3ba2d3d	384	if (!fCaloClusters)
	385	return 0;
	386
	387	Float_t sum = 0;
	388
	389	// Cluster loop
	390	Int_t nclusters = fCaloClusters->GetEntriesFast();
	391
	392	for (Int_t iClusters = 0; iClusters < nclusters; iClusters++) {
e44e8726	393	AliVCluster* cluster = static_cast<AliVCluster*>(fCaloClusters->At(iClusters));
c3ba2d3d	394	if (!cluster) {
	395	AliError(Form("Could not receive cluster %d", iClusters));
	396	continue;
	397	}
c3ba2d3d	398
e44e8726	399	if (!AcceptCluster(cluster, kTRUE))
e44e8726	400	continue;
c3ba2d3d	401
c3ba2d3d	402	sum += cluster->E();
c3ba2d3d	403
7cf4626b	404	TLorentzVector nPart;
7cf4626b	405	cluster->GetMomentum(nPart, fVertex);
c3ba2d3d	406
a487deae	407	fHistClusPhiEtaEnergy[fCentBin]->Fill(nPart.Eta(), nPart.Phi(), cluster->E());
7cf4626b	408	fHistNCellsEnergy->Fill(cluster->E(), cluster->GetNCells());
090a0c3e	409
090a0c3e	410	fHistClusTimeEnergy->Fill(cluster->E(), cluster->GetTOF());
7c1d624c	411
7c1d624c	412	fNclusters++;
c3604199	413	}
c3ba2d3d	414
	415	return sum;
	416	}
	417
	418	//________________________________________________________________________
	419	Float_t AliAnalysisTaskSAQA::DoTrackLoop()
	420	{
c3604199	421	// Do track loop.
c3604199	422
c3ba2d3d	423	if (!fTracks)
	424	return 0;
	425
	426	Float_t sum = 0;
	427
c3ba2d3d	428	Int_t ntracks = fTracks->GetEntriesFast();
	429	Int_t nclusters = 0;
	430	if (fCaloClusters)
	431	nclusters = fCaloClusters->GetEntriesFast();
	432
c3604199	433	for (Int_t i = 0; i < ntracks; i++) {
c3ba2d3d	434
e44e8726	435	AliVParticle* track = static_cast<AliVParticle*>(fTracks->At(i)); // pointer to reconstructed to track
1f6fff78	436
c3604199	437	if (!track) {
c3ba2d3d	438	AliError(Form("Could not retrieve track %d",i));
	439	continue;
	440	}
c3ba2d3d	441
1f6fff78	442	AliVTrack* vtrack = dynamic_cast<AliVTrack*>(track);
	443
	444	if (vtrack && !AcceptTrack(vtrack, kTRUE))
	445	continue;
7c1d624c	446
7c1d624c	447	fNtracks++;
c3ba2d3d	448
c3ba2d3d	449	sum += track->P();
a487deae	450
a487deae	451	fHistTrPhiEtaPt[fCentBin]->Fill(track->Eta(), track->Phi(), track->Pt());
c3ba2d3d	452
c3ba2d3d	453	Int_t label = track->GetLabel();
c3ba2d3d	454
	455	if (label >= 0 && label < 4) {
	456	fHistTrackEta[label]->Fill(track->Eta());
	457	fHistTrackPhi[label]->Fill(track->Phi());
	458	}
	459
a487deae	460	fHistTrackEta[4]->Fill(track->Eta());
a487deae	461	fHistTrackPhi[4]->Fill(track->Phi());
aa4d701c	462
1f6fff78	463	if (!vtrack)
c3ba2d3d	464	continue;
c3ba2d3d	465
079b4732	466	if (vtrack->GetTrackEtaOnEMCal() == -999 \|\| vtrack->GetTrackPhiOnEMCal() == -999)
	467	fHistTrPhiEtaNonProp->Fill(vtrack->Eta(), vtrack->Phi());
	468
1f6fff78	469	fHistTrEmcPhiEta->Fill(vtrack->GetTrackEtaOnEMCal(), vtrack->GetTrackPhiOnEMCal());
3e951e32	470	fHistDeltaEtaPt->Fill(vtrack->Pt(), vtrack->Eta() - vtrack->GetTrackEtaOnEMCal());
3e951e32	471	fHistDeltaPhiPt->Fill(vtrack->Pt(), vtrack->Phi() - vtrack->GetTrackPhiOnEMCal());
c3ba2d3d	472	}
	473
	474	return sum;
	475	}
	476
	477	//________________________________________________________________________
	478	void AliAnalysisTaskSAQA::DoJetLoop()
	479	{
c3604199	480	// Do jet loop.
c3604199	481
c3ba2d3d	482	if (!fJets)
	483	return;
	484
	485	Int_t njets = fJets->GetEntriesFast();
	486
	487	for (Int_t ij = 0; ij < njets; ij++) {
	488
e44e8726	489	AliEmcalJet* jet = static_cast<AliEmcalJet*>(fJets->At(ij));
c3ba2d3d	490
	491	if (!jet) {
	492	AliError(Form("Could not receive jet %d", ij));
	493	continue;
	494	}
	495
a487deae	496	if (!AcceptJet(jet))
c3ba2d3d	497	continue;
6fd5039f	498
7c1d624c	499	fNjets++;
7c1d624c	500
a487deae	501	fHistJetsPhiEtaPt[fCentBin]->Fill(jet->Eta(), jet->Phi(), jet->Pt());
99cfd012	502	fHistJetsPtArea[fCentBin]->Fill(jet->Pt(), jet->Area());
c3ba2d3d	503	}
	504	}
	505
c3ba2d3d	506
	507	//________________________________________________________________________
	508	void AliAnalysisTaskSAQA::Terminate(Option_t *)
	509	{
	510	// Called once at the end of the analysis.
	511	}