X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PWGLF%2FSTRANGENESS%2FLambdaK0PbPb%2FAddTaskV0ForRAA.C;h=6148819a1ccf34ac71f8a6da132f84eda33ba49c;hb=11e769d963c5c14c7b8421846850c1bb4bd73efb;hp=f22adc0f047d1fbf82b66790cf8ac74f4aab220b;hpb=6cdc2a74638e82f5a663c5a7dce39ab63d3f6fc0;p=u%2Fmrichter%2FAliRoot.git diff --git a/PWGLF/STRANGENESS/LambdaK0PbPb/AddTaskV0ForRAA.C b/PWGLF/STRANGENESS/LambdaK0PbPb/AddTaskV0ForRAA.C index f22adc0f047..6148819a1cc 100644 --- a/PWGLF/STRANGENESS/LambdaK0PbPb/AddTaskV0ForRAA.C +++ b/PWGLF/STRANGENESS/LambdaK0PbPb/AddTaskV0ForRAA.C @@ -1,4 +1,4 @@ -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){ +AliAnalysisTaskV0ForRAA *AddTaskV0ForRAA(Bool_t anaPP=kFALSE, Bool_t wSDD=kFALSE,Int_t cent=0,Int_t centDet=1,Int_t centRange=0, Bool_t mcMode=kFALSE, Bool_t mcTruthMode=kFALSE,Bool_t usePID=kFALSE,Double_t radCut=0.0,const Char_t * addname=""){ @@ -58,48 +58,66 @@ AliAnalysisTaskV0ForRAA *AddTaskV0ForRAA(Bool_t anaPP=kFALSE, Int_t cent=0,Int_t outname += cent; } + cutsname += addname; + taskname += addname; + outname += addname; AliAnalysisTaskV0ForRAA *task = new AliAnalysisTaskV0ForRAA(taskname); - Double_t minPt=0.0;//15; - - //--- esd track cuts V0 daughters ---// + Double_t minPt=0.0; + task->SetESDTrackCuts(70,4,kTRUE); + task->SetESDTrackCutsCharged(70,4,kTRUE); + task->SetESDTrackCutsLowPt(70,4,kTRUE); - AliESDtrackCuts* esdTrackCuts = new AliESDtrackCuts(cutsname); - // esdTrackCuts->SetMinNClustersTPC(70); - esdTrackCuts->SetMaxChi2PerClusterTPC(4); - esdTrackCuts->SetMinNCrossedRowsTPC(70); - esdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.8); - - esdTrackCuts->SetAcceptKinkDaughters(kFALSE); - esdTrackCuts->SetRequireTPCRefit(kTRUE); - esdTrackCuts->SetRequireSigmaToVertex(kFALSE); - - //esdTrackCuts->SetEtaRange(-0.9,+0.9); - task->SetESDTrackCuts(esdTrackCuts); + //Add cuts to task + //--- analysis modes ---// task->SetAnapp(anaPP); task->SetMCMode(mcMode); task->SetMCTruthMode(mcTruthMode); - - //--- general cuts ---// - task->SetUseOnthefly(onFly); - task->SetUsePID(usePID,4.0,2100.0); + task->SelectWithSDD(wSDD); + //---------- cuts -------------// + //general cuts + task->SetUseOnthefly(kTRUE); + task->SetUsePID(usePID,3.0,100.0); task->SetPrimVertexZCut(10.0,kTRUE); + + //rapidity + task->SetRapidityCutMother(kTRUE,0.5); + task->SetDoEtaOfMCDaughtersCut(kFALSE,0.8); + + //TPC cuts + task->SetCutMoreNclsThanRows(kTRUE); + task->SetCutMoreNclsThanFindable(kTRUE); + task->SetLowPtTPCCutAliESDTrackCut(-1.0); + task->SetRatioFoundOverFindable(0.5); + + //V0 specific cuts task->SetCosOfPointingAngleK(0.99,1000.0); - task->SetCosOfPointingAngleL(0.99,1000.0); - task->SetRapidityCutMother(kTRUE,0.6); - task->SetDoEtaOfMCDaughtersCut(mcEtaCut,etaCut); + task->SetCosOfPointingAngleL(0.998,1000.0); + + + task->SetArmenterosCutQt(-1.0,6.0,kTRUE,kFALSE); + + task->SetDCAV0ToVertexK0(0.4); + task->SetDCAV0ToVertexL(1.2); + + task->SetDCADaughtersK0(0.23); + task->SetDCADaughtersL(0.35); + task->SetDCADaughtersAL(0.35); + + task->SetDecayRadiusXYMinMax(radCut,1000.0); + //--- centrality ---// - task->SetUseCentrality(centDet); // 0=off, 1=VZERO, 2=SPD - task->SetUseCentralityBin(cent); // Bin to be used 0,5,10,20,30,40,50,60,70,80,90,(100=SPDonly) - task->SetUseCentralityRange(centRange); + task->SetUseCentrality(centDet); // 0=off, 1=VZERO, 2=SPD + task->SetUseCentralityBin(cent); // bin to be used 0,5,10,20,30,40,50,60,70,80,90,(100=SPDonly) + task->SetUseCentralityRange(centRange); // add centrality bin for increasing original bin range. For cent 60-80%: cent = 60 and centRange = 10 - task->SelectCollisionCandidates(AliVEvent::kMB | AliVEvent::kCentral | AliVEvent::kSemiCentral ); + task->SelectCollisionCandidates(AliVEvent::kMB);// | AliVEvent::kCentral | AliVEvent::kSemiCentral ); mgr->AddTask(task); @@ -113,7 +131,7 @@ AliAnalysisTaskV0ForRAA *AddTaskV0ForRAA(Bool_t anaPP=kFALSE, Int_t cent=0,Int_t AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(outname, TList::Class(), - AliAnalysisManager::kOutputContainer,Form("%s:foldername", AliAnalysisManager::GetCommonFileName())); + AliAnalysisManager::kOutputContainer,Form("%s:simones", AliAnalysisManager::GetCommonFileName())); //--- connect containers ---// mgr->ConnectInput (task, 0, cinput );