#include "AliAnalysisTaskEmcal.h"
-#include <TChain.h>
#include <TClonesArray.h>
#include <TList.h>
#include <TObject.h>
//
// Author: S.Aiola
-#include <TObject.h>
-#include <TChain.h>
#include <TClonesArray.h>
#include <TH1F.h>
#include <TH2F.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"
// Embedding
// _________________________________
- if (!fJets)
+ if (!fEmbJets)
return kTRUE;
AliEmcalJet *embJet = 0;
#include <TH2F.h>
#include <TClonesArray.h>
-#include "AliAnalysisManager.h"
#include "AliLog.h"
#include "AliRhoParameter.h"
#include "AliEmcalJet.h"
fCompareRhoScaledName(),
fRhoFunction(0),
fScaleFunction(0),
- fInEventSigmaRho(71.66),
+ fInEventSigmaRho(35.83),
fRhoScaled(0),
fCompareRho(0),
fCompareRhoScaled(0),
fHistRhoScaledvsNcluster(0)
{
// Constructor.
+
for (Int_t i = 0; i < 4; i++) {
- fHistNjUEvsNj[i] = 0;
+ fHistJetNconstVsPt[i] = 0;
+ }
+ for (Int_t i = 0; i < 12; i++) {
fHistNjUEoverNjVsNj[i] = 0;
}
}
fCompareRhoScaledName(),
fRhoFunction(0),
fScaleFunction(0),
- fInEventSigmaRho(71.66),
+ fInEventSigmaRho(35.83),
fRhoScaled(0),
fCompareRho(0),
fCompareRhoScaled(0),
// Constructor.
for (Int_t i = 0; i < 4; i++) {
- fHistNjUEvsNj[i] = 0;
+ fHistJetNconstVsPt[i] = 0;
+ }
+ for (Int_t i = 0; i < 12; i++) {
fHistNjUEoverNjVsNj[i] = 0;
}
fOutput->Add(fHistNjetvsNtrack);
}
+
+ TString name;
for (Int_t i = 0; i < 4; i++) {
- TString name1(Form("NjUEvsNj_%d",i));
- fHistNjUEvsNj[i] = new TH2F(name1, name1, 150, -0.5, 149.5, 150, -0.5, 149.5);
- fHistNjUEvsNj[i]->GetXaxis()->SetTitle("N_{jet}");
- fHistNjUEvsNj[i]->GetYaxis()->SetTitle("N_{jet_{UE}}");
- fOutput->Add(fHistNjUEvsNj[i]);
-
- TString name2(Form("NjUEoverNjVsNj_%d",i));
- fHistNjUEoverNjVsNj[i] = new TH2F(name2, name2, 150, -0.5, 149.5, 120, -0.01, 1.19);
- fHistNjUEoverNjVsNj[i]->GetXaxis()->SetTitle("N_{jet_{UE}} / N_{jet}");
- fHistNjUEoverNjVsNj[i]->GetYaxis()->SetTitle("N_{jet}");
- fOutput->Add(fHistNjUEoverNjVsNj[i]);
+ name = Form("fHistJetNconstVsPt_%d",i);
+ fHistJetNconstVsPt[i] = new TH2F(name, name, 150, -0.5, 149.5, fNbins, fMinBinPt, fMaxBinPt);
+ fHistJetNconstVsPt[i]->GetXaxis()->SetTitle("# of constituents");
+ fHistJetNconstVsPt[i]->GetYaxis()->SetTitle("p_{T,jet} (GeV/c)");
+ fOutput->Add(fHistJetNconstVsPt[i]);
+
+ for (Int_t j = 0; j < 3; j++) {
+ name = Form("NjUEoverNjVsNj_%d_%d",i,j+1);
+ fHistNjUEoverNjVsNj[i*3+j] = new TH2F(name, name, 150, -0.5, 149.5, 120, 0.01, 1.21);
+ fHistNjUEoverNjVsNj[i*3+j]->GetXaxis()->SetTitle("N_{jet}");
+ fHistNjUEoverNjVsNj[i*3+j]->GetYaxis()->SetTitle("N_{jet_{UE}} / N_{jet}");
+ fOutput->Add(fHistNjUEoverNjVsNj[i*3+j]);
+ }
}
}
Nclusters = fCaloClusters->GetEntriesFast();
if (fJets) {
- Int_t Njets = fJets->GetEntries();
- Int_t NjetAcc = 0;
- Int_t NjetUE = 0;
- Double_t rhoPlusSigma = fRho->GetVal() + fInEventSigmaRho;
+ Int_t Njets = fJets->GetEntries();
+ Int_t NjetAcc = 0;
+ Int_t NjetUE1Sigma = 0;
+ Int_t NjetUE2Sigma = 0;
+ Int_t NjetUE3Sigma = 0;
+ Double_t rhoPlus1Sigma = fRho->GetVal() + fInEventSigmaRho;
+ Double_t rhoPlus2Sigma = fRho->GetVal() + 2*fInEventSigmaRho;
+ Double_t rhoPlus3Sigma = fRho->GetVal() + 3*fInEventSigmaRho;
for (Int_t i = 0; i < Njets; ++i) {
fHistJetAreavsNtrack->Fill(Ntracks, jet->Area());
}
- if (jet->Pt() < rhoPlusSigma * jet->Area())
- NjetUE++;
+ fHistJetNconstVsPt[fCentBin]->Fill(jet->GetNumberOfConstituents(), jet->Pt());
+
+ if (jet->Pt() < rhoPlus1Sigma * jet->Area())
+ NjetUE1Sigma++;
+
+ if (jet->Pt() < rhoPlus2Sigma * jet->Area())
+ NjetUE2Sigma++;
+
+ if (jet->Pt() < rhoPlus3Sigma * jet->Area())
+ NjetUE3Sigma++;
NjetAcc++;
}
-
- fHistNjUEvsNj[fCentBin]->Fill(NjetAcc,NjetUE);
- if (NjetAcc>0)
- fHistNjUEoverNjVsNj[fCentBin]->Fill(NjetAcc,1.*NjetUE/NjetAcc);
+ if (NjetAcc>0) {
+ fHistNjUEoverNjVsNj[fCentBin*3 ]->Fill(NjetAcc,1.*NjetUE1Sigma/NjetAcc);
+ fHistNjUEoverNjVsNj[fCentBin*3+1]->Fill(NjetAcc,1.*NjetUE2Sigma/NjetAcc);
+ fHistNjUEoverNjVsNj[fCentBin*3+2]->Fill(NjetAcc,1.*NjetUE3Sigma/NjetAcc);
+ }
fHistNjetvsCent->Fill(fCent, NjetAcc);
if (fTracks)
TH2F *fHistJetPtvsNtrack; //!jet pt vs. no. of tracks
TH2F *fHistJetAreavsNtrack; //!jet area vs. no. of tracks
TH2F *fHistNjetvsNtrack; //!no. of jets vs. no. of tracks
- TH2F *fHistNjUEvsNj[4]; //!no. of jets below rho*A+sigma_rho vs. no. of jets
- TH2F *fHistNjUEoverNjVsNj[4]; //!ratio no. of jets below rho*A+sigma_rho over. no. of jets vs. no. of jets
+ TH2F *fHistNjUEoverNjVsNj[12]; //!ratio no. of jets below rho*A+sigma_rho over. no. of jets vs. no. of jets
+ TH2F *fHistJetNconstVsPt[4]; //!jet no. of constituents vs. pt
TH2F *fHistRhovsCent; //!rho vs. centrality
TH2F *fHistRhoScaledvsCent; //!rhoscaled vs. centrality
AliAnalysisTaskRhoBase(const AliAnalysisTaskRhoBase&); // not implemented
AliAnalysisTaskRhoBase& operator=(const AliAnalysisTaskRhoBase&); // not implemented
- ClassDef(AliAnalysisTaskRhoBase, 5); // Rho base task
+ ClassDef(AliAnalysisTaskRhoBase, 6); // Rho base task
};
#endif
AliVCluster *ClusterAt(Int_t idx, TClonesArray *ca) const { if (!ca) return 0; return dynamic_cast<AliVCluster*>(ca->At(ClusterAt(idx))); }
UShort_t GetNumberOfClusters() const { return fClusterIDs.GetSize(); }
UShort_t GetNumberOfTracks() const { return fTrackIDs.GetSize(); }
+ UShort_t GetNumberOfConstituents() const { return GetNumberOfClusters()+GetNumberOfTracks(); }
Double_t FracEmcalArea() const { return fAreaEmc/fArea; }
Bool_t IsInsideEmcal() const { return (fAreaEmc/fArea>0.999); }
Bool_t IsInEmcal() const { return (fAreaEmc>0); }
Double_t fPtSub; //! background subtracted pt (not stored set from outside)
Double_t fPtVectSub; //! background vector subtracted pt (not stored set from outside)
- ClassDef(AliEmcalJet,8) // Emcal jet class in cylindrical coordinates
+ ClassDef(AliEmcalJet,9) // Emcal jet class in cylindrical coordinates
};
#endif
//
// Author: S.Aiola
-#include <TObject.h>
-#include <TChain.h>
#include <TClonesArray.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TH3F.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),
fLeadingHadronType(0),
- fHistRhoVSleadJetPt(0)
+ fHistRhoVSleadJetPt(0),
+ fNjetsVsCent(0)
{
// Default constructor.
fHistConstituents[i] = 0;
fHistTracksJetPt[i] = 0;
fHistClustersJetPt[i] = 0;
+ fHistJetNconstVsPt[i] = 0;
}
SetMakeGeneralHistograms(kTRUE);
AliAnalysisTaskSAJF::AliAnalysisTaskSAJF(const char *name) :
AliAnalysisTaskEmcalJet(name, kTRUE),
fLeadingHadronType(0),
- fHistRhoVSleadJetPt(0)
+ fHistRhoVSleadJetPt(0),
+ fNjetsVsCent(0)
{
// Standard constructor.
fHistConstituents[i] = 0;
fHistTracksJetPt[i] = 0;
fHistClustersJetPt[i] = 0;
+ fHistJetNconstVsPt[i] = 0;
}
SetMakeGeneralHistograms(kTRUE);
AliAnalysisTaskEmcalJet::UserCreateOutputObjects();
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]");
+ fHistRhoVSleadJetPt->GetXaxis()->SetTitle("#rho * area (GeV/c)");
+ fHistRhoVSleadJetPt->GetYaxis()->SetTitle("Leading jet p_{T} (GeV/c)");
fOutput->Add(fHistRhoVSleadJetPt);
+ fNjetsVsCent = new TH2F("fNjetsVsCent","fNjetsVsCent", 100, 0, 100, 150, -0.5, 149.5);
+ fNjetsVsCent->GetXaxis()->SetTitle("Centrality (%)");
+ fNjetsVsCent->GetYaxis()->SetTitle("# of jets");
+ fOutput->Add(fNjetsVsCent);
+
const Int_t nbinsZ = 12;
Float_t binsZ[nbinsZ+1] = {0,1,2,3,4,5,6,7,8,9,10,20,1000};
Float_t *binsEta = GenerateFixedBinArray(50,-1, 1);
Float_t *binsPhi = GenerateFixedBinArray(101, 0, TMath::Pi() * 2.02);
Float_t *bins120 = GenerateFixedBinArray(120, 0, 1.2);
+ Float_t *bins150 = GenerateFixedBinArray(150, -0.5, 149.5);
TString histname;
histname = "fHistConstituents_";
histname += i;
- fHistConstituents[i] = new TH3F(histname.Data(), histname.Data(), 100, 1, 101, 100, -0.5, 99.5, fNbins, fMinBinPt, fMaxBinPt);
+ 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");
- fHistConstituents[i]->GetZaxis()->SetTitle("p_{T,jet} (GeV/c)");
+ fHistConstituents[i]->GetZaxis()->SetTitle("counts");
fOutput->Add(fHistConstituents[i]);
histname = "fHistTracksJetPt_";
fHistClustersJetPt[i]->GetZaxis()->SetTitle("counts");
fOutput->Add(fHistClustersJetPt[i]);
}
+
+ histname = "fHistJetNconstVsPt_";
+ histname += i;
+ fHistJetNconstVsPt[i] = new TH3F(histname.Data(), histname.Data(), 150, bins150, fNbins, binsPt, nbinsZ, binsZ);
+ fHistJetNconstVsPt[i]->GetXaxis()->SetTitle("# of constituents");
+ fHistJetNconstVsPt[i]->GetYaxis()->SetTitle("p_{T,jet} (GeV/c)");
+ fHistJetNconstVsPt[i]->GetZaxis()->SetTitle("p_{T,lead} (GeV/c)");
+ fOutput->Add(fHistJetNconstVsPt[i]);
}
PostData(1, fOutput); // Post data for ALL output slots >0 here, to get at least an empty histogram
if (jet2)
fHist2LeadingJetPt[fCentBin]->Fill(jet2->Pt());
- DoJetLoop();
+ Int_t njets = DoJetLoop();
+
+ fNjetsVsCent->Fill(fCent, njets);
return kTRUE;
}
//________________________________________________________________________
-void AliAnalysisTaskSAJF::DoJetLoop()
+Int_t AliAnalysisTaskSAJF::DoJetLoop()
{
// Do the jet loop.
if (!fJets)
- return;
+ return 0;
const Int_t njets = fJets->GetEntriesFast();
+ Int_t nAccJets = 0;
TH1F constituents("constituents", "constituents",
fHistConstituents[0]->GetNbinsX(), fHistConstituents[0]->GetXaxis()->GetXmin(), fHistConstituents[0]->GetXaxis()->GetXmax());
fHistJetPhiEta[fCentBin]->Fill(jet->Eta(), jet->Phi(), ptLeading);
fHistJetsPtArea[fCentBin]->Fill(jet->Pt(), jet->Area(), ptLeading);
fHistJetsCorrPtArea[fCentBin]->Fill(corrPt, jet->Area(), ptLeading);
+ fHistJetNconstVsPt[fCentBin]->Fill(jet->GetNumberOfConstituents(), jet->Pt(), ptLeading);
if (fCaloClusters)
fHistJetsNEFvsPt[fCentBin]->Fill(jet->NEF(), jet->Pt(), ptLeading);
}
for (Int_t i = 1; i <= constituents.GetNbinsX(); i++) {
- fHistConstituents[fCentBin]->Fill(constituents.GetBinCenter(i), constituents.GetBinContent(i), jet->Pt());
+ fHistConstituents[fCentBin]->Fill(constituents.GetBinCenter(i), constituents.GetBinContent(i));
}
constituents.Reset();
+ nAccJets++;
} //jet loop
+
+ return nAccJets;
}
//________________________________________________________________________
protected:
Bool_t FillHistograms() ;
- void DoJetLoop() ;
+ Int_t DoJetLoop() ;
Float_t* GenerateFixedBinArray(Int_t n, Float_t min, Float_t max) const;
Int_t fLeadingHadronType; // 0 = charged, 1 = neutral, 2 = both
TH1F *fHist2LeadingJetPt[4]; //!Second leading jet pt spectrum
TH1F *fHistLeadingJetCorrPt[4]; //!Corrected leading jet pt spectrum
TH2F *fHistRhoVSleadJetPt; //!Area(leadjet) * rho vs. leading jet pt
+ TH2F *fNjetsVsCent; //!No. of jets vs. centrality
// Inclusive jets histograms
TH3F *fHistJetPhiEta[4]; //!Phi-Eta distribution of jets
TH3F *fHistJetsCorrPtArea[4]; //!Jet corr pt vs. area
TH3F *fHistJetsNEFvsPt[4]; //!Jet neutral energy fraction vs. jet pt
TH3F *fHistJetsZvsPt[4]; //!Constituent Pt over Jet Pt ratio vs. jet pt
- TH3F *fHistConstituents[4]; //!x axis = constituents pt; y axis = no. of constituents; z axis = jet pt
+ TH2F *fHistConstituents[4]; //!x axis = constituents pt; y axis = no. of constituents
TH2F *fHistTracksJetPt[4]; //!Track pt vs. jet pt
TH2F *fHistClustersJetPt[4]; //!Cluster pt vs. jet pt
+ TH3F *fHistJetNconstVsPt[4]; //!Jet no. of constituents vs. pt
private:
AliAnalysisTaskSAJF(const AliAnalysisTaskSAJF&); // not implemented
AliAnalysisTaskSAJF &operator=(const AliAnalysisTaskSAJF&); // not implemented
- ClassDef(AliAnalysisTaskSAJF, 11) // jet analysis task
+ ClassDef(AliAnalysisTaskSAJF, 12) // jet analysis task
};
#endif
+// $Id$
AliAnalysisTaskDeltaPt* AddTaskDeltaPt(
const char *ntracks = "Tracks",
const char *nrho = "Rho",
Double_t jetradius = 0.2,
Double_t jetptcut = 1,
- Double_t jetareacut = 0.8,
- Double_t ptcut = 0.15,
+ Double_t jetareacut = 0.557,
+ Double_t trackptcut = 0.15,
+ Double_t clusptcut = 0.30,
UInt_t type = AliAnalysisTaskEmcal::kTPC,
const char *taskname = "AliAnalysisTaskDeltaPt"
)
jetTask->SetRandTracksName(nrandtracks);
jetTask->SetRandClusName(nrandclusters);
jetTask->SetRhoName(nrho);
- jetTask->SetPtCut(ptcut);
+ jetTask->SetClusPtCut(clusptcut);
+ jetTask->SetTrackPtCut(trackptcut);
jetTask->SetJetRadius(jetradius);
jetTask->SetJetPtCut(jetptcut);
jetTask->SetPercAreaCut(jetareacut);
const char *njets = "Jets",
const char *nmcjets = "MCJets",
const char *nmctracks = "MCParticles",
- Double_t jetradius = 0.4,
+ Double_t jetradius = 0.2,
Double_t jetptcut = 1,
- Double_t jetareacut = 0.8,
- Double_t ptcut = 0.15,
- Double_t jetBiasTrack = 10,
- Double_t jetBiasClus = 10,
+ Double_t jetareacut = 0.557,
+ Double_t jetBiasTrack = 5,
+ Double_t jetBiasClus = 5,
Double_t maxDistance = 0.25,
UInt_t type = AliAnalysisTaskEmcal::kTPC,
Int_t minPtBin = 1,
for (Int_t i = minPtBin; i <= maxPtBin; i++) {
AddTaskJetResponseMaker(ntracks, nclusters, njets, nmcjets, nmctracks,
- jetradius, jetptcut, jetareacut, ptcut, jetBiasTrack,
+ jetradius, jetptcut, jetareacut, jetBiasTrack,
jetBiasClus, maxDistance, type, i, taskname, jetTask + i - minPtBin);
jetTask[i - minPtBin].SetDoMatching(domatch);
if (biggerMatrix)
const char *ntracks = "Tracks",
const char *nclusters = "CaloClusters",
const char *njets = "Jets",
- const char *nmcjets = "MCJets",
const char *nmctracks = "MCParticles",
- Double_t jetradius = 0.4,
+ const char *nmcjets = "MCJets",
+ Double_t jetradius = 0.2,
Double_t jetptcut = 1,
- Double_t jetareacut = 0.8,
- Double_t ptcut = 0.15,
- Double_t jetBiasTrack = 10,
- Double_t jetBiasClus = 10,
+ Double_t jetareacut = 0.557,
+ Double_t jetBiasTrack = 5,
+ Double_t jetBiasClus = 5,
Double_t maxDistance = 0.25,
UInt_t type = AliAnalysisTaskEmcal::kTPC,
Int_t ptHardBin = -999,
// Init the task and do settings
//-------------------------------------------------------
- TString name(taskname);
- name += "_";
- name += njets;
- name += "_Track";
- name += jetBiasTrack;
- name += "_Clus";
- name += jetBiasClus;
- name += "_R0";
- name += floor(jetradius*10+0.5);
- name += "_PtCut";
- name += floor(ptcut*1000+0.5);
- name += "_";
- if (type == AliAnalysisTaskEmcal::kTPC)
+ TString name(Form("%s_%s_Track%f_Clus%f_R0%d_",taskname,njets,jetBiasTrack,jetBiasClus,(Int_t)floor(jetradius*100+0.5)));
+ if (type == AliAnalysisTaskEmcal::kTPC)
name += "TPC";
else if (type == AliAnalysisTaskEmcal::kEMCAL)
name += "EMCAL";
- else if (type == AliAnalysisTaskEmcal::kTPCSmall)
- name += "TPCSmall";
- else if (type == AliAnalysisTaskEmcal::kEMCALOnly)
- name += "EMCALOnly";
+ else if (type == AliAnalysisTaskEmcal::kUser)
+ name += "USER";
if (ptHardBin != -999) {
name += "_PtHard";
name += ptHardBin;
jetTask->SetJetsName(njets);
jetTask->SetMCJetsName(nmcjets);
jetTask->SetMCTracksName(nmctracks);
- jetTask->SetPtCut(ptcut);
jetTask->SetJetRadius(jetradius);
jetTask->SetJetPtCut(jetptcut);
jetTask->SetPercAreaCut(jetareacut);
UInt_t type = AliAnalysisTaskEmcal::kTPC,
Double_t jetareacut = 0.01,
Double_t emcareacut = 0,
- Double_t ptcut = 0.15,
TF1 *sfunc = 0,
const UInt_t exclJets = 2,
const Bool_t histo = kFALSE,
rhotask->SetRhoName(nRho);
rhotask->SetJetAreaCut(jetareacut);
rhotask->SetAreaEmcCut(emcareacut);
- rhotask->SetPtCut(ptcut);
rhotask->SetJetPtCut(0);
rhotask->SetJetRadius(jetradius);
rhotask->SetExcludeLeadJets(exclJets);
UInt_t type = AliAnalysisTaskEmcal::kTPC,
Double_t jetareacut = 0.01,
Double_t emcareacut = 0,
- Double_t ptcut = 0.15,
+ Double_t trackptcut = 0.15,
+ Double_t clusptcut = 0.30,
TF1 *sfunc = 0,
const UInt_t exclPart = 2,
const UInt_t rhotype = 1,
rhotask->SetRhoName(nRho);
rhotask->SetJetAreaCut(jetareacut);
rhotask->SetAreaEmcCut(emcareacut);
- rhotask->SetPtCut(ptcut);
+ rhotask->SetClusPtCut(clusptcut);
+ rhotask->SetTrackPtCut(trackptcut);
rhotask->SetJetPtCut(0);
rhotask->SetJetRadius(jetradius);
rhotask->SetExcludeLeadPart(exclPart);
UInt_t type = AliAnalysisTaskEmcal::kTPC,
Double_t jetareacut = 0.01,
Double_t emcareacut = 0,
- Double_t ptcut = 0.15,
TF1 *sfunc = 0,
TF1 *rfunc = 0,
const Bool_t histo = kFALSE,
rhotask->SetRhoName(nRho);
rhotask->SetJetAreaCut(jetareacut);
rhotask->SetAreaEmcCut(emcareacut);
- rhotask->SetPtCut(ptcut);
rhotask->SetJetPtCut(0);
rhotask->SetJetRadius(jetradius);
const char *nrho = "Rho",
Double_t jetradius = 0.2,
Double_t jetptcut = 1,
- Double_t jetareacut = 0.8,
- Double_t ptcut = 0.15,
+ Double_t jetareacut = 0.557,
UInt_t type = AliAnalysisTaskEmcal::kTPC,
Int_t leadhadtype = 0,
const char *taskname = "AliAnalysisTaskSAJF"
jetTask->SetClusName(nclusters);
jetTask->SetJetsName(njets);
jetTask->SetRhoName(nrho);
- jetTask->SetPtCut(ptcut);
jetTask->SetJetRadius(jetradius);
jetTask->SetJetPtCut(jetptcut);
jetTask->SetPercAreaCut(jetareacut);
const char *njets = "Jets",
Double_t jetradius = 0.2,
Double_t jetptcut = 1,
- Double_t jetareacut = 0.8,
- Double_t ptcut = 0.15,
+ Double_t jetareacut = 0.557,
+ Double_t trackptcut = 0.15,
+ Double_t clusptcut = 0.30,
UInt_t type = AliAnalysisTaskEmcal::kTPC,
const char *taskname = "AliAnalysisTaskSAQA"
)
//-------------------------------------------------------
// Init the task and do settings
//-------------------------------------------------------
- TString name(Form("%s_%s_%s_%s_R0%d_",taskname,njets,ntracks,nclusters,(Int_t)floor(jetradius*100+0.5)));
- if (type == AliAnalysisTaskEmcal::kTPC)
- name += "TPC";
- else if (type == AliAnalysisTaskEmcal::kEMCAL)
- name += "EMCAL";
- else if (type == AliAnalysisTaskEmcal::kUser)
- name += "USER";
+
+ TString name;
+ if (strcmp(njets,"")) {
+ name = Form("%s_%s_%s_%s_R0%d_",taskname,njets,ntracks,nclusters,(Int_t)floor(jetradius*100+0.5));
+ if (type == AliAnalysisTaskEmcal::kTPC)
+ name += "TPC";
+ else if (type == AliAnalysisTaskEmcal::kEMCAL)
+ name += "EMCAL";
+ else if (type == AliAnalysisTaskEmcal::kUser)
+ name += "USER";
+ }
+ else {
+ name = Form("%s_%s_%s_",taskname,ntracks,nclusters);
+ }
+
AliAnalysisTaskSAQA* qaTask = new AliAnalysisTaskSAQA(name);
qaTask->SetTracksName(ntracks);
qaTask->SetClusName(nclusters);
qaTask->SetJetRadius(jetradius);
qaTask->SetJetPtCut(jetptcut);
qaTask->SetPercAreaCut(jetareacut);
- qaTask->SetPtCut(ptcut);
+ qaTask->SetTrackPtCut(trackptcut);
+ qaTask->SetClusPtCut(clusptcut);
qaTask->SetAnaType(type);
//-------------------------------------------------------