Merge remote-tracking branch 'origin/master' into TPCdev

[u/mrichter/AliRoot.git] / PWGLF / STRANGENESS / Correlations / macros / AddTaskLambdaOverK0sJets.C

AliAnalysisTaskLambdaOverK0sJets *AddTaskLambdaOverK0sJets( TString  name      = "LambdaOverK0sRatio", 
							    TString  data      = "PbPb2010", 
							    Float_t  minCen    = 0.,
							    Float_t  maxCen    = 40.,
							    Bool_t   sepInjec  = kTRUE,
							    Bool_t   isMC      = kFALSE,
							    Bool_t   doQA      = kTRUE){




  Float_t  ptMinTrig   = 5.;
  Float_t  ptMaxTrig   = 10.;
  Float_t  etaMaxTrig  = 0.7;
  Float_t  checkIDTrig = kTRUE;
  Float_t  rapMaxV0    = 0.7;
  Bool_t   usePID      = kFALSE;
  Float_t  nSigmaPID   = 3.0;
  Float_t  dcaDaug     = 0.095;  // looser cut
  Float_t  nclsDaug    = 50;     // looser cut


  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  if (!mgr) {
    ::Error("AddTaskCheckCascade", "No analysis manager to connect to.");
    return NULL;
  }   
  
  
  
  // Create and configure the task
  AliAnalysisTaskLambdaOverK0sJets *task = new AliAnalysisTaskLambdaOverK0sJets(name.Data());
  task->SetCollisionType(data);
  task->SetCentrality(minCen,maxCen);
  task->SetTriggerPt(ptMinTrig,ptMaxTrig);
  task->SetTriggerEta(etaMaxTrig);
  task->SetCheckIDTrig(checkIDTrig);
  task->SetMaxY(rapMaxV0);
  task->SetMaxDCADaughter(dcaDaug); // Add to perform systematics
  task->SetNSigmaPID(nSigmaPID);
  task->SetNClsTPC(nclsDaug);
  task->SetSeparateInjectedPart(sepInjec);
  task->SetMC(isMC);
  task->SetPID(usePID);
  task->SetQA(doQA);
  mgr->AddTask(task);
  
  
  
  // 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
  //==============================================================================
  
  AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();
  mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer());
  
  if (isMC) name+="_mc";
  
  AliAnalysisDataContainer *coutput1 =  
    mgr->CreateContainer(name, TList::Class(), 
                         AliAnalysisManager::kOutputContainer, 
                         "AnalysisResults.root");

  AliAnalysisDataContainer *coutput2 =  
    mgr->CreateContainer(name+"_ME", TList::Class(), 
                         AliAnalysisManager::kOutputContainer, 
                         "AnalysisResults.root");

  AliAnalysisDataContainer *coutput3 =  
    mgr->CreateContainer(name+"_QA", TList::Class(), 
                         AliAnalysisManager::kOutputContainer, 
                         "AnalysisResults.root");
  
  mgr->ConnectOutput(task,1,coutput1);
  mgr->ConnectOutput(task,2,coutput2);
  mgr->ConnectOutput(task,3,coutput3);

  
  return task;
}