1 AliAnalysisTaskV0ForRAA *AddTaskV0ForRAA(Bool_t anaPP=kFALSE, Int_t cent=0,Int_t centDet=1,Int_t centRange=0, Bool_t mcMode=kFALSE, Bool_t mcTruthMode=kFALSE,Bool_t onFly=kTRUE,Bool_t usePID=kFALSE,Double_t etaCut,Bool_t mcEtaCut,const Char_t * addname){
5 //--- get the current analysis manager ---//
6 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
8 Error("AddTask_V0ForRAA", "No analysis manager found.");
13 // -- check for ESD and MC ---//
14 Bool_t hasESD=kFALSE,hasMC=kFALSE;
15 AliESDInputHandler *esdH =
16 static_cast<AliESDInputHandler*>(mgr->GetInputEventHandler());
17 if (esdH) hasESD=kTRUE;
18 cout<<"ESD: "<<hasESD<<endl;
19 if(!hasESD) return NULL;
21 if(mcMode || mcTruthMode){
22 AliMCEventHandler *mcH = dynamic_cast<AliMCEventHandler*>
23 (AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler());
25 cout<<"MC: "<<hasMC<<endl;
26 if(!hasMC) return NULL;
30 //========= Add task to the ANALYSIS manager =====
31 TString cutsname = "AliESDtrackCutsV0ForRAA";
32 TString taskname = "V0ForRAA";
33 TString outname = "V0ForRAA";
41 cutsname +="_MCTruth";
42 taskname +="_MCTruth";
66 AliAnalysisTaskV0ForRAA *task = new AliAnalysisTaskV0ForRAA(taskname);
68 Double_t minPt=0.0;//15;
70 //--- esd track cuts V0 daughters ---//
72 AliESDtrackCuts* esdTrackCuts = new AliESDtrackCuts(cutsname);
73 // esdTrackCuts->SetMinNClustersTPC(70);
74 esdTrackCuts->SetMaxChi2PerClusterTPC(4);
75 esdTrackCuts->SetMinNCrossedRowsTPC(70);
76 esdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.8);
78 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
79 esdTrackCuts->SetRequireTPCRefit(kTRUE);
80 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
82 //esdTrackCuts->SetEtaRange(-0.9,+0.9);
83 task->SetESDTrackCuts(esdTrackCuts);
86 //--- analysis modes ---//
87 task->SetAnapp(anaPP);
88 task->SetMCMode(mcMode);
89 task->SetMCTruthMode(mcTruthMode);
91 //--- general cuts ---//
92 task->SetUseOnthefly(onFly);
93 task->SetUsePID(usePID,4.0,2100.0);
94 task->SetPrimVertexZCut(10.0,kTRUE);
95 task->SetCosOfPointingAngleK(0.99,1000.0);
96 task->SetCosOfPointingAngleL(0.99,1000.0);
97 task->SetRapidityCutMother(kTRUE,0.5);
98 task->SetDoEtaOfMCDaughtersCut(mcEtaCut,etaCut);
100 //--- centrality ---//
101 task->SetUseCentrality(centDet); // 0=off, 1=VZERO, 2=SPD
102 task->SetUseCentralityBin(cent); // Bin to be used 0,5,10,20,30,40,50,60,70,80,90,(100=SPDonly)
103 task->SetUseCentralityRange(centRange);
105 task->SelectCollisionCandidates(AliVEvent::kMB | AliVEvent::kCentral | AliVEvent::kSemiCentral );
110 //================================================
112 //================================================
113 // find input container
115 AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
117 AliAnalysisDataContainer *coutput1 =
118 mgr->CreateContainer(outname, TList::Class(),
119 AliAnalysisManager::kOutputContainer,Form("%s:simones", AliAnalysisManager::GetCommonFileName()));
121 //--- connect containers ---//
122 mgr->ConnectInput (task, 0, cinput );
123 mgr->ConnectOutput (task, 1, coutput1);
125 AliLog::SetClassDebugLevel("AliAnalysisTaskV0ForRAA",2);