]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGHF/vertexingHF/ConfigVertexingHF_highmult.C
Merge branch 'feature-movesplit'
[u/mrichter/AliRoot.git] / PWGHF / vertexingHF / ConfigVertexingHF_highmult.C
index 35f9de9af57df02ddc8a7d21512efb8a6ff3be60..ba14cc25b5d7ca4ddf18e9d29a6cb71b4cf46982 100644 (file)
@@ -8,6 +8,7 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   vHF->SetJPSItoEleOff();
   //vHF->Set3ProngOff();
   vHF->SetLikeSignOn(); // like-sign pairs and triplets
+  //  vHF->SetLikeSign3prongOff();
   vHF->Set4ProngOff();
   //vHF->SetDstarOff();
   vHF->SetFindVertexForDstar(kFALSE);
@@ -16,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,
@@ -76,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();
@@ -96,14 +98,12 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   cutsD0toKpipipi->SetCuts(9,cutsArrayD0toKpipipi);
   cutsD0toKpipipi->AddTrackCuts(esdTrackCuts);
   vHF->SetCutsD0toKpipipi(cutsD0toKpipipi);
-  AliRDHFCutsDStartoKpipi *cutsDStartoKpipi = new AliRDHFCutsDStartoKpipi("CutsDStartoKpipi");
-  cutsDStartoKpipi->SetStandardCutsPbPb2010();
-  cutsDStartoKpipi->SetUsePID(kFALSE);
 
 
  // D* pt dependent cuts ------------------------------------------
 
   AliRDHFCutsDStartoKpipi *cutsDStartoKpipi = new AliRDHFCutsDStartoKpipi("CutsDStartoKpipi");
+  cutsDStartoKpipi->SetUsePID(kFALSE);
   
   const Int_t nvars=16;
   const Int_t nptbins=2;
@@ -111,7 +111,7 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   Float_t* ptbins;
   ptbins=new Float_t[nptbins+1];
   ptbins[0]=0.;
-  ptbins[1]=4.;
+  ptbins[1]=5.;
   ptbins[2]=999.;
   
   cutsDStartoKpipi->SetPtBins(nptbins+1,ptbins);
@@ -121,54 +121,64 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   for(Int_t iv=0;iv<nvars;iv++){
     rdcutsvalmine[iv]=new Float_t[nptbins];
   }
-  //0-4
-  rdcutsvalmine[0][0]=0.15;
-  rdcutsvalmine[1][0]=0.07;
-  rdcutsvalmine[2][0]=0.8;
-  rdcutsvalmine[3][0]=0.8;
-  rdcutsvalmine[4][0]=0.8;
-  rdcutsvalmine[5][0]=0.08;
-  rdcutsvalmine[6][0]=0.08;
-  rdcutsvalmine[7][0]=-0.00002;
-  rdcutsvalmine[8][0]=0.72;
-  rdcutsvalmine[9][0]=0.15;
-  rdcutsvalmine[10][0]=0.03;
-  rdcutsvalmine[11][0]=0.2;
-  rdcutsvalmine[12][0]=100.;
-  rdcutsvalmine[13][0]=0.5;
-  rdcutsvalmine[14][0]=-1.;
-  rdcutsvalmine[15][0]=0.;
 //4-999
-  rdcutsvalmine[0][1]=0.24;
-  rdcutsvalmine[1][1]=0.07;
-  rdcutsvalmine[2][1]=0.8;
-  rdcutsvalmine[3][1]=0.8;
-  rdcutsvalmine[4][1]=0.8;
-  rdcutsvalmine[5][1]=0.1;
-  rdcutsvalmine[6][1]=0.1;
-  rdcutsvalmine[7][1]=0.00001;
-  rdcutsvalmine[8][1]=0.72;
-  rdcutsvalmine[9][1]=0.15;
-  rdcutsvalmine[10][1]=0.03;
-  rdcutsvalmine[11][1]=0.2;
-  rdcutsvalmine[12][1]=100.;
-  rdcutsvalmine[13][1]=0.5;
-  rdcutsvalmine[14][1]=-1.;
-  rdcutsvalmine[15][1]=0.;
+  //0-5
+  rdcutsvalmine[0][0]=0.10;  //D0 inv mass window
+  rdcutsvalmine[1][0]=0.06;   // dca
+  rdcutsvalmine[2][0]=0.9;    // thetastar
+  rdcutsvalmine[3][0]=0.5;    // pt Pion
+  rdcutsvalmine[4][0]=0.5;    // Pt Kaon
+  rdcutsvalmine[5][0]=0.1;    // d0K
+  rdcutsvalmine[6][0]=0.1;    // d0Pi
+  rdcutsvalmine[7][0]=0.0001; // d0xd0
+  rdcutsvalmine[8][0]=0.8;    // costhetapoint
+  rdcutsvalmine[9][0]=0.15;   // Dstar inv mass window
+  rdcutsvalmine[10][0]=0.03;  // half width of (M_Kpipi-M_D0)
+  rdcutsvalmine[11][0]=0.1;   // Pt min of Pi soft
+  rdcutsvalmine[12][0]=100.;  // Pt max of pi soft
+  rdcutsvalmine[13][0]=9999.; // theta
+  rdcutsvalmine[14][0]=0.9;   // |cosThetaPointXY|
+  rdcutsvalmine[15][0]=1.;    // NormDecayLenghtXY
//5-999
+  rdcutsvalmine[0][1]=0.10;   //D0 inv mass window
+  rdcutsvalmine[1][1]=0.06;   // dca
+  rdcutsvalmine[2][1]=0.9;    // thetastar
+  rdcutsvalmine[3][1]=0.5;    // pt Pion
+  rdcutsvalmine[4][1]=0.5;    // Pt Kaon
+  rdcutsvalmine[5][1]=0.1;    // d0K
+  rdcutsvalmine[6][1]=0.1;    // d0Pi
+  rdcutsvalmine[7][1]=0.0001; // d0xd0
+  rdcutsvalmine[8][1]=0.7;    // costhetapoint
+  rdcutsvalmine[9][1]=0.15;   // Dstar inv mass window
+  rdcutsvalmine[10][1]=0.03;  // half width of (M_Kpipi-M_D0)
+  rdcutsvalmine[11][1]=0.1;   // Pt min of Pi soft
+  rdcutsvalmine[12][1]=100.;  // Pt max of pi soft
+  rdcutsvalmine[13][1]=9999.; // theta
+  rdcutsvalmine[14][1]=0.8;   // |cosThetaPointXY|
+  rdcutsvalmine[15][1]=0.;    // NormDecayLenghtXY
 
   cutsDStartoKpipi->SetCuts(nvars,nptbins,rdcutsvalmine);
  
   cutsDStartoKpipi->AddTrackCuts(esdTrackCuts);
   cutsDStartoKpipi->AddTrackCutsSoftPi(esdTrackCutsSoftPi);
-  cutsDStartoKpipi->SetMinPtCandidate(3.);
+  cutsDStartoKpipi->SetMinPtCandidate(2.);
   vHF->SetCutsDStartoKpipi(cutsDStartoKpipi);
 
   //--------------------------------------------------------
 
   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