//
// Author: S.Aiola
-#include <TObject.h>
-#include <TChain.h>
#include <TClonesArray.h>
-#include <TH1F.h>
#include <TH2F.h>
+#include <TH3F.h>
+#include <THnSparse.h>
#include <TList.h>
#include <TLorentzVector.h>
-#include <TRandom3.h>
-#include <TParameter.h>
-#include "AliAnalysisManager.h"
-#include "AliCentrality.h"
#include "AliVCluster.h"
#include "AliVParticle.h"
-#include "AliVTrack.h"
#include "AliEmcalJet.h"
-#include "AliVEventHandler.h"
#include "AliRhoParameter.h"
#include "AliLog.h"
//________________________________________________________________________
AliAnalysisTaskSAJF::AliAnalysisTaskSAJF() :
AliAnalysisTaskEmcalJet("AliAnalysisTaskSAJF", kTRUE),
- fMCAna(kFALSE),
- fMinRC2LJ(-1),
- fEmbJetsName("EmbJets"),
- fEmbTracksName(""),
- fEmbCaloName(""),
- fRandTracksName("TracksRandomized"),
- fRandCaloName("CaloClustersRandomized"),
- fRhoName("Rho"),
- fRCperEvent(-1),
- fEmbJets(0),
- fEmbTracks(0),
- fEmbCaloClusters(0),
- fRandTracks(0),
- fRandCaloClusters(0),
- fRho(0),
- fRhoVal(0),
- fEmbeddedClusterId(-1),
- fEmbeddedTrackId(-1),
- fHistCentrality(0),
- fHistDeltaVectorPt(0),
- fHistRhoVSleadJetPt(0),
- fHistRCPhiEta(0),
- fHistRCPtExLJVSDPhiLJ(0),
- fHistRhoVSRCPt(0),
- fHistEmbJetPhiEta(0),
- fHistEmbPartPhiEta(0),
- fHistRhoVSEmbBkg(0)
+ fHistoType(1),
+ fHistJetObservables(0)
{
// Default constructor.
for (Int_t i = 0; i < 4; i++) {
- fHistEvents[i] = 0;
- fHistTracksPt[i] = 0;
- fHistClustersPt[i] = 0;
- fHistJetPhiEta[i] = 0;
- fHistJetsPt[i] = 0;
- fHistJetsPtArea[i] = 0;
- fHistLeadingJetPt[i] = 0;
- fHist2LeadingJetPt[i] = 0;
- fHistJetsNEFvsPt[i] = 0;
- fHistJetsZvsPt[i] = 0;
- fHistMaxTrackPtvsJetPt[i] = 0;
- fHistMaxClusPtvsJetPt[i] = 0;
- fHistMaxPartPtvsJetPt[i] = 0;
- fHistMaxTrackPtvsJetCorrPt[i] = 0;
- fHistMaxClusPtvsJetCorrPt[i] = 0;
- fHistMaxPartPtvsJetCorrPt[i] = 0;
- fHistConstituents[i] = 0;
- fHistRho[i] = 0;
- fHistJetsCorrPt[i] = 0;
- fHistJetsCorrPtArea[i] = 0;
- fHistLeadingJetCorrPt[i] = 0;
- fHistRCPtRigid[i] = 0;
- fHistRCPt[i] = 0;
- fHistRCPtExLJ[i] = 0;
- fHistRCPtRand[i] = 0;
- fHistDeltaPtRCRigid[i] = 0;
- fHistDeltaPtRC[i] = 0;
- fHistDeltaPtRCExLJ[i] = 0;
- fHistDeltaPtRCRand[i] = 0;
- fHistEmbNotFoundPhiEta[i] = 0;
- fHistEmbJetsPt[i] = 0;
- fHistEmbJetsCorrPt[i] = 0;
- fHistEmbPartPt[i] = 0;
- fHistDistEmbPartJetAxis[i] = 0;
- fHistDeltaPtEmb[i] = 0;
+ fHistTracksJetPt[i] = 0;
+ fHistClustersJetPt[i] = 0;
+ fHistTracksPtDist[i] = 0;
+ fHistClustersPtDist[i] = 0;
+
+ fHistJetPtEtaPhi[i] = 0;
+ fHistJetPtArea[i] = 0;
+ fHistJetPtEP[i] = 0;
+ fHistJetPtNEF[i] = 0;
+ fHistJetPtZ[i] = 0;
+ fHistJetPtLeadingPartPt[i] = 0;
+ fHistJetCorrPtEtaPhi[i] = 0;
+ fHistJetCorrPtArea[i] = 0;
+ fHistJetCorrPtEP[i] = 0;
+ fHistJetCorrPtNEF[i] = 0;
+ fHistJetCorrPtZ[i] = 0;
+ fHistJetCorrPtLeadingPartPt[i] = 0;
+ fHistJetPtCorrPt[i] = 0;
+ fHistJetPtMCPt[i] = 0;
+ fHistJetMCPtCorrPt[i] = 0;
}
+
+ SetMakeGeneralHistograms(kTRUE);
}
//________________________________________________________________________
AliAnalysisTaskSAJF::AliAnalysisTaskSAJF(const char *name) :
AliAnalysisTaskEmcalJet(name, kTRUE),
- fMCAna(kFALSE),
- fMinRC2LJ(-1),
- fEmbJetsName("EmbJets"),
- fEmbTracksName(""),
- fEmbCaloName(""),
- fRandTracksName("TracksRandomized"),
- fRandCaloName("CaloClustersRandomized"),
- fRhoName("Rho"),
- fRCperEvent(-1),
- fEmbJets(0),
- fEmbTracks(0),
- fEmbCaloClusters(0),
- fRandTracks(0),
- fRandCaloClusters(0),
- fRho(0),
- fRhoVal(0),
- fEmbeddedClusterId(-1),
- fEmbeddedTrackId(-1),
- fHistCentrality(0),
- fHistDeltaVectorPt(0),
- fHistRhoVSleadJetPt(0),
- fHistRCPhiEta(0),
- fHistRCPtExLJVSDPhiLJ(0),
- fHistRhoVSRCPt(0),
- fHistEmbJetPhiEta(0),
- fHistEmbPartPhiEta(0),
- fHistRhoVSEmbBkg(0)
+ fHistoType(1),
+ fHistJetObservables(0)
{
// Standard constructor.
for (Int_t i = 0; i < 4; i++) {
- fHistEvents[i] = 0;
- fHistTracksPt[i] = 0;
- fHistClustersPt[i] = 0;
- fHistJetPhiEta[i] = 0;
- fHistJetsPt[i] = 0;
- fHistJetsPtArea[i] = 0;
- fHistLeadingJetPt[i] = 0;
- fHist2LeadingJetPt[i] = 0;
- fHistJetsNEFvsPt[i] = 0;
- fHistJetsZvsPt[i] = 0;
- fHistMaxTrackPtvsJetPt[i] = 0;
- fHistMaxClusPtvsJetPt[i] = 0;
- fHistMaxPartPtvsJetPt[i] = 0;
- fHistMaxTrackPtvsJetCorrPt[i] = 0;
- fHistMaxClusPtvsJetCorrPt[i] = 0;
- fHistMaxPartPtvsJetCorrPt[i] = 0;
- fHistConstituents[i] = 0;
- fHistRho[i] = 0;
- fHistJetsCorrPt[i] = 0;
- fHistJetsCorrPtArea[i] = 0;
- fHistLeadingJetCorrPt[i] = 0;
- fHistRCPtRigid[i] = 0;
- fHistRCPt[i] = 0;
- fHistRCPtExLJ[i] = 0;
- fHistRCPtRand[i] = 0;
- fHistDeltaPtRCRigid[i] = 0;
- fHistDeltaPtRC[i] = 0;
- fHistDeltaPtRCExLJ[i] = 0;
- fHistDeltaPtRCRand[i] = 0;
- fHistEmbNotFoundPhiEta[i] = 0;
- fHistEmbJetsPt[i] = 0;
- fHistEmbJetsCorrPt[i] = 0;
- fHistEmbPartPt[i] = 0;
- fHistDistEmbPartJetAxis[i] = 0;
- fHistDeltaPtEmb[i] = 0;
+ fHistTracksJetPt[i] = 0;
+ fHistClustersJetPt[i] = 0;
+ fHistTracksPtDist[i] = 0;
+ fHistClustersPtDist[i] = 0;
+
+ fHistJetPtEtaPhi[i] = 0;
+ fHistJetPtArea[i] = 0;
+ fHistJetPtEP[i] = 0;
+ fHistJetPtNEF[i] = 0;
+ fHistJetPtZ[i] = 0;
+ fHistJetPtLeadingPartPt[i] = 0;
+ fHistJetCorrPtEtaPhi[i] = 0;
+ fHistJetCorrPtArea[i] = 0;
+ fHistJetCorrPtEP[i] = 0;
+ fHistJetCorrPtNEF[i] = 0;
+ fHistJetCorrPtZ[i] = 0;
+ fHistJetCorrPtLeadingPartPt[i] = 0;
+ fHistJetPtCorrPt[i] = 0;
+ fHistJetPtMCPt[i] = 0;
+ fHistJetMCPtCorrPt[i] = 0;
}
-}
-//________________________________________________________________________
-AliAnalysisTaskSAJF::~AliAnalysisTaskSAJF()
-{
- // Destructor.
+ SetMakeGeneralHistograms(kTRUE);
}
//________________________________________________________________________
-void AliAnalysisTaskSAJF::UserCreateOutputObjects()
+void AliAnalysisTaskSAJF::AllocateTHnSparse()
{
- // Create user output.
-
- AliAnalysisTaskEmcalJet::UserCreateOutputObjects();
-
- OpenFile(1);
- fOutput = new TList();
- fOutput->SetOwner();
-
- fHistCentrality = new TH1F("fHistCentrality","fHistCentrality", fNbins, 0, 100);
- fHistCentrality->GetXaxis()->SetTitle("Centrality (%)");
- fHistCentrality->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistCentrality);
-
- fHistDeltaVectorPt = new TH1F("fHistDeltaVectorPt", "fHistDeltaVectorPt", fNbins, -50, 50);
- fHistDeltaVectorPt->GetXaxis()->SetTitle("#deltap_{T} [GeV/c]");
- fHistDeltaVectorPt->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistDeltaVectorPt);
-
- fHistRhoVSleadJetPt = new TH2F("fHistRhoVSleadJetPt","fHistRhoVSleadJetPt", fNbins, fMinBinPt, fMaxBinPt, fNbins, fMinBinPt, fMaxBinPt);
- fHistRhoVSleadJetPt->GetXaxis()->SetTitle("#rho * area [GeV/c]");
- fHistRhoVSleadJetPt->GetYaxis()->SetTitle("Leading jet p_{T} [GeV/c]");
- fOutput->Add(fHistRhoVSleadJetPt);
-
- fHistRCPhiEta = new TH2F("fHistRCPhiEta","Phi-Eta distribution of rigid cones", 50, -1, 1, 101, 0, TMath::Pi() * 2.02);
- fHistRCPhiEta->GetXaxis()->SetTitle("#eta");
- fHistRCPhiEta->GetYaxis()->SetTitle("#phi");
- fOutput->Add(fHistRCPhiEta);
-
- fHistRCPtExLJVSDPhiLJ = new TH2F("fHistRCPtExLJVSDPhiLJ","fHistRCPtExLJVSDPhiLJ", fNbins, fMinBinPt, fMaxBinPt, 128, -1.6, 4.8);
- fHistRCPtExLJVSDPhiLJ->GetXaxis()->SetTitle("rigid cone p_{T} [GeV/c]");
- fHistRCPtExLJVSDPhiLJ->GetYaxis()->SetTitle("#Delta#phi");
- fOutput->Add(fHistRCPtExLJVSDPhiLJ);
-
- fHistRhoVSRCPt = new TH2F("fHistRhoVSRCPt","fHistRhoVSRCPt", fNbins, fMinBinPt, fMaxBinPt, fNbins, fMinBinPt, fMaxBinPt);
- fHistRhoVSRCPt->GetXaxis()->SetTitle("A#rho [GeV/c]");
- fHistRhoVSRCPt->GetYaxis()->SetTitle("rigid cone p_{T} [GeV/c]");
- fOutput->Add(fHistRhoVSRCPt);
-
- if (!fEmbJetsName.IsNull()) {
- fHistEmbJetPhiEta = new TH2F("fHistEmbJetPhiEta","Phi-Eta distribution of embedded jets", 50, -1, 1, 101, 0, TMath::Pi() * 2.02);
- fHistEmbJetPhiEta->GetXaxis()->SetTitle("#eta");
- fHistEmbJetPhiEta->GetYaxis()->SetTitle("#phi");
- fOutput->Add(fHistEmbJetPhiEta);
+ 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++;
- fHistEmbPartPhiEta = new TH2F("fHistEmbPartPhiEta","Phi-Eta distribution of embedded particles", 50, -1, 1, 101, 0, TMath::Pi() * 2.02);
- fHistEmbPartPhiEta->GetXaxis()->SetTitle("#eta");
- fHistEmbPartPhiEta->GetYaxis()->SetTitle("#phi");
- fOutput->Add(fHistEmbPartPhiEta);
-
- fHistRhoVSEmbBkg = new TH2F("fHistRhoVSEmbBkg","fHistRhoVSEmbBkg", fNbins, fMinBinPt, fMaxBinPt, fNbins, fMinBinPt, fMaxBinPt);
- fHistRhoVSEmbBkg->GetXaxis()->SetTitle("A#rho [GeV/c]");
- fHistRhoVSEmbBkg->GetYaxis()->SetTitle("background of embedded track [GeV/c]");
- fOutput->Add(fHistRhoVSEmbBkg);
+ title[dim] = "#phi_{jet} - #psi_{RP}";
+ nbins[dim] = 100;
+ min[dim] = 0;
+ max[dim] = TMath::Pi();
+ dim++;
}
- TString histname;
+ title[dim] = "#eta";
+ 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++;
+ }
- for (Int_t i = 0; i < 4; i++) {
- histname = "fHistEvents_";
- histname += i;
- fHistEvents[i] = new TH1F(histname,histname, 6, 0, 6);
- fHistEvents[i]->GetXaxis()->SetTitle("Event state");
- fHistEvents[i]->GetYaxis()->SetTitle("counts");
- fHistEvents[i]->GetXaxis()->SetBinLabel(1, "No jets");
- fHistEvents[i]->GetXaxis()->SetBinLabel(2, "Max jet not found");
- fHistEvents[i]->GetXaxis()->SetBinLabel(3, "Rho == 0");
- fHistEvents[i]->GetXaxis()->SetBinLabel(4, "Max jet <= 0");
- fHistEvents[i]->GetXaxis()->SetBinLabel(5, "OK");
- fOutput->Add(fHistEvents[i]);
-
- if (fAnaType != kEMCALOnly) {
- histname = "fHistTracksPt_";
- histname += i;
- fHistTracksPt[i] = new TH1F(histname.Data(), histname.Data(), (Int_t)(fNbins / 2.5), fMinBinPt, fMaxBinPt / 2.5);
- fHistTracksPt[i]->GetXaxis()->SetTitle("p_{T} [GeV/c]");
- fHistTracksPt[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistTracksPt[i]);
- }
+ if (!GetRhoName().IsNull()) {
+ title[dim] = "p_{T}^{corr} (GeV/c)";
+ nbins[dim] = fNbins*2;
+ min[dim] = -fMaxBinPt;
+ max[dim] = fMaxBinPt;
+ dim++;
+ }
- if (fAnaType == kEMCAL || fAnaType == kEMCALOnly) {
- histname = "fHistClustersPt_";
- histname += i;
- fHistClustersPt[i] = new TH1F(histname.Data(), histname.Data(), (Int_t)(fNbins / 2.5), fMinBinPt, fMaxBinPt / 2.5);
- fHistClustersPt[i]->GetXaxis()->SetTitle("p_{T} [GeV/c]");
- fHistClustersPt[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistClustersPt[i]);
- }
+ 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]);
+}
- histname = "fHistJetPhiEta_";
- histname += i;
- fHistJetPhiEta[i] = new TH2F(histname.Data(), histname.Data(), 50, -1, 1, 101, 0, TMath::Pi() * 2.02);
- fHistJetPhiEta[i]->GetXaxis()->SetTitle("#eta");
- fHistJetPhiEta[i]->GetYaxis()->SetTitle("#phi");
- fOutput->Add(fHistJetPhiEta[i]);
+//________________________________________________________________________
+void AliAnalysisTaskSAJF::AllocateTHX()
+{
+ for (Int_t i = 0; i < 4; i++) {
+ TString histname;
- histname = "fHistJetsPt_";
+ histname = "fHistJetPtEtaPhi_";
histname += i;
- fHistJetsPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt);
- fHistJetsPt[i]->GetXaxis()->SetTitle("p_{T}^{raw} [GeV/c]");
- fHistJetsPt[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistJetsPt[i]);
-
- histname = "fHistJetsPtArea_";
+ 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;
- fHistJetsPtArea[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, 40, 0, fJetRadius * fJetRadius * TMath::Pi() * 3);
- fHistJetsPtArea[i]->GetXaxis()->SetTitle("p_{T}^{raw} [GeV/c]");
- fHistJetsPtArea[i]->GetYaxis()->SetTitle("area");
- fOutput->Add(fHistJetsPtArea[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 = "fHistLeadingJetPt_";
+ histname = "fHistJetPtEP_";
histname += i;
- fHistLeadingJetPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt);
- fHistLeadingJetPt[i]->GetXaxis()->SetTitle("p_{T}^{raw} [GeV]");
- fHistLeadingJetPt[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistLeadingJetPt[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 = "fHist2LeadingJetPt_";
+ histname = "fHistJetPtNEF_";
histname += i;
- fHist2LeadingJetPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt);
- fHist2LeadingJetPt[i]->GetXaxis()->SetTitle("p_{T}^{raw} [GeV]");
- fHist2LeadingJetPt[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHist2LeadingJetPt[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 = "fHistJetsZvsPt_";
+ histname = "fHistJetPtZ_";
histname += i;
- fHistJetsZvsPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins, 0, 1.2, fNbins, fMinBinPt, fMaxBinPt);
- fHistJetsZvsPt[i]->GetXaxis()->SetTitle("Z");
- fHistJetsZvsPt[i]->GetYaxis()->SetTitle("p_{T}^{raw} [GeV/c]");
- fOutput->Add(fHistJetsZvsPt[i]);
-
- if (fAnaType == kEMCAL || fAnaType == kEMCALOnly) {
- histname = "fHistJetsNEFvsPt_";
- histname += i;
- fHistJetsNEFvsPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins, 0, 1.2, fNbins, fMinBinPt, fMaxBinPt);
- fHistJetsNEFvsPt[i]->GetXaxis()->SetTitle("NEF");
- fHistJetsNEFvsPt[i]->GetYaxis()->SetTitle("p_{T}^{raw} [GeV/c]");
- fOutput->Add(fHistJetsNEFvsPt[i]);
- }
-
- if (fAnaType != kEMCALOnly) {
- histname = "fHistMaxTrackPtvsJetPt_";
- histname += i;
- fHistMaxTrackPtvsJetPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, (Int_t)(fNbins / 2.5), fMinBinPt, fMaxBinPt / 2.5);
- fHistMaxTrackPtvsJetPt[i]->GetXaxis()->SetTitle("p_{T}^{jet} [GeV/c]");
- fHistMaxTrackPtvsJetPt[i]->GetYaxis()->SetTitle("p_{T}^{track} [GeV/c]");
- fOutput->Add(fHistMaxTrackPtvsJetPt[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]);
- if (fAnaType == kEMCAL || fAnaType == kEMCALOnly) {
- histname = "fHistMaxClusPtvsJetPt_";
- histname += i;
- fHistMaxClusPtvsJetPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, (Int_t)(fNbins / 2.5), fMinBinPt, fMaxBinPt / 2.5);
- fHistMaxClusPtvsJetPt[i]->GetXaxis()->SetTitle("p_{T}^{jet} [GeV/c]");
- fHistMaxClusPtvsJetPt[i]->GetYaxis()->SetTitle("p_{T}^{clus} [GeV/c]");
- fOutput->Add(fHistMaxClusPtvsJetPt[i]);
- }
-
- histname = "fHistMaxPartPtvsJetPt_";
+ histname = "fHistJetPtLeadingPartPt_";
histname += i;
- fHistMaxPartPtvsJetPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, (Int_t)(fNbins / 2.5), fMinBinPt, fMaxBinPt / 2.5);
- fHistMaxPartPtvsJetPt[i]->GetXaxis()->SetTitle("p_{T}^{jet} [GeV/c]");
- fHistMaxPartPtvsJetPt[i]->GetYaxis()->SetTitle("p_{T}^{part} [GeV/c]");
- fOutput->Add(fHistMaxPartPtvsJetPt[i]);
-
- if (fAnaType != kEMCALOnly) {
- histname = "fHistMaxTrackPtvsJetCorrPt_";
+ 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;
- fHistMaxTrackPtvsJetCorrPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt, (Int_t)(fNbins / 2.5), fMinBinPt, fMaxBinPt / 2.5);
- fHistMaxTrackPtvsJetCorrPt[i]->GetXaxis()->SetTitle("p_{T}^{jet} [GeV/c]");
- fHistMaxTrackPtvsJetCorrPt[i]->GetYaxis()->SetTitle("p_{T}^{track} [GeV/c]");
- fOutput->Add(fHistMaxTrackPtvsJetCorrPt[i]);
- }
-
- if (fAnaType == kEMCAL || fAnaType == kEMCALOnly) {
- histname = "fHistMaxClusPtvsJetCorrPt_";
+ 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;
- fHistMaxClusPtvsJetCorrPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt, (Int_t)(fNbins / 2.5), fMinBinPt, fMaxBinPt / 2.5);
- fHistMaxClusPtvsJetCorrPt[i]->GetXaxis()->SetTitle("p_{T}^{jet} [GeV/c]");
- fHistMaxClusPtvsJetCorrPt[i]->GetYaxis()->SetTitle("p_{T}^{clus} [GeV/c]");
- fOutput->Add(fHistMaxClusPtvsJetCorrPt[i]);
- }
-
- histname = "fHistMaxPartPtvsJetCorrPt_";
- histname += i;
- fHistMaxPartPtvsJetCorrPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt, (Int_t)(fNbins / 2.5), fMinBinPt, fMaxBinPt / 2.5);
- fHistMaxPartPtvsJetCorrPt[i]->GetXaxis()->SetTitle("p_{T}^{jet} [GeV/c]");
- fHistMaxPartPtvsJetCorrPt[i]->GetYaxis()->SetTitle("p_{T}^{part} [GeV/c]");
- fOutput->Add(fHistMaxPartPtvsJetCorrPt[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 = "fHistConstituents_";
- histname += i;
- fHistConstituents[i] = new TH2F(histname.Data(), histname.Data(), 100, 1, 101, 100, -0.5, 99.5);
- fHistConstituents[i]->GetXaxis()->SetTitle("p_{T,part} [GeV/c]");
- fHistConstituents[i]->GetYaxis()->SetTitle("no. of particles");
- fOutput->Add(fHistConstituents[i]);
-
- histname = "fHistRho_";
- histname += i;
- fHistRho[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt * 2);
- fHistRho[i]->GetXaxis()->SetTitle("p_{T} [GeV/c]");
- fHistRho[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistRho[i]);
-
- histname = "fHistJetsCorrPt_";
- histname += i;
- fHistJetsCorrPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
- fHistJetsCorrPt[i]->GetXaxis()->SetTitle("p_{T}^{corr} [GeV/c]");
- fHistJetsCorrPt[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistJetsCorrPt[i]);
-
- histname = "fHistJetsCorrPtArea_";
- histname += i;
- fHistJetsCorrPtArea[i] = new TH2F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt, 40, 0, fJetRadius * fJetRadius * TMath::Pi() * 3);
- fHistJetsCorrPtArea[i]->GetXaxis()->SetTitle("p_{T}^{corr} [GeV/c]");
- fHistJetsCorrPtArea[i]->GetYaxis()->SetTitle("area");
- fOutput->Add(fHistJetsCorrPtArea[i]);
-
- histname = "fHistLeadingJetCorrPt_";
- histname += i;
- fHistLeadingJetCorrPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
- fHistLeadingJetCorrPt[i]->GetXaxis()->SetTitle("p_{T}^{RC} [GeV/c]");
- fHistLeadingJetCorrPt[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistLeadingJetCorrPt[i]);
-
- histname = "fHistRCPtRigid_";
- histname += i;
- fHistRCPtRigid[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt * 2);
- fHistRCPtRigid[i]->GetXaxis()->SetTitle("rigid cone p_{T} [GeV/c]");
- fHistRCPtRigid[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistRCPtRigid[i]);
-
- histname = "fHistRCPt_";
- histname += i;
- fHistRCPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt * 2);
- fHistRCPt[i]->GetXaxis()->SetTitle("rigid cone p_{T} [GeV/c]");
- fHistRCPt[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistRCPt[i]);
-
- histname = "fHistRCPtExLJ_";
- histname += i;
- fHistRCPtExLJ[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt * 2);
- fHistRCPtExLJ[i]->GetXaxis()->SetTitle("rigid cone p_{T}^{RC} [GeV/c]");
- fHistRCPtExLJ[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistRCPtExLJ[i]);
-
- histname = "fHistRCPtRand_";
- histname += i;
- fHistRCPtRand[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt * 2);
- fHistRCPtRand[i]->GetXaxis()->SetTitle("rigid cone p_{T}^{RC} [GeV/c]");
- fHistRCPtRand[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistRCPtRand[i]);
-
- histname = "fHistDeltaPtRCRigid_";
- histname += i;
- fHistDeltaPtRCRigid[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
- fHistDeltaPtRCRigid[i]->GetXaxis()->SetTitle("#deltap_{T}^{RC} [GeV/c]");
- fHistDeltaPtRCRigid[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistDeltaPtRCRigid[i]);
-
- histname = "fHistDeltaPtRC_";
- histname += i;
- fHistDeltaPtRC[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
- fHistDeltaPtRC[i]->GetXaxis()->SetTitle("#deltap_{T}^{RC} [GeV/c]");
- fHistDeltaPtRC[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistDeltaPtRC[i]);
-
- histname = "fHistDeltaPtRCExLJ_";
- histname += i;
- fHistDeltaPtRCExLJ[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
- fHistDeltaPtRCExLJ[i]->GetXaxis()->SetTitle("#deltap_{T}^{RC} [GeV/c]");
- fHistDeltaPtRCExLJ[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistDeltaPtRCExLJ[i]);
-
- histname = "fHistDeltaPtRCRand_";
- histname += i;
- fHistDeltaPtRCRand[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
- fHistDeltaPtRCRand[i]->GetXaxis()->SetTitle("#deltap_{T}^{RC} [GeV/c]");
- fHistDeltaPtRCRand[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistDeltaPtRCRand[i]);
-
- if (!fEmbJetsName.IsNull()) {
- histname = "fHistEmbJetsPt_";
+ histname = "fHistJetCorrPtEP_";
histname += i;
- fHistEmbJetsPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt);
- fHistEmbJetsPt[i]->GetXaxis()->SetTitle("embedded jet p_{T}^{raw} [GeV/c]");
- fHistEmbJetsPt[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistEmbJetsPt[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 = "fHistEmbJetsCorrPt_";
+ histname = "fHistJetCorrPtNEF_";
histname += i;
- fHistEmbJetsCorrPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
- fHistEmbJetsCorrPt[i]->GetXaxis()->SetTitle("embedded jet p_{T}^{corr} [GeV/c]");
- fHistEmbJetsCorrPt[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistEmbJetsCorrPt[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 = "fHistEmbJetsArea_";
+ histname = "fHistJetCorrPtZ_";
histname += i;
- fHistEmbJetsArea[i] = new TH1F(histname.Data(), histname.Data(), 40, 0, fJetRadius * fJetRadius * TMath::Pi() * 3);
- fHistEmbJetsArea[i]->GetXaxis()->SetTitle("area");
- fHistEmbJetsArea[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistEmbJetsArea[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 = "fHistEmbPartPt_";
+ histname = "fHistJetCorrPtLeadingPartPt_";
histname += i;
- fHistEmbPartPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt);
- fHistEmbPartPt[i]->GetXaxis()->SetTitle("embedded particle p_{T}^{emb} [GeV/c]");
- fHistEmbPartPt[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistEmbPartPt[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 = "fHistDistEmbPartJetAxis_";
+ histname = "fHistJetPtCorrPt_";
histname += i;
- fHistDistEmbPartJetAxis[i] = new TH1F(histname.Data(), histname.Data(), 50, 0, 0.5);
- fHistDistEmbPartJetAxis[i]->GetXaxis()->SetTitle("distance");
- fHistDistEmbPartJetAxis[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistDistEmbPartJetAxis[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]);
+ }
+ }
- histname = "fHistEmbNotFoundPhiEta_";
- histname += i;
- fHistEmbNotFoundPhiEta[i] = new TH2F(histname.Data(), histname.Data(), 40, -2, 2, 64, 0, 6.4);
- fHistEmbNotFoundPhiEta[i]->GetXaxis()->SetTitle("#eta");
- fHistEmbNotFoundPhiEta[i]->GetYaxis()->SetTitle("#phi");
- fOutput->Add(fHistEmbNotFoundPhiEta[i]);
-
- histname = "fHistDeltaPtEmb_";
+ if (fIsEmbedded) {
+ histname = "fHistJetPtMCPt_";
histname += i;
- fHistDeltaPtEmb[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
- fHistDeltaPtEmb[i]->GetXaxis()->SetTitle("#deltap_{T}^{emb} [GeV/c]");
- fHistDeltaPtEmb[i]->GetYaxis()->SetTitle("counts");
- fOutput->Add(fHistDeltaPtEmb[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]);
}
}
-
- PostData(1, fOutput); // Post data for ALL output slots >0 here, to get at least an empty histogram
-}
-
-//________________________________________________________________________
-Bool_t AliAnalysisTaskSAJF::RetrieveEventObjects()
-{
- // Retrieve event objects.
-
- if (!AliAnalysisTaskEmcalJet::RetrieveEventObjects())
- return kFALSE;
-
- if (fRho)
- fRhoVal = fRho->GetVal();
-
- return kTRUE;
}
//________________________________________________________________________
-Bool_t AliAnalysisTaskSAJF::FillHistograms()
+void AliAnalysisTaskSAJF::UserCreateOutputObjects()
{
- // Fill histograms.
-
- Int_t maxJetIndex = -1;
- Int_t max2JetIndex = -1;
-
- GetLeadingJets(maxJetIndex, max2JetIndex);
-
- if (maxJetIndex < 0) { // no accepted jet, skipping
- fHistEvents[fCentBin]->Fill("No jets", 1);
- return kTRUE;
- }
-
- AliEmcalJet* jet = static_cast<AliEmcalJet*>(fJets->At(maxJetIndex));
- if (!jet) { // error, I cannot get the lead jet from collection (should never happen), skipping
- fHistEvents[fCentBin]->Fill("Max jet not found", 1);
- return kTRUE;
- }
-
- // OK, event accepted
-
- if (fRhoVal == 0)
- fHistEvents[fCentBin]->Fill("Rho == 0", 1);
-
- Float_t maxJetCorrPt = jet->Pt() - fRhoVal * jet->Area();
-
- if (maxJetCorrPt <= 0)
- fHistEvents[fCentBin]->Fill("Max jet <= 0", 1);
-
- fHistEvents[fCentBin]->Fill("OK", 1);
+ // Create user output.
- // ************
- // General histograms
- // _________________________________
+ AliAnalysisTaskEmcalJet::UserCreateOutputObjects();
- DoJetLoop();
- DoTrackLoop();
- DoClusterLoop();
+ if (fHistoType == 0)
+ AllocateTHX();
+ else
+ AllocateTHnSparse();
- fHistCentrality->Fill(fCent);
- fHistRho[fCentBin]->Fill(fRhoVal);
+ for (Int_t i = 0; i < 4; i++) {
+ TString histname;
- if (jet) {
- fHistLeadingJetPt[fCentBin]->Fill(jet->Pt());
- fHistRhoVSleadJetPt->Fill(fRhoVal * jet->Area(), jet->Pt());
- fHistLeadingJetCorrPt[fCentBin]->Fill(maxJetCorrPt);
- }
-
- AliEmcalJet* jet2 = 0;
- if (max2JetIndex >= 0)
- jet2 = static_cast<AliEmcalJet*>(fJets->At(max2JetIndex));
-
- if (jet2)
- fHist2LeadingJetPt[fCentBin]->Fill(jet2->Pt());
-
- // ************
- // Random cones
- // _________________________________
-
- const Float_t rcArea = fJetRadius * fJetRadius * TMath::Pi();
-
- // Simple random cones
- for (Int_t i = 0; i < fRCperEvent; i++) {
- Float_t RCpt = 0;
- Float_t RCptRigid = 0;
- Float_t RCeta = 0;
- Float_t RCphi = 0;
- GetRigidCone(RCpt, RCptRigid, RCeta, RCphi, 0);
- if (RCpt > 0) {
- fHistRCPt[fCentBin]->Fill(RCpt / rcArea);
- fHistDeltaPtRC[fCentBin]->Fill(RCpt - rcArea * fRhoVal);
- }
- if (RCptRigid > 0) {
- fHistRCPtRigid[fCentBin]->Fill(RCptRigid / rcArea);
- fHistDeltaPtRCRigid[fCentBin]->Fill(RCptRigid - rcArea * fRhoVal);
- }
-
- // Random cones far from leading jet
- RCpt = 0;
- RCptRigid = 0;
- RCeta = 0;
- RCphi = 0;
- GetRigidCone(RCpt, RCptRigid, RCeta, RCphi, jet);
- if (RCpt > 0) {
- fHistRCPhiEta->Fill(RCeta, RCphi);
- fHistRhoVSRCPt->Fill(fRhoVal, RCpt / rcArea);
+ 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]);
- Float_t dphi = RCphi - jet->Phi();
- if (dphi > 4.8) dphi -= TMath::Pi() * 2;
- if (dphi < -1.6) dphi += TMath::Pi() * 2;
- fHistRCPtExLJVSDPhiLJ->Fill(RCpt, dphi);
-
- fHistRCPtExLJ[fCentBin]->Fill(RCpt / rcArea);
- fHistDeltaPtRCExLJ[fCentBin]->Fill(RCpt - rcArea * fRhoVal);
- }
-
- // Random cones with randomized particles
- RCpt = 0;
- RCptRigid = 0;
- RCeta = 0;
- RCphi = 0;
- GetRigidCone(RCpt, RCptRigid, RCeta, RCphi, 0, fRandTracks, fRandCaloClusters);
- if (RCpt > 0) {
- fHistRCPtRand[fCentBin]->Fill(RCpt / rcArea);
- fHistDeltaPtRCRand[fCentBin]->Fill(RCpt - rcArea * fRhoVal);
- }
- }
-
- // ************
- // Embedding
- // _________________________________
-
- if (!fEmbJets)
- return kTRUE;
-
- AliEmcalJet *embJet = 0;
- TObject *embPart = 0;
-
- DoEmbJetLoop(embJet, embPart);
-
- if (embJet) {
- Double_t probePt = 0;
- Double_t probeEta = 0;
- Double_t probePhi = 0;
-
- AliVCluster *cluster = dynamic_cast<AliVCluster*>(embPart);
- if (cluster) {
- TLorentzVector nPart;
- cluster->GetMomentum(nPart, fVertex);
-
- probeEta = nPart.Eta();
- probePhi = nPart.Phi();
- probePt = nPart.Pt();
- }
- else {
- AliVTrack *track = dynamic_cast<AliVTrack*>(embPart);
- if (track) {
- probeEta = track->Eta();
- probePhi = track->Phi();
- probePt = track->Pt();
- }
- else {
- AliWarning(Form("%s - Embedded jet found but embedded particle not found (wrong type?)!", GetName()));
- return kTRUE;
- }
- }
- Double_t distProbeJet = TMath::Sqrt((embJet->Eta() - probeEta) * (embJet->Eta() - probeEta) + (embJet->Phi() - probePhi) * (embJet->Phi() - probePhi));
-
- fHistEmbPartPt[fCentBin]->Fill(probePt);
- fHistEmbPartPhiEta->Fill(probeEta, probePhi);
-
- fHistDistEmbPartJetAxis[fCentBin]->Fill(distProbeJet);
-
- fHistEmbJetsPt[fCentBin]->Fill(embJet->Pt());
- fHistEmbJetsCorrPt[fCentBin]->Fill(embJet->Pt() - fRhoVal * embJet->Area());
- fHistEmbJetsArea[fCentBin]->Fill(embJet->Area());
- fHistEmbJetPhiEta->Fill(embJet->Eta(), embJet->Phi());
-
- fHistDeltaPtEmb[fCentBin]->Fill(embJet->Pt() - embJet->Area() * fRhoVal - probePt);
- fHistRhoVSEmbBkg->Fill(embJet->Area() * fRhoVal, embJet->Pt() - probePt);
- }
- else {
- if (fAnaType != kEMCALOnly)
- DoEmbTrackLoop();
- if (fAnaType == kEMCAL || fAnaType == kEMCALOnly)
- DoEmbClusterLoop();
- if (fEmbeddedTrackId >= 0) {
- AliVTrack *track2 = static_cast<AliVTrack*>(fEmbTracks->At(fEmbeddedTrackId));
- fHistEmbNotFoundPhiEta[fCentBin]->Fill(track2->Eta(), track2->Phi());
- }
- else if (fEmbeddedClusterId >= 0) {
- AliVCluster *cluster2 = static_cast<AliVCluster*>(fEmbCaloClusters->At(fEmbeddedClusterId));
- TLorentzVector nPart;
- cluster2->GetMomentum(nPart, fVertex);
- fHistEmbNotFoundPhiEta[fCentBin]->Fill(nPart.Eta(), nPart.Phi());
- }
- else {
- AliWarning(Form("%s - Embedded particle not found!", GetName()));
+ 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]);
}
- }
-
- return kTRUE;
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskSAJF::GetLeadingJets(Int_t &maxJetIndex, Int_t &max2JetIndex)
-{
- // Get the leading jets.
- if (!fJets)
- return;
-
- const Int_t njets = fJets->GetEntriesFast();
-
- Float_t maxJetPt = -999;
- Float_t max2JetPt = -999;
- 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;
-
- Float_t corrPt = jet->Pt() - fRhoVal * jet->Area();
+ 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]);
- if (maxJetIndex == -1 || maxJetPt < corrPt) {
- max2JetPt = maxJetPt;
- max2JetIndex = maxJetIndex;
- maxJetPt = corrPt;
- maxJetIndex = ij;
- }
- else if (max2JetIndex == -1 || max2JetPt < corrPt) {
- max2JetPt = corrPt;
- max2JetIndex = ij;
+ 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]);
}
}
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskSAJF::DoClusterLoop()
-{
- // Do cluster loop.
- if (!fCaloClusters)
- return;
-
- 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))
- continue;
+ PostData(1, fOutput); // Post data for ALL output slots >0 here, to get at least an empty histogram
- fHistClustersPt[fCentBin]->Fill(cluster->E());
- }
}
//________________________________________________________________________
-void AliAnalysisTaskSAJF::DoTrackLoop()
+Bool_t AliAnalysisTaskSAJF::FillHistograms()
{
- // Do track loop.
-
- if (!fTracks)
- return;
-
- Int_t ntracks = fTracks->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;
- }
+ // Fill histograms.
- AliVTrack* vtrack = dynamic_cast<AliVTrack*>(track);
-
- if (vtrack && !AcceptTrack(vtrack))
- continue;
-
- fHistTracksPt[fCentBin]->Fill(track->Pt());
+ if (!fJets) {
+ AliError(Form("%s - Jet array not provided, returning...", GetName()));
+ return kFALSE;
}
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskSAJF::DoJetLoop()
-{
- // Do the jet loop.
-
- if (!fJets)
- return;
- const Int_t njets = fJets->GetEntriesFast();
-
- TH1F constituents("constituents", "constituents",
- fHistConstituents[0]->GetNbinsX(), fHistConstituents[0]->GetXaxis()->GetXmin(), fHistConstituents[0]->GetXaxis()->GetXmax());
-
- for (Int_t ij = 0; ij < njets; ij++) {
+ for (Int_t ij = 0; ij < fJets->GetEntriesFast(); ij++) {
AliEmcalJet* jet = static_cast<AliEmcalJet*>(fJets->At(ij));
if (!jet) {
AliError(Form("Could not receive jet %d", ij));
continue;
- }
+ }
if (!AcceptJet(jet))
continue;
+ Float_t ptLeading = GetLeadingHadronPt(jet);
Float_t corrPt = jet->Pt() - fRhoVal * jet->Area();
- fHistJetsPt[fCentBin]->Fill(jet->Pt());
- fHistJetsPtArea[fCentBin]->Fill(jet->Pt(), jet->Area());
- fHistJetsCorrPt[fCentBin]->Fill(corrPt);
- fHistJetsCorrPtArea[fCentBin]->Fill(corrPt, jet->Area());
-
- fHistJetPhiEta[fCentBin]->Fill(jet->Eta(), jet->Phi());
-
- fHistMaxPartPtvsJetPt[fCentBin]->Fill(jet->Pt(), jet->MaxPartPt());
- fHistMaxPartPtvsJetCorrPt[fCentBin]->Fill(corrPt, jet->MaxPartPt());
+ // Fill THnSparse
+ Double_t ep = jet->Phi() - fEPV0;
+ while (ep < 0) ep += TMath::Pi();
+ while (ep >= TMath::Pi()) ep -= TMath::Pi();
- if (fAnaType != kEMCALOnly) {
- fHistMaxTrackPtvsJetPt[fCentBin]->Fill(jet->Pt(), jet->MaxTrackPt());
- fHistMaxTrackPtvsJetCorrPt[fCentBin]->Fill(corrPt, jet->MaxTrackPt());
- }
-
- if (fAnaType == kEMCAL || fAnaType == kEMCALOnly) {
- fHistMaxClusPtvsJetPt[fCentBin]->Fill(jet->Pt(), jet->MaxClusterPt());
- fHistMaxClusPtvsJetCorrPt[fCentBin]->Fill(corrPt, jet->MaxClusterPt());
- fHistJetsNEFvsPt[fCentBin]->Fill(jet->NEF(), jet->Pt());
- }
-
- Float_t scalarpt = 0;
+ 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) {
- fHistJetsZvsPt[fCentBin]->Fill(track->Pt() / jet->Pt(), jet->Pt());
- constituents.Fill(track->Pt());
- scalarpt += track->Pt();
+ 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 (cluster) {
TLorentzVector nPart;
cluster->GetMomentum(nPart, fVertex);
- fHistJetsZvsPt[fCentBin]->Fill(nPart.Et() / jet->Pt(), jet->Pt());
- scalarpt += nPart.Pt();
- constituents.Fill(nPart.Pt());
+
+ 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);
}
}
}
-
- fHistDeltaVectorPt->Fill(scalarpt - jet->Pt());
-
- for (Int_t i = 1; i <= constituents.GetNbinsX(); i++) {
- fHistConstituents[fCentBin]->Fill(constituents.GetBinCenter(i), constituents.GetBinContent(i));
- }
-
- constituents.Reset();
} //jet loop
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskSAJF::DoEmbTrackLoop()
-{
- // Do track loop.
-
- if (!fEmbTracks)
- return;
-
- fEmbeddedTrackId = -1;
-
- Int_t ntracks = fEmbTracks->GetEntriesFast();
-
- for (Int_t i = 0; i < ntracks; i++) {
-
- AliVParticle* track = static_cast<AliVParticle*>(fEmbTracks->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;
-
- if (track->GetLabel() == 100) {
- fEmbeddedTrackId = i;
- continue;
- }
- }
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskSAJF::DoEmbClusterLoop()
-{
- // Do cluster loop.
-
- if (!fEmbCaloClusters)
- return;
-
- fEmbeddedClusterId = -1;
-
- Int_t nclusters = fEmbCaloClusters->GetEntriesFast();
-
- for (Int_t iClusters = 0; iClusters < nclusters; iClusters++) {
- AliVCluster* cluster = static_cast<AliVCluster*>(fEmbCaloClusters->At(iClusters));
- if (!cluster) {
- AliError(Form("Could not receive cluster %d", iClusters));
- continue;
- }
-
- if (!AcceptCluster(cluster, kTRUE))
- continue;
-
- if (cluster->Chi2() == 100) {
- fEmbeddedClusterId = iClusters;
- continue;
- }
- }
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskSAJF::DoEmbJetLoop(AliEmcalJet* &embJet, TObject* &embPart)
-{
- // Do the embedded jet loop.
-
- if (!fEmbJets)
- return;
-
- TLorentzVector maxClusVect;
-
- Int_t nembjets = fEmbJets->GetEntriesFast();
-
- for (Int_t ij = 0; ij < nembjets; ij++) {
-
- AliEmcalJet* jet = static_cast<AliEmcalJet*>(fEmbJets->At(ij));
-
- if (!jet) {
- AliError(Form("Could not receive jet %d", ij));
- continue;
- }
-
- if (!AcceptJet(jet))
- continue;
-
- if (!jet->IsMC())
- continue;
- AliVParticle *maxTrack = 0;
-
- for (Int_t it = 0; it < jet->GetNumberOfTracks(); it++) {
- AliVParticle *track = jet->TrackAt(it, fEmbTracks);
-
- if (!track)
- continue;
-
- if (track->GetLabel() != 100)
- continue;
-
- if (!maxTrack || track->Pt() > maxTrack->Pt())
- maxTrack = track;
- }
-
- AliVCluster *maxClus = 0;
-
- for (Int_t ic = 0; ic < jet->GetNumberOfClusters(); ic++) {
- AliVCluster *cluster = jet->ClusterAt(ic, fEmbCaloClusters);
-
- if (!cluster)
- continue;
-
- if (cluster->Chi2() != 100)
- continue;
-
- TLorentzVector nPart;
- cluster->GetMomentum(nPart, fVertex);
-
- if (!maxClus || nPart.Et() > maxClusVect.Et()) {
- new (&maxClusVect) TLorentzVector(nPart);
- maxClus = cluster;
- }
- }
-
- if ((maxClus && maxTrack && maxClusVect.Et() > maxTrack->Pt()) || (maxClus && !maxTrack)) {
- embPart = maxClus;
- embJet = jet;
- }
- else if (maxTrack) {
- embPart = maxTrack;
- embJet = jet;
- }
-
- return; // MC jet found, exit
- }
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskSAJF::GetRigidCone(Float_t &pt, Float_t &ptrigid, Float_t &eta, Float_t &phi,
- AliEmcalJet *jet, TClonesArray* tracks, TClonesArray* clusters) const
-{
- // Get rigid cone.
-
- if (!tracks)
- tracks = fTracks;
-
- if (!clusters)
- clusters = fCaloClusters;
-
- eta = 0;
- phi = 0;
- pt = 0;
- ptrigid = 0;
-
- if (!tracks && !clusters)
- return;
-
- Float_t LJeta = 999;
- Float_t LJphi = 999;
-
- if (jet) {
- LJeta = jet->Eta();
- LJphi = jet->Phi();
- }
-
- Float_t maxEta = fMaxEta;
- Float_t minEta = fMinEta;
- Float_t maxPhi = fMaxPhi;
- Float_t minPhi = fMinPhi;
-
- if (maxPhi > TMath::Pi() * 2) maxPhi = TMath::Pi() * 2;
- if (minPhi < 0) minPhi = 0;
-
- Float_t dLJ = 0;
- Int_t repeats = 0;
- do {
- eta = gRandom->Rndm() * (maxEta - minEta) + minEta;
- phi = gRandom->Rndm() * (maxPhi - minPhi) + minPhi;
- dLJ = TMath::Sqrt((LJeta - eta) * (LJeta - eta) + (LJphi - phi) * (LJphi - phi));
- repeats++;
- } while (dLJ < fMinRC2LJ && repeats < 999);
-
- if (repeats == 999) {
- AliWarning(Form("%s: Could not get random cone!", GetName()));
- return;
- }
-
- TVector3 rigidAxis;
- rigidAxis.SetPtEtaPhi(1, eta, phi);
- rigidAxis = rigidAxis.Unit();
-
- if ((fAnaType == kEMCAL || fAnaType == kEMCALOnly) && clusters) {
- Int_t nclusters = clusters->GetEntriesFast();
- for (Int_t iClusters = 0; iClusters < nclusters; iClusters++) {
- AliVCluster* cluster = static_cast<AliVCluster*>(clusters->At(iClusters));
- if (!cluster) {
- AliError(Form("Could not receive cluster %d", iClusters));
- continue;
- }
-
- if (!AcceptCluster(cluster, fMCAna))
- continue;
-
- TLorentzVector nPart;
- cluster->GetMomentum(nPart, const_cast<Double_t*>(fVertex));
-
- Float_t d = TMath::Sqrt((nPart.Eta() - eta) * (nPart.Eta() - eta) + (nPart.Phi() - phi) * (nPart.Phi() - phi));
-
- if (d <= fJetRadius) {
- TVector3 vect = nPart.Vect();
- vect *= vect * rigidAxis / vect.Mag();
- ptrigid += vect.Pt();
-
- pt += nPart.Pt();
- }
- }
- }
-
- if (tracks) {
- Int_t ntracks = tracks->GetEntriesFast();
- for(Int_t iTracks = 0; iTracks < ntracks; iTracks++) {
- AliVTrack* track = static_cast<AliVTrack*>(tracks->At(iTracks));
- if(!track) {
- AliError(Form("Could not retrieve track %d",iTracks));
- continue;
- }
-
- if (!AcceptTrack(track, fMCAna))
- continue;
-
- Float_t tracketa = track->Eta();
- Float_t trackphi = track->Phi();
-
- if (TMath::Abs(trackphi - phi) > TMath::Abs(trackphi - phi + 2 * TMath::Pi()))
- trackphi += 2 * TMath::Pi();
- if (TMath::Abs(trackphi - phi) > TMath::Abs(trackphi - phi - 2 * TMath::Pi()))
- trackphi -= 2 * TMath::Pi();
-
- Float_t d = TMath::Sqrt((tracketa - eta) * (tracketa - eta) + (trackphi - phi) * (trackphi - phi));
- if (d <= fJetRadius){
- TVector3 vect(track->Px(), track->Py(), track->Pz());
- vect *= vect * rigidAxis / vect.Mag();
- ptrigid += vect.Pt();
-
- pt += track->Pt();
- }
- }
- }
+ return kTRUE;
}
//________________________________________________________________________
-void AliAnalysisTaskSAJF::ExecOnce()
+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)
{
- // Initialize the analysis.
-
- if (!fRhoName.IsNull() && !fRho) {
- fRho = dynamic_cast<AliRhoParameter*>(InputEvent()->FindListObject(fRhoName));
- if (!fRho) {
- AliError(Form("%s: Could not retrieve rho %s!", GetName(), fRhoName.Data()));
- return;
- }
+ 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);
}
-
- if (!fEmbJetsName.IsNull() && !fEmbJets) {
- fEmbJets = dynamic_cast<TClonesArray*>(InputEvent()->FindListObject(fEmbJetsName));
- if (!fEmbJets) {
- AliError(Form("%s: Could not retrieve emb jets %s!", GetName(), fEmbJetsName.Data()));
- return;
- }
- else if (!fEmbJets->GetClass()->GetBaseClass("AliEmcalJet")) {
- AliError(Form("%s: Collection %s does not contain AliEmcalJet objects!", GetName(), fEmbJetsName.Data()));
- fEmbJets = 0;
- return;
- }
- }
-
- if (!fEmbCaloName.IsNull() && (fAnaType == kEMCAL || fAnaType == kEMCALOnly) && !fEmbCaloClusters) {
- fEmbCaloClusters = dynamic_cast<TClonesArray*>(InputEvent()->FindListObject(fEmbCaloName));
- if (!fEmbCaloClusters) {
- AliError(Form("%s: Could not retrieve embedded clusters %s!", GetName(), fEmbCaloName.Data()));
- return;
- }
- else if (!fEmbCaloClusters->GetClass()->GetBaseClass("AliVCluster") && !fEmbCaloClusters->GetClass()->GetBaseClass("AliEmcalParticle")) {
- AliError(Form("%s: Collection %s does not contain AliVCluster nor AliEmcalParticle objects!", GetName(), fEmbCaloName.Data()));
- fEmbCaloClusters = 0;
- return;
- }
- }
-
- if (!fEmbTracksName.IsNull() && fAnaType != kEMCALOnly && !fEmbTracks) {
- fEmbTracks = dynamic_cast<TClonesArray*>(InputEvent()->FindListObject(fEmbTracksName));
- if (!fEmbTracks) {
- AliError(Form("%s: Could not retrieve embedded tracks %s!", GetName(), fEmbTracksName.Data()));
- return;
- }
- else if (!fEmbTracks->GetClass()->GetBaseClass("AliVParticle") && !fEmbTracks->GetClass()->GetBaseClass("AliEmcalParticle")) {
- AliError(Form("%s: Collection %s does not contain AliVParticle nor AliEmcalParticle objects!", GetName(), fEmbTracksName.Data()));
- fEmbTracks = 0;
- return;
- }
- }
-
- if (!fRandCaloName.IsNull() && (fAnaType == kEMCAL || fAnaType == kEMCALOnly) && !fRandCaloClusters) {
- fRandCaloClusters = dynamic_cast<TClonesArray*>(InputEvent()->FindListObject(fRandCaloName));
- if (!fRandCaloClusters) {
- AliError(Form("%s: Could not retrieve randomized clusters %s!", GetName(), fRandCaloName.Data()));
- return;
- }
- else if (!fRandCaloClusters->GetClass()->GetBaseClass("AliVCluster") && !fRandCaloClusters->GetClass()->GetBaseClass("AliEmcalParticle")) {
- AliError(Form("%s: Collection %s does not contain AliVCluster nor AliEmcalParticle objects!", GetName(), fRandCaloName.Data()));
- fRandCaloClusters = 0;
- return;
- }
- }
-
- if (!fRandTracksName.IsNull() && fAnaType != kEMCALOnly && !fRandTracks) {
- fRandTracks = dynamic_cast<TClonesArray*>(InputEvent()->FindListObject(fRandTracksName));
- if (!fRandTracks) {
- AliError(Form("%s: Could not retrieve randomized tracks %s!", GetName(), fRandTracksName.Data()));
- return;
- }
- else if (!fRandTracks->GetClass()->GetBaseClass("AliVParticle") && !fRandTracks->GetClass()->GetBaseClass("AliEmcalParticle")) {
- AliError(Form("%s: Collection %s does not contain AliVParticle nor AliEmcalParticle objects!", GetName(), fRandTracksName.Data()));
- fRandTracks = 0;
- return;
- }
- }
-
- AliAnalysisTaskEmcalJet::ExecOnce();
-
- if (fRCperEvent < 0) {
- Double_t area = (fMaxEta - fMinEta) * (fMaxPhi - fMinPhi);
- Double_t jetArea = TMath::Pi() * fJetRadius * fJetRadius;
- fRCperEvent = TMath::FloorNint(area / jetArea - 0.5);
- if (fRCperEvent == 0)
- fRCperEvent = 1;
- }
-
- if (fMinRC2LJ < 0)
- fMinRC2LJ = fJetRadius * 1.5;
-
- const Float_t maxDist = TMath::Max(fMaxPhi - fMinPhi, fMaxEta - fMinEta) / 2;
- if (fMinRC2LJ > maxDist) {
- AliWarning(Form("The parameter fMinRC2LJ = %f is too large for the considered acceptance. "
- "Will use fMinRC2LJ = %f", fMinRC2LJ, maxDist));
- fMinRC2LJ = maxDist;
+ 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")
+ 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);
}
}
-
-//________________________________________________________________________
-void AliAnalysisTaskSAJF::Terminate(Option_t *)
-{
- // Called once at the end of the analysis.
-}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff