prod/acrcaf/qa_pp/AddTaskQAsym.C

   1 void AddTaskQAsym(Int_t runNumber)\r
   2 {\r
   3   // Creates a QA task exploiting simple symmetries phi, eta +/-, charge ...\r
   4   \r
   5   // Get the pointer to the existing analysis manager via the static access method.\r
   6   //==============================================================================\r
   7   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
   8   if (!mgr) {\r
   9     ::Error("AddTaskQAsym", "No analysis manager to connect to.");\r
  10     return NULL;\r
  11   }  \r
  12   \r
  13   // Check the analysis type using the event handlers connected to the analysis manager.\r
  14   //==============================================================================\r
  15   if (!mgr->GetInputEventHandler()) {\r
  16     ::Error("AddTasQAsym", "This task requires an input event handler");\r
  17     return NULL;\r
  18   }\r
  19    TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"\r
  20   \r
  21    // Configure analysis\r
  22    //===========================================================================\r
  23    \r
  24  \r
  25  \r
  26    AliAnalysisTaskQASym *task = new AliAnalysisTaskQASym("AliAnalysisTaskQASym");\r
  27    AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("AliESDtrackCuts","Standard");\r
  28    esdTrackCutsL->SetMinNClustersTPC(50);\r
  29    esdTrackCutsL->SetRequireTPCRefit(kTRUE);\r
  30    esdTrackCutsL->SetMaxDCAToVertexXY(3.);\r
  31    esdTrackCutsL->SetMaxDCAToVertexZ(3.);\r
  32    esdTrackCutsL->SetAcceptKinkDaughters(kFALSE);\r
  33    \r
  34    task->SetCuts(esdTrackCutsL);\r
  35 \r
  36    mgr->AddTask(task);\r
  37   \r
  38    AliAnalysisDataContainer *cout  = mgr->CreateContainer("QAsymHists",TList::Class(),\r
  39                                                         AliAnalysisManager::kOutputContainer, Form("run%d.root",runNumber));\r
  40 \r
  41    mgr->ConnectInput  (task, 0, mgr->GetCommonInputContainer());\r
  42    mgr->ConnectOutput (task, 1, cout);\r
  43   \r
  44 }\r
  45 \r
  46 \r