Protection for pool number overflow

[u/mrichter/AliRoot.git] / PWGHF / vertexingHF / ConfigVertexingHF_highmult.C

diff --git a/PWGHF/vertexingHF/ConfigVertexingHF_highmult.C b/PWGHF/vertexingHF/ConfigVertexingHF_highmult.C
index 27a4f5e83afd5722ea5b9f664a4d97e37a06d55e..ba14cc25b5d7ca4ddf18e9d29a6cb71b4cf46982 100644 (file)
--- a/PWGHF/vertexingHF/ConfigVertexingHF_highmult.C
+++ b/PWGHF/vertexingHF/ConfigVertexingHF_highmult.C
@@ -17,12 +17,13 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   vHF->SetCascadesOff();
   vHF->SetFindVertexForCascades(kFALSE);
   vHF->SetMassCutBeforeVertexing(kTRUE); // PbPb
+  vHF->SetV0TypeForCascadeVertex(AliRDHFCuts::kAllV0s);
 
   //--- set cuts for single-track selection  
   //     displaced tracks
   AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts","default");
   esdTrackCuts->SetRequireTPCRefit(kTRUE);
-  esdTrackCuts->SetMinNClustersTPC(70);
+  esdTrackCuts->SetMinNClustersTPC(50);
   esdTrackCuts->SetRequireITSRefit(kTRUE);
   //esdTrackCuts->SetMinNClustersITS(4);
   esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
@@ -77,11 +78,11 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   AliRDHFCutsDstoKKpi *cutsDstoKKpi = new AliRDHFCutsDstoKKpi("CutsDstoKKpi");
   cutsDstoKKpi->SetStandardCutsPbPb2010();
   cutsDstoKKpi->SetUsePID(kFALSE);
-  Float_t cutsArrayDstoKKpi[16]={0.2,0.4,0.4,0.,0.,0.005,0.045,0.,0.,0.9,0.,100000.,0.035,0.05,-1.,1.};
+  Float_t cutsArrayDstoKKpi[20]={0.2,0.4,0.4,0.,0.,0.005,0.045,0.,0.,0.9,0.,100000.,0.035,0.05,-1.,1.,0.,0.,0.,-1.};
   cutsDstoKKpi->SetPtBins(nptbins,ptlimits);
-  cutsDstoKKpi->SetCuts(16,cutsArrayDstoKKpi);
+  cutsDstoKKpi->SetCuts(20,cutsArrayDstoKKpi);
   cutsDstoKKpi->AddTrackCuts(esdTrackCuts);
-  cutsDstoKKpi->SetMinPtCandidate(4.);
+  cutsDstoKKpi->SetMinPtCandidate(3.);
   vHF->SetCutsDstoKKpi(cutsDstoKKpi);
   AliRDHFCutsLctopKpi *cutsLctopKpi = new AliRDHFCutsLctopKpi("CutsLctopKpi");
   cutsLctopKpi->SetStandardCutsPbPb2010();
@@ -165,9 +166,19 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   //--------------------------------------------------------
 
   AliRDHFCutsLctoV0 *cutsLctoV0 = new AliRDHFCutsLctoV0("CutsLctoV0");
-  Float_t cutsArrayLctoV0[9]={4.0,4.0,2.0,2.0,0.0,0.0,0.0,1000.,1000.};
-  cutsLctoV0->SetCuts(9,cutsArrayLctoV0);
+  Float_t cutsArrayLctoV0[17]={1.0,1.0,0.05,0.05,0.0,0.0,0.0,1000.,1000.,0.99,3.,1000.,0.,0.,0.,0.,0.0};
+  cutsLctoV0->SetCuts(17,cutsArrayLctoV0);
   cutsLctoV0->AddTrackCuts(esdTrackCuts);
+  AliESDtrackCuts *esdV0daughterTrackCuts = new AliESDtrackCuts("AliESDtrackCutsForV0D","default cuts for V0 daughters");
+  esdV0daughterTrackCuts->SetRequireTPCRefit(kTRUE);
+  esdV0daughterTrackCuts->SetMinNClustersTPC(30);
+  esdV0daughterTrackCuts->SetRequireITSRefit(kFALSE);
+  esdV0daughterTrackCuts->SetMinDCAToVertexXY(0.);
+  esdV0daughterTrackCuts->SetPtRange(0.05,1.e10);
+  esdV0daughterTrackCuts->SetEtaRange(-1.1,+1.1);
+  esdV0daughterTrackCuts->SetAcceptKinkDaughters(kTRUE);
+  esdV0daughterTrackCuts->SetRequireSigmaToVertex(kFALSE);
+  cutsLctoV0->AddTrackCutsV0daughters(esdV0daughterTrackCuts);
   vHF->SetCutsLctoV0(cutsLctoV0);
   // 
   //--- set this if you want to reconstruct primary vertex candidate by