ANALYSIS/macros/AddTaskPhysicsSelection.C

   1 AliPhysicsSelectionTask* AddTaskPhysicsSelection(Bool_t mCAnalysisFlag = kFALSE, Bool_t deprecatedFlag = kTRUE, UInt_t computeBG = 0)
   2 {
   3   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
   4   if (!mgr) {
   5     ::Error("AddTaskPhysicsSelection", "No analysis manager to connect to.");
   6     return NULL;
   7   }
   8
   9   // Check the analysis type using the event handlers connected to the analysis manager.
  10   //==============================================================================
  11   if (!mgr->GetInputEventHandler()) {
  12     ::Error("AddTaskPhysicsSelection", "This task requires an input event handler");
  13     return NULL;
  14   }
  15   TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
  16
  17   // Configure analysis
  18   //===========================================================================
  19
  20
  21
  22   AliPhysicsSelectionTask *task = new AliPhysicsSelectionTask("");
  23   mgr->AddTask(task);
  24
  25   AliPhysicsSelection* physSel = task->GetPhysicsSelection();
  26   if (mCAnalysisFlag)
  27     physSel->SetAnalyzeMC();
  28   if (computeBG)
  29     physSel->SetComputeBG(computeBG);
  30
  31   if(!deprecatedFlag)
  32     AliFatal("The BG ID flag is deprecated. Please use the OADB to configure the cuts");
  33
  34   AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();
  35   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("cstatsout",
  36                 TList::Class(),
  37                 AliAnalysisManager::kOutputContainer,
  38                 "EventStat_temp.root");
  39
  40   mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer());
  41   mgr->ConnectOutput(task,1,coutput1);
  42
  43   mgr->RegisterExtraFile("event_stat.root");
  44
  45   return task;
  46 }