1 AliAnalysisTaskSE *AddTaskNuclei(Bool_t kAOD=kTRUE, Double_t CentralityMin=0.0, Double_t CentralityMax=100.0, Int_t filterBit=16, Double_t AbsEtaMin=0.0, Double_t AbsEtaMax=0.8, Double_t DCAxyCut=0.1, Double_t DCAzCut=1000.0, Double_t fNsigmaTpcCut=2.0, Int_t NminTpcCluster=0, Int_t iTRDslices=0, Int_t kSignalCheck=1, Int_t iMtof=1, Int_t kPvtxCorr=1){
3 //get the current analysis manager
4 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
7 Error("No manager found in AddTaskVZERO. Why?");
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!");
15 else if(1){ // check ESD
19 //========= Add tender to the ANALYSIS manager and set default storage =====
22 snprintf(mytaskName,100,"AliAnalysisNucleiMass");
24 AliAnalysisNucleiMass *task = new AliAnalysisNucleiMass(mytaskName);
25 task->SetCentrality(CentralityMin,CentralityMax);
26 task->SetFilterBit(filterBit);
27 task->SetAbsEtaLimit(AbsEtaMin,AbsEtaMax);
28 task->SetDCACut(DCAxyCut,DCAzCut);
29 task->SetNsigmaTPCCut(fNsigmaTpcCut);
30 task->SetNminTPCcluster(NminTpcCluster);
31 task->SetTrdCut(iTRDslices);
32 task->SetisSignalCheck(kSignalCheck);
33 task->SetMtofMethod(iMtof);
34 task->SetPvtxNucleiCorrection(kPvtxCorr);
38 //Attach input to my tasks
41 snprintf(name,1000,"cchain1%.0f_%.0f_FilterBit=%02i_EtaMin=%.1f_EtaMax=%.1f_DCAxyCUT=%.2f_DCAzCUT=%.1f_NsigTPCcut=%1.0f_NminTpcClusters=%03i_iTrdCut=%i_kSignCheck=%i_iMtof=%i_kPvtxCorr=%i",CentralityMin,CentralityMax,filterBit,AbsEtaMin,AbsEtaMax,DCAxyCut,DCAzCut,fNsigmaTpcCut,NminTpcCluster,iTRDslices,kSignalCheck,iMtof,kPvtxCorr);
43 AliAnalysisDataContainer *cinput = mgr->CreateContainer(name,TChain::Class(),AliAnalysisManager::kInputContainer);
44 mgr->ConnectInput(task,0,mgr->GetCommonInputContainer());
46 // Attach output to my tasks
48 snprintf(name,1000,"ResultsBmm_CC%.0f_%.0f_FilterBit=%02i_EtaMin=%.1f_EtaMax=%.1f_DCAxyCUT=%.2f_DCAzCUT=%.1f_NsigTPCcut=%1.0f_NminTpcClusters=%03i_iTrdCut=%i_kSignCheck=%i_iMtof=%i_kPvtxCorr=%i",CentralityMin,CentralityMax,filterBit,AbsEtaMin,AbsEtaMax,DCAxyCut,DCAzCut,fNsigmaTpcCut,NminTpcCluster,iTRDslices,kSignalCheck,iMtof,kPvtxCorr);
49 AliAnalysisDataContainer *cOutputL= mgr->CreateContainer(name,TList::Class(), AliAnalysisManager::kOutputContainer, AliAnalysisManager::GetCommonFileName());
50 mgr->ConnectOutput(task, 1, cOutputL);
52 snprintf(name,1000,"ResultsBpp_CC%.0f_%.0f_FilterBit=%02i_EtaMin=%.1f_EtaMax=%.1f_DCAxyCUT=%.2f_DCAzCUT=%.1f_NsigTPCcut=%1.0f_NminTpcClusters=%03i_iTrdCut=%i_kSignCheck=%i_iMtof=%i_kPvtxCorr=%i",CentralityMin,CentralityMax,filterBit,AbsEtaMin,AbsEtaMax,DCAxyCut,DCAzCut,fNsigmaTpcCut,NminTpcCluster,iTRDslices,kSignalCheck,iMtof,kPvtxCorr);
53 AliAnalysisDataContainer *cOutputL2= mgr->CreateContainer(name,TList::Class(), AliAnalysisManager::kOutputContainer, AliAnalysisManager::GetCommonFileName());
54 mgr->ConnectOutput(task, 2, cOutputL2);