end-of-line normalization
[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / PiKaPr / TPCTOF / AddTaskCombinedHadronSpectra.C
1
2
3 AliAnalysisTask *AddTaskCombinedHadronSpectra(Int_t identifier = 0, Bool_t isMC = kFALSE, Bool_t isTPConly = kFALSE, Bool_t writeOwnFile = kFALSE, Bool_t saveMotherPDG = kFALSE, Bool_t useEvenDcaBinning = kFALSE, Bool_t smallTHnSparse = kFALSE, Double_t nSigmaTPCLow= -3., Double_t nSigmaTPCHigh = 3., Double_t rapidityLow = -0.2, Double_t rapidityHigh = 0.2, Bool_t setTrackCuts = kFALSE, AliESDtrackCuts *ESDtrackCuts = 0){
4
5
6   //get the current analysis manager
7   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
8   if (!mgr) {
9     Error("AddTask_janielsk_CombinedHadron", "No analysis manager found.");
10     return 0;
11   }
12   //============= Set Task Name ===================
13   TString taskName=("AliAnalysisCombinedHadronSpectra.cxx+g");
14   //===============================================
15   //            Load the task
16   gROOT->LoadMacro(taskName.Data());
17
18
19   
20   //========= Add task to the ANALYSIS manager =====
21
22   //normal tracks
23   AliAnalysisCombinedHadronSpectra *task = new AliAnalysisCombinedHadronSpectra("janielskTaskCombinedHadron");
24   task->SelectCollisionCandidates(AliVEvent::kMB);
25
26   //switches
27   if (isMC) task->SetIsMCtrue(isMC);
28   if (isTPConly)task->SetUseTPConlyTracks(isTPConly);
29   if (saveMotherPDG) task->SetSaveMotherPDG(saveMotherPDG);
30   if (useEvenDcaBinning) task->SetEvenDCAbinning(kTRUE);
31   if (smallTHnSparse){
32     task->SetSmallTHnSparse(kTRUE);
33     task->SetTPCnSigmaCuts(nSigmaTPCLow,nSigmaTPCHigh);
34     task->SetRapidityCuts(rapidityLow,rapidityHigh);
35   }
36
37   //initialize task
38   task->Initialize();
39
40   //esd cuts need to be set after initialize or cuts will be replaced by standard cuts in initialize
41   if (setTrackCuts) task->SetESDtrackCuts(ESDtrackCuts);
42
43   //add task to manager
44   mgr->AddTask(task);
45
46
47   
48
49   //================================================
50   //              data containers
51   //================================================
52   //            find input container
53   //below the trunk version
54   AliAnalysisDataContainer *cinput  = mgr->GetCommonInputContainer();
55
56 /*
57   //dumm output container
58   AliAnalysisDataContainer *coutput0 =
59       mgr->CreateContainer(Form("akalweit_tree%i",identifier),
60                            TTree::Class(),
61                            AliAnalysisManager::kExchangeContainer,
62                            Form("akalweit_default%i",identifier));
63
64
65   //define output containers, please use 'username'_'somename'
66   AliAnalysisDataContainer *coutput1 = 
67       mgr->CreateContainer(Form("akalweit_CombinedHadron%i",identifier), TList::Class(),
68                            AliAnalysisManager::kOutputContainer,Form("akalweit_CombinedHadron%i.root",identifier));
69 */
70   if (!writeOwnFile) {
71         AliAnalysisDataContainer *coutput1 =  mgr->CreateContainer(Form("janielsk_CombinedHadron%i",identifier), TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:janielsk_CombinedHadron", AliAnalysisManager::GetCommonFileName())); 
72         }
73   else {
74         AliAnalysisDataContainer *coutput1 =  mgr->CreateContainer(Form("janielsk_CombinedHadron%i",identifier), TList::Class(), AliAnalysisManager::kOutputContainer, Form("janielsk_CombinedHadron.root"));
75         }
76
77
78
79   //connect containers
80
81   //
82   mgr->ConnectInput  (task,  0, cinput );
83   //mgr->ConnectOutput (task,  0, coutput0);
84   mgr->ConnectOutput (task,  1, coutput1);
85
86   return task;
87 }