plotted delta pT for different cents

[u/mrichter/AliRoot.git] / PWG4 / macros / AddTaskKMeans.C

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