PWGJE/macros/AddTaskJetBackgroundSubtract.C

   1 AliAnalysisTaskJetBackgroundSubtract *AddTaskJetBackgroundSubtract(TString sJetBranches,Int_t iSubtraction = 1,char* cToReplace = "B0",char* cReplaceMask = "B%d", char* cSuffix = ""){
   2
   3
   4
   5
   6   //get the current analysis manager
   7   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
   8   if (!mgr) {
   9     Error("AddTask_obusch_jets", "No analysis manager found.");
  10     return 0;
  11   }
  12
  13
  14
  15
  16   //========= Add task to the ANALYSIS manager =====
  17
  18   // TOkenize...
  19   TObjArray *objArr = sJetBranches.Tokenize(" ");
  20
  21
  22   TString cAdd;
  23   cAdd += Form(cReplaceMask,iSubtraction);
  24   if(cSuffix!="") cAdd += Form("_%s",cSuffix);
  25
  26   AliAnalysisTaskJetBackgroundSubtract *task = new AliAnalysisTaskJetBackgroundSubtract(Form("JetSubtract_%s",cAdd.Data()));
  27   for(int iJB = 0;iJB<objArr->GetEntries();iJB++){
  28     TObjString *ostr = (TObjString*)objArr->At(iJB);
  29     task->AddJetBranch(ostr->GetString().Data());
  30   }
  31   task->SetToReplace(cToReplace);
  32   task->SetReplacementMask(cReplaceMask);
  33   task->SetSubtractionMethod(iSubtraction);
  34   mgr->AddTask(task);
  35
  36
  37   //================================================
  38   //              data containers
  39   //================================================
  40   //            find input container
  41
  42   AliAnalysisDataContainer *cinput  = mgr->GetCommonInputContainer();
  43   AliAnalysisDataContainer *coutput = mgr->CreateContainer(
  44                                                            Form("pwgjeJetSubtract_%s",cAdd.Data()),
  45                                                            TList::Class(),
  46                                                            AliAnalysisManager::kOutputContainer,
  47                                                            Form("%s:PWGJE_JetSubtract_%s",AliAnalysisManager::GetCommonFileName(),cAdd.Data()));
  48
  49   mgr->ConnectInput(task,0,cinput );
  50   mgr->ConnectOutput(task,1,coutput);
  51
  52   return task;
  53 }