1 // $Id: AddTaskSAQA.C 56113 2012-05-01 21:39:27Z loizides $
3 AliAnalysisTaskSAQA* AddTaskSAQA(
4 const char *ntracks = "Tracks",
5 const char *nclusters = "CaloClusters",
6 const char *njets = "Jets",
7 const char *ntrgclusters = "",
8 Double_t jetradius = 0.4,
10 Double_t jetareacut = 0.4,
11 Double_t ptcut = 0.15,
12 Double_t jetBiasTrack = 5,
13 Double_t jetBiasClus = 5,
14 UInt_t type = AliAnalysisTaskEmcal::kTPC,
15 const char *taskname = "AliAnalysisTaskSAQA"
18 // Get the pointer to the existing analysis manager via the static access method.
19 //==============================================================================
20 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
23 ::Error("AddTaskSAQA", "No analysis manager to connect to.");
27 // Check the analysis type using the event handlers connected to the analysis manager.
28 //==============================================================================
29 if (!mgr->GetInputEventHandler())
31 ::Error("AddTaskSAQA", "This task requires an input event handler");
35 //-------------------------------------------------------
36 // Init the task and do settings
37 //-------------------------------------------------------
38 TString name(taskname);
46 name += floor(jetradius*10+0.5);
48 name += floor(ptcut*1000+0.5);
50 if (type == AliAnalysisTaskEmcal::kTPC)
52 else if (type == AliAnalysisTaskEmcal::kEMCAL)
54 else if (type == AliAnalysisTaskEmcal::kTPCSmall)
56 AliAnalysisTaskSAQA* qaTask = new AliAnalysisTaskSAQA(name);
57 qaTask->SetTracksName(ntracks);
58 qaTask->SetClusName(nclusters);
59 qaTask->SetJetsName(njets);
60 qaTask->SetTrgClusName(ntrgclusters);
61 qaTask->SetJetRadius(jetradius);
62 qaTask->SetJetPtCut(jetptcut);
63 qaTask->SetJetAreaCut(jetareacut);
64 qaTask->SetPtCut(ptcut);
65 qaTask->SetPtBiasJetTrack(jetBiasTrack);
66 qaTask->SetPtBiasJetClus(jetBiasClus);
67 qaTask->SetAnaType(type);
69 //-------------------------------------------------------
70 // Final settings, pass to manager and set the containers
71 //-------------------------------------------------------
75 // Create containers for input/output
76 AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer() ;
78 TString contName(name);
79 contName += "_histos";
80 AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(contName.Data(),
81 TList::Class(),AliAnalysisManager::kOutputContainer,
82 Form("%s", AliAnalysisManager::GetCommonFileName()));
83 mgr->ConnectInput (qaTask, 0, cinput1 );
84 mgr->ConnectOutput (qaTask, 1, coutput1 );