1 AliAnalysisTaskFullpAJets *AddTaskFullpAJets(const char* proj_name, const Double_t jetRadius=0.4, Bool_t IsMC=kFALSE, const char* track_name="PicoTracks", const char* clus_name="caloClusters", const char* corrclus_name="caloClustersCorr", const char* mcpart_name="MCParticles", const char* Centrality_name="V0A", Double_t scaleFactor = 1.45)
3 char *usedTracks = track_name;
4 char *usedClusters = clus_name;
5 char *outClusName = corrclus_name;
6 char *usedMCParticles = mcpart_name;
7 char *centEst = Centrality_name;
8 char *projName = proj_name;
9 const Double_t minTrackPt=0.15;
10 const Double_t minClusterPt=0.30;
11 const Double_t minMCPartPt=0.00;
12 Double_t scaleFac = scaleFactor; // Obtained from previous runs...
14 // Some constants for the jet finders
16 const Int_t cANTIKT=1;
17 const Int_t cFULLJETS=0;
18 const Int_t cCHARGEDJETS=1;
19 const Int_t cNEUTRALJETS=2;
21 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
24 Error("AddTaskJetCommon","No analysis manager found.");
28 // Determine the int of the jet radius for naming purposes
29 Int_t drjet=Int_t(100*jetRadius);
35 TString taskName = Form("AnalysisFullpAJetsR%d",drjet);
36 TString listName = Form("List%sR%d",projName,drjet);
37 TString fileName = Form("%s:FullpAJets", AliAnalysisManager::GetCommonFileName());
39 // Jet finders (RECONSTRUCTED DATA)
40 TString tmpTaskName("");
41 AliEmcalJetTask* jetFinderTask = NULL;
43 gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C");
46 AliAnalysisTaskFullpAJets *task = new AliAnalysisTaskFullpAJets(taskName);
50 task->SetTrackName(usedMCParticles);
51 task->SetClusterName("");
52 task->SetTrackPtCut(minMCPartPt);
53 task->SetClusterPtCut(minMCPartPt);
55 // ########## CHARGED JETS ##########
56 jetFinderTask = AddTaskEmcalJet(usedMCParticles,"",cKT,jetRadius,cCHARGEDJETS,minMCPartPt,minMCPartPt);
57 task->SetkTChargedJetName(jetFinderTask->GetName());
59 jetFinderTask = AddTaskEmcalJet(usedMCParticles,"",cANTIKT,jetRadius,cCHARGEDJETS,minMCPartPt,minMCPartPt);
60 task->SetAkTChargedJetName(jetFinderTask->GetName());
62 // ########## FULL JETS ##########
63 // No Full jets or clusters are used if run over MCParticles!
64 task->SetkTFullJetName("");
65 task->SetAkTFullJetName("");
69 task->SetTrackName(usedTracks);
70 task->SetClusterName(outClusName);
71 task->SetTrackPtCut(minTrackPt);
72 task->SetClusterPtCut(minClusterPt);
74 // ########## CHARGED JETS ##########
75 jetFinderTask = AddTaskEmcalJet(usedTracks,"",cKT,jetRadius,cCHARGEDJETS,minTrackPt,minClusterPt);
76 task->SetkTChargedJetName(jetFinderTask->GetName());
78 jetFinderTask = AddTaskEmcalJet(usedTracks,"",cANTIKT,jetRadius,cCHARGEDJETS,minTrackPt,minClusterPt);
79 task->SetAkTChargedJetName(jetFinderTask->GetName());
81 // ########## FULL JETS ##########
82 jetFinderTask = AddTaskEmcalJet(usedTracks,outClusName,cKT,jetRadius,cFULLJETS,minTrackPt,minClusterPt);
83 task->SetkTFullJetName(jetFinderTask->GetName());
85 jetFinderTask = AddTaskEmcalJet(usedTracks,outClusName,cANTIKT,jetRadius,cFULLJETS,minTrackPt,minClusterPt);
86 task->SetAkTFullJetName(jetFinderTask->GetName());
90 task->SetCentralityTag(centEst);
91 task->SetScaleFactor(scaleFac);
92 task->SelectCollisionCandidates(AliVEvent::kINT7);
97 AliAnalysisDataContainer *coutput = mgr->CreateContainer(listName,TList::Class(),AliAnalysisManager::kOutputContainer,fileName);
98 mgr->ConnectInput(task,0,mgr->GetCommonInputContainer());
99 mgr->ConnectOutput(task,1,coutput);