adding histogram of multiplicity vs centrality wo vertex selection (G.Luparello)

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

diff --git a/PWGHF/vertexingHF/ConfigVertexingHF_Pb_AllCent_NoLS_NoDstar_16var.C b/PWGHF/vertexingHF/ConfigVertexingHF_Pb_AllCent_NoLS_NoDstar_16var.C
index 1952efa3fa08eeec2c3534d782d1f8c1871d1c2f..23558285f9162eac3cd617fae0d8036da7f3051d 100644 (file)
--- a/PWGHF/vertexingHF/ConfigVertexingHF_Pb_AllCent_NoLS_NoDstar_16var.C
+++ b/PWGHF/vertexingHF/ConfigVertexingHF_Pb_AllCent_NoLS_NoDstar_16var.C
@@ -2,7 +2,7 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
 
   printf("Call to AliAnalysisVertexingHF parameters setting :\n");
   vHF = new AliAnalysisVertexingHF();
- 
+
   //--- switch-off candidates finding (default: all on)
   //vHF->SetD0toKpiOff();
   vHF->SetJPSItoEleOff();
@@ -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 
+  //--- 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,
@@ -40,14 +41,14 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   //     displaced tracks for 20% mot central events 2 prongs
   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,
                                         AliESDtrackCuts::kAny);
   // |d0|>75 micron for pt<2GeV, no cut above 2
   esdTrackCuts->SetMinDCAToVertexXYPtDep("0.0050*TMath::Max(0.,(1-TMath::Floor(TMath::Abs(pt)/2.)))");
-  esdTrackCuts->SetMaxDCAToVertexXY(1.);  
+  esdTrackCuts->SetMaxDCAToVertexXY(1.);
   esdTrackCuts->SetMaxDCAToVertexZ(1.);
   esdTrackCuts->SetPtRange(0.5,1.e10);
   esdTrackCuts->SetEtaRange(-0.8,+0.8);
@@ -58,14 +59,14 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   //     displaced tracks for 20% mot central events 3 prongs
   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,
                                         AliESDtrackCuts::kAny);
   // |d0|>75 micron for pt<2GeV, no cut above 2
   esdTrackCuts->SetMinDCAToVertexXYPtDep("0.0075*TMath::Max(0.,(1-TMath::Floor(TMath::Abs(pt)/2.)))");
-  esdTrackCuts->SetMaxDCAToVertexXY(1.);  
+  esdTrackCuts->SetMaxDCAToVertexXY(1.);
   esdTrackCuts->SetMaxDCAToVertexZ(1.);
   esdTrackCuts->SetPtRange(0.6,1.e10);
   esdTrackCuts->SetEtaRange(-0.8,+0.8);
@@ -76,7 +77,7 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
  //     D* soft pion tracks
   AliESDtrackCuts *esdTrackCutsSoftPi = new AliESDtrackCuts("AliESDtrackCuts","default");
   esdTrackCutsSoftPi->SetMinNClustersITS(3);
-  esdTrackCutsSoftPi->SetMaxDCAToVertexXY(1.);  
+  esdTrackCutsSoftPi->SetMaxDCAToVertexXY(1.);
   esdTrackCutsSoftPi->SetMaxDCAToVertexZ(1.);
   esdTrackCutsSoftPi->SetPtRange(0.1,1.e10);
   esdTrackCutsSoftPi->SetEtaRange(-0.8,+0.8);
@@ -91,7 +92,8 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   cutsD0toKpi->SetUsePID(kFALSE);
   cutsD0toKpi->SetUseSpecialCuts(kFALSE);
   cutsD0toKpi->SetUsePhysicsSelection(kFALSE);
-  cutsD0toKpi->SetMaxCentrality(80.);
+  cutsD0toKpi->SetMinCentrality(-10.);
+  cutsD0toKpi->SetMaxCentrality(110.);
   cutsD0toKpi->SetMaxVtxZ(1.e6);
   cutsD0toKpi->SetTriggerClass("");
   Float_t cutsArrayD0toKpi[11]={0.4,999999.,1.1,0.,0.,999999.,999999.,0.,0.5,-1,0.};
@@ -144,15 +146,15 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   cutsDStartoKpipi->SetUsePID(kFALSE);
   const Int_t nvars=16;
   const Int_t nptbins=2;
-  
+
   Float_t* ptbins;
   ptbins=new Float_t[nptbins+1];
   ptbins[0]=0.;
   ptbins[1]=5.;
   ptbins[2]=999.;
-  
+
   cutsDStartoKpipi->SetPtBins(nptbins+1,ptbins);
-  
+
   Float_t** rdcutsvalmine;
   rdcutsvalmine=new Float_t*[nvars];
   for(Int_t iv=0;iv<nvars;iv++){
@@ -194,7 +196,7 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   rdcutsvalmine[15][1]=0.;    // NormDecayLenghtXY
 
   cutsDStartoKpipi->SetCuts(nvars,nptbins,rdcutsvalmine);
- 
+
   cutsDStartoKpipi->AddTrackCuts(esdTrackCuts);
   cutsDStartoKpipi->AddTrackCutsSoftPi(esdTrackCutsSoftPi);
   cutsDStartoKpipi->SetMinPtCandidate(1.);
@@ -203,16 +205,27 @@ 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
-  //    candidate using other tracks in the event (for pp, broad 
+  //    candidate using other tracks in the event (for pp, broad
   //    interaction region)
   //vHF->SetRecoPrimVtxSkippingTrks();
-  //--- OR set this if you want to remove the candidate daughters from 
+  //--- OR set this if you want to remove the candidate daughters from
   //    the primary vertex, without recostructing it from scratch
   //vHF->SetRmTrksFromPrimVtx();
 
@@ -221,7 +234,7 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   //--- verbose
   //  AliLog::SetClassDebugLevel("AliAnalysisVertexingHF",2);
 
- 
+
   return vHF;
 }