PWGCF/Correlations/DPhi/MuonHadron/AddAnalysisTaskMuonHadronCorrelations.C

   1 AliAnalysisTaskMuonHadronCorrelations *AddAnalysisTaskMuonHadronCorrelations(const char *centMethod = "V0M") {
   2
   3   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
   4   if (!mgr) {
   5     printf("Error in adding AnalysisTaskMuonHadronCorrelations: no Analysis Manager found!\n");
   6     return NULL;
   7   }
   8
   9   AliAnalysisTaskMuonHadronCorrelations *task = new AliAnalysisTaskMuonHadronCorrelations(Form("AliAnalysisTaskMuonHadronCorrelations_%s",centMethod));
  10
  11   // Set analysis cuts
  12   task->SetFilterBitCentralBarrel(7);  // -> 128
  13   task->SetMaxEtaCentralBarrel(1.0);
  14   task->SetMinEtaCentralBarrel(-1.0);
  15   task->SetTriggerMatchLevelMuon(1);
  16
  17   const Int_t nBinCent = 4;
  18   Double_t centLimits[nBinCent+1] = {0., 20., 40, 60., 100.};
  19   task->SetCentBinning(nBinCent, centLimits);
  20
  21   task->SetCentMethod(centMethod);
  22
  23   const Int_t nBinPt = 3;
  24   Double_t ptLimits[nBinPt+1] = {0.5, 1., 2., 4.};
  25   task->SetPtBinning(nBinPt, ptLimits);
  26
  27   const Int_t nBinEta = 3;
  28   Double_t etaLimits[nBinEta+1] = {-4., -3.5, -3.0, -2.5};
  29   task->SetEtaBinning(nBinEta, etaLimits);
  30
  31   mgr->AddTask(task);
  32
  33   // create output container
  34   AliAnalysisDataContainer *output = mgr->CreateContainer(Form("MuonHadronCorrHistos_%s",centMethod), TList::Class(), AliAnalysisManager::kOutputContainer,
  35                                                           Form("%s:MuonHadronCorrelations_%s", AliAnalysisManager::GetCommonFileName(), centMethod));
  36
  37   // finaly connect input and output
  38   mgr->ConnectInput(task, 0,  mgr->GetCommonInputContainer());
  39   mgr->ConnectOutput(task, 1, output);
  40
  41   return task;
  42 }
  43