PWGLF/STRANGENESS/Correlations/macros/AddTaskLambdaOverK0sJets.C

   1 AliAnalysisTaskLambdaOverK0sJets *AddTaskLambdaOverK0sJets( TString  name      = "LambdaOverK0sRatio",
   2                                                             TString  data      = "PbPb2010",
   3                                                             Float_t  minCen    = 0.,
   4                                                             Float_t  maxCen    = 40.,
   5                                                             Bool_t   sepInjec  = kTRUE,
   6                                                             Bool_t   isMC      = kFALSE,
   7                                                             Bool_t   doQA      = kTRUE){
   8
   9
  10
  11
  12   Float_t  ptMinTrig   = 5.;
  13   Float_t  ptMaxTrig   = 10.;
  14   Float_t  etaMaxTrig  = 0.7;
  15   Float_t  checkIDTrig = kTRUE;
  16   Float_t  rapMaxV0    = 0.7;
  17   Bool_t   usePID      = kFALSE;
  18   Float_t  nSigmaPID   = 3.0;
  19   Float_t  dcaDaug     = 1.0;
  20   Float_t  dca2PrmVtx  = 0.095;  // tighter cut
  21   Float_t  nclsDaug    = 50;     // looser cut
  22
  23
  24   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  25   if (!mgr) {
  26     ::Error("AddTaskCheckCascade", "No analysis manager to connect to.");
  27     return NULL;
  28   }
  29
  30
  31
  32   // Create and configure the task
  33   AliAnalysisTaskLambdaOverK0sJets *task = new AliAnalysisTaskLambdaOverK0sJets(name.Data());
  34   // collision type
  35   task->SetCollisionType(data);
  36   task->SetCentrality(minCen,maxCen);
  37   // trigger particle
  38   task->SetTriggerPt(ptMinTrig,ptMaxTrig);
  39   task->SetTriggerEta(etaMaxTrig);
  40   task->SetCheckIDTrig(checkIDTrig);
  41   // V0 candidates
  42   task->SetMaxY(rapMaxV0);
  43   task->SetMaxDCADaughter(dcaDaug); // Added to perform systematics
  44   task->SetDCAToPrimVtx(dca2PrmVtx); // Added to perform systematics
  45   task->SetNSigmaPID(nSigmaPID);
  46   task->SetNClsTPC(nclsDaug);  // Added to perform systematics
  47   // PID
  48   task->SetSeparateInjectedPart(sepInjec);
  49   // MC
  50   task->SetMC(isMC);
  51   task->SetPID(usePID);
  52   // QA
  53   task->SetQA(doQA);
  54   // Add task
  55   mgr->AddTask(task);
  56
  57
  58
  59   // Create ONLY the output containers for the data produced by the task.
  60   // Get and connect other common input/output containers via the manager as below
  61   //==============================================================================
  62
  63   AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();
  64   mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer());
  65
  66   if (isMC) name+="_mc";
  67
  68   AliAnalysisDataContainer *coutput1 =
  69     mgr->CreateContainer(name, TList::Class(),
  70                          AliAnalysisManager::kOutputContainer,
  71                          "AnalysisResults.root");
  72
  73   AliAnalysisDataContainer *coutput2 =
  74     mgr->CreateContainer(name+"_ME", TList::Class(),
  75                          AliAnalysisManager::kOutputContainer,
  76                          "AnalysisResults.root");
  77
  78   AliAnalysisDataContainer *coutput3 =
  79     mgr->CreateContainer(name+"_QA", TList::Class(),
  80                          AliAnalysisManager::kOutputContainer,
  81                          "AnalysisResults.root");
  82
  83   mgr->ConnectOutput(task,1,coutput1);
  84   mgr->ConnectOutput(task,2,coutput2);
  85   mgr->ConnectOutput(task,3,coutput3);
  86
  87
  88   return task;
  89 }