Coverity

[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / ConfigVertexingHF_Pb_Cent2080.C

diff --git a/PWG3/vertexingHF/ConfigVertexingHF_Pb_Cent2080.C b/PWG3/vertexingHF/ConfigVertexingHF_Pb_Cent2080.C
index 1fc1a933f732ef4984c76778ee9281212ee32eec..8ae2e447a156f0a461cb1c1002ab7559c5bccaab 100644 (file)
--- a/PWG3/vertexingHF/ConfigVertexingHF_Pb_Cent2080.C
+++ b/PWG3/vertexingHF/ConfigVertexingHF_Pb_Cent2080.C
@@ -40,7 +40,7 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   esdTrackCutsSoftPi->SetMinNClustersITS(3);
   esdTrackCutsSoftPi->SetMaxDCAToVertexXY(1.);  
   esdTrackCutsSoftPi->SetMaxDCAToVertexZ(1.);
-  esdTrackCutsSoftPi->SetPtRange(0.2,1.e10);
+  esdTrackCutsSoftPi->SetPtRange(0.1,1.e10);
   esdTrackCutsSoftPi->SetEtaRange(-0.8,+0.8);
   AliAnalysisFilter *trkFilterSoftPi = new AliAnalysisFilter("trackFilterSoftPi");
   trkFilterSoftPi->AddCuts(esdTrackCutsSoftPi);
@@ -50,10 +50,11 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   AliRDHFCutsD0toKpi *cutsD0toKpi = new AliRDHFCutsD0toKpi("CutsD0toKpi");
   cutsD0toKpi->SetStandardCutsPbPb2010();
   cutsD0toKpi->SetMinPtCandidate(0.);
-  cutsD0toKpi->SetMinCentrality(20.);
-  cutsD0toKpi->SetMaxCentrality(80.);
   cutsD0toKpi->SetUsePID(kFALSE);
   cutsD0toKpi->SetUseSpecialCuts(kFALSE);
+  cutsD0toKpi->SetUsePhysicsSelection(kFALSE);
+  cutsD0toKpi->SetMinCentrality(20.);
+  cutsD0toKpi->SetMaxCentrality(80.);
   cutsD0toKpi->SetMaxVtxZ(1.e6);
   cutsD0toKpi->SetTriggerClass("");
   Float_t cutsArrayD0toKpi[11]={0.2,999999.,1.1,0.,0.,999999.,999999.,0.,0.5,-1,0.};
@@ -98,14 +99,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;
@@ -113,7 +112,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);
@@ -123,46 +122,46 @@ 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(0.);
+  cutsDStartoKpipi->SetMinPtCandidate(1.);
   vHF->SetCutsDStartoKpipi(cutsDStartoKpipi);
 
   //--------------------------------------------------------