PWGLF/SPECTRA/Nuclei/masses/AddTaskNucleiInfo.C

   1 AliAnalysisTaskSE *AddTaskNucleiInfo(Bool_t kAOD=kTRUE, Int_t filterBit=16, Double_t EtaMin=-0.8, Double_t EtaMax=0.8, Double_t DCAxyCut=1000.0, Double_t DCAzCut=1000.0, Double_t fNsigmaTpcCut=2.0, Double_t StartTimeTofRes=9999.9){
   2
   3   //get the current analysis manager
   4   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
   5
   6   if (!mgr) {
   7     Error("No manager found in AddTaskVZERO. Why?");
   8     return 0;
   9   }
  10   // currently don't accept AOD input
  11   if (kAOD && !mgr->GetInputEventHandler()->InheritsFrom(AliAODInputHandler::Class())) { // check AOD
  12     Error("AddTaskVZERO","This task works only with AOD input!");
  13     return 0;
  14   }
  15   else if(1){ // check ESD
  16
  17   }
  18
  19   //========= Add tender to the ANALYSIS manager and set default storage =====
  20
  21   char mytaskName[100];
  22   snprintf(mytaskName,100,"AliAnalysisNucleiInfo");
  23
  24   AliAnalysisNucleiInfo *task = new AliAnalysisNucleiInfo(mytaskName);
  25   task->SetFilterBit(filterBit);
  26   task->SetEtaLimit(EtaMin,EtaMax);
  27   task->SetDCACut(DCAxyCut,DCAzCut);
  28   task->SetNsigmaTPCCut(fNsigmaTpcCut);
  29   task->SetStartTimeTofRes(StartTimeTofRes);
  30
  31   mgr->AddTask(task);
  32
  33   //Attach input to my tasks
  34   char name[1000];
  35
  36   snprintf(name,1000,"cchain1_FilterBit=%02i_EtaMin=%.1f_EtaMax=%.1f_DCAxyCUT=%.1f_DCAzCUT=%.1f_NsigTPCcut=%1.0f_StartTimeRes=%.1f",filterBit,EtaMin,EtaMax,DCAxyCut,DCAzCut,fNsigmaTpcCut,StartTimeTofRes);
  37
  38   AliAnalysisDataContainer *cinput = mgr->CreateContainer(name,TChain::Class(),AliAnalysisManager::kInputContainer);
  39   mgr->ConnectInput(task,0,mgr->GetCommonInputContainer());
  40
  41   // Attach output to my tasks
  42
  43   snprintf(name,1000,"ResultsBmm_FilterBit=%02i_EtaMin=%.1f_EtaMax=%.1f_DCAxyCUT=%.1f_DCAzCUT=%.1f_NsigTPCcut=%1.0f_StartTimeRes=%.1f",filterBit,EtaMin,EtaMax,DCAxyCut,DCAzCut,fNsigmaTpcCut,StartTimeTofRes);
  44   AliAnalysisDataContainer *cOutputL= mgr->CreateContainer(name,TList::Class(), AliAnalysisManager::kOutputContainer, AliAnalysisManager::GetCommonFileName());
  45   mgr->ConnectOutput(task, 1, cOutputL);
  46
  47   snprintf(name,1000,"ResultsBpp_FilterBit=%02i_EtaMin=%.1f_EtaMax=%.1f_DCAxyCUT=%.1f_DCAzCUT=%.1f_NsigTPCcut=%1.0f_StartTimeRes=%.1f",filterBit,EtaMin,EtaMax,DCAxyCut,DCAzCut,fNsigmaTpcCut,StartTimeTofRes);
  48   AliAnalysisDataContainer *cOutputL2= mgr->CreateContainer(name,TList::Class(), AliAnalysisManager::kOutputContainer, AliAnalysisManager::GetCommonFileName());
  49   mgr->ConnectOutput(task, 2, cOutputL2);
  50
  51   return task;
  52 }