]>
Commit | Line | Data |
---|---|---|
624bef5b | 1 | // $Id$ |
a487deae | 2 | |
3 | AliAnalysisTaskDeltaPt* AddTaskDeltaPt( | |
4 | const char *ntracks = "Tracks", | |
5 | const char *nclusters = "CaloClusters", | |
6 | const char *njets = "Jets", | |
7 | const char *nembtracks = "TracksEmbedded", | |
8 | const char *nembclusters = "CaloClustersEmbedded", | |
9 | const char *nembjets = "EmbJets", | |
10 | const char *nrandtracks = "TracksRandomized", | |
11 | const char *nrandclusters = "CaloClustersRandomized", | |
12 | const char *nrho = "Rho", | |
13 | Double_t jetradius = 0.2, | |
14 | Double_t jetptcut = 1, | |
624bef5b | 15 | Double_t jetareacut = 0.557, |
16 | Double_t trackptcut = 0.15, | |
17 | Double_t clusptcut = 0.30, | |
a487deae | 18 | UInt_t type = AliAnalysisTaskEmcal::kTPC, |
19 | const char *taskname = "AliAnalysisTaskDeltaPt" | |
20 | ) | |
21 | { | |
22 | // Get the pointer to the existing analysis manager via the static access method. | |
23 | //============================================================================== | |
24 | AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); | |
25 | if (!mgr) | |
26 | { | |
27 | ::Error("AddTaskSAJF", "No analysis manager to connect to."); | |
28 | return NULL; | |
29 | } | |
30 | ||
31 | // Check the analysis type using the event handlers connected to the analysis manager. | |
32 | //============================================================================== | |
33 | if (!mgr->GetInputEventHandler()) | |
34 | { | |
35 | ::Error("AddTaskSAJF", "This task requires an input event handler"); | |
36 | return NULL; | |
37 | } | |
38 | ||
39 | //------------------------------------------------------- | |
40 | // Init the task and do settings | |
41 | //------------------------------------------------------- | |
ca5c29fa | 42 | TString name; |
43 | if (strcmp(ntracks, "") == 0 && strcmp(nclusters, "") == 0) | |
44 | name = Form("%s_%s_R0%d",taskname,nrho,(Int_t)floor(jetradius*100+0.5)); | |
45 | else if (strcmp(ntracks, "") == 0) | |
46 | name = Form("%s_%s_%s_R0%d",taskname,nclusters,nrho,(Int_t)floor(jetradius*100+0.5)); | |
47 | else if (strcmp(nclusters, "") == 0) | |
48 | name = Form("%s_%s_%s_R0%d",taskname,ntracks,nrho,(Int_t)floor(jetradius*100+0.5)); | |
49 | else | |
50 | name = Form("%s_%s_%s_%s_R0%d",taskname,ntracks,nclusters,nrho,(Int_t)floor(jetradius*100+0.5)); | |
51 | ||
a487deae | 52 | if (type == AliAnalysisTaskEmcal::kTPC) |
ca5c29fa | 53 | name += "_TPC"; |
a487deae | 54 | else if (type == AliAnalysisTaskEmcal::kEMCAL) |
ca5c29fa | 55 | name += "_EMCAL"; |
a487deae | 56 | else if (type == AliAnalysisTaskEmcal::kUser) |
ca5c29fa | 57 | name += "_USER"; |
a487deae | 58 | |
59 | AliAnalysisTaskDeltaPt* jetTask = new AliAnalysisTaskDeltaPt(name); | |
60 | jetTask->SetAnaType(type); | |
61 | jetTask->SetTracksName(ntracks); | |
62 | jetTask->SetClusName(nclusters); | |
63 | jetTask->SetJetsName(njets); | |
64 | jetTask->SetEmbTracksName(nembtracks); | |
65 | jetTask->SetEmbClusName(nembclusters); | |
66 | jetTask->SetEmbJetsName(nembjets); | |
67 | jetTask->SetRandTracksName(nrandtracks); | |
68 | jetTask->SetRandClusName(nrandclusters); | |
69 | jetTask->SetRhoName(nrho); | |
624bef5b | 70 | jetTask->SetClusPtCut(clusptcut); |
71 | jetTask->SetTrackPtCut(trackptcut); | |
a487deae | 72 | jetTask->SetJetRadius(jetradius); |
73 | jetTask->SetJetPtCut(jetptcut); | |
74 | jetTask->SetPercAreaCut(jetareacut); | |
75 | ||
76 | //------------------------------------------------------- | |
77 | // Final settings, pass to manager and set the containers | |
78 | //------------------------------------------------------- | |
79 | ||
80 | mgr->AddTask(jetTask); | |
81 | ||
82 | // Create containers for input/output | |
83 | AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer(); | |
84 | TString contname(name); | |
85 | contname += "_histos"; | |
86 | AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(contname.Data(), | |
87 | TList::Class(),AliAnalysisManager::kOutputContainer, | |
88 | Form("%s", AliAnalysisManager::GetCommonFileName())); | |
ca5c29fa | 89 | mgr->ConnectInput(jetTask, 0, cinput1); |
90 | mgr->ConnectOutput(jetTask, 1, coutput1); | |
a487deae | 91 | |
92 | return jetTask; | |
93 | } |