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