// Default constructor.
for (Int_t i = 0;i<7;++i){
fHistJetPtvsTrackPt[i] = 0;
+ fHistJetPtScvsTrackPt[i] = 0;
fHistRawJetPtvsTrackPt[i] = 0;
fHistTrackPt[i] = 0;
fHistEP0[i] = 0;
fHistEPAvsC[i] = 0;
fHistJetPtvsdEP[i] = 0;
fHistJetPtvsdEPBias[i] = 0;
+ fHistJetPtvsEP[i] = 0;
+ fHistJetPtvsEPBias[i] = 0;
fHistRhovsEP[i] = 0;
fHistJetPtEtaPhi[i] = 0;
{
for (Int_t i = 0;i<7;++i){
fHistJetPtvsTrackPt[i] = 0;
+ fHistJetPtScvsTrackPt[i] = 0;
fHistRawJetPtvsTrackPt[i] = 0;
fHistTrackPt[i] = 0;
fHistEP0[i] = 0;
fHistEPAvsC[i] = 0;
fHistJetPtvsdEP[i] = 0;
fHistJetPtvsdEPBias[i] = 0;
+ fHistJetPtvsEP[i] = 0;
+ fHistJetPtvsEPBias[i] = 0;
fHistRhovsEP[i] = 0;
fHistJetPtEtaPhi[i] = 0;
}
deltapt->SetCentralityEstimator(CentEst);
TString emcdeltaname(Form("DeltaPt_%s", nRhosChEm));
- AliAnalysisTaskDeltaPt* deltaptEMC = AddTaskDeltaPt(usedTracks,outClusName,nJets,"","","","","",nRhosChEm,radius,1,AreaCut,minTrackPt,minClusterPt,AliAnalysisTaskEmcal::kEMCAL,"DeltaPtTask");
+ AliAnalysisTaskDeltaPt* deltaptEMC = AddTaskDeltaPt(usedTracks,outClusName,nJets,"","","","","",nRhosChEm,radius,1,AreaCut,minTrackPt,minClusterPt,AliAnalysisTaskEmcal::kEMCAL,emcdeltaname);
deltapt->SetCentralityEstimator(CentEst);
gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskScale.C");
- AliAnalysisTaskScale* scaletask = AddTaskScale(usedTracks,outClusName,0.15,0.3);
+
+ Int_t radlabel=(Int_t)floor(radius*100+0.5);
+ Int_t mincluslabel=(Int_t)floor(minClusterPt*1000+0.5);
+ Int_t mintracklabel=(Int_t)floor(minTrackPt*1000+0.5);
+ TString scalename(Form("Scale_R0%d", radlabel));
+
+ AliAnalysisTaskScale* scaletask = AddTaskScale(usedTracks,outClusName,minTrackPt,minClusterPt,scalename);
scaletask->SetCentralityEstimator(CentEst);
scaletask->SetScaleFunction(sfunc);
-
+
}
+
cout << "Running non charged jet finders..." <<endl;
AliEmcalJetTask* jetFinderTaskEm = AddTaskEmcalJet("",outClusName,cANTIKT,radius,cNEUTRALJETS,minTrackPt,minClusterPt);