correcting conenction of output container
[u/mrichter/AliRoot.git] / PWG4 / macros / AddTaskJetChem.C
1 AliAnalysisTask *AddTaskJetChem(){
2
3
4
5   //cout<<" OB : add JetTasks inlcude path ! "<<endl;
6   //gSystem->AddIncludePath("-I$ALICE_ROOT/PWG4/JetTasks");
7
8   //get the current analysis manager
9   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
10   if (!mgr) {
11     Error("AddTask_obusch_jets", "No analysis manager found.");
12     return 0;
13   }
14
15   // physics event selection task - not for AOD analysis
16
17   //gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
18   //AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection(); 
19   //AliPhysicsSelection* physSel = physSelTask->GetPhysicsSelection();
20   //physSel->AddBackgroundIdentification(new AliBackgroundSelection()); 
21   //physSel->SetAnalyzeMC();
22
23   /*
24   //============= Set Task Name ===================
25   TString taskName=("AliAnalysisTaskJetChem.cxx+");
26   //===============================================
27   //            Load the task
28   gROOT->LoadMacro(taskName.Data());
29   if (gProof){
30     TString taskSO=gSystem->pwd();
31     taskSO+="/";
32     taskSO+=taskName(0,taskName.First('.'))+"_cxx.so";
33     gProof->Exec(Form("gSystem->Load(\"%s\")",taskSO.Data()),kTRUE);
34   }
35   */
36   //========= Add task to the ANALYSIS manager =====
37   AliAnalysisTaskJetChem *task = new AliAnalysisTaskJetChem;
38   
39   // configure task 
40
41
42   task->SetUseLOConeJets();
43   task->SetUseLOConeMCJets();
44   
45   //task->SetUsePythiaJets();
46   task->SetConeRadius(0.4);
47   task->SetTrackPtCutJF(0.150); //
48   task->SetFilterBitJF(0x01);   // official PWG4 high pt filter bit 0x10, but not all AliAnalysisTaskESDFilter had configured this ESDTrackCut
49   task->SetRequireITSRefitJF(); // 0x01 + ITS refit = 0x10 
50   task->SetRejectK0TracksJF();  // uncomment for K0 analysis running jet finder in task - modifies jet spectrum
51
52   task->SetJetPtCut(2.0);
53   //task->SetJetPtCut(0.150); // lower pt cut: for plot of diffractive contribution to jet spectrum (goes up to 2 GeV ...) 
54   task->SetJetEtaCut(0.5);
55
56   task->SetFilterBit((UInt_t) 0X01); // std AOD track cuts
57   task->SetTrackPtCut(0.150); 
58   task->SetTrackEtaCut(0.9);
59
60   task->SetUseOnFlyV0s(); 
61   task->SetCutnSigdEdx(2); 
62  
63   //task->ReadDeltaAOD(); // uncomment for DeltaAODs
64   //task->SelectDeltaAODBranch("bla"); 
65   //task->SelectAODBranch("jetsAOD_FASTKT04");
66   //task->SelectAODBranch("jetsAOD_UA107");
67   //task->SelectAODBranch("jets");
68
69   //task->SelectCollisionCandidates(); // either here or in userExec of task - but not for AODs ...
70
71   //  AliLog::SetGlobalLogLevel(AliLog::kInfo);    // kInfo // kDebug // kFatal
72   //  task->SetDebugLevel(10);
73   //  mgr->SetDebugLevel(10);
74
75   mgr->AddTask(task);
76
77
78   //================================================
79   //              data containers
80   //================================================
81   //            find input container
82   //below the trunk version
83   AliAnalysisDataContainer *cinput  = mgr->GetCommonInputContainer();
84   AliAnalysisDataContainer *coutput = mgr->CreateContainer("PWG4_JetChem", 
85                                                            TList::Class(), AliAnalysisManager::kOutputContainer,
86                                                            Form("%s:PWG4_JetChem",AliAnalysisManager::GetCommonFileName()));
87
88   mgr->ConnectInput(task,0,cinput );
89   mgr->ConnectOutput(task,1,coutput);
90  
91   return task;
92 }