-AliAnalysisTaskKMeans *AddTaskKMeans()\r
-{\r
-// Creates a dijet task, configures it and adds it to the analysis manager.\r
-\r
- // Get the pointer to the existing analysis manager via the static access method.\r
- //==============================================================================\r
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
- if (!mgr) {\r
- ::Error("AddTaskJets", "No analysis manager to connect to.");\r
- return NULL;\r
- }\r
-\r
- // Check the analysis type using the event handlers connected to the analysis manager.\r
- //==============================================================================\r
- if (!mgr->GetInputEventHandler()) {\r
- ::Error("AddTaskKMeans", "This task requires an input event handler");\r
- return NULL;\r
- }\r
-\r
- // Create the task and configure it.\r
- //===========================================================================\r
-\r
- AliAnalysisTaskKMeans *taskKMeans = new AliAnalysisTaskKMeans("K-Means Analysis");\r
- taskKMeans->SetDebugLevel(0);\r
- AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("AliESDtrackCuts", "Standard");\r
- esdTrackCutsL->SetMinNClustersTPC(50);\r
- esdTrackCutsL->SetRequireTPCRefit(kTRUE);\r
- esdTrackCutsL->SetRequireITSRefit(kTRUE);\r
- esdTrackCutsL->SetMaxDCAToVertexXY(3.);\r
- esdTrackCutsL->SetMaxDCAToVertexZ(3.);\r
- esdTrackCutsL->SetAcceptKinkDaughters(kFALSE);\r
- taskKMeans->SetCuts(esdTrackCutsL);\r
- taskKMeans->SetK(4);\r
- taskKMeans->SetMinimumMultiplicity(10);\r
- AliKMeansClustering::SetBeta(1.);\r
- mgr->AddTask(taskKMeans);\r
-\r
- AliAnalysisDataContainer* cout_kmeans = mgr->CreateContainer("KMeans", TList::Class(),AliAnalysisManager::kOutputContainer,\r
- Form("%s:PWG4_KMeans", AliAnalysisManager::GetCommonFileName()));\r
-\r
- // Create ONLY the output containers for the data produced by the task.\r
- // Get and connect other common input/output containers via the manager as below\r
- //==============================================================================\r
- mgr->ConnectInput (taskKMeans, 0, mgr->GetCommonInputContainer());\r
- mgr->ConnectOutput (taskKMeans, 1, cout_kmeans);\r
-\r
- return taskKMeans;\r
-}\r
+AliAnalysisTaskKMeans *AddTaskKMeans()
+{
+// Creates a dijet task, configures it and adds it to the analysis manager.
+
+ // Get the pointer to the existing analysis manager via the static access method.
+ //==============================================================================
+ AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+ if (!mgr) {
+ ::Error("AddTaskJets", "No analysis manager to connect to.");
+ return NULL;
+ }
+
+ // Check the analysis type using the event handlers connected to the analysis manager.
+ //==============================================================================
+ if (!mgr->GetInputEventHandler()) {
+ ::Error("AddTaskKMeans", "This task requires an input event handler");
+ return NULL;
+ }
+
+ // Create the task and configure it.
+ //===========================================================================
+
+ AliAnalysisTaskKMeans *taskKMeans = new AliAnalysisTaskKMeans("K-Means Analysis");
+ taskKMeans->SetDebugLevel(0);
+ AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("AliESDtrackCuts", "Standard");
+ esdTrackCutsL->SetMinNClustersTPC(50);
+ esdTrackCutsL->SetRequireTPCRefit(kTRUE);
+ esdTrackCutsL->SetRequireITSRefit(kTRUE);
+ esdTrackCutsL->SetMaxDCAToVertexXY(3.);
+ esdTrackCutsL->SetMaxDCAToVertexZ(3.);
+ esdTrackCutsL->SetAcceptKinkDaughters(kFALSE);
+ taskKMeans->SetCuts(esdTrackCutsL);
+ taskKMeans->SetK(4);
+ taskKMeans->SetMinimumMultiplicity(10);
+ AliKMeansClustering::SetBeta(1.);
+ mgr->AddTask(taskKMeans);
+
+ AliAnalysisDataContainer* cout_kmeans = mgr->CreateContainer("KMeans", TList::Class(),AliAnalysisManager::kOutputContainer,
+ Form("%s:PWG4_KMeans", AliAnalysisManager::GetCommonFileName()));
+
+ // Create ONLY the output containers for the data produced by the task.
+ // Get and connect other common input/output containers via the manager as below
+ //==============================================================================
+ mgr->ConnectInput (taskKMeans, 0, mgr->GetCommonInputContainer());
+ mgr->ConnectOutput (taskKMeans, 1, cout_kmeans);
+
+ return taskKMeans;
+}