-\r
-\r
-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){\r
-\r
-\r
- //get the current analysis manager\r
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
- if (!mgr) {\r
- Error("AddTask_janielsk_CombinedHadron", "No analysis manager found.");\r
- return 0;\r
- }\r
- //============= Set Task Name ===================\r
- TString taskName=("AliAnalysisCombinedHadronSpectra.cxx+g");\r
- //===============================================\r
- // Load the task\r
- gROOT->LoadMacro(taskName.Data());\r
-\r
-\r
- \r
- //========= Add task to the ANALYSIS manager =====\r
-\r
- //normal tracks\r
- AliAnalysisCombinedHadronSpectra *task = new AliAnalysisCombinedHadronSpectra("janielskTaskCombinedHadron");\r
- task->SelectCollisionCandidates(AliVEvent::kMB);\r
-\r
- //switches\r
- if (isMC) task->SetIsMCtrue(isMC);\r
- if (isTPConly)task->SetUseTPConlyTracks(isTPConly);\r
- if (saveMotherPDG) task->SetSaveMotherPDG(saveMotherPDG);\r
- if (useEvenDcaBinning) task->SetEvenDCAbinning(kTRUE);\r
- if (smallTHnSparse){\r
- task->SetSmallTHnSparse(kTRUE);\r
- task->SetTPCnSigmaCuts(nSigmaTPCLow,nSigmaTPCHigh);\r
- task->SetRapidityCuts(rapidityLow,rapidityHigh);\r
- }\r
-\r
- //initialize task\r
- task->Initialize();\r
-\r
- //esd cuts need to be set after initialize or cuts will be replaced by standard cuts in initialize\r
- if (setTrackCuts) task->SetESDtrackCuts(ESDtrackCuts);\r
-\r
- //add task to manager\r
- mgr->AddTask(task);\r
-\r
-\r
- \r
-\r
- //================================================\r
- // data containers\r
- //================================================\r
- // find input container\r
- //below the trunk version\r
- AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();\r
-\r
-/*\r
- //dumm output container\r
- AliAnalysisDataContainer *coutput0 =\r
- mgr->CreateContainer(Form("akalweit_tree%i",identifier),\r
- TTree::Class(),\r
- AliAnalysisManager::kExchangeContainer,\r
- Form("akalweit_default%i",identifier));\r
-\r
-\r
- //define output containers, please use 'username'_'somename'\r
- AliAnalysisDataContainer *coutput1 = \r
- mgr->CreateContainer(Form("akalweit_CombinedHadron%i",identifier), TList::Class(),\r
- AliAnalysisManager::kOutputContainer,Form("akalweit_CombinedHadron%i.root",identifier));\r
-*/\r
- if (!writeOwnFile) {\r
- AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(Form("janielsk_CombinedHadron%i",identifier), TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:janielsk_CombinedHadron", AliAnalysisManager::GetCommonFileName())); \r
- }\r
- else {\r
- AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(Form("janielsk_CombinedHadron%i",identifier), TList::Class(), AliAnalysisManager::kOutputContainer, Form("janielsk_CombinedHadron.root"));\r
- }\r
-\r
-\r
-\r
- //connect containers\r
-\r
- //\r
- mgr->ConnectInput (task, 0, cinput );\r
- //mgr->ConnectOutput (task, 0, coutput0);\r
- mgr->ConnectOutput (task, 1, coutput1);\r
-\r
- return task;\r
-}\r
+
+
+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){
+
+
+ //get the current analysis manager
+ AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+ if (!mgr) {
+ Error("AddTask_janielsk_CombinedHadron", "No analysis manager found.");
+ return 0;
+ }
+ //============= Set Task Name ===================
+ TString taskName=("AliAnalysisCombinedHadronSpectra.cxx+g");
+ //===============================================
+ // Load the task
+ gROOT->LoadMacro(taskName.Data());
+
+
+
+ //========= Add task to the ANALYSIS manager =====
+
+ //normal tracks
+ AliAnalysisCombinedHadronSpectra *task = new AliAnalysisCombinedHadronSpectra("janielskTaskCombinedHadron");
+ task->SelectCollisionCandidates(AliVEvent::kMB);
+
+ //switches
+ if (isMC) task->SetIsMCtrue(isMC);
+ if (isTPConly)task->SetUseTPConlyTracks(isTPConly);
+ if (saveMotherPDG) task->SetSaveMotherPDG(saveMotherPDG);
+ if (useEvenDcaBinning) task->SetEvenDCAbinning(kTRUE);
+ if (smallTHnSparse){
+ task->SetSmallTHnSparse(kTRUE);
+ task->SetTPCnSigmaCuts(nSigmaTPCLow,nSigmaTPCHigh);
+ task->SetRapidityCuts(rapidityLow,rapidityHigh);
+ }
+
+ //initialize task
+ task->Initialize();
+
+ //esd cuts need to be set after initialize or cuts will be replaced by standard cuts in initialize
+ if (setTrackCuts) task->SetESDtrackCuts(ESDtrackCuts);
+
+ //add task to manager
+ mgr->AddTask(task);
+
+
+
+
+ //================================================
+ // data containers
+ //================================================
+ // find input container
+ //below the trunk version
+ AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
+
+/*
+ //dumm output container
+ AliAnalysisDataContainer *coutput0 =
+ mgr->CreateContainer(Form("akalweit_tree%i",identifier),
+ TTree::Class(),
+ AliAnalysisManager::kExchangeContainer,
+ Form("akalweit_default%i",identifier));
+
+
+ //define output containers, please use 'username'_'somename'
+ AliAnalysisDataContainer *coutput1 =
+ mgr->CreateContainer(Form("akalweit_CombinedHadron%i",identifier), TList::Class(),
+ AliAnalysisManager::kOutputContainer,Form("akalweit_CombinedHadron%i.root",identifier));
+*/
+ if (!writeOwnFile) {
+ AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(Form("janielsk_CombinedHadron%i",identifier), TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:janielsk_CombinedHadron", AliAnalysisManager::GetCommonFileName()));
+ }
+ else {
+ AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(Form("janielsk_CombinedHadron%i",identifier), TList::Class(), AliAnalysisManager::kOutputContainer, Form("janielsk_CombinedHadron.root"));
+ }
+
+
+
+ //connect containers
+
+ //
+ mgr->ConnectInput (task, 0, cinput );
+ //mgr->ConnectOutput (task, 0, coutput0);
+ mgr->ConnectOutput (task, 1, coutput1);
+
+ return task;
+}